本篇内容主要讲解“如何使用TCGAbiolinks进行甲基化和转录组数据的联合”，感兴趣的朋友不妨来看看。本文介绍的方法操作简单快捷，实用性强。下面就让小编来带大家学习“如何使用TCGAbiolinks进行甲基化和转录组数据的联合”吧!
DNA甲基化作为表观遗传的一种重要标记，在基因表达调控中发挥重要作用，已有研究表明，位于启动子区的甲基化会抑制基因表达。结合甲基化数据和基因表达谱数据，可以更好的分析甲基化的调控功能。
通过TCGAbiolinks不仅仅可以对某种类型的数据进行分析，还可以进行多组学数据的联合分析，DNA甲基化和基因表达谱数据的联合分析过程如下
以肾上腺皮质癌为例，首先从TCGA数据库下载450K甲基化芯片的数据，分析后得到差异甲基化位点
差异甲基化位点的火山图结果如下
首先下载数据
然后进行差异分析，代码如下
差异基因火山图的绘制代码如下
效果图如下所示
starburst plot是联合分析最常用的可视化方法，代码如下
生成的效果图如下
这里采用了差异分析的FDR值来绘制图片，其他文章中也有采用fold change来绘制该图的，因为根据fold change的大小更好判断差异表达的趋势，根据FDR值只能判断是否显著差异。
在DNA甲基化和基因表达谱的联合分析中，鉴于DNA甲基化抑制基因表达的思想，通常关注二者间为负相关的位点，即甲基化下调而基因上调，甲基化下调基因上调的相关位点。在上图中，显著差异且存在负相关的位点用圆圈进行了标注。
通过TCGAbiolinks可以轻松实现DNA甲基化和基因表达谱数据的联合分析，美中不足的是它的starburst plot没有采用常用的fold change来绘制。
到此，相信大家对“如何使用TCGAbiolinks进行甲基化和转录组数据的联合”有了更深的了解，不妨来实际操作一番吧！这里是亿速云网站，更多相关内容可以进入相关频道进行查询，关注我们，继续学习！
免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
