Occurrence and genotype of E. bieneusi in red pandas
Among 198 fecal samples detected by nested PCR, 24 tested positives for E. bieneusi, the overall infection rate reached 12.1% (24/198) (Table 1). The infection rates of E. bieneusi were 18.0% (18/100), 7.1% (3/42), 7.1% (1/14), 9.1% (1/11),0% (0/18),7.7% (1/13) in Chengdu Research Base of Giant Panda Breeding(CDBGPB), Sanlang Resort(SLR), Bifengxia Zoo(BFXZ), Chengdu Zoo(CDZ), Wolong China Giant Panda Garden(WCGPG), Panda Valley(PV), respectively. The differences among different regions were not significant (χ2 = 7.364, df = 5, P = 0.1949 > 0.05). In groups, Group i was 5.9% (5/85), and Group ii reached up to 16.8% (19/113). There was a significant difference between in different groups (χ2 = 5.442, df = 1, P = 0.0197 < 0.05).
Table 1
Sampling information and epidemic results of E. bieneusi in red pandas in this study.
Group
|
Region
|
Site
|
Breeding mode
|
No. of samples
|
No. of positive (%)
|
ITS
genotypes
|
Ⅰ
|
Chenghua
|
CDZ
|
Captive
|
11
|
1 (9.1%)
|
D
|
|
Wenchuan
|
WLCGPG
|
Captive
|
18
|
0 (0.0%)
|
|
|
Chongzhou
|
SLR
|
Captive
|
42
|
3 (7.1%)
|
SCR1
|
|
Ya'an
|
BFXZ
|
Captive
|
14
|
1 (7.1%)
|
D
|
Ⅱ
|
Chenghua
|
CDRBGPB
|
Semi-free-range
|
100
|
18 (18.0%)
|
PL2
|
|
Dujangyan
|
PV
|
Semi-free-range
|
13
|
1 (7.7%)
|
SC02
|
Note: CDBGPB: Chengdu Research Base of Giant Panda Breeding; SLR: Sanlang Resort; BFXZ: Bifengxia Zoo; CDZ: Chengdu Zoo; WCGPG: Wolong China Giant Panda Garden; PV: Panda Valley; |
Phylogenetic analysis and MLST genotyping of E. bieneusi
All 24 positive isolates were successfully sequencing. observed Four genotypes, including 3 described genotypes (D, PL2, and SC02) and one novel genotype SCR1. The obtained ITS gene sequences of genotypes D, PL2, and SC02 were identical with reference sequences KY950534, MT497891, and KU852476 from GenBank, respectively. The homology was 98.6% between novel genotype SCR1 and genotype FJL (MK357781), with one single nucleotide substitution and one insertion. Among them, genotype PL2 (n = 18) was dominant and genotype SCR1 (n = 3) was ranked second. Both SC02 and D only identified one sample.
Table 2
The primers and annealing temperatures of PCR amplification.
Gene
|
Primer sequences (5'-3')
|
AT (℃)
|
Size (bp)
|
Reference
|
ITS1
|
F1: GATGGTCATAGGGATGAAGAGCTT
|
55
|
392
|
[31]
|
R1: AATACAGGATCACTTGGATCCGT
|
|
F2: AGGGATGAAGAGCTTCGGCTCTG
|
55
|
|
R2: AATATCCCTAATACAGGATCACT
|
|
ITS2
|
F1: GGTCATAGGGATGAAGAG
|
57
|
390
|
[14]
|
|
R1: TTCGAGTTCTTTCGCGCTC
|
|
|
|
|
F2: GCTCTGAATATCTATGGCT
|
55
|
|
|
|
R2: ATCGCCGACGGATCCAAGTG
|
|
|
|
MS1
|
F1: AAGTTGCAAGTTCAGTGTTTGAA
|
58
|
676
|
[14]
|
R1: GATGAATATGCATCCATTGATGTT
|
|
F2: TTGTAAATCGACCAAATGTGCTAT
|
58
|
|
R2: ACATAAACCACTAATTAATGTAAC
|
|
MS3
|
F1: CAAGCACTGTGGTTACTGTT
|
55
|
537
|
[14]
|
R1:AAGTTA GGGCATTTAATAAAATTA
|
|
F2: GTTCAAGTAATTGATACCAGTCT
|
55
|
|
R2: CTCATTGAATCTAAATGTGTATAA
|
|
MS4
|
F1: GCATATCGTCTCATAGGAACA
|
55
|
885
|
[30]
|
R1: GTTCATGGTTATTAATTCCAGAA
|
|
F2: CGA AGTGTACTACATGTCTCT
|
55
|
|
R2: GGACTTTAATAAGTTACCTATAGT
|
|
MS7
|
F1: GTTGATCGTCCAGATGGAATT
|
55
|
471
|
[14]
|
R1: GACTATCAGTATTACTGATTATAT
|
|
F2: CAATAGTAAAGGAAGATGGTCA
|
55
|
|
R2: CGTCGCTTTGTTTCATAATCTT
|
|
Note: AT: Annealing temperatures; F1:Forward primer of primary PCR; R1: Reverse primer of the primary PCR; F2:Forward primer of the second PCR; R2:Reverse primer of the second PCR. |
Phylogenetic analysis results (Fig. 2) showed that genotype D and SCR1 belonged to subgroup1a in Group 1. SC02 is clustered in subgroup 1b of Group 1, and PL2 formed a separate clade related to Group 2 [33]. All genotypes were human-pathogenic [31]. To further analyze the gene polymorphism and the sub-genotype by using the MLST (Table 3). The amplification efficiencies in loci MS1, MS3, MS4, and MS7 were 25.0%(6/24), 45.8%(11/24),16.7%(4/24), and 50.0%(12/24), respectively. Only 3 genotype PL2 specimens tested positive in these five loci. Correspondingly generated 1, 2, 4, and 1 types in MS1, MS3, MS4, and MS7 locus based on different nucleotide repeats and SNPs. 3 distinct MLGs were formed in genotype PL2, and some degree of sequence polymorphism was shown.
Table 3
Multilocus genotypes of positive E.bieneusi isolates from red pandas in this study.
Sampling sites
|
Sampling
no.
|
ITS genotype
|
Multilocus genotypes
|
MS1
|
MS3
|
MS4
|
MS7
|
MLGs
|
CDRBGPB
|
3
|
PL2
|
|
|
|
|
|
|
5
|
PL2
|
|
type1
|
|
type1
|
|
|
10
|
PL2
|
|
type1
|
|
type1
|
|
|
11
|
PL2
|
type1
|
type1
|
type1
|
type1
|
MLG1
|
|
12
|
PL2
|
type1
|
type1
|
|
|
|
|
15
|
PL2
|
|
|
|
|
|
|
16
|
PL2
|
|
|
|
type1
|
|
|
18
|
PL2
|
|
|
|
|
|
|
19
|
PL2
|
type1
|
|
type1
|
type1
|
|
|
21
|
PL2
|
|
|
|
|
|
|
22
|
PL2
|
|
type1
|
|
|
|
|
40
|
PL2
|
|
|
|
|
|
|
43
|
PL2
|
|
|
|
|
|
|
47
|
PL2
|
|
type1
|
|
type1
|
|
|
63
|
PL2
|
|
type1
|
|
type1
|
|
|
86
|
PL2
|
type1
|
type1
|
type2
|
type1
|
MLG2
|
|
93
|
PL2
|
type1
|
|
|
type1
|
|
|
95
|
PL2
|
type1
|
type1
|
type3
|
type1
|
MLG3
|
SLR
|
122
|
SCR1
|
|
|
|
|
|
SLR
|
136
|
SCR1
|
|
type1
|
|
type1
|
|
SLR
|
145
|
SCR1
|
|
|
|
|
|
BFXZ
|
174
|
D
|
|
type2
|
|
|
|
PV
|
185
|
SC02
|
|
|
|
type1
|
|
CDZ
|
211
|
D
|
|
|
|
|
|
Note: CDBGPB: Chengdu Research Base of Giant Panda Breeding; SLR: Sanlang Resort; BFXZ: Bifengxia Zoo; CDZ: Chengdu Zoo; WCGPG: Wolong China Giant Panda Garden; PV: Panda Valley; |
This study sampled 6 zoos in 5 different regions and demonstrated the prevalence and genotype characterization of E. bieneusi among captive red pandas in Sichuan province, southwest China. In this study, red pandas in Group ii have a significantly higher prevalence of infection than Group i. The difference may be related to the larger range of red pandas and the more frequent contact between different animal populations and between animals and humans under the semi-free-range breeding mode. Meanwhile, higher prevalence indicated that the safety of semi-free-range breeding mode needs to be carefully regulated, due to zoonotic transmission potential.
The overall prevalence of this study was 12.1%, which is close to the previous report of red pandas in Shaanxi province, northwestern China, 11.1%-13.9% [13, 28]. In the present study, the prevalence of Chengdu Zoo and Bifengxia Zoo was 7.1% (1/14) and 9.1 % (1/11), both were lower than former prevalence data (10.6% and 29.7%)[19]. The differences in E. bieneusi prevalence between studies may be influenced by geographical region, sample size, sampling time, animal health condition, breeding mode, and population density[34–36]. The role of those factors can be further examined. As the closest phylogenetic relationship animal with red panda [37], the American raccoon (Procyon lotor), also possessed a higher prevalence (27.3%, 15/55) [27, 31]. Red pandas have a low overall prevalence compared to other common zoo species with a sampling size of more than 20 isolates (Table 4). Ignoring the interference of other factors, the adaptation of E.bieneusi to the red panda may be in the middle level.
Table 4
The prevalence of E.bieneusi in some species of Chinese zoos.
Family
|
Species
|
No. of tested
|
No.of positive (%)
|
Reference
|
Ailuridae
|
Red panda (Ailurus fulgens)
|
198
|
24(12.1%)
|
This study
|
Cercopithecidae
|
Golden snub-nosed monkey (Cercopithecus kandti)
|
160
|
74(46.2%)
|
[26]
|
|
Rhesus macaque (Macaca mulatta)
|
304
|
60(19.7%)
|
[16, 33]
|
|
Hamadryas baboon (Papio hamadryas)
|
21
|
6 (28.6%)
|
[42]
|
|
Cynomolgus monkey (Macaca fascicularis)
|
80
|
47(58.8%)
|
[16, 18]
|
Cebidae
|
Squirrel monkey (Saimiri sp.)
|
43
|
17(39.5%)
|
[16]
|
|
Black-capped capuchin (Cebus apella)
|
22
|
6 (27.3%)
|
[42]
|
Ursidae
|
Giant panda (Ailuropoda melanoleuca)
|
200
|
69(34.5%)
|
[16]
|
|
Asiatic black bears (Ursus thibetanus)
|
106
|
29(27.4%)
|
[42]
|
Bovidae
|
Golden takin (Budorcas taxicolor bedfordi)
|
191
|
28(14.7%)
|
[42]
|
Anatidae
|
Whooper swans (Cygnus cygnus)
|
467
|
35(7.5%)
|
[42]
|
Moschidae
|
Musk deer (Moschus berezovskii)
|
223
|
38(17.0%)
|
[42]
|
Lemuridae
|
Ring-tailed lemur (Lemur catta)
|
45
|
11(24.0%)
|
[42]
|
Hominidae
|
Bornean orangutan (Pongo pygmaeus)
|
23
|
4 (17.4%)
|
[16]
|
Macropodidae
|
Red kangaroo (Macropus rufus)
|
38
|
14(36.8%)
|
[38]
|
|
Grey kangaroo (Macropus fuliginosu)
|
23
|
0 (0.0%)
|
[38]
|
Four observed genotypes in this study are completely different from genotypes previously identified in red panda (EBPC, CHB1)[13, 19, 28]. As the dominant genotype, PL2 was only found in masked palm civets (Paguma larvata) and formed a clade related to Group 2 in phylogenetic analysis. [33, 39]. All positive isolates from CDRBGPB were PL2, which revealed red panda is a new host of genotype PL2, which expands the host range of this genotype. The high prevalence in CDRBGPB indicates that red pandas may occur inter-species transmission and be an adaptive host of PL2, but the infection route is unknown and needs more future research. Long-term monitoring and further study are also required to ascertain whether the semi-free-range breeding mode has increased the transmission of E. bieneusi within the red panda populations and the potential risk of cross-specific transmission. Both Chengdu Zoo and Bifengxia Zoo observed only 1 genotype D isolates. Considering that genotype D has previously been discovered in the above-mentioned zoos, pathogens originated from other animal hosts are possible [26, 40]. Genotype SC02 was once detected in humans and various animals such as Asian black bears, Tibetan blue bears, sun bears, raccoons, horses, and squirrels [19, 41, 42]. The exchange and recombination in ITS genes of E.bieneusi may have resulted in genetic diversity and the new genotype SCR1 [43].
MLST analysis showed that amplification efficiency is different in various ITS genotypes, which is consistent with previous studies [11]. We found 3 MLGs, which showed higher genetic diversity in PB isolates than the other 5 sampling sites; however, the influence of sample size cannot be excluded. MLST tools revealed higher resolution and genetic diversity than ITS sequence genotyping, which is according to the conclusions of previous studies[18, 41, 42, 44]. The amplification efficiency of locus MS4 was higher than MS1, MS4, and MS7, which was inconsistent with previous reports[14, 45].