转录组测序多少生物重复合适?2个?3个?48个?

2016年英国邓迪大学的Geoffrey J Barton教授在RNA发表一篇文章专门评估这一问题。作者对野生型和snf2突变型酵母样品分别测序了48个生物学重复;质控后,野生型样品保留42个生物学重复,突变型样品剩余44个生物学重复。

在控制假阳性率为0.05的标准下,用所有可用生物重复计算出的差异基因定义为该差异基因鉴定工具的金标准差异基因集

作者评估了11个常用的差异基因分析工具,性能最好的是edgeR, DESeq2limma

下面以edgeR为例展示了不同生物学重复数目对鉴定差异基因的影响。

A展示了,在控制相同的假阳性率水平下,不同的生物学重复鉴定出的差异基因数目 (nr: number of biological replicates)。作者从所有生物重复中随机抽取2组、3组、4组…生物学重复,分别计算差异基因,发现:

  1. 差异基因的数目整体与生物重复数量正相关。

  2. 差异基因数目的稳定性与生物重复数量负相关;

    生物重复较少时,不同的抽样导致的差异基因数目波动较大;

    生物重复较多时,检测出的差异基因数目受抽样影响较小,体现在柱状图数据分布更集中 (可视化之为什么要使用箱线图?)。

image

B展示了,不同生物学重复与鉴定的差异基因的真阳性率的关系。不同的实线代表不同的差异基因筛选倍数变化(T=|Log2(FC)|)条件下的真阳性率。虚线代表假阳性率,近乎一条直线,说明edgeR的假阳性率控制的还是比较好的,比较低,且不受生物重复数影响。如果筛选阈值比较高,比如4倍差异(T=2)时,较低的重复数即可获得较高的真阳性率。而筛选阈值较低(T=0)时,真阳性率受生物学重复影响较大;生物学重复越少,真阳性率越低。常规筛选标准2倍差异(T=1)时需要20个生物重复才能达到与4倍差异相同的真阳性率。

image

C则是图B的另一种展现,横轴是筛选倍数阈值 (T=|Log2(FC)|)。蓝色虚线代表3个生物重复条件下的假阳性率,在常规筛选标准2倍差异(T=1)时,假阳性率已趋近于0。不同颜色的实现代表不同生物重复下的真阳性率随筛选阈值差异倍数的变化,整体呈现正相关;且生物重复越多,真阳性率越高,并受筛选阈值影响越少。

D展示了真阳性、真阴性 (非金标准差异基因定义为真阴性(无差异)基因)、假阳性、假阴性基因数目随生物重复数的变化。生物重复越多,漏掉的差异基因(假阴性基因)越少。

image

原文:https://rnajournal.cshlp.org/content/22/6/839.long

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

推荐阅读更多精彩内容