一只半野生病死川金絲猴腸道微生物多樣性分析

劉燕鄭常明李林海張成林*劉學鋒姚輝余輝亮李瑩楊明海賈婷盧巖夏茂華

(1.北京動物園，北京，100044；2.圈養(yǎng)野生動物技術北京市重點實驗室，北京，100044；3.國家林業(yè)局林產(chǎn)工業(yè)規(guī)劃設計院，北京，100010；4.湖北省金絲猴保護研究中心，神農(nóng)架，442421)

一只半野生病死川金絲猴腸道微生物多樣性分析

劉 燕1，2鄭常明1，2李林海張成林1，2*劉學鋒1，2姚 輝4余輝亮李 瑩1，2楊明海1，2賈 婷1，2盧 巖1，2夏茂華1，2

(1.北京動物園，北京，100044；2.圈養(yǎng)野生動物技術北京市重點實驗室，北京，100044；3.國家林業(yè)局林產(chǎn)工業(yè)規(guī)劃設計院，北京，100010；4.湖北省金絲猴保護研究中心，神農(nóng)架，442421)

稿件運行過程

川金絲猴； 腸道菌群； 多樣性分析； 投食群； 神農(nóng)架

川金絲猴(Rhinopithecusroxellana)是中國特有的珍稀動物，具有極高的生態(tài)地位和研究價值，僅次于大熊貓，被譽為中國的第二國寶。神農(nóng)架保護區(qū)內(nèi)的川金絲猴種群相對孤立，主要包括大龍?zhí)锻妒橙汉托↓執(zhí)度︷B(yǎng)種群。大龍?zhí)兜暮锶鹤?006年開始進行人工投食以來，種群數(shù)量正在逐年增加，種群內(nèi)和種群間存在著個體遷移，近交水平較低，遺傳多樣性得以維持[1]。

本次死亡川金絲猴是1只老年雄性金絲猴，年齡在20歲左右。自2006～2015年間一直是大龍?zhí)侗Ｗo站投食種群中一個16只金絲猴家庭的家長。其在家庭中地位較為穩(wěn)固，對家庭中雌性成員較好，與它們分享食物，幫它們梳理毛發(fā)，經(jīng)常協(xié)助照顧幼猴，并干預家庭中雌性之間的沖突等。2015年6月9日該川金絲猴出現(xiàn)咳嗽，16日服用乙酰螺旋素片2片，治療無效，25日因肺炎死亡。隨后對其進行剖檢，采集4個部位的腸內(nèi)容物利用高通量測序技術進行腸道菌群分析[2-4]。

1 材料與方法

1.1樣品采集

川金絲猴死后，立即剖檢，采集十二指腸段、小腸(包括空腸和回腸)段、盲腸段及大腸(包括結腸和直腸)段共4份樣品進行-80℃保存。

1.2腸道菌群基因組DNA的提取

使用E.Z.N.A.@Stool DNA Kit(Omega Bio-tek，Norcross，GA，U.S.)，均按試劑盒要求操作。提取DNA，用紫外分光光度計測定其A260值，計算其濃度，通過A260/A280值鑒定純度，進行0.8%的瓊脂電泳驗證其完整性。

1.3 16S rRNA的PCR擴增和Illumina MiSeq測序

以樣品中菌群總DNA 為模板，用引物338F(5′-ACTCCTACGGGAGGCAGCA-3′)和806R(5′-GGACTACCAGGGTATCTAAT-3′)，對細菌16S rRNA 的V3 ～ V4 區(qū)域進行PCR擴增。PCR采用20 μL 反應體系，包括0.4 μL FastPfu DNA Polymerase(購買于北京全式金生物技術有限公司)，濃度為5U/μL，4 μL 5×FastPfu Buffer，引物濃度均為5 μmol / L，各0.4 μL，2.5 mmol/LdNTP 2 μL，總DNA 10 ng，補ddH2O 至20 μL。

PCR 過程在ABI Gene Amp 9700 上進行，反應條件為：95 ℃ 預變性3 min；隨后變性循環(huán)27次(95 ℃，30 s；55 ℃，30 s；72 ℃，45 s 為1 個循環(huán))；72 ℃ 延伸10 min，10 ℃ 冷卻。PCR 產(chǎn)物用含溴化乙錠的2%瓊脂糖凝膠電泳鑒定并用PCR純化試劑盒Qiagen純化和回收，回收的PCR產(chǎn)物用QuantiFlourTM-ST fluorometer定量，并按測序要求進行混合。然后送至上海美吉生物醫(yī)藥科技有限公司。

1.4 ITS rRNA的PCR擴增和Illumina MiSeq測序

以樣品中菌群總DNA為模板，用真菌通用引物1737F(5′-GGAAGTAAAAGTCGTAACAAGG-3′)和2043R(5′-GCTGCGTTCATCGATGC-3′)，針對真菌ITS區(qū)域進行PCR擴增。PCR采用反應體系同上。PCR 過程在ABI Gene Amp 9700 上進行，反應條件為：95 ℃ 預變性3 min；隨后變性循環(huán)32 次(95 ℃，30 s；55 ℃，30 s；72 ℃，45 s 為1 個循環(huán))；72 ℃ 延伸10 min，10 ℃ 冷卻。其他同上。

1.5數(shù)據(jù)分析

首先對原始數(shù)據(jù)進行過濾處理，得到優(yōu)化序列。然后在去除嵌合體序列后進行OTU聚類分析，對OTU的代表序列作分類學分析。基于OTU聚類分析結果，對OTU進行多種多樣性指數(shù)分析，以及對測序深度的檢測；基于分類學信息，在各個分類水平上進行群落結構的統(tǒng)計分析。

2 結果

2.1樣品信息

將采集的十二指腸段、小腸段、盲腸段及大腸段4份腸內(nèi)容物，分別命名為SJ1、SJ2、SJ3和SJ4。提取微生物DNA，進行測序和分析。

2.2樣品測序深度

用Illumina MiSeq平臺對細菌16S rRNA的V3～V4區(qū)進行測序，共獲得184 412條有效序列(cleandata)，平均每個樣品含有(46 103)條序列，每條序列的平均長度是439 bp?；谙嗨贫却笥?7%的原則，將獲得的有效序列進行OTU聚類，共獲得240個OTU，其中SJ1含181個、SJ2含208個、SJ3含220個和SJ4含213個。根據(jù)樣品的稀釋曲線(圖1)，在本試驗的測序深度下，全部樣品的稀疏曲線均趨于平緩，獲得的OTU數(shù)量變化不大，說明本試驗的測序量可以覆蓋各樣品的大多數(shù)微生物，測序數(shù)據(jù)量合理。

用Illumina MiSeq平臺對真菌的ITS片段進行測序，共獲得103 316條有效序列，平均每個樣品含有(25 829)條序列，每條序列的平均長度是241 bp。基于相似度大于97%的原則，將獲得的有效序列進行OTU聚類，共獲得341個OTU，其中SJ1含150個、SJ2含166個、SJ3含129個和SJ4含104個。根據(jù)樣品的稀釋曲線(圖2)，在本試驗的測序深度下，全部樣品的稀疏曲線均趨于平緩，獲得的OTU數(shù)量變化不大，說明本試驗的測序量可以覆蓋各樣品的大多數(shù)微生物，測序數(shù)據(jù)量合理。

2.3 α多樣性分析

應用軟件mothur version v.1.30.1進行ace、chao、coverage、shannon、simpson 5種常用指數(shù)的測定(表1)，4份樣品的序列覆蓋率coverage均達到99%以上。根據(jù)反應菌群多樣性的指數(shù)Simpson和Shannon和菌群豐度指數(shù)Chao和Ace可以看出，細菌豐度最高的是盲腸段，細菌多樣性最高的是大腸段，其中十二指腸是細菌多樣性較高但細菌豐度較低。

表1 4份樣品16S rRNA基因庫多樣性指數(shù)和覆蓋率

Tab.1 Phylotype coverage and diversity estimation of the 16S rRNA gene of four samples

應用軟件mothur version v.1.30.1進行ace、chao、coverage、shannon、simpson 5種常用指數(shù)的測定(表2)，4份樣品的序列覆蓋率均達到99%以上。根據(jù)反應菌群多樣性的指數(shù)Simpson和Shannon和菌群豐度指數(shù)Chao和Ace可以看出，真菌豐度最高的是小腸段，真菌多樣性最高的是大腸段，其中大腸段是真菌多樣性較低同時真菌豐度也較低的腸段。

2.4 菌群分類組成

2.4.1 門分類水平上細菌分類組成

從門的水平上看，金絲猴腸道內(nèi)的細菌主要分屬于7個門類，包括厚壁菌門(Firmicutes)、疣微菌門(Verrucomicrobia)、變形菌門(Proteobacteria)、梭桿菌門(Fusobacteria)、放線菌門(Actinobacteria)、擬桿菌門(Bacteroidetes)、互養(yǎng)菌門(Synergistetes)及其他。其中占有優(yōu)勢的是厚壁菌門、疣微菌門和變形菌門。十二指腸中厚壁菌門占比可高達69%，其他門類細菌豐度很低。從十二指腸到最后的直腸，厚壁菌門所占比例逐漸降低，疣微菌門、變形菌門及梭桿菌門所占比例逐漸升高。

表2 4份樣品ITS基因庫多樣性指數(shù)和覆蓋率

Tab.2 Phylotype coverage and diversity estimation of the ITS gene of four samples

圖1 4份樣品的稀釋性曲線Fig.1 Rarefaction curve of four samples

圖2 4份樣品的稀釋性曲線Fig.2 Rarefaction curve of four samples

圖3 各腸段內(nèi)容物在屬水平上的細菌分布Fig.3 Bacterial distribution of the intestinal bacteria at the level of genus

圖4 各腸段內(nèi)容物在屬水平上的真菌分布Fig.4 Bacterial distribution of the intestinal fungus at the level of genus

2.4.2 屬分類水平上細菌分類組成

從屬的水平上看，金絲猴腸道內(nèi)的細菌主要分屬于31個屬(圖3)，其中Akkermansia、瘤胃菌科未培養(yǎng)菌屬(Ruminococcaceae_uncultured)、大腸-志賀氏菌屬(Escherichia_Shigella)、韋榮氏菌科未分類菌屬(Veillonellaceae_unclassified)以及毛螺旋菌屬(Lachnospiraceae_unclassified)是主要的屬。其中厚壁菌門的瘤胃菌科(Ruminococcaceae)是十二指腸中最主要的細菌。

2.4.3 門分類水平上真菌分類組成

從門的水平上看，金絲猴腸道內(nèi)的真菌主要分屬于4類，子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)、接合菌門(Zygomycota)以及未分類真菌(Fungi-unclassified)。從所占比例上看，子囊菌門具有絕對的優(yōu)勢地位，是金絲猴腸道內(nèi)真菌的最主要的門類。從變化趨勢上看，從十二指腸到直腸，接合菌門和未分類真菌所占比例逐漸升高，而子囊菌門和擔子菌門所占比例略微降低。子囊菌門真菌在腸道內(nèi)比例可達67%～74%，是最主要的真菌。

2.4.4 屬分類水平上真菌分類組成

從屬的水平上看，金絲猴腸道內(nèi)的真菌主要分屬于28個屬(圖4)。從十二指腸到大腸，最主要的4個屬未分類Montagnulaceae、松蘿屬(Usnea)、鬼筆屬(Phallus)、槽枝屬(Sulcaria)所占比例逐漸變大，從31%升至82%，而其他24個屬的真菌所占比例逐漸降低，從54%降至8%。

3 討論

湖北神農(nóng)架國家級自然保護區(qū)的大龍?zhí)侗Ｗo站多年來一直對川金絲猴進行人工投食。以往對投食群的研究主要集中于行為學觀察、社群行為等方面[5-6]。對腸道內(nèi)的菌群的研究，傳統(tǒng)方法較多[7]，通過高通量測序的方法進行研究較少[8]。本次死亡的投食群的川金絲猴屬于老年半野生金絲猴，為我們研究提供了難得的材料。其從開始出現(xiàn)咳嗽現(xiàn)象到最后死亡大概經(jīng)歷半月左右時間。期間食欲逐步下降，藥物治療是在16日投喂了2片乙酰螺旋素片。本研究中川金絲猴腸道內(nèi)菌群多樣性特征代表了老年死亡金絲猴腸道特點。另外，不可忽略的是乙酰螺旋素會對腸道菌群多樣性有一定的影響。乙酰螺旋素是大環(huán)內(nèi)酯類抗生素，對金黃色葡萄球菌、鏈球菌、肺炎球菌以及腦膜炎球菌等有抑菌作用，這將降低陽性球菌所占比例。

根據(jù)本次研究川金絲猴腸道內(nèi)菌群多樣性有如下幾個特點：(1)十二指腸中厚壁菌門細菌占比可高達69%，其他門類細菌豐度很低。厚壁菌門[9]通常是人類和靈長類動物腸道內(nèi)的主要細菌，但是能達到70%左右的還是非常少見。可能因為十二指腸離胃較近，pH較低，不太適宜其他門類細菌的存在。(2)該金絲猴腸道內(nèi)細菌組成的另一個特點是疣微菌門在整個腸道內(nèi)所占比例較高，僅次于厚壁菌門，位居第二，比例在4%～35%之間。金絲猴和人同屬于靈長類，食性都是雜食，腸道菌群結構具有一定相似性。人類的腸道細菌主要是屬于厚壁菌門和擬桿菌門，疣微菌門所占比例甚小。疣微菌門是一門被劃出不久的門類，包括的細菌種類不是很多，具有胞質(zhì)突出形成的兩個到多個突起，因此得名，在人類糞便中常被發(fā)現(xiàn)。疣微菌門的Akkermansia在川金絲猴的腸道比在人類腸道內(nèi)占有豐度要高。(3)子囊菌門[10]是該只金絲猴腸道內(nèi)真菌的最主要的門類，具有絕對的比例優(yōu)勢。從變化趨勢上看，從十二指腸到直腸，接合菌門和未分類真菌所占比例逐漸升高，而子囊菌門和擔子菌門所占比例略微降低。子囊菌門真菌在腸道內(nèi)比例可達67%～74%，是最主要的真菌。以往對金絲猴腸道內(nèi)真菌分布情況研究較少，本研究為老年金絲猴腸道內(nèi)真菌多樣性提供參考。

腸道菌群對動物健康有著重要影響和作用，近年來研究結果很多。人的優(yōu)勢菌門為厚壁菌門、擬桿菌門、變形菌門[11]；大熊貓(Ailuropodamelanoleuca)腸道優(yōu)勢菌門為厚壁菌門、變形菌門、梭桿菌門[12]；老虎腸菌優(yōu)勢菌門為厚壁菌門、梭桿菌門、擬桿菌門[13]。已有研究顯示健康成年川金絲猴糞樣菌群以厚壁菌門和擬桿菌門為主[14]。另外，不同年齡段多樣性也不同，變化趨勢是隨著年齡的增大，多樣性也變高[7]。不同年齡階段的金絲猴腸道菌群在種屬上的類別上差別不大，不同主要體現(xiàn)在比例上。本次研究的川金絲猴與以往測定的腸道菌群多樣性有一定差異，可能是由于健康狀況不同、年齡不同、個體不同以及用過1次乙酰螺旋素片有關系。

[1] 周素容.神農(nóng)架川金絲猴投食群的遺傳多樣性及種群動態(tài)[D].武漢：華中師范大學，2015.

[2] 蔡言安，李冬，畢學軍，等.基于不同測序技術的生物群落結構及功能菌分析[J].中國環(huán)境科學，2016，36(6)：1830-1834.

[3] 曹露芬，殷婷婷，袁振亞，等.基于高通量測序的陜西花馬池和茍池微生物群落多樣性分析[J].生物技術通訊，2016，27(3)：374-380.

[4] 胡海兵，崔立，郭靚驊，等.基于高通量測序技術的冠心病患者腸道菌群多樣性研究[J].上海交通大學學報：農(nóng)業(yè)科學版，2016，34(2)：1-11，19.

[5] 周蕓蕓，楊萬吉，張于光，等.神農(nóng)架川金絲猴人工補食群的遺傳多樣性和親緣關系[J].獸類學報，2015，35(3)：229-240.

[6] 蔚培龍，楊敬元，鮑偉東，等.神農(nóng)架川金絲猴投食群的攻擊行為及等級序列[J].獸類學報，2009，29(1)：7-11.

[7] 簡平，王強，王劍，等.不同年齡段川金絲猴腸道菌群結構差異分析[J].動物營養(yǎng)學報，2015，27(4)：1302-1309.

[8] 許文濤，郭星，羅云波，等.微生物菌群多樣性分析方法的研究進展[J].食品科學，2009，30(7)：258-260.

[9] 韓崗，馬婧，陶士珩，等.厚壁菌門細菌及其祖先蛋白的氨基酸偏好性研究[J].安徽農(nóng)業(yè)科學，2011，39(30)：18408-18410.

[10] 孟慶峰，魏江春.中國霜盤衣屬(文字衣科，子囊菌門)地衣的研究[J].菌物學報，2008，27(4)：525-531.

[11] Wang Tingting，Cai Guoxiang，Qiu Yunping，et al.Structural segregation of gut microbiota between colorectal cancer patients and healthy volunteers[J].The ISME Journal，2012，6(2)：320-329.

[12] Zhu Lifeng，Wu Qi，Dai Jiayin，et al.Evidence of cellulose metabolism by the giant panda gut microbiome[J].Proceedings of the National Academy of Sciences，2011，108(43)：17714-17719.

[13] 田銀平.人工飼養(yǎng)條件下虎腸道菌群結構研究[D].哈爾濱：東北林業(yè)大學，2013.

[14] 王劍，王強，曾東，等.健康和腹瀉川金絲猴糞樣菌群比較分析[J].中國獸醫(yī)學報，2015，35(8)：1232-1238.

Golden snub-nosed monkey； Intestinal flora； Diversity analysis； Provisioned group； Shennongjia

對神農(nóng)架保護區(qū)大龍?zhí)侗Ｗo站1只半野生病死川金絲猴的不同部位的腸道內(nèi)容物進行了微生物多樣性研究。在無菌條件下取十二指腸段、小腸段、盲腸段及大腸段4個部位的腸內(nèi)容物，提取微生物基因組DNA，應用高通量測序技術，對細菌和真菌多樣性進行了詳細的分析，比較了不同部位腸道的細菌和真菌的豐度、分布、多樣性及變化趨勢等。發(fā)現(xiàn)該川金絲猴腸道細菌主要分屬于厚壁菌門(Firmicutes)、疣微菌門(Verrucomicrobia)、變形菌門(Proteobacteria)、梭桿菌門(Fusobacteria)、放線菌門(Actinobacteria)、擬桿菌門(Bacteroidetes)和Synergistetes 7個門，真菌主要分屬于子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)、接合菌門(Zygomycota)以及未分類真菌(Fungi_unclassified)。細菌豐度最高的是盲腸段，細菌多樣性最高的是大腸段，真菌豐度最高的是小腸段，真菌多樣性最高的是大腸段。

Gut Microbiota of a Dead Semi-Wild Golden Snub-Nosed Monkey in a Shennongjia Nature Reserve

Liu Yan1，2Zheng Changming1，2Li Linhai3Zhang Chenglin1，2*Liu Xuefeng1，2Yao Hui4Yu Huiliang4Li Ying1，2Yang Minghai1，2Jia Ting1，2Lu Yan1，2Xia Maohua1，2

(1.Beijing zoo，Beijing，100044，China；2.Beijing Key Laboratory of Captive Related Technologies in Wildlife，Beijing，100044，China；3.Planning and Design Academy of Forest Products Industry，Beijing，100010，China；4.Hubei Golden Monkey Research Center，Shennongjia，442421，China)

We analyzed the gut microbiota from different parts of the intestine of a dead semi-wild golden snub-nosed monkey in Dalongtan Station of Shennongjia National Nature Reserve.Intestinal content was sampled from the duodenum，small intestine(jejunum and ileum)，cecum，and large intestine(colon and rectum)under aseptic condition.Microbial genome DNA were extracted and sequenced by high-output sequencing.Then we analyzed the intestinal bacterial and fungal diversity and compared bacterial and fungal abundance，distribution，diversity and change trends between different parts of the intestinal.The intestinal bacteria of the golden snub-nosed monkey mainly belong to Firmicutes，Verrucomicrobia，Proteobacteria，F(xiàn)usobacteria，Actinobacteria，Bacteroidetes，and Synergistetes.The intestinal fungi mainly belong to Ascomycota，Basidiomycota，Zygomycota，and unclassified fungi.Bacteria abundance was highest in the cecum while bacterial diversity was highest in the large intestine.Fungi were most abundant in the small intestine while fungus diversity was highest in the large intestine.

北京市公園管理中心科研項目：高通量測序技術分析北京動物園大熊貓腸道菌群多樣性，項目編號：zx2015013；國家科技支撐計劃課題(2013BAD03B01)

劉燕，女，35歲，碩士，畜牧師；主要從事野生動物致病細菌分離等方面研究。E-mail：liuyanllh@sohu.com

*通訊作者：張成林，E-mail：zhch6465@263.net

2016-09-02

S852.61

A

修回日期：2016-10-25

發(fā)表日期：2017-05-10

2310-1490(2017)02-194-06