动植物重测序--全基因组关联分析GWAS

GWAS简介

GWAS(Genome-wide association study)是对遗传多样性丰富的自然群体的每个个体进行基因组测序,结合目标性状的表型数据,基于一定的统计方法进行全基因组关联分析,可以快速获得影响目标性状表型变异的染色体区段或基因位点。

当然,GWAS可以应用于人的表型分析,这里暂时先说动植物的。

GWAS已经发表的物种:玉米水稻拟南芥大豆毛果杨番茄果蝇白虎疟原虫等,物种很丰富。因为不需要构建家系群体,像BSA遗传图谱还得杂交自交,作物植物之类的只杂交养大就得小半年时间,GWAS就直接采集自然群体就好也可以多地采样,这样变异范围广一些,做出来的结果适用性也比较强;其次就是一次性可以定位多个性状,基本是唯一的一个用于自然样本性状定位的分析方法;定位精度也高,可以直接获得与目标性状相关的基因。

GWAS方案策略

①样本选择:自然群体,可以多品种多地域取样,原则上样本间不能有明显的亚群分化(例如生殖隔离等)即遗传背景差异小,如果样本间存在亚群分化,则每个亚群建议样本量在200个以上;

②样本数量:≥300个个体,每种性状能记录到表型的植株数尽量保证在200株以上;

③测序策略:有参选择WGS,无参选择简化基因组,测序数据量10x,如果样本数量很多500个样本以上,测5x即可。

※※另外非常重要的一点,前期对采集样本的表型记录一定要准确详细,因为分析中要用到表型数据所以表型的记录也会对分析结果产生影响,原因:每次计算为单表型 GWAS 分析,在当前表型中,缺失样本不参与此次分析,导致有效个体数目减少;GWAS 通过统计模型找出影响表型差异的 SNP 及基因,记录不准确的话,可能会导致找到的SNP 或基因有很大的假阳性。

表型选择:生育期 叶子果实大小 枝茎长短弯曲程度 花瓣果实品质含糖量 抗病抗虫抗逆性,除此之外,基因的表达量也可作为表型分析叫做eWAS,代谢物类型也可叫做mWAS,只有想不到没有做不到。

样本数量对GWAS结果的影响

GWAS分析


GWAS分析流程

分析流程其他软文有很多讲过的,我就不班门弄斧了,我把链接放在文末供参考。

GWAS常见问题

1、家系样本是否可以进行 GWAS 分析?

通常禽类研究选择半同胞家系(由同父异母或异父同母所生子女的集合体成为半同胞家系)或全同胞家系(由同父同母所生子女的集合体称为全同胞家系)。理论上,其他动物也可以选择表型多样性丰富的家系样本进行GWAS 研究, 优先推荐选用 F2代群体。如果老师的 F2 代群体样本数少,也可以将 F1 代和 F2 代作为一个研究群体,在后续GWAS 分析过程中会有 KINSHIP(亲缘关系矩阵)去校正模型。

2、是否 GWAS 样本的表型分布必须呈现正态分布?

对于数量性状所选择的个体,在表型上尽量呈正态分布;若无法取到完全正态分布的样本的性状, 针对该类性状,在后续分析过程中采用非正态分布性状的分析方法对其进行校正,以保证分析结果的可靠性。

3、是否一定要超过 200 个样本才能进行 GWAS 分析?

基于 GWAS 文章中的样本数和我们的项目经验,推荐选取 200 个以上的样本,以尽可能保证分析结果的可靠性。但是对于一些珍稀的物种研究,材料相对稀少且难以获得,如老虎,熊猫,金丝猴等,可适当减少样本个数。 如果研究内容是由单基因控制的性状,其实材料少也可以得到不错的结果。但如果是复杂性状,少量样本的关联分析结果不可信。

4、如何降低分析结果的假阳性?

在全基因组关联分析中,前期对样本的采集情况(表型分布均匀,环境一致)会对后续分析的假阳性结果存在最大的影响,在分析过程中会采用如下方法降低分析结果的假阳性:a)结合群体分层信息,利用混合线性模型,对结果进行校正; b)必要时还会采取多种线性模型进行分析降低假阳性;c)采用 Bonferroni 校正法来校正 GWAS 分析中多重假设检验后的 P 值可以降低假阳性的概率。因此,前期对样品的采集及表型的记录需要准确详细。

5、得到的结果如何验证?

a)已发表的本物种或相近物种的QTL或基因比较分析;b)若没有已发表的,可通过家系群体QTL分析进行验证;c)结合转录组或者表达谱分析,验证相关区域或者位点调控的通路信息。

参考学习:

1、全基因组关联分析 (GWAS) - 简介

2、全基因组关联分析学习资料(GWAS tutorial)

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

推荐阅读更多精彩内容