Abstract
It is critical to characterize changes in the structure and composition of the host fungal community in natural Cryptosporidium infection, because it gives the possible overview of gut microbiome in host homeostasis and disease progression. A total of 168 rectal fecal samples were collected and examined using nPCR. The positive samples were double-checked using 18S rDNA high-throughput sequencing. After confirmation, ITS high-throughput sequencing was utilized to investigate the fungal community’s response to natural Cryptosporidium infection. Results showed that a total of three positive samples (1.79%) were identified with an increased abundance of fungi associated with health hazards, such as class Dothideomycetes, families, i.e., Cladosporiaceae, Glomerellaceae, and genera, i.e., Wickerhamomyces, Talaromyces, Cladosporium, Dactylonectria, and Colletotrichum. On the contrary, taxa associated with favorable physiological effects on the host were shown to have the reverse impact, such as families, i.e., Psathyrellaceae, Pseudeurotiaceae and genera (Beauveria, Nigrospora, and Diversispora). For the first time, we evaluated the condition of natural Cryptosporidium infection in horses in Wuhan, China, and discovered distinct variations in the fungal microbiome in response to natural infection. It might prompt a therapy or prevention strategy to apply specific fungal microorganisms that are probably responsible for decreased susceptibility or increased resistance to infection.
Introduction
The gut is a diverse and complex ecosystem of microbes, encompassing viruses, bacteria, fungi, parasites, and other microorganisms that have been demonstrated to play a fundamental role in maintaining the health or disease state of the host (). The intestinal immune system has evolved to adequately meet the presence of various microorganisms, including beneficial, detrimental, and harmless microbes (). The host intestinal immune system has to make accurate assessments to accept or reject the new microbes in each condition. Meanwhile, the incoming organism has also evolved to maximize opportunities for acceptance (). Hence, investigating the interplay between various diseases and gut microbiota remains the key focus. Evidence shows that the symbiotic relationship between host and resident microbial taxa in the gut contributes to extracting essential nutrients as a first-line for gut pathogens and shaping the immune system []. There is an increasing realization that the influences between the parasites and gut microbiota are intertwined. reported that parasitic infection changed the structure and composition of the gut microbiome, and conversely, the stable microbiota affects intestinal parasite colonization and persistence in the gut. Infection with parasites can disrupt or restructure the host’s gut communities, both in vertebrates and invertebrates (; ). On the contrary, a disrupted intestinal microbiome can increase the susceptibility to some parasitic diseases, e.g., malaria and Giardia infections (; ).
Interaction between intestinal organisms in host individuals is one of the key determinants of dynamics in intestinal communities. These multitudes of interactions are both detrimental and beneficial to their hosts, affecting symbiont communities’ ecological niche and evolution (; ). Recently, few observational studies reported that the interactions between parasites and microbiota have become more feasible on the base of new sequencing technologies. Some have reported no bacterial community changes in infection, whereas others did (; ). In all conditions with significant alterations, a higher diversity of bacteria has been associated with infections. At the same time, the impact on operational taxonomic units (OTUs) have also been identified to have parasitic specificity (). Even though most research is focused on the shifts of gut bacterial communities in parasites infections. Growing evidence confirmed the undisputable effects of fungi in driving various gut-associated and metabolic diseases (). Compared with bacterial communities, fungal presence is relatively insignificant that limits the in-depth understanding of this “rare biosphere.”
Nonetheless, deep-sequencing technologies have revealed the complexity of the fungi microflora that existed on the gut mucosal surfaces. In addition, accumulated evidence showed that specific fungi played an important role in modulating host immune response and maintaining host dynamic gut microflora changes (), which interacted with the host immune system through the innate immune receptor Dectin-1. More profoundly, gut fungal communities may be the reservoir for opportunistic pathogens when a host is immunocompromised (). For example, previous research reported that the fungal burden increased in the gut of the Crohn’s Disease patients (). The evidence clearly indicated that fungi microflora was essential in host homeostasis and disease development, suggesting a potential role of gut fungi in host homeostasis and disease development in mammals. So, it is essential to characterize the shifts of structure and composition of host fungi community in natural parasitic infection.
Cryptosporidium spp. are opportunistic apicomplexan pathogenic parasites, infecting gut epithelial cells and causing debilitating gastrointestinal and diarrheal ailment in the host (). Cryptosporidiosis can be life-threatening in immunocompromised individuals. High morbidity and mortality have been reported in patients with AIDS/HIV (). One single Cryptosporidium oocyst can trigger an infection in immunosuppressed individuals (). The infection can be spread directly from the infective oocysts in the feces or indirectly from vegetables or other food items (contaminated by the feces of infected hosts) (). The findings of Cryptosporidium spp. in the infected horses and donkeys, reared as companion animals suggest that these could be a reason for transmission to humans (). In equines, most reports on cryptosporidiosis were related to Immune Deficiency Syndrome, e.g., Arabian foal (). Moreover, it has been reported that the infection by Cryptosporidium in healthy foals occurred in the United States, Spain, etc. (). A remarkable phenomenon is the difference in the severity and susceptibility of Cryptosporidium infection (). These conditions may be multifactorial, including host species, immunological status, and differences in the host gut microflora.
Therefore, the study aimed to test the hypothesis that gut fungal microbiota as a holobiont could respond to natural Cryptosporidium infection in horses. Here, we used specific polymerase chain reaction (PCR) amplification and 18S rDNA high-throughput sequencing to reveal the occurrence of natural Cryptosporidium infection. In addition, we screened positively infected feces to identify changes in the structure and composition of fungi to explore any correlation to Cryptosporidium infection. This study was the first report that revealed natural Cryptosporidium infection in horses in Wuhan, China, and response of gut fungal microbiota to infection.
Materials and Methods
Collection of Fecal Samples
The sterile fecal collection tubes were used to collect 168 rectal fecal samples from horses in Wuhan, China. Each fecal sample (about 2 g) obtained was immediately transferred into a fecal collection tube and rapidly frozen in dry ice before storing at −80°C in the clinical lab of Huazhong Agricultural University. The study procedures were approved by the Ethics Committee of Huazhong Agricultural University and Animals Research Centre of Hubei Province.
DNA Extraction
Before DNA extraction, samples were vortexed to homogenization. The genomic DNA (gDNA) was extracted from feces using the commercially available product E.Z.N.A® Stool DNA Kit (Omega Bio-tek, Inc., United States) following the manufacturer’s recommended procedure. Then, the quality and purity of gDNA were checked using spectrophotometry and 1% agarose gel electrophoresis before storage at –20°C.
PCR Amplification and DNA Agarose Gel Electrophoresis
Following the reports described previously, we amplified the 18S rDNA gene of Cryptosporidium using two-step nested PCR (). In the primary PCR, a pair of primers 18SiCF2: 5′-GACATATCATTCAAGTTTCTGACC-3′ and 18SiCR2: 5′-CTGAAGGAGTAAGGAACAACC-3′ were used to obtain a PCR product of 763 bp. Total 25 μl PCR reaction mixture contained 1 μl Taq, total 2 μl forward and reverse primers, 2.5 μl PCR Buffer (10×), 2.5 μl DNA, 5 μl dNTPs (2.5 mM), and 12 μl sterile distilled water. The annealing temperature was 58°C for 30 s during 45 PCR cycles. During the second PCR, a fragment of 587 bp was amplified using forward primer 18SiCF1 (5′-CCTATCAGCTTTAGACGGTAGG-3′) and reverse primer 18SiCR1 (5′-TCTAAGAATTTCACCTCTGACTG-3′). The total 50 μl reaction volume contained a total of 4 μl forward and reverse primers, 2 μl Taq, 5 μl primary PCR product, 5 μl PCR Buffer (10 × ), 10 μl dNTPs (2.5 mM), and 24 μl sterile distilled water. The conditions for two-step PCR were identical. Ultimately, all PCR products were analyzed using 1% agarose gels (w/v).
Double-Test to Verify Cryptosporidium Spp. Infection Using 18S rDNA High-Throughput Sequencing
The Cryptosporidium-positive feces verified by the nested PCR were double-checked to confirm infection using 18S rDNA high-throughput sequencing. Specific primers based on the V4 region were synthesized for amplification (F: 5′-CCAGCASCYGCGGTAATTCC-3′, and R: 5′-ACTTTCGTTCTTGATYRA-3′). The annealing temperature was 56°C for 30 s during 35 PCR cycles. All the PCR amplification products were performed 1% agarose gel electrophoresis, then the gel recovery kit (Axygen, United States) was used to recover targeted fragments. Subsequently, purified all PCR products to remove the unspecific products for library preparation. The average molecular length and quantity of amplifications were validated using the Agilent bioanalyzer instrument (Agilent, United States) and real-time quantitative PCR (EvaGreen™). The final validated libraries were loaded on Illumina Novaseq 6000 platform (Illumina, San Diego, CA, United States) for high-throughput sequencing using a paired-end configuration.
ITS Genes Amplification and Sequencing
Both Cryptosporidium-positive (YangS group) and Cryptosporidium-negative (YinC group) feces were used to perform ITS genes amplification. The ITS1 variable region of the fungal ITS gene was amplified using forward primers containing the sequence 5′-CTTGGTCATTTAGAGGAAGTAA-3′ and reverse primers containing the sequence 5′-GCTGCGTTCTTCATCGATGC-3′. Both the primers were tagged with an Illumina adapter for PCR amplification. The annealing temperature was 58°C during 35 PCR cycles. A total of 30 ng verified gDNA was used to perform PCR per reaction. The conditions of purifying PCR products and preparing sequencing libraries were identical with 18S rDNA high-throughput sequencing. Ultimately, paired-end sequencing was performed on qualified libraries using the Illumina Novaseq platform.
Bioinformatics and Statistical Analysis
Performed the following pre-processed steps on the raw FASTQ format sequences obtained by high-throughput sequencing to eliminate the low quality and adapter pollution for acquiring more accurate and reliable Clean Reads: (1) Quality filtering: those reads (<200 bp) and low quality were filtered using Trimmomatic v0.33 software; (2) To obtain Clean Sequences that do not contain primer sequences, the primer sequences were identified and removed using Cutadapt 1.9.1 software; (3) Both the ends of Clean Reads were stitched through overlap using Usearch v10 software; (4) The UCHIME v4.2 software was used to identify and remove chimera sequences to obtain the final Effective Reads. The obtained high-quality Effective Reads were clustered into operational taxonomic units (OTUs) against the SILVA 132 (18S rDNA high-throughput sequencing) or UNITE (ITS genes high-throughput sequencing) databases at over 97% sequence similarity. A representative sequence per OTU has executed classification and phylogenetic analysis. The Rarefaction, Shannon, and Species accumulation curves evaluated the quality and depth of current sequencing, illustrated by R (v3.0.3) software. Alpha-diversity indexes contained the Chao1 and ACE (calculates species richness), Shannon and Simpson indexes (calculates species diversity), which were processed and calculated by QIIME2 software (v1.7.0). Beta diversity was performed by principal coordinates analysis (PCoA) and combination of (unweighted pair-group method with arithmetic mean) UPGMA clustering tree and histogram, which reflected the difference and similarity of gut fungal community structure in positive and negative infection groups. Metastatic analysis was utilized to assess the fungal community difference between groups. GraphPad Prism (v7.0) and SPSS software (v17.0) were applied to statistical analysis. The values were expressed as the mean ± SD, with P-value < 0.05 identified as statistically significant.
Results
Cryptosporidium Spp. Prevalence in Horses
The positive samples of Cryptosporidium infection were preliminarily identified by nested PCR amplification of 18S rDNA genes. Positive samples exhibited a target band of about 587 bp on the agarose gel after staining with ethidium bromide (Figure 1A). The 18S rDNA PCR analysis of 168 rectal fecal samples showed the occurrence of natural Cryptosporidium infection in horses of Wuhan, China, to be 1.79% found in three samples (Figure 1B).
FIGURE 1
Analysis of 18S rDNA High-Throughput Sequencing Data
The gDNA of feces initially identified as positive for infection was further performed 18S rDNA high-throughput sequencing to confirm Cryptosporidium spp. infection. After optimizing the original data, 236,026 effective reads were acquired, and 78,911, 77,938, and 79,177 high-quality valid sequences in the WS3, SS2, and SS7 samples, respectively (Table 1). Based on 97% similarity, 59, 61, and 95 OTUs were separately acquired from WS3, SS2, and SS7, respectively (Figure 1C). Figures 1D–F showed the taxonomic information of the top 10 abundant species identified at the phylum and genus levels. Results showed phylum Apicomplexa and genus Cryptosporidium, confirming natural Cryptosporidium infection in horses.
TABLE 1
| Sample ID | Raw reads | Clean reads | Effective reads | Effective (%) |
| WS3 | 79,869 | 79,596 | 78,911 | 98.8 |
| SS2 | 79,954 | 79,679 | 77,938 | 97.48 |
| SS7 | 79,935 | 79,624 | 79,177 | 99.05 |
The data of 18S rDNA high-throughput sequencing.
ITS High-Throughput Sequencing Data Analysis
The ITS amplicon sequence was performed to investigate the alterations in gut fungal community driven by the natural Cryptosporidium spp. infection. After the quality control, a total of 478,486 clean reads were obtained, with an average value of 79,748 counts per sample (Table 2). The rarefaction curve, Shannon curve, and species accumulation curve were wide and tended to be flat, suggesting that the current sequencing quantity and depth were reasonable and sufficient to reflect the fungal richness and diversity included in the samples (Figures 2A–C). We recognized a total of 1,519 OTUs based on 97% sequence similarity. Among them, 784 OTUs were acquired for group YangS, and 735 OTUs were acquired for group YinC; among them, 434 OTUs were core OTUs, accounting for 55.36% and 59.05% of the total OTUs in group YangS and group YinC, respectively (Figures 3A,B).
TABLE 2
| Sample ID | Raw reads | Clean reads | Effective reads | Effective (%) |
| SS2 | 79,834 | 79,618 | 78,791 | 98.69 |
| SS7 | 80,182 | 79,992 | 78,801 | 98.28 |
| MC10 | 80,083 | 79,887 | 79,046 | 98.71 |
| MC15 | 79,673 | 79,453 | 78,475 | 98.5 |
| WS3 | 80,042 | 79,799 | 79,150 | 98.89 |
| MC30 | 79,978 | 79,737 | 78,495 | 98.15 |
The data of ITS genes high-throughput sequencing.
FIGURE 2
FIGURE 3
Shifts in Gut Fungal Community Diversity and Structure With the Infection of Cryptosporidium Spp.
The Chao1 and ACE indices showed that the difference in species richness between YangS and YinC groups was insignificant (Figure 3C). A significantly lower Shannon diversity index was found in the YinC group than in the YangS group (p < 0.05), indicating that the fungal community in the YangS group had higher diversity than the YinC group. However, the other diversity index, Simpson index had no significant difference between these two groups (p > 0.05). In addition, the values of coverage index estimates showed excellent coverage in both the groups, ranging from 99.97 to 99.98%. The PCoA based on the fungal OTUs showed that all the fecal samples from the YangS group were clustered together and separated from the YinC group and indicated the variation in the gut fungal community with natural Cryptosporidium spp. infection (Figure 2D) was supported by visualizing UPGMA clustering tree associations with horses with or without infected Cryptosporidium spp. (Figure 2E).
Mapping different taxonomical levels (class, family, and genus levels) relative abundances onto the species distribution histogram revealed gut fungal community structure in horses with or without infected Cryptosporidium spp. (Figures 4A–C). Significantly, the same species hosts generally had similar phylum-level abundances, e.g., phyla Mortierellomycota and Ascomycota were dominant in fungal communities regardless of Cryptosporidium spp. infection, which consisted of over 25% of total tags on average. Besides these two dominant phyla, the abundantly present phylum was Basidiomycota (16.1%) in the YangS group, while Neocallimastigomycota (15%) was observed as the predominant phylum in YinC group. Classes Eurotiomycetes (8.44, 13.40%), Mortierellomycetes (28, 24.9%), and Sordariomycetes (10.8, 10.3%) were the most abundant in both YangS and YinC groups, whereas Agaricomycetes (9.6%) and Neocallimastigomycetes (15%) were enriched in the YangS and YinC groups, respectively. The outcomes of the taxonomic classification level of the gut fungal community identified one dominating family at its level, Mortierellaceae, accounting for 24.6% of identified taxonomic sequences. Besides, Aspergillaceae (3.99%) and Neocallimastigaceae (15%) dominated the YangS and YinC groups, respectively. Using genus-level cluster analysis, Mortierella was the predominant genus of total sequences in both groups, consisting of over 24.6% of total tags, whereas Aspergillus was secondary (over 3.24%). Moreover, the most abundant genera in the YangS group were Archaeorhizomyces (2.82%), Wallemia (2.29%), Anaeromyces (2.14%), Fusarium (1.92%), and Inocybe (1.33%). By comparison, Caecomyces (2.67%), Fusarium (2.62%), Zygosaccharomyces (2.53%), Archaeorhizomyces (2.06%), and Piromyces (1.37%) were enriched in the YinC group. Furthermore, the heatmap results exhibited relative species abundance at different levels (class, family, and genus) were further validated by analytical strategy (Figures 4D–F).
FIGURE 4
Changes in Gut Fungal Composition of Horses With the Infection of Cryptosporidium Spp.
Metastatic analysis was used to explore the significant changes in the gut fungal community at class, family, and genus levels to undertake exploratory shifts in the gut fungal microbiota of natural Cryptosporidium spp. infection (Figure 5). At the class level, 2 taxa were significantly different when comparing YangS and YinC in the exploration panel (p < 0.05), which were Dothideomycetes and Archaeorhizomycetes. In comparing all the taxa between YangS and YinC, four families had different abundances in the combined panel; families Glomerellaceae and Cladosporiaceae were enriched in YangS, while families Psathyrellaceae and Pseudeurotiaceae were enriched in YinC (all: p < 0.05). Results of the genus-level cluster showed that 304 genera were recognized across all the sequences. Several taxa at the genus level were found to have significant infection-related variations. The high signals of the YangS group of fungi were enriched at these genera levels (Dactylonectria, Trichocladium, Colletotrichum, Wickerhamomyces, Talaromyces, Cladosporium, Serendipita, and Strelitziana). Conversely, Beauveria, Hannaella, Subulicystidium, Nigrospora, Pseudeurotium, Symmetrospora, and Diversispora were represented significantly in YinC group compared with YangS group (p < 0.05 or p < 0.01).
FIGURE 5
Discussion
Microbiota typically exceeds intestinal parasites in biomass, while commensal fungi and parasites occupy different taxonomic spaces; both have responded to evolutionary capability by developing some strategies for modulating host immunity (). It is evident that they have established dialog with a common theme of creating new homeostasis in the host gut. Specifically, parasitic infection changes the structure and composition of the gut microbiota, and conversely, the inherent microbiota affects invasions of new parasite species within hosts (). The bacterial microbiota patterns associated with susceptibility to Cryptosporidium infection have been identified in humans (), while the theme of the fungal–parasite interaction was relatively unexplored. In the current study, we first performed an epidemiological survey to investigate the condition of natural Cryptosporidium infection in horses of Wuhan, P. R. China. By doing so, we obtained the positively infected feces to investigate the gut microbiome in which fungal communities have been changed by exposure to Cryptosporidium spp.
Cryptosporidium spp. are protozoan parasites of medical importance that causes moderate-to-severe diarrhea and death in varieties of vertebrate hosts (). Cryptosporidium spp. are now considered important waterborne parasites worldwide. The positively infected horses pose a serious public health threat in many countries, for example, the United States and France (; ). The main symptoms includes diarrhea, abdominal pain, vomiting and low-grade fever (). In 1976, the first case of human cryptosporidiosis was reported in patient with severe watery diarrhea (). The severity of Cryptosporidium infection depends usually on the host factors e.g., immunocompetent individuals suffer a transient self-limiting illness (). When Cryptosporidium infection was considered as a major inducement of diarrhea in immunocompetent individuals, the public health meaning of it became apparent (). Evidence from epidemiology and microbiology revealed the direct and indirect transmission routes of Cryptosporidium spp.: (1) direct infection occurs by the fecal–oral route, including animal-to-human (zoonotic), animal-to-animal, human-to-human, and human-to-animal transmissions (); (2) Indirect infection occurs through environmental contamination, involving contact with Cryptosporidium fecal contaminated material (). The infection rate in adult equines ranges from 0.33 to 9% () with 1.79% prevalence of natural Cryptosporidium spp. The current study was consistent with previous reports in Algeria (2.3%) and China’s southwest (1.8%) (; ). It was inconsistent with previous reports of prevalence in cattle and donkeys (; ). We inferred that the potential reason might due to the different animal species that had different immunity against Cryptosporidium spp. infection.
The composition and structure of the intestinal microbiome can be influenced by diet, age, environmental factors, and the role of protozoal parasites (; ; ; ), but few studies focusing on gut fungi have been performed with Cryptosporidium spp. infection. In neonatal mice, the innate responses requiring microbiota played a key role in protection against Cryptosporidium infection (). Although the impacts of microbiome parasite on the host were predictable in the long term, the potential competition or facilitation in the intestinal micro-ecosystem might happen in a much shorter timeframe (). After 3 days of Cryptosporidium infection, as seen in protein-malnourished mice, intestinal injury and inflammation were recorded (). The composition of gut microflora in monkeys also shifted significantly following Cryptosporidium infection (). According to a previous report, in all conditions with significant alteration, a higher diversity of bacteria seems to associate with infections (), which was consistent with the current results. More Shannon diversity index was found in natural Cryptosporidium infection horses. Moreover, six branches in the UPGMA evolutionary tree were sequestered, showing the segregation of fungal community structure and composition in Cryptosporidium infections. Such findings supported the hypothesis about the response of fungal microbiome to Cryptosporidium infection in horses ().
As mentioned previously, resident microbiota composition could affect the success rate of invading parasites; meanwhile, parasites actively secrete molecules affecting microbiota composition. This information conveyed that the changes in some functional microflora abundance exacerbate infection susceptivity or that particular gut condition, which infected Cryptosporidium spp., drives the shifts of the functional microorganisms. The significant enrichment of the classes Dothideomycetes and Archaeorhizomycetes contributed to the unique gut fungal signature found in YangS group in this study. Although research on the fungal microbiota is still in its early stages, various reports have revealed its potential effect on maintaining hosts’ homeostasis. It coexists with bacterial microorganisms and substantially expands gut microecology interacting with the intestinal immune system. Mainly fungal microorganisms were contributed to decomposing the lignocellulose in the gut and implementing this functionality through secretion of cell wall degradation enzymes and physical penetration (). The final metabolites of this procedure were mainly formate, hydrogen, and acetate (). Class Dothideomycetes is the largest taxa in the kingdom of fungi, many of which have evolved to an incredible diversity of lifestyles. The major ecological niche of this class is plant pathogens, and most food crops and feedstocks for biofuel production are its’ infected subjects (). While the ecology of class Archaeorhizomycetes is poorly resolved, their abundance in the soil is currently disputed. Even so, reported that these root-associated ascomycetes probably promote high Carbon turnover. At the family level, we found one taxa, affiliated with degraded hemicellulose, significantly enriched in the YinC group, i.e., Pseudeurotiaceae (). In contrast, classes Glomerellaceae and Cladosporiaceae had the opposite trend (enrichment in YangS group). Cladosporium species is an opportunistic human pathogen belonging to the class Cladosporiaceae, reported from cases of various superficial and invasive infections worldwide (). Coincidentally, leaf anthracnose is caused by genus Colletotrichum, a notorious necrotrophic fungus that belongs to the family Glomerellaceae (). These variations in the abundances of involved fungal microflora might be associated with infection development. The robust connection warrants further study.
At the genus level, Wickerhamomyces, Talaromyces, and Cladosporium were recorded as most abundant in positively infected horses (all: p < 0.05), probably associating with their physiological effects in the gut. For example, genus Talaromyces marneffei, an essential pathogenic thermally dimorphic fungus, have been reported of causing clinical systemic mycosis in patients with human immunodeficiency virus (HIV) or non-HIV infected with damaged cell-mediated immunity (). Although the genus Wickerhamomyces has been used in various biotechnology applications of agriculture, some studies have recognized this specie as an opportunistic human pathogen (). Remarkably, several microbial components were found connected to plant pathogen, i.e., Dactylonectria, Colletotrichum, and Cladosporium, in YangS group (). There is no knowledge about the physiological effects of plant pathogens in vertebrate intestines. From current data, we speculated that its enrichment in the gut might not be an accidental phenomenon, which is worthy of further advancement in the field of fungi. The substantial rise in taxa connected with health threats in Cryptosporidium infected horses was followed by a significant decrease in several beneficial fungi classifications. Beauveria was reported for its eminent function as a mycoinsecticide and tool in the integrated pest management (), which might consist with down-regulation of fungi diversity in negative infected horses because its highest abundance was recorded in YinC group (p < 0.01). The secondary metabolites produced from a piezotolerant fungus, Nigrospora spp., exhibited great antimicrobial and anticancer activities and has been applicated in pharmaceutical industries (), while Diversispora as an arbuscular mycorrhizal symbiosis, can stimulate growth and nutrient uptake in host plants (). Cryptosporidium spp. survival in the hosts’ gut must interact with intestinal constituents to compensate for deficiencies, i.e., lack of numerous metabolic systems (). Gut microbiome–parasite relationship was demonstrated to be a medium in modulating significant biochemical changes during cryptosporidiosis (). identified that specific changes in the gut microbiome allow for more favorable survival of Cryptosporidium parasite in mice (), hence this conclusion supported the current results. A question, which could only be speculated upon, is whether changes in the abundance of fungal microflora is inducer in Cryptosporidium infection; do the microflora changes represent an actual connection between the intestinal conditions and Cryptosporidium infection, or are they indirectly involved in infection? Anyway, it was speculated that the marked fungal microbiota alterations in YangS were involved in infection.
Conclusion
The present research focused on the prevalence of Cryptosporidium spp. infection in horses. According to our findings, the positive samples exhibited an increased abundance of fungi associated with health issues, while taxa associated with beneficial physiological effects on the host showed the opposite effect. Consequently, our findings suggest that certain modifications in the host fungal microbiome might enable Cryptosporidium parasites to survive. These findings might pave the way for creating a therapeutic or preventative approach that explicitly targets fungal microorganisms that produce infection susceptibility or resistance.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The raw data of current study had been submitted to NCBI Sequence Read Archive (SRA) with accession no: PRJNA797580.
Ethics statement
The animal study was reviewed and approved by Ethics Committee of Huazhong Agricultural University and Animals Research Centre of Hubei province.
Author contributions
YaW, KL, and YiW: research idea and methodology. XL, KD, HL, and XC: reagents, materials, and analysis tools. YaW: writing – original draft and preparation. MK and ZB: writing – review and editing. KL: visualization and supervision. All authors known and approved the final manuscript.
Funding
This study was supported by the start-up fund of the Nanjing Agricultural University (804131) and the start-up fund for distinguished scholars of the Nanjing Agricultural University (80900219).
Conflict of interest
KD was employed by China Tobacco Henan Industrial Co., Ltd., China. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AbdoliA.MaspiN. (2018). Commentary: estimates of global, regional, and national morbidity, mortality, and aetiologies of diarrhoeal diseases: a systematic analysis for the global burden of disease study 2015.Front. Med.5:11. 10.3389/fmed.2018.00011
2
AiveloT.NorbergA. (2018). Parasite-microbiota interactions potentially affect intestinal communities in wild mammals.J. Anim. Ecol.87438–447. 10.1111/1365-2656.12708
3
AlcockF.WebbC. T.DolezalP.HewittV.Shingu-VasquezM.LikicV. A.et al (2012). A small Tim homohexamer in the relict mitochondrion of Cryptosporidium.Mol. Biol. Evol.29113–122. 10.1093/molbev/msr165
4
ArumugamG. K.SrinivasanS. K.JoshiG.GopalD.RamalingamK. (2015). Production and characterization of bioactive metabolites from piezotolerant deep sea fungus Nigrospora sp. In submerged fermentation.J. Appl. Microbiol.11899–111. 10.1111/jam.12693
5
BarashN. R.MaloneyJ. G.SingerS. M.DawsonS. C. (2017). Giardia Alters Commensal Microbial Diversity throughout the Murine Gut.Infect. Immun.85e00948–16. 10.1128/IAI.00948-16
6
BatraN.KaurH.MohindraS.SinghS.ShamanthA. S.RudramurthyS. M. (2019). Cladosporium sphaerospermum causing brain abscess, a saprophyte turning pathogen: case and review of published reports.J. Mycol. Med.29180–184. 10.1016/j.mycmed.2019.04.005
7
BerrilliF.Di CaveD.CavalleroS.D’AmelioS. (2012). Interactions between parasites and microbial communities in the human gut.Front. Cell Infect. Microbiol.2:141. 10.3389/fcimb.2012.00141
8
BolickD. T.Mayneris-PerxachsJ.MedlockG. L.KollingG. L.PapinJ. A.SwannJ. R.et al (2017). Increased urinary trimethylamine N-Oxide following cryptosporidium infection and protein malnutrition independent of microbiome effects.J. Infect. Dis.21664–71. 10.1093/infdis/jix234
9
BouzidM.HunterP. R.ChalmersR. M.TylerK. M. (2013). Cryptosporidium pathogenicity and virulence.Clin. Microbiol. Rev.26115–134. 10.1128/CMR.00076-12
10
ChanJ. F.LauS. K.YuenK. Y.WooP. C. (2016). Talaromyces (Penicillium) marneffei infection in non-HIV-infected patients.Emerg. Microbes. Infect.5:e19. 10.1038/emi.2016.18
11
ChappellC. L.DarkohC.ShimminL.FarhanaN.KimD. K.OkhuysenP. C.et al (2016). Fecal indole as a biomarker of susceptibility to cryptosporidium infection.Infect. Immun.842299–2306. 10.1128/IAI.00336-16
12
CharaniaR.WadeB. E.McNairN. N.MeadJ. R. (2020). Changes in the microbiome of Cryptosporidium-Infected mice correlate to differences in susceptibility and infection levels.Microorganisms8:879. 10.3390/microorganisms8060879
13
ChenS.ZhaoH.ZouC.LiY.ChenY.WangZ.et al (2017). Combined inoculation with multiple arbuscular mycorrhizal fungi improves growth, nutrient uptake and photosynthesis in cucumber seedlings.Front. Microbiol.8:2516. 10.3389/fmicb.2017.02516
14
ChinV. K.YongV. C.ChongP. P.AminN. S.BasirR.AbdullahM. (2020). Mycobiome in the gut: a multiperspective review.Mediators Inflamm.2020:9560684. 10.1155/2020/9560684
15
ColemanS. U.KleiT. R.FrenchD. D.ChapmanM. R.CorstvetR. E. (1989). Prevalence of Cryptosporidium sp in equids in Louisiana.Am. J. Vet. Res.50575–577.
16
CooperP.WalkerA. W.ReyesJ.ChicoM.SalterS. J.VacaM.et al (2013). Patent human infections with the whipworm, Trichuris trichiura, are not associated with alterations in the faecal microbiota.PLoS One8:e76573. 10.1371/journal.pone.0076573
17
CurrentW. L.GarciaL. S. (1991). Cryptosporidiosis.Clin. Microbiol. Rev.4325–358. 10.1128/CMR.4.3.325
18
DengL.LiW.ZhongZ.GongC.CaoX.SongY.et al (2017). Occurrence and Genetic Characteristics of Cryptosporidium hominis and Cryptosporidium andersoni in Horses from Southwestern China.J. Eukaryot. Microbiol.64716–720. 10.1111/jeu.12399
19
FanW.HuoG.LiX.YangL.DuanC. (2014). Impact of diet in shaping gut microbiota revealed by a comparative study in infants during the six months of life.J. Microbiol. Biotechnol.24133–143. 10.4014/jmb.1309.09029
20
FayerR.SantinM.MacarisinD.BauchanG. (2013). Adhesive-tape recovery combined with molecular and microscopic testing for the detection of Cryptosporidium oocysts on experimentally contaminated fresh produce and a food preparation surface.Parasitol. Res.1121567–1574. 10.1007/s00436-013-3305-7
21
GilchristC. A.PetriS. E.SchneiderB. N.ReichmanD. J.JiangN.BegumS.et al (2016). Role of the gut microbiota of children in diarrhea due to the protozoan parasite entamoeba histolytica.J. Infect. Dis.2131579–1585. 10.1093/infdis/jiv772
22
GraczykT. K.FayerR.CranfieldM. R. (1997). Zoonotic transmission of Cryptosporidium parvum: implications for water-borne cryptosporidiosis.Parasitol. Today13348–351. 10.1016/s0169-4758(97)01076-4
23
GruningerR. J.PuniyaA. K.CallaghanT. M.EdwardsJ. E.YoussefN.DagarS. S.et al (2014). Anaerobic fungi (phylum Neocallimastigomycota): advances in understanding their taxonomy, life cycle, ecology, role and biotechnological potential.FEMS Microbiol. Ecol.901–17. 10.1111/1574-6941.12383
24
HaridasS.AlbertR.BinderM.BloemJ.LaButtiK.SalamovA.et al (2020). 101 Dothideomycetes genomes: a test case for predicting lifestyles and emergence of pathogens.Stud. Mycol.96141–153. 10.1016/j.simyco.2020.01.003
25
HayesK. S.BancroftA. J.GoldrickM.PortsmouthC.RobertsI. S.GrencisR. K. (2010). Exploitation of the intestinal microflora by the parasitic nematode Trichuris muris.Science3281391–1394. 10.1126/science.1187703
26
HooperL. V.LittmanD. R.MacphersonA. J. (2012). Interactions between the microbiota and the immune system.Science3361268–1273. 10.1126/science.1223490
27
Human Microbiome Project Consortium [HMPC] (2012). Structure, function and diversity of the healthy human microbiome.Nature486207–214. 10.1038/nature11234
28
HunterP. R.NicholsG. (2002). Epidemiology and clinical features of Cryptosporidium infection in immunocompromised patients.Clin. Microbiol. Rev.15145–154. 10.1128/CMR.15.1.145-154.2002
29
HuseyinC. E.O’TooleP. W.CotterP. D.ScanlanP. D. (2017). Forgotten fungi-the gut mycobiome in human health and disease.FEMS Microbiol. Rev.41479–511. 10.1093/femsre/fuw047
30
IlievI. D.FunariV. A.TaylorK. D.NguyenQ.ReyesC. N.StromS. P.et al (2012). Interactions between commensal fungi and the C-type lectin receptor Dectin-1 influence colitis.Science3361314–1317. 10.1126/science.1221789
31
JianF.LiuA.WangR.ZhangS.QiM.ZhaoW.et al (2016). Common occurrence of Cryptosporidium hominis in horses and donkeys.Infect. Genet. Evol.43261–266. 10.1016/j.meegid.2016.06.004
32
JiangJ.AlderisioK. A.XiaoL. (2005). Distribution of cryptosporidium genotypes in storm event water samples from three watersheds in New York.Appl. Environ. Microbiol.714446–4454. 10.1128/AEM.71.8.4446-4454.2005
33
KarpeA. V.HuttonM. L.MiletoS. J.JamesM. L.EvansC.ShahR. M.et al (2021). Cryptosporidiosis modulates the gut microbiome and metabolism in a murine infection model.Metabolites11:380. 10.3390/metabo11060380
34
LantierL.DrouetF.GuesdonW.MancassolaR.MettonC.Lo-ManR.et al (2014). Poly(I:C)-induced protection of neonatal mice against intestinal Cryptosporidium parvum infection requires an additional TLR5 signal provided by the gut flora.J. Infect. Dis.209457–467. 10.1093/infdis/jit432
35
LeeS. C.TangM. S.LimY.ChoyS. H.KurtzZ. D.CoxL. M.et al (2021). Correction: helminth colonization is associated with increased diversity of the gut microbiota.PLoS Negl. Trop. Dis.15:e9325. 10.1371/journal.pntd.0009325
36
LeungH. T.MaasK. R.WilhelmR. C.MohnW. W. (2016). Long-term effects of timber harvesting on hemicellulolytic microbial populations in coniferous forest soils.ISME J.10363–375. 10.1038/ismej.2015.118
37
LiF.SuJ.ChahanB.GuoQ.WangT.YuZ.et al (2019). Different distribution of Cryptosporidium species between horses and donkeys.Infect. Genet. Evol.75:103954. 10.1016/j.meegid.2019.103954
38
MaJ.LiP.ZhaoX.XuH.WuW.WangY.et al (2015). Occurrence and molecular characterization of Cryptosporidium spp. And Enterocytozoon bieneusi in dairy cattle, beef cattle and water buffaloes in China.Vet. Parasitol.207220–227. 10.1016/j.vetpar.2014.10.011
39
MaL.ZhangX.JianY.LiX.WangG.HuY.et al (2019). Detection of Cryptosporidium and Giardia in the slaughterhouse, sewage and river waters of the Qinghai Tibetan plateau area (QTPA), China.Parasitol. Res.1182041–2051. 10.1007/s00436-019-06330-w
40
MairT. S.TaylorF. G.HarbourD. A.PearsonG. R. (1990). Concurrent cryptosporidium and coronavirus infections in an Arabian foal with combined immunodeficiency syndrome.Vet. Rec.126127–130.
41
MartinezI.StegenJ. C.Maldonado-GomezM. X.ErenA. M.SibaP. M.GreenhillA. R.et al (2015). The gut microbiota of rural papua new guineans: composition, diversity patterns, and ecological processes.Cell Rep.11527–538. 10.1016/j.celrep.2015.03.049
42
MascarinG. M.JaronskiS. T. (2016). The production and uses of Beauveria bassiana as a microbial insecticide.World J. Microbiol. Biotechnol.32:177. 10.1007/s11274-016-2131-3
43
McKenneyE. A.GreeneL. K.DreaC. M.YoderA. D. (2017). Down for the count: cryptosporidium infection depletes the gut microbiome in Coquerel’s sifakas.Microb. Ecol. Health Dis.28:1335165. 10.1080/16512235.2017.1335165
44
McKenneyE. A.WilliamsonL.YoderA. D.RawlsJ. F.BilboS. D.ParkerW. (2015). Alteration of the rat cecal microbiome during colonization with the helminth Hymenolepis diminuta.Gut. Microbes.6182–193. 10.1080/19490976.2015.1047128
45
MishraP. K.PalmaM.BleichD.LokeP.GauseW. C. (2014). Systemic impact of intestinal helminth infections.Mucosal. Immunol.7753–762. 10.1038/mi.2014.23
46
NimeF. A.BurekJ. D.PageD. L.HolscherM. A.YardleyJ. H. (1976). Acute enterocolitis in a human being infected with the protozoan Cryptosporidium.Gastroenterology70592–598.
47
NundaengS.SuwannarachN.LimtongS.KhunaS.KumlaJ.LumyongS. (2021). An updated global species diversity and phylogeny in the genus wickerhamomyces with addition of two new species from thailand.J Fungi7:957. 10.3390/jof7110957
48
OttS. J.KühbacherT.MusfeldtM.RosenstielP.HellmigS.RehmanA.et al (2008). Fungi and inflammatory bowel diseases: alterations of composition and diversity.Scand J Gastroenterol.43831–841. 10.1080/00365520801935434
49
ParkinsonL. E.LeD. P.DannE. K. (2019). Development of three Loop-Mediated isothermal amplification (LAMP) assays for the rapid detection of calonectria ilicicola, dactylonectria macrodidyma, and the dactylonectria genus in avocado roots.Plant Dis.1031865–1875. 10.1094/PDIS-11-18-2005-RE
50
PedersenA. B.FentonA. (2007). Emphasizing the ecology in parasite community ecology.Trends Ecol. Evol.22133–139. 10.1016/j.tree.2006.11.005
51
PolviE. J.LiX.O’MearaT. R.LeachM. D.CowenL. E. (2015). Opportunistic yeast pathogens: reservoirs, virulence mechanisms, and therapeutic strategies.Cell. Mol. Life Sci.722261–2287. 10.1007/s00018-015-1860-z
52
ReynoldsL. A.FinlayB. B.MaizelsR. M. (2015). Cohabitation in the intestine: interactions among helminth parasites, bacterial microbiota, and host immunity.J. Immunol.1954059–4066. 10.4049/jimmunol.1501432
53
RigaudT.Perrot-MinnotM. J.BrownM. J. (2010). Parasite and host assemblages: embracing the reality will improve our knowledge of parasite transmission and virulence.Proc. Biol. Sci.2773693–3702. 10.1098/rspb.2010.1163
54
RyanU.XiaoL.ReadC.ZhouL.LalA. A.PavlasekI. (2003). Identification of novel Cryptosporidium genotypes from the Czech Republic.Appl. Environ. Microbiol.69, 4302–4307. 10.1128/AEM.69.7.4302-4307.2003
55
SenderR.FuchsS.MiloR. (2016). Revised estimates for the number of human and bacteria cells in the body.PLoS Biol.14:e1002533. 10.1371/journal.pbio.1002533
56
SowS. O.MuhsenK.NasrinD.BlackwelderW. C.WuY.FaragT. H.et al (2016). The burden of cryptosporidium diarrheal disease among children < 24 months of age in Moderate/High mortality regions of Sub-Saharan africa and south asia, utilizing data from the global enteric multicenter study (GEMS).PLoS Negl. Trop. Dis.10:e4729. 10.1371/journal.pntd.0004729
57
SturdeeA. P.Bodley-TickellA. T.ArcherA.ChalmersR. M. (2003). Long-term study of Cryptosporidium prevalence on a lowland farm in the United Kingdom.Vet. Parasitol.11697–113. 10.1016/s0304-4017(03)00261-9
58
TonjerL. R.ThoenE.MorgadoL.BotnenS.MundraS.NybakkenL.et al (2021). Fungal community dynamics across a forest-alpine ecotone.Mol. Ecol.304926–4938. 10.1111/mec.16095
59
TziporiS.SmithM.BirchC.BarnesG.BishopR. (1983). Cryptosporidiosis in hospital patients with gastroenteritis.Am. J. Trop. Med. Hyg.32931–934. 10.4269/ajtmh.1983.32.931
60
VillarinoN. F.LeCleirG. R.DennyJ. E.DearthS. P.HardingC. L.SloanS. S.et al (2016). Composition of the gut microbiota modulates the severity of malaria.Proc Natl Acad Sci U S A.1132235–2240. 10.1073/pnas.1504887113
61
WangY.FuY.HeY.KulyarM. F.IqbalM.LiK.et al (2021a). Longitudinal Characterization of the Gut Bacterial and Fungal Communities in Yaks.J Fungi7:559. 10.3390/jof7070559
62
WangM.LuoX.WangH. (2021b). First Draft Genome Sequence Resource of Colletotrichum liriopes Causing Leaf Anthracnose on Ophiopogon japonicus.Plant Dis.1051179–1182. 10.1094/PDIS-11-20-2326-A
63
YatsunenkoT.ReyF. E.ManaryM. J.TrehanI.Dominguez-BelloM. G.ContrerasM.et al (2012). Human gut microbiome viewed across age and geography.Nature486222–227. 10.1038/nature11053
Summary
Keywords
Cryptosporidium infection, prevalence, fungi microbiota, horse, protozoa
Citation
Wang Y, Li X, Chen X, Kulyar MF-A, Duan K, Li H, Bhutta ZA, Wu Y and Li K (2022) Gut Fungal Microbiome Responses to Natural Cryptosporidium Infection in Horses. Front. Microbiol. 13:877280. doi: 10.3389/fmicb.2022.877280
Received
22 February 2022
Accepted
18 May 2022
Published
06 July 2022
Volume
13 - 2022
Edited by
Joyce Siwila, University of Zambia, Zambia
Reviewed by
Guanghui Zhao, Northwest A&F University, China; Prerona Boruah, DY Patil Deemed to be University, India
Updates
Copyright
© 2022 Wang, Li, Chen, Kulyar, Duan, Li, Bhutta, Wu and Li.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kun Li, lk3005@njau.edu.cn
This article was submitted to Microorganisms in Vertebrate Digestive Systems, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.