其他
科研 | Gut:多组学分析揭示了胎羊肠道中存在微生物组(国人佳作)
编译:微科盟草重木雪,编辑:微科盟Tracy、江舜尧。
微科盟原创微文,欢迎转发转载。
论文ID
原名:Multiomics analysis reveals the presence of a microbiome in the gut of fetal lambs译名:多组学分析揭示了胎羊肠道中存在的微生物组
期刊:Gut
IF:19.819发表时间:2021.05通讯作者:刁其玉
通讯作者单位:中国农业科学院饲料研究所
实验设计
实验结果
(A)胎儿盲肠内容物样本数据去污前后微生物宏基因组(MG)和元转录组(MT)总基因计数。Before表示数据去污前的基因计数;After数据去污后的基因计数。(B)数据去污前后每个盲肠内容物样品中微生物MG和MT的平均基因计数。在(B)中,方框代表第一至第三四分位数之间的IQRs,方框内的水平线代表中位数;须表示距离第一或第三四分位数1.5 IQR内的最小值或最大值。 2. 宏基因组分析确定的分类组成 为了在宏基因组水平上鉴定胎儿肠道的微生物群落结构,我们对6只无菌子宫切除术分娩羔羊的盲肠内容物样本进行了宏基因组测序,将宏基因组基因与国家生物技术信息中心(NCBI)的nr核苷酸数据库进行比对,以进行分类分配。微生物类群被定义为至少存在于六分之三的样本中。我们发现,胎儿羔羊肠道显示出较低的α多样性(在线补充图S2),表明胎儿肠道中的细菌丰富度较低。盲肠内容物样本中鉴定出四个主要门:变形杆菌(95.30%±2.19%)、厚壁菌(4.85%±1.71%)、放线杆菌(0.53%±0.22%)和奇古菌门(0.02%±0.01%)(图2A,在线补充表S7)。在属和种水平上,盲肠内容物样品中检出33属50种,大肠杆菌(88.76%±2.04%)和链球菌(4.19%±1.65%)是主要的细菌属(图2B,在线补充表S7),而大肠杆菌(86.89%±2.21%)和海洋链球菌(4.11%±1.61%)是其中最丰富的物种(图2C,在线补充表S7)。我们还在所有盲肠内容物样本(在线补充图S3A)和胎儿肠道中检测到高水平的噬菌体phiX174(52.29%±5.55%)和Orf病毒(0.03%±0.01%)(在线补充表S8)。
(A-C)分别代表胎儿肠道微生物组的门级、属级和种级细菌组成(相对丰度为0.1%)。微生物组组成显示所有样本中细菌的平均相对丰度。n=6。(D)平均微生物相对丰度与MG (MT/MG)的比值。 3. 宏转录组分析确定的分类组成
为了检查检测到的微生物的潜在活性,我们对盲肠内容样品中提取的RNA进行了转录组分析。宏转录组数据的分类结果显示总共26个细菌属和32个细菌种。四个菌门包括变形菌门(94.73%±3.66%),放线菌(2.75%±2.52%),厚壁菌(1.69%±1.13%)和拟杆菌(0.67%±0.62%)(图2A,在线补充表S9)。在属和种水平上,大肠埃希菌(64.23%±19.30%)和沙门氏菌(15.51%±6.84%)是主要细菌属(图2B,在线补充表S9),而大肠杆菌(64.14%±19.17%),肠沙门氏菌和沙门氏菌(15.60%±6.87%)是最丰富的物种(图2C,在线补充表S9)。盲肠内容物样品的宏转录组分析还表明,噬菌体phiX174的水平很高(占所有映射reads的58.88%±14.73%)(在线补充图S3B)。为了推断样品收集时的微生物相对活性,我们基于微生物功能基因转录与基因的比率分析了样品中微生物基因的表达。MT/MG比率定义为宏转录组中平均相对微生物丰度的比率,放线菌门的成员,例如分枝杆菌属和角质杆菌属,显示出增加的MT/MG比(图2D);而厚壁菌门的成员,例如链球菌属,表现出降低的MT/MG比率,但是梭菌的MT/MG比率增加。变形杆菌的MT/MG比率接近1,但属于该门的成员差异很大,例如,大肠埃希氏菌和柠檬酸杆菌属的成员显示降低的MT/MG比,而属于变形杆菌的其他微生物显示增加的MT/MG比。 4. 盲肠内容物样品和对照样品中的微生物定量
我们进行了绝对定量定量PCR(qPCR),以定量总细菌和五个选定微生物(图3)的拷贝数,这些微生物在我们宏基因组学数据中显示出较高的相对丰度。盲肠含量样品和阴性对照中每克总细菌的拷贝数(表示为平均值±SD)分别为4.6×107±3.4×107和1.6×107±1.1×106。盲肠内容物样品显示每克总细菌的拷贝数显著高于阴性对照(p<0.05),而大肠杆菌,海洋杆菌,肺炎链球菌和酿酒酵母的拷贝数显著高于阴性对照(p <0.05),而这两类样本的慢生根瘤菌每克拷贝数没有显著差异。
方框表示第一和第三四分位数之间的IQRs,方框内的水平线表示中位数;须表示从第一或第三四分位数开始的1.5× IQR内的最小值或最大值。胡须上方不同字母的盒数有显著差异(p<0.05)。 5. 盲肠内容物微生物组的功能特征 为了在功能上表征盲肠微生物组,我们针对《京都议定书》的基因和基因组百科全书(KEGG)和基因的进化谱系以及非监督直系同源群体(eggNOG)数据库对宏基因组和宏转录组的nr基因目录进行了比对。在宏基因组基因中,我们用KEGG和eggNOG分别对5458和8749个基因进行了注释。隶属于信号转导、碳水化合物代谢、氨基酸代谢和能量代谢的KEGG通路,在胎儿肠道微生物组中高度富集(图4A)。同样地,在eggNOG功能类中,涉及氨基酸运输与代谢、碳水化合物运输与代谢、能量生产与转换、信号转导机制途径的微生物基因也高度富集。此外,涉及转录、翻译后修饰、蛋白转换、伴侣、细胞内转运、分泌和泡囊运输的通路在eggNOG功能类别中也高度富集(图4C)。我们使用KEGG和eggNOG分别对宏转录组数据中912和1168个基因进行了注释。涉及碳水化合物代谢、能量代谢、氨基酸代谢、膜转运和信号转导的KEGG通路在胎儿肠道微生物组中高度富集(图4B)。eggNOG在胎儿肠道中最活跃的功能类别是能量产生与转换、氨基酸运输与代谢、碳水化合物运输与代谢,其次是翻译、核糖体结构与生物发生、细胞壁/膜/包膜生物发生(图4D)。
(A,B)分别为胎儿肠道元基因组和元转录组的KEGG分类。(C,D)分别为胎儿肠道元基因组和元转录组的eggNOG功能分类。eggNOG,非监督直系同源群体;KEGG,《京都基因与基因组百科全书》。 6. 盲肠内容物和脐血的代谢组学分析
为了进一步分析胎儿肠道中的微生物活性,我们进行了代谢组学分析,以检测样本中的微生物代谢产物。我们还测定了短链脂肪酸(SCFA)的浓度,并将其作为肠道微生物代谢活指标。在本研究中,我们使用的乙酸,丙酸和异丁酸的检出限为5×10-5 µg/mg,丁酸,异戊酸,戊酸和己酸的检出限为5×10-6 µg/mg。此范围允许检测盲肠内容物样品中除丙酸,异丁酸和异戊酸之外的所有SCFA。结果表明乙酸,丁酸,戊酸和己酸的浓度分别为1.88×10-3±4.05×10-4、1.01×10-4±2.43×10-5、2.71×10-5±4.56×10-6、1.02×10−4±3.01×10−5µg/mg。为了追踪肠道中这些SCFA的可能来源,我们检查了脐带血中的SCFA。脐带血中乙酸,丙酸,丁酸,异丁酸,戊酸,异戊酸和己酸的浓度分别为6.23±2.94、1.69±0.55、0.72±0.24、0.56±0.19、0.24±0.09、0.04±0.03和0.41±0.24 µg/mL(在线补充表S10)。我们还使用高效液相色谱-四级杆飞行时间质谱(UHPLC-QTOF-MS)进行了非靶向代谢组学分析,以检测样品中的其他代谢物。我们从盲肠内容物样品中共检测到2808个峰,其中正离子模式(ESI+)和负离子模式(ESI-)的峰分别为1271个和1537个。使用人类代谢组数据库(HMDB)、PubChem和KEGG数据库对这些峰进行比对,我们能够注释1990个总峰(ESI+中937个峰和ESI-中1053个峰)(在线补充表S11和S12),在这1990个代谢物中,有151个被注释到180个KEGG途径,而其他1839个代谢物没有被指定为特定的途径,然后将这些KEGG途径分为寄主相关(即绵羊)途径和非寄主相关途径(即所有与绵羊不共享的KEGG途径)。在这151种代谢物中,143种参与了127种与宿主相关的KEGG途径(在线补充表S13),77种参与了53种与非宿主微生物相关的KEGG途径(在线补充表S14)。53条非宿主微生物相关代谢途径包含8种独特的代谢产物,包括脱氧野霉素、丝裂霉素、妥布霉素、对苯醌、氨基磺酸环己酯、大豆苷元、氯仿和橙皮素。我们注意到这8种代谢物都与已知的微生物代谢途径有关。具体来说,脱氧野霉素、丝裂霉素和妥布霉素参与“抗生素生物合成”,妥布霉素参与“新霉素、卡那霉素和庆大霉素生物合成”,对苯醌和环己基氨基磺酸参与“不同环境中的微生物代谢”,大豆苷元和氯代物参与“次级代谢物的生物合成”,橙皮苷参与“类黄酮生物合成”(在线补充表S14)。为了追踪肠道中这些微生物代谢物的可能来源,我们进行了非靶向代谢组学分析以检测脐血中的代谢物,但没有发现脱氧野霉素、丝裂霉素或妥布霉素的存在(在线补充表S15和S16)。 7. 盲肠微生物组和盲肠内容物代谢物之间的关系 为了确定盲肠微生物组与盲肠含量代谢物之间的相关性,我们使用Spearman秩相关系数进行了共现网络和热图分析。在共存网络和热图中,丁酸和己酸的浓度与肺炎链球菌和科氏菌属的相对丰度呈正相关(Spearman的相关值> 0.6,p <0.05;图5A,B)。此外,己酸浓度与肺炎克雷伯菌的相对丰度呈正相关,炔诺酮浓度与肺炎克雷伯菌和肺炎链球菌的相对丰度呈正相关,但与肠沙门氏菌的丰度呈负相关。除了确定盲肠微生物组和盲肠内容物代谢物之间的相关性外,我们还发现盲肠内容物宏基因组,宏转录组和代谢组共享许多KEGG途径。在宏基因组学和宏转录组数据中,分别有376条和198条KEGG途径以及代谢组共享152条和95条KEGG途径(在线补充表S17和S18)。
(A)胎儿肠道转录组中盲肠含量代谢物和微生物物种之间共现网络的分析。节点代表代谢物和微生物种类。通过节点的线数增加会增加节点的大小。红色边缘表示代谢物与微生物种类之间呈正相关;Spearman等级相关系数> 0.6,p <0.05。蓝色边缘表示代谢产物与微生物种类之间呈负相关;Spearman等级相关系数<-0.6,p <0.05。(B)胎儿肠道转录组中盲肠含量代谢物和微生物种类之间Spearman等级相关系数的热图。Spearman等级相关系数> 0.6或<-0.6,* p <0.05,** p <0.01。红色和蓝色分别表示正相关和负相关。颜色强度与Spearman等级相关值成正比。
讨论
结论
综上所述,以无菌子宫切除术分娩的羔羊为动物模型,通过多组学分析,此研究提供了强有力的证据,证明产前肠道具有低多样性和低生物量的代谢活性微生物群,胎儿肠道微生物群是在产前定植的。这些发现促进了我们对胎儿肠道微生物群的理解,这有助于设计旨在改善人类健康和治疗与肠道微生物群失调相关的非传染性疾病的临床疗法。然而,还需要进一步的研究来深入了解胎儿肠道微生物群的起源、组成、功能、动力学和定植时间,以及阐明胎儿肠道微生物群及其代谢物在整个生命早期对胎儿发育和健康的影响。
----------微科盟更多推荐----------
免费生信作图平台——生科云 | |
长按左侧二维码 进入生科云 | |
生科云所有分析工具可以免费使用,不收取任何直接或间接费用;您还可以在微信上联系微生态老师,随时获取免费的指导,帮助您解决分析时遇到的问题;专业的生信分析团队,持续添加、更新、优化生信云上的分析工具,集成多种生信分析流程,一键批量生成主流科研图,帮您节省时间,有更多的时间探究生物学意义。 |
----------微科盟精彩文章----------
科研 | Hort. Res.:bHLH3的异常表达会破坏类黄酮稳态网络,从而导致桑葚果实色素成分差异(国人佳作)
科研 | Aging:基于代谢组学和网络药理学的衰老及便秘相关性研究(国人佳作)
如果需要原文pdf,请扫描文末二维码,加助理获取
代谢组科研学术群期待与您交流更多代谢组学科研问题
(联系代谢组学老师即可申请入群)。
了解更多代谢组学知识,请关注下列二维码。