First detection and genotyping of Enterocytozoon bieneusi in reindeers (Rangifer tarandus): a zoonotic potential of ITS genotypes

Abstract

Background

Enterocytozoon bieneusi is the most common pathogen of 14 microsporidian species infecting humans worldwide. In China, E. bieneusi has been reported in some common livestock and environmental specimens. However, no information is available on occurrence of E. bieneusi in reindeers. The objective of the present study was to detect and genotype E. bieneusi in reindeers in China, and assess the zoonotic potential.

Findings

125 fecal specimens were collected from wild reindeers in the northeast forest region of Great Hinggan Mountains of China. By PCR and sequencing of the internal transcribed spacer (ITS) region of the ribosomal RNA (rRNA) gene of E. bieneusi, an average infection rate of 16.8 % (21/125) was observed in reindeers. E. bieneusi was detected in two age groups: 7.7 % (3/39) in the youths (aged 1 to 2 years) and 22.2 % (18/81) in the adults (aged 3 to 8 years). Five genotypes were identified: one known genotype Peru6 (n = 6) and four novel genotypes named as CHN-RD1 (n = 12), and CHN-RD 2 to CHN-RD4 (one each). In phylogenetic analysis, all the novel genotypes together with known genotype Peru 6 were clustered into group 1.

Conclusions

This is the first report of E. bieneusi infection in reindeers, expanding the host range of E. bieneusi. The fact of genotype Peru 6 previously reported in humans and the result of all the novel genotypes falling into zoonotic group 1 suggest the possibility of E. bieneusi transmitted from reindeers to humans.

Findings

Background

Microsporidia are obligate intracellular eukaryotic pathogens composed of about 1300 species in 160 genera, and they have the ability to infect almost all animal phyla [1]. To date, at least 14 microsporidian species in eight genera have been described as human pathogens. Enterocytozoon bieneusi is the most frequently diagnosed species of microsporidia in humans [2]. Microsporidiosis caused by E. bieneusi is mainly characterized by chronic diarrhea and wasting in HIV-infected patients, but it appears to be asymptomatic or self-limited diarrhea in immunocompetent persons [3]. E. bieneusi is also common inhabitants of the gastrointestinal tract of a wide range of animal hosts, including mammals, birds and reptiles [4, 5].

Application of PCR-based molecular tools for genotyping E. bieneusi has contributed to a better understanding of the characteristics of this pathogen about its host specificity and transmission patterns. Due to the fact of a hypervariable sequence (243 bp) in the internal transcribed spacer (ITS) region of the ribosomal RNA(rRNA) gene within E. bieneusi, sequencing of the ITS gene is the standard method for genotyping E. bieneusi isolates [6]. Molecular data has shown that E. bieneusi is a complex species with multiple genotypes [2]. To date, at least 220 genotypes of E. bieneusi have been described, 64 of which have been detected in humans [2, 7, 8]. 51.56 % (33/64) of human-pathogenic genotypes are also found in animals, supporting a presumption of a zoonotic possibility of E. bieneusi [4, 9]. Phylogenetic analysis indicates that 94 % of the identified ITS genotypes of E. bieneus are in a large group named as human-pathogenic group or group 1, and the remaining ones are clustered into several potentially host-adapted groups named as groups 2 to 8 [10, 11].

At present, molecular epidemiological studies of animal microsporidiosis are mainly from some common livestock, and limited reports involve wild animals [11–15]. In China, reindeers are one of vulnerable animal species and there are at most 700 reindeers alive. The aims of the present study were to determine the natural infection rate of E. bieneusi in reindeers and genotype E. bieneusi isolates by PCR and sequencing of ITS gene as well as assess the potential of zoonotic transmission by phylogenetic analysis.

Methods

Ethics statement

The present study was carried out in accordance with the Law of the People’s Republic of China on the Protection of Wildlife of 1989. The research protocol was reviewed and approved by the Research Ethics Committee of Henan Agriculture University. Before collecting fecal specimens, we contacted the managers of reindeers and obtained their permission to have their animals involved. No animals were injured during this procedure.

Collection of fecal specimens

In October 2014, approximately 10 g of fecal specimens was collected from 125 wild reindeers in captivity from three farms in the northeast forest region of Great Hinggan Mountains (51°10′N, 121°74′E) (Table 1). All the fecal specimens were collected from the ground immediately after defecation by using a sterile disposable latex glove and then placed in a labeled sterile bag individually. To avoid duplicate sampling of animals, each individual was identified according to their ear tag, neck rope, body characteristics such as color and size. All the specimens were transported to our laboratory in a cooler with ice packs within 24 h and stored at 4 °C prior to being used in molecular biological characterizations. Meanwhile, the source, age, gender and health status of each animal was recorded at the time of sampling. The ages of the reindeers ranged from one to 11 years, with all of them showing no clinical signs of illness.

Table 1 Prevalence and distribution of E. bieneusi genotypes by geography in China

DNA extraction

Each fecal specimen was homogenized in distilled water, filtered through gauze and centrifuged at 1500 g for 10 min at room temperature, followed by a wash in distilled water three times. Genomic DNA was extracted from 200 mg of each processed specimen using a QIAamp DNA Mini Stool Kit (Qiagen, Hilden, Germany) according to the manufacturer’s recommended procedures. The eluted DNA was stored at −20 °C until its analysis with PCR.

PCR amplification

All the DNA preparations were detected for the presence of E. bieneusi by nested PCR amplification of a 389 bp nucleotide fragment of the rRNA gene of E. bieneusi and the primers and the cycling parameters in nested PCR analysis were used as previously described by Buckholt et al. [16]. TaKaRa Taq DNA Polymerase (TaKaRa Bio Inc., Tokyo, Japan) was used for all the PCR amplifications. A negative control with no DNA added was included in all PCR tests. All the secondary PCR products were subjected to electrophoresis in a 1.5 % agarose gel and visualized by staining the gel with ethidium bromide.

Nucleotide sequencing and analysis

All the secondary PCR products of expected size were directly sequenced with a set of primers used for the secondary PCR on an ABI PRISM 3730 XL DNA Analyzer by Sinogeno-max Biotechnology Co., Ltd. (Beijing, China), using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA). Sequence accuracy was confirmed by two-directional sequencing and by sequencing additional PCR products if necessary for some DNA preparations.

ITS gene sequences obtained in the present study were aligned with each other and reference sequences downloaded from GenBank database using the Basic Local Alignment Search Tool (BLAST) (http://blast.ncbi.nlm.nih.gov/Blast.cgi) and Clustal X 1.83 (http://www.clustal.org/) to determine E. bieneusi genotypes. If the sequences obtained were identical to those published in GenBank, they were considered to be known genotypes and given the first published name. If not, they were considered to be novel genotypes. All the genotypes were named based on 243 base pairs of the ITS gene region of E. bieneusi according to the established nomenclature system [6].

Phylogenetic analysis

To present the diversity of all the genotypes obtained in the present study and to assess the genetic relationship of novel ones here to known ones, a comparison of the ITS region of all the sequences obtained here and reference sequences published in previous studies was made using the software Mega 5 (http://www.megasoftware.net/) by constructing a neighbor-joining tree, based on the evolutionary distances calculated by a Kimura 2-parameter model. The reliability of these trees was assessed using bootstrap analysis with 1000 replicates.

Results and discussions

In the present study, 16.8 % of 125 reindeers were found to be infected with E. bieneusi, with the highest infection rate in Alongshan Farm (19.5 %; 8/41), followed by 16.9 % (10/59) in Genhe Farm and 12.0 % (3/25) in Jinhe Farm (Table 1). E. bieneusi was only found in two age groups, youths aged 1–2 years (7.7 %, 3/39) and adults aged 3–8 years (22.2 %, 18/81), while infection rate of 9.6 % (5/52) in males and 21.9 % (16/73) in females (Table 2). In general, the present average infection rate was higher than that of the white-tail deer in New York City, the USA (12.2 %, 6/49) [11], but was lower than those in the white-tail deer (32.5 %, 26/80) in Maryland, the USA [15], and in sika deer 32.6 % (28/86) and in red deer 20 % (1/5) in China [9].

Table 2 Prevalence and distribution of E. bieneusi genotypes by age and gender

Five ITS genotypes of E. bieneusi were obtained in the three farms, including one known genotype Peru 6 (syn. PtEbI, PtEbVII) (n = 6) and four novel genotypes (GenBank: KR632538-KR632541) named as CHN-RD1 (n = 12) and CHN-RD2 to CHN-RD4 (one each) (Table 2). 13 nucleotide polymorphic sites were observed among them (Fig. 1). In fact, a high degree of genetic variability of deer-derived E. bieneusi isolates has been described in the ITS gene region of the rRNA gene [9, 11, 15]. To date, the number of genotypes of E. bieneusi in deer has increased to 29, including 19 genotypes (WL4, WL18, WL19, I, J, LW1 and DeerEb1-DeerEb13) in white-tail deer [11, 15], five genotypes (BEB6 and HLJD-I to HLJD-IV) in sika deer and one genotype (HLJD-V) in a red deer in China [9].

Fig. 1

Sequence variation in the ITS region of the rRNA gene of E. bieneusi isolates. The ITS sequences of the known genotype Peru6 and four novel genotypes (CHN-RD1 to CHN-RD4) identified in this study were aligned with one another. Dots indicate the same base identity as the ITS gene sequence of genotype CHN-RD4

Genotype Peru 6 has been found in humans in Peru and Portugal [17–20], in mammals (cattle, sheep, goats, dogs) in USA, China and Portugal [16, 21, 22] and in birds (pigeons and a lovebird) in Portugal [23] as well as in wastewater in China [24]. The findings above suggest the reindeers infected with genotype Peru 6 may pose a threat to humans and other susceptible animal hosts. In phylogenetic analysis, all the novel genotypes belonged to group 1 previously described as a zoonotic group. Genotype CHN-RD1 was sub-clustered into 1e while genotypes CHN-RD2 to CHN-RD4 fell into 1b together with Peru 6 (Fig. 2), suggesting the potential for zoonotic transmission.

Fig. 2

Phylogenetic relationship of Enterocytozoon bieneusi genotype groups. The relationship of E. bieneusi genotypes identified in the present study and other known genotypes deposited in GenBank was inferred by a neighbor-joining analysis of ITS sequences based on genetic distance by the Kimura two-parameter model. The numbers on the branches are percent bootstrapping values from 1000 replicates. Each sequence is identified by its accession number, host origin, and genotype designation. The group terminology for the clusters is based on the work of Karim et al. [25]. The solid and open circles indicate novel and known genotypes identified in this study, respectively

Conclusion

This is the first report of E. bieneusi infection in reindeers, expanding the host range of E. bieneusi. The fact of genotype Peru 6 reported previously in humans and the finding of all the novel genotypes falling into zoonotic group 1 suggest the possibility of zoonotic transmission of E. bieneusi from reindeers to humans. However, due to the limited geographical distribution and the small population of reindeers as well as the fewer infected animals in China, public health significance of infected reindeers is relatively low. In spite of this, advice should be given to those people having close contact with reindeers.