利用连锁不平衡进行畜禽QTL定位的研究进展
- 格式:pdf
- 大小:317.88 KB
- 文档页数:4
利用连锁不平衡进行畜禽QT L定位的研究进展谢水华1,李加琪1,张豪1,陈瑶生2(1.华南农业大学动物科学学院,广州 510642;2.中山大学生命科学学院,广州 510275)摘要:一般家养畜禽群体规模小、选择强度高,连锁不平衡(linkage disequilibrium,LD)现象普遍存在,尤其在杂种或新培育品种群体中更为突出。
随着分子遗传学的发展,我们可以直接选择影响性状的基因或包含基因的染色体区域(Q TL),可分型的遗传标记和Q TL之间是否存在连锁不平衡是畜禽群体中Q TL定位和应用的关键。
目前,畜禽群体中用基因组范围内SN Ps的连锁不平衡精细定位影响畜禽重要经济性状的Q TL正受到广泛关注。
作者从LD定位Q TL的原理、畜禽群体中LD程度、利用LD进行Q TL定位的试验设计及统计分析方法、LD在畜禽Q TL精细定位中的应用等几个方面进行了综述。
并进一步阐述了畜禽群体中利用LD进行标记辅助选择(MAS)的策略,对LD在畜禽Q TL定位及标记辅助选择方面的应用进行了新的探讨。
关键词:连锁不平衡;Q TL定位;关联分析;标记辅助选择中图分类号:Q348 文献标识码:A 文章编号:167127236(2007)0620036204 连锁不平衡(linkage disequilibrium,LD)又称等位基因关联(allelic association),它是相邻基因座位上等位基因的非随机性相关。
在多基因疾病中若应用连锁分析定位微效基因,其所需要的家系数目将大得惊人,故提出了关联分析(test of associa2 tion)的方法,其中基于LD的关联分析在定位复杂疾病基因上显示出强大的功能(Risch等,1996)。
LD分析在人类基因精细定位、群体遗传学参数估计及人类的进化和迁移研究中已经得到了广泛应用。
近年来随着家鸡、牛、猪全基因组高精度序列图的完成,在畜禽群体中利用LD分析定位重要经济性状的Q TL已经引起了畜牧工作者的高度关注,其中用基因组范围内SN Ps的LD来精细定位重要经济性状的Q TL更是受到广泛关注。
可分型的遗传标记和Q TL之间是否存在LD是畜禽群体中Q TL定位和应用的关键,其LD水平直接决定关联分析的统计效力。
1 LD定位QT L的原理群体中产生LD的主要因素有突变、选择、漂变、迁移和杂交,LD的概念可用一个简单的等式来描述和衡量:假设存在相邻基因座位1和2,座位1的等位基因为A、a,其频率为P A、P a;座位2的等位收稿日期:2006212208作者简介:谢水华(1984-),男,湖南人,硕士生,研究方向:分子数量遗传与动物育种。
通讯作者:李加琪,教授,博士生导师。
E2mail:jqli@.cn,电话:020*********基因为B、b,其频率为P B、P b。
等位基因A和B在同一染色体上同时出现的频率,即它们组成的单倍型频率为P AB,则LD值D=P AB-P A×P B。
如果位点2上的等位基因B与畜禽优良经济性状有关,那么等位基因A的频率在经济性状较好的群体中将高于对照群体。
换句话说,等位基因A与优良经济性状相关,事实上可以通过检测遍布基因组中的大量遗传标记位点或者候选基因附近的遗传标记来寻找与影响性状的Q TL距离足够近而表现出与性状相关的标记位点,这就是LD定位Q TL的基本原理。
利用LD分析定位Q TL时,依赖于是否能发现标记座位的基因型和表型之间的关联,通过比较不同基因型个体表型平均值的差异,可以发现是否存在这种关联,如果不同基因型个体间的表型值存在差异,则说明该标记和一个Q TL连锁。
在畜禽杂交群体中群体范围内的LD广泛存在,这些广泛存在的LD可以用来检测与标记有一定距离的Q TL,但是不能用来精确定位Q TL;在远交群中家系内LD也是广泛存在的,我们可以用覆盖整个基因组有限的与Q TL有一定距离的显著标记进行Q TL 定位,远交群中群体范围内有足够多的LD可用于遗传改良,除了紧密连锁的基因座以外,经过多个世代闭锁的群体一般处于连锁平衡(linkage equilibri2 um,L E)状态,在这些群体中,只有与Q TL紧密连锁的标记才有可能与表型相关联,即使如此,也由于存在随机抽样效应,不能保证一定会存在关联(J ack,2005)。
2 LD值的估算2.1 双等位基因LD值的估算 对于LD的度量已有多种不同的方法,其中大多数都是用于双等位基因的配对检验。
目前常用的2种配对检验方法为LD系数(coefficient of linkage disequilibrium,D′)和r2(Pritchard等,2001)。
当D≥0时,D′=D/min (P A P b,P a P B);当D≤0时,D′=D/min(P A P B, P a P b)。
r2的计算公式为:r2=D2/(P A P a P B P b)。
这2个参数的取值范围都在0~1之间,D′等于1时说明两位点等位基因没有分离。
r2的数值表示一个位点可反映另一位点信息量的程度,在某种程度上可看作D′的补充,r2等于1称为完全连锁不平衡(perfect LD),这时只观察一个标记即可提供另一标记的全部信息。
当Q TL和一个遗传标记间的连锁不平衡用r2度量时,要达到与直接检测这个Q TL相似的显著性水平,样本容量需达到原来的1/r2倍(Risch等,1996)。
2.2 多等位基因LD值的估算 在研究群体遗传标记LD时,标准的χ2检验比r2和D′更有效,r2在多等位基因时会低估LD值,D′在广泛的分离和非同线性标记配对时会给予高的LD值(Heifetz等, 2005)。
利用多等位基因标记LD进行Q TL定位时,标准化的χ2检验是一个非常不错的方法(Zhao 等,2005)。
Farnir等(2000)、McRae等(2002)、Vallejo 等(2003)定义了D′和χ2‘:D′=∑i ∑jP(A i)P(B j)|D ij/D max ij|χ2′=χ2/[2N(n-1)]D ij=P(A i B j)-P(A i)P(B j)D ij<0时,D max ij=min[P(A i)P(B j),(1-P (A i))(1-P(B j))]D ij>0时,D max ij=min[P(A i)(1-P(B j)),(1-P(A i))P(B j)]χ2=N∑i ∑jD2ij/[P(A i)P(B j)]其中P(A i)为座位A上等位基因i的频率,P (B j)为座位B上等位基因j的频率,N为群体大小,n 为标记位点具有最少等位基因的等位基因个数。
3 利用LD进行QT L定位的试验设计和统计分析畜禽群体处于LD状态是LD分析定位畜禽Q TL的前提和基础,其定位方法从标记应用的数量可以分为单标记、两标记和多标记等,并延伸出所谓的区间定位(interval mapping)。
利用LD精细定位Q TL主效基因时,影响其成功的因素有:试验设计、表型测定的精确性、基因效应大小、遗传杂合度、LD 水平以及适当的统计方法。
3.1 试验设计 利用LD进行Q TL定位的试验群体一般包括近交系杂交群(crosses between inbred lines)和分离群体(segregating pop ulations)。
当使用近交系杂交群时,通常采用回交设计(BC)、F2设计、测交设计(TC)和重组近交系设计(recombinant inbred lines);当使用分离群体时,通常采用全同胞(f ull sibs)设计、半同胞(half sibs)设计、祖父2孙女设计(grand2daughter design)、动物模型(animal model),对于小的全同胞家系一般采用同胞配对分析(sib2pair analysis),大的全同胞家系采用全同胞设计,大的半同胞家系采用半同胞和祖父2孙女设计(Joel,2001)。
当使用不同的试验设计进行Q TL定位时,区间定位要优于单标记定位;有显性效应存在时,单亲回交的检测效率比F2设计要高,无显性效应存在时F2代设计比回交设计效率高(Soller等, 1976),因为与单个亲本回交只能检测出杂合子效应,在F2代检测出的Q TL加性效应是纯合子效应的一半,为了估计纯合子效应,需要与双亲进行回交,但他的效率也比F2设计低。
3.2 统计分析 利用t检验或方差分析来定位与标记连锁的Q TL基本上是可靠的,除非严重偏离正态分布,如由少量的、效应大的Q TL明显区分的2个亲本系。
最大似然估计法(maximum2likelihood estimation)考虑了Q TL分离的标记基因型准确分布特性,可以提供更适合的参数估计和显著性检验,因此在利用LD定位Q TL时使用最大似然估计更为理想(盛志廉等,1999)。
Meuwissen等(2000, 2002)对LD定位Q TL时用最大似然法估计Q TL 的位置进行了大量研究,接下来我们将对Meuwis2 sen在这方面的报道进行简单介绍。
畜禽群体中用最大似然法估计Q TL位置时,首先必须对最后一代的表型、标记数据及父母代的标记数据进行记录,分清哪些标记等位基因来自父亲哪些来自母亲,弄清楚标记等位基因的连锁相后就可以构建标记单倍型,例如:个体基因型为M1M2/N1N2,父母分别为M1M3/N1N3、M2M4/ N2N4,则个体的连锁相为M1N1/M2N2。
如果不知道标记等位基因连锁相,将无法得到单倍型信息,因此在分析时这种数据不能用。
计算最后一代表型记录的模型为y=X b+Zh+e(y为记录值,b为固定效应,h为单倍型随机效应,e为残差,X、Z分别为b、h的概率矩阵)。
Var(e)=σe2R(R为单位矩阵);Var(h)=σh2H p(H p为单倍型效应协方差矩阵)。
接着通过下面的模型可以计算得到假设位点存在Q TL的概率:L(H p,σ2h,σ2e)∝-0.5ln|V|+ln|X′V-1X|+y-X^b′V-1y-X^b 式中V=Var(y)=Z H p Z‘σh2+σe2R,^b为b的估计值,从公式中不难看出要计算出Q TL位置概率就必须获得单倍型效应方差矩阵H p。
H p为单倍型效应协方差矩阵,协方差在两个单倍型效应h i、h j 的计算公式为:CoV(h i,h j)=Prob(IBD/marker haplotypes)×σh2 Pro b(IBD/marker haplotypes)为在给定标记单倍型时Q TL座位IBD概率。
通过上面的方法同样可以算出假设位点不存在Q TL时的最大似然值。
用遗传距离作为横坐标,似然值的常用对数作为纵坐标,在坐标轴上就可清晰的看到Q TL峰,而Q TL峰所对应的横坐标就是Q TL在染色体上的位置。