r语言 转录本结构及丰度_生信实操|一个生信素人的上道经验分享转录组测序(绘图篇)...

c4807bec9fd78d0a5dd0f2ab02a9e978.gif

a5de2cc26d0ca637f791656d4df1bbee.png

转录组测序技术(RNA-seq)作为目前二代测序领域最普遍的技术手段,自从转录组测序问世以来,已经开发了数百种分析工具。根据转录组分析内容可大致将其分析流程分为比对,转录本组装,差异表达分析和差异基因注释。目前,分析的每一步都有很多软件,其软件的性能和分析效率不尽相同。

上篇文章小编为各位小伙伴介绍了转录组分析的第五步——新转录本预测,新转录本预测主要是对参考基因组信息进行补充,下面小编来介绍转录组分析的第六步——各种图形的绘制吧~我们在进行转录组数据分析时,除了会得到各种结果文件外还需要将结果文件绘制成各种图形,以展示分析结果。本文主要介绍转录组分析中常见的几种图形的绘制。

b9aca10ab4aaae2d00027ace9e28ddcf.png042738ecfc414765f34f7137e52e338c.gif1 密度图

密度图主要用于展示各个样本表达量水平的分布。输入文件为各个样本的各个基因的FPKM值。

输入文件格式如下: 83ef91d1abc9374b690c35dab43ad6fa.png

密度图绘制代码如下: ef18d769c0523dc7f3b4aabc1a832d3d.png

密度图绘制如下:

 1cb37655d4effe40448b4358723836fa.png

042738ecfc414765f34f7137e52e338c.gif2 柱状图

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

输入文件格式如下:

 d3f397c533211a951f69f27828c5ef50.png

柱状图绘制代码如下: 53046555a33ae76425520ed298b27618.png

柱状图绘制如下:

 110d80a14252fd481d41a425d9dce399.png

042738ecfc414765f34f7137e52e338c.gif3聚类热图

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

输入文件格式如下: 83ef91d1abc9374b690c35dab43ad6fa.png

聚类热图绘制代码如下:

1b77e7a82163eb693ad8f62fcf8f219e.png

聚类热图绘制如下:

 df646c01c16829665b0faaaef7b95069.png

042738ecfc414765f34f7137e52e338c.gif4韦恩图

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

输入文件格式为:

 7d199a1baf0b3d248740d9ee997f19dc.png

韦恩图绘制代码如下: 16881bd6cb63c3e1b79d7846f10c212c.png

韦恩图绘制如下:

 cf1faf64de131366877a247f917829c9.png

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

  • 一个生信素人的上道经验分享-转录组测序(差异分析篇)

  • 一个生信素人的上道经验分享-转录组测序(基因定量篇)

  • 一个生信菜鸟的上道经验分享-转录组测序(转录本组装篇)

  • 一个生信菜鸟的上道经验分享-转录组测序(比对篇)

6a607f3385bb4d8665372ff1f5ca9894.gif