動物園圈養大熊貓畢氏腸微孢子蟲分子流行病學調查

曹鈺瑩，李 威，鐘志軍，肖 潔，弓 超，王承東，李才武，楊海迪，黃祥明，吳孔菊，*，彭廣能，*

(1.四川農業大學 動物醫學院 動物疫病與人類健康四川省重點實驗室，四川 成都 611130； 2.中國大熊貓保護研究中心臥龍基地，四川 阿壩 624000； 3.成都大熊貓繁育研究基地，四川 成都 610081)

畢氏腸微孢子蟲(Enterocytozoonbieneusi)是一種專性細胞內寄生的單核真核生物，是一種重要的人畜共患病原體。它可以導致免疫缺陷者發生脫水性致命性的腹瀉，也會導致免疫力正常動物的自限性腹瀉。畢氏腸微孢子蟲屬真菌界(Fungi)，微孢子蟲門(Microsporidia)，無膜泛成孢子蟲亞目(Apansporoblastina)，腸孢子蟲科(Enterocy-tozoonidae)，腸孢子蟲屬(Enterocytozoon)[1-2]。畢氏腸微孢子蟲(E.bieneusi)在自然界分布廣泛，從無脊椎動物到脊椎動物乃至人中均有發現[3]，其孢子也廣泛存在于水源中。美國環境保護協會(EPA)將其列為能夠導致水源性暴發流行的微生物類污染物，美國國立衛生研究院(NIH)將其列為B類生物防范病原體[4]。

畢氏腸微孢子蟲核糖體RNA(rRNA)基因的內部轉錄間隔區(ITS)具有高度的多樣性，可被用于E.bieneusi的基因分型[5]。目前，根據畢氏腸微孢子蟲(E.bieneusi)的ITS基因型，將其劃分為8個組(group 1～group 8)，其中group 1中的基因型能夠感染人和動物[6]。

大熊貓是我國重要的保護動物。據第4次大熊貓普查的統計結果，目前，我國存活的大熊貓僅有2 239只，其中，野生大熊貓1 864只，圈養大熊貓375只[7]。有研究發現，繁殖力低下、棲息地減少及食物缺乏等因素是造成大熊貓數量稀少的主要原因，除此之外，寄生蟲的高感染率也嚴重威脅著大熊貓的生存，現已成為大熊貓保護工作不可忽視的重要內容[8]。本研究采用巢氏PCR基于ITS基因擴增樣品DNA，對國內10個省(市、自治區)動物園內的大熊貓進行畢氏腸微孢子蟲流行病學調查，為圈養大熊貓畢氏腸微孢子蟲病的防控奠定一定的理論基礎。

1 材料與方法

1.1 樣品采集與保存

本研究于2017年采集國內江蘇、福建、浙江等10個省(市、自治區)的動物園內圈養大熊貓新鮮糞便共計31份。將采集的新鮮糞便分別置于密封袋中，相應標記并編號，并放置于4 ℃冰箱保存備用。

1.2 主要試劑

糞便DNA提取試劑盒購自Omega公司。2×TaqPCR Master Mix 購自天根生化科技(北京)有限公司。DL2000 DNA Marker、BSA均購自寶生物工程(大連)有限公司。

1.3 糞便DNA提取

取200 mg糞便于2 mL離心管中，按照DNA提取試劑盒的步驟進行大熊貓糞便的提取。將提取的DNA分裝，一部分保存于4 ℃待檢，另一部分保存于-20 ℃備用。

1.4 巢氏PCR

將提取的DNA進行巢氏PCR擴增，擴增其ITS基因，引物及退火溫度參考Sulaiman等[9]的文獻，詳細見表1。擴增體系均為25 μL，其中2×TaqPCR MasterMix 12.5 μL，去離子水7.5 μL，0.1% BSA 1.0 μL，上、下游引物各1.0 μL，模板DNA 2.0 μL。

1.5 產物檢測及測序

取5 μL擴增產物，將其點入1%瓊脂糖凝膠孔中，在130 V電壓下電泳1 h，然后于電泳凝膠成像系統儀中拍照。將出現的與目的條帶一致的擴增產物送至英濰捷基(上海)貿易有限公司進行雙向測序。測序結果在GenBank中經BLAST比對校正后提交至NCBI，獲得登錄號。

1.6 序列分析及種系發育分析

從NCBI分別下載ITS位點相關序列，用Clustal X(1.83)按照默認設置進行比對，然后用MEGA 6.06軟件NJ法構建系統發育進化樹。

F1， 第一套上游引物； R1， 第一套下游引物； F2， 第二套上游引物； R2， 第二套下游引物。

F1， The first set of upstream primers; R1， The first set of downstream primers; F2， The second set of upstream primers; R2， The second set of downstream primers.

2 結果與分析

2.1 大熊貓畢氏腸微孢子蟲感染情況

PCR結果顯示，實驗的31份大熊貓新鮮糞便中，畢氏腸微孢子蟲陽性共14份，總感染率45.2%，其中雌性大熊貓感染率為50%(6/12)，雄性大熊貓感染率為42.1%(8/19)(表2)，不同性別的大熊貓畢氏腸微孢子蟲感染率無顯著差異(P>0.05)。亞成體(年齡在1.5～5.5歲)大熊貓畢氏腸微孢子蟲感染率為100%(1/1)；成體(年齡>5.5歲)的大熊貓畢氏腸微孢子蟲感染率為43.3%(13/30)(表2)。在此次的實驗中還發現，有2個省份的動物園圈養大熊貓畢氏腸微孢子蟲感染率為100%(表3)，另有3個省份的感染率為0，有2個省份的感染率為50%，剩余的感染率分別為67.7%和33.3%(表3)。

2.2 ITS基因型及種系發育分析

對畢氏腸微孢子蟲陽性分離株進行PCR擴增、測序，發現了2種感染大熊貓的畢氏腸微孢子蟲基因型，分別為SC02(n=13)(GenBank登錄號：KY950533)和Peru6(n=1)(GenBank登錄號：KY950540)，其中，基因型Peru6來自于成都動物園內的感染大熊貓?；贗TS序列構建的種系發育進化樹(圖1)表明，SC02與Peru6均屬于group 1b，具有人畜共患潛能。

3 討論

近年來，由于對畢氏腸微孢子蟲的研究逐漸增多，在器官移植受者、兒童、老年人和惡性疾病與糖尿病患者等人群中發現了若干例感染的情況[10]。隨著對畢氏腸微孢子蟲研究的逐漸深入，畢氏腸微孢子蟲也越來越多地在家畜以及野生動物上被發現。

本次調查的圈養大熊貓畢氏腸微孢子蟲感染率為45.2%，是此前Tian等[11]在陜西秦嶺野生動物園等的大熊貓上發現的畢氏腸微孢子蟲感基因型括號內的數字為樣品數。

表2 不同年齡及性別的圈養大熊貓畢氏腸微孢子蟲感染率及基因型

The values in brackets are sample numbers of the corresponding genotype

圖1 大熊貓源畢氏腸微孢子蟲ITS基因的系統發育分析Fig.1 Phylogenetic relationships of ITS nucleotide sequences of the Enterocytozoon bieneusi in giant panda

染率(8.7%，4/46)的5倍；而Deng等[12]在中國西南地區5個動物園內亞洲黑熊上發現的畢氏腸微孢子蟲感染率為27.4%(29/106)，低于此次實驗所得感染率。造成感染率差距如此大的原因可能與動物的種類、身體狀況、年齡大小，以及其飼養密度等有關。在對比各地動物園大熊貓的畢氏腸微孢子蟲感染率后發現，此次實驗中沿海地區的動物園感染率較高，如廣東、福建、上海等(表3)，這可能與地理位置及氣候有關。從表3還可以看出，在感染率100%的動物園內發現的畢氏腸微孢子蟲陽性分離株的基因型均為SC02(如廣東、湖南)，而大多數只存在一只大熊貓的動物園內畢氏腸微孢子蟲感染率為0，因此，也可能存在動物園內大熊貓交叉感染的可能性。另外，高感染率還可能與大熊貓“搬家”造成的應激有關：這些都需要再進一步研究。此外，在本研究中不同性別及年齡的大熊貓的畢氏腸微孢子蟲感染率無顯著差異。

本實驗檢測到的圈養大熊貓畢氏腸微孢子蟲優勢基因型SC02(92.6%，13/14)，SC02基因型是2016年Li等[13]從圈養西藏棕熊、北方浣熊等熊科動物上發現并命名的畢氏腸微孢子蟲新基因型，根據Li等[13]的研究發現，基因型SC02與感染人的基因型R同屬于畢氏腸微孢子蟲的ITS基因分組中group 1的1b，具有潛在的人畜共患可能性，目前已經在中國人群發現(GenBank登錄號KY465443)。另外，該基因型也在亞洲黑熊、月熊等熊科動物和馬中被發現[12,14]。SC02基因型首次在大熊貓中鑒定到，拓寬了該基因型的宿主范圍。在本次實驗中鑒定到的畢氏腸微孢子蟲基因型Peru6之前也在成都動物園中的一只大熊貓中鑒定到，該基因型除了在大熊貓上被發現，之前也在鳥類、犬、牛等動物，以及人上有過報道[1,5,13,15]。有研究證明，基因型Peru6屬于group 1中的1b[16]，具有人畜共患特性。

本實驗利用ITS基因位點對圈養大熊貓畢氏腸微孢子蟲的感染情況進行調查，并開展了遺傳進化分析，發現感染大熊貓的2種畢氏腸微孢子蟲基因型(SC02和Peru6)皆屬于group 1，具有潛在的人畜共患可能性。SC02基因型是第一次在大熊貓上被發現，擴大了SC02的宿主范圍。由于圈養大熊貓與人類接觸頻繁，增加了大熊貓源畢氏腸微孢子蟲感染人的可能性。