多变量阈值 分子遗传数据 Gibbs

  1. Stock KF, Distl O, Hoeschele I. Bayesian estimation of genetic parameters for multivariate threshold and continuous phenotypes and molecular genetic data in simulated horse populations using Gibbs sampling. BMC Genet. [Internet]. 2007;8:19. Available from: http://dx.doi.org/10.1186/1471-2156-8-19

使用Gibbs采样的模拟马群中多变量阈值和连续表型和分子遗传数据的遗传参数的贝叶斯估计

抽象

背景
成功实施包括表型和基因型信息的多变量动物阈值模型的要求尚不清楚。这里使用模拟马数据来调查遗传参数的遗传参数的多变量估计器的性质在分类,连续和分子遗传数据在重要的放射性健康性状的背景下使用混合线性阈值动物模型通过Gibbs采样。模拟谱系包括7代和40000动物每代。模拟了四个二元性状的一个连续性状和负债liabilities的加性遗传值,残差​​和固定效应,类似于Warmblood马中遇到的情况。对一种负债模拟了数量性状基因座(QTL)效应和遗传标记信息。研究了遗传标记与QTL之间的重组率和遗传标记多态性信息含量的不同情况。对于每种情况,从模拟群体中抽取10个重复,并且在每个重复中产生具有性状记录和遗传标志信息可用性的动物的数量和分布不同的6个不同数据集。 (Co)方差分量使用​​贝叶斯混合线性阈值动物模型通过Gibbs采样进行估计。残差方差固定为零,适当的先验用于遗传协方差矩阵。

结果
有效样本大小(ESS)和遗传参数的偏差在数据集之间显着不同。遗传力估计的偏差对于连续性状为-6%至+ 6%,对于中度遗传性的二元性状为-6%至+ 10%,对于低遗传力的二元性状为-21%至+ 25%。加性遗传相关性大多低估了低遗传力的连续性状和二元性状之间,在连续性状和中度遗传性的二元性状之间低估或高估,并且在两个二元性状之间高估。使用性状信息对后续两代动物增加了ESS和减少参数估计的偏差,而不仅仅是增加来自一代的信息性动物的数量。将基因型信息考虑为模型中的固定效应导致高估QTL性状的多基因遗传性,但增加QTL性状的加性遗传相关性的准确性**。

结论
联合使用表型和基因型信息对父母和后代将有助于确定感兴趣的性状之间的激动和拮抗遗传相关性,促进有效的多特征选择方案的设计


背景

使用线性模型估计分类性状的遗传参数违反了混合线性模型方法的基本假设。已经开发了算法以将线性模型估计转换为潜在的责任尺度,以补偿由线性模型中的非线性性状的分析所导致的估计偏差,但是变换的遗传参数估计可能仍然是显着偏差的[例如, 2,3,4]]。阈值模型用于估计遗传参数的使用直接说明了分类性状的非线性性质,并且阈值模型估计应当比线性模型估计或变换的线性模型估计更可靠[例如[5,6,7]]。

马尔科夫链蒙特卡罗(MCMC)方法如吉布斯取样(GS)使得实现多变量阈值模型或多变量混合线性阈值模型成为可能。动物模型充分利用可用的谱系信息,但遗传方差和协方差估计的准确性以及GS链的收敛可能在低性状流行率和每个个体的少量观察的情况下是一个问题[例如[5,8,9,10] ]。因此,在实际情况下,多变量动物阈值模型的实现并不总是直接的。连续性状的包含,即使用多元线性阈值模型,预期增加遗传参数估计的可靠性[例如,[2,11,12,13]]。

在动物育种中,健康数据通常使用离散类别记录,而大多数性能特征是连续的。在马中,二进制编码已经用于对马on的放射学检查的研究,并且在年轻马的肢体中已经确定了放射学可见的改变的高流行性,主要在10-25%的范围内。 [14,15,16,17]]。这促进了寻求预防而不是治疗措施。由于马运动系统的强度和健壮性在马业的所有部门都非常重要,因此已经提出在Warmblood骑乘马的当前育种计划中包括放射健康性状[18,19]。这些分类性状的可靠估计的遗传参数提供了这样做的基础。对于分析的射线照相健康数据,根据Dempster和Lerner [20]和Vinson等人的转换因子的适用性[21]通过残差最大似然(REML)获得的线性动物模型估计已经通过模拟证明[22]。然而,由线性模型中的二元性状的分析引起的并且必须通过变换进行补偿的过估计或低估速率取决于数据结构。因此,需要重新评估变换过程,以分析关于可用信息的分布和种类的不同结构的数据。用于(共)方差分量估计的贝叶斯多变量动物阈值或混合线性阈值模型的实现可以提供有价值的替代方案

已经为不同物种的生产和健康性状鉴定了定量性状基因座(QTL),即包括影响个体关于特定性状的表型的基因的基因组区域[23]。增加对遗传决定放射性健康特性的知识[24,25]意味着在马中改进遗传评价和选择计划的进一步机会有效利用标记辅助选择的条件已经被描述[26],但是结合使用表型和基因型数据对分类性状的遗传参数估计的准确性的影响还没有被研究成功实施包括表型和基因型信息的多变量动物阈值模型的要求是未知的

本研究的目的是使用线性阈值动物模型和Gibbs采样来表征分类,连续和分子遗传数据的遗传参数的多变量估计的性质。在Warmblood马的重要放射健康性状的上下文中研究了数据结构和分子遗传数据的质量对遗传参数估计的准确性的影响。

统计分析
使用用于动物模型的多特征Gibbs采样器(MTGSAM)的阈值版本的Gibbs采样来评估遗传参数[35],该软件支持连续和分类性状的任何组合的多变量遗传分析。随机和残差效应被假定为正态分布,平均先验被用于固定效应。用户可以为加性遗传和残差(共)方差矩阵指定起始值和先验。对于我们的分析,为所有加性遗传变异选择起始值1,对所有加性遗传协方差选择起始值0,并且所有性状之间的残余协方差固定为零。在单变量和多变量二进制阈值模型中,通过将阈值和残差方差分别固定为零和一的值来确保模型的可识别性[36]。用于限制固定在一个的对角线元素的残余协方差矩阵的有效抽样的方法仍在开发中。然而,本研究的残差协方差的固定由以前的实际数据分析的结果证明[14],并拟合我们的模拟数据。因为残差协方差在实际数据中可以忽略,因此在数据模拟中设置为零,进一步信息的可能增益似乎与残差协方差的采样的额外成本不成比例。因此,在本研究中,残余协方差被固定为零。对于遗传协方差矩阵,使用具有最小形状参数(即υIW= 7)的逆Wishart分布的适当的先验,以确保后验适当性。在所有分析中考虑了性别和当代组的固定效应。在仅对数据集C1和C2的分析中考虑标记基因型的固定效应,区分所有遗传标记的不利等位基因纯合阴性的个体,至少一个遗传标记的不利等位基因杂合的个体和纯合的个体对于四种遗传标记中的至少一种的不利等位基因。

在所有分析中,吉布斯链的总长度设定为205000,并且保存5000轮的老化后的所有样品。通过样品图的目视检查来检查Gibbs链的收缩以及对要丢弃的加性轮次的老化的需要。通过在后Gibbs分析程序POSTGIBBSF90 [37]中实施的时间序列方法,以10的稀释率计算所有(共)方差估计值的有效样本大小(ESS)和蒙特卡罗误差(MCE)。未稀释的链用于计算加性遗传(共)方差,遗传性和加性遗传相关性估计的后验均值。选择后验手段而不是模式作为点估计,因为在初步分析中,手段在大多数情况下比模式更接近于真实的,即模拟的值,这与先前的研究一致[8,38]。遗传性和加性遗传相关性估计的偏差计算为估计值(parest)与真实(即模拟)值(partrue)的平均相对偏差。

bias =(par est-par true)/ par true

使用统计分析系统(SAS),版本9.1.3(SAS Institute,Cary,NC,USA,2005)的程序GLM,通过方差分析测试数据结构和遗传标记信息的质量对ESS和偏差的影响。 。有效样本量或偏倚被认为是因变量,遗传标记信息(r0p9,r1p9,r0p7)的数据集(A1,A2,B1,B2,C1,C2)和质量被认为是固定效应。

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

推荐阅读更多精彩内容