英文题目: associations between host gene expression, the mucosal microbiome, and clinical outcome in the pelvic pouch of patients with inflammatory bowel disease
期刊名:genome biology 影响因子:11.313
研究背景:
囊袋炎是溃疡性结肠炎(uc)回肠袋肛管吻合手术(ipaa)后常见的并发症,与炎症性肠病(ibd)相似,宿主遗传学和微生物也参与到其发病机制中。本研究对患有溃疡性结肠炎(uc)和家族性腺瘤性息肉病(fap)的ipaa患者配对活检样本进行宿主转录组和16s rrna扩增子测序数据分析。
研究流程:
研究结果:
宿主上皮组织基因表达、肠组织菌群结构、队列特征和临床表型的多因素联合分析模型
为了更好地了解ipaa手术(ipaa是一种能有效治疗溃疡性结肠炎的手术方式。它通过切除全部大肠以及部分肛管,再用回肠末端做成一个储袋-当成新的直肠用,并与保留下来的肛门括约肌进行吻合,该术式在移除病灶的同时,又能保留正常的排便功能)后宿主和微生物之间的关系,我们在一个大型元数据丰富的横断面队列中,利用微阵列和16s rrna扩增子测序分别检测了宿主基因表达和菌群结构,队列包括265名患者(51%是女性),年龄在18-78岁之间(表1),手术治疗的uc或fap患者被纳入其中,所有患者在活检样本收集前至少1年进行ipaa手术。根据临床结果患者分为家族性腺瘤性息肉(fap)、无结肠袋炎(np)、急性结肠袋炎(ap)、慢性结肠袋炎(cp)或克罗恩病类炎症(cdl)。大多数患者在囊袋(pouch,p)和前袋回肠(ppi)处进行了活检,共包括196个ppi 样本和59个pouch样本。
表1 ipaa队列的人口统计学和临床特征
组织间宿主基因表达变异度高而菌群变异度低
我们观察到囊袋(pouch,p)和前袋回肠(ppi)之间的菌群结构类似,但是它们之间的转录水平差异很大。
队列多组学整合的降维方法
为了提高菌群组成与宿主转录活性关联的效率,我们降低了宿主和微生物特征的维数。我们首先计算出,如果微生物丰度和基因表达数据之间的协方差是0.5,使用bonferroni校正最多执行10 4配对检验,并保留90%的power和α等于0.05,因此非常有必要将19908个宿主转录本和6999个观察到的otus减少到10 4次测试,大约是100个转录本和100个微生物clades。为了实现这一目标,首先我们对otus分析仅限于那些存在于多个个体并且丰富的otus(平均丰度>0.005);其次在下游分析之前,我们采用了进一步的无监督和监督策略来进行数据降维,包括多元线性建模(旨在将微生物与宿主相关联)和线性判别分析(确定哪些微生物和转录物能用于临床结果判别)(图1)。
通过这些数据降维方法,我们将19908个宿主转录物和6999个观察到的otus转化为138个特征:9个转录本主成分和9个微生物主成分分支,此外还有75个基因medoids和45个微生物分支。
图1 数据分析流程图。
组织位置和抗生素使用分别导致宿主基因表达和微生物组成最大变化
为了使gpcs、cpcs、感兴趣的基因中心medoids、炎症、抗生素使用和临床结果之间关系的可视化,我们生成了一个双标图(图2),结果发现抗生素使用与慢性口腔炎(chronic pouchitis)表型和肠球菌(enterococcus)丰度高度相关。gpc8与抗生素使用呈负相关。克罗恩病类炎症(crohn’s disease-like)与肠杆菌科细菌(enterobacteriaceae)增多有关。gpc9高表达与丰富的sutterella和有益梭菌(clostridia)相关。gpc1、gpc9、gpc6与fap和无结肠袋炎(np)最密切相关。
图2 微生物分支、宿主转录本和临床数据的双标图。紫色代表微生物,棕色和蓝色箭头分别代表宿主转录本和临床数据(包括抗生素使用和临床结果,例如fap、ap、np),样品根据炎症进行颜色区分,其范围从无(绿色)到高(红色),这个图示仅基于ppi样本。
接着我们量化了被样本位置(囊袋和前袋回肠)、临床结果、抗生素使用和炎症影响的微生物和总宿主转录本的比例(表2),结果发现宿主转录与样本位置密切相关,其次是炎症,而与抗生素使用几乎没有关联,与此相反,微生物的差异与抗生素密切相关,但与炎症或样本位置关系不大。
表2 炎症、抗生素、临床结果和样本位置对转录组和微生物的影响
为了进一步研究抗生素、样本位置、临床结果和炎症对特定微生物分支的影响,以及观察这些受影响的微生物分支系统发育关系,我们进行了独立的单变量分析(lefse)(图3)。与抗生素使用有关的bacteroides, firmicutes, tenericutes的丰度广泛减少,与抗生素使用有关的bacilli,gammaproteo的丰度增加。埃希氏菌属(escherichia)与炎症呈正相关,放线菌(actinobacteria)与炎症呈负相关。考虑使用抗生素后,sutterella属和bacteroidetes属与fap临床结果也强烈相关。actinomycetales,flavobacteria与ppi弱相关。
图3单变量分析中的微生物和临床数据之间的关系。队列中的主要临床数据是抗生素使用,炎症,组织(pouch或ppi),和临床结果(ap,np,cp,fap或cdl)。在抗生素使用分层之后,对炎症(环2)、组织类型(环3)和临床结果(环4, 5和6)的影响lda值分别计算。环的颜色强度对应于lda效应显著的分类水平,从门(最小强度)到属(最大强度)。
寄主基因表达不是囊袋(pouch)菌群组成的主要决定因素
在数据降维之后,为了检测基因-微生物的关联,使用多元线性模型来控制炎症、组织位置和抗生素使用的影响。监督(curated gene)和无监督(gpc/cpc)列表分别独立地通过maasl运行,结果只有无监督列表的结果是显著的(图4)。与cpcs显著相关的gpcs只有gpc8 和gpc9(图4)。gpc9中补体级联(cfi、c2和cfb)、干扰素调节因子1、干扰素诱导鸟苷酸结合蛋白和白细胞趋化因子ccl2的表达减少,预示gpc9高表达可能是对应于较低的炎症状态,事实上,当样品按临床结果分层时,gpc9在克罗恩病炎症患者中表达最低,在fap患者中表达最高。gpc8中脂多糖激活p38map激酶map2k6和pla2g10的表达减少,这两个元件涉及钙和脂肪介导的炎症信号,因此gpc8也可能与炎症有关。
只有四个微生物cpcs与gpc8,gpc9相关:分别是cpc1, cpc3, cpc6, cpc8。cpc1显示了明显与抗生素使用相关的几个特征:肠杆菌科细菌数量增加,拟杆菌属和厚壁菌属广泛减少,肠球菌数量最多(图4)。cpc3具有最低水平的双歧杆菌。cpc1和cpc3与gpc8和gpc9呈负相关, 这些模式表明这组具有抗生素特征的微生物与潜在的炎症基因高表达相关。
总之,宿主转录物和微生物之间的线性关系一般是适度的,占总变异的大约25%,因为其变异主要是由组织位置和抗生素使用分别驱动的。
图4 多元线性模型结果。主成分分析将数据减少到9个gpcs(左侧)和9个cpcs(中间),橙色和蓝色分别表示表达的增加或减少。右侧是同时控制炎症、组织位置和抗生素使用影响后cpcs 和gpcs之间的关联分析。
用宿主-微生物混合模型分鉴定pouch结果
接下来想探讨宿主转录本、微生物或某些组合是否可用于临床结果诊断。因为抗生素使用在临床结果中高度不均衡,并且对慢性口腔炎的预后具有高度预测性,因此我们将分析限制在没有抗生素使用的样本中。慢性结肠袋炎(cp)和克罗恩病类炎症(cdl)可以通过这个模型有效区分,特别是fap(图5)。然而交叉验证的准确性较低(auc 0.57),主要是由于模型对ap和np结果的低区分度造成的。
图5针对临床结果的线性鉴别分析。使用线性鉴别分析来确定在控制抗生素使用后,哪些基因和微生物在临床结果上最具鉴别力。在分析之前除去抗生素使用的所有样品,并使用lda混合模型进行一次性交叉验证。(a,b)使用ld1和ld2鉴别临床结果。
研究结果: