In order to describe the molecular epidemiology of Blastocystis infection in Italy, 189 isolates, which had been collected during the years 2012–2014 from mildly symptomatic patients, or those affected by inflammatory bowel disease (IBD), irritable bowel syndrome (IBS) or chronic diarrhoea, or otherwise immunosuppressed, were subtyped by sequence analysis of the SSU rRNA gene (536 bp). Six subtypes (STs) were detected: ST1 (15·3%), ST2 (13·8%), ST3 (46·0%), ST4 (21·7%), ST6 (3·2%) and ST8 (0·5%). They clustered in distinct clades, as inferred from Bayesian inference phylogenetic and median joining network analyses. A high genetic differentiation was found at the inter-subtype level; it ranged from Jukes–Cantor (JC) distance = 0·02 (between ST1 and ST4) to JC = 0·11 (between ST6 and ST2). At the intra-ST level, a high genetic homogeneity was registered in ST4, whereas higher genetic variation was found in isolates corresponding to ST1 and ST2. Accordingly, high values of haplotype and nucleotide diversity were observed in ST1, ST2 and ST3. No association was found between patient gender and ST, whereas ST3 and ST1 were significantly more prevalent in patients aged 15–50 years. A significant occurrence of Blastocystis ST4 in patients suffering from IBS, IBD or chronic diarrhoea was observed; in addition, a slight significant association between ST1 and ST3 and IBS patients was found. Multiple correspondence analysis showed some significant contribution of different variables (subtypes, haplotypes, age) in the observed pattern of ordination of the 189 patients in the symptom categories.