WGCNA使用非常广泛,网上流传的示例代码也非常多,但还是推荐大家直接参考官方代码,最全面且权威:https://horvath.genetics.ucla.edu/html/CoexpressionNetwork/Rpackages/WGCNA/Tutorials/index.html

今天不讲代码,只讲讲这些图怎么看,能表达什么意义。

图解|WGCNA那么多图,都啥意思? - 知乎

WGCNA文章[1]给出了工作流程图,主要包括:构建共表达网络→模块识别→模块与性状关联→模块相关性→关键模块中识别核心基因

在构建共表达网络之前,还有一步比较关键的步骤——寻找最优软阈值(soft thresholding或power),使构建的网络更符合无标度拓扑结构。这时,我们会看到如下图:

确定最优软阈值主要参考左图,即在不同软阈值(x轴)情况下的无标度拟合指数(scale-free fit index,y轴)。其中红线表示主观选择的无标度拟合指数取值,如图为0.9,一般我们选择0.8~0.95之间较为稳妥。从左图看,当无标度拟合指数在0.9时,符合构建无标度网络的最小软阈值为6,因此可以选择6作为最优软阈值用于后续分析。右图为不同软阈值情况下的网络连通度。

基于最优软阈值构建共表达网络,将基因划分到不同模块后,可以绘制基因聚类树:

这个图可以分为两部分看:上半部分是基因的层次聚类树状图,下半部分是基因模块,也就是网络模块。上下对应,可以看到距离较近的基因(聚类到同一条分支)被划分到了同一模块。

接着,计算模块与感兴趣的临床特征之间的相关性与显著性,绘制相关性热图:

上图中,最左侧的颜色块代表模块,最右侧的颜色条代表相关性范围。中间部分的热图中,颜色越深相关性越高,红色表示正相关,绿色表示负相关;每个单元格中的数字表示相关性和显著性。如上图,brown模块与weight性状表现为正相关且相关性最高。此时,我们可以选择相关性最高的brown模块作为关键模块。一般,我们会按相关性的绝对值筛选最相关模块,即负相关模块也应该考虑在内。需要注意的是,grey模块中包含了所有未参与聚类的基因,因此是无效模块,不应用于后续分析。

基于拓扑重叠矩阵,可以绘制基因之间的相关性热图。其中颜色越深,说明基因之间的相互作用越强。由于对角线表示模块内部基因之间的相互作用,所以自然是对角线上的颜色最深:

也可以看看模块之间的聚类树和相关性热图,探索模块之间的互作关系:

可以将感兴趣的性状,如weight也加入到图中。如上图,可以发现red、brown、blue模块之间有很强的相关性,甚至比他们与weight之间的相关性更强。

之后针对关键模块brown和感兴趣的性状weight进一步挖掘,看看基因与模块的相关性(Module Membership, MM)和基因与性状的相关性(Gene Significance, GS)之间是否有某种关联。 `

通过以上散点图,可以发现MM和GS呈正相关,说明这些与性状高度相关的基因,在关键模块中也扮演着举足轻重的角色。

我们也可以基于关键模块中基因之间的相似性,构建关键模块网络,以Cytoscape或VisANT软件可视化展示网络中基因之间的互作关系:

最后,也是最重要的——基于关键模块筛选核心基因。这一步的可用方法比较多,没有固定限制。比如,可以根据MM>0.8且GS>0.3筛选核心基因[2];可以根据关键模块的网络degree筛选top 30作为核心基因[3];也可以计算kME值,依据|kME|≥0.7筛选[4];或者将网络导入Cytoscape后,利用插件cytoHubba筛选[5]。

1. Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559. Published 2008 Dec 29. doi:10.1186/1471-2105-9-559

2. Song ZY, Chao F, Zhuo Z, Ma Z, Li W, Chen G. Identification of hub genes in prostate cancer using robust rank aggregation and weighted gene co-expression network analysis. Aging (Albany NY). 2019;11(13):4736-4756. doi:10.18632/aging.102087

3. Zhang X, Feng H, Li Z, et al. Application of weighted gene co-expression network analysis to identify key modules and hub genes in oral squamous cell carcinoma tumorigenesis. Onco Targets Ther. 2018;11:6001-6021. Published 2018 Sep 19. doi:10.2147/OTT.S171791

4. Panahi B, Hejazi MA. Weighted gene co-expression network analysis of the salt-responsive transcriptomes reveals novel hub genes in green halophytic microalgae Dunaliella salina. Sci Rep. 2021;11(1):1607. Published 2021 Jan 15. doi:10.1038/s41598-020-80945-3

5. Li CY, Cai JH, Tsai JJP, Wang CCN. Identification of Hub Genes Associated With Development of Head and Neck Squamous Cell Carcinoma by Integrated Bioinformatics Analysis. Front Oncol. 2020;10:681. Published 2020 May 22. doi:10.3389/fonc.2020.00681

WGCNA那么多图,都啥意思? 官网相关推荐

  1. 为什么现在越来越多的人买手机都喜欢从官网购买,而不从实体店买?

    准确的说,现在越来越多的人喜欢在淘宝网.天猫网或是京东网等购物平台购买手机(官网是企业公司为宣传自己的产品而创办的独立网页).这些网站的旗舰店或者是工厂直营店大多数都支持七天无理由退换货,还可享受无人 ...

  2. 南都娱乐周刊官网突然崩溃 疑遭黑客攻击

    南都娱乐周刊近期爆出的"文章出轨恋姚笛"迅速火遍网络,并创造"周一见"."且行且珍惜"等热点词语,当事人文章的道歉微博甚至超越了王菲离婚微博 ...

  3. 左右全屏banner焦点图 代码特效+苹果官网首页左右全屏banner焦点图效果+包括JS图片CSS样式等

    介绍 源码名称:[左右全屏banner焦点图]代码特效+苹果官网首页左右全屏banner焦点图效果+包括JS图片CSS样式等 源码大小:16.6KB 开发语言:PHP+Mysql 操作系统:Windo ...

  4. 华硕FA506IV重装系统以后驱动都打过官网的驱动,无线网还是不能用的解决办法。

    华硕天选FA506IV,重装过系统无线网还是没有信号,无线网卡驱动也在官网更新过了还是不能用,原来bios版本是315,在官网更新BIOS版本到316重启下电脑直接就可以了!更新BIOS版本记着插上原 ...

  5. 表单内如何直接贴图而不用上传图片_重磅更新|偷偷告诉你,表单大师官网改版啦啦啦啦...

    你们心心念念的 版本 更新 终(ke)于(suan) 来了(!) 昨天表单大师新版本悄咪咪的上线啦,相信眼尖的小伙伴们已经发现官网首页大变样~偷偷告诉你们哦,本次不仅仅只有更新,还上线了许多便宜好用的 ...

  6. 基于django和vue的xdh官网设计

    前言 本项目是使用三段分离的设计 前台 使用materialize框架搭建的前台页面,后端使用的django写的接口 后台 使用Amazon UI 模板搭建的界面,管理各个部分的内容 项目环境 pyt ...

  7. NVIDIA Jetson官网资料整理

    可看到,对比树莓派相关方面的资料,NVIDIA Jetson的资料不多,特别是一些国内的培训机构,这方面系统性的介绍也几乎没有,所以我们应该围绕官网转.博主之前学习TensorRT以及OpenVINO ...

  8. Gstore官网学习二:安装(笔者自带填坑)

    文章目录 一.系统要求 二.安装环境准备 2.1 安装gcc 和 g++ 安装 2.1.1官网指导 2.1.2 笔者方法 第一步:安装centos-release-scl 第二步:安装devtools ...

  9. Bootstrap框架: 模拟小风车科技官网

    插件地址:http://www.bootcss.com Bootstrap是一个简洁.直观.强悍的前端开发框架,让web开发更迅速.简单. 本次我们利用它来模拟小风车科技的官网,因为它的官网开发也是利 ...

  10. 同样的电器,为什么官网能比实体店的价格便宜那么多?

    所谓羊毛出在羊身上,实体店的价格之所以比网店更贵,因为实体店的成本更高. 对于官网上所卖的电器,它的价格构成包括生产成本+运营成本(宣传.人工等费用).物流成本+税费+企业利润. 而对于实体店来说,它 ...

最新文章

  1. 【Sql Server】DateBase-SQL调整
  2. 二级域名怎么设置_怎么建设自己的网站?看完这四个步骤你就明白了
  3. java组件代码_[java]常用组件
  4. [react] 为什么属性使用className而不是class呢?
  5. [学习笔记]半平面交
  6. 曾为浮云遮望眼,平沙莽莽黄入天。 终有一缘得一道,拨开云雾见青天! ——论如何快速删除上亿级数据集思路
  7. 宅在家里写数据库中DML(增删改)
  8. jeecms 数据被引用,不能删除。解决办法
  9. 自走棋投降代码_自走棋太火惨遭针对!代码哥之后又是锁钱挂?玩家:还有拒绝员!...
  10. 函数指针以及在类成员函数中应用函数指针
  11. dart 语言中的列表(List)
  12. 通过Spring Boot中的手动Bean定义提高启动性能
  13. SmartSvn遗漏.a文件的解决办法
  14. 如何恢复删除的文件?wps文件恢复,4种方法教你找回来
  15. python音频 降噪_一种基于深度神经网络的音频降噪方法技术
  16. 小米笔记本android,全球唯一 独家定制:「安卓中国版」小米笔记本电脑Air图赏...
  17. 【TP5】安装Guzzle并简单使用
  18. 违反协议 oracle,Oracle增加表字段导致“违反协议”异常
  19. 企业业务架构设计方法论及实践(二)
  20. 虚拟机2012搭建DNS服务器,Windows Server2012 安装配置DNS服务器方法详解

热门文章

  1. 李洪超 硬件工程师_工作六年,如何从 layout 转硬件工程师?丨版主采访(一)...
  2. 在 iPad和 iPhone的浏览器上查看网页源代码
  3. 生产任务分配问题 matlab+lingo
  4. 项目启动报错No appropriate protocol (protocol is disabled or cipher suites are inappropriate) 解决办法
  5. Android获取本机号码(双卡双待无法获取两个号码)
  6. make出错,编译通不过的解决办法
  7. JqGrid 表格基本使用(一)
  8. canvas二次、三次贝赛尔曲线 绘制气泡圆角矩形立体桃心
  9. ai如何做倒角和圆角_AI怎么将直角矩形设置成圆角-AI将直角矩形设为圆角的方法 - 河东软件园...
  10. 关键路径例题图表_关键路径法典型范例