转录组分析结果阅读指南

测序容易,解读困难

现在大部分同学做转录组实验都是让测序公司去分析,经过了那么多年的发展,公司提供的分析结果也是越来越丰富了,光结题报告就有几十页,结果文件有几百兆甚至超过1G。

但是内容多没有指导的话,反而适得其反,很多同学拿到结果研究了几个月,却还是搞不懂各个结果文件夹里到底存放了什么东西。

测序公司分析结果概述:

今天我们就拿一份典型的有参物种转录组结果文件,来给大家概述下各个文件夹里存放了什么,有什么用。建议大家收藏,需要时拿出来看看。

先来看看公司交付结果的组成:

一、前面的10个文件夹储存了本次实验所有的分析结果,写文章时需要用到的统计数据、表格、高清图片都得从这里面找;

二、最下面几个就是结题报告了,主要的用途是比较清楚的描述了整个项目的实验、分析流程,也对部分分析结果进行了概述,强调一下,结题报告里只有部分结果!!想找详细的分析结果还得从上面10个文件夹里找。

结题报告组成

公司一般会同时提供网页版和PDF版的,建议直接看网页版结题报告,当然别用IE打开,不然格式不对,看瞎眼,推荐用Chrome、火狐、360极速浏览器(不是360安全浏览器)等打开。如果只有PDF版的,建议您换家测序公司,都9102年了,连个网页版结题报告都没有,还开啥公司?

看到这里,可能就有同学要喊停了,“等下,我这里怎么只有结题报告,那10个文件夹呢???”赶紧联系那个不负责任的销售吧,十有八九是忘了把结果给你了。这个情况还真不是小编杜撰,我起码碰到过两次,老师项目都做完2个多月了,钱也结清了,结果手里就一个结题报告。

结果文件夹内容概述

下面我们具体讲讲这10个文件夹都存放了什么:

1.data_assess:下机数据的质量评估。每个样品的数据量、reads数、GC含量、Q20、Q30指标,以及碱基含量分布图和碱基测序错误率分布图。

2.Tophat_Cufflinks:测序数据与参考基因组比对结果。有参物种做分析时第一步就是把clean data与参考基因组做比对,这个文件夹里就储存了测序数据在基因组上的覆盖比例、在基因不同区域的比例、插入片段的分布等等和比对相关的文件和图表。

3.gene_expression:基因表达量文件夹。这个好理解,基因定量的一堆相关指标,比如每个样品的表达量、表达密度、表达箱型图、样品间表达量相关性等等都在这里了。

4.Saturation:基因定量饱和度结果。这个文件夹里东西很少,就是为了评估样品测序数据量是否足够而做了一个饱和度图,没啥别的了,我简直怀疑这个文件夹是测序公司拿来凑整的。

5.Gene_Func_Anno:基因注释结果。虽然有参物种的基因组本身有功能注释,但往往不够全面。测序公司还会对本次实验检测到的基因在NR,Swiss-Prot,GO,COG,KOG,Pfam,KEGG等数据库里进行功能注释,这样后续才能全面了解感兴趣基因的作用。这个文件夹里又分成AllGene和NewGene两部分,分别是所有的基因和新注释的基因,大家可以根据自己的研究目的来进行选择查看。

6.DEG_Analysis:差异基因分析结果。大家要注意了,这个文件夹里储存的往往就是大部分同学最关注的部分了。各比对组的差异基因火山图、聚类分析、venn图、注释结果等等都在这里了。

7.DEG_PPI:差异基因蛋白质互作网络结果。应用STRING蛋白质互作数据库(http://string-db.org/)中的互作关系进行差异基因蛋白互作网络的分析。

8.Altsplice_Analysis:可变剪切结果。储存了各个样品可变剪切事件的数目、样品间差异可变剪切事件统计。

9.SNP_Analysis:SNP/InDel结果文件夹。转录组不仅提供了基因表达量的分析,还同时提供了各个样品的SNP和InDel信息,比如SNP、InDel在基因组各个区域的分布情况,InDel长度统计,SNP突变类型统计等。

10.Gene_Structure_Optimize:基因结构优化文件夹。这个文件夹里内容也很少,根据此次转录组测序结果对原参考基因组的基因区范围进行了优化。

深入学习方法推荐

好了,看到这里相信同学们对公司交付的转录组结果有了一个概况性的认识了,当然结果文件有几百个,彻底理解、搞懂全部的分析内容还是需要再花费功夫的。想深入学习的同学可以选择观看我们制作的《转录组分析结果的解读》视频课程,来更好的利用手中的数据!

转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接转录组(有参)结果解读转录组(无参)结果解读

更多生物信息课程:

1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4分文章简单快速,学习链接:基因家族分析实操课程基因家族文献思路解读

2. 转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接:转录组(有参)结果解读转录组(无参)结果解读

3. 转录组数据深入挖掘技能-WGCNA,提升你的文章档次,学习链接:WGCNA-加权基因共表达网络分析

4. 转录组数据怎么挖掘?学习链接:转录组标准分析后的数据挖掘转录组文献解读

5. 微生物16S/ITS/18S分析原理及结果解读OTU网络图绘制cytoscape与网络图绘制课程

6. 生物信息入门到精通必修基础课,学习链接:linux系统使用perl入门到精通perl语言高级R语言画图

7. 医学相关数据挖掘课程,不用做实验也能发文章,学习链接:TCGA-差异基因分析GEO芯片数据挖掘GSEA富集分析课程TCGA临床数据生存分析TCGA-转录因子分析TCGA-ceRNA调控网络分析

8.其他课程链接:二代测序转录组数据自主分析NCBI数据上传二代测序数据解读

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

推荐阅读更多精彩内容