生命科學與工程學院李英研究員團隊在大熊貓腸道菌群研究中取得重要進展
2023年8月14日,，Microbiome（Top 一區(qū),，IF= 15.5）在線發(fā)表了題為“The unique gut microbiome of giant pandas involved in protein metabolism contributes to the host’s dietary adaption to bamboo”的學術論文。我校生命科學與工程學院鄧飛龍博士為本文的第一作者,，我校為第一署名單位,，中國國家大熊貓中心王承東、李德生正高級工程師為共同一作,，我校李英研究員和美國阿肯色大學趙江潮教授為共同通訊作者,。
Deng, F., Wang, C., Li, D. et al. The unique gut microbiome of giant pandas involved in protein metabolism contributes to the host’s dietary adaption to bamboo. Microbiome 11, 180 (2023). 原文鏈接：
大熊貓是全球野生動物保護的旗艦物種。盡管大熊貓屬于食肉目動物,，并擁有典型的食肉類動物的胃腸道特征,，但其飲食主要由高纖維的竹子組成。大熊貓基因組缺乏消化纖維素和半纖維素所需的酶,，因此研究人員普遍認為腸道微生物對大熊貓食性變遷的適應發(fā)揮了重要作用,。本課題組既往研究發(fā)現(xiàn)大熊貓腸道微生物包含少量降解纖維食物/半纖維素的酶，但其含量顯著低于典型的草食動物,。因此,，大熊貓腸道微生物對宿主食性改變的適應性貢獻仍然有待深入探究。
本論文通過深度的illumina二代測序和三代測序數據,，構建了較為完整的大熊貓腸道微生物基因組（MAGs）,；利用宏基因組和宏轉錄組測序技術,，將大熊貓腸道微生物與肉食動物、草食動物和雜食動物的消化道微生物進行比較,，發(fā)現(xiàn)大熊貓腸道微生物的氨基酸代謝顯著區(qū)分于草食動物,；進一步分析發(fā)現(xiàn)非解乳糖鏈球菌(Streptococcus alactolyticus)可能是大熊貓腸道微生物氨基酸代謝的關鍵微生物，小鼠灌胃實驗證明非解乳糖鏈球菌顯著提高小鼠空腸內容物中必需氨基酸水平,。
圖1. 實驗設計概述示意圖
結果解讀
1,、大熊貓腸道微生物宏基因組組裝基因組（MAGs）構建
盡可能完整的重構腸道微生物基因組對于研究腸道微生物功能，尤其是獨特的微生物和功能至關重要,。本論文結合6.3T Illumina 二代宏基因組測序數據和0.3T Nanopore三代宏基因組數據,，構建出610個中高質量的微生物基因組，非冗余基因組408個,，高質量基因組148個,，平均覆蓋76.47%的宏基因組reads。主要以厚壁菌門 (204 MAGs, 50.00%) ,、變形菌門 (139 MAGs, 34.07%), 擬桿菌門 (33 MAGs, 8.09%) and 和放線菌門(21 MAGs, 5.15%)為主,。
圖2. 408個宏基因組組裝基因組(MAGs)的分類注釋和系統(tǒng)發(fā)育樹
2. 大熊貓腸道微生物宏轉錄組分析
既往研究基于宏轉錄組測序技術解析大熊貓腸道微生物功能，但是對微生物的表達情況知之甚少,。論文利用14個圈養(yǎng)大熊貓糞便樣本,，結合宏基因組和宏轉錄組數據分析了大熊貓腸道微生物的表達活躍度。
在種水平,，宏基因組數據顯示的大熊貓腸道微生物的優(yōu)勢種為E. flexneri,，平均相對豐度為27.10%，其次是S. alactolyticus（24.34%）和E. coli（6.08%）,。然而,，通過宏轉錄組數據分析發(fā)現(xiàn)表達豐度最高的物種是S. alactolyticus，平均豐度為45.45%,，其次是E. flexneri（10.49%）和Leuconostoc lactis A（3.66%）,。值得注意的是，盡管根據宏基因組數據,，E. coli是第3豐富的物種,，但根據宏轉錄組數據，其平均豐度只有2.38%,，排名第六。
圖3. 基于宏基因組和宏轉錄組的相對豐度（表達量）Top10物種
3,、宏轉錄組數據分析為大熊貓營養(yǎng)代謝基因表達概況提供了新的見解
課題研究了主要營養(yǎng)物（例如,，碳水化合物，氨基酸）代謝中的微生物代謝途徑,，并基于宏基因組豐度和宏轉錄組表達與草食性（牛和羊）,、雜食性（豬,，老鼠和人）和肉食性（貓）進行了比較。對27個纖維/半纖維/木質素降解相關的基因家族的相對豐度和表達豐度比較,，發(fā)現(xiàn)大熊貓中涉及植物細胞壁降解的微生物基因家族的代謝途徑與貓的更為相似,。肉食動物腸道中這些基因家族的相對豐度與草食動物（ANOSIM，R=0.81,，P=0.001）和雜食動物（ANOSIM,，R=0.93，P=0.001）（圖4A）顯著不同,。而基于宏轉錄組數據的分析,，也顯示出類似的模式，即肉食動物顯著地與雜食動物（ANOSIM,，R=0.37,，P=0.001）和草食動物（ANOSIM，R=0.75,，P=0.001）分離（圖4B）,。包括endo-1,4-beta-xylanase (EC ), cellobiose phosphorylase (EC ), and endo-cellulase (EC ) 等關鍵酶在內都與肉食動物貓相似，低于草食動物和雜食動物,，這表明大熊貓的腸道微生物不具備幫助宿主利用竹子纖維獲得額外能量的潛能,。
圖4. 不同宿主種類腸道微生物組中膳食纖維代謝酶相關基因的組成與表達模式
既往研究表明，竹子中的蛋白質可能是大熊貓的主要能量來源,。因此,，本課題分析了涉及氨基酸代謝的基因家族的豐度和表達，包括61個編碼氨基酸降解酶的基因家族和98個涉及氨基酸合成的基因家族,?；诤昊蚪M數據中的基因豐度的Beta多樣性（Bray-Curtis）顯示，不同的宿主物種在涉及氨基酸降解的基因上有顯著的差異（ANOSIM,，P≤0.05）（圖5A）,。肉食動物（大熊貓和貓）與雜食動物（R=0.616，P=0.001）和草食動物（R=0.576,，P=0.001）差異顯著,，而雜食動物和草食動物差異顯著但R值較低（R=0.203，P=0.002）,。在氨基酸生物合成方面,，基于宏基因組（圖5C）和宏轉錄組（圖5D）數據集的PCoA圖上，大熊貓與其他宿主物種顯著不同（P<0.01）,。
4. 大熊貓腸道微生物的必需氨基酸合成
植物蛋白往往缺乏某些動物所需的必需氨基酸,，為了探索腸道微生物是否幫助宿主合成必需氨基酸，本研究分析了 S. alactolyticus（大熊貓腸道中表達活躍度最高的微生物）中涉及必需氨基酸生物合成的基因的豐度,。根據S. alactolyticus的相對豐度（RA）對圈養(yǎng)大熊貓重新分組為三組：高豐度組[RA ≥?50%,，n =?28],，中豐度組(50% >?RA?≥?25%，n =?29)和低豐度組(RA <?25%,，n =?33),，比較分析顯示，高豐度組的map00300 KEGG途徑(圖6A)和map00400 KEGG途徑(圖6B)的豐度顯著高于中等(map00300: P =?0.027, map00400: P?=?0.0056)和低豐度組(map00300: P =?0.0049, map00400: P?=?0.0014)組,。對于map00290 KEGG途徑(圖6C),，其在高組的豐度顯著高于中等組(P =?0.040)，高豐度組與低豐度組之間差異不顯著,。結果顯示S. alactolyticus可能對大熊貓氨基酸合成有重要的貢獻,。
為了驗證非解乳糖鏈球菌（S. alactolyticus）在腸道中合成必需氨基酸的作用，首先從大熊貓中分離出S. alactolyticus,，然后灌胃給低蛋白和正常蛋白日糧的BALB/c小鼠,。灌胃結束后3周取空腸內容測定代謝組。結果表明,，對于低蛋白飲食的小鼠,，灌胃組的非必需氨基酸（8.64 ×?1010 vs. 9.89?×?1010, P?=?0.034）和必需氨基酸（3.71 ×?1010 vs. 4.83?×?1010, P?=?0.042）的總豐度均顯著高對照組。而正常日糧下,，灌胃組的非必需氨基酸（8.01 ×?1010 vs. 8.72?×?1010, P?=?0.21）和必需氨基酸（3.37 ×?1010 vs. 3.99?×?1010, P?=?0.087）豐度更高,，但沒有統(tǒng)計學意義。對于特定的氨基酸種類,，在低蛋白日糧下,，纈氨酸（P =?0.039）和酪氨酸（P =?0.027）在灌胃組中的豐度顯著更高；而當給予正常蛋白飲食時,，谷氨酸（P =?0.033）在灌胃組中顯著增加,。
本研究重建了總共408個非冗余的宏基因組組裝基因組（MAG），發(fā)現(xiàn)大熊貓的腸道微生物群在蛋白質代謝而非碳水化合物代謝中起到了重要作用,。進一步證明非解乳糖鏈球菌（S. alactolyticus）作為大熊貓腸道微生物最豐富的細菌種類,，參與了大熊貓腸道的蛋白質代謝。本項研究為提高我們對大熊貓腸道微生物多樣性和代謝潛力的認識提供了基礎,，從微生物氨基酸代謝的角度解釋大熊貓腸道微生物對大熊貓食性轉變的作用,，為我們解析大熊貓的適應性進化提供了一個新的視角。李英課題組近10年以來,，對大熊貓腸道菌群的組成,、潛在功能、對宿主適應竹子的作用及其在大熊貓野化訓練及遷地保護中的潛在影響進行了較為系統(tǒng)的研究,，成果在Molecular Ecology,，Frontiers in Microbiology等期刊發(fā)表SCI論文10余篇。
