网络图 Network
网络图虽然给人高大上的感觉,但是由于信息太多,无法给读者提供读有效的可读信息或是读者不知道该理解什么,总是让人望尔却步。那是因为大家太不了解网络,自己读不懂网络想表达的意思及其重要性。
 
因此我要举一个网络分析重要的例子,大家都知道2012年诺贝尔奖得主Yamanaka提出的诱导干细胞的四因子,其时在之前已经有多篇网络分析文章提出了这四因子,只是Yamanaka是第一实验验证的。值得一提的是这篇获得诺奖的Cell文章只有一个一作和一个通讯,据説是当时这课题没人看好,也没人帮助,全是一作自己干的。
 
网络绘制工具
 
以OTU共丰度网络为例,常用的工具有MENA, LSA, SparCC和 CoNet四种网络构建方法,可视化的部分常用Cytoscape。详细的翻译见文章:微生物相关网络构建教程中文http://blog.csdn.net/woodcorpse/article/details/73148650
 
网络图的应用
 
1、显示群体内个体间的关系;正/负/无相关
2、显示群体内间的调控关系,具体可分类几个高度相关的子网络/模块;
3、寻找中心结点物种,用于人工重组微生物组实验,验证网络的结构和功能;
4、显示不同群体间共有或特有的OTU;
 
看图实战(Result)
 
示例1. OTU共丰度网络
这篇文章分析了水稻根不同区域的细菌组成,16S分析文章较系统的作品,两年被引用147次,推荐阅读
图5. 共丰度OTU网络,展示与甲烷循环相关的OTU模块
- A. 展示网络分析中11个与甲烷循环相关的模块,图中所有线为Pearson相关系数大于0.6的两个OTU间的连线;
- B. 单独展示与产甲烷菌、营养共生、甲烷氧化菌和其它甲烷循环功能高度相关的模块119,其中与以上四类相关的OTU添加了字母标签,对图例对应;
- C. 展示模块119在不同时间和生态位间的丰度变化
 
1、图中元素解释:以图A为例说明
图中每个点代表一个OTU
图中的连线/边代表相连的两个OTU间Pearson相关系数大于0.6,即有明显的正相关。
图中点的颜色代表该OTU处于不同的模块,图中的图注标示不同模块对应的颜色。
2、图表结果:图中关注了水稻根细菌组相关网络中与甲烷相关的模块;并展示了其中一个模板示例存在多个与甲烷相关的OTU,并发现这类菌在不同生态位和地区的相对丰度存在明显变化。
3、经验和技巧:每个模块一般为相关联接很紧密的一组OTU形成的子网络。在同一套数据中,一般会形成数量从大到小的若干子网络。本图只关注了正相关的关系,其实还存在许多菌/OTU间负相关的关系
 
示例2. 网络图展示不同海拔区域真菌的共有和特有情况
这篇文章是南土所褚海燕老师组杨腾博士2016年发表在Enviromental microbiology上的文章,主要研究了长白山地区不同海拔分布下树木叶片内生真菌的分布,研究真菌的推荐阅读
- 图4. OTU网络图展示不同海拔样品真菌属水平的OTU互作。
1. 图表元素
- 图中每个点代表一个OTU的真菌属,一共有242个属。
- 左侧的点表示这些属只出现在某一特定海拔
- 右侧则表示这些属出现在多个海拔位置。
- 图中间为不同海位置的类,与对应的各组相连;
2. 图表结果:图中展示了不同海拔特有和共有的真菌属数量。其中有62个属出现在所有的海拔位置,而且这些属所占的丰度达到总测序量的98%。
3. 经验和技巧:此图为使用QIIME的make_otu_network.py程序分析获得的点和边文件结果,虽然叫网络,个人认为只是Venn图的变形,还是展示OTU共有和特有的情况,只是图片感觉更高端大气。而且分类的组可自由定义。比如右下角的OTU shared 2-5 组,Venn是做不到的。还有图片的配色,黑底绿线白字是不是很B格满满,有骇客帝国的感觉。

扩增子图表解读8网络图:节点OTU或类Venn比较的更多相关文章

  1. 扩增子分析解读5物种注释 OTU表操作

    本节课程,需要先完成<扩增子分析解读>系列之前的操作 1质控 实验设计 双端序列合并 2提取barcode 质控及样品拆分 切除扩增引物 3格式转换 去冗余 聚类 4去嵌合体 非细菌序列 ...

  2. 扩增子图表解读4曼哈顿图:差异分类级别Taxonomy

    曼哈顿图 Manhattan Plot 曼哈顿图本质上是一个散点图,用于显示大量非零大范围波动数值,最早应用于全基因组关联分析(GWAS)研究展示高度相关位点.它得名源于样式与曼哈顿天际线相似(如下图 ...

  3. 扩增子图表解读5火山图:差异OTU的数量及变化规律

    火山图 Volcano plot 在统计学上,火山图是一种类型的散点图,被用于在大数据中快速鉴定变化.由于它的形成像火山喷发的样子,所以被称为火山图.和上文讲的曼哈顿图类似.   火山图基本元素 火山 ...

  4. 扩增子图表解读3热图:差异菌、OTU及功能

    热图是使用颜色来展示数值矩阵的图形.通常还会结合行.列的聚类分析,以表达实验数据多方面的结果.  热图在生物学领域应用广泛,尤其在高通量测序的结果展示中很流行,如样品-基因表达,样品-OTU相对丰度矩 ...

  5. 扩增子图表解读6韦恩图:比较组间共有和特有OTU或分类单元

    韦恩图 Venn Diagram Venn Diagram,也称韦恩图.维恩图.文氏图,用于显示元素集合重叠区域的图示.   韦图绘制工具 常用R语言的VennDiagram包绘制,输出PDF格式方便 ...

  6. 扩增子图表解读1箱线图:Alpha多样性

    箱线图 箱形图(Box-plot)又称为盒须图.盒式图或箱线图,是一种用作显示一组数据分散情况资料的统计图.因形状如箱子而得名.在宏基因组领域,常用于展示样品组中各样品Alpha多样性的分布 第一种情 ...

  7. 宏基因组扩增子图表解读2散点图:组间整体差异分析(Beta多样性)

    散点图 数据点在直角坐标系平面上的分布图.在宏基因组领域,散点图常用于展示样品组间的Beta多样性,常用的分析方法有主成分分析(PCA),主坐标轴分析(PCoA/MDS)和限制条件的主坐标轴分析(CP ...

  8. 扩增子分析解读6进化树 Alpha Beta多样性

    分析前准备 # 进入工作目录 cd example_PE250 上一节回顾:我们的OTU获得了物种注释,并学习OTU表的各种操作————添加信息,格式转换,筛选信息.   接下来我们学习对OTU序列的 ...

  9. 扩增子分析解读4去嵌合体 非细菌序列 生成代表性序列和OTU表

    本节课程,需要先完成 扩增子分析解读1质控 实验设计 双端序列合并 2提取barcode 质控及样品拆分 切除扩增引物 3格式转换 去冗余 聚类   先看一下扩增子分析的整体流程,从下向上逐层分析 分 ...

随机推荐

  1. Java虚拟机深入JVM内核—原理、诊断与优化视频教程

    http://www.eimhe.com/forum.php?mod=viewthread&tid=142832&highlight=%C4%DA%BA%CB

  2. Oracle RAC load blance

    首先声明 本文基本是阅读 大话RAC 后的笔记.OK, 进入正题. Oracle 10g RAC中采取两种方式提供负载均衡.第一种是connection blance.在用户连接的时候,根据随机算法把 ...

  3. 友盟 个推 微信sdk spring boot

    友盟 个推 微信sdk spring boot ngix mongodb memory cache

  4. HDU 1973

    Prime Path Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others) Total ...

  5. Spark SQL CLI 实现分析

    背景 本文主要介绍了Spark SQL里眼下的CLI实现,代码之后肯定会有不少变动,所以我关注的是比較核心的逻辑.主要是对照了Hive CLI的实现方式,比較Spark SQL在哪块地方做了改动,哪些 ...

  6. 解决华为手机不出现logcat日志的问题

    问题描写叙述:公司一部华为手机在连接Eclipse时在Logcat中看不到相关日志 解决方法:1 进入手机拨号界面2 输入*#*#2846579#*#*3 输入完成后自己主动跳转到測试界面4 依次选择 ...

  7. 自然常数 e 的理解与应用

    某彩票中奖率是百万分之一,则一个人买一百万张彩票仍不中奖的概率是: (1−1106)106≈1e e 往往出现在: 许多微小事件带来的总体变化 随机性和无穷多:

  8. 腾讯云linux服务器安装lnmp一键包

    这边域名已经实名了. 然后修改DNS服务器 然后备案吧 还是先不备案,直接云解析DNS 哦,想起来了,阿里云自己都可以生成SSL证书.重新弄一次吧,其实腾讯云也可以申请域名型免费版DV

  9. PCB genesis SET取中心点--算法实现

    最新ICS工厂有一项incam脚本新需求,这里介绍5种解决方法解决 需求如下图所示:绿色所圈处是是需求出的中心点(图形间距一致归为一类并计算中心点坐标) 前题条件:1.一个SET里面可能有多个CAM, ...

  10. struts2基础学习----struts.xml配置文件

    <?xml version="1.0" encoding="UTF-8" ?> <!DOCTYPE struts PUBLIC "- ...