微生物分析用图的看图说明1.docx

1、微生物分析用图的看图说明11. 稀释性曲线（Rarefaction Curve）采用对测序序列进行随机抽样的方法，以抽到的序列数与它们所能代表OTU的数目构建曲线，即稀释性曲线。当曲线趋于平坦时，说明测序数据量合理，更多的数据量对发现新OTU的边际贡献很小；反之则表明继续测序还可能产生较多新的OTU。横轴：从某个样品中随机抽取的测序条数；Label 0.03 表示该分析是基于OTU 序列差异水平在0.03，即相似度为97% 的水平上进行运算的，客户可以选取其他不同的相似度水平。纵轴：基于该测序条数能构建的OTU数量。曲线解读：图1中每条曲线代表一个样品，用不同颜色标记； 随测序深度增加，被发现

2、OTU的数量增加。当曲线趋于平缓时表示此时的测序数据量较为合理。2. Shannon-Wiener 曲线反映样品中微生物多样性的指数，利用各样品的测序量在不同测序深度时的微生物多样性指数构建曲线，以此反映各样本在不同测序数量时的微生物多样性。当曲线趋向平坦时，说明测序数据量足够大，可以反映样品中绝大多数的微生物物种信息。横轴：从某个样品中随机抽取的测序条数。纵轴：Shannon-Wiener指数，用来估算群落多样性的高低。Shannon 指数计算公式：其中，Sobs=实际测量出的OTU数目；ni=含有i 条序列的OTU数目；N=所有的序列数。曲线解读：图2每条曲线代表一个样品，用不同颜色标记，

3、末端数字为实际测序条数；起初曲线直线上升，是由于测序条数远不足覆盖样品导致； 数值升高直至平滑说明测序条数足以覆盖样品中的大部分微生物。3.Rank-Abundance 曲线用于同时解释样品多样性的两个方面，即样品所含物种的丰富程度和均匀程度。物种的丰富程度由曲线在横轴上的长度来反映，曲线越宽，表示物种的组成越丰富；物种组成的均匀程度由曲线的形状来反映，曲线越平坦，表示物种组成的均匀程度越高。横轴：OTU 相对丰度含量等级降序排列。纵轴：相对丰度比例。曲线解读：图3与图4中每条曲线对应一个样本（参考右上角图标）；图3与图4中横坐标表示的是OTU（物种）丰度排列顺序，纵坐标对应的是OTU（物种）

4、所占相对丰度比例（图3为相对百分比例，图4为换算后Log值），曲线趋于水平则表示样品中各物种所占比例相似；曲线整体斜率越大则表示样品中各物种所占比例差异较大。4. 样本群落组成分析：多样本柱状图/ 单样本饼状图根据分类学分析结果，可以得知一个或多个样品在各分类水平上的物种组成比例情况，反映样品在不同分类学水平上的群落结构。柱状图（图5）横轴：各样品的编号。纵轴：相对丰度比例。图标解读：颜色对应此分类学水平下各物种名称，不同色块宽度表示不同物种相对丰度比例；可以在不同分类学水平下作图分析。饼状图（图6）在某一分类学水平上，不同菌群所占的相对丰度比例。不同颜色代表不同的物种。5. 样品OTU 分布

5、Venn 图用于统计多个样品中共有或独有的OTU数目，可以比较直观地表现各环境样品之间的OTU组成相似程度。不同样品用不同颜色标记，各个数字代表了某个样品独有或几种样品共有的OTU数量，对应的OTU编号会以EXCEL表的形式在结题报告中呈现。分析要求单张分析图，样本分组至少两个，最多5个。默认设置为97%相似度水平下以OTU为单位进行分析作图。6.Heatmap图用颜色变化来反映二维矩阵或表格中的数据信息，它可以直观地将数据值的大小以定义的颜色深浅表示出来。将高丰度和低丰度的物种分块聚集，通过颜色梯度及相似程度来反映多个样品在各分类水平上群落组成的相似性和差异性。相对丰度比例：热图（图8）中每

6、小格代表其所在样品中某个OTU的相对丰度。以图8为例，红框高亮的小格所对应的信息为：样本（R11-1Z）中OTU（OTU128）的相对丰度比例大概为0.2%。丰度比例计算公式（BrayCurtis算法）：其中，SA,i=表示A样品中第i个OTU所含的序列数SB,i=表示B样品中第i个OTU所含的序列数样品间聚类关系树：进化树表示在选用成图数据中，样本与样本间序列的进化关系（差异关系）。处于同一分支内的样品序列进化关系相近。物种/OTU丰度相似性树：丰度相似性树表示选用成图的数据中样品与样品中的OTU或序列在丰度上的相似程度。丰度最相近的会分配到同一分支上。客户自定义分组：根据研究需求对菌群物种

7、/OTU研究样本进行二级分组二级物种/OTU分组：将下级分类学水平物种或OTU分配到对应的上级分类学水平，以不同颜色区分；二级样品分组：根据研究需要，对样品进行人为的分组，以不同颜色区分。7.主成分分析PCA(PrincipalComponentAnalysis)在多元统计分析中，主成分分析是一种简化数据集的技术。主成分分析经常用于减少数据集的维数，同时保持数据集中对方差贡献最大的特征，从而有效地找出数据中最“主要”的元素和结构，去除噪音和冗余，将原有的复杂数据降维，揭示隐藏在复杂数据背后的简单结构。通过分析不同样品的OTU组成可以反映样品间的差异和距离，PCA运用方差分解，将多组数据的差异反

8、映在二维坐标图上，坐标轴为能够最大程度反映方差的两个特征值。如样品组成越相似，反映在PCA图中的距离越近。横轴和纵轴：以百分数的形式体现主成分主要影响程度。以图9为例，主成分1（PC1）和主成分2（PC2）是造成四组样品（红色，蓝色，黄色和绿色）的两个最大差异特征，贡献率分别为41.1%和27.1%。十字交叉线：在图9中作为0点基线存在，起到辅助分析的作用，本身没有意义。图例解读：PCA分析图是基于每个样品中所含有的全部OTU完成的；图9中每个点代表了一个样本；颜色则代表不同的样品分组；两点之间在横、纵坐标上的距离，代表了样品受主成分（PC1或PC2）影响下的相似性距离；样本数量越多，该分析意

9、义越大；反之样本数量过少，会产生个体差异，导致PCA分析成图后形成较大距离的分开，建议多组样品时，每组不少于5个，不分组时样品不少于10个；图10中的圆圈为聚类分析结果，圆圈内的样品，其相似距离比较接近。8. RDA/ CCA 分析图基于对应分析发展的一种排序方法，将对应分析与多元回归分析相结合，每一步计算均与环境因子进行回归，又称多元直接梯度分析。主要用来反映菌群与环境因子之间的关系。RDA是基于线性模型，CCA是基于单峰模型。分析可以检测环境因子、样品、菌群三者之间的关系或者两两之间的关系。横轴和纵轴：RDA和CCA分析，模型不同，横纵坐标上的刻度为每个样品或者物种在与环境因子进行回归分析

10、计算时产生的值，可以绘制于二维图形中。图例解读：冗余分析可以基于所有样品的OTU作图，也可以基于样品中优势物种作图；箭头射线：图11中的箭头分别代表不同的环境因子（即图中的碳酸氢根离子HCO3-，醋酸根离子AC-等，图中的其它环境因子因研究不同代表的意义不同，因此不再赘述）；夹角：环境因子之间的夹角为锐角时表示两个环境因子之间呈正相关关系，钝角时呈负相关关系。环境因子的射线越长，说明该影响因子的影响程度越大；图11中不同颜色的点表示不同组别的样品或者同一组别不同时期的样品，图中的拉丁文代表物种名称，可以将关注的优势物种也纳入图中；环境因子数量要少于样本数量，同时在分析时，需要提供环境因子的数据

11、，比如pH值，测定的温度值等。9.单样品/多样品分类学系统组成树根据NCBI提供的已有微生物物种的分类学信息数据库，将测序得到的物种丰度信息回归至数据库的分类学系统关系树中，从整个分类系统上全面了解样品中所有微生物的进化关系和丰度差异。单样品图（图12）：可以了解单样品中的序列在各个分类学水平上的分布情况。图例解读：图12中不同的层次反映不同的分类学水平；分支处的圆面积说明了分布在该分类学水平，且无法继续往下级水平比对的序列数量，面积越大，说明此类序列越多；每个分支上的名词后面的两组数字分别表示比对到该分支上的序列数和驻留在该节点上的序列数；图13中为某单一水平物种分布情况，并非是序列分布。多

12、样品图（图14）：比对多个样品在不同分类学分支上序列数量差异。图例解读：比对不同样品在某分支上的序列数量差异，通过带颜色的饼状图呈现，饼状图的面积越大，说明在分支处的序列数量越多，不同的颜色代表不同的样品。某颜色的扇形面积越大，说明在该分支上，其对应样品的序列数比其他样品多。多样品在做该分析时，建议样品数量控制在10个以内，或者将重复样本数据合并成一个样本后，总样品数在10个以内。10.系统发生进化树在分子进化研究中，基于系统发生的推断来揭示某一分类水平上序列间碱基的差异，进而构建进化树。图例解读：图15中体现的是序列进化差异情况，处在同一分支上的物种说明进化关系较近。图15左下角的图例为距离

13、标尺，分支距离越长，进化关系越远。11.(un)WeightedUniFracPCoA/Tree分析利用各样品序列间的进化信息来计算样品间距离，反映环境样品在进化树中是否有显著的微生物群落差异。PCoA（principalco-ordinatesanalysis）是一种研究数据相似性或差异性的可视化方法，通过一系列的特征值和特征向量进行排序后，选择主要排在前几位的特征值，PCoA可以找到距离矩阵中最主要的坐标，结果是数据矩阵的一个旋转，它没有改变样品点之间的相互位置关系，只是改变了坐标系统。通过PCoA可以观察个体或群体间的差异。图例解读：图16和图17中不同颜色代表不同分组；PCoA分析建议

14、不分组时，样本数量不少于10个；多组样本时，每组样本数量不少于5个；对于某一功能基因，进行进化树分析时，建议采用OTU数目控制在10,000以内，或者由客户指定分析优势OTU个数。12.NMDS分析NMDS（NonmetricMultidimensionalScaling）常用于比对样本组之间的差异，可以基于进化关系或数量距离矩阵。横轴和纵轴：表示基于进化或者数量距离矩阵的数值在二维表中成图。图例解读：图18中不同的颜色代表不同的分组；建议不分组时，样本数量不少于10个；多组样本时，每组样本数量不少于5个；图18中的点代表样本，点与点之间的距离表示差异程度。13.含相似性树柱状图根据样品中相似

15、程度进行排布，并绘制对应样本树状图反映样本中群落结构。图例解读：图19中左侧是相似度树状图，样本之间的差异越小，样本便会处在相近的同一分支上；右侧柱状图，展示样本中微生物的群落结构。不同颜色代表不同物种。14.Unifrac显著性差异分析比较样品间进化差异的显著性分析。图例解读：图20横坐标为两组样品；纵轴坐标为unifrac进化距离（序列差异）。15.单因素unweightedunifracPCA分析在某单一因素上，进行unweightedunifracPCA分析。图例解读：图21横轴为不同变量（本例为不同年龄阶段）下的样本；纵坐标为主成分，图21中显示同一年龄阶段内和不同年龄阶段间的由主成

16、分导致的差异情况。16. 个性化分析案例展示案例描述一Community composition of root-associated fungi in a Quercus-dominated temperate forest:codominance of mycorrhizal and root-endophytic fungi. Ecol Evol. 2013 May; 3(5):1281-93.样本来源样本数量高变区测序平台测序量（reads）植物根系/泥土159/38ITSRocheGS70,495样本来自于以日本橡木为主的温带森林，使用454测序分析多个生态系统的真菌多样性。样品来源

17、为植物根系样本，从12株植物中提取了159份根系。分析结果表明，外生真菌群落和根系内生真菌群落互相作用，维持了这种以日本橡木为主的温带森林的生态环境。图例解读： 图22构建了地下植物与真菌的相互关系网络；图22中灰色圆点表示与植物共生的微生物物种，菱形代表真菌OTU，它们之间的关系用灰色连接线表示；它们之间的密集程度越高表示它们之间相互作用被观察的次数越多；互利共生真菌OTU用粉色菱形表示；寄生微生物OTU用橘色菱形表示；未知功能OTU用蓝色菱形表示；A：子囊真菌门；B：担子菌门；G：球囊菌门；U：表示门水平未知真菌。案例描述二Theignoreddiversity:complexbacter

18、ialcommunitiesinintensivecareunitsrevealedby16Spyrosequencing.SciRep.2013;3:1413.室内微生物群落对日常生活中人类健康起着重要作用，尤其是医院的重症监护室。采用扩增焦磷酸测序研究ICU中微生物群落可以检测多种微生物序列，与现有的传统标准培养技术相比，有极大的优越性。传统培养方式只能检测总细菌多样性的2.5。结合外部环境与物种系统发育谱分析发现，许多微生物与潜在的人类病原菌相关，当然也包含有益菌，一共7个门76个属。此外，丙酸杆菌属，假单孢菌属和伯克霍尔德氏菌被确定为感染的重要来源。在地板，医疗器械和工作间微生物组成有

19、显著差异，但网络分析和一致性分子指纹印记分析发现该三个地点微生物组成也有一定的相似性。这些信息将帮助加护病房进行新的公共卫生风险评估，帮助建立新的卫生协议，帮助深入了解医院获得性感染的情况。样本来源样本数量高变区测序平台测序量（reads）重症监护室（ICU）中地板，医疗设备及工作间表面样品3416SrRNARoche454FLX+5000/样本图23为丰度最高的40种菌(OTU)在不同采样区域的分布情况。不同颜色代表不同样本来源，绿色为ICU地板，红色为医疗设备，蓝色为工作间。结点面积表示该OTU在三个采样区域之间的相对丰度。案例描述三SkinMicrobiomeImbalanceinPat

20、ientswithSTAT1/STAT3DefectsImpairsInnateHostDefenseResponses.JournalofInnateImmunity.2013,DOI:10.1159/000351912.CMC和HIES均属于罕见的原发性免疫缺陷病。本文研究了CMC和HIES病人皮肤及口腔微生物与正常人群的微生物菌群差异。结果表明，相比正常人体的微生物组成（主要包括棒杆菌属和放线菌科等），病人的所携带的微生物中革兰氏阴性菌比例增加。而革兰氏阴性菌（如不动杆菌属）可以抑制机体对假丝酵母和金色葡萄球菌的免疫反应，结果很可能导致患者对这些菌感染的敏感性增加。本文的研究表明，此类免

21、疫缺陷病人所携带的微生物菌群可以影响宿主的免疫防御系统，进而有可能通过基于微生物的辅助疗法来治疗免疫缺陷的患者。样本来源样本数量高变区测序平台测序量（reads）慢性皮肤黏膜念珠菌病（CMC）、高IgE综合症（HIES）患者皮肤及口腔微生物6016SrRNAV4区MiSeq6389/样本图24为免疫缺陷疾病各亚型患者皮肤和口腔微生物群落结构的协方差。（a）图表示35个皮肤样本的微生物多样性的差异；（b）图表示21个口腔样本微生物多样性差异。其中，绿色方块表示健康对照组，红色圆圈表示CMC，蓝色三角形表示HIES。健康人体皮肤表面的微生物主要包括葡萄球菌和棒状杆菌，而疾病组多为莫拉菌科。培训文档

22、境微生物多样性分析-基础境微生物多样性分析-高级常见问题Q1.能对哪些环境下的样品进行分析？A 目前已经成功运用Roche454技术发表的文章涵盖农业、土壤、林业、海洋、矿井（石油等）、人体医学等诸多领域，共计近1400篇，Paper目录可向本公司索取。其中有大量文献发表于国际顶级杂志上，包括Science、Nature、ISME、PNAS、AEM等。Q2.DGGE技术与Roche454高通量测序技术在环境微生物群落多样性研究中有什么区别？A DGGE等分子指纹图谱技术，在其实验结果中往往只含有数十条条带，只能反映出样品中的优势种群，也无法得到细菌种类及其绝对含量。而Roche454高通量测序

23、技术能同时对样品中的优势种群及微量菌进行检测，获得样品中的微生物群落组成，并将其含量进行数字化。Q3.在医学领域有哪些应用？A 在人类的皮肤、口腔、呼吸道、胃肠道和尿道等处存在着大量与人体健康密切相关的正常菌群。它们能够合成并辅助机体吸收一些必需氨基酸和维生素；加工诸如植物多糖等人类饮食中一些难以消化的组分；占据人体的不同粘膜表面并产生天然的抗生素，抑制有害菌的着落与生长。目前，人们已经采用传统方法对人人体微生物作了许多研究，但是这种方法需要对微生物进行纯培养，从而大大限制了我们对机体正常菌群，尤其是许多不可培养微生物的认识。第二代高通量测序技术能对人体菌群进行深度测序，从而精确检测到千分之一

24、机体微生物群落的数量和组成结构的变化，突破了传统方法基于纯培养的限制，能使我们从整体上认识微生物群落与宿主之间的相互关系及其对人体健康的影响。例如，利用第二代高通量测序技术可以对处于疾病状态下的人体微生境进行研究，比较分析正常和疾病状态下或疾病不同进程人体微生物群落的结构和功能变化；可以对正常人群与某些疾病患者体内的微生物群体多样性进行比较分析，研究微生物群落变化同疾病之间的关系；通过深度测序还可以快速地发现和检测常见病原及新发传染病病原微生物。Q4.高通量测序是否需要做平行样？A 随着高通量测序的发展，提出研究过程中设置重复样的要求，也日益被大家所接受。在高通量研究中，设置重复样不仅仅是体现

25、了科研的一种严谨态度，同时也体现了结果的真实性，避免了个体差异造成的影响。目前的发展趋势：普通样品（水体，土壤，处理样本内部设置重复等）一般会要求设置重复样；对于大面积研究水体和土壤等样本时，可以采用多点采样混样研究。稀有样本（稀有动物粪便，冰川极地冰样，过多采样会影响整体环境构成的样本等）对于重复样设置的要求相对比较宽松。Q5.454一块PTP板最多可以放多少样本？1个相同的样本分别研究细菌，真菌，古菌是否可以放在同一块PTP板上测序？A理论上，知道每个样本的测序量和一个454PTP板能够容纳的测序量，便能得知可以测序的样本量。实际上，考虑到测序质量的因素，不可如此计算上样量。一块PTP能够

26、容纳多少样本往往应该考虑以下两个因素：1. barcode数量：barcode是用于区分上样样品的标签序列，标签的个数限制着一块PTP能够容纳多少样本；2. PTP板分区情况：PTP板分区越多可以增加上样的样品个数，但分区越多，会影响测序总量。一块PTP板，在不分区的情况下，一般可以产出80万条序列，综合考虑barcode数量，分区情况以及单样本的测序量，便可计算一块PTP板容纳的样本数量。对于同一个样本，分别研究细菌，真菌，古菌一般是可以放到同一个板里进行测序，但要考虑到不同类型PCR产物的长度对测序结果的影响，如果PCR产物长度差异较大，则不可以放入同一块PTP板进行测序。参考文献Alek

27、sandarD.Kostic,DirkGevers,etal.GenomicanalysisidentifiesassociationofFusobacteriumwithcolorectalcarcinoma.GenomeResearch,2012,22:292298.ChiachiHwang,FangqiongLing,etal.Microbialcommunitydynamicsofanurbandrinkingwaterdistributionsystemsubjectedtophasesofchloraminationandchlorinationtreatments.Microbi

28、ol,2012,78(22):7856.CraigE.Nelson,StuartJ.Goldberg,etal.Coralandmacroalgalexudatesvaryinneutralsugarcompositionanddifferentiallyenrichreefbacterioplanktonlineages.TheISMEJournal,2013,7:962979.ElieJami,AdiIsrael,etal.Exploringthebovinerumenbacterialcommunityfrombirthtoadulthood.TheISMEJournal,2013.2:

29、111.ErikCvonRosenvinge,YangSong,etal.Immunestatus,antibioticmedicationandpHareassociatedwithchangesinthestomachfluidmicrobiota.TheISMEJournal,2013.33:113.KrissiM.Hewitt,CharlesP.Gerba,etal.OfficespacebacterialabundanceanddiversityinthreemetropolitanAreas.PLoSONE,2012,7(5):e37849.SizhongYang,XiWen,et

30、al.PyrosequencinginvestigationintothebacterialcommunityinpermafrostsoilsalongtheChina-Russiacrudeoilpipeline(CRCOP).PLoSONE,2012,7(12):e52730.VasviChaudhry,AteequrRehman,etal.Changesinbacterialcommunitystructureofagriculturallandduetolong-termorganicandchemicalamendments.MicrobEcol,2012,64:450-460.Weig