Genotypic analysis of zoonotic Enterocytozoon bieneusi in wild deer in Korea
Article information
Abstract
Enterocytozoon bieneusi is an important microsporidian protozoa that causes intestinal disorders in humans. We collected 191 fecal samples from roadkill deer carcasses, among which 13 (6.8%) showed positive reaction for E. bieneusi by polymerase chain reaction assay. Phylogenetic analysis revealed 6 distinct genotypes, 1 of which was novel. All genotypes belonged to Group 1, which has low host specificity, indicating possible transmission through sylvatic cycle. E. bieneusi infection was predominant in female deer (P<0.05).
Enterocytozoon bieneusi, previously classified in the kingdom Fungi, is a zoonotic microsporidia currently grouped into Cryptomycota, a sister group to fungi [1,2]. This opportunistic pathogen causes intestinal disorders, such as chronic diarrhea and intestinal malabsorption, in immunocompromised or immunocompetent individuals, while most infections are asymptomatic [3,4].
Polymerase chain reaction (PCR)-based sequencing analysis of the polymorphic internal transcribed spacer (ITS) region of E. bieneusi has resulted in the identification of over 500 E. bieneusi genotypes in a wide range of host animals, including mammals and birds [2,5]. To date, 11 major groups have been established using phylogenetic analysis: zoonotic genotypes Groups 1 and 2, and host-adapted genotypes Groups 3–11. Thus, E. bieneusi is a critical public health concern [6]. E. bieneusi transmission mainly occurs via the fecal–oral route, beginning with the intake of environmentally-resistant spores that invade host cells via a polar tube that penetrates the cell. There are also reports of airborne and hospital-associated transmission [7–12].
Because of its zoonotic potential, E. bieneusi has been globally studied. However, compared with other animals, studies investigating the prevalence of E. bieneusi infection in deer are scarce, particularly in Korea. Therefore, we investigated the prevalence of E. bieneusi infection in deer in Korea, including the Korean water deer (Hydropotes inermis argyropus), which is a dominant species in Korea, and roe deer (Capreolus capreolus).
We collected 191 fresh fecal samples from Korean water deer or roe deer that were accidentally killed on highways or roadsides from March 2018 to June 2021. When a roadkill carcass is found in Korea, it is usually reported to the government office by citizens passing by on the road, leading to immediate collection by a field dispatch unit. Thus, the carcass is kept in a fresh condition. The fecal samples were directly collected per rectum by trained veterinarians, placed in plastic tubes, and delivered to the laboratory for analysis. Data on the collection region, season, sex, weight, and host species were recorded, with missing information marked as “unknown.” The samples were classified into three regions according to provincial boundaries: northern (Gangwon, Gyeonggi, Incheon, and Seoul), central (Chungbuk, Chungnam, Gyeongbuk, Jeonbuk, Sejong, and Daejeon), and southern (Gyeongnam, Jeonnam, Jeju, Gwangju, Busan, and Ulsan). The study population was categorized into five groups according to weight: <10 kg, ≥10 to <15 kg, ≥15 to <20 kg, ≥20 kg, and unknown. Statistical analysis was performed using the χ2 test in SPSS v26.0 (IBM Corporation, Armonk, NY, USA). P-values <0.05 were considered statistically significant. Approval number of dead wild animals was not needed.
DNA was extracted from the fecal samples using a QIAamp Fast DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. The quality and quantity of the extracted DNA were determined spectrophotometrically using an Infinite 200 PRO NanoQuant plate reader (Tecan, Männedorf, Switzerland) before storage at −20°C until use. The presence of E. bieneusi was determined using nested PCR targeting the ITS region, with the primer pair ITSF1 (forward, 5′-GGTCATAGGGATGAAGAG-3′) and ITSR1 (reverse, 5′-TTCGAGTTCTTTCGCGCTC-3′) and the primer pair ITSF2 (forward, 5′-GCTCTGAATATCTATGGCT-3′) and ITSR2 (reverse, 5′-ATCGCCGACGGATCCAAGTG-3′) used in the first and second PCR rounds, respectively. The expected amplicon size was 389 bp [12]. Each PCR assay was performed in a final reaction volume of 20 μl comprising 1 μl of each primer, 2 μl template DNA, and 16 μl distilled H2O in the lyophilized AccuPower HotStart PCR Premix (Bioneer, Daejeon, Korea) using a MastercyclerPro thermal cycler (Eppendorf, Hamburg, Germany). The cycling conditions were as follows: initial denaturation at 95°C for 5 min; 35 cycles of denaturation at 95°C for 30 sec, annealing at 55°C for 30 sec, and extension at 72°C for 30 sec; and final extension at 75°C for 5 min. The resultant PCR amplicons were electrophoresed on a 1% agarose gel and visualized with ethidium bromide staining. All PCR-positive samples underwent bidirectional sequencing (Macrogen, Daejeon, Korea). The obtained sequences were aligned using BioEdit v7.2.5 (https://thalljiscience.github.io/page2.html, accessed 10 Oct. 2024), revealing 13 consensus sequences, of which 6 were identical and consequently removed. The remaining 7 sequences were submitted to GenBank (accession numbers: OR335887–OR335893) and compared with an E. bieneusi reference sequence (DQ885585) in NCBI BLAST. A maximum likelihood phylogenetic tree with 1,000 bootstrap replicates was generated using MEGA 7 software [13].
The E. bieneusi infection rate in our study was 6.8% (13/191 samples) (Table 1). E. bieneusi infection has been studied in various species, including bats, wild boars, domestic pigs, dairy cattle, primates in zoos, lambs, and other wildlife species [2,5,7,14–20]. Studies on E. bieneusi infection in deer have been performed in several countries, mainly in China [16,20]. However, studies on E. bieneusi infection in water deer and roe deer are scarce. This study compared the E. bieneusi infection rate with other Cervidae species, including sambar deer (Rusa unicolor), red deer (Cervus elaphus), fallow deer (Dama dama), and sika deer (Cervus nippon). The rate of E. bieneusi infection was reported as 4.1% (25/610 samples) in an Australian population including sambar deer, red deer, and fallow deer [17], which was slightly lower than in our study. However, the infection rate of E. bieneusi in wild red deer in Spain was reported as much lower (1.5%, 5/329 samples) [20] and the infection rate in sika deer and red deer in China was reported to be much higher (35.9%, 221/615 samples) [16] than in our study. Studies on E. bieneusi infection in bats, wild boars, calves, Korean water deer, raccoon dogs, and dairy cattle have been conducted in Korea, and the infection rate in Korean water deer was reported as 53.6% (52/97) [15], considerably higher compared with that in our study. This result demonstrated that there may be various differences depending on the methodology, study region, and sampling time.
We also analyzed the risk factors predisposing to E. bieneusi infection, such as locality, season, sex, weight, and host species. Female deer had a significantly higher infection rate than males (19.2% vs. 0%, P=0.008), although the sex of eight of the infected deer was unknown. Deer thrives in the southern region had a higher infection rate (13.2%) and in the spring season (9.7%), but not statistically significant (P=0.054 and 0.583, respectively). Other studies have not reported significant differences in the infection prevalence between females and males [16,18,19]. In Chinese deer, the prevalence was 30.1% in females and 42.1% in males [16]; in Ethiopian lambs, the prevalence was 11.0% in females and 9.0% in males [19]; and in Korean wild boars, the prevalence was 3.0% in females and 2.6% in males [18]. Thus, further studies with a larger sample size are necessary to clarify this association.
Regarding weight and species, the infection rate was high in the ≥15 to <20 kg group (15.0%) and in roe deer (10.5%), albeit not statistically significant (P=0.483 and P=0.239, respectively).
Phylogenic analysis revealed 6 distinct E. bieneusi genotypes: Korea-WL2 (OR335891), Korea-WL6 (OR335887 and OR335892), H (OR335888), KWB1 (OR335889), D (OR335893), and RDK (OR335890). RDK (OR335890) was a novel genotype with 3 nucleotide substitutions at positions 200 (T/A), 207 (A/G), and 260 (A/G) compared with the reference genotype PigHN-II sequence (MN699293). Among the 7 samples positive for E. bieneusi, 1 sample was identical to genotype Korea-WL2 (subgroup 1a, identified in raccoon dogs and Korean water deer) and two were identical to genotype Korea-WL6 (subgroup 1e, identified in Korean water deer) [15]; 1 sample each was identical to genotype H (subgroup 1e), KWB1 (subgroup 1e, identified in wild boar), and D (subgroup 1a, identified in wild boar) [18], respectively. Finally, 1 sample was classified as the novel RDK genotype, which belonged to subgroup 1e. All 6 genotypes identified in our study belonged to Group 1 (Fig. 1), the largest genetic group of E. bieneusi. This group is a significant public health concern because of its low host specificity [5,6].
Of the 191 fecal samples tested in our study, 13 (6.8%) were infected with E. bieneusi. There are few studies on E. bieneusi infection in roe deer. Only 1 study on Korean water deer has been reported [15]. Although we found a significant association between sex and E. bieneusi infection, further studies are necessary to verify this relationship. We identified 6 distinct E. bieneusi genotypes, all belonging to Group 1 and 1 of which was novel. E. bieneusi is genetically diverse and infects a wide range of host species, including humans. Korea has large populations of Korean water deer and roe deer. Given its zoonotic potential and wide geographic distribution, E. bieneusi infection poses a significant public health concern in wild life. Further research on the infection rates and the pathogenicity of the novel RDK genotype is warranted.
Notes
Author contributions
Conceptualization: Kwak D
Project administration: Kwak D
Methodology: Noh G, Lee H, Lee SH, Seo MG, Kim KT, Lee J, Nazim K, Park SJ, Rhee MH
Investigation: Noh G, Lee H, Lee SH, Seo MG, Kim KT
Data curation: Lee H, Lee SH, Seo MG, Kim KT, Lee J, Nazim K
Formal analysis: Noh G, Lee H, Lee SH, Seo MG, Kim KT
Visualization: Nazim K, Park SJ, Rhee MH
Validation: Nazim K, Park SJ, Rhee MH
Supervision: Kwak D
Writing – original draft: Noh G, Lee H, Lee SH, Seo MG, Kim KT, Lee J, Nazim K
Writing– review and editing: Noh G, Lee H, Lee SH, Seo MG, Kim KT, Lee J, Nazim K, Park SJ, Rhee MH, Kwak D
The authors declare no conflict of interest related to this study.