Extensive intra-host genetic diversity uncovered in Cryptosporidium parvum using Next Generation Sequencing

First identification of Cryptosporidium parvum virus 1 (CSpV1) in various subtypes of Cryptosporidium parvum from diarrheic calves, lambs and goat kids from France

Cryptosporidium spp. remain a major cause of waterborne diarrhea and illness in developing countries and represent a significant burden to farmers worldwide. Cryptosporidium parvum virus 1 (CSpV1), of the genus Cryspovirus, was first reported to be present in the cytoplasm of C. parvum in 1997. Full-length genome sequences have been obtained from C. parvum from Iowa (Iowa), Kansas (KSU) and China. We aimed at characterizing the genome of CSpV1 from France and used sequence analysis from Cryptosporidium isolates to explore whether CSpV1 genome diversity varies over time, with geographical sampling location, C. parvum genetic diversity, or ruminant host species. A total of 123 fecal samples of cattle, sheep and goats were collected from 17 different French departments (57 diseased animal fecal samples and 66 healthy animal fecal samples). Subtyping analysis of the C. parvum isolates revealed the presence of two zoonotic subtype families IIa and IId. Sequence analysis of CSpV1 revealed that all CSpV1 from France, regardless of the subtype of C. parvum (IIaA15G2R1, IIaA17G2R1 and IIdA18G1R1) are more closely related to CSpV1 from Turkey, and cluster on a distinct branch from CSpV1 collected from C. parvum subtype IIaA15G2R1 from Asia and North America. We also found that samples collected on a given year or successive years in a given location are more likely to host the same subtype of C. parvum and the same CSpV1 strain. Yet, there is no distinct clustering of CSpV1 per French department or ruminants, probably due to trade, and transmission of C. parvum among host species. Our results point towards (i) a close association between CSpV1 movement and C. parvum movement, (ii) recent migrations of C. parvum among distantly located departments and (iii) incidental transmission of C. parvum between ruminants. All together, these results provide insightful information regarding CSpV1 evolution and suggest the virus might be used as an epidemiological tracer for C. parvum. Future studies need to investigate CSpV1's role in C. parvum virulence and on subtype ability to infect different species.

Cross-Border Investigations on the Prevalence and Transmission Dynamics of Cryptosporidium Species in Dairy Cattle Farms in Western Mainland Europe

Cryptosporidium is an apicomplexan parasitic protist, which infects a wide range of hosts, causing cryptosporidiosis disease. In farms, the incidence of this disease is high in animals such as cows, leading to extensive economic loss in the livestock industry. Infected cows may also act as a major reservoir of Cryptosporidium spp., in particular C. parvum, the most common cause of cryptosporidiosis in these animals. This poses a risk to the trading of livestock, to other farms via breeding centres, and to human health. This study is a part of a global project aimed at strategies to tackle cryptosporidiosis. To reach this target, it was essential to determine whether prevalence was dependent on the studied countries or if the issue was borderless. Indeed, C. parvum occurrence was assessed across dairy farms in certain regions of Belgium, France, and the Netherlands. At the same time, the animal-to-animal transmission of the circulating C. parvum subtypes was studied. To accomplish this, we analysed 1084 faecal samples, corresponding to 57 dairy farms from all three countries. To this end, 18S rRNA and gp60 genes fragments were amplified, followed by DNA sequencing, which was subsequently used for detection and subtyping C. parvum. Bioinformatic and phylogenetic methods were integrated to analyse and characterise the obtained DNA sequences. Our results show 25.7%, 24.9% and 20.8% prevalence of Cryptosporidium spp. in Belgium, France, and the Netherlands respectively. Overall, 93% of the farms were Cryptosporidium positive. The gp60 subtyping demonstrated a significant number of the C. parvum positives belonged to the IIa allelic family, which has been also identified in humans. Therefore, this study highlights how prevalent C. parvum is in dairy farms and further suggests cattle as a possible carrier of zoonotic C. parvum subtypes, which could pose a threat to human health.

Taxonomy and molecular epidemiology of Cryptosporidium and Giardia - a 50 year perspective (1971-2021)

The protozoan parasites Cryptosporidium and Giardia are significant causes of diarrhoea worldwide and are responsible for numerous waterborne and foodborne outbreaks of diseases. Over the last 50 years, the development of improved detection and typing tools has facilitated the expanding range of named species. Currently at least 44 Cryptosporidium spp. and >120 genotypes, and nine Giardia spp., are recognised. Many of these Cryptosporidium genotypes will likely be described as species in the future. The phylogenetic placement of Cryptosporidium at the genus level is still unclear and further research is required to better understand its evolutionary origins. Zoonotic transmission has long been known to play an important role in the epidemiology of cryptosporidiosis and giardiasis, and the development and application of next generation sequencing tools is providing evidence for this. Comparative whole genome sequencing is also providing key information on the genetic mechanisms for host specificity and human infectivity, and will enable One Health management of these zoonotic parasites in the future.

Cryptosporidium spp. diagnosis and research in the 21st century

The protozoan parasite Cryptosporidium has emerged as a leading cause of diarrhoeal illness worldwide, posing a significant threat to young children and immunocompromised patients. While endemic in the vast majority of developing countries, Cryptosporidium also has the potential to cause waterborne epidemics and large scale outbreaks in both developing and developed nations. Anthroponontic and zoonotic transmission routes are well defined, with the ingestion of faecally contaminated food and water supplies a common source of infection. Microscopy, the current diagnostic mainstay, is considered by many to be suboptimal. This has prompted a shift towards alternative diagnostic techniques in the advent of the molecular era. Molecular methods, particularly PCR, are gaining traction in a diagnostic capacity over microscopy in the diagnosis of cryptosporidiosis, given the laborious and often tedious nature of the latter. Until now, developments in the field of Cryptosporidium detection and research have been somewhat hampered by the intractable nature of this parasite. However, recent advances in the field have taken the tentative first steps towards bringing Cryptosporidium research into the 21st century. Herein, we provide a review of these advances.

TIDE analysis of Cryptosporidium infections by gp60 typing reveals obscured mixed infections

BACKGROUND

Cryptosporidiosis is a parasitic disease associated with potentially fatal diarrhea. The most used method in Cryptosporidium subtyping is based on the glycoprotein gene gp60. Each infection can represent a parasite population, and it is important to investigate the influence on transmission and virulence, as well as any impact on public health investigations. However, an easy-to-use method for detection is lacking.

METHODS

Here we report on the use of the bioinformatic program TIDE for deconvolution of gp60 chromatograms. A combination of single oocyst analysis and cloning successfully confirmed the within-sample parasite population diversity. Retrospective sample analysis was conducted on archived chromatograms.

RESULTS

For C. parvum, 8.6% multi-strain infections (13 out of 152) obscured by currently used consensus base calling were detected. Importantly, we show that single oocysts can harbor a mixed population of sporozoites. We also identified a striking dominance of unappreciated polymerase stutter artefacts in all 218 chromatograms analyzed, challenging the uncritical use of gp60 typing.

DISCUSSION

We demonstrate the value of a new easy-to-use analytical procedure for critical characterization of C. parvum and C. hominis in epidemiological investigations, also applicable in retrospect. Our findings illuminate the hidden parasite diversity with important implications for tracing zoonotic and person-to-person transmissions.

An Illumina MiSeq-Based Amplicon Sequencing Method for the Detection of Mixed Parasite Infections Using the Blastocystis SSU rRNA Gene as an Example

Parasite mixed infections remain a relatively unexplored field in part due to the difficulties of unraveling complex mixtures of parasite DNA using classical methods of sequencing. Next-generation amplicon sequencing (NGS) is a powerful tool for exploring mixed infections of multiple genetic variants of the same parasite in clinical, environmental (water or soil), or food samples. Here, we provide a method for NGS-based detection of mixed parasite infections which uses the Blastocystis SSU rRNA gene as an example and includes steps for parasite concentration, DNA extraction, sequencing library preparation, and bioinformatic analysis.

Molecular analysis of cryptosporidiosis cases in Western Australia in 2019 and 2020 supports the occurrence of two swimming pool associated outbreaks and reveals the emergence of a rare C. hominis IbA12G3 subtype

Cryptosporidium is an important protozoan parasite and due to its resistance to chlorine is a major cause of swimming pool-associated gastroenteritis outbreaks. The present study combined contact tracing and molecular techniques to analyse cryptosporidiosis cases and outbreaks in Western Australia in 2019 and 2020. In the 2019 outbreak, subtyping at the 60 kDa glycoprotein (gp60) gene identified 89.0% (16/18) of samples were caused by the C. hominis IdA15G1 subtype. Amplicon next generation sequencing (NGS) at the gp60 locus identified five C. hominis IdA15G1 subtype samples that also had C. hominis IdA14 subtype DNA, while multi locus sequence typing (MLST) analysis on a subset (n = 14) of C. hominis samples identified three IdA15G1 samples with a 6 bp insertion at the end of the trinucleotide repeat region of the cp47 gene. In 2020, 88.0% (73/83) of samples typed were caused by the relatively rare C. hominis subtype IbA12G3. Four mixed infections were observed by NGS with three IdA15G1/ IdA14 mixtures and one C. parvum IIaA18G3R1 sample mixed with IIaA16G3R1. No genetic diversity using MLST was detected. Epidemiological and molecular data indicates that the outbreaks in 2019 and 2020 were each potentially from swimming pool point sources and a new C. hominis subtype IbA12G3 is emerging in Australia. The findings of the present study are important for understanding the introduction and transmission of rare Cryptosporidium subtypes to vulnerable populations.

Assessment of next generation amplicon sequencing of the beta-giardin gene for the detection of Giardia duodenalis assemblages and mixed infections

Giardia duodenalis is an enteric protozoan parasite commonly found in humans and many other animals around the world. The parasite is grouped into genetically related strains called assemblages which display differing degrees of host specificity. Although mixed assemblage infections have been documented the full extent of the occurrence and importance of mixed infections remains to be characterized as current sequencing technologies lack the sensitivity to readily detect mixed infections. Here we have developed a next generation amplicon sequencing (NGS) protocol and analysis pipeline for detecting Giardia assemblages using the beta-giardin gene. NGS was validated using 37 isolates that included Giardia muris and six assemblages (A-F) of Giardia duodenalis obtained from seven different hosts. NGS was compared to traditional PCR and direct Sanger sequencing for its ability to detect Giardia species, assemblages, and mixed assemblage infections. We demonstrate that NGS works as well as PCR and Sanger sequencing for assemblage detection as the same assemblage was observed in all samples by both methods. NGS has the further benefit of detecting mixed assemblage infections, low abundance assemblages, and intra-assemblage variation in samples which would have been missed using direct Sanger sequencing alone. NGS represents a powerful new tool for exploring Giardia infections not only in infected hosts but also in environmental specimens which may aide in understanding Giardia epidemiology.

•

Developed a next generation sequencing (NGS) protocol to detect Giardia

•

Comparison of NGS to Sanger sequencing for sensitivity and accuracy

•

NGS improved detection of mixed assemblages and detects low abundance assemblages

•

NGS allowed detection of intra-assemblage variabilities

•

Mixed assemblages may be far more common than previously thought

Direct Sequencing of Cryptosporidium in Stool Samples for Public Health

The protozoan parasite Cryptosporidium is an important cause of diarrheal disease (cryptosporidiosis) in humans and animals, with significant morbidity and mortality especially in severely immunocompromised people and in young children in low-resource settings. Due to the sexual life cycle of the parasite, transmission is complex. There are no restrictions on sexual recombination between sub-populations, meaning that large-scale genetic recombination may occur within a host, potentially confounding epidemiological analysis. To clarify the relationships between infections in different hosts, it is first necessary to correctly identify species and genotypes, but these differentiations are not made by standard diagnostic tests and more sophisticated molecular methods have been developed. For instance, multilocus genotyping has been utilized to differentiate isolates within the major human pathogens, Cryptosporidium parvum and Cryptosporidium hominis. This has allowed mixed populations with multiple alleles to be identified: recombination events are considered to be the driving force of increased variation and the emergence of new subtypes. As yet, whole genome sequencing (WGS) is having limited impact on public health investigations, due in part to insufficient numbers of oocysts and purity of DNA derived from clinical samples. Moreover, because public health agencies have not prioritized parasites, validation has not been performed on user-friendly data analysis pipelines suitable for public health practitioners. Nonetheless, since the first whole genome assembly in 2004 there are now numerous genomes of human and animal-derived cryptosporidia publically available, spanning nine species. It has also been demonstrated that WGS from very low numbers of oocysts is possible, through the use of amplification procedures. These data and approaches are providing new insights into host-adapted infectivity, the presence and frequency of multiple sub-populations of Cryptosporidium spp. within single clinical samples, and transmission of infection. Analyses show that although whole genome sequences do indeed contain many alleles, they are invariably dominated by a single highly abundant allele. These insights are helping to better understand population structures within hosts, which will be important to develop novel prevention strategies in the fight against cryptosporidiosis.

Shedding of Cryptosporidium in calves and dams: evidence of re-infection and shedding of different gp60 subtypes

One of the most common causes of calf diarrhoea is the parasite Cryptosporidium parvum. Two longitudinal studies were carried out on a dairy farm Scotland to determine the prevalence of Cryptosporidium species and subtypes in a group of calves and to determine whether dams were a possible source of calfhood infection. Fecal samples were collected from 25 calves from birth to 12 months in the first year. In the second year, fecal samples were collected from pregnant cows (n = 29) and their calves (n = 30) from birth to 6 months. The samples were tested for Cryptosporidium and speciated. Cryptosporidium parvum-positive samples were subtyped by GP60 fragment analysis. All calves in both studies shed Cryptosporidium during the study period. Cryptosporidium parvum was the predominant species detected in calves ⩽6 weeks of age and at 6 months of age, C. bovis and C. ryanae were detected in calves older than 4 weeks of age but ⩽6 months of age. The prevalence of Cryptosporidium was higher in younger animals than in older animals. GP60 subtyping revealed two subtypes in calves on this farm (IIaA15G2R1 and IIaA19G2R1) that differed in frequency by age. Adult cattle also shed C. parvum, of four gp60 genotypes.

Isolation, genotyping and subtyping of single Cryptosporidium oocysts from calves with special reference to zoonotic significance

The ability of the small-subunit ribosomal RNA (SSU rRNA) based nested PCR and Restriction Fragment Length Polymorphism (PCR-RFLP) to identify and genotype a single Cryptosporidium oocyst isolated from bovine faecal samples was evaluated in this study. In addition, subtyping was carried out by sequencing the 60 kDa glycoprotein (gp60) gene from the same single oocyst. Faecal samples were collected from 40 pre-weaned calves (5-20 days old) from 7 dairy farms located in 3 different counties within the Finger Lakes region of Upstate New York. All the samples were microscopically positive for Cryptosporidium spp. A total of 400 Cryptosporidium oocysts (10 single oocysts from each calf sample) were individually isolated and analyzed using a nested PCR targeting the SSU rRNA gene. The SSU rRNA gene was amplified in 324 (81%) individual oocysts. All SSU rRNA amplified individual oocysts DNA was genotyped using PCR-RFLP. C. parvum was the only identified species; 107 single oocysts generated PCR products from the A gene, 18 generated PCR products from the B gene and 199 generated PCR products from both. Sequence analysis of the gp60 gene in 99 individual oocysts revealed the presence of only subtype IIaA15G2R1 with 99.4-100% and 99.1-100% identity of nucleotides and amino acids, respectively. These sequences were identical (100%) in oocysts from 35 calves and exhibited mutations in the non-repeat region of the gp60 gene in those of 5 other calves. The examination of DNA from individual oocysts with genotyping and subtyping tools provides methodology to more clearly define the genetic characteristics of Cryptosporidium spp. on farms and within individual animals.

Genomics and molecular epidemiology of Cryptosporidium species

Cryptosporidium is one of the most widespread protozoan parasites that infects domestic and wild animals and is considered the second major cause of diarrhea and death in children after rotavirus. So far, around 20 distinct species are known to cause severe to moderate infections in humans, of which Cryptosporidium hominis and Cryptosporidium parvum are the major causative agents. Currently, ssurRNA and gp60 are used as the optimal markers for differentiating species and subtypes respectively. Over the last decade, diagnostic tools to detect and differentiate Cryptosporidium species at the genotype and subtype level have improved, but our understanding of the zoonotic and anthroponotic transmission potential of each species is less clear, largely because of the paucity of high resolution whole genome sequencing data for the different species. Defining which species possess an anthroponotic vs. zoonotic transmission cycle is critical if we are to limit the spread of disease between animals and humans. Likewise, it is unclear to what extent genetic hybridization impacts disease potential or the emergence of outbreak strains. The development of high resolution genetic markers and whole genome sequencing of different species should provide new insights into these knowledge gaps. The aim of this review is to outline currently available molecular epidemiology and genomics data for different species of Cryptosporidium.

Identification of pig-specific Cryptosporidium species in mixed infections using Illumina sequencing technology

Nowadays molecular methods are widely used in epidemiological studies of Cryptosporidium infections in humans and animals. However to gain better understanding of parasite species or genotypes, especially when mixed infections are noticed, highly sensitive tools with adequate resolution power need to be employed. In this article, we report an application of the next generation sequencing method (NGS) for detection and characterisation of Cryptosporidium species concurrently present in pig faeces. A mixture of Cryptosporidium DNA obtained from two faecal samples was amplified at the 18 SSU rRNA gene locus and the resulting amplicons were subsequently used for MiSeq sequencing. Although initial molecular analyses indicated the possible presence of another Cryptosporidium species other than Cryptosporidium scrofarum and Cryptosporidium suis, deep sequencing only confirmed the presence of pig-specific Cryptosporidium.

Increased genetic diversity and prevalence of co-infection with Trypanosoma spp. in koalas (Phascolarctos cinereus) and their ticks identified using next-generation sequencing (NGS)

Infections with Trypanosoma spp. have been associated with poor health and decreased survival of koalas (Phascolarctos cinereus), particularly in the presence of concurrent pathogens such as Chlamydia and koala retrovirus. The present study describes the application of a next-generation sequencing (NGS)-based assay to characterise the prevalence and genetic diversity of trypanosome communities in koalas and two native species of ticks (Ixodes holocyclus and I. tasmani) removed from koala hosts. Among 168 koalas tested, 32.2% (95% CI: 25.2–39.8%) were positive for at least one Trypanosoma sp. Previously described Trypanosoma spp. from koalas were identified, including T. irwini (32.1%, 95% CI: 25.2–39.8%), T. gilletti (25%, 95% CI: 18.7–32.3%), T. copemani (27.4%, 95% CI: 20.8–34.8%) and T. vegrandis (10.1%, 95% CI: 6.0–15.7%). Trypanosoma noyesi was detected for the first time in koalas, although at a low prevalence (0.6% 95% CI: 0–3.3%), and a novel species (Trypanosoma sp. AB-2017) was identified at a prevalence of 4.8% (95% CI: 2.1–9.2%). Mixed infections with up to five species were present in 27.4% (95% CI: 21–35%) of the koalas, which was significantly higher than the prevalence of single infections 4.8% (95% CI: 2–9%). Overall, a considerably higher proportion (79.7%) of the Trypanosoma sequences isolated from koala blood samples were identified as T. irwini, suggesting this is the dominant species. Co-infections involving T. gilletti, T. irwini, T. copemani, T. vegrandis and Trypanosoma sp. AB-2017 were also detected in ticks, with T. gilletti and T. copemani being the dominant species within the invertebrate hosts. Direct Sanger sequencing of Trypanosoma 18S rRNA gene amplicons was also performed and results revealed that this method was only able to identify the genotypes with greater amount of reads (according to NGS) within koala samples, which highlights the advantages of NGS in detecting mixed infections. The present study provides new insights on the natural genetic diversity of Trypanosoma communities infecting koalas and constitutes a benchmark for future clinical and epidemiological studies required to quantify the contribution of trypanosome infections on koala survival rates.

Cryptosporidium within-host genetic diversity: systematic bibliographical search and narrative overview

Knowledge of the within-host genetic diversity of a pathogen often has broad implications for disease management. Cryptosporidium protozoan parasites are among the most common causative agents of infectious diarrhoea. Current limitations of in vitro culture impose the use of uncultured isolates obtained directly from the hosts as operational units of Cryptosporidium genotyping. The validity of this practice is centred on the assumption of genetic homogeneity of the parasite within the host, and genetic studies often take little account of the within-host genetic diversity of Cryptosporidium. Yet, theory and experimental evidence contemplate genetic diversity of Cryptosporidium at the within-host scale, but this diversity is not easily identified by genotyping methods ill-suited for the resolution of DNA mixtures. We performed a systematic bibliographical search of the occurrence of within-host genetic diversity of Cryptosporidium parasites in epidemiological samples, between 2005 and 2015. Our results indicate that genetic diversity at the within-host scale, in the form of mixed species or intra-species diversity, has been identified in a large number (n=55) of epidemiological surveys of cryptosporidiosis in variable proportions, but has often been treated as a secondary finding and not analysed. As in malaria, there are indications that the scale of this diversity varies between geographical regions, perhaps depending on the prevailing transmission pathways. These results provide a significant knowledge base from which to draw alternative population genetic structure models, some of which are discussed in this paper.

Evaluation of next generation sequencing for the analysis of Eimeria communities in wildlife

Next-generation sequencing (NGS) techniques are well-established for studying bacterial communities but not yet for microbial eukaryotes. Parasite communities remain poorly studied, due in part to the lack of reliable and accessible molecular methods to analyse eukaryotic communities. We aimed to develop and evaluate a methodology to analyse communities of the protozoan parasite Eimeria from populations of the Australian marsupial Petrogale penicillata (brush-tailed rock-wallaby) using NGS. An oocyst purification method for small sample sizes and polymerase chain reaction (PCR) protocol for the 18S rRNA locus targeting Eimeria was developed and optimised prior to sequencing on the Illumina MiSeq platform. A data analysis approach was developed by modifying methods from bacterial metagenomics and utilising existing Eimeria sequences in GenBank. Operational taxonomic unit (OTU) assignment at a high similarity threshold (97%) was more accurate at assigning Eimeria contigs into Eimeria OTUs but at a lower threshold (95%) there was greater resolution between OTU consensus sequences. The assessment of two amplification PCR methods prior to Illumina MiSeq, single and nested PCR, determined that single PCR was more sensitive to Eimeria as more Eimeria OTUs were detected in single amplicons. We have developed a simple and cost-effective approach to a data analysis pipeline for community analysis of eukaryotic organisms using Eimeria communities as a model. The pipeline provides a basis for evaluation using other eukaryotic organisms and potential for diverse community analysis studies.

Intra-Species Diversity and Panmictic Structure of Cryptosporidium parvum Populations in Cattle Farms in Northern Spain

The intra-herd and intra-host genetic variability of 123 Cryptosporidium parvum isolates was investigated using a multilocus fragment typing approach with eleven variable-number tandem-repeat (VNTR) loci and the GP60 gene. Isolates were collected from intensively farmed diarrheic pre-weaned calves originating from 31 dairy farms in three adjoining regions in northern Spain (País Vasco, Cantabria and Asturias). The multilocus tool demonstrated an acceptable typeability, with 104/123 samples amplifying at all twelve loci. The ML2, TP14, GP60 and the previously un-described minisatellite at locus cgd2_3850 were the most discriminatory markers, while others may be dismissed as monomorphic (MSB) or less informative (CP47, ML1 and the novel minisatellites at loci Cgd1_3670 and Cgd6_3940). The 12-satellite typing tool provided a Hunter-Gaston index (HGDI) of 0.987 (95% CI, 0.982-0.992), and differentiated a total of 70 multilocus subtypes (MLTs). The inclusion of only the four most discriminatory markers dramatically reduced the number of MLTs (n: 44) but hardly reduced the HGDI value. A total of 54 MLTs were distinctive for individual farms, indicating that cryptosporidiosis is an endemic condition on most cattle farms. However, a high rate of mixed infections was detected, suggesting frequent meiotic recombination. Namely, multiple MLTs were seen in most farms where several specimens were analyzed (90.5%), with up to 9 MLTs being found on one farm, and individual specimens with mixed populations being reported on 11/29 farms. Bayesian Structure analysis showed that over 35% of isolates had mixed ancestry and analysis of evolutionary descent using the eBURST algorithm detected a high rate (21.4%) of MLTs appearing as singletons, indicating a high degree of genetic divergence. Linkage analysis found evidence of linkage equilibrium and an overall panmictic structure within the C. parvum population in this discrete geographical area.

Cross sectional study of prevalence, genetic diversity and zoonotic potential of Cryptosporidium parvum cycling in New Zealand dairy farms

BACKGROUND

The estimation of the prevalence and zoonotic potential of Cryptosporidium parvum cycling in bovine populations requires the use of genotyping, as several morphologically similar non-parvum genetic variants of unproven clinical and public health impact are found in cattle. However, robust C. parvum prevalence estimates in cattle are lacking and comparative data of bovine and human isolates collected from the same regions are scarce. Thus, the relative contribution of the C. parvum oocysts released by farmed animals to animal and human cryptosporidiosis burden is, in general, poorly understood.

METHODS

The New Zealand farm-level C. parvum prevalence was estimated using a cross-sectional sample of 1283 faecal specimens collected from newborn calves on 97 dairy farms. Faeces were analysed by immunofluorescence and the Cryptosporidium parasites were genetically identified. Finally, bovine C. parvum were genetically compared with historical human clinical isolates using a bilocus subtyping scheme.

RESULTS

Immunofluoresence-positive faeces were found in 63/97 (65%) farms. C. parvum was identified in 49 (50.5%) farms, C. bovis in 6 (6.1%) farms, and on 8 (8.2%) farms the species could not be identified. The dominant C. parvum genetic variants were geographically widespread and found in both host populations, but several variants were found in humans only.

CONCLUSIONS

Phenotypic tests offered by New Zealand veterinary diagnostic laboratories for the diagnosis of C. parvum may have moderate to high positive predictive values for this species. The genetic similarities observed between the human and bovine parasites support a model considering calves as significant amplifiers of zoonotic C. parvum in New Zealand. However, data suggest that transmission routes not associated with dairy cattle should also be taken into account in future source-attribution studies of human cryptosporidiosis.

Comparative genomic analysis reveals occurrence of genetic recombination in virulent Cryptosporidium hominis subtypes and telomeric gene duplications in Cryptosporidium parvum

Background

Cryptosporidium hominis is a dominant species for human cryptosporidiosis. Within the species, IbA10G2 is the most virulent subtype responsible for all C. hominis–associated outbreaks in Europe and Australia, and is a dominant outbreak subtype in the United States. In recent yearsIaA28R4 is becoming a major new subtype in the United States. In this study, we sequenced the genomes of two field specimens from each of the two subtypes and conducted a comparative genomic analysis of the obtained sequences with those from the only fully sequenced Cryptosporidium parvum genome.

Results

Altogether, 8.59-9.05 Mb of Cryptosporidium sequences in 45–767 assembled contigs were obtained from the four specimens, representing 94.36-99.47% coverage of the expected genome. These genomes had complete synteny in gene organization and 96.86-97.0% and 99.72-99.83% nucleotide sequence similarities to the published genomes of C. parvum and C. hominis, respectively. Several major insertions and deletions were seen between C. hominis and C. parvum genomes, involving mostly members of multicopy gene families near telomeres. The four C. hominis genomes were highly similar to each other and divergent from the reference IaA25R3 genome in some highly polymorphic regions. Major sequence differences among the four specimens sequenced in this study were in the 5′ and 3′ ends of chromosome 6 and the gp60 region, largely the result of genetic recombination.

Conclusions

The sequence similarity among specimens of the two dominant outbreak subtypes and genetic recombination in chromosome 6, especially around the putative virulence determinant gp60 region, suggest that genetic recombination plays a potential role in the emergence of hyper-transmissible C. hominis subtypes. The high sequence conservation between C. parvum and C. hominis genomes and significant differences in copy numbers of MEDLE family secreted proteins and insulinase-like proteases indicate that telomeric gene duplications could potentially contribute to host expansion in C. parvum.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1517-1) contains supplementary material, which is available to authorized users.

Multiple Cryptosporidium parvum subtypes detected in a unique isolate of a Chilean neonatal calf with diarrhea

To further understand the composition of population of parasite in a single host, we analyzed the GP60 gene of Cryptosporidium parvum amplified from DNA of a randomly selected isolate found in the feces of a diarrheic calf from a dairy farm in Central Chile. Direct sequencing of the amplicon yield the IIaA17G4R1 C. parvum subtype. The same amplicon was cloned in Escherichia coli (22 clones) and sequenced, yielding three different GP60 subtypes, IIaA17G4R1 (16/22), IIaA16G4R1 (1/22), and IIaA15G4R1 (1/22), and four sequences with nucleotide substitutions in the serine repeats, which subtype would be otherwise IIaA17G4R1. It is thus possible to determine allelic polymorphism using Sanger sequencing with an additional step of bacterial cloning. The results also indicate the necessity to further characterize parasite populations in a single host to better understand the dynamics of Cryptosporidium epidemiology.

Cryptosporidiumspecies in humans and animals: current understanding and research needs

Cryptosporidiumis increasingly recognized as one of the major causes of moderate to severe diarrhoea in developing countries. With treatment options limited, control relies on knowledge of the biology and transmission of the members of the genus responsible for disease. Currently, 26 species are recognized as valid on the basis of morphological, biological and molecular data. Of the nearly 20Cryptosporidiumspecies and genotypes that have been reported in humans,Cryptosporidium hominisandCryptosporidium parvumare responsible for the majority of infections. Livestock, particularly cattle, are one of the most important reservoirs of zoonotic infections. Domesticated and wild animals can each be infected with severalCryptosporidiumspecies or genotypes that have only a narrow host range and therefore have no major public health significance. Recent advances in next-generation sequencing techniques will significantly improve our understanding of the taxonomy and transmission ofCryptosporidiumspecies, and the investigation of outbreaks and monitoring of emerging and virulent subtypes. Important research gaps remain including a lack of subtyping tools for manyCryptosporidiumspecies of public and veterinary health importance, and poor understanding of the genetic determinants of host specificity ofCryptosporidiumspecies and impact of climate change on the transmission ofCryptosporidium.

Next-generation sequencing in veterinary medicine: how can the massive amount of information arising from high-throughput technologies improve diagnosis, control, and management of infectious diseases?

The development of high-throughput molecular technologies and associated bioinformatics has dramatically changed the capacities of scientists to produce, handle, and analyze large amounts of genomic, transcriptomic, and proteomic data. A clear example of this step-change is represented by the amount of DNA sequence data that can be now produced using next-generation sequencing (NGS) platforms. Similarly, recent improvements in protein and peptide separation efficiencies and highly accurate mass spectrometry have promoted the identification and quantification of proteins in a given sample. These advancements in biotechnology have increasingly been applied to the study of animal infectious diseases and are beginning to revolutionize the way that biological and evolutionary processes can be studied at the molecular level. Studies have demonstrated the value of NGS technologies for molecular characterization, ranging from metagenomic characterization of unknown pathogens or microbial communities to molecular epidemiology and evolution of viral quasispecies. Moreover, high-throughput technologies now allow detailed studies of host-pathogen interactions at the level of their genomes (genomics), transcriptomes (transcriptomics), or proteomes (proteomics). Ultimately, the interaction between pathogen and host biological networks can be questioned by analytically integrating these levels (integrative OMICS and systems biology). The application of high-throughput biotechnology platforms in these fields and their typical low-cost per information content has revolutionized the resolution with which these processes can now be studied. The aim of this chapter is to provide a current and prospective view on the opportunities and challenges associated with the application of massive parallel sequencing technologies to veterinary medicine, with particular focus on applications that have a potential impact on disease control and management.

Population structure of natural and propagated isolates of Cryptosporidium parvum, C. hominis and C. meleagridis

The three protozoan species Cryptosporidium parvum, C. meleagridis and C. hominis (phylum Apicomplexa) are enteric pathogens of humans. The former two species are zoonotic and the latter is thought to infect only humans. To better characterize the structure and transmission of natural and laboratory-propagated isolates, we analyzed a collection of archived human and animal isolates of these three species by deep-sequencing polymerase chain reaction products amplified from a polymorphic sequence on chromosome 1. Thousands of screened 200-nucleotide sequences were analyzed to compare the diversity among samples, to assess the impact of laboratory propagation on population complexity and to identify taxonomically mixed isolates. Contrary to our expectation, repeated propagation in animals did not reduce intra-isolate diversity nor was diversity associated with host species. Significantly, in most samples, sequences characteristic of a different species were identified. The presence of C. hominis alleles in C. parvum and C. meleagridis isolates confirms earlier reports of mixed isolates and raises the possibility that the host range of C. hominis is broader than typically assumed. In a genetically divergent isolate of C. parvum, a majority of sequences was found to be recombinant, suggesting that this genotype originated from a C. parvum × C. hominis recombination event.

Validation of fragment analysis by capillary electrophoresis to resolve mixed infections by Cryptosporidium parvum subpopulations

The potential of capillary electrophoresis (CE)-based DNA fragment analysis to identify mixed infections by Cryptosporidium parvum subpopulations was validated using high-resolution slab-gel electrophoresis. A selection of genomic DNA samples from C. parvum isolates with CE electropherogram profiles indicative of two concurrent alleles at one or more of six mini and microsatellite loci (MSB, MS5, ML1, ML2, TP14, 5B12) were analysed. These loci were PCR-amplified and products separated on precast Spreadex EL600 slab gels. ML1 PCR products differing by as little as 3 bp in length were visible after Spreadex gel electrophoresis and fragments were clearly separated for all but the ML2 and 5B12 loci, which generated stutter bands. No stuttering was seen for the remaining markers, having three or more nucleotide motifs in the repeat region. For each sample, the two bands of interest were excised separately, DNA extracted and re-amplified by PCR. Sequencing of these PCR products revealed the expected sequences for both alleles at most samples, except for the longest ML2 and 5B12 alleles which generated indeterminate sequences. Two novel MS5 alleles were successfully sequenced after PCR re-amplification. These findings demonstrate the utility of high-resolution Spreadex gels for analysing the polymorphism of satellite markers of Cryptosporidium isolates and support the validity of CE as a reliable and sensitive tool for detecting mixed Cryptosporidium subpopulations in a single-host infection.

New Cryptosporidium parvum subtypes of IIa subfamily in dairy calves from Brazil

Bovine cryptosporidiosis is mainly caused by four distinct species: Cryptosporidium parvum, C. bovis, C. ryanae and C. andersoni. The first, C. parvum, is a major concern in livestock causing economic losses, in addition to public health impact because of its zoonotic characteristics. The present study aimed to determine the occurrence of different species and subtypes of Cryptosporidium using molecular techniques. A total of 143 fecal samples were collected from calves from three dairy farms located in the state of Rio de Janeiro, Brazil. Saturated sugar centrifugal flotation method was used for the microscopic evaluation of the samples. Among these samples, 19.6% (28) were positive by microscopy, and 82.1% (23) of these 28 samples had their diagnosis confirmed by PCR using 18S as gene target. After sequencing, three species of Cryptosporidium were found to infect calves in different age groups. In pre-weaning phase (<2 months), 10% (3/30) of the calves were infected with C. parvum, whereas 14.2% (16/113) of post-weaning calves (≥2 months) were observed to be infected with C. andersoni and 1.8% (2/113) by C. ryanae with the latter diagnosed for the first time in the state of Rio de Janeiro. Those samples identified as C. parvum were further characterized at the GP60 locus, and PCR products were cloned. Eight different subtypes (IIaA20G2R1, IIaA20G2R2, IIaA19G2R1, IIaA19G2R2, IIaA18G1R1, IIaA18G2R2, IIaA16G3R2 and IIaA14G2R2) of C. parvum were identified, all belonging to the IIa family subtype, which is considered of high zoonotic potential. The subtypes mentioned above have not yet been detected in Brazilian cattle, and four of these subtypes (IIaA20G2R2, IIaA19G2R2, IIaA18G2R2 and IIaA14G2R2) had not been diagnosed elsewhere in calves until this study.

Cryptosporidiosis and Cryptosporidium species in animals and humans: a thirty colour rainbow?

Parasites of the genus Cryptosporidium (Apicomplexa) cause cryptosporidiosis in humans and animals worldwide. The species names used for Cryptosporidium spp. are confusing for parasitologists and even more so for non-specialists. Here, 30 named species of the genus Cryptosporidium are reviewed and proposed as valid. Molecular and experimental evidence suggests that humans and cattle are the hosts for 14 and 13 out of 30 named species, respectively. Two, four and eight named species are considered of major, moderate and minor public health significance, respectively. There are at least nine named species that are shared between humans and cattle. The aim of this review is to outline available species information together with the most commonly used genetic markers enabling the identification of named Cryptosporidium spp. Currently, 28 of 30 named species can be identified using the complete or partial ssrRNA, serving as a retrospective 'barcode'. Currently, the ssrRNA satisfies the implicit assumption that the reference databases used for comparison are sufficiently complete and applicable across the whole genus. However, due to unreliable annotation in public DNA repositories, the reference nucleotide entries and alignment of named Cryptosporidium spp. has been compiled. Despite its known limitations, ssrRNA remains the optimal marker for species identification.