微生物多样研究—β多样性分析

一、β-多样性分析

1. 样品间距离计算

样品间的物种丰度分布差异程度可通过统计学中的距离进行量化分析,使用统计算法Euclidean,Bray-Curtis,Unweighted_unifrac,weighted_unifrac等,计算两两样品间距离,获得距离矩阵,可用于后续进一步的beta多样性分析和可视化统计分析。

例如:将距离矩阵使用热图表示可直观观察样品间的差异高低分布。


基于Bray-Curtis距离heatmap图示意图

2. PCA 分析

主成分分析(PCA,PrincipalComponent Analysis),是一种应用方差分解,对多维数据进行降维,从而提取出数据中最主要的元素和结构的方法。  

应用PCA分析,能够提取出最大程度反映样品间差异的两个坐标轴,从而将多维数据的差异反映在二维坐标图上,进而揭示复杂数据背景下的简单规律。  

如果样品的群落组成越相似,则它们在PCA图中的距离越接近。

3.  PCoA分析

主坐标分析(PCoA,PrincipalCo-ordinates Analysis),是一种与PCA类似的降维排序方法,通过一系列的特征值和特征向量排序从多维数据中提取出最主要的元素和结构。  

可以基于bray_curtis、WeightedUnifrac距离和UnweightedUnifrac距离分别来进行PCoA分析,并选取贡献率最大的主坐标组合进行作图展示。  

如果样品距离越接近,表示物种组成结构越相似,因此群落结构相似度高的样品倾向于聚集在一起,群落差异很大的样品则会远远分开。


PCA & PCoA分析示意图

※   当PCA或PCoA分析的前两个成分(解释度)较小(如pc1与pc2之和小于50%)时,可尝试将前三个成分用于对假设因素进行验证,并作三维图来反应样品间群落组成的关系。

4.  NMDS分析

非度量多维尺度分析(NMDS分析)是一种将多维空间的研究对象(样品或变量)简化到低维空间进行定位、分析和归类,同时又保留对象间原始关系的数据分析方法。 

适用于无法获得研究对象间精确的相似性或相异性数据,仅能得到他们之间等级关系数据的情形。

基本特征是将对象间的相似性或相异性数据看成点间距离的单调函数,在保持原始数据次序关系的基础上,用新的相同次序的数据列替换原始数据进行度量型多维尺度分析。换句话说,当资料不适合直接进行变量型多维尺度分析时,对其进行变量变换,再采用变量型多维尺度分析,对原始资料而言,就称之为非度量型多维尺度分析。

特点是根据样品中包含的物种信息,以点的形式反映在多维空间上,而对不同样品间的差异程度,则是通过点与点间的距离体现的,最终获得样品的空间定位点图。


NMDS分析示意图

5. 多样品相似度树状图

利用树枝结构描述和比较多个样品间的相似性和差异关系。 

首先使用描述群落组成关系和结构的算法计算样品间的距离,即根据beta多样性距离矩阵进行层次聚类(Hierarchicalcluatering)分析,使用非加权组平均法UPGMA(Unweightedpair group method with arithmetic mean)算法构建树状结构,得到树状关系形式用于可视化分析。


多样品相似度树状图示意图

6.  PLS-DA分析

PLS-DA(PartialLeast Squares Discriminant Analysis)分析是以偏最小二乘回归模型为基础,作为一种有监督的模式识别方法,根据给定的样品分布/分组信息,对群落结构数据进行判别分析。 

PLS-DA通过寻找物种丰度矩阵和给定的样品分布/分组信息的最大协方差,从而在新的低维坐标系中对样品重新排序。

PLS-DA可以减少变量间多重共线性产生的影响,因此,比较适合用于微生物群落数据的研究。


PLD-DA分析示意图

分析时,会计算每个物种的VIP(Variableimportance in projection)系数(VIP值需>1,值越大,说明该物种对于组间差异的贡献越大)

7.组合(变换)分析图

特点: 

集多种分析结果于一身组合成图,即一整图表解释多种生物学意义。 

展现形式、分析名称发生变化并进行重新调整,但所表述的生物学意义未变化。 

具有一定的观赏性。

  分析形式多种多样,但万变不离其宗。

例如:样本聚类树与柱状图组合分析


样本聚类树与柱状图组合分析示意图
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 195,783评论 5 462
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 82,360评论 2 373
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 142,942评论 0 325
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,507评论 1 267
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 61,324评论 5 358
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,299评论 1 273
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,685评论 3 386
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,358评论 0 254
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,652评论 1 293
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,704评论 2 312
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,465评论 1 326
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,318评论 3 313
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,711评论 3 299
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 28,991评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,265评论 1 251
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,661评论 2 342
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,864评论 2 335

推荐阅读更多精彩内容