Multilocus sequence typing of Enterocytozoon bieneusi in crab-eating macaques (Macaca fascicularis) in Hainan, China
Parasites & Vectors volume 13, Article number: 182 (2020)
Enterocytozoon bieneusi is one of common intestinal pathogens in humans and animals including non-human primates (NHPs). Many zoonotic pathogens including E. bieneusi have been found in these animals. However, there are few studies on the population structure of E. bieneusi in NHPs. To infer the gene diversity and population genetics of E. bieneusi, we selected 88 E. bieneusi-positive samples from crab-eating macaques for multilocus characterizations in this study.
The E. bieneusi isolates examined belonged to three common genotypes with different host ranges by sequence analysis of the ribosomal internal transcribed spacer (ITS): Type IV (n = 44), Macaque3 (n = 24) and Peru8 (n = 20). They were further characterized by sequence analysis at four microsatellite and minisatellite loci (MS1, MS3, MS4 and MS7). DnaSP, Arlequin and LIAN were used to analyze the sequence data together with those from the ITS locus to infer the population genetics. Subpopulation structure was inferred using phylogenetic and STRUCTURE analyses.
Seventy-two (81.8%), 71 (80.7%), 76 (86.4%) and 79 (89.8%) samples were amplified and sequenced successfully at the MS1, MS3, MS4 and MS7 loci, respectively, with 53 having sequence data at all five MLST loci including ITS. Altogether, 33 multilocus genotypes (MLGs) were produced based on concatenated sequences from the 53 samples. In phylogenetic analyses of sequences and allelic data, four major subpopulations (SPs) were observed with different ITS genotypes in each of them: Type IV and Peru8 in SP1 and SP2; Type IV, Macaque3 and Peru8 in SP3; and Type IV and Macaque3 in SP4. SP3 and SP4 were phylogenetically related and might be NHP-specific based on the fact that Macaque3 is mostly found in NHPs. A strong linkage disequilibrium (LD) was observed among the multilocus sequences and allelic data.
The significant LD in the multilocus sequence analysis indicated the presence of an overall clonal population structure of E. bieneusi in crab-eating macaques. The inconsistent segregation of MLGs among ITS genotypes suggested some occurrence of genetic recombination. These observations should improve our understanding of the population genetics of E. bieneusi in NHPs.
Enterocytozoon bieneusi is an obligate intracellular pathogen infecting humans, domestic animals and wildlife . Since the detection of E. bieneusi in AIDS patients in 1985 for the first time, this pathogen has been found in other immunocompromised/immunodeficient individuals, such as organ transplant recipients, children, and life-threatening and persistent diarrhea often occur in these populations [2,3,4,5,6,7]. It has been recently identified in several species of non-human primates (NHPs) in China and elsewhere [8,9,10,11,12,13,14].
Approximately 500 E. bieneusi genotypes have been identified based on sequence analysis of the internal transcribed spacer (ITS) region of the rRNA gene [15, 16]. They belong to 11 phylogenetic groups. Among them, ITS genotypes in Group 1 and Group 2 have been found in a broad range of hosts including humans and are probably responsible for most zoonotic or cross-species E. bieneusi infections whereas host adaptation seems to be more common in genotypes of Groups 3 to 11 according to a recent review of E. bieneusi . To date, more than 100 E. bieneusi genotypes have been detected in NHPs, mostly belonging to Group 1 [6, 9, 10, 13, 14, 17,18,19,20,21]. Among them, Type IV and Peru8 are common ITS genotypes in humans and various animals in many countries, including China . Conversely, Macaque3 has been mostly found in monkeys and occasionally in dogs in China [1, 21, 23,24,25,26].
Multilocus sequence typing (MLST) based on sequence analysis of four markers (MS1, MS3, MS4 and MS7) with simple tandem repeats and the ITS has been used in investigating the transmission dynamics of E. bieneusi in humans and some animals, including pigs, sheep, rabbits, calves, black bears, giant pandas and horses [27,28,29,30,31,32,33,34,35]. Thus far, only three studies have used the MLST analysis to characterize E. bieneusi isolates in NHPs, and only one study has found the presence of an overall clonal population structure and an epidemic population structure in subpopulations of this pathogen in NHPs [11, 36, 37].
In this study, 88 E. bieneusi-positive samples from crab-eating macaques (Macaca fascicularis) in Hainan Province, China were used for MLST analysis and population genetic characterization of E. bieneusi. The present study aims at evaluating the genetic diversity and population structure of E. bieneusi in NHPs to better understand the transmission and epidemiology of E. bieneusi in these animals.
Source of E. bieneusi-positive samples
The 88 E. bieneusi-positive samples used in this study were from crab-eating macaques in Hainan, China. They consisted of ITS genotypes Type IV (n = 44), Macaque3 (n = 24) and Peru8 (n = 20), which were identified by sequence analysis of the ITS region of the rRNA gene in a previous study . They were collected from different animals and used for further molecular characterization in this study. Other minor E. bieneusi genotypes from that study were not included in the present study, as the small numbers of samples for each of them made examinations of intra-genotypic variations of them difficult.
MLST analysis of E. bieneusi
Four established MLST loci were targeted in this study, including three microsatellite loci (MS1, MS3 and MS7) and one minisatellite locus (MS4). These loci were amplified by nested PCR using published primers and conditions . A positive control (DNA of genotype BEB4 which usually infects cattle rather than NHPs) and a negative control (reagent water without DNA) were used in all PCR tests. PCR was performed twice on each DNA sample at each genetic locus. The secondary PCR products generated were bi-directionally sequenced on an ABI 3730 Genetic Analyzer (Applied Biosystems, Foster City, CA). The nucleotide sequences obtained were assembled using software ChromasPro 1.32 (http://technelysium.com.au/ChromasPro.html), and aligned with each other and homologous sequences downloaded from GenBank using ClustalX (http://clustal.org). The sequence alignments at the four genetic loci (MS1, MS3, MS4 and MS7) were used in the assessment of sequence polymorphism within the ITS genotypes, determination of MLGs and phylogenetic analyses of the MLST data.
Population genetic analysis of sequence data
The multilocus sequences data (concatenated sequences from all five markers, including the ITS) generated were used in the analysis of the population genetics of E. bieneusi using the software DnaSP version 5.10.01 (http://www.ub.edu/dnasp/) . The population structure of E. bieneusi was assessed by intragenic and intergenic linkage disequilibrium (LD), neutrality (Fu’s statistics or Fs values on selective neutrality based on genotype frequency), and recombination events (Rms).
Population genetic analysis of allelic profile data
Pairwise intergenic LD based on allelic profiles at the five genetic loci was evaluated using the exact test and Markov chain parameters implemented in Arlequin version 22.214.171.124 . The standardized index of association (ISA) was measured using LIAN 3.7 .
Subpopulation structure analysis
The phylogenetic relationship among MLGs including both SNPs and INDELs was assessed using the maximum likelihood (ML) analysis and general time reversible model implemented in MEGA 6 (https://www.megasoftware.net). The reliability of clusters formed was assessed using bootstrap analysis with 1000 replicates. An un-weighted pair group mean average (UPGMA) tree of the allelic data was generated using MEGA 6. Wright’s fixation index (FST) of E. bieneusi isolates was calculated using software Arlequin version 126.96.36.199.
STRUCTURE version 2.3.1 was used to identify distinct sub-populations, with K-values of 2–4 . FST was calculated to assess the robustness of the sub-structuring and to determine genetic distances between subpopulations.
PCR amplification efficiency at the MLST loci
The PCR amplification efficiencies were 81.8%, 80.7%, 86.4% and 89.8% at the MS1, MS3, MS4 and MS7 loci, respectively (Table 1). For ITS genotypes Type IV and Peru8, the amplification efficiencies were similar at the MS1 (88.7% and 80.0%), MS3 (81.8% and 85.0%), MS4 (86.4% and 85.0%) and MS7 (91.2% and 95.0%) loci. For Macaque3, however, the amplification efficiencies at the MS1 (70.8%), MS3 (75.0%), MS4 (87.5%) and MS7 (79.2%) loci were correspondingly lower than those of ITS genotypes Type IV and Peru8 at the three loci except MS4.
Genetic diversity at the MLST loci was in the form of both SNPs and indels, with 24, 5, 25, 4 and 3 segregation sites at the MS1, MS3, MS4, MS7 and ITS loci, respectively (Table 2). Altogether, 6, 4, 6, 2 genotypes were detected at the MS1, MS3, MS4 and MS7 loci, respectively. The intragenic LD among segregating sites for each locus was calculated using the linear regression equation. The MS3, MS7 and ITS loci were in complete linkage disequilibrium (LD = 1). In contrast, the remaining two loci (MS1 and MS4) had incomplete LD, indicating the potential presence of genetic recombination. This was confirmed by further assessment of intragenic recombination, which showed the presence of two and five Rms at the MS1 and MS4 loci, respectively.
Of the 88 samples analyzed by MLST, 53 were positive at all four genetic loci. To determine the MLGs, sequences of the five genetic loci (MS1, MS3, MS4, MS7 and ITS) from each of the 53 samples were concatenated. Genetic diversity was assessed based on both total polymorphic sites (finite population variance estimates) and segregating sites only (infinite population variance estimates). Together they formed 33 MLGs with an Hd of 0.98 (Table 3). The mean numbers of pairwise differences in polymorphic and segregating sites were 16 and 14, respectively. The frequency of the MLGs ranged from 9.4% (one MLG with five samples), 7.5% (four MLGs each with three samples), 3.8% (eight MLGs each with two samples), to 1.9% (20 MLGs each with one sample) (Table 4).
Inter-locus LD over all segregating sites within the multilocus contig was assessed for all 53 samples. This produced a Zns value of 0.185; the probability P-values for expected Zns ≤ 0.185 was 0.976. The incomplete LD among these 53 samples (|D′| Y = 0.8679 − 0.2781X) indicated a possible occurrence of genetic recombination, which was confirmed by results of the recombination analysis (10 Rms) (Table 3). Thereafter, we used Fu’s statistic test to evaluate the genotype frequency in this population. In this test, we obtained Fs values of -5.56 (P (Fs < − 5.56) = 0.158) and − 6.91 (P (Fs < − 6.91) = 0.122) based on polymorphic sites and segregating sites, respectively (Table 3). The results of Fu’s statistic tests indicated the presence of an excessive number of alleles.
We calculated ISA to assess the LD between alleles in pairwise combinations of genetic loci . Results of this analysis of allelic profile data showed a positive ISA value together with VD > L. Furthermore, a significant PMC value (< 0.001) was obtained using the Monte Carlo method based on the null hypothesis of panmixia (Table 5). All these results indicated that the population was in LD. However, values of VD < L and PMC > 0.001 were obtained for subpopulations identified in this study or ITS genotype. Thus, these subpopulations and ITS genotypes were in linkage equilibrium (LE) and had the epidemic population structure.
Maximum likelihood analysis of the concatenated sequences of the 33 MLGs showed four main subpopulations (SPs). SP1 contained samples of three MLGs of ITS genotype Type IV and two MLGs of ITS genotype Peru8, SP2 contained six MLGs of ITS genotype Type IV and two MLGs of ITS genotype Peru8, SP3 contained seven MLGs of ITS genotype Type IV, four MLGs of ITS genotype Peru8, and five MLGs of ITS genotype Macaque3, while SP4 contained two MLGs of ITS genotype Type IV and two MLGs of ITS genotype Macaque3 (Fig. 1). Similarly, four major clades were formed in UPGMA analysis of 33 MLGs in agreement with results of ML analysis of the concatenated sequence data (Fig. 2).
In STRUCTURE analysis of the MLST data, we used the smallest value of K (K = 2) to measure population differentiation within the E. bieneusi population. Two subpopulations (SPs) with alpha (α) value of 0.104 were produced (1, 2; Fig. 3a). SP1 contained 14 MLGs, while SP2 contained the remaining 19 MLGs (Fig. 3a). When the K value was increased to 3, three clear SPs were formed (α = 0.074; Fig. 3b). SP1 included 8 MLGs belonging to SP1 at K = 2, SP2 included 3 MLGs belonging to SP1 and 7 MLGs belonging to SP2 at K = 2, while SP3 included 3 MLGs belonging to SP1 and 12 MLGs belonging to SP2 at K = 2 (Fig. 3a, b). When using K = 4, four clear SPs formed, with individuals in SP1 being similar to SP1 at K = 2 and 3 (Fig. 3). In addition, the α value of 0.061 was similar to the values generated at K = 2 and 3, suggesting the substructures at K = 4 could be realistic.
FST analysis was used to assess the sub-structuring within the E. bieneusi population. For SP1 and SP2 generated by STRUCTURE at K = 2, the differentiation was significant (FST = 0.161, P < 0.0001). At K = 3, the differentiation between SP1 and SP2 (FST = 0.362, P < 0.0001), SP1 and SP3 (FST = 0.274, P < 0.0001), and SP2 and SP3 (FST = 0.234, P < 0.0001) were also significant. Similarly at K = 4, the differentiation between SP1 and SP2 (FST = 0.349, P < 0.0001), SP1 and SP3 (FST = 0.345, P < 0.0001), SP1 and SP4 (FST = 0.313, P < 0.0001), SP2 and SP3 (FST = 0.293, P < 0.0001), SP2 and SP4 (FST = 0.233, P < 0.0001), and SP3 and SP4 (FST = 0.254, P < 0.0001) were all significant.
Currently, multilocus sequence typing has been used to analyze the genetic diversity of E. bieneusi around the world [15, 42, 43]. In Peru, India and Nigeria, 66 MLGs were identified in HIV-infected individuals [33, 37]. In China, 59 and 44 MLGs were identified in NHPs and pigs, respectively, while some genetic diversity based on MLST data was also observed in E. bieneusi isolates from domestic, wild, and zoo animals [29, 36, 42,43,44,45,46,47,48]. However, few studies analyzed the MLGs of this pathogen in NHPs in China. In the present study, 33 MLGs were identified from 53 E. bieneusi-positive samples from macaque monkeys, indicating a high genetic diversity within some common ITS genotypes (Type IV, Macaques3 and Peru8), and supporting high-resolution of the multilocus sequencing tool established by Feng et al. .
A strong and significant intragenic LD was observed in the analysis of the MLST data, indicating an overall clonal population structure of E. bieneusi in this study. In previous studies, the strong and significant intragenic LD obtained was used as an indication for the presence of a clonal population structure of E. bieneusi in the AIDS patients, giant pandas, pigs, fur animals and NHPs [36, 37, 47, 49, 50]. Therefore, there was a clonal endemicity of E. bieneusi in the study population.
Nevertheless, the MLGs found in crab-eating macaques were not segregated by ITS genotypes, indicating the existence of genetic recombination in the E. bieneusi population studied. In previous investigations, E. bieneusi MLGs were segregated by ITS genotypes for some genotypes with narrower host ranges, such as A, EbpA, SCO2 and Cs-4, but not for zoonotic genotypes such as Type IV, D, WL11 and Peru11 [15, 16]. Similarly, in this study, MLGs of E. bieneusi in crab-eating macaques were also not segregated by ITS genotypes, especially for the zoonotic genotypes Type IV and Peru8. Genetic recombination appears to be common in some microsporidia of insects, including Kneallhazia solenopsae in fire ants and Andreanna caspii in mosquitoes. In a previous study in China, the MLGs of some common ITS genotypes (Type IV, Peru8, D, Macaque3, CM2, Peru11, Henan V, PigEBITS7) in NHPs were not segregated by ITS genotypes in phylogenetic analysis, suggesting a possible occurrence of sexual recombination in E. bieneusi . Therefore, data from both studies have indicated that there could be some sexual recombination among common E. bieneusi genotypes in NHPs.
Four SPs were obtained in this study among MLGs isolates within Group 1, according to results of phylogenetic and subpopulation genetics analyses. The four SPs all appear to have an epidemic population structure. In previous studies, the ITS genotype Macaque3 was mostly found in NHPs, and only occasionally in dogs [8,9,10, 26, 51]. Therefore, the clustering of the host-adapted ITS genotype Macaque3 in SP3 and SP4 suggests the NHP-specific nature of these two subpopulations. Both potentially zoonotic and host-adapted subpopulations have been identified in Group 1 in previous studies. The former included MLGs of ITS genotypes Type IV, D, Peru8, WL11, and Peru11 in humans, NHPs, pigs and companion animals, whereas the latter included MLGs of ITS genotypes A in humans, EbpA and CS-4 in pigs, and Macaque3 in NHPs [1, 15, 21, 36, 52, 53]. Therefore, SP3 and SP4 including MLGs of Macaque3 were potentially NHPs-adapted subpopulations, while SP1 and SP2 only including MLGs of ITS genotype Type IV and Peru8 were potentially zoonotic subpopulations.
A much higher genetic diversity was revealed in E. bieneusi isolates from crab-eating macaques by multilocus sequence typing (33 MLGs) than by ITS genotyping (three ITS genotypes). This appears to have resulted from genetic recombination among common ITS genotypes despite an overall clonal population structure in E. bieneusi. Nevertheless, the MLGs generated from the analysis have formed four subpopulations with different host ranges and an epidemic population structure. These results provide data on the population and subpopulation structure of E. bieneusi in NHPs, and would improve our perception of E. bieneusi epidemiology in these animals.
Availability of data and materials
The data supporting the conclusions of this article are included within the article. Representative nucleotide sequences generated in this study were deposited in the GenBank database under the accession numbers MT137078-MT137095.
internal transcribed spacer
polymerase chain reaction
un-weighted pair group mean average
Santin 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.
Tavalla M, Mardani-Kateki M, Abdizadeh R, Nashibi R, Rafie A, Khademvatan S. Molecular identification of Enterocytozoon bieneusi and Encephalitozoon spp. in immunodeficient patients in Ahvaz, Southwest of Iran. Acta Trop. 2017;172:107–12.
Didier ES, Weiss LM. Microsporidiosis: not just in AIDS patients. Curr Opin Infect Dis. 2011;24:490–5.
Lobo ML, Xiao L, Antunes F, Matos O. Microsporidia as emerging pathogens and the implication for public health: a 10-year study on HIV-positive and -negative patients. Int J Parasitol. 2012;42:197–205.
Li W, Feng Y, Santin M. Host specificity of Enterocytozoon bieneusi and public health implications. Trends Parasitol. 2019;35:436–51.
Matos O, Lobo ML, Xiao L. Epidemiology of Enterocytozoon bieneusi infection in humans. J Parasitol Res. 2012;2012:981424.
Chen L, Zhao J, Li N, Guo Y, Feng Y, Feng Y, et al. Genotypes and public health potential of Enterocytozoon bieneusi and Giardia duodenalis in crab-eating macaques. Parasites Vectors. 2019;12:254.
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.
Yang H, Lin Y, Li Y, Song M, Lu Y, Li W. Molecular characterization of Enterocytozoon bieneusi isolates in laboratory macaques in north China: zoonotic concerns. Parasitol Res. 2017;116:2877–82.
Zhong Z, Li W, Deng L, Song Y, Wu K, Tian Y, et al. Multilocus genotyping of Enterocytozoon bieneusi derived from nonhuman primates in southwest China. PLoS ONE. 2017;12:e0176926.
Mynarova A, Foitova I, Kvac M, Kvetonova D, Rost M, Morrogh-Bernard H, et al. Prevalence of Cryptosporidium spp., Enterocytozoon bieneusi, Encephalitozoon spp. and Giardia intestinalis in wild, semi-wild and captive orangutans (Pongo abelii and Pongo pygmaeus) on Sumatra and Borneo, Indonesia. PLoS ONE. 2016;11:e0152771.
Sak B, Petrzelkova KJ, Kvetonova D, Mynarova A, Pomajbikova K, Modry D, et al. Diversity of microsporidia, Cryptosporidium and Giardia in mountain gorillas (Gorilla beringei beringei) in Volcanoes National Park, Rwanda. PLoS ONE. 2014;9:e109751.
Li W, Kiulia NM, Mwenda JM, Nyachieo A, Taylor MB, Zhang X, et al. Cyclospora papionis, Cryptosporidium hominis, and human-pathogenic Enterocytozoon bieneusi in captive baboons in Kenya. J Clin Microbiol. 2011;49:4326–9.
Li W, Xiao L. Multilocus sequence typing and population genetic analysis of Enterocytozoon bieneusi: host specificity and its impacts on public health. Front Genet. 2019;10:307.
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.
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.
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.
Du S, Zhao G, Shao J, Fang Y, Tian G, Zhang L, et al. Cryptosporidium spp., Giardia intestinalis, and Enterocytozoon bieneusi in captive non-human primates in Qinling mountains. Korean J Parasitol. 2015;53:395–402.
Sak B, Petrzelkova KJ, Kvetonova D, Mynarova A, Shutt KA, Pomajbikova K, et al. Long-term monitoring of microsporidia, Cryptosporidium and Giardia infections in western lowland gorillas (Gorilla gorilla gorilla) at different stages of habituation in Dzanga Sangha Protected Areas, Central African Republic. PLoS ONE. 2013;8:e71840.
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.
Wang S, Wang R, Fan X, Liu T, Zhang L, Zhao G. Prevalence and genotypes of Enterocytozoon bieneusi in China. Acta Trop. 2018;183:142–52.
Lores B, Lopez-Miragaya I, Arias C, Fenoy S, Torres J, del Aguila C. Intestinal microsporidiosis due to Enterocytozoon bieneusi in elderly human immunodeficiency virus-negative patients from Vigo, Spain. Clin Infect Dis. 2002;34:918–21.
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.
Santin M, Fayer R. Enterocytozoon bieneusi genotype nomenclature based on the internal transcribed spacer sequence: a consensus. J Eukaryot Microbiol. 2009;56:34–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.
Luo R, Xiang L, Liu H, Zhong Z, Liu L, Deng L, et al. First report and multilocus genotyping of Enterocytozoon bieneusi from Tibetan pigs in southwestern China. Parasite. 2019;26:24.
Zou Y, Hou J, Li F, Zou F, Lin R, Ma J, et al. Prevalence and genotypes of Enterocytozoon bieneusi in pigs in southern China. Infect Genet Evol. 2018;66:52–6.
Tang C, Cai M, Wang L, Guo Y, Li N, Feng Y, et al. Genetic diversity within dominant Enterocytozoon bieneusi genotypes in pre-weaned calves. Parasites Vectors. 2018;11:170.
Tian G, Zhao G, Du S, Hu X, Wang H, Zhang L, 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.
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.
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.
Li W, Cama V, Feng Y, Gilman RH, Bern C, Zhang X, et al. Population genetic analysis of Enterocytozoon bieneusi in humans. Int J Parasitol. 2012;42:287–93.
Desoubeaux G, Nourrisson C, Moniot M, De Kyvon MA, Bonnin V, De La Bretonniere ME, et al. Genotyping approach for potential common source of Enterocytozoon bieneusi infection in hematology unit. Emerg Infect Dis. 2019;25:1625–31.
Widmer G, Dilo J, Tumwine JK, Tzipori S, Akiyoshi DE. Frequent occurrence of mixed Enterocytozoon bieneusi infections in humans. Appl Environ Microbiol. 2013;79:5357–62.
Karim MR, Wang R, He X, Zhang L, Li J, Rume FI, et al. Multilocus sequence typing of Enterocytozoon bieneusi in nonhuman primates in China. Vet Parasitol. 2014;200:13–23.
Li W, Cama V, Akinbo FO, Ganguly S, Kiulia NM, Zhang X, et al. Multilocus sequence typing of Enterocytozoon bieneusi: lack of geographic segregation and existence of genetically isolated sub-populations. Infect Genet Evol. 2013;14:111–9.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003;19:2496–7.
Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.
Haubold B, Hudson RR. LIAN 3.0: detecting linkage disequilibrium in multilocus data. Linkage Analysis. Bioinformatics. 2000;16:847–8.
Falush D, Stephens M, Pritchard JK. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics. 2003;164:1567–87.
Gui B, Zou Y, Chen Y, Li F, Jin Y, Liu M, et al. Novel genotypes and multilocus genotypes of Enterocytozoon bieneusi in two wild rat species in China: potential for zoonotic transmission. Parasitol Res. 2020;119:283–90.
Zhang N, Wu R, Ji T, Cui L, Cao H, Li D, et al. Molecular detection, multilocus genotyping, and population genetics of Enterocytozoon bieneusi in pigs in southeastern China. J Eukaryot Microbiol. 2020;67:107–14.
Ma Y, Zou Y, Liu Q, Xie S, Li R, Zhu X, et al. Prevalence and multilocus genotypes of Enterocytozoon bieneusi in alpacas (Vicugna pacos) in Shanxi Province, northern China. Parasitol Res. 2019;118:3371–5.
Zhao G, Du S, Wang H, Hu X, Deng M, Yu S, 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.
Deng L, Li W, Zhong Z, Gong C, Liu X, Huang X, et al. Molecular characterization and multilocus genotypes of Enterocytozoon bieneusi among horses in southwestern China. Parasites Vectors. 2016;9:561.
Li W, Wan Q, Yu Q, Yang Y, Tao W, Jiang Y, et al. Genetic variation of mini- and microsatellites and a clonal structure in Enterocytozoon bieneusi population in foxes and raccoon dogs and population differentiation of the parasite between fur animals and humans. Parasitol Res. 2016;115:2899–904.
Wang X, Wang R, Ren G, Yu Z, Zhang L, Zhang S, et al. Multilocus genotyping of Giardia duodenalis and Enterocytozoon bieneusi in dairy and native beef (Qinchuan) calves in Shaanxi Province, northwestern China. Parasitol Res. 2016;115:1355–61.
Li W, Song Y, Zhong Z, Huang X, Wang C, Li C, et al. Population genetics of Enterocytozoon bieneusi in captive giant pandas of China. Parasites Vectors. 2017;10:499.
Wan Q, Xiao L, Zhang X, Li Y, Lu Y, Song M, et al. Clonal Evolution of Enterocytozoon bieneusi populations in swine and genetic differentiation in subpopulations between isolates from swine and humans. PLoS Negl Trop Dis. 2016;10:e0004966.
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.
Thellier M, Breton J. Enterocytozoon bieneusi in human and animals, focus on laboratory identification and molecular epidemiology. Parasite. 2008;15:349–58.
Li D, Zheng S, Zhou C, Karim MR, Wang L, Wang H, et al. Multilocus typing of Enterocytozoon bieneusi in pig reveals the high prevalence, zoonotic potential, host adaptation and geographical segregation in China. J Eukaryot Microbiol. 2019;66:707–18.
We thank the farm owner and staff for their assistance in sample collection during this study.
This study was supported by the National Key R&D Program of China (2017YFC1601206) and the National Natural Science Foundation of China (31820103014).
Ethics approval and consent to participate
This study was approved by the Research Ethics Committee of the East China University of Science and Technology, with the approval number of 2015018 on 3rd January 2015. Permission for collecting fecal specimens was obtained from the managers of the farms. Animals were kept according to the Animal Ethics Procedures and Guidelines of the People’s Republic of China.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, L., Li, N., Guo, Y. et al. Multilocus sequence typing of Enterocytozoon bieneusi in crab-eating macaques (Macaca fascicularis) in Hainan, China. Parasites Vectors 13, 182 (2020). https://doi.org/10.1186/s13071-020-04046-w