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

SKSR1 identified as key virulence factor in Cryptosporidium by genetic crossing

Cryptosporidium is a major cause of severe diarrhea. Although Cryptosporidium isolates exhibit significant differences in infectivity and virulence, the genetic determinants for these traits are not clear. In this study, we use classical genetics to cross two Cryptosporidium parvum isolates of different virulence and use bulk segregant analysis of whole-genome sequences from the progeny to identify quantitative trait loci (QTL) associated with Cryptosporidium infectivity and virulence. Of the 23 genes in three QTL, two have loss-of-function mutations in the low-virulence isolates, including the SKSR1 gene encoding a variant secretory protein. Deletion of the SKSR1 gene or expression of the frame-shifted sequence reduces the pathogenicity of the virulent isolate. SKSR1 is expressed in small granules and secreted into the parasite-host interface during invasion. These results demonstrate that SKSR1 is an important virulence factor in Cryptosporidium, and suggest that the extended SKSR protein family, encoded by clusters of subtelomeric genes, may contribute to pathogenesis.

Stage-specific expression and divergent functions of two insulinase-like proteases associated with host infectivity in Cryptosporidium

BACKGROUND

The determinants of differences in host infectivity among Cryptosporidium species and subtypes are poorly understood. Results from recent comparative genomic studies suggest that gains and losses of multicopy subtelomeric genes encoding insulinase-like proteases (INS-19 and INS-20 in Cryptosporidium parvum and their orthologs in closely related species) may potentially contribute to these differences.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we investigated the expression and biological function of the INS-19 and INS-20 of C. parvum. CRISPR/Cas9 was used to endogenously tag both genes with the hemagglutinin epitope. Immunofluorescence analysis revealed that INS-19 and INS-20 are expressed at different developmental stages of the pathogen. Although knockout of either had no detectable effect on the in vitro growth of C. parvum, knockout of INS-20, deletion of its multiple domains, or mutation of the active motif in the functional domain reduced the intensity of C. parvum infection in IFN-γ knockout mice. Consistent with this, mice infected with the INS-20-deleted mutant had reduced intestinal damage and parasite burden.

CONCLUSIONS/SIGNIFICANCE

These results suggest that INS-19 and INS-20 have stage-specific expression with distinct biological functions, and that the presence of the INS-20 in zoonotic C. parvum contributes to its infectivity and fitness in mice.

Multicopy subtelomeric genes underlie animal infectivity of divergent Cryptosporidium hominis subtypes

The anthroponotic Cryptosporidium hominis differs from the zoonotic C. parvum in its lack of infectivity to animals, but several divergent subtypes have recently been found in nonhuman primates and equines. Here, we sequence 17 animal C. hominis isolates and generate a new IbA12G3 genome at the chromosome level. Comparative analysis with 222 human isolates shows significant genetic divergence of the animal isolates, with genetic recombination among them. They have additional subtelomeric insulinase and MEDLE genes. In interferon-γ knockout mice, three monkey isolates show differences in infectivity and induce higher and longer oocyst shedding than a reference C. parvum isolate. Deletion of the MEDLE genes significantly reduces the growth and pathogenicity of a virulent strain in mice. Co-infection of two fluorescence-tagged C. hominis subtypes produces bicolored oocysts, supporting the conclusion that mixed subtype infections can lead to genetic recombination. These data provide insight into potential determinants of host infectivity in Cryptosporidium, and a convenient animal model for biological studies of C. hominis.

Characterization of NFDQ1 in Cryptosporidium parvum

BACKGROUND

Cryptosporidium spp. are important zoonotic parasites that can cause moderate to severe diarrhea in humans and animals. Among the three Cryptosporidium species infecting the intestines of calves, Cryptosporidium parvum has a broad host range and causes severe diarrhea in calves, while Cryptosporidium bovis and Cryptosporidium ryanae mainly infect calves without obvious clinical symptoms. Comparative genomic analysis revealed differences in the copy number of genes encoding the nonfinancial disclosure quality (NFDQ) secretory protein family among the three species, suggesting that this protein family may be associated with the host range or pathogenicity of Cryptosporidium spp. To understand the function of cgd8_10 encoded NFDQ1, tagged and knockout strains were constructed and characterized in this study.

METHODS

To determine the localization of NFDQ1, we used clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) technology to tag the C-terminus of NFDQ1 with three hemagglutinin epitopes (3 × HA). The tagged strain was constructed, and the genomic insertion was confirmed by polymerase chain reaction (PCR). Immunofluorescence assays were performed to observe the localization of NFDQ1 both in extracellular sporozoites and at various intracellular developmental stages. Immunoelectron microscopy was used to study the ultrastructural localization of NFDQ1. Then, the ΔNFDQ1 strain was generated by CRISPR/Cas9 and the in vitro growth assay on HCT-8 cells was used to analyze of phenotypic changes after knockout NFDQ1 in parasites.

RESULTS

The NFDQ1 tagging and knockout stains were successfully constructed by CRISPR/Cas9 technology and the insertions of transgenic strains were validated by PCR. The expression of NFDQ1 was validated in parasite by western blot. Immunofluorescence and immune-electron microscopy assay showed that NFDQ1 expressed in both asexual and sexual stages of C. parvum, where it was localized to the cytoplasm of the parasite. Upon ablation of NFDQ1, the ΔNFDQ1 strain showed an apparent growth retardation during sexual replication in vitro.

CONCLUSIONS

NFDQ1 is a cytoplasmic protein without specific localization to secretory organelles, and it may participate in C. parvum growth during sexual reproduction. Future study should determine the role of NFDQ1 following C. parvum infection in vivo.

Exploring the impact of digestive physicochemical parameters of adults and infants on the pathophysiology of Cryptosporidium parvum using the dynamic TIM-1 gastrointestinal model

BACKGROUND

Human cryptosporidiosis is distributed worldwide, and it is recognised as a leading cause of acute diarrhoea and death in infants in low- and middle-income countries. Besides immune status, the higher incidence and severity of this gastrointestinal disease in young children could also be attributed to the digestive environment. For instance, human gastrointestinal physiology undergoes significant changes with age, however the role this variability plays in Cryptosporidium parvum pathogenesis is not known. In this study, we analysed for the first time the impact of digestive physicochemical parameters on C. parvum infection in a human and age-dependent context using a dynamic in vitro gastrointestinal model.

RESULTS

Our results showed that the parasite excystation, releasing sporozoites from oocysts, occurs in the duodenum compartment after one hour of digestion in both child (from 6 months to 2 years) and adult experimental conditions. In the child small intestine, slightly less sporozoites were released from excystation compared to adult, however they exhibited a higher luciferase activity, suggesting a better physiological state. Sporozoites collected from the child jejunum compartment also showed a higher ability to invade human intestinal epithelial cells compared to the adult condition. Global analysis of the parasite transcriptome through RNA-sequencing demonstrated a more pronounced modulation in ileal effluents compared to gastric ones, albeit showing less susceptibility to age-related digestive condition. Further analysis of gene expression and enriched pathways showed that oocysts are highly active in protein synthesis in the stomach compartment, whereas sporozoites released in the ileum showed downregulation of glycolysis as well as strong modulation of genes potentially related to gliding motility and secreted effectors.

CONCLUSIONS

Digestion in a sophisticated in vitro gastrointestinal model revealed that invasive sporozoite stages are released in the small intestine, and are highly abundant and active in the ileum compartment, supporting reported C. parvum tissue tropism. Our comparative analysis suggests that physicochemical parameters encountered in the child digestive environment can influence the amount, physiological state and possibly invasiveness of sporozoites released in the small intestine, thus potentially contributing to the higher susceptibility of young individuals to cryptosporidiosis.

Novel strategy to quantify the viability of oocysts of Cryptosporidium parvum and C. hominis, a risk factor of the waterborne protozoan pathogens of public health concern

While waters might be contaminated by oocysts from >40 Cryptosporidium species, only viable oocysts of C. parvum and C. hominis truly pose the main health risk to the immunocompetent population. Oocyst viability is also an important but often neglected risk factor in monitoring waterborne parasites. However, commonly used methods in water monitoring and surveys cannot distinguish species (microscopic observation) or oocyst viability (PCR), as dead oocysts in water could retain gross structure and DNA content for weeks to months. Here, we report new TaqMan qRT-PCR/qPCR assays for quantitative detection of viable C. parvum and C. hominis oocysts. By targeting a hypothetical protein-encoding gene cgd6_3920 that is highly expressed in oocysts and variable between species, the qRT-PCR/qPCR assays achieve excellent analytical specificity and sensitivity (limit of quantification [LOQ] = 0.25 and 1.0 oocyst/reaction). Using calibration curves, the number and ratio of viable oocysts in specimens could be calculated. Additionally, we also establish a TaqMan-18S qPCR for cost-effective screening of pan-Cryptosporidium-positive specimens (LOQ = 0.1 oocyst/reaction). The assay feasibility is validated using field water (N = 43) and soil (79) specimens from 17 locations in Changchun, China, which detects four Cryptosporidium species from seven locations, including three gp60-subtypes (i.e., IIdA19G1, IIdA17G1 and IIdA24G2) of C. parvum oocysts showing varied viability ratios. These new TaqMan q(RT)-PCR assays supplement current methods in the survey of waters and other samples (e.g., surfaces, foods and beverages), and are applicable to assessing the efficiency of oocyst deactivation protocols.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

WHOLE GENOME TARGETED ENRICHMENT AND SEQUENCING OF HUMAN-INFECTING CRYPTOSPORIDIUM spp

Cryptosporidium spp. are protozoan parasites that cause severe illness in vulnerable human populations. Obtaining pure Cryptosporidium DNA from clinical and environmental samples is challenging because the oocysts shed in contaminated feces are limited in quantity, difficult to purify efficiently, may derive from multiple species, and yield limited DNA (<40 fg/oocyst). Here, we develop and validate a set of 100,000 RNA baits (CryptoCap_100k) based on six human-infecting Cryptosporidium spp. (C. cuniculus, C. hominis, C. meleagridis, C. parvum, C. tyzzeri, and C. viatorum) to enrich Cryptosporidium spp. DNA from a wide array of samples. We demonstrate that CryptoCap_100k increases the percentage of reads mapping to target Cryptosporidium references in a wide variety of scenarios, increasing the depth and breadth of genome coverage, facilitating increased accuracy of detecting and analyzing species within a given sample, while simultaneously decreasing costs, thereby opening new opportunities to understand the complex biology of these important pathogens.

Sequence introgression from exogenous lineages underlies genomic and biological differences among Cryptosporidium parvum IOWA lines

The IOWA strain of Cryptosporidium parvum is widely used in studies of the biology and detection of the waterborne pathogens Cryptosporidium spp. While several lines of the strain have been sequenced, IOWA-II, the only reference of the original subtype (IIaA15G2R1), exhibits significant assembly errors. Here we generated a fully assembled genome of IOWA-CDC of this subtype using PacBio and Illumina technologies. In comparative analyses of seven IOWA lines maintained in different laboratories (including two sequenced in this study) and 56 field isolates, IOWA lines (IIaA17G2R1) with less virulence had mixed genomes closely related to IOWA-CDC but with multiple sequence introgressions from IOWA-II and unknown lineages. In addition, the IOWA-IIaA17G2R1 lines showed unique nucleotide substitutions and loss of a gene associated with host infectivity, which were not observed in other isolates analyzed. These genomic differences among IOWA lines could be the genetic determinants of phenotypic traits in C. parvum. These data provide a new reference for comparative genomic analyses of Cryptosporidium spp. and rich targets for the development of advanced source tracking tools.

Cultivation, cryopreservation, and transcriptomic studies of host-adapted Cryptosporidium parvum and Cryptosporidium hominis using enteroids

Cryptosporidium hominis and Cryptosporidium parvum are major causes of severe diarrhea. Comparative studies of them are hampered by the lack of effective cultivation and cryopreservation methods, especially for C. hominis. Here, we describe adapted murine enteroids for the cultivation and complete development of host-adapted C. parvum and C. hominis subtypes, producing oocysts infectious to mice. Using the system, we developed a cryopreservation method for Cryptosporidium isolates. In comparative RNA-seq analyses of C. hominis cultures, the enteroid system generated significantly more host and pathogen responses than the conventional HCT-8 cell system. In particular, the infection was shown to upregulate PI3K-Akt, Ras, TNF, NF-κB, IL-17, MAPK, and innate immunity signaling pathways and downregulate host cell metabolism, and had significantly higher expression of parasite genes involved in oocyst formation. Therefore, the enteroid system provides a valuable tool for comparative studies of the biology of divergent Cryptosporidium species and isolates.

WHOLE GENOME TARGETED ENRICHMENT AND SEQUENCING OF HUMAN-INFECTING CRYPTOSPORIDIUM spp

Cryptosporidium spp. are protozoan parasites that cause severe illness in vulnerable human populations. Obtaining pure Cryptosporidium DNA from clinical and environmental samples is challenging because the oocysts shed in contaminated feces are limited in quantity, difficult to purify efficiently, may derive from multiple species, and yield limited DNA (<40 fg/oocyst). Here, we develop and validate a set of 100,000 RNA baits (CryptoCap_100k) based on six human-infecting Cryptosporidium spp. ( C. cuniculus , C. hominis , C. meleagridis , C. parvum , C. tyzzeri , and C. viatorum ) to enrich Cryptosporidium spp. DNA from a wide array of samples. We demonstrate that CryptoCap_100k increases the percentage of reads mapping to target Cryptosporidium references in a wide variety of scenarios, increasing the depth and breadth of genome coverage, facilitating increased accuracy of detecting and analyzing species within a given sample, while simultaneously decreasing costs, thereby opening new opportunities to understand the complex biology of these important pathogens.

Identification of host protein ENO1 (alpha-enolase) interacting with Cryptosporidium parvum sporozoite surface protein, Cpgp40

BACKGROUND

Cryptosporidium parvum is an apicomplexan zoonotic parasite causing the diarrheal illness cryptosporidiosis in humans and animals. To invade the host intestinal epithelial cells, parasitic proteins expressed on the surface of sporozoites interact with host cells to facilitate the formation of parasitophorous vacuole for the parasite to reside and develop. The gp40 of C. parvum, named Cpgp40 and located on the surface of sporozoites, was proven to participate in the process of host cell invasion.

METHODS

We utilized the purified Cpgp40 as a bait to obtain host cell proteins interacting with Cpgp40 through the glutathione S-transferase (GST) pull-down method. In vitro analysis, through bimolecular fluorescence complementation assay (BiFC) and coimmunoprecipitation (Co-IP), confirmed the solid interaction between Cpgp40 and ENO1. In addition, by using protein mutation and parasite infection rate analysis, it was demonstrated that ENO1 plays an important role in the C. parvum invasion of HCT-8 cells.

RESULTS

To illustrate the functional activity of Cpgp40 interacting with host cells, we identified the alpha-enolase protein (ENO1) from HCT-8 cells, which showed direct interaction with Cpgp40. The mRNA level of ENO1 gene was significantly decreased at 3 and 24 h after C. parvum infection. Antibodies and siRNA specific to ENO1 showed the ability to neutralize C. parvum infection in vitro, which indicated the participation of ENO1 during the parasite invasion of HCT-8 cells. In addition, we further demonstrated that ENO1 protein was involved in the regulation of cytoplasmic matrix of HCT-8 cells during C. parvum invasion. Functional study of the protein mutation illustrated that ENO1 was also required for the endogenous development of C. parvum.

CONCLUSIONS

In this study, we utilized the purified Cpgp40 as a bait to obtain host cell proteins ENO1 interacting with Cpgp40. Functional studies illustrated that the host cell protein ENO1 was involved in the regulation of tight junction and adherent junction proteins during C. parvum invasion and was required for endogenous development of C. parvum.

Genetic crosses within and between species of Cryptosporidium

Parasites and their hosts are engaged in reciprocal coevolution that balances competing mechanisms of virulence, resistance, and evasion. This often leads to host specificity, but genomic reassortment between different strains can enable parasites to jump host barriers and conquer new niches. In the apicomplexan parasite Cryptosporidium, genetic exchange has been hypothesized to play a prominent role in adaptation to humans. The sexual lifecycle of the parasite provides a potential mechanism for such exchange; however, the boundaries of Cryptosporidium sex are currently undefined. To explore this experimentally, we established a model for genetic crosses. Drug resistance was engineered using a mutated phenylalanyl tRNA synthetase gene and marking strains with this and the previously used Neo transgene enabled selection of recombinant progeny. This is highly efficient, and genomic recombination is evident and can be continuously monitored in real time by drug resistance, flow cytometry, and PCR mapping. Using this approach, multiple loci can now be modified with ease. We demonstrate that essential genes can be ablated by crossing a Cre recombinase driver strain with floxed strains. We further find that genetic crosses are also feasible between species. Crossing Cryptosporidium parvum, a parasite of cattle and humans, and Cryptosporidium tyzzeri a mouse parasite resulted in progeny with a recombinant genome derived from both species that continues to vigorously replicate sexually. These experiments have important fundamental and translational implications for the evolution of Cryptosporidium and open the door to reverse- and forward-genetic analysis of parasite biology and host specificity.

Comparison of in vitro growth characteristics of Cryptosporidium hominis (IdA15G1) and Cryptosporidium parvum (Iowa-IIaA17G2R1 and IIaA18G3R1)

Cryptosporidium is a major cause of diarrhoeal disease and mortality in young children in resource-poor countries, for which no vaccines or adequate therapeutic options are available. Infection in humans is primarily caused by two species: C. hominis and C. parvum. Despite C. hominis being the dominant species infecting humans in most countries, very little is known about its growth characteristics and life cycle in vitro, given that the majority of our knowledge of the in vitro development of Cryptosporidium has been based on C. parvum. In the present study, the growth and development of two C. parvum isolates (subtypes Iowa-IIaA17G2R1 and IIaA18G3R1) and one C. hominis isolate (subtype IdA15G1) in HCT-8 cells were examined and compared at 24 h and 48 h using morphological data acquired with scanning electron microscopy. Our data indicated no significant differences in the proportion of meronts or merozoites between species or subtypes at either time-point. Sexual development was observed at the 48-h time-point across both species through observations of both microgamonts and macrogamonts, with a higher frequency of macrogamont observations in C. hominis (IdA15G1) cultures at 48-h post-infection compared to both C. parvum subtypes. This corresponded to differences in the proportion of trophozoites observed at the same time point. No differences in proportion of microgamonts were observed between the three subtypes, which were rarely observed across all cultures. In summary, our data indicate that asexual development of C. hominis is similar to that of C. parvum, while sexual development is accelerated in C. hominis. This study provides new insights into differences in the in vitro growth characteristics of C. hominis when compared to C. parvum, which will facilitate our understanding of the sexual development of both species.

Investigating Cryptosporidium spp. Using Genomic, Proteomic and Transcriptomic Techniques: Current Progress and Future Directions

Cryptosporidiosis is a widespread disease caused by the parasitic protozoan Cryptosporidium spp., which infects various vertebrate species, including humans. Once unknown as a gastroenteritis-causing agent, Cryptosporidium spp. is now recognized as a pathogen causing life-threatening disease, especially in immunocompromised individuals such as AIDS patients. Advances in diagnostic methods and increased awareness have led to a significant shift in the perception of Cryptosporidium spp. as a pathogen. Currently, genomic and proteomic studies play a main role in understanding the molecular biology of this complex-life-cycle parasite. Genomics has enabled the identification of numerous genes involved in the parasite's development and interaction with hosts. Proteomics has allowed for the identification of protein interactions, their function, structure, and cellular activity. The combination of these two approaches has significantly contributed to the development of new diagnostic tools, vaccines, and drugs for cryptosporidiosis. This review presents an overview of the significant achievements in Cryptosporidium research by utilizing genomics, proteomics, and transcriptomics approaches.

Genetic crosses within and between species of Cryptosporidium

Parasites and their hosts are engaged in rapid coevolution that balances competing mechanisms of virulence, resistance, and evasion. This often leads to host specificity, but genomic reassortment between different strains can enable parasites to jump host barriers and conquer new niches. In the apicomplexan parasite Cryptosporidium genetic exchange has been hypothesized to play a prominent role in adaptation to humans. The sexual lifecycle of the parasite provides a potential mechanism for such exchange; however, the boundaries of Cryptosporidium sex are currently undefined. To explore this experimentally, we established a model for genetic crosses. Drug resistance was engineered using a mutated phenylalanyl tRNA synthetase gene and marking strains with this and the previously used Neo transgene enabled selection of recombinant progeny. This is highly efficient, and genomic recombination is evident and can be continuously monitored in real time by drug resistance, flow cytometry, and PCR mapping. Using this approach multiple loci can now be modified with ease. We demonstrate that essential genes can be ablated by crossing a Cre recombinase driver strain with floxed strains. We further find that genetic crosses are also feasible between species. Crossing C. parvum, a parasite of cattle and humans, and C. tyzzeri a mouse parasite resulted in progeny with a recombinant genome derived from both species that continues to vigorously replicate sexually. These experiments have important fundamental and translational implications for the evolution of Cryptosporidium and open the door to reverse- and forward- genetic analysis of parasite biology and host specificity.

High subtelomeric GC content in the genome of a zoonotic Cryptosporidium species

Cryptosporidium canis is a zoonotic species causing cryptosporidiosis in humans in addition to its natural hosts dogs and other fur animals. To understand the genetic basis for host adaptation, we sequenced the genomes of C. canis from dogs, minks, and foxes and conducted a comparative genomics analysis. While the genomes of C. canis have similar gene contents and organisations, they (~41.0 %) and C. felis (39.6 %) have GC content much higher than other Cryptosporidium spp. (24.3-32.9 %) sequenced to date. The high GC content is mostly restricted to subtelomeric regions of the eight chromosomes. Most of these GC-balanced genes encode Cryptosporidium-specific proteins that have intrinsically disordered regions and are involved in host-parasite interactions. Natural selection appears to play a more important role in the evolution of codon usage in GC-balanced C. canis, and most of the GC-balanced genes have undergone positive selection. While the identity in whole genome sequences between the mink- and dog-derived isolates is 99.9 % (9365 SNVs), it is only 96.0 % (362 894 SNVs) between them and the fox-derived isolate. In agreement with this, the fox-derived isolate possesses more subtelomeric genes encoding invasion-related protein families. Therefore, the change in subtelomeric GC content appears to be responsible for the more GC-balanced C. canis genomes, and the fox-derived isolate could represent a new Cryptosporidium species.

Cryptosporidium equi n. sp. (Apicomplexa: Cryptosporidiidae): biological and genetic characterisations

The horse genotype is one of three common Cryptosporidium spp. in equine animals and has been identified in some human cases. The species status of Cryptosporidium horse genotype remains unclear due to the lack of extensive morphological, biological, and genetic data. In the present study, we have conducted biological and whole genome sequence analyses of an isolate of the genotype from hedgehogs and proposed to name it Cryptosporidium equi n. sp. to reflect its common occurrence in equine animals. Oocysts of C. equi measured 5.12 ± 0.36 μm × 4.46 ± 0.21 μm with a shape index of 1.15 ± 0.08 (n = 50). Cryptosporidium equi was infectious to 3-week-old four-toed hedgehogs (Atelerix albiventris) and mice, with a prepatent period of 2-9 days and a patent period of 30-40 days in hedgehogs. It was not infectious to rats and rabbits. Phylogenetic analyses of small subunit rRNA, 70 kDa heat shock protein, actin, 60 kDa glycoprotein and 100 other orthologous genes revealed that C. equi is genetically distinct from other known Cryptosporidium species and genotypes. The sequence identity between C. equi and Cryptosporidium parvum genomes is 97.9%. Compared with C. parvum, C. equi has lost two MEDLE genes and one insulinase-like protease gene and gained one SKSR gene. In addition, 60 genes have highly divergent sequences (sequence differences ≥ 5.0%), including those encoding mucin-like glycoproteins, insulinase-like peptidases, and MEDLE and SKSR proteins. The genetic uniqueness of C. equi supports its increasing host range and the naming of it as a valid Cryptosporidium species. This is the first known use of whole genome sequence data in delineating new Cryptosporidium species.

Multiple introductions and recombination events underlie the emergence of a hyper-transmissible Cryptosporidium hominis subtype in the USA

The parasite Cryptosporidium hominis is a leading cause of the diarrheal disease cryptosporidiosis, whose incidence in the United States has increased since 2005. Here, we show that the newly emerged and hyper-transmissible subtype IfA12G1R5 is now dominant in the United States. In a comparative analysis of 127 newly sequenced and 95 published C. hominis genomes, IfA12G1R5 isolates from the United States place into three of the 14 clusters (Pop6, Pop13, and Pop14), indicating that this subtype has multiple ancestral origins. Pop6 (IfA12G1R5a) has an East Africa origin and has recombined with autochthonous subtypes after its arrival. Pop13 (IfA12G1R5b) is imported from Europe, where it has recombined with the prevalent local subtype, whereas Pop14 (IfA12G1R5c) is a progeny of secondary recombination between Pop6 and Pop13. Selective sweeps in invasion-associated genes have accompanied the emergence of the dominant Pop14. These observations offer insights into the emergence and evolution of hyper-transmissible pathogens.

High infectivity and unique genomic sequence characteristics of Cryptosporidium parvum in China

Zoonotic Cryptosporidium parvum infections are mainly caused by IIa and IId subtypes. As most biological characterizations have been performed on IIa subtypes, the biological and genetic characteristics of IId subtypes in China are not clear. We evaluated the infection and genetic characteristics of IId isolates in interferon-γ-knockout mice using qPCR to quantify oocyst shedding, histological examination to monitor pathological changes and comparative genomic analyses to identify infectivity and virulence-associated differences. Compared with the reference IIa isolate, mice infected with the IId isolates had significantly higher and longer oocyst shedding and lower body weight gain. In addition, the four IId isolates examined differed significantly in infectivity (as indicated by the median infective dose), oocyst shedding duration, and pathogenicity. Comparative genomic analysis indicated that the IId isolates had three more subtelomeric genes than the reference IIa isolate and 5385–5548 nucleotide substitutions, with the hypervariable genes mostly in two blocks on chromosome 1. In contrast, the four IId isolates differed from each other by 77–1,452 nucleotides, with virulence-associated sequence differences mainly in nine genes within a 28-kb block on chromosome 6. These data indicate the newly emerged C. parvum IId subtypes in China have high animal infectivity and unique genomic characteristics.

Cryptosporidiosis is the most important waterborne disease in industrialized nations and a primary cause of severe diarrhea in children in low- and middle-income countries. While the IIa subtype family of Cryptosporidium parvum is responsible for most zoonotic cryptosporidiosis, its IId subtype family has emerged in China in recent years. To understand the biological differences between the two major zoonotic subtype families, we have compared the infection patterns, virulence, and genetic characteristics of IIa and IId isolates using a newly established mouse model and whole genome sequencing. We have shown that IId isolates induce significantly higher infection intensity, longer infection duration, and more severe pathogenicity than the reference IIa isolate. They also have three more invasion-associated genes and substantial nucleotide sequence differences. In contrast, the four IId isolates with different virulence differ from each other mainly in sequences of nine genes within a small area on chromosome 6. We conclude that C. parvum isolates in China have high infectivity and unique genomic characteristics, and the productive infection model developed in the study should be useful in evaluations of potential therapeutics and studies of pathogenesis of C. parvum.

Characterization of Dense Granule Metalloproteinase INS-16 in Cryptosporidium parvum

The protozoan pathogen Cryptosporidium parvum infects intestinal epithelial cells and causes diarrhea in humans and young animals. Among the more than 20 genes encoding insulinase-like metalloproteinases (INS), two are paralogs with high sequence identity. In this study, one of them, INS-16 encoded by the cgd3_4270 gene, was expressed and characterized in a comparative study of its sibling, INS-15 encoded by the cgd3_4260 gene. A full-length INS-16 protein and its active domain I were expressed in Escherichia coli, and antibodies against the domain I and an INS-16-specific peptide were produced in rabbits. In the analysis of the crude extract of oocysts, a ~60 kDa fragment of INS-16 rather than the full protein was recognized by polyclonal antibodies against the specific peptide, indicating that INS-16 undergoes proteolytic cleavage before maturation. The expression of the ins-16 gene peaked at the invasion phase of in vitro C. parvum culture, with the documented expression of the protein in both sporozoites and merozoites. Localization studies with antibodies showed significant differences in the distribution of the native INS-15 and INS-16 proteins in sporozoites and merozoites. INS-16 was identified as a dense granule protein in sporozoites and macrogamonts but was mostly expressed at the apical end of merozoites. We screened 48 candidate INS-16 inhibitors from the molecular docking of INS-16. Among them, two inhibited the growth of C. parvum in vitro (EC₅₀ = 1.058 µM and 2.089 µM). The results of this study suggest that INS-16 may have important roles in the development of C. parvum and could be a valid target for the development of effective treatments.

Marine gregarine genomes reveal the breadth of apicomplexan diversity with a partially conserved glideosome machinery

Our current view of the evolutionary history, coding and adaptive capacities of Apicomplexa, protozoan parasites of a wide range of metazoan, is currently strongly biased toward species infecting humans, as data on early diverging apicomplexan lineages infecting invertebrates is extremely limited. Here, we characterized the genome of the marine eugregarine Porospora gigantea, intestinal parasite of Lobsters, remarkable for the macroscopic size of its vegetative feeding forms (trophozoites) and its gliding speed, the fastest so far recorded for Apicomplexa. Two highly syntenic genomes named A and B were assembled. Similar in size (~ 9 Mb) and coding capacity (~ 5300 genes), A and B genomes are 10.8% divergent at the nucleotide level, corresponding to 16-38 My in divergent time. Orthogroup analysis across 25 (proto)Apicomplexa species, including Gregarina niphandrodes, showed that A and B are highly divergent from all other known apicomplexan species, revealing an unexpected breadth of diversity. Phylogenetically these two species branch sisters to Cephaloidophoroidea, and thus expand the known crustacean gregarine superfamily. The genomes were mined for genes encoding proteins necessary for gliding, a key feature of apicomplexans parasites, currently studied through the molecular model called glideosome. Sequence analysis shows that actin-related proteins and regulatory factors are strongly conserved within apicomplexans. In contrast, the predicted protein sequences of core glideosome proteins and adhesion proteins are highly variable among apicomplexan lineages, especially in gregarines. These results confirm the importance of studying gregarines to widen our biological and evolutionary view of apicomplexan species diversity, and to deepen our understanding of the molecular bases of key functions such as gliding, well known to allow access to the intracellular parasitic lifestyle in Apicomplexa.

Sympatric Recombination in Zoonotic Cryptosporidium Leads to Emergence of Populations with Modified Host Preference

Genetic recombination plays a critical role in the emergence of pathogens with phenotypes such as drug resistance, virulence, and host adaptation. Here, we tested the hypothesis that recombination between sympatric ancestral populations leads to the emergence of divergent variants of the zoonotic parasite Cryptosporidium parvum with modified host ranges. Comparative genomic analyses of 101 isolates have identified seven subpopulations isolated by distance. They appear to be descendants of two ancestral populations, IIa in northwestern Europe and IId from southwestern Asia. Sympatric recombination in areas with both ancestral subtypes and subsequent selective sweeps have led to the emergence of new subpopulations with mosaic genomes and modified host preference. Subtelomeric genes could be involved in the adaptive selection of subpopulations, while copy number variations of genes encoding invasion-associated proteins are potentially associated with modified host ranges. These observations reveal ancestral origins of zoonotic C. parvum and suggest that pathogen import through modern animal farming might promote the emergence of divergent subpopulations of C. parvum with modified host preference.

An update on Cryptosporidium biology and therapeutic avenues

Cryptosporidium species has been identified as an important pediatric diarrheal pathogen in resource-limited countries, particularly in very young children (0–24 months). However, the only available drug (nitazoxanide) has limited efficacy and can only be prescribed in a medical setting to children older than one year. Many drug development projects have started to investigate new therapeutic avenues. Cryptosporidium’s unique biology is challenging for the traditional drug discovery pipeline and requires novel drug screening approaches. Notably, in recent years, new methods of oocyst generation, in vitro processing, and continuous three-dimensional cultivation capacities have been developed. This has enabled more physiologically pertinent research assays for inhibitor discovery. In a short time, many great strides have been made in the development of anti-Cryptosporidium drugs. These are expected to eventually turn into clinical candidates for cryptosporidiosis treatment in the future. This review describes the latest development in Cryptosporidium biology, genomics, transcriptomics of the parasite, assay development, and new drug discovery.

Global population genomics of two subspecies of Cryptosporidium hominis during 500 years of evolution

Cryptosporidiosis is a major global health problem and a primary cause of diarrhea, particularly in young children in low- and middle-income countries (LMICs). The zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis cause most human infections. Here, we present a comprehensive whole-genome study of C. hominis, comprising 114 isolates from 16 countries within five continents. We detect two lineages with distinct biology and demography, which diverged circa 500 years ago. We consider these lineages two subspecies and propose the names C. hominis hominis and C. hominis aquapotentis (gp60 subtype IbA10G2). In our study, C. h. hominis is almost exclusively represented by isolates from LMICs in Africa and Asia and appears to have undergone recent population contraction. In contrast, C. h. aquapotentis was found in high-income countries, mainly in Europe, North America, and Oceania, and appears to be expanding. Notably, C. h. aquapotentis is associated with high rates of direct human-to-human transmission, which may explain its success in countries with well-developed environmental sanitation infrastructure. Intriguingly, we detected genomic regions of introgression following secondary contact between the subspecies. This resulted in high diversity and divergence in genomic islands of putative virulence genes, including muc5 (CHUDEA2_430) and a hypothetical protein (CHUDEA6_5270). This diversity is maintained by balancing selection, suggesting a co-evolutionary arms race with the host. Finally, we find that recent gene flow from C. h. aquapotentis to C. h. hominis, likely associated with increased human migration, maybe driving the evolution of more virulent C. hominis variants.

A productive immunocompetent mouse model of cryptosporidiosis with long oocyst shedding duration for immunological studies

OBJECTIVES

Studies on the pathogenesis and immune responses of Cryptosporidium infection and development of drugs and vaccines use mostly immunocompromised mouse models. In this study, we establish an immunocompetent mouse model of cryptosporidiosis with high intensity and long duration of infection.

METHODS

We have obtained a Cryptosporidium tyzzeri isolate from laboratory mice, and infect adult C57BL/6J mice experimentally with the isolate for determinations of infectivity, infection patterns, pathological changes, and transcriptomic responses.

RESULTS

The isolate has an ID₅₀ of 5.2 oocysts, with oocyst shedding lasting at high levels for >2 months. The oocyst shedding is boosted by immunosuppression of animals and suppressed by paromomycin treatment. The isolate induces strong inflammatory and acquired immune responses, but down-regulates the expression of α-defensins in epithelium. Comparative genomics analysis has revealed significant sequence differences from other isolates in subtelomeric genes. The down-regulation of the expression of α-defensins may be responsible for the high-intensity and long-lasting infection in this animal model.

CONCLUSIONS

The immunocompetent mouse model of cryptosporidiosis developed has the advantages of high oocyst shedding intensity and long oocyst shedding duration. It provides an effective mechanism for the propagation of Cryptosporidium, evaluations of potential therapeutics, and studies of pathogen biology and immune responses.

Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveal expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design, and interpretation. We have generated a new C. parvum IOWA genome assembly supported by Pacific Biosciences (PacBio) and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species: C. parvum, Cryptosporidium hominis, and Cryptosporidium tyzzeri We made 1926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns, and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis, and C. tyzzeri revealed that most "missing" orthologs are found, suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation, and single-nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free, and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

Cryptosporidium felis differs from other Cryptosporidium spp. in codon usage

Cryptosporidium spp. are important enteric pathogens in a wide range of vertebrates including humans. Previous comparative analysis revealed conservation in genome composition, gene content, and gene organization among Cryptosporidium spp., with a progressive reductive evolution in metabolic pathways and invasion-related proteins. In this study, we sequenced the genome of zoonotic pathogen Cryptosporidium felis and conducted a comparative genomic analysis. While most intestinal Cryptosporidium species have similar genomic characteristics and almost complete genome synteny, fewer protein-coding genes and some sequence inversions and translocations were found in the C. felis genome. The C. felis genome exhibits much higher GC content (39.6 %) than other Cryptosporidium species (24.3–32.9 %), especially at the third codon position (GC3) of protein-coding genes. Thus, C. felis has a different codon usage, which increases the use of less energy costly amino acids (Gly and Ala) encoded by GC-rich codons. While the tRNA usage is conserved among Cryptosporidium species, consistent with its higher GC content, C. felis uses a unique tRNA for GTG for valine instead of GTA in other Cryptosporidium species. Both mutational pressures and natural selection are associated with the evolution of the codon usage in Cryptosporidium spp., while natural selection seems to drive the codon usage in C. felis. Other unique features of the C. felis genome include the loss of the entire traditional and alternative electron transport systems and several invasion-related proteins. Thus, the preference for the use of some less energy costly amino acids in C. felis may lead to a more harmonious parasite–host interaction, and the strengthened host-adaptation is reflected by the further reductive evolution of metabolism and host invasion-related proteins.

The enteric pathogen Cryptosporidium parvum exports proteins into the cytosol of the infected host cell

The parasite Cryptosporidium is responsible for diarrheal disease in young children causing death, malnutrition, and growth delay. Cryptosporidium invades enterocytes where it develops in a unique intracellular niche. Infected cells exhibit profound changes in morphology, physiology, and transcriptional activity. How the parasite effects these changes is poorly understood. We explored the localization of highly polymorphic proteins and found members of the Cryptosporidium parvum MEDLE protein family to be translocated into the cytosol of infected cells. All intracellular life stages engage in this export, which occurs after completion of invasion. Mutational studies defined an N-terminal host-targeting motif and demonstrated proteolytic processing at a specific leucine residue. Direct expression of MEDLE2 in mammalian cells triggered an ER stress response, which was also observed during infection. Taken together, our studies reveal the presence of a Cryptosporidium secretion system capable of delivering parasite proteins into the infected enterocyte.

An Update on Zoonotic Cryptosporidium Species and Genotypes in Humans

The enteric parasite, Cryptosporidium is a major cause of diarrhoeal illness in humans and animals worldwide. No effective therapeutics or vaccines are available and therefore control is dependent on understanding transmission dynamics. The development of molecular detection and typing tools has resulted in the identification of a large number of cryptic species and genotypes and facilitated our understanding of their potential for zoonotic transmission. Of the 44 recognised Cryptosporidium species and >120 genotypes, 19 species, and four genotypes have been reported in humans with C. hominis, C. parvum, C. meleagridis, C. canis and C. felis being the most prevalent. The development of typing tools that are still lacking some zoonotic species and genotypes and more extensive molecular epidemiological studies in countries where the potential for transmission is highest are required to further our understanding of this important zoonotic pathogen. Similarly, whole-genome sequencing (WGS) and amplicon next-generation sequencing (NGS) are important for more accurately tracking transmission and understanding the mechanisms behind host specificity.

Cryptosporidium hominis Phylogenomic Analysis Reveals Separate Lineages With Continental Segregation

Cryptosporidium is a leading cause of waterborne outbreaks globally, and Cryptosporidium hominis and C. parvum are the principal cause of human cryptosporidiosis on the planet. Thanks to the advances in Next-Generation Sequencing (NGS) sequencing and bioinformatic software development, more than 100 genomes have been generated in the last decade using a metagenomic-like strategy. This procedure involves the parasite oocyst enrichment from stool samples of infected individuals, NGS sequencing, metagenomic assembly, parasite genome computational filtering, and comparative genomic analysis. Following this approach, genomes of infected individuals of all continents have been generated, although with striking different quality results. In this study, we performed a thorough comparison, in terms of assembly quality and purity, of 100+ de novo assembled genomes of C. hominis. Remarkably, after quality genome filtering, a comprehensive phylogenomic analysis allowed us to discover that C. hominis encompasses two lineages with continental segregation. These lineages were named based on the observed continental distribution bias as C. hominis Euro-American (EA) and the C. hominis Afro-Asian (AA) lineages.

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.

Challenges for Cryptosporidium Population Studies

Cryptosporidiosis is ranked sixth in the list of the most important food-borne parasites globally, and it is an important contributor to mortality in infants and the immunosuppressed. Recently, the number of genome sequences available for this parasite has increased drastically. The majority of the sequences are derived from population studies of Cryptosporidium parvum and Cryptosporidium hominis, the most important species causing disease in humans. Work with this parasite is challenging since it lacks an optimal, prolonged, in vitro culture system, which accurately reproduces the in vivo life cycle. This obstacle makes the cloning of isolates nearly impossible. Thus, patient isolates that are sequenced represent a population or, at times, mixed infections. Oocysts, the lifecycle stage currently used for sequencing, must be considered a population even if the sequence is derived from single-cell sequencing of a single oocyst because each oocyst contains four haploid meiotic progeny (sporozoites). Additionally, the community does not yet have a set of universal markers for strain typing that are distributed across all chromosomes. These variables pose challenges for population studies and require careful analyses to avoid biased interpretation. This review presents an overview of existing population studies, challenges, and potential solutions to facilitate future population analyses.

Prevalence and Genetic Characterization of Cryptosporidium and Microsporidia Infecting Hematological Malignancy Patients

PURPOSE

To evaluate the prevalence of Cryptosporidium and Microsporidia, associated risk factors and species identification in patients with haematological malignancies (HM).

METHODS

A total of 148 consecutive patients with HM and 101 healthy subjects were evaluated for Cryptosporidium and Microsporidia using modified Kinyoun and modified Trichrome staining. Clinical, demographic and laboratory parameters were studied. The species of Cryptosporidium and Microsporidia were studied using PCR-RFLP.

RESULTS

Of 148 HM patients initially screened, 47 were excluded from the final analysis due to inadequate clinical records. Patients with HM [n = 101, 63 (62.4%) male] more often had Cryptosporidium than healthy subjects [n = 101, 65 (74.4%) male] [3/101 (3%) vs. 0/101 (0%), p = 0.02]. Two of 101 (2%) HM patients and none of the healthy subjects had Microsporidia (p = 0.155). Diarrhea was more prevalent in HM patients with Cryptosporidium than those without [3, 100% vs. 39/96, 40.62%; p = 0.04). Both patients infected with Microsporidia presented with persistent diarrhea and fever. Cryptosporidium hominis was identified in all the three HM patients. Enterocytozoon bieneusi was identified in one HM patient infected with Microsporidia, which was classified as genotype Ind2.

CONCLUSION

Cryptosporidium and Microsporidia may infect HM patients leading to overwhelming diarrhea. The commonest species of Cryptosporidium and Microsporidia found to infect HM patients are C. hominis and E. bieneusi.

Zoonotic parasites in farmed exotic animals in China: Implications to public health

Several species of wild mammals are farmed in China as part of the rural development and poverty alleviation, including fur animals, bamboo rats, and macaque monkeys. Concerns have been raised on the potential dispersal of pathogens to humans and other farm animals brought in from native habitats. Numerous studies have been conducted on the genetic identity and public health potential of Cryptosporidium spp., Giardia duodenalis, and Enterocytozoon bieneusi in these newly farmed exotic animals. The data generated have shown a high prevalence of the pathogens in farmed wildlife, probably due to the stress from the short captivity and congregation of large numbers of susceptible animals. Host adaptation at species/genotype and subtype levels has reduced the potential for cross-species and zoonotic transmission of pathogens, but the farm environment appears to favor the transmission of some species, genotypes, and subtypes, with reduced pathogen diversity compared with their wild relatives. Most genotypes and subtypes of the pathogens detected appear to be brought in from their native habitats. A few of the subtypes have emerged as human pathogens. One Health measures should be developed to slow the dispersal of indigenous pathogens among farmed exotic animals and prevent their spillover to other farm animals and humans.

Molecular Epidemiology of Human Cryptosporidiosis in Low- and Middle-Income Countries

Cryptosporidiosis is one of the most important causes of moderate to severe diarrhea and diarrhea-related mortality in children under 2 years of age in low- and middle-income countries. In recent decades, genotyping and subtyping tools have been used in epidemiological studies of human cryptosporidiosis. Results of these studies suggest that higher genetic diversity of Cryptosporidium spp. is present in humans in these countries at both species and subtype levels and that anthroponotic transmission plays a major role in human cryptosporidiosis. Cryptosporidium hominis is the most common Cryptosporidium species in humans in almost all the low- and middle-income countries examined, with five subtype families (namely, Ia, Ib, Id, Ie, and If) being commonly found in most regions. In addition, most Cryptosporidium parvum infections in these areas are caused by the anthroponotic IIc subtype family rather than the zoonotic IIa subtype family. There is geographic segregation in Cryptosporidium hominis subtypes, as revealed by multilocus subtyping. Concurrent and sequential infections with different Cryptosporidium species and subtypes are common, as immunity against reinfection and cross protection against different Cryptosporidium species are partial. Differences in clinical presentations have been observed among Cryptosporidium species and C. hominis subtypes. These observations suggest that WASH (water, sanitation, and hygiene)-based interventions should be implemented to prevent and control human cryptosporidiosis in low- and middle-income countries.

Comparative genomics revealed adaptive admixture in Cryptosporidium hominis in Africa

Cryptosporidiosis is a major cause of diarrhoeal illness among African children, and is associated with childhood mortality, malnutrition, cognitive development and growth retardation. Cryptosporidium hominis is the dominant pathogen in Africa, and genotyping at the glycoprotein 60 (gp60) gene has revealed a complex distribution of different subtypes across this continent. However, a comprehensive exploration of the metapopulation structure and evolution based on whole-genome data has yet to be performed. Here, we sequenced and analysed the genomes of 26 C. hominis isolates, representing different gp60 subtypes, collected at rural sites in Gabon, Ghana, Madagascar and Tanzania. Phylogenetic and cluster analyses based on single-nucleotide polymorphisms showed that isolates predominantly clustered by their country of origin, irrespective of their gp60 subtype. We found a significant isolation-by-distance signature that shows the importance of local transmission, but we also detected evidence of hybridization between isolates of different geographical regions. We identified 37 outlier genes with exceptionally high nucleotide diversity, and this group is significantly enriched for genes encoding extracellular proteins and signal peptides. Furthermore, these genes are found more often than expected in recombinant regions, and they show a distinct signature of positive or balancing selection. We conclude that: (1) the metapopulation structure of C. hominis can only be accurately captured by whole-genome analyses; (2) local anthroponotic transmission underpins the spread of this pathogen in Africa; (3) hybridization occurs between distinct geographical lineages; and (4) genetic introgression provides novel substrate for positive or balancing selection in genes involved in host-parasite coevolution.

Comparative genomic analysis of the principal Cryptosporidium species that infect humans

Cryptosporidium parasites are ubiquitous and can infect a broad range of vertebrates and are considered the most frequent protozoa associated with waterborne parasitic outbreaks. The intestine is the target of three of the species most frequently found in humans: C. hominis, C. parvum, and. C. meleagridis. Despite the recent advance in genome sequencing projects for this apicomplexan, a broad genomic comparison including the three species most prevalent in humans have not been published so far. In this work, we downloaded raw NGS data, assembled it under normalized conditions, and compared 23 publicly available genomes of C. hominis, C. parvum, and C. meleagridis. Although few genomes showed highly fragmented assemblies, most of them had less than 500 scaffolds and mean coverage that ranged between 35X and 511X. Synonymous single nucleotide variants were the most common in C. hominis and C. meleagridis, while in C. parvum, they accounted for around 50% of the SNV observed. Furthermore, deleterious nucleotide substitutions common to all three species were more common in genes associated with DNA repair, recombination, and chromosome-associated proteins. Indel events were observed in the 23 studied isolates that spanned up to 500 bases. The highest number of deletions was observed in C. meleagridis, followed by C. hominis, with more than 60 species-specific deletions found in some isolates of these two species. Although several genes with indel events have been partially annotated, most of them remain to encode uncharacterized proteins.

Comparative analyses and structural insights of new class glutathione transferases in Cryptosporidium species

Cryptosporidiosis, caused by protozoan parasites of the genus Cryptosporidium, is estimated to rank as a leading cause in the global burden of neglected zoonotic parasitic diseases. This diarrheal disease is the second leading cause of death in children under 5 years of age. Based on the C. parvum transcriptome data, glutathione transferase (GST) has been suggested as a drug target against this pathogen. GSTs are diverse multifunctional proteins involved in cellular defense and detoxification in organisms and help pathogens to alleviate chemical and environmental stress. In this study, we performed genome-wide data mining, identification, classification and in silico structural analysis of GSTs in fifteen Cryptosporidium species. The study revealed the presence three GSTs in each of the Cryptosporidium species analyzed in the study. Based on the percentage identity and comprehensive comparative phylogenetic analysis, we assigned Cryptosporidium species GSTs to three new GST classes, named Vega (ϑ), Gamma (γ) and Psi (ψ). The study also revealed an atypical thioredoxin-like fold in the C. parvum GST1 of the Vega class, whereas C. parvum GST2 of the Gamma class and C. melagridis GST3 of the Psi class has a typical thioredoxin-like fold in the N-terminal region. This study reports the first comparative analysis of GSTs in Cryptosporidium species.

Molecular characterization of Cryptosporidium spp. from patients with diarrhoea in Lusaka, Zambia

Cryptosporidium is a major etiological agent of diarrhoeal diseases among children and immune-compromised individuals in sub-Saharan African countries. We conducted a study to determine the prevalence and genetic characteristics of Cryptosporidium spp. in stool samples from patients with diarrhoea who presented at the University Teaching Hospital in Lusaka, Zambia. Cryptosporidium species and subtypes from 71 microscopically confirmed cryptosporidiosis stool samples collected between 2017 and 2019 were determined by polymerase chain reaction followed by partial sequencing of the small subunit rRNA and 60-kDa glycoprotein (gp60) gene. Additionally, data for the period between 2014 and 2019 were reviewed and analysed for cryptosporidiosis seasonal and age distribution. Cryptosporidium was more prevalent in the rainy season. The highest number of cases was reported among the 1–4 year age group. By sequence analysis of the 71 positive isolates, Cryptosporidium hominis (n = 42; 59.2%), C. parvum (n = 27; 38%), C. felis (n = 1; 1.4%), and C. meleagridis (n = 1; 1.4%) were identified. Four C. hominis subtype families (Ia, Ib, Id, and Ie) and three C. parvum subtype families (IIc, IIe, and IIs) were identified. The most frequent subtypes were IeA11G3T3 (n = 20; 28.2%), IIcA5G3 (n = 12; 16.9%), IIeA12G1 (n = 11; 15.5%) and IaA30R3 (n = 10; 14.1%). The observed species/subtypes of C. hominis and C. parvum indicated that the infection was mainly transmitted through the anthroponotic route. The identification of C. felis and C. meleagridis suggests that an atypical zoonotic transmission cycle also exists.

Comparative genomic analysis of three intestinal species reveals reductions in secreted pathogenesis determinants in bovine-specific and non-pathogenic Cryptosporidium species

The three common intestinal Cryptosporidium species in cattle differ significantly in host range, pathogenicity and public health significance. While Cryptosporidium parvum is pathogenic in pre-weaned calves and has a broad host range, C. bovis and C. ryanae are largely non-pathogenic and bovine-specific species in post-weaned calves. Thus far, only the genome of C. parvum has been sequenced. To improve our understanding of the genetic determinants of biological differences among Cryptosporidium spcies, we sequenced the genomes of C. bovis and C. ryanae and conducted a comparative genomics analysis. The genome of C. bovis has a gene content and organization more similar to C. ryanae than to other Cryptosporidium species sequenced to date; the level of similarity in amino acid and nucleotide sequences between the two species is 75.2 and 69.4 %, respectively. A total of 3723 and 3711 putative protein-encoding genes were identified in the genomes of C. bovis and C. ryanae, respectively, which are fewer than the 3981 in C. parvum. Metabolism is similar among the three species, although energy production pathways are further reduced in C. bovis and C. ryanae. Compared with C. parvum, C. bovis and C. ryanae have lost 14 genes encoding mucin-type glycoproteins and three for insulinase-like proteases. Other gene gains and losses in the two bovine-specific and non-pathogenic species also involve the secretory pathogenesis determinants (SPDs); they have lost all genes encoding MEDLE, FLGN and SKSR proteins, and two of the three genes for NFDQ proteins, but have more genes encoding secreted WYLE proteins, secreted leucine-rich proteins and GPI-anchored adhesin PGA18. The only major difference between C. bovis and C. ryanae is in nucleotide metabolism. In addition, half of the highly divergent genes between C. bovis and C. ryanae encode secreted or membrane-bound proteins. Therefore, C. bovis and C. ryanae have gene organization and metabolic pathways similar to C. parvum, but have lost some invasion-associated mucin glycoproteins, insulinase-like proteases, MEDLE secretory proteins and other SPDs. The multiple gene families under positive selection, such as helicase-associated domains, AMP-binding domains, protein kinases, mucins, insulinases and TRAPs could contribute to differences in host specificity and pathogenicity between C. parvum and C. bovis. Biological studies should be conducted to assess the contribution of these copy number variations to the narrow host range and reduced pathogenicity of C. bovis and C. ryanae.

Expression and Functional Studies of INS-5, an Insulinase-Like Protein in Cryptosporidium parvum

The small Cryptosporidium genome (∼9 Mb) has over 20 copies of genes encoding insulinase-like proteases (INS), suggesting that these enzymes may have important biological functions in the pathogen and could be developmentally regulated. In this study, INS-5, a unique member of the INS family in Cryptosporidium parvum, was cloned and expressed in Escherichia coli BL21 (DE3). In addition to the predicted INS-5 of ∼78 kDa, smaller fragments of ∼70, ∼55, and ∼30 kDa were simultaneously generated. After purification through a nickel-nitrilotriacetic acid affinity column, the full recombinant protein obtained was used to prepare polyclonal antibodies. Antibodies raised against INS-5 recognized the recombinant protein and native protein in sporozoite extracts. Further characterization of INS-5 included qRT-PCR assessment of gene expression; immunofluorescence localization of the protein expression in sporozoites, merozoites, and other developmental stages; and neutralization of invasion of C. parvum in vitro. The results obtained indicated that although INS-5 was expressed in sporozoites and merozoites, the high gene expression was from 36 to 48 h of the in vitro culture after invasion. Anti-INS-5 antibodies partially neutralized the invasion (inhibition rate = 38.5%). Results of this study suggest that INS-5 plays some role in the invasion and growth of C. parvum.

Genetic basis for virulence differences of various Cryptosporidium parvum carcinogenic isolates

Cryptosporidium parvum is known to cause life-threatening diarrhea in immunocompromised hosts and was also reported to be capable of inducing digestive adenocarcinoma in a rodent model. Interestingly, three carcinogenic isolates of C. parvum, called DID, TUM1 and CHR, obtained from fecal samples of naturally infected animals or humans, showed higher virulence than the commercially available C. parvum IOWA isolate in our animal model in terms of clinical manifestations, mortality rate and time of onset of neoplastic lesions. In order to discover the potential genetic basis of the differential virulence observed between C. parvum isolates and to contribute to the understanding of Cryptosporidium virulence, entire genomes of the isolates DID, TUM1 and CHR were sequenced then compared to the C. parvum IOWA reference genome. 125 common SNVs corresponding to 90 CDSs were found in the C. parvum genome that could explain this differential virulence. In particular variants in several membrane and secreted proteins were identified. Besides the genes already known to be involved in parasite virulence, this study identified potential new virulence factors whose functional characterization can be achieved through CRISPR/Cas9 technology applied to this parasite.

Update on Cryptosporidium spp.: highlights from the Seventh International Giardia and Cryptosporidium Conference

While cryptosporidiosis is recognized as being among the most common causes of human parasitic diarrhea in the world, there is currently limited knowledge on Cryptosporidium infection mechanisms, incomplete codification of diagnostic methods, and a need for additional therapeutic options. In response, the Seventh International Giardia and Cryptosporidium Conference (IGCC 2019) was hosted from 23 to 26 June 2019, at the Rouen Normandy University, France. This trusted event brought together an international delegation of researchers to synthesize recent advances and identify key research questions and knowledge gaps. The program of the interdisciplinary conference included all aspects of host-parasite relationships from basic research to applications to human and veterinary medicine, and environmental issues associated with waterborne parasites and their epidemiological consequences. In relation to Cryptosporidium and cryptosporidiosis, the primary research areas for which novel findings and the most impressive communications were presented and discussed included: Cryptosporidium in environmental waters, seafood, and fresh produce; Animal epidemiology; Human cryptosporidiosis and epidemiology; Genomes and genomic evolution encompassing: Comparative genomics of Cryptosporidium spp., Genomic insights into biology, Acquiring and utilizing genome sequences, Genetic manipulation; Host-parasite interaction (immunology, microbiome); and Diagnosis and treatment. High quality presentations discussed at the conference reflected decisive progress and identified new opportunities that will engage investigators and funding agencies to spur future research in a “one health” approach to improve basic knowledge and the clinical and public health management of zoonotic cryptosporidiosis.

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.

Divergent Copies of a Cryptosporidium parvum-Specific Subtelomeric Gene

Subtype families of Cryptosporidium parvum differ in host range, with IIa and IId being found in a broad range of animals, IIc in humans, and IIo and IIp in some rodents. Previous studies indicated that the subtelomeric cgd6_5520-5510 gene in C. parvum is lost in many Cryptosporidium species, and could potentially contribute to the broad host range of the former. In this study, we identified the presence of a second copy of the gene in some C. parvum subtype families with a broad host range, and showed sequence differences among them. The sequence differences in the cgd6_5520-5510 gene were not segregated by the sequence type of the 60 kDa glycoprotein gene. Genetic recombination appeared to have played a role in generating divergent nucleotide sequences between copies and among subtype families. These data support the previous conclusion on the potential involvement of the insulinase-like protease encoded by the subtelomeric cgd6_5520-5510 gene in the broad host range of C. parvum IIa and IId subtypes.

Genetic Diversity of Cryptosporidium hominis in a Bangladeshi Community as Revealed by Whole-Genome Sequencing

We studied the genetic diversity of Cryptosporidium hominis infections in slum-dwelling infants from Dhaka over a 2-year period. Cryptosporidium hominis infections were common during the monsoon, and were genetically diverse as measured by gp60 genotyping and whole-genome resequencing. Recombination in the parasite was evidenced by the decay of linkage disequilibrium in the genome over <300 bp. Regions of the genome with high levels of polymorphism were also identified. Yet to be determined is if genomic diversity is responsible in part for the high rate of reinfection, seasonality, and varied clinical presentations of cryptosporidiosis in this population.

Molecular typing of Cryptosporidium in Israel

Cryptosporidium is a protozoan parasite associated with gastrointestinal illness. In immune-compromised individuals, the infection may become life-threatening. Cryptosporidiosis is a mandatory-reported disease but little was known about its prevalence and associated morbidity in Israel. Currently, laboratory diagnosis is based on microscopy or copro-antigen tests and the disease is underreported. Molecular assays, which are more sensitive and specific, are now increasingly used for identification and screening. Here, the molecular epidemiology of cryptosporidiosis is explored for the first time. Samples from 33 patients infected during an outbreak of 146 laboratory confirmed cases that occurred in Haifa and Western Galilee in 2015 were genotyped, as well as samples from 36 patients sporadically infected during 2014–2018 in different regions. The results suggest that Cryptosporidium subtypes found in Israel are more similar to those reported in the neighboring countries Jordan and Egypt than in European countries. C. hominis was the predominant species in the center and the north of Israel, implicating human-to-human transmission. C. hominis IeA11G3T3 was the most prevalent subtype contributing to morbidity.

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.