Skip to main content

Identification and genotyping of Enterocytozoon bieneusi in wild Himalayan marmots (Marmota himalayana) and Alashan ground squirrels (Spermophilus alashanicus) in the Qinghai-Tibetan Plateau area (QTPA) of Gansu Province, China



Enterocytozoon bieneusi is the most frequently detected microsporidian species in humans and animals. Currently, to the best of our knowledge, no information on E. bieneusi infection in Himalayan marmots (Marmota himalayana) and Alashan ground squirrels (Spermophilus alashanicus) is available worldwide. The aim of the present study was to understand the occurrence and genetic characterizations of E. bieneusi in Himalayan marmots and Alashan ground squirrels in the Qinghai-Tibetan Plateau area (QTPA), Gansu Province, China.


A total of 498 intestinal contents were collected from 399 Himalayan marmots and 99 Alashan ground squirrels in QTPA. These samples were screened for the presence of E. bieneusi by using nested polymerase chain reaction and sequencing of the internal transcribed spacer (ITS) region of the ribosomal RNA (rRNA) gene. The ITS-positive sequences were aligned and phylogenetically analyzed to determine the genotypes of E. bieneusi.


The average infection rate of E. bieneusi was 10.0% (50/498), with 11.8% (47/399) in Himalayan marmots and 3.0% (3/99) in Alashan ground squirrels. A total of 7 distinct E. bieneusi genotypes were confirmed: 1 known genotype, YAK1 (n = 18) and 6 novel genotypes, named as ZY37 (n = 27), HN39 (n = 1), HN96 (n = 1), SN45 (n = 1), XH47 (n = 1) and ZY83 (n = 1). All the genotypes obtained in the present study were classified into group 1.


To our knowledge, this is the first report of E. bieneusi in Himalayan marmots and Alashan ground squirrels in China. The identification of genotype YAK1 in the two rodent species expanded the host range of this genotype. All the seven genotypes were clustered into zoonotic group 1, suggesting that these animal species can be potential epidemiological vectors of zoonotic microsporidiosis caused by E. bieneusi and pose a threat to ecological security. It is necessary to strengthen management practices and surveillance in the investigated areas to reduce the risk of E. bieneusi infection from the two rodent species to humans.


Microsporidia are a large phylum of fungus-related pathogens with an obligate intracellular lifestyle, which are composed of at least 1500 species within 200 genera, and 17 of these species have been reported to cause human infections [1, 2]. Enterocytozoon bieneusi is the most frequently detected species (more than 90%) in microsporidiosis cases [3]. Microsporidiosis caused by E. bieneusi is characterized by asymptomatic infection or self-limiting diarrhea in immunocompetent and healthy humans, and life-threatening and persistent diarrhea in immunocompromised individuals, such as HIV-positive patients, cancer patients, organ transplant recipients, and children [3]. Enterocytozoon bieneusi has also been widely reported in numerous mammalian and avian hosts [4]. In fact, wild animals can play a role in spreading microsporidium spores to humans and farmlands [5]. Consumption of E. bieneusi spores from human and animal feces via contaminated food or water is the most probable transmission route of E. bieneusi to humans or animals [3]. Enterocytozoon bieneusi has been detected in foods such as milk, raspberries, cucumbers, beans and lettuce, which can act as intermediate carriers of this parasite [6,7,8]. In 2009, a food-borne microsporidiosis outbreak caused by E. bieneusi was reported in Sweden [7]. No water-borne microsporidiosis outbreak has been reported to date; however, E. bieneusi spores have been found in multiple water bodies, for example drinking-source watersheds [9], recreational watersheds [10], and wastewater resources used for irrigation of farmlands [11], indicating the possibility of water-borne transmission.

Significant progress has been made in understanding the host specificity and transmission patterns of E. bieneusi by using genotyping tools and phylogenetic analyses. Currently, the standard method for identifying and genotyping E. bieneusi isolates is the analysis of genetic polymorphisms in the ribosomal internal transcribed spacer (ITS) region [12]. To date, almost 500 genotypes have been identified in a wide range of host species and divided into 11 different genetic groups [13]. At least 106 genotypes have been found in humans, with 46 of them detected in animals [14]. Groups 1 and 2 are zoonotic in nature and mostly composed of genotypes that spread in human and animal hosts [15]. To better understand the epidemiology of E. bieneusi and assess the role of animals in its transmission to humans, molecular epidemiological studies should be conducted for genotyping E. bieneusi isolates from under-sampled animal hosts.

Rodents are well-known hosts or reservoirs for a number of zoonotic infectious diseases (e.g. plague, leptospirosis, leishmaniasis, salmonellosis and viral hemorrhagic fevers) and may act as possible mediators of these diseases [16]. Previous epidemiological studies of E. bieneusi in rodents have identified 56 valid genotypes [9, 17,18,19,20,21,22,23,24,25,26,27,28], with 31 and 8 genotypes belonging to group 1 and 2, respectively (Additional file 1: Table S1). However, there is limited genetic information on E. bieneusi in rodent species belonging to genera Marmota and Spermophilus. To date, only one study has reported the occurrence of E. bieneusi in a species of the genus Marmota, i.e. wild woodchucks (Marmota monax) [9]. In China, Himalayan marmots (Marmota himalayana) and Alashan ground squirrels (Spermophilus alashanicus) are common wild rodents in plateau areas, and they often share the same pasture with humans, herbivorous animals, and other wild animals. However, to the best of our knowledge, no reports are available on the occurrence and genetic characterizations of E. bieneusi in the two rodent species in these areas. The aims of the present study were to determine the epidemiological status and genotypes of E. bieneusi in the two wild animals in the Qinghai-Tibetan Plateau Area (QTPA) in China by nested polymerase chain reaction (PCR) and sequence analysis of the ITS region of the ribosomal RNA (rRNA) gene, and assess the potential zoonotic transmission.


Sample collection

A total of 498 intestinal contents were collected from 399 Himalayan marmots and 99 Alashan ground squirrels from June to September 2017 in QTPA in Gansu Province (geographical coordinates: 32°31′–42°57′N latitude, 92°13′–108°46′E longitude), China. The Himalayan marmots were obtained from 4 cities/counties, namely Luqu (n = 98), Sunan (n = 100), Xiahe (n = 102) and Zhangye (n = 99), whereas the Alashan ground squirrels were from Huining County (n = 99) (Table 1). The animals were captured alive with mousetraps in the field, and the intestinal contents were collected after euthanizing the animals with a high-dose of CO2 in the laboratory of the local Center for Disease Control and Prevention (CDC) with security measures. All the collected intestinal content samples were stored in 2 ml sterile tubes at − 20 °C, and then transported in a cooler with ice packs to our laboratory prior to DNA extraction in Shanghai.

Table 1 Prevalence and distribution of E. bieneusi genotypes in wild Himalayan marmots and Alashan ground squirrels

DNA extraction

Before DNA extraction, each intestinal content sample was transferred (up to 200 mg) to a 2 ml microcentrifuge tube containing 200 mg of magnetic beads. The samples were then vortexed (Voterx-5 Qilinbeier, Scientific Industries, New York, USA) at maximum speed and incubated at 56 °C until the intestinal content samples were completely lysed. Then DNA was extracted using the DNeasy Blood & Tissue Kit (Cat. #69506; Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The extracted DNA was stored at − 20 °C in a freezer until further use.

PCR amplification

Enterocytozoon bieneusi was identified and genotyped using nested PCR, in which an approximately 390 bp nucleotide fragment, including 243 bp of the ITS region, was amplified [29]. Two pairs of primers, EBITS3/EBITS4 and EBITS1/EBITS2, were used in the first and second PCR amplifications, respectively. The two cycling parameters were as follows: 35 cycles at 94 °C for 30 s, 57 °C for 30 s and 72 °C for 40 s, and 30 cycles at 94 °C for 30 s, 55 °C for 30 s and 72 °C for 40 s, with both having a final extension step at 72 °C for 10 min [30]. A positive control (DNA of human-derived genotype D) and a negative control (DNase-free water) were used for the primary and secondary PCR tests to ensure accuracy. TaKaRa Taq DNA polymerase (TaKaRa Bio Inc., Tokyo, Japan) was used and the secondary PCR products were visualized using agarose gel electrophoresis and staining the 1.5% gel with GelRed (Biotium Inc., Hayward, CA, USA).

Nucleotide sequencing and analysis

MEGA 7 ( was used to align and compare the nucleotide sequences obtained in the present study with each other and the referenced E. bieneusi sequences obtained from GenBank ( If the genotypes obtained in this study were identical to the known genotypes in GenBank, they would be given the first published names. If the ITS sequences had nucleotide deletions, insertions, and substitutions that were different from those of known genotypes, they were considered as novel genotypes. All the genotypes were identified on the basis of the 243 bp of the ITS region of the rRNA gene of E. bieneusi according to the established nomenclature system [29].

Phylogenetic analysis

A Maximum Likelihood tree was constructed for E. bieneusi by using MEGA 7 ( to evaluate the genetic relationship between the ITS genotypes of E. bieneusi obtained in this study and known genotypes deposited in GenBank. The Tamura 3-parameter model was used, and bootstrap analysis with 1000 replicates was performed to assess the reliability of the tree. Nucleotide sequences of Enterocytozoon hepatopenaei from shrimp (GenBank: MH260592) and E. bieneusi from the red kangaroo (GenBank: KY706128) were used as the outgroup for the phylogenetic analysis.


Occurrence of E. bieneusi

A total of 498 intestinal contents of Himalayan marmots and Alashan ground squirrels were screened for the presence of E. bieneusi by amplifying the ITS region with nested PCR. Fifty samples were successfully amplified and sequenced; the overall infection rate of E. bieneusi was 10.0% (50/498), with 11.8% (47/399) for Himalayan marmots and 3.0% (3/99) for Alashan ground squirrels. Enterocytozoon bieneusi was detected in all five investigated areas. For Himalayan marmots, the infection rates varied from 7.8% to 17.2% in 4 areas (Table 1).

Genetic characterizations and distribution of E. bieneusi genotypes

Based on sequence analysis of the ITS region, 7 genotypes were identified from the 50 E. bieneusi isolates: 1 known genotype, YAK1 (GenBank: MK843246, n = 18), and 6 novel genotypes, named as ZY37 (GenBank: MN378367, n = 27), HN39 (GenBank: MN378365, n = 1), HN96 (GenBank: MN378366, n = 1), SN45 (GenBank: MN378369, n = 1), XH47 (GenBank: MN378370, n = 1) and ZY83 (GenBank: MN378371, n = 1). Genotypes ZY37, YAK1, SN45, XH47 and ZY83 were detected in the Himalayan marmots, with genotype ZY37 showing a predominance (57.4%, 27/47). Genotypes HN39, HN96 and YAK1 were detected in the Alashan ground squirrels (Table 1). A total of 14 nucleotide polymorphic sites were observed, and each genotype differed by 1–9 single-nucleotide polymorphisms (Table 2).

Table 2 Nucleotide variation in ITS gene region of E. bieneusi isolates obtained in this study

The novel genotype HN39 had 99.59% similarity with genotype D (syn. CEbC, NCF7, Peru9, PigEBITS9, PtEb VI, SHW1 and WL8), with one base substitution (A to G). The other five novel genotypes (ZY37, HN96, SN45, XH47 and ZY83) showed the highest similarity with genotype YAK1: 99.59% for genotypes HN96 (C to T), SN45 (A to G) and XH47 (G to T), with one base substitution, and 99.18% for genotypes ZY37 (T to -, G to -) and ZY83 (G to -, T to -), with two base deletions.

Phylogenetic relationship of E. bieneusi genotypes

Phylogenetic analysis of the Maximum Likelihood tree for the ITS sequences of E. bieneusi showed that all the genotypes obtained in the present study were classified into group 1. Genotypes ZY37, YAK1, HN96, SN45, XH47 and ZY83 were clustered into one branch, and genotype HN39 was clustered into another branch (Fig. 1).

Fig. 1

Phylogenetic relationships of E. bieneusi genotypes identified in this study and known genotypes published on GenBank inferred using the Maximum-Likelihood method. Genetic distances of the ITS sequences were calculated using the Tamura 3-parameter model. The numbers above the branches are bootstrap values from 1000 replicates. The group terminology for the clusters is based on the Li W et al. [14]. The black triangles and squares indicate the known genotype YAK1 and six novel genotypes identified in the present study, respectively


Rodents, especially wild ones, can act as vectors for numerous pathogens, and some of the pathogens, including E. bieneusi are zoonotic. Enterocytozoon bieneusi has been reported in some wild rodent species: with a prevalence of 38.9% (121/311) in wild bank voles and wild mice from Poland [25]; 26.8% (38/142) in wild beavers, chipmunks, deer mice, muskrats, porcupines, squirrels, voles, and woodchucks from New York, USA [9]; 10.7% (31/289) in wild mice from Germany/Czech Republic [24]; 10.2% (33/324) in wild beavers and muskrats from Maryland, USA [23]; 7.9% (19/242) in brown rats from Heilongjiang Province, China [19]; 4.0% (8/199) in brown rats and wild mice from Henan Province, China [28]; and 1.1% (3/280) in wild mice from Slovakia [27]. In the present study, E. bieneusi was detected in Himalayan marmots (11.8%, 47/399) and Alashan ground squirrels (3.0%, 3/99). The prevalence of E. bieneusi is difficult to compare because it may be influenced by many factors, such as the immune status of the hosts, number of samples, experimental methods, and climatic and geographical differences.

Because of the high degree of genetic polymorphisms in the ITS region of the rRNA gene, sequence analysis of the ITS region is regarded as the standard method for identifying and genotyping E. bieneusi isolates [29]. In 2002, Buckholt et al. [30] proposed that if enough isolates are sequenced, all 243 nucleotides in the ITS region can be found to be polymorphic. In fact, numerous studies on genotyping of the ITS region of E. bieneusi isolates provide further information on the transmission routes and genetic characterization of E. bieneusi. In the present study, seven genotypes were identified from 50 E. bieneusi isolates, with genotypes ZY37, YAK1, SN45, XH47 and ZY83 in the Himalayan marmots and genotypes YAK1, HN39 and HN96 in the Alashan ground squirrels. To date, 63 genotypes of E. bieneusi including those obtained in this study, have been identified in seven rodent families from seven countries: six genotypes (D, EbpC, WL7, WL9, WL12 and WL15) in the family Castoridae; one genotype (Peru16) in the family Caviidae; two genotypes (BEB6 and D) in the family Chinchillidae; 16 genotypes (D, EbpC, Peru11, Type IV, WL4, WL6, WL10, WL14, WL15, WL20, WL21, WL23, WL25, WR2, WR6, WR10) in the family Cricetidae; 22 genotypes (BEB6, C, CD6, CHG2, CZ3, D, EbpA, gorilla 1, H, Peru6, Peru8, Peru16, PigEBITS5, S6, WR1, WR3 to WR9) in the family Muridae; and 27 genotypes (CE01, CHG9, D, EbpC, HN39, HN96, Horse2, Nig7, PtEb V, S7, SC02, SCC-1 to SCC-4, SN45, Type IV, WL4, WL20 to WL23, WW6, XH47, YAK1, ZY37 and ZY83) in the family Sciuridae; 6 genotypes (BR1, BR2, D, EbpA, J, PigEBITS7) in the family Spalacidae. In humans, 34.9% (22/63) of the 63 genotypes have been identified: C, CZ3, D, EbpA, EbpC, H, Peru 6, Peru8, Peru11, Peru16, PigEBITS5, PigEBITS7, S6, SC02, Type IV, WL7, WL12, and WL15 in zoonotic group 1; BEB6 and J in group 2; Nig 7 in group 6; and S7 in group 10 (Additional file 1: Table S1). Zoonotic E. bieneusi genotypes have been identified in 12 of 13 epidemiological studies. The prevalence of zoonotic genotypes was observed to vary in different host species ranging between 0–100% (Fig. 2). All the E. bieneusi genotypes have been found in humans in some epidemiological studies, including those conducted in chinchillas and brown rats in China [18, 19], guinea pigs in Peru [26], and house mice in Germany/Czech Republic and Slovakia [24, 27] (Additional file 1: Table S1). These results indicate that rodents infected with E. bieneusi may play an important role in its transmission to humans and become an important source of water contamination in the areas investigated in this study.

Fig. 2

Percentage of zoonotic E. bieneusi genotypes in rodents worldwide. 1, beaver (Castor fiber); 2, guinea pig (Cavia porcellus); 3, chinchilla (Chinchilla lanigera); 4, bank vole (Myodes glareolus); 5, boreal red-backed vole (Myodes gapperi) and meadow vole (Microtus pennsylvanicus) (No separate genotypes information provided in the original reference); 6, deer mouse (Peromyscus sp.); 7, muskrat (Ondatra zibethicus); 8, brown rat (Rattus norvegicus); 9, brown rat (Rattus norvegicus) and house mouse (Mus musculus) (No separate genotypes information provided in the original reference); 10, house mouse (Mus musculus); 11, striped field mouse (Apodemus agrarius); 12, yellow-necked mouse (Apodemus flavicollis); 13, Alashan ground squirrel (Spermophilus alashanicus); 14, chipmunk (Eutamias asiaticus); 15, eastern chipmunk (Tamias striatus); 16, eastern gray squirrel (Sciurus carolinensis), red squirrel (Sciurus vulgaris) and southern flying squirrel (Glaucomys volans) (No separate genotypes information provided in the original reference)); 17, Himalayan marmot (Marmota himalayana); 18, red-bellied tree squirrel (Callosciurus erythraeus); 19, woodchuck (Marmota monax). 20, bamboo rat (Rhizomys sinensis). The bar chart was produced from the data presented in Additional file 1: Table S1

To date, the known genotype YAK1 has been detected in only domestic yaks in QTPA, Qinghai Province, China [31]. Although no zoonotic E. bieneusi genotypes were obtained in this study, genotype YAK1 and six novel genotypes were clustered into group 1 in the phylogenetic analysis, indicating the zoonotic ability and public health significance of these genotypes. In fact, the genotype HN39 shows only one base difference when compared with D (pathogenic in humans). In the investigated areas, wild Himalayan marmots and Alashan ground squirrels share the same pasture with humans, herbivorous animals (cattle, sheep and yaks, etc.), and other wild animals (plateau pikas and Tibetan antelopes, etc.). Wild animals usually have strong migration habits and a wide range of activity. Infection occurs through oral ingestion of the environmentally resistant spores of E. beneusi from contaminated feces, food, drink and pasture (for grazing animals), which provides feasible opportunities for the transmission of the spores from wild Himalayan marmots and Alashan ground squirrels to humans and animals. Therefore, these animals infected with E. bieneusi can be potential sources of microsporidiosis in humans and other animals and pose a threat to public health and ecological security.

Generally, the ITS region is 243-bp long in natural E. bieneusi genotypes. However, in the present study, amplicons of 241 bp were observed in the ITS region of the novel genotypes ZY37 and ZY83. In fact, length variations in the ITS region of E. bieneusi have been found in at least 15 E. bieneusi genotypes. Amplicons of less than 243 bp in the ITS region have been reported for a few genotypes: 241 bp in genotype MAY1 from a renal transplant recipient in France [32]; CHN-RR1 from rex rabbits in China [33]; CHN3, CHN4 and CHN5 from diarrheal children in China [34]; YN249 from a Yao person in China [35]; and PigEb12 and PigEb15 from pigs in Brazil [38] and 242 bp in genotype CAF4 from HIV-positive patients in Gabon and villagers in Cameroon [36]. In contrast, amplicons of more than 243 bp in the ITS region have been found in some genotypes: 244 bp in genotype SCBB1 from captive black bears in China [37], as well as genotypes HNM-V, HNM-VI, and HNM-VII from captive long-tailed macaques in China [15]; 245 bp in genotype PigEb16 from pigs in Brazil [38], and genotype Q from patients with diarrhea in Switzerland [39]. These observations indicate that base deletions and insertions can occur in the ITS region of the rRNA gene of E. bieneusi.


To the best of our knowledge, this is the first report of the identification and genotyping of E. bieneusi in wild Himalayan marmots and Alashan ground squirrels in QTPA, Gansu Province, China. All seven genotypes (ZY37, YAK1, HN39, HN96, SN45, XH47 and ZY83) were classified into group 1, suggesting the possibility of zoonotic and cross-species transmission. Further molecular epidemiological surveys are required to assess the importance and transmission dynamics of E. bieneusi in other animals and humans, including local residents and tourists in close contact with pasture. Currently, no effective treatment or vaccine is available for microsporidiosis caused by E. bieneusi. Therefore, it is necessary to make people aware of the potential risk of zoonotic transmission of this disease and importance of strengthening management practices and surveillance for E. bieneusi in wild animals.

Availability of data and materials

The data supporting the conclusions of this article are included within the article and its additional file. Nucleotide sequences were deposited in the GenBank database under the Accession numbers MN378365-MN378367, MN378369-MN378371.



Qinghai-Tibetan Plateau area


ribosomal internal transcribed spacer


polymerase chain reaction


  1. 1.

    Stentiford GD, Becnel JJ, Weiss LM, Keeling PJ, Didier ES, Williams BP, et al. Microsporidia—emergent pathogens in the global food chain. Trends Parasitol. 2016;32:336–48.

    CAS  Article  Google Scholar 

  2. 2.

    Fayer R, Santín M. Epidemiology of microsporidia in human infections. In: Weiss LM, Becnel JJ, editors. Microsporidia: pathogens of opportunity. Chichester: Wiley; 2014. p. 135–64.

    Google Scholar 

  3. 3.

    Qiu L, Xia W, Li W, Ping J, Ding S, Liu H. The prevalence of microsporidia in China: a systematic review and meta-analysis. Sci Rep. 2019;9:3174.

    Article  Google Scholar 

  4. 4.

    Santín M, Fayer R. Microsporidiosis: Enterocytozoon bieneusi in domesticated and wild animals. Res Vet Sci. 2011;90:363–71.

    Article  Google Scholar 

  5. 5.

    Javanmard E, Nemati S, Sharifdini M, Rostami A, Mirjalali H, Zali MR. The First report and molecular analysis of Enterocytozoon bieneusi from Raccoon (Procyon lotor) in North of Iran. J Eukaryot Microbiol. 2020;67:359–68.

    CAS  Article  Google Scholar 

  6. 6.

    Lee JH. Molecular detection of Enterocytozoon bieneusi and identification of a potentially human-pathogenic genotype in milk. Appl Environ Microbiol. 2008;74:1664–6.

    CAS  Article  Google Scholar 

  7. 7.

    Decraene V, Lebbad M, Botero-Kleiven S, Gustavsson AM, Löfdahl M. First reported foodborne outbreak associated with microsporidia, Sweden, October 2009. Epidemiol Infect. 2012;140:519–27.

    CAS  Article  Google Scholar 

  8. 8.

    Jedrzejewski S, Graczyk TK, Slodkowicz-Kowalska A, Tamang L, Majewska AC. Quantitative assessment of contamination of fresh food produce of various retail types by human-virulent microsporidian spores. Appl Environ Microbiol. 2007;73:4071–3.

    CAS  Article  Google Scholar 

  9. 9.

    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.

    CAS  Article  Google Scholar 

  10. 10.

    Hu Y, Feng Y, Huang C, Xiao L. Occurrence, source, and human infection potential of Cryptosporidium and Enterocytozoon bieneusi in drinking source water in Shanghai, China, during a pig carcass disposal incident. Environ Sci Technol. 2014;48:14219–27.

    CAS  Article  Google Scholar 

  11. 11.

    Javanmard E, Mirjalali H, Niyyati M, Jalilzadeh E, Seyed Tabaei SJ, Asadzadeh Aghdaei H. Molecular and phylogenetic evidences of dispersion of human-infecting microsporidia to vegetable farms via irrigation with treated wastewater: one-year follow up. Int J Hyg Environ Health. 2018;221:642–51.

    Article  Google Scholar 

  12. 12.

    Santín M. Enterocytozoon bieneusi. In: Xiao L, Ryan U, Feng Y, editors. Biology of foodborne parasites. Boca Raton: CRC Press; 2015. p. 149–74.

    Google Scholar 

  13. 13.

    Li W, Feng Y, Zhang L, Xiao L. Potential impacts of host specificity on zoonotic or interspecies transmission of Enterocytozoon bieneusi. Infect Genet Evol. 2019;75:104033.

    CAS  Article  Google Scholar 

  14. 14.

    Li W, Feng Y, Santín M. Host specificity of Enterocytozoon bieneusi and public health implications. Trends Parasitol. 2019;35:436–51.

    CAS  Article  Google Scholar 

  15. 15.

    Zhao W, Zhou H, Jin H, Sun L, Li P, Liu M, et al. Genotyping of Enterocytozoon bieneusi among captive long-tailed macaques (Macaca fascicularis) in Hainan Province: high genetic diversity and zoonotic potential. Acta Trop. 2020;201:105211.

    CAS  Article  Google Scholar 

  16. 16.

    Rabiee MH, Mahmoudi A, Siahsarvie R, Kryštufek B, Mostafavi E. Rodent-borne diseases and their public health importance in Iran. PLoS Negl Trop Dis. 2018;12:e0006256.

    Article  Google Scholar 

  17. 17.

    Deng L, Li W, Yu X, Gong C, Liu X, Zhong Z, et al. First report of the human-pathogenic Enterocytozoon bieneusi from red-bellied tree squirrels (Callosciurus erythraeus) in Sichuan, China. PLoS ONE. 2016;11:e0163605.

    Article  Google Scholar 

  18. 18.

    Qi M, Luo N, Wang H, Yu F, Wang R, Huang J, et al. Zoonotic Cryptosporidium spp. and Enterocytozoon bieneusi in pet chinchillas (Chinchilla lanigera) in China. Parasitol Int. 2015;64:339–41.

    Article  Google Scholar 

  19. 19.

    Zhao W, Wang J, Ren G, Yang Z, Yang F, Zhang W, et al. Molecular characterizations of Cryptosporidium spp and Enterocytozoon bieneusi in brown rats (Rattus norvegicus) from Heilongjiang Province, China. Parasites Vectors. 2018;11:313.

    Article  Google Scholar 

  20. 20.

    Deng L, Li W, Zhong Z, Chai Y, Yang L, Zheng H, et al. Molecular characterization and new genotypes of Enterocytozoon bieneusi in pet chipmunks (Eutamias asiaticus) in Sichuan Province, China. BMC Microbiol. 2018;18:37.

    Article  Google Scholar 

  21. 21.

    Wang H, Liu Q, Jiang X, Zhang Y, Zhao A, Cui Z, et al. Dominance of zoonotic genotype D of Enterocytozoon bieneusi in bamboo rats (Rhizomys sinensis). Infect Genet Evol. 2019;73:113–8.

    Article  Google Scholar 

  22. 22.

    Roellig DM, Salzer JS, Carroll DS, Ritter JM, Drew C, Gallardo-Romero N, et al. Identification of Giardia duodenalis and Enterocytozoon bieneusi in an epizoological investigation of a laboratory colony of prairie dogs, Cynomys ludovicianus. Vet Parasitol. 2015;210:91–7.

    Article  Google Scholar 

  23. 23.

    Sulaiman IM, Fayer R, Lal AA, Trout JM, Schaefer 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.

    CAS  Article  Google Scholar 

  24. 24.

    Sak B, Kváč M, Květoňová D, Albrecht T, Piálek J. The first report on natural Enterocytozoon bieneusi and Encephalitozoon spp. infections in wild east-European house mice (Mus musculus musculus) and west-European house mice (M. m. domesticus) in a hybrid zone across the Czech Republic-Germany border. Vet Parasitol. 2011;178:246–50.

    Article  Google Scholar 

  25. 25.

    Perec-Matysiak A, Buńkowska-Gawlik K, Kváč M, Sak B, Hildebrand J, Leśniańska K. Diversity of Enterocytozoon bieneusi genotypes among small rodents in southwestern Poland. Vet Parasitol. 2015;214:242–6.

    Article  Google Scholar 

  26. 26.

    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.

    CAS  Article  Google Scholar 

  27. 27.

    Danišová O, Valenčáková A, Stanko M, Luptáková L, Hasajová A. First report of Enterocytozoon bieneusi and Encephalitozoon intestinalis infection of wild mice in Slovakia. Ann Agric Environ Med. 2015;22:251–2.

    Article  Google Scholar 

  28. 28.

    Yu F, Qi M, Zhao Z, Lv C, Wang Y, Wang R, et al. The potential role of synanthropic rodents and flies in the transmission of Enterocytozoon bieneusi on a dairy cattle farm in China. J Eukaryot Microbiol. 2019;66:435–41.

    Article  Google Scholar 

  29. 29.

    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.

    Article  Google Scholar 

  30. 30.

    Buckholt MA, Lee JH, Tzipori S. Prevalence of Enterocytozoon bieneusi in swine: an 18-month survey at a slaughterhouse in Massachusetts. Appl Environ Microbiol. 2002;68:2595–9.

    CAS  Article  Google Scholar 

  31. 31.

    Zhang Q, Zhang Z, Ai S, Wang X, Zhang R, Duan Z. Cryptosporidium spp., Enterocytozoon bieneusi, and Giardia duodenalis from animal sources in the Qinghai-Tibetan Plateau Area (QTPA) in China. Comp Immunol Microbiol Infect Dis. 2019;67:101346.

    Article  Google Scholar 

  32. 32.

    Pomares C, Santín M, Miegeville M, Espern A, Albano L, Marty P, et al. A new and highly divergent Enterocytozoon bieneusi genotype isolated from a renal transplant recipient. J Clin Microbiol. 2012;50:2176–8.

    Article  Google Scholar 

  33. 33.

    Yang Z, Zhao W, Shen Y, Zhang W, Shi Y, Ren G, et al. Subtyping of Cryptosporidium cuniculus and genotyping of Enterocytozoon bieneusi in rabbits in two farms in Heilongjiang Province, China. Parasite. 2016;23:52.

    Article  Google Scholar 

  34. 34.

    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.

    CAS  Article  Google Scholar 

  35. 35.

    Gong B, Yang Y, Liu X, Cao J, Xu M, Xu N, et al. First survey of Enterocytozoon bieneusi and dominant genotype Peru6 among ethnic minority groups in southwestern China’s Yunnan Province and assessment of risk factors. PLoS Negl Trop Dis. 2019;13:e0007356.

    CAS  Article  Google Scholar 

  36. 36.

    Breton J, Bart-Delabesse E, Biligui S, Carbone A, Seiller X, Okome-Nkoumou M, et al. New highly divergent rRNA sequence among biodiverse genotypes of Enterocytozoon bieneusi strains isolated from humans in Gabon and Cameroon. J Clin Microbiol. 2007;45:2580–9.

    CAS  Article  Google Scholar 

  37. 37.

    Huang X, Zhou Z, Liu H, Deng L, Bi B, Chai Y, et al. New genotypes and molecular characterization of Enterocytozoon bieneusi in captive black bears in China. Int J Parasitol Parasites Wildl. 2019;10:1–5.

    Article  Google Scholar 

  38. 38.

    Fiuza VR, Oliveira FC, Fayer R, Santín M. First report of Enterocytozoon bieneusi in pigs in Brazil. Parasitol Int. 2015;64:18–23.

    CAS  Article  Google Scholar 

  39. 39.

    Rinder H, Dengjel B, Gothe R, Loscher T. Microsporidosis of man: where is the reservoir? Mitt Osterr Ges Trop Med Parasitol. 2000;22:1–6.

    Google Scholar 

Download references


We are grateful to Dr Daqin Xu at Gansu Provincial Centre for Disease Control and Prevention for her help in collecting the samples.


This work was supported by the National Science and Technology Major Program of China (Nos. 2018ZX10713001-004 and 2018ZX107134-404 to YS, No. 2018ZX10201002-009 to JY), the National Key Research and Development Program of China (2017YFD0501300 to YJ) and the Chinese Special Program for Scientific Research of Public Health (No. 201502021 to JC). The funders had no role in the study design, data collection and analysis, the decision to publish, or preparation of the manuscript.

Author information




YS and JC designed the study. JX, XZ, XW and SC participated in the sample collection and methodology. JX, HJ and JY contributed to data analysis. YS and JC contributed reagents and materials. JX wrote the manuscript. YS and JC revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Jianping Cao or Yujuan Shen.

Ethics declarations

Ethics approval and consent to participate

This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention. The protocol was approved by the Laboratory Animal Welfare & Ethics Committee (LAWEC), National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention (Permit Number: NIPD 2015-011).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Genotypes of E. bieneusi in rodents worldwide.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Xu, J., Wang, X., Jing, H. et al. Identification and genotyping of Enterocytozoon bieneusi in wild Himalayan marmots (Marmota himalayana) and Alashan ground squirrels (Spermophilus alashanicus) in the Qinghai-Tibetan Plateau area (QTPA) of Gansu Province, China. Parasites Vectors 13, 367 (2020).

Download citation


  • Enterocytozoon bieneusi
  • Himalayan marmots
  • Alashan ground squirrels
  • Prevalence
  • Qinghai-Tibetan Plateau area