- Open Access
Molecular survey of Enterocytozoon bieneusi in sheep and goats in China
Parasites & Vectorsvolume 9, Article number: 23 (2016)
Over recent years, several studies have conducted genotyping of Enterocytozoon bieneusi in various hosts worldwide using sequence analysis of the ribosomal internal transcribed spacer (ITS), however, relatively little is known about E. bieneusi in sheep and goats in China. Therefore this research was conducted to understand the prevalence and genotype distribution of E. bieneusi in farmed sheep and goats in China.
A total of 1025 fecal specimens from farmed animals in various geographic areas were collected. Overall, PCR and sequence analysis of the ITS detected E. bieneusi in 34.4 % (353/1025) of isolates; of which the prevalence in goats was 28.8 % (176/611) and in sheep was 42.8 % (177/414). Phylogenetic analysis of ITS sequences identified 42 genotypes (nine known and 33 novel ones). These consisted of four known genotypes (D, KIN-1, EbpC, and F) and 10 novel genotypes (CHG6, CHG7, CHG9, CHG19, CHG23, CHG25, CHS5 and CHS10–CHS12) which all belonged to the so-called zoonotic group 1. A further four known genotypes (CD6, COS-I, BEB6, and J) and 22 novel genotypes (CHG1–CHG3, CHG5, CHG8, CHG10–CHG14, CHG16–CHG18, CHG20, CHG22, CHG24, CHS3, CHS4 and CHS6–CHS9) formed a clade within the group 2. One novel genotype (CHG21) was clustered in the group 9 with the genotype CM4.
E. bieneusi is highly prevalent, widely distributed, and genetically diverse in Chinese farmed goats and sheep. Some of the genotypes identified are potentially zoonotic.
Microsporidia are a diverse group of obligate intracellular pathogens related to fungi . To date, approximately 1300 species belonging to 160 microsporidian genera are known to infect invertebrate and vertebrate hosts, including humans [2, 3]. Among the 14 microsporidian species infectious to humans, Enterocytozoon bieneusi is regarded as the most commonly identified microsporidian species in humans [4, 5]. Thus far, in addition to humans, E. bieneusi has been detected in domestic animals such as cattle, sheep, goats, rabbits, pigs, dogs, cats, and horses [6–14], and in wild animals such as raccoons, muskrats, beavers, foxes, otters, nonhuman primates (NHPs), llamas, guinea pigs, deer and snakes [15–21] and birds .
E. bieneusi spores are morphologically indistinguishable and in vitro cultivation methods for this pathogen have not been developed. Therefore, molecular methods are required to identify genotypes within the species that infect humans, animals or both . Until now, the internal transcribed spacer (ITS) of the rRNA gene has been used widely for species identification and genotyping of E. bieneusi isolates [11, 23, 24] and more than 200 E. bieneusi ITS genotypes have been identified [3–12]. Previous studies have reported that some genotypes are host-specific, whereas others are zoonotic because of their capability to infect both humans and animals [15, 16, 25–28]. Currently, over 30 countries have been surveyed to gather data on the molecular epidemiology of E. bieneusi in humans and other animals . In China, this pathogen has been detected in HIV-positive and -negative patients, children, non-human primates (NHPs), dogs, cats, pigs, cows, sheep, goats, deer and snakes [8, 9, 11, 12, 20, 21, 30–37]. Moreover, this parasite has been identified in wastewater in several Chinese cities . However, few studies have been conducted on E. bieneusi infections in sheep and goats in China . In the present study, a large number of farms were surveyed to gain a better understanding of the prevalence and genotype distribution of E. bieneusi in sheep and goats in various parts of China.
This research was conducted according to the Chinese Laboratory Animal Administration Act (1988) after review and approval of its protocol by the Research Ethics Committee of Henan Agricultural University. Appropriate permission was obtained from the farm owners before collection of fecal specimens from their sheep and goats.
Fecal specimen collection
A total of 1025 fresh fecal samples were collected from sheep and goats in Henan, Heilongjiang, Liaoning, Yunnan, Anhui, and Shaanxi Provinces and Chongqing city between March 2011 and October 2013 (Table 1 and Additional file 1). Specimens were collected either directly from the rectum or were collected in sterile plastic containers during defecation (without the fecal sample having contact with the ground). One fecal specimen per sheep or goat was collected in a container marked with the breed, gender, age, and feeding habits of the animal. No obvious clinical symptoms were observed during specimen collection, which were then transported in an icebox to the laboratory for examination.
The fecal specimen from each container was transferred to a sterile beaker immediately on arrival at the laboratory, soaked for thirty minutes at room temperature in 30 ml of tap water, and then stirred with a clean stick to ensure thorough mixing. Each mixture was then put through a 7.62-cm-diameter sieve (pore size 45 μm) and concentrated by centrifugation at 3500 r · min-1 for 5 min. The concentrated fecal material was stored in 2.5 % potassium dichromate solution at 4 °C prior to DNA extraction.
DNA extraction and PCR amplification
The stored fecal specimens were washed three times with distilled water and genomic DNA was extracted by using an E.Z.N.A.® Stool DNA kit (Omega Biotek Inc., Norcross, GA, USA) according to the procedure recommended by the manufacturer. The extracted DNA was stored at −20 °C before it was used for PCR amplification. DNA from each specimen was tested at least three times for E. bieneusi DNA using a nested PCR assay that amplified a 392 bp product that encompassed the terminal part of 18S rRNA gene, the complete internal transcribed spacer (ITS), and a partial region of 5.8S rRNA gene . The KOD-Plus-Neo amplification enzyme (Toyobo Co. Ltd., Osaka, Japan) was used for PCR amplification. Nonacetylated bovine serum albumin (Solarbio Co. Ltd., Beijing, China) (final concentration 400 ng/L) was used to neutralize PCR inhibitors during the primary PCR. The secondary PCR products were examined by agarose gel electrophoresis and visualized after GelRed (Biotium Inc., Hayward, CA) staining.
Sequencing and phylogenetic analysis
Since some fecal samples were collected from the same animals, positive specimens from different farms, animals, species and genders were selected for sequencing. The selected secondary PCR products were purified using Montage PCR filters (Millipore, Bedford, MA) and sequenced using a BigDye Terminator v3.1 cycle sequencing kit (Applied Biosystems) on an ABI 3730 DNA analyzer (Applied Biosystems, Foster City, CA). The nucleotide sequences were confirmed by bidirectional sequencing and by sequencing a new PCR product if necessary. The sequences obtained were aligned with reference sequences downloaded from GenBank using the ClustalX 2.0 (http://www.clustal.org/) program to determine the genotypes. Genotypes from this study were compared with known E. bieneusi ITS genotypes by neighbor-joining analysis of the aligned E. bieneusi sequences using Mega 5.0 (http://www.megasoftware.net/). Bootstrap analysis was used to assess the robustness of clusters using 1000 replicates. A previously established nomenclature system was used for naming the E. bieneusi ITS genotypes .
Infection rates were compared using the chi-square test and a difference was considered statistically significant when the P value was <0.05. The analysis was done using QuickCalcs software (GraphPad Software Inc., La Jolla, CA).
Nucleotide sequence accession numbers
Representative nucleotide sequences from this study have been deposited in GenBank under accession numbers KP262357–KP262398.
Prevalence of E. bieneusi in goats and sheep
Of the 1025 fecal specimens from goats and sheep, overall 353 specimens (34.4 %) were E. bieneusi-positive by PCR amplification of the ITS locus. The prevalence in goats was 28.8 % (176/611) and in sheep was 42.8 % (177/414). The differences in the infection rates between sheep and goats were statistically significant different (χ2 = 21.27, p < 0.01). For goats, E. bieneusi was found on 17 of the studied farms from four provinces and one city (Table 1 and Additional file 1). The prevalence ranged from 7.5 % to 62.5 % for each goat farm, with significant differences in prevalence between different goat farms(χ2 = 73.21, p < 0.01). The prevalence rates were 62.5 % in goats from Chongqing City, 47.8 % in Shaanxi Province, 32.9 % in Henan Province, 22.4 % in Yunnan Province, and 7.5 % in goats from Anhui Province (χ2 = 35.81, p < 0.01). For sheep, E. bieneusi was detected in six farms from three provinces. The prevalence varied from 9.4 % to 67.4 % in the studied sheep farms, and the differences in infection rates were statistically significant (χ2 = 67.84, p < 0.01) (Table 1 and Additional file 1). The prevalence were 51.9 %, 25.0 %, and 9.4 % in sheep from Henan Province, Heilongjiang Province, and Liaoning Province, respectively (χ2 = 44.96, p < 0.01) (Table 1).
There were no statistically significant differences in the prevalence among the age groups in goats (χ2 = 4.50, p > 0.05) and sheep (χ2 = 5.61, p > 0.05) (Table 3). Likewise, there were no statistically significant differences for gender in goats (χ2 = 0.01, p > 0.05) and sheep (χ2 = 2.05, p > 0.05) (Table 2). There were some differences in prevalence between grazing and household goats, with a higher prevalence in household goats (χ2 = 5.75, 0.01 < p < 0.05) (Table 2), but there was no significant difference in prevalence for grazing and household sheep (χ2 = 0.65, p > 0.05) (Table 2).
Distribution of E. bieneusi genotypes in goats and sheep
A total of 42 ITS genotypes were obtained from 230 successfully sequenced specimens from goats and sheep. Of them, nine genotypes have been previously reported, including BEB6 (EU153584), D (AF101200, syn. CEbC; PigEBITS9; VI; WL18; Peru 9; PtEb), EbpC (AF135832, syn. E; WL13; Peru4; WL17), F (AF135833, syn. EbpA), KIN-1(JQ437573), J(AF135837), COS-I(KJ850432, syn. CM7), CM4(KF543866) and CD6(KJ668733), while other 33 genotypes were novel (CHG1–CHG3, CHG5–CHG14, CHG16–CHG25 and CHS3–CHS12). A higher sequence polymorphism was identified in goat isolates compared to sheep specimens. Of the 139 goat isolates that were successfully sequenced, 31 genotypes were identified, eight of which were known (BEB6, D, EbpC, F, KIN-1, J, COS-I and CD6), while other 23 were new genotypes (CHG1–CHG3, CHG5–CHG14, CHG16–CHG25). For sheep, 91 specimens were sequenced successfully and 14 genotypes were identified, including three known genotypes (BEB6, CM4 and COS-I) and 11 novel genotypes (CHG3 and CHS3–CHS12). Genotypes BEB6, COS-I and CHG3 were found in both sheep and goat specimens (Tables 1, 2, 3).
The genotype BEB6 was the dominant genotype in goat (41/139) and sheep isolates (60/91) in this study. In goats, in addition to BEB6, the other dominant genotypes, namely CHG1 and CHG3, were detected in 19 and 17 specimens, respectively. The genotypes CD6, CHG5, EbpC, CHG2, F, and D were identified in nine, seven, eight, six, four, and three specimens, respectively. Genotypes KIN-1, COS-I and CHG16 were observed in two specimens each, while all other genotypes we identified were found in one specimen each. In sheep, the genotype COS-I was the second dominant genotype, being identified in 14 specimens. Genotypes CHG3 and CHS2 were identified in five and two specimens, respectively, while all other genotypes we obtained were identified in one specimen each. The genotype distribution according to the age, gender and the feeding habits of the animals is shown in Tables 2 and 3.
Phylogenetic analysis revealed that some novel genotypes from this study belonged to the human-pathogenic group , recently renamed group 1 [39, 40]. We identified genotypes D, KIN-1, EbpC and F reported previously and 10 novel genotypes (CHG6, CHG7, CHG9, CHG19, CHG23, CHG25, CHS5, CHS10, CHS11 and CHS12) belonged to this group. The novel genotypes CHG23, CHG25, CHS10 and CHG19 had one (C to A substitution at nucleotide position 362), two (T to C and T to A substitutions at positions 243 and 245), two (T to G substitutions at positions 223 and 252) and three (A to G, T to G and G to A substitutions at positions 137, 201 and 374) single nucleotide polymorphisms (SNPs) compared with the established genotype, EbpC (AF135832). Therefore, genotypes CHG23, CHG25, NESH1, CHS10 and CHG19 were clustered into subgroup 1d with the EbpC genotype. The following new genotypes CHS5, CHS11 and CHS12 were grouped with subgroup 1e with known genotypes COG-1 (KJ850434) and F (AF135833). Compared with the known genotype F (AF135833), one or several SNPs were present (data not shown). The new genotype CHG9 differed from the genotype C (AF101199) at nucleotide positions 123 (G to A substitution), 230 (G to A substitution), 255(G to T substitution) and 284 (G to A substitution), respectively, thus they formed subgroup 1 f. However, the new genotypes CHG6 and CHG7 formed a cluster distinct from any known subgroup 1. We found that 22 new genotypes (CHG1–CHG3, CHG5, CHG8, CHG10–CHG14, CHG16–CHG18, CHG20, CHG22, CHG24, CHS3, CHS4, and CHS6–CHS9) and 7 known genotypes (CD6, COS-I, BEB6, J, COS-III, NESH4 and CHN2) belonged to the host-specific group 2, but the known genotypes J (AF135837) and CHN2 (HM992510) formed one branch, and the 22 new genotypes with the reported genotypes CD6(KJ668733), COS-I(KJ850432),COS-III(KP063053), NESH4(KP732478) and BEB6(EU153584) formed another clade. Compared with genotype BEB6 (EU153584), new genotypes CHG2, CHG13, CHG16, CHG17, CHS3, CHS4, CHS7, CHS8 and CHS9 contained one SNP, genotypes CHG3, CHG5, CHG12, CHG22 and CHS6 contained two SNPs, genotypes CHG1, CHG14, CHG18 and CHG20 contained three SNPs, and genotypes CHG11, CHG10 and CHG24 contained 4, 11 and 13 SNPs, respectively, the genotype CHG8 contained over 15 SNPs (Additional file 2). However, the new genotype CHG21 was clustered with the CM4 (KF543866) genotype [20, this study] in the group 9  ( Additional file 3: Figure S1).
In the present study, E. bieneusi was found to be a common parasite in sheep and goats. The prevalence of this organism in healthy sheep and goats was 42.8 % and 28.8 %, respectively. A similar prevalence (21.8 %) was reported for E. bieneusi in goats in Heilongjiang Province , while lower prevalence rates (22.5 % and 13.9 %) were detected in Heilongjiang Province  and northeast China  in sheep. A higher E. bieneusi prevalence (69.3 %) was found in sheep in Inner Mongolia . The differences in the prevalence between this study and others [6, 9, 34, 35, 37] are likely due to differences in the age groups studied, as well as the number of study animals, detection methods, and the geographic and ecological environments. Here we described the infection rates of E. bieneusi in different age groups in sheep and goats and showed that there were no significant differences among different age groups, which was consistent with the conclusion obtained by Jiang and others  but inconsistent with the results gained by Stensvold . In this study, we have described the prevalence of E. bieneusi by the gender and different feeding habits of the host animals for the first time in detail, although there were no significant differences between them. Nevertheless, these results provide an useful reference data for researchers in this field.
All nine known genotypes identified in the present study have been previously reported in China: the genotype BEB6 has been reported in one child (reported as SH5, ), sheep , NHPs [20, 42], deer , cats  and golden takins , even in urban wastewater . This genotype was first detected in cattle in the United States  and is considered to be a cattle-adopted genotype. However, the genotype BEB6 was detected in 43.9 % (101/230) sequenced specimens and found in almost every farm (Additional file 1), which showed that this genotype was observed as a dominant genotype with a wide geographic distribution, extensive host range and zoonotic potential both in the present study and in other reports [9, 10, 23, 34–37]. The genotype D has been found in AIDS patients , children , NHPs [33, 42], pigs , dogs [12, 44], cats [12, 44], foxes , raccoon dogs , cattle , golden takins , urban wastewater  and drinking source watershed . In this report, this genotype was only found in goats, in contrast to the results reported in Heilongjiang Province . Though the genotype D only infected goats in this study, the transmission source and transmission dynamics of this genotype is still unclear and further research is required to understand the role of this genotype in microsporidiosis. The genotype EbpC has been found in HIV-positive and -negative patients , children , NHPs [20, 33], dogs [12, 44], pigs [8, 9], lake water , red pandas  and urban wastewater . Here, the genotype EbpC was detected only in goat specimens, which was similar to the results of Zhao et al , but was in contrast with Li et al . The genotype F has been found in children , pigs , cattle  and wastewater . Genotypes J has been reported in humans , cattle  and yaks , KIN-1 in humans (JQ437573, unpublished), CD6(KJ668733), COS-I(KJ850432) (syn. CM7) in sheep [35, 37], deer  and NHPs , CM4 in NHPs  and CD6 in dogs .
Among the nine established genotypes found in this study, six (BEB6, D, EbpC, F, KIN-1, J) genotypes have been reported in humans [11, 30, 31]. Thus, some of the genotypes in sheep and goats have zoonotic potential, and sheep and goats may be reservoirs of human microsporidiosis.
Phylogenetic analysis showed that 10 new E. bieneusi genotypes (CHG6, CHG7, CHG9, CHG19, CHG23, CHG25, CHS5, CHS10, CHS11, and CHS12) belonged to the so-called group 1 with zoonotic potential . This finding merits close attention because of the ability of these genotypes to cause E. bieneusi-related microsporidiosis in humans. The other 22 new genotypes (CHG1–CHG3, CHG5, CHG8, CHG10–CHG14, CHG16–CHG18, CHG20, CHG22, CHG24, CHS3, CHS4, and CHS6–CHS9) were classified into the cattle-specific group 2 . Some genotypes from this group initially considered to be cattle-specific (such as BEB6 and J), were later found in other animals including humans [11, 12, 20, 30, 50, this research], therefore some genotypes in the group 2 may have limited zoonotic potential and should be examined in more detail. The remaining new genotype CHG21 with only one SNP compared with the genotype CM4 formed a cluster, in the previously described group 9 .
This study to infer the prevalence and genotype distribution of E. bieneusi in sheep and goats in China. The data showed there were a widespread distribution, high prevalence, and considerable genetic diversity in E. bieneusi genotypes from farmed goats and sheep, which may cause microsporidiosis among sheep, goats, other animals and humans resulting in economic losses in animal husbandry. Moreover, many genotypes from this study have been reported in humans and urban wastewater in China, which may result in zoonotic transmission. Therefore, contact between susceptible human populations and goats or sheep should be avoided to reduce zoonotic or anthroponotic parasite transmission. Further studies are needed to fully elucidate the importance of goats and sheep in the epidemiology of human microsporidiosis.
Keeling PJ, Fast NM. Microsporidia: biology and evolution of highly reduced intracellular parasites. Annu Rev Microbiol. 2002;56:93–116.
Keeling P. Five questions about microsporidia. PLoS Pathog. 2009;5:e1000489.
Santín M, Fayer R. Microsporidiosis: Enterocytozoon bieneusi in domesticated and wild animals. Res Vet Sci. 2011;90:363–71.
Desportes I, Le Charpentier Y, Galian A, Bernard F, Cochand-Priollet B, Lavergne A, et al. Occurrence of a new microsporidan: Enterocytozoon bieneusi n. g., n. sp., in the enterocytes of a human patient with AIDS. J Protozool. 1985;32:250–4.
Didier ES. Microsporidiosis: an emerging and opportunistic infection in humans and animals. Acta Trop. 2005;94:61–76.
Lores B, del Aguila C, Arias C. Enterocytozoon bieneusi (Microsporidia) in faecal specimen from domestic animals from Galicia, Spain. Mem Inst Oswaldo Cruz. 2002;97:941–5.
Del Aguila C, Izquierdo F, Navajas R, Pieniazek NJ, Miró G, Alonso AI, et al. Enterocytozoon bieneusi in animals: rabbits and dogs as new hosts. J Eukaryot Microbiol. 1999;46:8S–9S.
Li W, Diao R, Yang J, Xiao L, Lu Y, Li Y, et al. High diversity of human-pathogenic Enterocytozoon bieneusi genotypes in swine in northeast China. Parasitol Res. 2014;113:1147–53.
Li W, Li Y, Li W, Yang J, Song M, Diao R, et al. Genotypes of Enterocytozoon bieneusi in livestock in China: high prevalence and zoonotic potential. PLoS One. 2014;9:e97623.
Stensvold CR, Beser J, Ljungström B, Troell K, Lebbad M. Low host-specific Enterocytozoon bieneusi genotype BEB6 is common in Swedish lambs. Vet Parasitol. 2014;205:371–4.
Zhang X, Wang Z, Su Y, Liang X, Sun X, Peng S, et al. Identification and genotyping of Enterocytozoon bieneusi in China. J Clin Microbiol. 2011;49:2006–8.
Karim MR, Dong H, Yu F, Jian F, Zhang L, Wang R, et al. Genetic diversity in Enterocytozoon bieneusi isolates from dogs and cats in China: host specificity and public health implications. J Clin Microbiol. 2014;52:3297–302.
Fayer R, Santin M, Trout JM. Enterocytozoon bieneusi in mature dairy cattle on farms in the eastern United States. Parasitol Res. 2007;102:15–20.
Wagnerová P, Sak B, Květoňová D, Buňatová Z, Civišová H, Maršálek M, et al. Enterocytozoon bieneusi and Encephalitozoon cuniculi in horses kept under different management systems in the Czech Republic. Vet Parasitol. 2012;190:573–7.
Sulaiman IM, Fayer R, Lal AA, Trout JM, Schaefer 3rd FW, Xiao L. Molecular characterization of microsporidia indicates that wild mammals harbor host-adapted Enterocytozoon spp. as well as human-pathogenic Enterocytozoon bieneusi. Appl Environ Microbiol. 2003;69:4495–501.
Dengjel B, Zahler M, Hermanns W, Heinritzi K, Spillmann T, Thomschke A, et al. Zoonotic potential of Enterocytozoon bieneusi. J Clin Microbiol. 2001;39:4495–9.
Drosten C, Laabs J, Kuhn EM, Schottelius J. Interspecies transmission of Enterozytozoon bieneusi supported by observations in laboratory animals and phylogeny. Med Microbiol Immunol. 2005;194:207–9.
Cama VA, Pearson J, Cabrera L, Pacheco L, Gilman R, Meyer S, et al. Transmission of Enterocytozoon bieneusi between a child and guinea pigs. J Clin Microbiol. 2007;45:2708–10.
Santín M, Fayer R. Enterocytozoon bieneusi, Giardia, and Cryptosporidium infecting white-tailed deer. J Eukaryot Microbiol. 2014;194:207–9.
Karim MR, Wang R, Dong H, Zhang L, Li J, Zhang S, et al. Genetic polymorphism and zoonotic potential of Enterocytozoon bieneusi from nonhuman primates in China. Appl Environ Microbiol. 2014;80:1893–8.
Karim MR, Yu F, Li J, Li J, Zhang L, Wang R, et al. First molecular characterization of enteric protozoa and the human pathogenic microsporidian, Enterocytozoon bieneusi, in captive snakes in China. Parasitol Res. 2014;113:3041–8.
Haro M, Izquierdo F, Henriques-Gil N, Andrés I, Alonso F, Fenoy S, et al. First detection and genotyping of human-associated microsporidia in pigeons from urban parks. Appl Environ Microbiol. 2005;71:3153–7.
Feng Y, Li N, Dearen T, Lobo ML, Matos O, Cama V, et al. Development of a multilocus sequence typing tool for high-resolution genotyping of Enterocytozoon bieneusi. Appl Environ Microbiol. 2011;77:4822–8.
Santín M, Fayer R. Enterocytozoon bieneusi genotype nomenclature based on the internal transcribed spacer sequence: a consensus. J Eukaryot Microbiol. 2009;56:34–8.
Lobo ML, Xiao L, Cama V, Stevens T, Antunes F, Matos O. Genotypes of Enterocytozoon bieneusi in mammals in Portugal. J Eukaryot Microbiol. 2006;53:S61–4.
Santín M, Trout JM, Fayer R. Enterocytozoon bieneusi genotypes in dairy cattle in the eastern United States. Parasitol Res. 2005;97:535–8.
Santín M, Trout JM, Vecino JA, Dubey JP, Fayer R. Cryptosporidium, Giardia and Enterocytozoon bieneusi in cats from Bogota (Colombia) and genotyping of isolates. Vet Parasitol. 2006;141:334–9.
Sulaiman IM, Bern C, Gilman R, Cama V, Kawai V, Vargas D, et al. A molecular biologic study of Enterocytozoon bieneusi in HIV infected patients in Lima, Peru. J Eukaryot Microbiol. 2003;50:591–6.
Abe N, Kimata I, Iseki M. Molecular evidence of Enterocytozoon bieneusi in Japan. J Vet Med Sci. 2009;71:217–9.
Wang L, Zhang H, Zhao X, Zhang L, Zhang G, Guo M, et al. Zoonotic Cryptosporidium species and Enterocytozoon bieneusi genotypes in HIV-positive patients on antiretroviral therapy. J Clin Microbiol. 2013;51:557–63.
Wang L, Xiao L, Duan L, Ye J, Guo Y, Guo M, et al. Concurrent infections of Giardia duodenalis, Enterocytozoon bieneusi, and Clostridium difficile in children during a cryptosporidiosis outbreak in a pediatric hospital in China. PLoS Negl Trop Dis. 2013;7:e2437.
Ye J, Xiao L, Ma J, Guo M, Liu L, Feng Y. Anthroponotic enteric parasites in monkeys in public park, China. Emerg Infect Dis. 2012;18:1640–3.
Ye J, Xiao L, Li J, Huang W, Amer SE, Guo Y, et al. Occurrence of human-pathogenic Enterocytozoon bieneusi, Giardia duodenalis and Cryptosporidium genotypes in laboratory macaques in Guangxi, China. Parasitol Int. 2014;63:132–7.
Ye J, Xiao L, Wang Y, Guo Y, Roellig DM, Feng Y. Dominance of Giardia duodenalis assemblage A and Enterocytozoon bieneusi genotype BEB6 in sheep in Inner Mongolia, China. Vet Parasitol. 2015;210:235–9.
Jiang Y, Tao W, Wan Q, Li Q, Yang Y, Lin Y, et al. Zoonotic and potentially host-adapted Enterocytozoon bieneusi genotypes in sheep and cattle in Northeast China and an increasing concern about the zoonotic importance of previously considered ruminant-adapted genotypes. Appl Environ Microbiol. 2015;81:3326–35.
Zhao W, Zhang W, Wang R, Liu W, Liu A, Yang D, et al. Enterocytozoon bieneusi in sika deer (cervus nippon) and red deer (cervus elaphus): deer specificity and zoonotic potential of ITS genotypes. Parasitol Res. 2014;113:4243–50.
Zhao W, Zhang W, Yang D, Zhang L, Wang R, Liu A. Prevalence of Enterocytozoon bieneusi and genetic diversity of ITS genotypes in sheep and goats in China. Infect Genet Evol. 2015;32:265–70.
Li N, Xiao L, Wang L, Zhao S, Zhao X, Duan L, et al. Molecular surveillance of Cryptosporidium spp., Giardia duodenalis, and Enterocytozoon bieneusi by genotyping and subtyping parasites in wastewater. PLoS Negl Trop Dis. 2012;6:e1809.
Thellier M, Breton J. Enterocytozoon bieneusi in human and animals, focus on laboratory identification and molecular epidemiology. Parasite. 2008;15:349–58.
Henriques-Gil N, Haro M, Izquierdo F, Fenoy S, del Aguila C. Phylogenetic approach to the variability of the microsporidian Enterocytozoon bieneusi and its implications for inter- and intrahost transmission. Appl Environ Microbiol. 2010;76:3333–42.
Karim MR, Dong H, Li T, Yu F, Li D, Zhang L, et al. Predomination and new genotypes of Enterocytozoon bieneusi in captive nonhuman primates in zoos in China: high genetic diversity and zoonotic significance. PLoS One. 2015;10:e0117991.
Du SZ, Zhao GH, Shao JF, Fang YQ, Tian GR, Zhang LX, et al. Cryptosporidium spp. Giardia intestinalis, and Enterocytozoon bieneusi in Captive Non-Human Primates in Qinling Mountains. Korean J Parasitol. 2015;53:395–402.
Zhao GH, Du SZ, Wang HB, Hu XF, Deng MJ, Yu SK, et al. First report of zoonotic Cryptosporidium spp., Giardia intestinalis and Enterocytozoon bieneusi in golden takins (Budorcas taxicolor bedfordi). Infect Genet Evol. 2015;34:394–401.
Li W, Li Y, Song M, Lu Y, Yang J, Tao W, et al. Prevalence and genetic characteristics of Cryptosporidium, Enterocytozoon bieneusi and Giardia duodenalis in cats and dogs in Heilongjiang province, China. Vet Parasitol. 2015;208:125–34.
Yang Y, Lin Y, Li Q, Zhang S, Tao W, Wan Q, et al. Widespread presence of human-pathogenic Enterocytozoon bieneusi genotype D in farmed foxes (Vulpes vulpes) and raccoon dogs (Nyctereutes procyonoides) in China: first identification and zoonotic concern. Parasitol Res. 2015;114:4341–8.
Zhao W, Zhang W, Yang F, Zhang L, Wang R, Cao J, et al. Enterocytozoon bieneusi in dairy cattle in the Northeast of China: genetic diversity of ITS gene and evaluation of zoonotic transmission potential. J Eukaryot Microbiol. 2015;62:553–60.
Guo Y, Alderisio KA, Yang W, Cama V, Feng Y, Xiao L. Host specificity and source of Enterocytozoon bieneusi genotypes in a drinking source watershed. Appl Environ Microbiol. 2014;80:218–25.
Tian GR, Zhao GH, Du SZ, Hu XF, Wang HB, Zhang LX, et al. First report of Enterocytozoon bieneusi from giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens) in China. Infect Genet Evol. 2015;34:32–5.
Ma J, Li P, Zhao X, Xu H, Wu W, Wang Y, et al. Occurrence and molecular characterization of Cryptosporidium spp. and Enterocytozoon bieneusi in dairy cattle, beef cattle and water buffaloes in China. Vet Parasitol. 2015;207:220–7.
Ma J, Cai J, Ma J, Feng Y, Xiao L. Enterocytozoon bieneusi genotypes in yaks (Bos grunniens) and their public health potential. J Eukaryot Microbiol. 2015;62:21–5.
Zhang Z, Huang J, Karim MR, Zhao J, Dong H, Ai W, et al. Zoonotic Enterocytozoon bieneusi genotypes in Pere David’s deer (Elaphurus davidianus) in Henan, China. Exp Parasitol. 2015;155:46–8.
This work was supported by the Earmarked Fund for China Modern Agro-industry Technology Research System (nycytx-39) and the Collaborative Innovation Center of Modern Animal Husbandry, Henan Province, China. We thanks for Professor Una Ryan from School of Veterinary and Life Sciences, Murdoch University, Murdoch, Western Australia, to copy-edit our MS.
The authors declare that they have no competing interests.
CSN conceived and designed the experiments; KS, MJL, and XXW performed the experiments; KS, JQL and MRK analyzed the data; KS, RJW, LXZ, FCJ and CSN wrote the manuscript. All the authors have read and approved the final version of the manuscript.
Prevalence and genotype distribution of E. bieneusi in goats and sheep on each farm. (DOC 117 kb)
Nucleotide substitutions among some new genotypes of E. bieneusi from this study versus the reported genotype BEB6 (EU 153584). (DOC 132 kb)
FIG 1 Phylogenetic relationships of the E. bieneusi genotypes identified in this study and other reported genotypes, as inferred by neighbor-joining analysis of ITS gene sequences based on the distances calculated using the Kimura 2-parameter model. Bootstrap values >50 % from 1000 replicates are shown on the nodes. The tree was rooted with GenBank sequence AF059610 from a dog. Known genotypes observed in goats and sheep are marked with open squares and triangles, and the new genotypes in goats and sheep from this study are indicated by filled squares and triangles. (TIF 2078 kb)