ImageVerifierCode 换一换
格式:PPT , 页数:128 ,大小:3.66MB ,
资源ID:10247666      下载积分:5 金币
验证码下载
登录下载
邮箱/手机:
验证码: 获取验证码
温馨提示:
支付成功后,系统会自动生成账号(用户名为邮箱或者手机号,密码是验证码),方便下次登录下载和查询订单;
特别说明:
请自助下载,系统不会自动发送文件的哦; 如果您已付费,想二次下载,请登录后访问:我的下载记录
支付方式: 支付宝    微信支付   
验证码:   换一换

开通VIP
 

温馨提示:由于个人手机设置不同,如果发现不能下载,请复制以下地址【https://www.zixin.com.cn/docdown/10247666.html】到电脑端继续下载(重复下载【60天内】不扣币)。

已注册用户请登录:
账号:
密码:
验证码:   换一换
  忘记密码?
三方登录: 微信登录   QQ登录  

开通VIP折扣优惠下载文档

            查看会员权益                  [ 下载后找不到文档?]

填表反馈(24小时):  下载求助     关注领币    退款申请

开具发票请登录PC端进行申请。


权利声明

1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前可先查看【教您几个在下载文档中可以更好的避免被坑】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时联系平台进行协调解决,联系【微信客服】、【QQ客服】,若有其他问题请点击或扫码反馈【服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【版权申诉】”,意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4009-655-100;投诉/维权电话:18658249818。

注意事项

本文(QTL定位的原理和方法.ppt)为本站上传会员【a199****6536】主动上传,咨信网仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知咨信网(发送邮件至1219186828@qq.com、拔打电话4009-655-100或【 微信客服】、【 QQ客服】),核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载【60天内】不扣币。 服务填表

QTL定位的原理和方法.ppt

1、单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,*,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,*,第三章,QTL,定位的原理和方法,2025/4/28 周一,1,2025/4/28 周一,2,QTL,是什么?,数量性状位点(,QTL,)是影响数量性状的一个染色体片段;,QTL,定位是确定数量性状基因在染色体上位置的一种方法;,QTL,和,QTLs,。,2025/4/28 周一,3,为什么要定位它?,它为了解个体数量性状基因之间的行为和交互作用等基础知识提供了一条路径,允许建立更加真实的表型变异、选择反

2、应和进化过程模型;,将标记信息综合到遗传评估中,辅助人工选择程序,主要方式有,MAS,和,MAI,;,能进行基因的位置克隆,允许对当前存在的数量变异进行分子机制的研究,并通过直接的分子干预,进一步增加增效等位基因频率。,2025/4/28 周一,4,QTL,定位的基本原则,QTL,定位的基本原则是关联度量的遗传变异和表型变异;,群体的选择、用于度量表型个体选择和基因型判型个体的选择是所有,QTL,定位设计要重点考虑的因素;,对于所有的,QTL,定位设计,标记等位基因和,QTL,等位基因之间的,LD,是必须的。,2025/4/28 周一,5,QTL,定位的关键,2025/4/28 周一,6,20

3、25/4/28 周一,7,2025/4/28 周一,8,第一节,LA,定位(连锁分析定位),2025/4/28 周一,9,linkage analysis,only considers the linkage disequilibrium that exists within families,which can extend for 10s of cM,and is broken down by recombination after only a few generations.,Such as BC and F2 design,2025/4/28 周一,10,单标记分析,2025/4/2

4、8 周一,11,是总平均;和 是加性和显性效应;是标记和,QTL,之间的重组率。,是给定个体标记位点基因型为,Aa,的条件下的,QTL,基因型,Qq,的,条件概率,;,是标记和,QTL,基因型的,联合概率,;,是标记基因型的,边际概率,。,2025/4/28 周一,12,来自近交系的回交群体的标记和,QTL,概率,标记基因型之间的表型值平均差异:,2025/4/28 周一,13,单标记分析的缺点,单标记使用标记平均值,不能获得,QTL,效应单独的估计值和,QTL,与标记的重组频率;因此,不能区分是一个大的,QTL,效应松散地与标记连锁,或是小效应紧密地与标记连锁。,2025/4/28 周一,1

5、4,区间定位,Lander and Botstein(1989),提出使用所有连续的标记进行,QTL,定位的方法;,该方法原则上能够区分,QTL,的效应和位置;,该方法需要一张带有一定数目的遗传图谱,相邻标记间的距离是已知的。,2025/4/28 周一,15,Haldane,作图函数,为遗传距离();,假设减数分裂期间的遗传物质交换沿着染色体是随机和独立发生的。,2025/4/28 周一,16,标记和,QTL,概率,2025/4/28 周一,17,数据分析,为具有,QTL,基因型 的个体 的性状记录;,为具有,QTL,基因型 的个体的期望效应(如 或 );,为随机误差,并且 ,因此有:,202

6、5/4/28 周一,18,最大似然法分析,前面回交例子的似然函数为:,为,QTL,位点的基因型;,和 为个体 在标记位点,A,和,B,的基因型;,为回交个体数。,2025/4/28 周一,19,似然率检验(,LRT,):,为零假设没有分离,QTL,条件下的似然值;,为有一个,QTL,分离条件下的似然值。,LOD,检验:,2025/4/28 周一,20,最小二乘分析,前面回交例子的最小二乘分析模型为:,需要估计的参数:一种为两个,QTL,基因型的平均值;另外一种为总平均值和两个基因型之间的效应差;,显著性检验:,MSQ,为拟合模型由,QTL,基因型解释的方差;,RMS,为拟合模型的残余均方。,2

7、025/4/28 周一,21,LS,和,ML,的比较,LS,只使用了标记平均值信息,标记基因型组内的方差变异没有被使用;而,ML,使用了所有可能的信息,这包括标记基因型和性状分布。,LS,的计算比较简单易行,能够使用标准的软件(,SAS,)进行分析;而,ML,计算非常困难,需要专门的软件将其扩展到非常复杂的模型。,2025/4/28 周一,22,似然率检验和,F,检验的比较:,对一个,QTL,,如果残差呈正态分布,则,LS,和,ML,估计是相同的;,对一般情形,关系变为:,大部分,QTL,定位分析结果显示,LS,获得与,ML,极端近似的结果。,2025/4/28 周一,23,基因组扫描,区间定

8、位的优势在于能对整个标记的基因组进行扫描;,QTL,定位是在整个基因组内进行,某一个区间内,QTL,基因型的条件概率根据侧翼标记信息进行计算,然后一个区间接着一个区间,使用最小二乘或最大似然法进行分析,同时每个区间的检验统计量(,F-ratio,或,LRT,)也被计算,具有最大检验统计量的位置就是,QTL,最可能存在的位置,而该位置的,QTL,效应就是最好的,QTL,估计效应。,2025/4/28 周一,24,2025/4/28 周一,25,多次检测问题,如果有许多独立的零假设被检验,而且事先知道所有的零假设都为真,则,至少出现一次假显著(false positive)的概率为,2025/4/

9、28 周一,26,伯努利校正,2025/4/28 周一,27,Permutation test,对表型和标记基因型数据进行随机重排,它消除了标记基因型和表型之间的关联;,每次重排数据,都要重新在整个基因组中进行,QTL,定位分析;,通过多次重排,可获得每次检验,LRT,统计量在没有,QTL,的零假设条件下的分布;,2025/4/28 周一,28,Permutation test,的具体步骤:,2025/4/28 周一,29,FDR(false discovery rate),is declared FDR(such as 0.05),j is the largest order that me

10、t formula(1),m is the number of marker,2025/4/28 周一,30,FDR(false discovery rate),方法,Sort p values of all marker interval based on ascending order,2025/4/28 周一,31,LOD,下降支撑区间(,LOD drop support interval,),如果某一特定位置检测到一个,QTL,,需要对,QTL,所在的位置执行检验;,零假设是该,QTL,位于估计的峰值位置,备择假设为,QTL,位于距峰值距离为 的位置,,检验统计量为全,QTL,模型在峰

11、值位置和距离峰值位置 图距单位位置的似然函数的差值的两倍,当样本为大样本时,它近似呈自由度为,1,的 分布;,因此可以通过偏离峰值位置,使检验统计量降到一个给定的数值来对,QTL,位置置信区间进行检验。,2025/4/28 周一,32,例如:,95%,的,QTL,置信区间对应的检验统计量下降,3.84,;,1 LOD,下降对应,97%,的,QTL,置信区间;,2 LOD,下降对应,99.8%,的,QTL,置信区间;,2025/4/28 周一,33,2025/4/28 周一,34,Bootstrap,置信区间,对于一个大小为 的群体,抽取 个带有覆盖性质的记录(有些记录被抽取多次,而有些记录没被

12、抽取);,分析并估计,QTL,位置;,重复上面的,1,和,2,两个过程,如,200,次或更多;,在分布的两尾去掉,2.5%,的极端的,QTL,位置估计值;,剩余的,95%,表示置信区间的估计值。,2025/4/28 周一,35,QTL,位置估计的置信区间,2025/4/28 周一,36,预测置信区间,置信区间的长度受样本大小、,QTL,效应和标记密度的影响,对一个高密度标记图谱,,Darvasi and Soller(1997),给出了一个预测的近似,95%,的置信区间(单位,cM),:,为样本大小;,和 为标准的加性和显性效应(以基因型标准差为单位)。,2025/4/28 周一,37,统计能

13、力(,Statistical power,),2025/4/28 周一,38,为什么要计算检测能力?,给定样本大小,计算能够检测到的,QTL,效应;,给定,QTL,效应,估计检测到该,QTL,需要的群体大小;,检测特定的,QTL,时,比较不同的群体设计。,2025/4/28 周一,39,完全连锁标记统计能力的计算理论,型错误():当零假设为真,拒绝零假设所犯错误的概率;,型错误():当零假设为假,接受零假设所犯错误的概率;,统计能力被定义为:,2025/4/28 周一,40,P,(,T,),T,Critical value,H,A,H,0,Statistical errors,2025/4/2

14、8 周一,41,Rejection of H,0,Nonrejection of H,0,H,0,true,H,A,true,Type I error,at rate,Type II error,at rate,Significant result,Nonsignificant result,STATISTICS,R E A L I T Y,POWER,=(1-,),2025/4/28 周一,42,Impact of,alpha,P,(,T,),T,Critical value,2025/4/28 周一,43,Impact of,effect size,N,P,(,T,),T,Critical

15、 value,2025/4/28 周一,44,影响检测能力的重要因素,群体类型;,样本大小;,QTL,效应;,基因组大小;,标记密度;,显著性阈值;,分析类型。,2025/4/28 周一,45,完全连锁标记统计能力的计算,近交系杂交情形下的,QTL,定位检测能力计算基于单标记的,t-,检验和,F-,检验。,2025/4/28 周一,46,F,2,设计:,BC,设计:,2025/4/28 周一,47,对于合理的样本大小和小的,QTL,效应,要求的,t,值为,:,2025/4/28 周一,48,Sample size BC 672 128 42 11 6,2025/4/28 周一,49,BC,和,

16、F2,设计的合理样本大小之比为:,BC,比,F2,的基因组扫描所需的显著性阈值要低;,BC,:,F2,:,BC,比,F2,的 可能要低。,2025/4/28 周一,50,考虑两种设计阈值的变化:,2025/4/28 周一,51,2025/4/28 周一,52,如果连锁不完全(),且使用单标记分析:,如果连锁不完全(),且使用区间定位分析:,2025/4/28 周一,53,为了增加,QTL,检测能力,可以增加判型的个体数目或标记密度;两者之间花费依赖于标记的成本与获得个体表型成本之间的比率。,2025/4/28 周一,54,增加检测能力的方式,增加样本大小;,增加效应大小。,后者可以通过选择一个

17、具有丰富分离,QTL,的群体结构或样本;,如后裔检验。,2025/4/28 周一,55,精细定位,QTL,的群体设计,2025/4/28 周一,56,Fine Mapping Strategies,Genomewide-based strategies:,Large scale BC,F2,half sibs,etc.,Recombinant inbred lines(RIL),Advanced Intercross Lines(AIL),Locus-based strategies:,Selective phenotyping,Recombinant progeny testing,Inte

18、rval specific congenic strains(ISCS),Recombinant inbred segregation test(RIST),2025/4/28 周一,57,重组近交系(,Recombinant inbred lines,RIL,),重组近交系来源于,F2,群体的近交;,RIL,只需要被判型一次,却能很好地度量多个性状(,clonal Lines,);,RIL,关键的特性是比,F2,发生更多的重组,数量性状通过使用系平均值能被准确度量;,RIL,只能定位加性,QTL,;,RIL,的产生慢而困难。,2025/4/28 周一,58,深度杂交系(,Advanced i

19、ntercross lines AIL,),AIL,开始于,F2,群体,杂交后裔继续杂交一定数目的世代(与,RIL,近似,但是远交,而不是近交);,AIL,是在,F2,群体,QTL,定位的基础上进一步提高,QTL,的定位精度;,AIL,的任何性状都能被度量,但基因型判型只着眼于感兴趣的区域;,AIL,的关键特性是在目标区域创造了附加的重组事件,类似于扩大了,F2,群体。,2025/4/28 周一,59,Advanced intercross lines(AIL),Semi-random intercrossing,P,F,1,F,2,F,3,F,t,CI=CI,F2,/(t/2),2025/4

20、/28 周一,60,AIL,要保持一定的群体大小;,AIL,相对于,F2,使重组近似增加 ,置信区间为:,AIL,能定位几个,QTL,或多个,QTL,到,15cM,。,2025/4/28 周一,61,Locus-based strategies:,Selective phenotyping(SPh),Theoretical basis:,Only recombinants increase mapping accuracy for a detected QTL.,Procedure:,A large F2 or BC population,Only individuals recombinan

21、t at a QTL-containing interval are subsequently phenotyped.,2025/4/28 周一,62,Requires only 2 generations,.,Requires very large samples,.,2025/4/28 周一,63,SPh-Experimental results,Lesions density,Paigen et al.,BC,SPh-BC,2025/4/28 周一,64,重组后裔检验(,Recombinant progeny testing,),2025/4/28 周一,65,Recombinant p

22、rogeny testing,QTL,Males,recombinant at an interval of interest,are progeny,tested to check which QTL allele was retained.,Requires only 3 generations.Efficient for dominant effects,Requires large sample,2025/4/28 周一,66,Interval-specific congenic strains,区间特异,同源异基因品系,2025/4/28 周一,67,Interval specifi

23、c congenic strains(ISCS),QTL,ISCS,are produced by a series of backcrosses and intercrosses,-,Requires very few individuals.Useful fro further studies-Complicated and lengthy development process.,2025/4/28 周一,68,Recombinant inbred segregation test(RIST),P,1,RI P,2,x,x,F,1,1,F,1,2,F,2,1,F,2,2,QTL,2025

24、/4/28 周一,69,Each selected RIL is backcrossed to each parent and then the BC1 is selfed and grown out for phenotypiing and genotyping in the QTL region.Because the QTL was previously mapped to this region,the BC to one of the parents will segregate while the other will not;thus,indicating whether the

25、 gene controlling the QTL is above or below the breakpoint.The overlapping results of the various RILs will narrow the QTL interval.,-,Requires only 2 generations.Few individuals required;-Requires RILs with recombinations in region of interest.,2025/4/28 周一,70,RIST-Experimental results,F,21,F,22,C5

26、7L,AKR,AKXL-16,P=0.41,D2MIT64,D2MIT200,P=0.02,B.Taylor,A.Darvasi,Obesity QTL,2025/4/28 周一,71,第二节 特定结构的远交群体,QTL,定位,2025/4/28 周一,72,近交和远交群体的差别,远交群体也存在部分的近交;,远交群体的主要特征是群体内部没有故意尽力让亲属之间进行配种而创造近交(随机交配);,远交群体与近交群体的主要差别是远交群体内有遗传变异正在分离;,远交群体,QTL,与标记的关联是特定家系的关联,而不是群体范围的关联;,远交群体家系间存在附加的遗传方差。,2025/4/28 周一,73,使用

27、远交群体进行,QTL,定位的基本策略,在存在差异的远交群体之间寻找,QTL,;,在一个群体内寻找正在分离的,QTL,;,具体策略为:,使用遗传标记追踪从父母亲到后裔的遗传,在基因组的所有位置获得不同基因型可能的概率;,在家系内关联表型数据和基因型概率数据。,2025/4/28 周一,74,远交系杂交,2025/4/28 周一,75,远交系杂交与,F2,情形比较类似,但现在是两个经济性状存在差异的远交系或品种进行异型杂交;,所有三代,包括祖代、,F1,和,F2,在多个标记位点都要进行标记基因型判型,但只有,F2,个体获得表型;,F2,个体的,QTL,基因型(,QQ,、,Qq,、,qQ,和,qq,

28、的两个等位基因不能区分那一个来自父亲,那一个来自母亲;,在远交群体杂交,要考虑加性效应、显性效应和父母亲来源效应(印记效应)。,2025/4/28 周一,76,半同胞群体,2025/4/28 周一,77,家畜群体存在很大的半同胞家系;,半同胞家系,QTL,定位的原则是关联半同胞后裔的表型和它们遗传自共同祖先的等位基因的概率;,父母亲和半同胞后裔都要判定基因型,只有后裔度量表型。,2025/4/28 周一,78,单个半同胞家系,单标记,对于单个半同胞家系,唯一的要求就是共同祖先在一个标记位点是杂合子;,因此,能看到后代在一个位点的两个等位基因的表型值平均差异;,该差异能使用,t,检验进行显著性

29、检验;,单标记的单个半同胞家系类似于,BC,设计。,2025/4/28 周一,79,多标记,多标记条件下,共同父母亲标记之间的连锁相未知;,需要重新构建父母亲的“单倍型”;,根据父母亲的配子和后代的基因型,获得最可能的父母亲“单倍型”;,参考,BC,的计算方法,在父母亲“单倍型”已知的条件下计算每个,HS,后裔遗传自父母亲某一个配子的条件概率;,表型对估计的条件概率进行回归获得,QTL,等位基因之间的差异,利用,t,检验进行差异显著性检验。,2025/4/28 周一,80,公畜单倍型重构,确定每个,HS,后代的信息标记,即确定公共父母亲的那些标记是杂合的,且等位基因的传递是清楚的;,考虑某一个

30、公共父母亲的那一些相邻标记对是杂合的;,计算出两个相邻位点等位基因能确定遗传自公共父母亲的后代数目;,利用期望最大法(,EM,)在最小化重组数目的基础上建立标记位点的连锁相,2025/4/28 周一,81,多个半同胞家系,对上面四个父母亲的后代进行标记对比分析将不会检测到,QTL,,因为,M,和,m,的差异为零,所以应该考虑家系内的嵌套分析。,2025/4/28 周一,82,单标记,对于单标记的多个家系,可以使用嵌套的,ANOVA,,考虑嵌套在家系内的标记效应:,2025/4/28 周一,83,孙女设计(,Granddaughter design,GDD,),Weller et al.(199

31、0),介绍了在半同胞家系中一个孙女设计被应用来进行,QTL,定位;,该设计要求公畜、儿子和女儿的三代系谱,公畜和儿子被判定基因型,孙女获得表型;,利用在某一位点遗传了公畜两个可选等位基因的儿子的女儿的表型平均值来比较定位,QTL,。,GDD,的优势在于获得相同检测能力的条件下比较少的个体需要被判定基因型;,GDD,比较容易收集数据,因为公牛的,AI,体系。,2025/4/28 周一,84,该分析通常采用儿子的女儿离差(,daughter yield deviations DYD,)来进行;,因此能使用女儿设计模型应用,ANOVA,和回归进行分析;,儿子女儿数目如果变化很大,这时需要对,DYD,

32、进行加权。,2025/4/28 周一,85,NCP for the daughter design as:,NCP for the granddaughter design as:,Once the NCP parameters is calculated,power is derived as the probability that a non-central variate exceeds the threshold from a central distribution.,GDD is generally much more powerful than a daughter desig

33、n,2025/4/28 周一,86,全同胞家系,单个或多个大的全同胞家系在绝大部分物种内都是不可能的,但检测到,QTL,的能力很强。可能的原因有;,全同胞家系存在两个标记差异,一个是父亲,另一个是母亲;,全同胞相对于半同胞,期望的标记差异包含加性和显性方差;,2025/4/28 周一,87,同胞对和核心家系(配对设计),大部分物种都不大可能获得大的全同胞或半同胞家系;,怎样在那样的群体内定位,QTL,呢?,一种设计是收集没有亲缘关系的同胞对或核心家系;,这时要将,QTL,效应作为随机效应,在同胞对之间关联类似的表型和它们类似的等位基因。,2025/4/28 周一,88,2025/4/28 周一

34、89,第三节 方差组分,QTL,定位,2025/4/28 周一,90,模型和检验统计量,An example of a linear mixed model for a single QTL analysis is,:,2025/4/28 周一,91,Variance components can be estimated using maximum likelihood or restricted maximum likelihood(REML),The log-likelihood function is:,The assumed mean and variance structure

35、of the observations:,Q,is the IBD matrix:,2025/4/28 周一,92,The distribution of the test statistics are,asymptotically,a mixture of zero(with probability)and a with 1 degree of freedom(also with probability of).,2025/4/28 周一,93,The advantage of this likelihood-based approach.,The full maximum likeliho

36、od approach simultaneously estimates the IBD probabilities and the variance components,in a combined segregation analysis and linkage analysis framework.,“distribution method”,“expectation method”,2025/4/28 周一,94,So why is QTL mapping in general pedigrees not used more frequently,in particular in la

37、rge,deep pedigrees?,IBD estimation in large pedigrees.,the unavailability of(user-friendly)software for the variance component estimation part of the analysis.,a finite budget.,the unavailability of DNA samples from most ancestors,2025/4/28 周一,95,IBD,估计,2025/4/28 周一,96,Perfect marker,As in the case

38、of sibpairs,IBD sharing using a fully informative marker is straightforward,because we can simply count the number of alleles that two relatives share by descent.,At a location linked to a perfect marker,IBD probabilities can be calculated from the observed,IBD probability at the marker,the average

39、relationship between individuals,and,the recombination rate between the marker and putative QTL position,.,2025/4/28 周一,97,The general case:missing data and non-informative markers,The marker information in complex pedigrees is often incomplete.,Unknown linkage phases,non-informative markers and/or

40、missing marker genotypes complicate the calculation of,Q,.,The calculation methods of,Q,are:,recursive algorithms,correlation based algorithms,simulation based algorithms.,2025/4/28 周一,98,Implementation in Loki,The multiple-site segregation sampler in Loki is a cleverly designed Gibbs sampler with b

41、atch updating.,is the probability of the segregation indicators across,n,loci at the,i,th segregation conditional on all other segregation indicators and observed marker data.,2025/4/28 周一,99,A two step strategy to sample,The first step involves moving through the genome,calculating locus by locus,c

42、umulative probabilities for S,ij,.,the second step involves moving back down the genome,sampling S,ij,from a univariate density that is a function of the associated cumulative probability,the previous sampled segregation indicator(S,i j,+1,)and the recombination rate between loci,j,and,j,+1.,2025/4/

43、28 周一,100,Introduction to Loki,Loki was originally designed for multipoint linkage analysis in general pedigrees using MCMC methods.,Then,it has since been modified for IBD probability calculation.,The user supplies Loki with the pedigree structure,marker genotypes,marker positions and QTL positions

44、 for which the IBD matrices are to be calculated.,Dependent chains of IBD probabilities are then obtained for each QTL position.,Convergence is determined by monitoring the IBD probabilities over the iteration number.,Once the probabilities stabilize,the sampler is deemed to have reached convergence

45、2025/4/28 周一,101,Variance component estimation,After having calculated IBD probabilities,there are two difficulties in estimating variance components by ML(REML).,Firstly,the IBD matrix is a completely general symmetrical matrix and does not have an obvious inverse.,Secondly,the IBD matrix is like

46、ly to be singular.,2025/4/28 周一,102,why the IBD matrices are often singular?,The reason is that two related relatives can share 0 or 100%of their alleles IBD,which can cause a dependency in the matrix of IBD probabilities.,The genotypes of the parents are M1M2 and M3M4.If the progeny have genotypes

47、M1M3 and M2M4(a),or M1M3 and M1M3(b),then the resulting IBD matrix is:,a,b,2025/4/28 周一,103,If the maximisation algorithm is based upon the complete matrix,V,(or,V,-1,),then there should not be a problem.,If the maximisation is based upon an algorithm that requires,Q,-1,then using genomic positions

48、which are slightly distant from the markers will give a positive-definite,Q,2025/4/28 周一,104,Implementation example,Visscher et al.(1999)used the combination of an MCMC sampling approach and REML variance component estimation to map a QTL for bipolar disorder(manic depression)in a human pedigree.,Th

49、e pedigree size was 168,over 4 generations,and 143 individuals had a phenotypic score.,The incidence of major recurrent depression(unipolar disorder)and bipolar disorder was 17/143 and 11/143.,A small segment of chromosome 4 was considered because this region had previously shown linkage to bipolar

50、disorder using a parametric linkage analysis,and 11 microsatellite markers were scored spanning 26 cM.,2025/4/28 周一,105,IBD probabilities were estimated using Loki,using 10,000 samples.,REML was used to estimate 81 variance components,with an algorithm based upon the complete(co)variance matrix,V,to

移动网页_全站_页脚广告1

关于我们      便捷服务       自信AI       AI导航        抽奖活动

©2010-2025 宁波自信网络信息技术有限公司  版权所有

客服电话:4009-655-100  投诉/维权电话:18658249818

gongan.png浙公网安备33021202000488号   

icp.png浙ICP备2021020529号-1  |  浙B2-20240490  

关注我们 :gzh.png    weibo.png    LOFTER.png 

客服