青海毛驢腸道微生物優勢菌群分析

郭 榮，曹雁行，艾德強，董建寶，張國梁，黃嬌嬌，官卻扎西，索佳佳，孫玉江,6*，劉書琴*

（1.青島農業大學動物科技學院，山東青島 266109；2.東營市現代畜牧業發展服務中心，山東東營 257091；3.青海省海北州剛察縣青海省種羊繁育推廣服務中心，青海海北藏族自治州 812200；4.山東畜牧獸醫職業學院，山東濰坊 261061；5.青海省畜牧總站，青海西寧 810008；6.東營職業學院，山東東營 257091）

馬屬動物是食草性非反芻動物，具有高度分隔的胃腸道（Gastrointestinal Tract，GIT），可利用多種植物纖維[1-2]。驢作為馬屬動物，其與反芻動物最大的區別是屬于后腸發酵，具有發育良好的后腸結構，驢后腸的潛在長度和體積超過4.5 m 和110 L，約為其前腸體積的15.95 倍[3]。

動物腸道中的微生物生態系統相對比較復雜，它們對宿主的生長代謝、營養消化與吸收、免疫力以及抵抗病菌入侵的能力等有著重要影響[4-5]。區域[6]及環境[7]是影響腸道微生物的主要因素。趙金香等[8]研究了干旱沙漠地帶新疆驢盲腸中的細菌多樣性，發現乳球菌科和松毛蟲科在消化粗飼料過程中起關鍵作用。已有研究成果證實地理環境條件對鬣蜥、蝠屬動物和歐洲熒火蟲等腸道微生物的生物多樣性有著非常重要的塑造效果[9-11]。

我國多樣化的生態類型造就了品種高達30 多個豐富的驢種資源，青海毛驢主要分布在青海省海東、海南、海北以及黃南等地的農區和半農半牧區，屬于高原大陸性氣候[12]。不同地區往往有著不同環境條件，而宿主與共生菌如何適應不同地區環境且體內不同部位微生物差異一直是研究的熱點之一[13]?，F階段隨著對腸道微生物研究的日益重視，而對于驢腸道微生物的研究還為之甚少。本文利用高通量16S rRNA 測序的方法，與山東德州驢相比較，分析不同地區（青海和山東）和不同部位（胃、小腸、盲腸和直腸）對青海毛驢腸道菌群的組成和結構差異，對改善青海毛驢胃腸道微生態環境、預防與控制疾病的發生有重要的理論指導。

1 材料和方法

1.1 樣本采集 各選取來自青海?。ㄇ嗪ＣH，海拔＞3 000 m，試驗組）和山東?。ǖ轮蒹H，海拔＜30 m，對照組）的6 個個體（3 公3 母），采樣對象均為成年健康驢（5～10 歲）。青海毛驢為舍飼，以麥草、青稞和豆秸等飼料為主，每日飼喂2～3 kg；德州驢為舍飼，以花生秧、玉米秸和麥糠麥秸等飼料為主，每日飼喂5～6 kg。

驢屠宰后，正確分離胃、小腸和大腸，具體部位確定如下：胃-賁門部至幽門部；十二指腸-胃結合部后10 cm；空腸-十二指腸結合部后10 cm；回腸-回盲口向前10 cm 處；盲腸-呈“逗號”狀，為一盲端；直腸-靠近肛門端。用繩子扎緊各段間狹小接口，水平放置，從各腸道段中部取樣，取等量的十二指腸、空腸和回腸內容物，混合后定義為小腸內容物。對包含內容物的胃、小腸（十二指腸、空腸和回腸）和大腸（盲腸、結腸和直腸）的重量進行稱重記錄，之后立即無菌采集胃、小腸、盲腸和直腸4 個部分的新鮮內容物樣品，無菌收集于2 mL 無菌康寧冷凍管中（Corning Incorporated，New York，USA），做好標記后立即放于液氮中進行速凍，然后放于-80℃低溫冰箱保存；收集完畢后排出內容物，對不含內容物的胃、小腸和大腸的重量再次進行稱重記錄，并測量記錄小腸和大腸的長度。收集的樣本送往Majorbio 生物公司進行16S rRNA 基因組分析。其中青海毛驢的小腸和大腸長度及其比例，胃、小腸和大腸內容物重量（內容物重量=含內容物重量-不含內容物重量）及比例如表1 所示。

表1 青海毛驢胃腸道長度及其內容物重量統計

1.2 16S rRNA 基因測序

1.2.1 總DNA 提取和PCR 擴增 使用 TransGen AP221-02 DNA 試劑盒從腸道內容物樣品中提取微生物群落基因組總DNA。使用1%瓊脂糖凝膠來檢測 DNA 提取物，并使用NanoDrop 2000 紫外可見分光光度計（Thermo Scientific，Wilmington，USA）測定DNA濃度和純度。通過ABI GeneAmP? 9700 PCR 熱循環儀（ABI，CA，USA）對細菌16S rRNA 基因的高變區V3-V4 進行擴增用，用的引物對為338F（5'-ACTCCTACGGGA GGCAGCAG-3'）和806R（5'-GGACTACHVGGGTW TCTAAT-3'）。16S rRNA 基因的PCR 擴增流程為：95℃3 min；95℃30 s，55℃30 s，72℃45 s，共27 個循環；72 ℃10 min，4 ℃終止。PCR 混合物包含5 ×TransStart FastPfu 緩沖液4 μL、2.5 mmol/L dNTPs 2 μL、正反向引物（5 μmol/L）各0.8 μL、TransStart FastPfu DNA 聚合酶0.4 μL、模板DNA 10 ng，最后ddH2O 補充至20 μL。PCR 反應一式三份進行。PCR 產物從2%瓊脂糖凝膠中提取，并使用AxyPreP DNA 凝膠提取試劑盒（Axygen Biosciences，Union City，CA，USA）根據使用說明進行純化，并使用QuantusTM熒光計（Promega，USA）進行定量。

1.2.2 Illumina MiSeq 測序 根據美吉生物公司（中國上海）的標準方案，在Illumina MiSeq PE300 平臺（Illumina，San Diego，USA）上以等摩爾和雙端方式匯集純化的擴增子。

1.2.3 統計分析 原始16S rRNA 基因測序讀數被解復用，由FASTQ 版本0.20.0[14]對其質量進行過濾，用FLASH 版本1.2.7[15]對其合并，標準如下：在任何位點的300 bp 讀數在50 bp 滑動窗口上被截斷，從而獲得＜20 的平均質量分數，并且小于50 bp 的截斷讀數被棄去，包含歧義字符的讀數同樣被舍棄；只有超過10 bp的重疊序列才會根據它們的重疊序列進行組裝。重疊區域的最大失配率為0.2。棄去無法組裝的讀數；區分樣品的方法是根據條形碼和引物，調整序列方向，條形碼精確匹配。

使用 UPARSE 7.1 版[16]對具有97%相似性截止值的操作分類單元（OTU）[16-17]進行聚類，并識別并刪除嵌合序列。每個OTU 代表序列的分類通過RDP 分類器版本2.2[18]針對16S rRNA 數據庫（Silva v138）使用0.7的置信閾值進行分析。

2 結果

2.1 稀釋曲線分析 圖1 顯示，在8 個分組中，Shannon指數曲線上升后均逐步趨于平緩，說明測序數據達到飽和，能夠覆蓋腸道微生物組群落的絕大部分物種，結果可用于后續數據分析。

圖1 Shannon 指數曲線

2.2α多樣性分析 在OTU 水平，青海毛驢的胃、小腸、盲腸和直腸的微生物物種豐富度和多樣性均高于德州驢。其中，小腸的豐富度和多樣性最低，盲腸和直腸的豐富度和多樣性均高于胃和小腸。

Shannon 指數表明（圖2-A），青海毛驢胃的物種多樣性低于其盲腸和直腸（P≤0.05），青海毛驢胃的物種多樣性高于德州驢小腸（P≤0.01）；青海毛驢小腸的物種多樣性低于其盲腸和直腸（P≤0.01），而青海毛驢小腸的物種多樣性低于德州驢盲腸和直腸的物種多樣性（P≤0.05）；青海毛驢盲腸和直腸的物種多樣性高于德州驢各組（P≤0.05）；德州驢小腸的物種多樣性顯著低于其盲腸和直腸（P≤0.01）。

Ace 指數表明（圖2-B），青海毛驢胃的物種豐富度低于其直腸（P≤0.05）；青海毛驢盲腸的物種豐富度顯著低于青海毛驢直腸，青海毛驢盲腸的物種豐富度高于德州驢胃（P≤0.05）和小腸（P≤0.01）的物種豐富度；青海毛驢直腸的物種豐富度高于德州驢胃和盲腸（P≤0.01）以及小腸（P≤0.001）的物種豐富度；德州驢胃的物種豐富度低于其直腸（P≤0.05）；德州驢小腸的物種豐富度低于其盲腸（P≤0.05）和其直腸（P≤0.01）；德州驢盲腸的物種豐富度低于其直腸的物種豐富度（P≤0.05）。

圖2 α 多樣性分析

2.3 基于門和屬水平的微生物群落結構組成 在門水平上，共鑒定到20 個門水平菌，其中共有門18 個（90%），青海毛驢小腸有2 個獨有菌門為P__Caldisericota（50%）和氫化戊烯菌（P__Hydrogenedentes，50%）（圖3-A）。

青海毛驢前腸的優勢菌門為變形菌門（Proteobacte ria，36.01%）和厚壁菌門（Fimicutes，32.68%），青海毛驢后腸的優勢菌門為厚壁菌門（Fimicutes，56.89%）和擬桿菌門（Bacteroidota，29.19%）。德州驢前腸的優勢菌門為厚壁菌門（Fimicutes，89.20%）和擬桿菌門（Bacteroidota，6.37%），德州驢后腸的優勢菌門為厚壁菌門（Fimicutes，54.87%）、擬桿菌門（Bacteroidota，36.28%）和螺旋菌門（Spirochaetota，5.12%）（圖3-C）。

在屬水平上，共鑒定到276 個屬水平菌，其中共有屬125 個（45.29%），青海毛驢胃獨有77 個菌門，占比前兩個的為黃桿菌屬和噬菌體屬（Flavobacterium24.19%，Algoriphagus9.72%）；青海毛驢小腸獨有19個菌門，占比前2 個的為普羅威登斯菌屬和原衣原體屬（Providencia65.97%，_Candidatus_Protochlamydia4.17%）；青海毛驢盲腸獨有2 個菌門，占比前兩個的屬為琥珀弧菌屬和副桿菌屬（Succinivibrionacia_UCG-00250%，Paracaedibacteracia50%）；青海毛驢直腸獨有5 個菌屬，占比前兩個的為DMER64 和脫硫孢屬（DMER6434.78%，DesulfurisPora30.43%）；德州驢胃獨有37 個菌屬，占比前兩個的為巨型單胞菌屬和桿菌屬（Megamonas19.52%，Coriobacteriacia8.5%）；德州驢小腸獨有9 個菌屬，占比前兩個的為脫甲基菌屬和乳酸菌屬（Demequinacia29.99%，Lacticigenium16.67%）；德州驢盲腸獨有1 個菌屬為g__norank_f__NS11-12_marine_group；德州驢直腸獨有1 個菌屬為奧秘桿菌屬（Arcanobacterium）（圖3-B）。

青海毛驢前腸的優勢菌屬為鏈球菌屬（Streptococc us，14.12%）、放線菌屬（Actinobacteria，16.62%）和乳酸桿菌屬（Lactobacillus，6.31%），青海毛驢后腸的優勢菌屬為理研菌屬（Rikenellacia，10.11%）密螺旋體（Treponema，6.60%）和乳酸桿菌屬（Lactobac illus，5.64%）；德州驢前腸的優勢菌屬為乳酸桿菌屬（Lactobacillus，61.13%）和鏈球菌屬（Streptococcus，14.22%），德州驢后腸的優勢菌屬為鏈球菌屬（Strepto coccus，15.77%）和理研菌屬（Rikenellacia，8.73%）（圖3-D）。

圖3 腸道微生物群落結構組成

2.4 Beta 多樣性分析

2.4.1 樣本層級聚類分析 試驗樣品的菌群組成關系和結構利用相似度樹狀圖描述。首先使用描述群落組成關系和結構的Weighted-unifrac 算法計算在OTU 水平各個樣本間的距離，即根據 beta 多樣性距離矩陣進行層次聚類分析，得到相似度樹狀圖，結果如圖4 所示，青海毛驢胃和小腸的菌群組成結構相似，成簇聚類到同一樹杈分支點上；德州驢胃和小腸菌群具有不同的聚落結構，成簇聚類到另一樹杈分支點上。青海毛驢直腸和德州驢直腸的菌群分別聚類在樹杈分支點上，青海毛驢盲腸和德州驢盲腸的菌群分別聚類在樹杈分支點上，從而形成鮮明的菌群結構派別。

圖4 樣本層級聚類分析

2.4.2 PCoA 分析 利用主坐標（Principal co-ordinates analysis，PCoA）比較微生物的組成，結果表明，青海毛驢和德州驢的地區差異顯著，且來自相同/相鄰區域（胃和小腸，盲腸和直腸）的微生物群落比來自其他區域的微生物群落更相似（圖5）。

圖5 PCoA 分析

2.4.3 組間差異顯著性檢驗 在門水平上，厚壁菌門（Fimicutes）的豐度在德州驢中高于青海毛驢（P=0.002），在德州驢小腸中的豐度最高，在德州驢胃中的豐度高于青海毛驢小腸（P=0.037），在德州驢胃中的豐度高于青海毛驢直腸（P=0.04），在德州驢胃中的豐度高于青海毛驢胃（P=0.0019）。擬桿菌門（Bacteroidota）在青海毛驢和德州驢中，盲腸和直腸的豐度均高于胃和小腸（P=0.0002）。變形菌門（Proteobacteria）在青海毛驢的豐度高于德州驢（P=0.003），且在青海毛驢胃和小腸的豐度高于盲腸和直腸。藍細菌（Cyanobacteria）只出現在青海毛驢中（P=0.0001），且只存在胃中。螺旋菌門（Spirochaetota）在青海毛驢和德州驢的盲腸和直腸的豐度高于胃和小腸（P=0.0006），在德州驢中，在盲腸中的豐度高于直腸，且在青海毛驢小腸的豐度高于德州驢小腸（P=0.048），在青海毛驢盲腸的豐度高于德州驢直腸（P=0.04）。放線菌門（Actinobacteria）在青海毛驢的豐度高于德州驢（P=0.001），且在青海毛驢中胃和小腸的豐度高于盲腸和直腸；在青海毛驢盲腸的豐度高于德州驢盲腸（P=0.0007）。乳頭菌門（Patescibacteria）在青海毛驢的豐度高于德州驢（P=0.008）（圖6）。

圖6 在門水平進行組間差異顯著性檢驗

在屬水平上，乳酸菌屬（Lactobacillus）在德州驢中的豐度高于青海毛驢（P=0.000006），且德州驢的平均相對豐度在胃和小腸中高于盲腸和直腸且高于青海毛驢的胃和小腸，青海毛驢的盲腸和直腸的乳酸菌豐度較低。鏈球菌屬（Streptococcus）在德州驢的豐度高于青海毛驢（P=0.014），且在德州驢盲腸的豐度顯著高于青海毛驢直腸（P=0.03），在青海毛驢胃的豐度顯著高于青海毛驢直腸（P=0.04）。理研菌屬（Rikenellacia）在盲腸和直腸的豐度高于胃和小腸，在青海毛驢直腸中的豐度高于青海毛驢盲腸（P=0.04）。放線菌屬（Actinobacteria）在青海毛驢胃和小腸的豐度較高，在德州驢的豐度較低。葉綠體（Chloroplast）只存在于青海毛驢的胃中（P=0.001）。密螺旋體屬（Treponema）在兩組中盲腸和直腸的豐度高于胃和小腸，在德州驢小腸的豐度較低。毛螺菌屬（Lachnospiracia）在青海毛驢的豐度高于德州驢，在盲腸和直腸中的豐度高于胃和小腸，在青海毛驢的盲腸豐度較高，在青海毛驢盲腸中的豐度高于青海毛驢直腸（P=0.0017）。產氣莢膜桿菌（Clostridium）在德州驢的豐度高于青海毛驢（P=0.006）（圖7）。

圖7 在屬水平進行組間差異顯著性檢驗

2.4.4 LEfSe 多級物種差異判別分析 通過LEfSe 分析進行組別間差異比較，找到與豐度有顯著性差異的類群，如圖8 所示，對樣本劃分產生顯著差異影響的群落或物種含量前三的細菌在青海毛驢胃中是藍細菌（Cyanobacteria）、葉綠體（Chloroplast）和放線桿菌門（Actinobacteria）；青海毛驢小腸中是變形菌門（Proteobacteria）、γ-變形菌綱（Gammaproteobacteria）、巴氏桿菌科（Pasteurellaceae）；青海毛驢盲腸中是梭狀芽孢桿菌綱（Clostridia）、乳螺目（Lachnospirales）、振螺目（Oscillospirales）；青海毛驢直腸是理研菌科（Rikenellaceae）、毛螺菌科_AC2044（Lachnospiraceae_AC2044）、疣菌群屬（Verrucomicrobiota）。

圖8 LEfSe 多級物種層級樹圖

德州驢胃中是伯克氏菌目（Burkholderiales）、伯克氏菌科（Burkholderiaceae）、糞桿菌屬（Faecalibacte rium）；德州驢小腸是乳桿菌目（Lactobacillales）、桿菌（Bacilli）、厚壁菌門（Firmicutes）；德州驢盲腸是普雷沃菌科（Prevotellaceae）、螺旋體屬（Spiroch aetota）、密螺旋體屬（Treponema）；德州驢直腸是擬桿菌目（Bacteroidales）、擬桿菌屬（Bacteroidia）、擬桿菌屬（Bacteroidota）。

2.4.5 功能預測分析 通過PICRUSt 對OTU 豐度表進行標準化，然后對OTU 進行COG 功能注釋，獲得OTU 在COG 功能水平的注釋信息及各功能在不同樣本中的豐度信息。圖9 顯示，除去COG 的功能中比重最高的功能未知基因（Function unknown）和一般功能預測基因（General Function Prediction only）這兩大類別，其中碳水化合物轉運與代謝基因（Carbohydrate Transport and Metabolism）除了青海毛驢胃外，在青海毛驢的小腸、盲腸和直腸以及德州驢4 個部位均呈現高豐度分布狀態；氨基酸轉運與代謝基因（Amino Acid Transport and Metabolism）除了德州驢小腸外，在青海毛驢的4 個部位以及德州驢的胃、盲腸和直腸均呈現高豐度分布狀態其余樣品中均呈現高豐度分布狀態，說明腸道中碳源物質代表碳水化合物和氮源物質代表氨基酸兩種物質代謝旺盛。

圖9 COG 功能分類統計圖

3 討 論

3.1 在門水平的微生物群落結構組成 微生物群落結構組成在門水平上的結果顯示，青海毛驢前腸的優勢菌門為厚壁菌門和變形菌門，青海毛驢后腸的優勢菌門為厚壁菌門和擬桿菌門。厚壁菌門的物種豐度在德州驢中高于青海毛驢，在青海毛驢中的豐度后腸高于前腸，而在德州驢中的豐度前腸高于后腸。劉桂琴等[19]基于16S rRNA 測序對德州公驢的研究表明，驢前腸優勢菌門順序為厚壁菌門、變形菌門、擬桿菌門和放線菌門；后腸為厚壁菌門、擬桿菌門、螺旋體門和變形菌門。劉宇等[20]發現驢盲腸菌群優勢順序為厚壁菌門、擬桿菌門、疣微菌門和螺旋體門。有報道稱厚壁菌門在草食動物中主要是利用碳水化合物[21]，有提高動物免疫力和增強腸道功能的作用。

在馬屬動物中，前腸主要負責食物的消化和吸收，后腸與微生物發酵有關[22]。雖然在馬的胃和小腸中已證實有輕微的發酵，但大腸仍是飼料發酵的最重要部分[23-24]。已知擬桿菌擁有豐富的編碼碳水化合物活性酶的基因，能夠根據可用性在腸道中不同類型的能量源之間輕松切換[25]。擬桿菌在大腸的消化相關微生物群中顯著富集，并且不受粘膜相關微生物群中腸道位置的影響。李志鵬等[26]在梅花鹿瘤胃中菌群的研究發現，其中87.9%的16S rRNA 基因序列都屬于擬桿菌門降解纖維的菌種序列。因此，擬桿菌在大腸中富集可能有助于宿主適應復雜的內部環境和大腸中纖維性碳水化合物的降解。之前已有研究證明，厚壁菌門和擬桿菌門是山羊[27]、牛[28]和馬[29]等食草動物中最具優勢的2 個門。同樣，Chen 等[30]研究發現，肉牛瘤胃中厚壁菌門和擬桿菌門的含量隨飼糧干草含量的增加而增加，特別是厚壁菌門。以上研究成果和本試驗結果均說明，厚壁菌門和擬桿菌門可能與動物個體飲食特征和消化大量纖維性碳水化合物有關。瘤胃疣微菌屬于典型的纖維降解菌，與纖維降解密切相關，可產生纖維酶降解纖維二糖等纖維類物質[31]。Patra 等[32]和Zhao 等[33]均證明瘤胃菌科微生物的數量與纖維素消化利用存在密切關系。由此可以看出，驢后腸具有極強的纖維素分解微生物基礎，可以很好地解釋其耐粗飼的原因。

3.2 在屬水平的微生物群落結構組成 在屬的水平上觀察細菌組成時，青海毛驢前腸的優勢菌屬為鏈球菌屬、放線菌屬和乳酸桿菌屬，青海毛驢后腸的優勢菌屬為理研菌屬、密螺旋體和乳酸桿菌屬。在青海毛驢中，乳酸桿菌在胃內和小腸內的平均相對豐度高于盲內和直內，表明乳酸桿菌是青海毛驢前腸的主要屬。據報道，乳酸桿菌具有非纖維性碳水化合物降解能力（如戊糖、己糖和淀粉）[34]，參與膽汁鹽的解偶聯[35]，并通過產生抗微生物物質（如細菌素和乳酸鹽）或與病原菌競爭粘膜粘附位點和營養物來抑制病原菌的增殖[36-37]。劉桂琴等[19]在驢腸道微生物多樣性及代謝功能差異分析的研究表明，前腸中乳酸桿菌屬為優勢菌群，約占50.1%，其次為鏈球菌屬，占21.2%。Liu 等[38]研究表明，在屬水平上，乳酸桿菌在前腸中占優勢，而鏈球菌在后腸中占優勢。淀粉分解細菌-鏈球菌總是存在于飼喂高濃度日糧的反芻動物的瘤胃中[39]。與乳酸桿菌屬相似，其發酵產物乳酸鹽使腸粘膜成為酸性環境，因此，乳酸桿菌和鏈球菌都被認為是對馬屬動物有益的益生菌[40-41]。

3.3α多樣性分析 在細菌門水平上，變形菌門、放線菌門和疣菌群門在青海毛驢的占比高于德州驢；在屬水平上，乳酸菌屬在德州驢高于青海毛驢；放線菌屬在青海毛驢的胃內和小腸內，德州驢幾乎沒有；毛螺菌屬在青海毛驢中的豐富度高于德州驢，葉綠體只在青海毛驢組出現，且只存在胃內中。葉綠體在青海毛驢的胃里分布的比較多，可能與青海地區海拔高，植被光合作用較強以及飼料因素有關，所以飼料中的含量相對較高。

青海地區比德州地區海拔高，隨著海拔高度的增加溫度會隨之下降，而腸道菌群的物種組成易受環境溫度的影響，環境溫度降低會導致菌群α多樣性降低[42-43]。動物實驗中也發現了隨著溫度的降低，小鼠腸道菌群的α多樣性會隨之下降。Zietak 等[44]研究發現，長期生活在12℃的小鼠腸道菌群的α多樣性明顯低于生活在29℃的小鼠。但也有研究發現環境溫度對宿主腸道菌群結構的影響也較小，如Li 等[45]研究發現生活在溫暖環境（22℃）的小鼠腸道菌群結構與生活在寒冷環境（4℃）的小鼠腸道菌群結構相似。隨著海拔高度增加，宿主生活環境中的氧濃度不斷下降，但低氧并不是影響腸道菌群結構的關鍵因素[46]；如Lucking 等[46]研究發現，長期處于低氧狀態下小鼠腸道菌群的群落結構與正常小鼠相似；Khanna 等[47]的研究也證實宿主處于低氧環境時腸道菌群結構未發生改變。梁田等[48]研究也發現生活在不同海拔高度的藏族人群腸道菌群結構相似，無明顯差異，因此推測低溫和低氧等環境因素可能對藏族人腸道菌群結構的影響較小。在趙俊松等[49]的研究中西藏恒河猴種群的腸道微生物組成多樣性和種群特有菌群種類均高于其他5 個地理種群（云南大理、貴州貴陽、四川西昌、廣西南寧、海南?？冢?，即西藏種群中腸道微生物菌群的多樣性最高，且擁有最多種類的特有菌群；同時，青藏高原地區的低溫、低氧和強紫外線的高海拔環境也有可能對恒河猴腸道微生物的組成多樣性產生影響[50]，幫助恒河猴適應青藏高原地區的高海拔環境[51]。Liu 等[52]通過高通量測序V3-V4 區16s rRNA 基因來分析藏野驢（中國青海省海北州高寒草甸）和非洲野驢（中國青海省玉樹州沼澤草甸）的腸道微生物群，結果表明，藏野驢的干物質消化率顯著高于非洲野驢，兩組之間的α多樣性沒有顯著差異，但β多樣性差異顯著，藏野驢比非洲野驢擁有相對更復雜的細菌網絡，總體結果表明，藏野驢在干物質消化、腸道微生物群落組成和功能方面均優于非洲野驢。本文中青海毛驢地處高原，青海毛驢的菌群多樣性和豐富度均高于德州驢，說明青海毛驢有其獨特的腸道菌群以適應高原的低溫、低氧以及高海拔。但是目前，低溫和低氧對于腸道菌群結構的具體作用機制尚不清楚，需要更多的研究與數據支持。

3.4 功能預測 為探索微生物群的潛在功能，本實驗使用PICRUSt 來確定不同腸道位置相關微生物群的潛在功能，結果顯示，碳水化合物轉運與代謝基因（Carbo hydrate Transport and Metabolism）除了在青海毛驢胃外，在其余各組均呈現高豐度分布狀態；氨基酸轉運與代謝基因（Amino Acid Transport and Metabolism）除了德州驢小腸外，其余樣品中均呈現高豐度分布狀態，說明腸道中碳源物質代表碳水化合物和氮源物質代表氨基酸兩種物質代謝旺盛。這與其他動物的腸道微生物群的功能一致[25-26]，也對應于微生物群的基礎代謝功能及其生長和增殖所需的營養物質（如碳源和氮源）[53]。

4 結 論

本研究主要基于16S rRNA 基因高通量測序研究了青海毛驢的不同部位腸道微生物的差異。結果表明，青海毛驢的菌群多樣性和豐富度均高于山東德州驢，且后腸的菌群多樣性和豐富度均高于前腸。這些發現對青海毛驢的腸道微生物群有了更多的了解，為進一步的研究奠定基礎。