華 彥 郎鵬飛 王洪堂 王 姣 曹和琴 姜廣順*

(1.東北林業(yè)大學(xué)野生動(dòng)物與自然保護(hù)地學(xué)院，哈爾濱，150040；2.廣東省林業(yè)科學(xué)研究院，廣州，510520)

華北豹(Pantherapardusjaponensis)是中國特有的豹亞種，曾廣布于中國的華北地區(qū)，由于棲息地喪失和過度捕獵等因素，近年來，華北豹的分布急劇退縮，現(xiàn)僅分布于中國的山西、河南、河北等省份[1]，在動(dòng)物園飼養(yǎng)有少量野外救護(hù)個(gè)體。胃腸道慢性炎癥[2-4]等病癥與腸道微生態(tài)紊亂有關(guān)，尤其是對區(qū)域性的極小種群而言，動(dòng)物的腸道疾病可能成為物種受脅的潛在重要因素。腸道微生物主要由細(xì)菌和少量真菌及其他微生物組成[5]，真菌雖然占腸道微生物比例低，但其群落結(jié)構(gòu)與物種、食性、疾病等因素密切相關(guān)[6]。近年來，腸道微生物中真菌的研究越來越受關(guān)注，Albenberg等[7]曾報(bào)道了炎癥性腸病(inflammatory bowel disease，IBD)腸道微生物群中真菌特征，發(fā)現(xiàn)患病樣本中特異性念珠菌分類群的豐富度有明顯的增加，腸道真菌紊亂是導(dǎo)致IBD發(fā)生、發(fā)展的重要因素[8]?？股厥谦F醫(yī)常用的藥物，研究表明長期使用抗生素會導(dǎo)致腸道菌群失衡，具有保護(hù)作用的細(xì)菌減少，真菌數(shù)量增加，從而導(dǎo)致腹瀉[9]。Li等[10]發(fā)現(xiàn)腹瀉牦牛(Bosgrunniens)腸道新麗鞭毛菌門(Neocallimastigomycota)顯著高于健康牦牛，新麗鞭毛菌門豐度過高可能是導(dǎo)致腹瀉的主要原因。

因腸道菌群失調(diào)造成其微生態(tài)失衡而引起一系列腸道疾病是人工圈養(yǎng)瀕危野生動(dòng)物主要原因之一。彭智偉[11]對野生東北虎(Pantheratigrisaltaica)、東北豹(Pantherapardusorientalis)等大型貓科(Felidae)動(dòng)物腸道微生物中的細(xì)菌群落進(jìn)行了分析，但對腸道真菌的研究顯見報(bào)到。本研究通過高通量測序技術(shù)對圈養(yǎng)和野生華北豹糞便中真菌類型進(jìn)行研究，以探究華北豹腸道真菌多樣性和優(yōu)勢菌群，比較人工圈養(yǎng)與野生兩種條件下華北豹腸道真菌菌群的差異，為華北豹的遷地保護(hù)提供科學(xué)管理的依據(jù)。

1 材料與方法

1.1 材料

在山西鐵橋山省級自然保護(hù)區(qū)采集24份野生華北豹新鮮糞便樣品，在山西太原動(dòng)物園和河南森林公園采集野外救助的人工圈養(yǎng)華北豹糞便樣品7份。每份樣品記錄了采集地坐標(biāo)及日期等信息，裝于樣品袋中，-80 ℃保存。

1.2 華北豹腸道真菌DNA提取

真菌胞壁比原核生物要厚10倍左右，所以先將樣本在超低溫冰箱和37 ℃水浴鍋中反復(fù)凍融5次，以提高真菌DNA的提取效果。通過E.Z.N.ATM Mag-Bind Soil DNA Kit的試劑盒進(jìn)行糞便樣品總DNA的提取。

1.3 華北豹腸道真菌18S rRNA的PCR擴(kuò)增

利用引物547F(5′-CCAGCASCYGCGGTAATTCC-3′)，引物V4R(5′-ACTTT-CGTTCTTGATYRA-3′)靶向18S SSU rRNA基因的V4區(qū)[12]。PCR混合物包含16 μL高壓滅菌蒸餾水，10 μL PCR Buffer(5x)，4 μL dNTPs(2.5 mM)，10 μL GC Buffer(5x)，5.5 μL DNA模板，0.5 μLTaqE，正向和反向引物(工作濃度：10 μM)各2 μL，每個(gè)反應(yīng)總體積為50 μL。35個(gè)PCR循環(huán)包括在98 ℃下持續(xù)25 s，55 ℃持續(xù)30 s和72 ℃持續(xù)50 s，在此之前在98 ℃下持續(xù)5 min的預(yù)變性，并在35個(gè)PCR循環(huán)后以72 ℃持續(xù)7 min后結(jié)束整個(gè)反應(yīng)。在Gold-View染色的2%瓊脂糖凝膠上電泳后分析PCR產(chǎn)物。

1.4 18S rRNA的文庫制備和Miseq測序

采用Illumina公司的TruSeq Nano DNA LT Library Prep Kit制備測序文庫，之后進(jìn)行文庫質(zhì)檢與定量，對合格的文庫，在MiSeq機(jī)器上利用Mi Seq Reagent Kit V3(600個(gè)循環(huán))進(jìn)行 2×300 bp 的雙端測序。通過對雙端序列數(shù)據(jù)的數(shù)據(jù)過濾和雙末端拼接操作，得到有效數(shù)據(jù)，再利用 FLASH 軟件將有效 reads 拼接組裝成一條序列。

為得到華北豹腸道中較好的真菌物種多樣性信息，對獲得的原始序列使用 QIIME 軟件對測序數(shù)據(jù)進(jìn)行優(yōu)化處理。去除接頭和無重復(fù)單序列，以便降低冗余計(jì)算量。以97%相似度，將所有樣品優(yōu)化后的序列聚類成為 OTU(operational taxonomic units)。

1.5 華北豹腸道真菌多樣性的分析

基于OTU聚類結(jié)果，利用Mothur軟件計(jì)算華北豹腸道真菌的Alpha多樣性。Alpha多樣性指數(shù)可用來對各樣品的豐富度和多樣性進(jìn)行比較分析，常用的Alpha多樣性指數(shù)包括Chao1指數(shù)、Shannon指數(shù)以及Simpson指數(shù)。Chao1指數(shù)表明樣品中群落的豐富度，Shannon 指數(shù)呈現(xiàn)群落的多樣性，而Simpson指數(shù)則是群落中優(yōu)勢種集中程度的體現(xiàn)。

利用 Mothur 軟件做豐富度稀釋性曲線及Shannon指數(shù)分析，作為測序數(shù)據(jù)量不同樣本中物種豐富度比較的依據(jù)，也可作為樣本測序量合理與否(是否足以覆蓋所有類群)的評價(jià)標(biāo)準(zhǔn)，得出樣本測序深度，同時(shí)反映華北豹腸道真菌多樣性指數(shù)，顯示華北豹腸道中絕大多數(shù)的真菌信息。

2 結(jié)果與分析

2.1 OTU數(shù)目與Alpha多樣性指數(shù)分析

31份華北豹糞便樣品原始序列共計(jì)1 352 804條，平均每份樣本43 638.84條；質(zhì)控后的序列共計(jì)1 209 271條，平均每份樣本39 008.74條。樣品質(zhì)控后的序列讀長主要位于300—390 bp，所有樣品優(yōu)化聚類后，共得到1 797個(gè) OTU(圖1)。

圖1 測序序列質(zhì)量控制前后讀長分布Fig.1 The reading lengths of the sequences

本實(shí)驗(yàn)中的Chao1指數(shù)為3—218，表明各樣品的豐富度不高。Shannon指數(shù)與Simpson指數(shù)的數(shù)值越大，則表示樣品中物種的多樣性越高，物種分布的均勻度越好。覆蓋率的大小則反映了樣品測序的深度，本實(shí)驗(yàn)31份樣品的覆蓋率均在99.23%以上，說明該實(shí)驗(yàn)測序結(jié)果能很好地反應(yīng)樣品中真菌的情況。

從圖2與圖3分別是野生與圈養(yǎng)華北豹樣品中真菌的豐富度稀疏曲線和Shannon指數(shù)稀疏曲線，從圖中可以看出當(dāng)測序量達(dá)到3 000時(shí)，各樣品的OTU數(shù)目增加基本趨向平緩，這說明更多的數(shù)據(jù)量只會產(chǎn)生少量OTU數(shù)目，這也表示該實(shí)驗(yàn)的測序數(shù)據(jù)量合理。從圖中也可明顯看出，野外與圈養(yǎng)華北豹所有樣品Shannon曲線，均已經(jīng)達(dá)到平緩的走勢，這也反映了樣品量較為合理。

表1 野生和圈養(yǎng)華北豹腸道真菌多樣性指數(shù)

Tab.1 The diversity of Alpha in wild leopard and captive leopard

2.2 物種組成及優(yōu)勢菌群分析

2.2.1 門水平上腸道真菌菌群組成及豐富度

對24份野生華北豹糞便樣品、7份圈養(yǎng)華北豹糞便樣品在門水平和屬水平上進(jìn)行物種注釋分析，可確定其腸道真菌菌群中的 6個(gè)菌門，無法鑒別菌群的序列這里統(tǒng)一歸并為其他(unidentified)。其中子囊菌門(Ascomycota)在華北豹腸道真菌菌群中含量都是最高的，各樣品的平均含量為61%、擔(dān)子菌門(Basidiomycota)的平均含量為19%、毛霉門(Mucoromycota)的平均含量為9%、壺菌門(Chytridiomycota)的平均含量為3%、芽枝霉門(Blastocladiomycota)的平均含量為0.05%，以及未確定(unidentified)的平均含量為6%，Zoopagomycota菌門只在野生華北豹個(gè)體樣品中檢測到平均含量為2%。

2.2.2 野生和圈養(yǎng)華北豹腸道真菌群落比較

圖5是野生和圈養(yǎng)華北豹這2組樣品共有和特有 OTU 數(shù)目的 Venn 圖。2組樣品共有631個(gè)OTU，重疊部分表示2組共有的 OTU 數(shù)目為 97個(gè)，野外華北豹組(YW)特有的OTU為 519個(gè)，圈養(yǎng)華北豹組(JY)特有的 OTU 個(gè)數(shù)為14個(gè)。野外特有的真菌有柄孢殼菌屬(Podospora)、念珠菌屬(Candida)、Cystofilobasidium、隱球菌屬(Cryptococcus)、Rhodosporidiobolus、隱囊菌屬(Aphanoascus)、Saitozyma、橫節(jié)霉菌屬(Chalara)、Helicodendron、子囊殼菌屬(Coniochaeta)；圈養(yǎng)特有的有復(fù)膜孢酵母菌屬(Saccharomycopsis)、Cyniclomyces、Lodderomyces、Gamsiella、Diutina、克魯維酵母屬(Kluyveromyces)、孢圓酵母(Torulaspora)。

圖2 豐富度稀疏曲線Fig.2 The richness rarefaction curve

圖3 香農(nóng)指數(shù)稀疏曲線Fig.3 The Shannon index rarefaction curve

圖4 野生與圈養(yǎng)華北豹腸道真菌門水平下的組成結(jié)構(gòu)Fig.4 The bar chart of fungal on phylum in wild leopard and captive leopard

圖5 野生與圈養(yǎng)華北豹樣本中獨(dú)有及共享操作分類單元(OTU)Venn圖及顯示各屬水平的比例餅形圖Fig.5 Venn diagram showing unique and shared operational taxonomic units(OTUs)between wild leopard and captive leopard samples.A pie chart of the respective proportions of genus is shown

2.2.3 屬水平上的物種豐度聚類熱圖

圖6所示野生華北豹與圈養(yǎng)華北豹腸道真菌屬水平上的豐度，選用豐度排名前50的屬繪制的聚類分析熱圖，圖中色塊越靠近紅色表示某份樣品對應(yīng)的某一物種比例越高，反之，色塊越靠近綠色表示某份樣品對應(yīng)的某一物種比例越低。總體來看，野生華北豹組較圈養(yǎng)華北豹組的物種豐度有變化，說明野生華北豹較圈養(yǎng)華北豹組間的物種多樣性存在差異。圖中可見，樣本19獨(dú)自為一組，其余各組腸道內(nèi)真菌菌群的組成及相對豐度差異較大，特有相對豐度顯著高于各組平均值的真菌有Polarella、Phaeocystis、Pelagophyceae等。

圖6 野生華北豹與圈養(yǎng)華北豹屬水平上的物種豐度熱圖Fig.6 The abundance of fungal heatmap in wild leopard and captive leopard on genus

圖7 LDA 值分布柱狀圖Fig.7 Distribution histogram based on LDA value

圖8 進(jìn)化分支圖LichtheimiaceaeFig.8 Evolutionary branch graph

2.2.4 基于LefSe(LDA Effect Size)分析野外與圈養(yǎng)華北豹2組間的物種差異

采用LefSe分析野外與圈養(yǎng)華北豹2組間具有顯著差異的物種(圖7，圖8)。LDA評分直方圖(圖7)顯示各組間具有統(tǒng)計(jì)學(xué)差異的生物學(xué)標(biāo)志，柱形圖的長度表示LDA評分，紅色代表野外組樣品，綠色代表圈養(yǎng)樣品。該圖顯示了野外與圈養(yǎng)組之間具有顯著差異的微生物類群(LDA score>3.3)。在野外組中具有顯著豐度差異的有：Coniochaetaceae科、銀耳綱(Tremellomycetes)、小子囊菌目(Microascales)、子囊殼菌屬；在圈養(yǎng)組中具有顯著豐富度差異的有：毛霉門、酵母綱(Saccharomycetes)、畢赤酵母科(Pichiaceae)、酵母菌屬(Saccharomyces)、毛霉目(Mucorales)、酵母科(Saccharomycetaceae)、環(huán)孢菌科(Inocybaceae)、根毛霉科(Rhizomucor)、德巴利(氏)酵母屬(Debaryomyces)、酵母菌目(Saccharomycetales)、畢赤酵母屬(Pichia)、絲蓋傘屬(Inocybe)、復(fù)膜孢酵母屬、金絲酵母科(Debaryomycetaceae)、Lichtheimiaceae、復(fù)膜孢酵母科(Saccharomycopsidaceae)。進(jìn)化分支圖(圖8)由內(nèi)至外輻射的圓圈，顯示了野外與圈養(yǎng)組在門至屬(或種)水平的分類群相對豐度差異。每個(gè)不同分類的小圓圈表示一個(gè)分類單元，而每個(gè)圓圈的直徑的大小反映了相對分度的高低。紅色和綠色的圓圈分別代表了野生和圈養(yǎng)華北豹相對豐度差異，未標(biāo)色的為無顯著差異的物種。從圖中結(jié)果顯示，LDA值分布柱狀圖與進(jìn)化分支圖的結(jié)果是相符合的。

3 討論

3.1 華北豹腸道真菌與其他動(dòng)物的比較

本研究發(fā)現(xiàn)子囊菌門、擔(dān)子菌門、毛霉門、壺菌門、芽枝霉門是華北豹腸道真菌的優(yōu)勢菌群，艾生權(quán)[5]等用RLFE法測得亞成體大熊貓腸道真菌主要為子囊菌門、擔(dān)子菌門2個(gè)門，周本翔等[12]通過實(shí)驗(yàn)證明草魚腸道中真菌同樣以子囊菌門數(shù)量最多，其次是擔(dān)子菌門。這些研究結(jié)果表明，至少在以上的幾種動(dòng)物類群間，其腸道真菌構(gòu)成大致是相似的。

3.2 野生華北豹與圈養(yǎng)華北豹腸道真菌多樣性的差異

真菌菌群占腸道總菌群的比例較小，但是對腸道的穩(wěn)態(tài)平衡維持具有重要的作用[13]，腸道中的真菌與細(xì)菌關(guān)系密切，兩者存在著穩(wěn)定的拮抗、互利共生關(guān)系[14-15]，它們?yōu)楸舜说纳L發(fā)育提供支持和保護(hù)作用[16]。本研究發(fā)現(xiàn)，圈養(yǎng)華北豹酵母綱較野生華北豹豐富。酵母綱作用動(dòng)物腸道的機(jī)理十分復(fù)雜，現(xiàn)僅見白色念珠菌屬(Candidaalbicans)對腸道黏膜損傷之后致炎作用，但酵母菌屬中的布拉酵母菌(Saccharomycesboulardii)是目前確認(rèn)人類腸道益真菌，還未見到其對于貓科食肉動(dòng)物的相關(guān)研究報(bào)道[17]。假單胞菌的生長、發(fā)育等活動(dòng)與白酵母菌密切相關(guān)[18]。Kohl等[19]曾研究報(bào)道指出飲食并不是圈養(yǎng)動(dòng)物微生物多樣性減少的唯一原因，并指出圈養(yǎng)動(dòng)物后代的腸道菌群與它們母親的腸道菌群結(jié)構(gòu)更相近。酵母菌屬中的白色念珠菌是哺乳動(dòng)物胃腸道中含量最多的一種真菌，有小鼠實(shí)驗(yàn)提示白色念珠菌的致炎作用更可能發(fā)生在腸道黏膜損傷之后[20]。Sokol等[21]通過觀察235名IBD患者與38名健康者的糞便微生物群的細(xì)菌與真菌組成，發(fā)現(xiàn)IBD患者擔(dān)子菌門與子囊菌門的比例增加，其中促炎真菌如白色念珠菌占比增加，抗炎真菌如酵母菌占比減少。

通過本研究積累了相關(guān)數(shù)據(jù)及資料，發(fā)現(xiàn)野生和圈養(yǎng)華北豹腸道真菌組成結(jié)構(gòu)中毛霉門存在明顯差異，這可能是由于野生與圈養(yǎng)條件下食物豐富度不同而導(dǎo)致的。因此飼養(yǎng)人員應(yīng)當(dāng)對圈養(yǎng)華北豹的食物構(gòu)成加以改善，以免引起華北豹因食物單一而引發(fā)疾病。

本實(shí)驗(yàn)首次就圈養(yǎng)與野外華北豹腸道真菌菌群進(jìn)行研究，目前國內(nèi)外對本實(shí)驗(yàn)中檢測到的一些真菌研究甚少，它們與機(jī)體及腸道微生態(tài)的相關(guān)性也不甚明確，無法就相關(guān)菌群在宿主腸道微生態(tài)中扮演的具體作用進(jìn)行詳細(xì)的分析，還要相關(guān)的后續(xù)實(shí)驗(yàn)進(jìn)一步研究。