Population genetic characterisation of dominant Cryptosporidium parvum subtype IIaA15G2R1

The subtype IIaA15G2R1 at the 60 kDa glycoprotein (gp60) gene locus is the most dominant Cryptosporidium parvum infecting dairy cattle and humans in industrialised nations. The reasons for its high transmissibility are not clear, and it remains to be determined whether this subtype represents a homogeneous parasite population. In this study, we sequence-characterised 26 IIaA15G2R subtype specimens and 26 non-IIaA15G2R subtype specimens from the United States, Canada, United Kingdom and Spain at seven other known polymorphic loci, including CP47, CP56, DZ-HRGP, MSC6-5, MSC6-7, RPGR and ZPT. Extensive heterogeneity within IIaA15G2R1 and discordance in typing results between gp60 and other genetic markers were observed. Results of inter-locus and intra-ZPT linkage disequilibrium and recombination analyses indicated that the heterogeneity within IIaA15G2R1 and discordance in typing results among genetic loci were largely due to the occurrence of genetic recombination, mostly within the gp60 subtype IIaA15G2R1. Although there was no clear population diversion between IIaA15G2R and non-IIaA15G2R subtypes, results of STRUCTURE and F_ST analyses suggested the presence of at least two subpopulations; subpopulation 1 had an epidemic population structure and was widely distributed, whereas subpopulation 2 had a clonal population structure and consisted of geographically segregated multilocus subtypes. Genetic recombination between epidemic and geographically segregated C. parvum populations appeared to be a driving force in the emergence of a hyper-transmissible IIaA15G2R1 subtype. Genetic recombination was observed even between the zoonotic IIa subtype family and anthroponotic subtype family IIc at CP56, MSC6-7 and ZPT. Thus, the IIaA15G2R1 subtype at gp60 is likely a fitness marker for C. parvum and the wide spread of IIaA15G2R1 subtype around the world is probably independent of the sequence characteristics at other genetic loci.