科学网—246.三元图在微生物组研究中的应用与绘图实战

在微生物组领域,我们主要利用的是基于相位网格的三元图,其特点是每一种成分在相应等边三角形的顶点的比例为100%,在其对面的线为0%,将零点线和顶点按比例划分,用以估计各成分的含量。如图1:y对面的底线代表落在底线上的所有点在y中占比为0,随着平行线逐渐靠近顶点,落在相应平行线上的点在y中的占比越来越高。

在解读三元图时,我们可以通过点的位置快速获得其在三个分组中的相对比例信息:目标点越靠近一角的顶点,说明他在相应的分组中比例越高;反之其相对比例就越低。根据图1,经过点的平行线在等边三角形两边的截距代表该点在对应顶点分组的占比,因此可以得出该点的在三个分组中的占比情况。如图2:过绿色点分别向三条边做平行线,然后按逆时针方向依次读取平行线在三条边的截距约为(0.33,0.33,0.33),因此该点在x,y,z三个组分中的占比为0.33,0.33,0.33;依此原则可得出,红点占比为0.1,0.2,0.7;蓝点占比为0.2,0.7,0.1;该结果跟我们的作图代码是相符的。

在微生物多样性实际分析中,由于OUTs数目繁多,我们并不需要获得如此精确的占比情况。一般来说,三元图不同的点代表不同的OTUs(或其他分类水平),点的大小代表该OTUs的平均丰度(一般需要log2或log10等转换)。不仅如此,还可以对OTUs在各微环境中的丰度数据进行统计检验后,得出各OTU分别在哪种微环境中显著富集,此时三元图不仅表现出OTUs或者物种在微环境中的相对比例,还包含显著性统计结果,它打破了火山图或韦恩图两两比较的结果,总共展示了6次两两比较的结果,即每个组的富集情况是相对于其他两组的。

三元图的应用有非常大的局限性,就是必须有3个实验组。只有两组的课题就不要考虑此类图型的应用了。值得一提的是,实验设计3组也是有着最高性价比和最广泛的应用,因为2组只有1种比较组合,而3组仅成本上升50%,却有3种比较组合,结果的图表量瞬间上升为3倍,性价比极高。而3组以上组合太多,反而不利于作者把故事讲清楚,读者阅读和理解也存在一定困难。因此强烈推荐在尽量可能的情况下,设计为3组起的实验,可以有更多的分析和讨论空间,也可以应用有监督的排序分析、三元图等常用结果样式。下面分别选取来自FrontierinMicrobiology、PNAS和NatureCommunications杂志中近年的三元图结果进行解读和点评。

图片描述:三元图展示不同区域(土体土、根际土、根内)样本的群落结构在门水平上的差异情况,每个圆点代表一个OTU,OTU的大小、颜色和位置分别代表其相对丰度、菌门水平名称和分组情况。

DistributionofOTUsandphylaacrossdifferentcompartments.(A)TernaryplotofallOTUsplottedbasedonthecompartment(Bulksoil,Rhizospheresoil,Endosphere)specificity.EachcirclerepresentsoneOTU.Thesize,colorandpositionofeachOTUrepresentsitrelativeabundance,bacterialphylaandaffiliationoftheOTUwithdifferentcompartments,respectively.

图注描述注意事项:

文中的图片解读:

主要描述了比较突出的情况,如丰度较高、规律性变化或明显聚集在某区域的门。在门水平的不同部位中,细菌群落结构明显不同。这些差异的主要原因是由于厚壁菌门在根内相对丰度较高引起,而在土体土和根际土中它们的丰度却很低。变形杆菌和拟杆菌的相对丰度从土体土到根际土再到根内逐渐增加,同时候选门AD3、芽单胞菌门和绿弯菌门的相对减少,它们共同构成根内微生物群落丰度<4%。

Bacterialcommunitystructureswereclearlydifferentinthedifferentcompartmentsatthephylumlevel.ThesedifferencesweremainlydrivenbystrongrelativeenrichmentofFirmicutesintheendosphere-derivedsequencedatasets,comparedtotheirverylowabundancesinthebulkandrhizospheresoils.TherelativeabundancesofProteobacteriaandBacteroidetesincreasedprogressivelyfrombulktorhizospheresoiltotheendosphere,withaconcomitantdecreaseinthoseofcandidatedivisionAD3,GemmatimonadetesandChloroflexi,whichcollectivelyconstituted<4%ofendospherecommunities.

本例选自2016年发表于PNAS杂志上一文的图4(Zgadzajetal.,2016),介绍了豆科模式植物百脉根的结瘤突变体中根际微生物组变化的研究。

图片描述:

TernaryplotsdepictingcompartmentRAofallOTUs(>5‰)forWTSampleID(A;WT;n=73)andmutantSampleID(B;nfr5-2,nfr5-3,nin-2,andlhk1-1;n=118)acrossthreesoilbatches(CAS8–CAS10).EachpointcorrespondstoanOTU.ItspositionrepresentsitsRAwithrespecttoeachcompartment,anditssizerepresentstheaverageacrossallthreecompartments.ColoredcirclesrepresentOTUsenrichedinonecompartmentcomparedwiththeothers(greeninroot,orangeinrhizosphere,andbrowninrootSampleID).AggregatedRAsofeachgroupofenrichedOTUs(root-,rhizosphere-andsoil-enrichedOTUs)ineachcompartmentfortheWTSampleID(C;WT;n=73)andmutantSampleID(D;nfr5-2,nfr5-3,nin-2,lhk1-1;n=118)areshown.Ineachcompartment,thedifferencefrom100%RAisexplainedbyOTUsthatarenotsignificantlyenrichedinaspecificcompartment.

总结:

该图在例1的基础上添加了组间显著性差异比较的结果,包含了6次两两比较和三次韦恩图比较的结果,信息高度概括。同时,作者还进一步结合相对丰度的箱线图突出组间差异;此外,通过使用相同的色系进行A与B图,C与D图的横向比较,突出明显的分布差异,使用两类材料微生物组不同的规律一目了然。

图5.土壤、叶片和毛虫共享的微生物。

a.在土壤,食用离体叶片,食用植株的毛虫中共有OTU的分布。b.在植物,食用离体叶片,食用植株的毛虫中共有OTU的分布。c.食用离体叶片毛虫和植株毛虫肠道共有和特有OTU的数量。图中颜色含义和a,b图中的一致。

结果

两种喂养方式毛虫共有的核心微生物群落也存在于植株根系(19.1%)和叶片(20.3%)中(图2a–c),同时也有特有的微生物群,大约16.7%微生物仅仅发现在毛虫中。毛虫核心微生物群组成细菌主要有变形菌门、放线菌门和厚壁菌门,真菌为未注释的OTU。值得注意的是以完整植株为食的毛虫大部分OTU也存在与土壤中(75%)。以完整植株为食的毛虫中有超过离体叶片为食的毛虫三倍的OTU数量。

Caterpillarsfedonintactplantsanddetachedleavessharedacommoncoremicrobiomewhichwasalsopresentintheleaves(20.3%oftheirmicrobiome)andintheroots(19.1%)(Fig.2a–c),butalsoharboreduniquemicrobes;16.7%ofthecaterpillarmicrobiomewasfoundonlyincaterpillars.ThiscoremicrobiomeofcaterpillarsconsistedpredominantlyofProteobacteria,Acidobacteria,Firmicutes,andunclassifiedfungi(SupplementaryFigs6,7).Remarkably,forcaterpillarsfedonintactplants,alargeproportionoftheOTUsfoundincaterpillars,wasalsodetectedinthesoil(75%;representedasnumbers1and4inFig.2a).MicrobiomesofcaterpillarsfeddetachedleaveshadvirtuallynoadditionalOTUsthatwerenotalsofoundincaterpillarskeptonintactplants(Fig.2c),butthemicrobiomesofthelattercontainedthreetimesmoreOTUs.

总结

由于前期数据处理是三元图主要的难点,所以在这里将数据处理和可视化分开,使用时便于检查异常、调整分析细节。

按需求安装,没必要每次都运行该安装代码,一般运行一次即可。

函数data_clean参数介绍:

三元图展示不同处理组(KO,OE,WT)的样本的群落结构在OTUs上的丰度差异情况,每个圆点代表一个OTU,OTU的大小和位置分别代表其相对丰度和分组情况。

由于该图中并未全部展示所有OTUs,所以在描述的时候,最好对数据筛选的参数进行描述:

从图中可以看到KO(基因敲除,knock-out)组与OE(过表达,over-expression)和WT(野生型,wild-type)组存在丰度差异,即基因的有无可对微生物群落的丰富度引起变化。

上图能够展示的信息和按主丰度着色后的可视化方案是一致的,但叙述时可以像图3一样着重介绍突出情况,比如丰度较高的OTUs有那些,富集情况很特殊的有哪些。

函数top_OTUs参数介绍:ggtern输入文件筛选丰度前N(10)的OTUs

图7.三元图展示OTUs的实验分组特异性

三元图展示不同处理组(KO,OE,WT)的样本的群落结构在OTUs上的丰度差异情况,每个圆点代表一个OTU,OTU的大小、颜色和位置分别代表其相对丰度、主丰度OTUs名称和分组情况。值得注意的是,红色原点所代表的ASV_5在KO组中丰度几乎为0,但在WT和OE均有相对较高的丰度。具体可以再结合物种注释进一步描述。

除了上面直接展示OTUs在不同分组中的丰度富集情况,还可以通过显著性富集分析,获得在不同分组中显著富集的OTUs,最后进行可视化(例2)。

图8.三元图展示各组特异富集的OTUs

除了图6中的信息外,该图还展示了在KO,OE,WT三个分组中,分别相对于另外两个分组显著性富集的OTUs。WT组一个特异的都没有。

由于该图中并未全部展示所有OTUs,所以在描述的时候,最好对数据统计分析过程的参数进行描述:

责编:刘永鑫中科院遗传发育所

版本更新历史

1.0.0,2020/6/6,吴一磊,初稿;文涛审阅;刘永鑫审阅

1.0.1,2020/6/9,吴一磊,小修,添加模拟数据,可独立运行的代码;刘永鑫审阅

1.0.2,2020/6/22,吴一磊,大修,添加实例讲解,代码封闭函数;刘永鑫审阅

1.0.3,2020/6/25,席娇,文字修改

1.0.4,2020/6/27,刘永鑫,整合修改意见;添加一篇NC实例讲解

THE END
1.「东方6+1基本走势图」东方6+1走势图彩经网基本走势图频道提供专业的基本走势图带连线,并免费提供上百种东方6+1走势图指标,为还在中奖路上的彩民服务。https://m.cjcp.cn/zst/cjwdf/df61_haoma.html
2.东方6+1第2024139期11月30日开奖号码结果查询乐彩网东方6+1 热点导读: 乐彩论坛乐彩网乐币排行榜太湖字谜3D走势图3D试机号分析3D和值和尾3D和值跨距排列3走势图双色球号码分布 关于乐彩|联系方式|意见反馈|免责声明|用户协议|隐私政策|网站地图|手机版 本站开奖数据仅供参考,最终结果请以福彩和体彩官方公告为准。本站展示的模拟开机号和模拟试机号来源于网络,与https://www.17500.cn/kj/detail-df61-2024139.html
3.浙江体彩网>>历史开奖查询>>6+1浙江体彩网、足彩竞猜、开奖分析、6+1、超级大乐透、中国足球彩票、20选5、大奖分布、进球游戏https://www.zjlottery.com/win/SResult.asp?flag=1
4.浙江6+1基本走势图2024088 4 2 7 56 50 1 22 8 17 9 9 11 5 3 2 23 17 3 10 8 5 14 36 3 1 3 7 4 8 6 0 15 3 1 7 9 36 4 6 13 2 13 34 19 1 5 10 9 4 3 0 7 10 8 3 6 14 1 5 9 9 2 4 26 15 5 3 10 14 1 2024089 5 3 8 57 51 2 23 7 18 1 10 12 6 3 3 24 https://lotto.sina.cn/trend/qxc_qlc_proxy.d.html?actionType=chzs&lottoType=zj61
5.浙江6+1走势图总和走势综合版带连线首页>浙江6+1总和走势 号码走势 历史开奖 基本走势 定位走势 五行走势 两面走势 升平降走势 总和走势 012路走势 大小比走势 单双比走势 质合比走势 期 至期 查看 近30期 近50期 近100期 期数 开奖号码 总和 单 双 大≧27 小< 27 24137期 787760 https://www.cpzyrj.com/zj61zhzs
6.浙江体彩6+1走势图标注形式选择: 折线 不带遗漏数据 遗漏分层 期号第一位第二位第三位第四位第五位第六位第七位 0123456789012345678901234567890123456789012345678901234567890123456789 24087 8 10 4 15 1 5 16 2 9 7 4 13 35 3 10 2 6 3 7 5 7 14 2 1 4 8 35 3 5 12 1 12 33 3 24 5 9 8 3 2 3 6 9 https://www.55128.cn/zs/63_356.htm?record=50
7.证券投资实务实训实验课精品资源共享课1、主菜单栏 如果对股票软件还不熟悉,那么通过菜单栏可以找到几乎所有它能够实现的功能。(图3) 2、综合信息栏和状态栏 在画面最低端有两行信息,即是综合信息栏和状态栏。 ①综合信息栏 图4 上证指数、涨跌、成交金额(亿):鼠标单击此处即可直接切换到上证指数分时走势图。 https://www.wdu.edu.cn/gljg/jjxy/xkjs/jpzygxk/201605/t20160527_25181.shtml
8.USDCNYUSDCNH配对交易4.4 最优协整货币对热力图 4.5 滑动时间窗口寻找最优协整货币对 5. CNY-CNH协整分析 5.1 整体数据回归的残差直方图 5.2 逐年回归的残差直方图 6. CNY-CNH配对策略交易代码 6.1 提取交易信号 6.2 计算回测收益并画图 7. 结语 1. 策略思路: 设 代表USDCNY序列, https://blog.csdn.net/lijunscut/article/details/126954619
9.福彩3D选号秘籍单选图里0、1、2路位置不固定,只在3D的百、十、个位填写当期开出的号码路数,如2006006期7 1 3在百位填1、十位填1、个位填0表示路数复式110组合,表示百、十、位的号码1、4、7个位号码是0、3、6、9,投资72元,中奖100元。 012路综合走势图里的组选图可以帮助我们判断哪路号码将缺失,一般情况我们多考虑http://www.360doc.com/content/15/0927/19/27902029_501865318.shtml
10.股票技术指标分析详解单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,股票技术分析详解,技术分析,分时图的基础知识,K,线图的基础知识,第一节,分时图的基础知识,分时走势图也叫即时走势图,它是把股票市场的交易信息实时地用曲线在http://m.zhuangpeitu.com/article/252979418.html
11.实践阿里QuickBI智能图表入门为了分析毛利额下滑原因,您需要先查看月度毛利额的走势,确认毛利额在哪个月份出现下滑。同时需要分析销售额的月度走势,确认近几个月的销售情况。本文以线图来展示月度毛利额和月度销售额的走势图。 查看月度毛利额统计 您可以通过线图查看毛利额的月度走势数据。 https://www.jianshu.com/p/f498b1ef7aae
12.煤炭价格最新行情走势图详解及指导手册,最新动态与操作指南1、确定数据来源 通过搜索引擎查找各大煤炭行业网站、财经网站等,关注其发布的煤炭价格信息,确保数据来源的可靠性,以确保获取的数据准确。 2、访问相关网站 打开浏览器,输入网址,访问相关网站,寻找煤炭价格信息,一般会有专门的煤炭价格行情页面。 3、查看煤炭价格走势图 https://zhuoyanxinli.com/post/8006.html
13.ICS2023国际尿控协会年会华西医院沈宏罗德毅教授团队16项研究亮相,展望结果:共收集中文社交媒体UI相关帖子24,770篇,LDA将其分为六个UI内容主题:1. UI对生活质量的影响(5,908);2. 传统中医治疗UI(4,793);3. 女性压力性尿失禁(4,651);4. 神经源性下尿路功能障碍(4,408);5. UI的病因(3,348);和6. 老年患者并发尿失禁(1,662)。与主题1、2、4和6相关的情感都是负面https://www.cd120.com/public/department/Urology/dynamics/76136.html
14.Python学习笔记StatsModels统计回归——可视化1、如何认识可视化? 图形总是比数据更加醒目、直观。解决统计回归问题,无论在分析问题的过程中,还是在结果的呈现和发表时,都需要可视化工具的帮助和支持。 需要指出的是,虽然不同绘图工具包的功能、效果会有差异,但在常用功能上相差并不是很大。与选择哪种绘图工具包相比,更重要的是针对不同的问题,需要思考选择什么https://www.flyai.com/article/893