

转录组测序技术(RNA-seq)作为目前二代测序领域最普遍的技术手段,自从转录组测序问世以来,已经开发了数百种分析工具。根据转录组分析内容可大致将其分析流程分为比对,转录本组装,差异表达分析和差异基因注释。目前,分析的每一步都有很多软件,其软件的性能和分析效率不尽相同。
上篇文章小编为各位小伙伴介绍了转录组分析的第五步——新转录本预测,新转录本预测主要是对参考基因组信息进行补充,下面小编来介绍转录组分析的第六步——各种图形的绘制吧~我们在进行转录组数据分析时,除了会得到各种结果文件外还需要将结果文件绘制成各种图形,以展示分析结果。本文主要介绍转录组分析中常见的几种图形的绘制。

1 密度图密度图主要用于展示各个样本表达量水平的分布。输入文件为各个样本的各个基因的FPKM值。
输入文件格式如下: 
密度图绘制代码如下: 
密度图绘制如下:

2 柱状图柱状图主要展示上调差异表达基因和下调差异表达基因数量的差异。输入文件为上下调基因的数量。
输入文件格式如下:

柱状图绘制代码如下: 
柱状图绘制如下:

3聚类热图聚类热图可用于可视化各个样本基因表达量的差异,表达模式相近的样本归为一类。输入文件为各个样本的各个基因的FPKM值。
输入文件格式如下: 
聚类热图绘制代码如下:

聚类热图绘制如下:

4韦恩图韦恩图可用于展示两组至多组共有和特有的基因,输入文件为各样本鉴定到的基因ID。
输入文件格式为:

韦恩图绘制代码如下: 
韦恩图绘制如下:

5经验总结
本文主要介绍了几种常见的转录组分析的图形绘制。绘制图形主要使用R进行绘制,R的下载和安装方法小伙伴可参考https://www.r-project.org/。建议安装R4.0版本以上。绘制图形时需要安装一些所需的R包,可使用install.packages(“R包名称”)进行安装。至此,转录本分析的第六步——绘图就介绍完了,希望对各位小伙伴有所帮助,我们下一期再见吧~
注:此推文未经许可禁止转载! 阅读推荐:生信实操|一个生信素人的上道经验分享-转录组测序(新转录本鉴定篇)
一个生信素人的上道经验分享-转录组测序(差异分析篇)
一个生信素人的上道经验分享-转录组测序(基因定量篇)
一个生信菜鸟的上道经验分享-转录组测序(转录本组装篇)
一个生信菜鸟的上道经验分享-转录组测序(比对篇)
