转录组分析流程

概念及类型:

  • 转录组(transcriptome)广义上指某一生理条件下,细胞内所有转录产物的集合,包括信使RNA核糖体RNA转运RNA非编码RNA;狭义上指所有mRNA的集合。

  • 转录组测序分析可以分为有参转录组分析和无参转录组分析。有参无参的意思是,有/无参考基因组。

1、有参转录组分析流程:

有参转录组分析流程.png

2、无参转录组分析流程:
Snipaste_2019-09-11_09-32-27.png

clean reads :为确保Reads有足够高的质量,将下机原始测序数据(raw reads)去掉含有带接头的、低质量的reads,得到clean reads,保证后续分析的准确性。
unigene:针对无参考基因组的物种,将小片段拼接出unigene,构建参考序列,以便后续分析,是研究无参物种分子机制与调控网络的有效手段。

三种转录组分析模式

基于研究物种是否有参考基因组,实验目的上是否需要分析新的转录本,转录组测序的分析模式大致可以分成3种类型,如下图:


553.jpg

一、有参,需要分析新转录本
部分有参考基因组的物种,由于注释信息不够完善,或需要分析一些非编码RNA,这时需要基于Reads与基因组比对信息对转录组进行组装,以期获得新的转录本来让分析结果更加完备。

这也是有参物种做转录组最常用的分析模式,其分析步骤如下:

  1. Reads与基因组比对

  2. 基于比对信息组装转录本

  3. 基因或转录本表达定量

  4. 差异分析和功能富集分析

二、有参,只分析已知转录本
针对参考基因组注释信息较为详细的物种,比如人,小鼠,拟南芥等模式生物,同时您的实验目的很明确,就是分析已知的基因或转录本,那就可以直接基于基因组注释信息中提取出的转录本序列来进行后续分析。该分析模式分析流程简单、速度快,其具体分析步骤如下:

  1. Reads与转录本序列进行比对

  2. 转录本表达定量

  3. 差异分析和功能富集分析

三、无参考基因组的转录组
而对于没有参考基因组的物种,或者基因组组装不好的物种,必须先使用测序数据组装一套转录本,再基于转录本进行后续分析。其分析步骤如下:

  1. Reads De novo组装转录本序列

  2. Reads 回比组装好的转录本序列

  3. 转录本表达定量

  4. 差异表达分析和功能分析


other:


image.png

image.png

image.png

怎么做无参转录组分析?


如下图1所示,在做转录组分析时,除了需要reads等原始数据外,还需要一系列的注释文件。然而,当我们的研究物种没有参考基因组但需要做转录组分析时,便不能通过以上流程获得我们需要的分析结果。


图1-有参转录组分析流程

此时,我们可以通过拼接reads数据获得转录本,该转录本用于后续reads比对。再将该拼接好的转录本与各大蛋白数据库注释,其中常用的数据库有NR、NT、Swiss-Prot、KEGG、COG、GO等。获得该注释结果后,我们便对该物种的蛋白信息有一个比较清楚的了解,可以进行后续的转录组分析流程。

最常使用的组装软件为Trinity,运行该软件时,主要分三部分运行,依次为Inchworm、Chrysalis和Butterfly(下图2)。


图 2 Trinity组装流程简图

运行Inchworm时,其目的是将每条短读长的reads打断成固定长度的Kmer(默认为25 bp)形成一个Kmer库,再将所有Kmer按其出现的次数进行排序,选择出现次数最高的Kmer进行3端延伸一个碱基(A、T、G、C),延伸后统计靠近3端25bp的Kmer在库中出现的次数,选择次数最高的那条路径(如有多条路径出现的次数一样,则几条路径保留),继续向3端延伸至不能延伸为止。同理进行5端延伸,最终形成一定长度的Contigs。

每形成一个Contigs时,将形成该Contigs的kmer从Kmer库中去除掉,继续选择次数最高的Kmer进行延伸。最终用完Kmer库中的Kmer,将Contigs放在一起形成一个contigs库。运行Chrysalis时,对Contigs库中的Contigs按照一定条件进行聚类,之后对每类Contigs构建de bruijin graphs。

最后运行Butterfly,解析上一步构建的de bruijin graphs,形成转录本。将所有转录本输出到一个文件,便是最Trinity终组装的文件,其格式如图3所示。


图 3 Trinity组装的转录本格式

获得拼接后的转录本后,便可以按照Trinity官网中的说明,进行注释和下游分析,其流程如图4所示。


图 4 Trinity官网中的分析流程

参考文章:

易基因-无参转录组分析流程

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 202,980评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,178评论 2 380
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 149,868评论 0 336
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,498评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,492评论 5 364
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,521评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,910评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,569评论 0 256
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,793评论 1 296
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,559评论 2 319
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,639评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,342评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,931评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,904评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,144评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,833评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,350评论 2 342

推荐阅读更多精彩内容