人类群体遗传学—第一讲

1、人类群体遗传学基本原理和分析方法,中科院-马普学会计算生物学伙伴研究所,中国科学院上海生命科学研究院研究生课程 人类群体遗传学,徐书华 金 力,20082009学年第二学期人类群体遗传学分析方法课程表 上课时间：每周四上午10:00-11:50 上课地点：中科大厦4楼403室第7教室,第一讲 Hardy-Weinberg平衡检验及其应用,Hardy-Weinberg Equilibrium (HWE)定律 HWE的偏离（固定指数） HWE的偏离（ Wahlund定律） HWE在单基因遗传疾病研究中的应用 HWE在多基因复杂疾病研究中的应用 检验HWE的常用方法和软件 全基因组数据中的HWE检验

2、问题,Why do students need to know Hardy-Weinberg Equilibrium？,As Thomas Merten (1992) states: If you have ever been asked questions such as the ones that follow, you begin to see why studying population genetics might be useful: 1. Im confused! How can O be the most common of the blood types if it is

3、a recessive trait? 2. If Huntingtons disease is a dominant trait, shouldnt three-fourths of the population have Huntingtons while one-fourth have the normal phenotype? 3. Shouldnt recessive traits be gradually swamped out so they disappear from the population?,基本概念,突变（mutation） 遗传多态性（genetic polymor

4、phism） 遗传座位（locus） 等位基因（allele） 核苷酸位点（nucleotide） 序列（sequence） 位点（site） 分离位点（segregating site）,基本概念,单倍体（haploid）、二倍体（diploid） 基因型（genotype） 杂合子（heterozygote）、纯合子（homozygotye） 等位基因频率（allele frequency） 基因型频率（genotype frequency）,一个群体由N个个体组成，其中有一对常染色体等位基因A, a，其频率分别以p、q表示 。 可能的基因型为AA，Aa，aa三种，其频率分别以D、H、R表

5、示，其中D+H+R=1。,若群体中有,n1 AA n2 Aa n3 aa N=n1+n2+n3,合计,1,计算基因型频率,等位基因A的频率,等位基因a的频率,其中，p+q=D+H/2+H/2+R=1,计算等位基因频率,精 子 A（p） a（q） 卵 A（p） AA（p2） Aa（pq） 子 a（q） Aa（pq） aa（q2） 由表可见子代基因型组成：p2+2pq+q2=1 这里基因型AA的频率为p2，基因型aa的频率为 q2，基因型Aa的频率为2pq。 AA：Aa：aa= p2：2pq：q2,子一代向下一代提供的配子中两种基因频率分别是： A=p2+1/2（2pq）=p2+pq=p（p+q）

6、=p a=q2+ 1/2（2pq）=q2+pq=q（p+q）=q 由此可见，子代基因A的频率仍然是p，基因a的频率仍然是q，而且将以这种频率在所有世代传递下去，这就是遗传平衡。,p2 2pq q2,Hardy-Weinberg平衡定律,Hardy-Weinberg平衡定律,群体中的等位基因频率以及基因型频率并不随世代的推移而变化。,保持群体HWE的理论条件,Random mating No differential fertility of the genotypes Equal genotype frequencies in the two sexes No mutation No immi

7、gration No differential emigration No differential viability Infinite population size,人类进化中的5个基本动力元素,突变（mutation） 婚配 （random mating） 选择 （selection） 漂变 （drift） 迁移 （migration）,现代人群是一个HWE群体吗？,有趣的是，对于一个特定的地理群体，绝大部分基因组区域处于HWE状态。,HWE平衡群体的性质（1）,杂合体的比例为H=2pq，H的最大值为1/2,当q=1/2，p=1/2时，H取最大值1/2,HWE平衡群体的性质（2）,Aa

8、Aa的交配频率永远为AAaa交配频率的2倍 AaAa=H2=4p2q2 AAaa=2DR=2p2q2,因此，隐性基因“a”在群体中不易消失，总是保持一定频率。当一般群体中隐性个体稀少的时候，其中大多数为AaAa交配的后代。,例如 人类白化病被认为是一个简单的隐性性状，某地白化病的婴儿的发病率q2 (aa)=1/10000 , q(a)=1/100 但群体中，杂合子的比例要大得多,其中由AaAa得到的aa个体占,当q很小时，p=1，则H-Ws law取得一种极限形式,H=2q 差不多所有的隐性基因都处于杂合状态中，杂合子个体的比例约为隐性基因频率的两倍。上述白化病 q=0.01,一般q2可看作发

9、病率，杂合子的频率=2,HWE平衡群体的性质（3）,HWE在单基因遗传疾病研究中的应用,计算基因型频率 随机交配时纯合子（基因型）的频率可用配子（基因）频率平方求得 ； 杂合子频率=2x隐性等位基因频率（如前所述）. 计算隐性性状的发病率（出生前死亡） 计算等位基因频率,计算等位基因频率（隐性）,如已知某种常染色体隐性遗传病（白化病）在一特定人群中频率，就能计算这个异常基因的携带者和基因频率。白化病发病（aa）(q2)的频率为110000，即其基因型频率，则致病基因（a）频率1/100= 0.01；基因A的频率10.01=0.99；而杂合子携带者的频率为2pq=2991001100150。因此

10、，在比例中，每个受累的个体将有200个左右在临床上无症状的携带者。,计算等位基因频率（显性）,常染色体显性遗传病，如并指症，在一个群体多为杂合子（Aa）发病。 杂合子（H）的频率为2pq，由于q值大，近于1，故H2p,p=1/2H。因此，只要知道杂合子发病率，就易求得基因A的频率。如并指症的发病率为 11000，H12000，p=1/2H=1/4000,即致病基因A的频率为0.000025。,在已知某一性状频率的群体中，确定等位基因频率和杂合子携带者频率。,例如：已知白化病的发病率为1/20000，求白化病致病基因频率q和携带者频率。 白化病为AR遗传病，患者为致病基因的纯合子，因此： 发病率

11、（aa） = q2 = 1/20000 a= q =0.007 A=p = 1 q =0.993 携带者频率(Aa) = 2pq = 2 0.993 0.007 = 0.0139 这提示人群中有1.4%为白化病致病基因携带者，对于遗传咨询很重要。,对于一种罕见的AD遗传病 致病基因频率p很低，p2可以忽略，因此： 2pq p2 + 2pq 患者几乎都以杂合子形式出现,=1, 对于一种罕见的XD遗传病 致病基因频率p很低，正常基因频率q1 男性患者 = 致病基因频率 = p 女性患者 = p2 + 2pq = 2p 男性患者 女性患者 女性患病率是男性患病率的2倍,= 1/2, 对于一种罕见的A

12、R遗传病 致病基因频率q很低，正常基因频率p1 ： 杂合携带者频率 = 2pq =2q 即杂合携带者频率约为致病基因频率的2倍。 杂合携带者 2pq 2 纯合患者 q2 q 提示: 致病基因频率q越低，致病基因在群体中主要以杂合携带者形式存在。,=,=, 对于一种罕见的XR遗传病 致病基因频率q很低， 男性患者 = 致病基因频率 = q 女性患者 = 致病基因频率的平方 = q2 男性患者 q 1 女性患者 q2 q 提示：疾病越罕见，q越低，女性患者越罕见。,=,=,HWE的偏离（固定指数）,固定指数（F）: 如果一个座位上有两个等位基因，Hardy-Weinberg比率的任何偏差可以由参量

13、F来度量，F称为固定指数，则基因型频率可以由下式给出：,由以上第二式可得：,固定指数F可正可负，视情况而定。 可以看出，当h0小于h时，F取正值；当h0大于h时，F取负值。在近亲交配时，杂合子频率的观察值减小，F就取正值。,随机交配（h）情况下杂合子的预期频率,群体（h0）中下杂合子的观察频率,上式可写成,亚群体（sub-population）,以上考虑的是一个简单的群体，不论其是否近亲交配。 然而，实际上大多数的自然群体可被再分为许多不同的繁殖单位或亚群体（sub-population），尽管这些群体并不是完全隔离的。这种情况下，研究群体内和群体间的遗传变异就显得十分重要。,可再分群体中的基

14、因型频率,假定一个群体可分为s个亚群体，每一个亚群体都满足Hardy-Weiberg平衡。设xk为第k个亚群体中等位基因A1的频率，则基因型A1A1，A1A2，A2A2的频率分别为,我们用wk来表示第k个亚群体的相对大小，且总和为1。则A1A1，A1A2，A2A2在整个群体中的频率为：,其中,和,是亚群体中等位基因频率的均值和方差。,可再分群体中的固定指数,比较,我们知道,，因此,表明如果一个群体被分为多个交配单位，纯合子的频率要高于Hardy-Weinberg比率。这个性质首先由Wahlund（1928）发现，被称为Wahlund定律，也称Wahlund现象。 当等位基因频率在所有亚群体中一

15、致时，F为0；而当每个亚群体都被固定为某一个等位基因时，F为1。,HWE的偏离（Wahlund定律）,Wahlund现象的启示,群体结构（population structure）的存在！ 反之，当F为负值的时候，,杂合子频率比Hardy-Weinberg平衡时预期的要高，意味着杂合优势，某种程度的自然选择发生。,杂合优势与平衡选择（后面“自然选择”章节细谈）,检验HWE的常用方法,X2-test,观察值,期望值,实例,例 MN血型基因频率的计算,结论：三个基因型频率的观察值与HWE定律所期望的频率没有显著差异。,精确检验,Am.J.Hum.Genet 2005, 76:887-893,HWE

16、在复杂疾病遗传研究中的应用,数据分析中，检测可能的数据质量问题（typing error）。 Case-control study 一般要求所有检测位点在control中处于HWE。 多个位点处于HWD暗示近亲婚配或群体结构的存在。 在排除其他可能的因素以后，筛选可能受到自然选择的基因。,全基因组数据中的HWE检验问题,Affymetrix 500K, Illumina 650K 多重检验（multiple testing）问题 If you have 10,000 genes in your genome, and perform a statistical analysis, a p-va

17、lue cutoff of 0.05 allows a 5% chance of error. That means that 500 genes out of 10,000 could be found to be significant by chance alone.,多重检验校正方法,Bonferroni correction Bonferroni Step-down (Holm) correction Westfall and Young Permutation Benjamini and Hochberg False Discovery Rate,Bonferroni correc

18、tion,The p-value of each gene is multiplied by the number of genes in the gene list. If the corrected p-value is still below the error rate, the gene will be significant: Corrected P-value= p-value * n (number of genes in test) 0.05 As a consequence, if testing 1000 genes at a time, the highest acce

19、pted individual p-value is 0.00005, making the correction very stringent. The expected number of false positives will be 0.05.,Bonferroni Step-down (Holm) correction,This correction is very similar to the Bonferroni, but a little less stringent: 1) The p-value of each gene is ranked from the smalles

20、t to the largest. 2) The first p-value is multiplied by the number of genes present in the gene list; if the end value is less than 0.05, the gene is significant; Corrected P-value= p-value * n 0.05 3) The second p-value is multiplied by the number of genes less 1. Corrected P-value= p-value * n-1 0

21、.05 4) The third p-value is multiplied by the number of genes less 2. Corrected P-value= p-value * n-2 0.05 It follows that sequence until no gene is found to be significant.,Westfall and Young Permutation,The Westfall and Young permutation follows a step-down procedure similar to the Holm method, c

22、ombined with a bootstrapping method to compute the p-value distribution: 1) P-values are calculated for each gene based on the original data set and ranked. 2) The permutation method creates a pseudo-data set by dividing the data into artificial treatment and control groups. 3) P-values for all gene

23、s are computed on the pseudo-data set. 4) The successive minima of the new p-values are retained and compared to the original ones. 5) This process is repeated a large number of times, and the roportion of resampled data sets where the minimum pseudo-p-value is less than the original p-value is the adjusted p-value. Because of the permutations, the method is very slow.,Benjamini and Hochberg False Discovery Rate,This correction is the least stringent of all 4 options, and therefo