X染色体印记效应检验及合并X染色体失活与印记效应的关联分析研究

摘要

背景:
　　关联分析(association analysis)是用来检验等位基因频率或基因型频率在两组个体（患病人群和正常人群）之间是否存在差异的方法。其主要目的是为了定位疾病位点，利用疾病位点和标记位点的关系来找到遗传变异。如果标记位点和疾病位点之间的距离很接近，它们可能存在连锁不平衡(linkage disequilibrium)，那么这两个位点很有可能是共同遗传给后代的。如果在患病的人群中这个标记位点的频率比正常人群的高，利用关联分析方法就可能将该疾病位点进行定位。
　　遗传印记效应(genomic imprinting)是熟知的一种病因学因素。遗传印记基因的表达取决于等位基因是来自父亲还是来自母亲，而这种现象可能会导致疾病和基因之间的关联。遗传印记现象会影响后代在子宫内的发育以及出生后个体的成长发育。与常染色体上变异的遗传印记作用有关的疾病有:Beckwith-Wiedemann综合征、Prader-Will综合征、Angelman综合征、糖尿病和精神分裂症等。另外，X染色体上的印记基因在复杂疾病和性状的形成上也许起到了关键的作用，如Turner综合征和孤独症等。目前已经有很多方法可以用来检验基于核心家庭和一般家系常染色体上的遗传印记效应。有关X染色体上印记效应的检验，仅仅有最近新提出的基于核心家庭数据的X染色体质量性状位点上的亲代不对称检验方法(parental-asymmetry test on X chromosome，XPAT)。然而该方法并不适用于一般的家系数据。
　　X染色体失活(X chromosome inactivation，XCI)是一种剂量补偿(dose compensation)的机制，通常发生在胚胎发育早期，女性的两条X染色体中的一条发生了失活现象，从而导致女性个体的每一个细胞中，只有一条来自父亲或者母亲的染色体表达。通常，这种X染色体失活的现象是随机的，即女性的两条X染色体失活的比例是50％∶50％。但是，如果失活现象不是随机的，那就称之为偏倚失活(skewed X chromosome inactivation)。偏倚失活通常意味着有75％或以上的细胞失活同一条X染色体，而另外25％的细胞失活另一条X染色体。另外，跟常染色体一样，X染色体还存在着一种特殊现象就是逃逸失活(escape from X chromosome inactivation)的情况，即两条X染色体均表达，没有失活。在一些女性个体中的某一条X染色体上，大概有四分之三的基因是沉默了的，其余四分之一的基因属于逃逸失活区域。然而，现有的X染色体上的关联分析方法中，仅考虑了X染色体失活，并没有合并遗传印记的效应。
　　目的：
　　本文的研究目的主要有以下两点:(1)基于家系数据，提出X染色体质量性状位点上遗传印记效应检验方法;(2)基于两代人的病例对照数据，提出同时合并X染色体失活和印记效应的关联分析方法。
　　方法:
　　(1)考虑X染色体上的位点，参考常染色体上的PPAT方法，利用家系中所有有信息的核心家庭，当感兴趣的SNP位点与疾病之间存在关联的情况下，提出基于家系数据的X染色体亲代不对称检验方法(pedigree parental-asymmetry test on X chromosome，XPPAT)。当某些家系中存在缺失基因型个体时，根据观察到的个体基因型，采用Monte Carlo重抽样方法对缺失个体的基因型进行反推，进一步提出基于家系数据的X染色体Monte Carlo亲代不对称检验方法(Monte Carlo pedigree parental-asymmetry test on X chromosome,XMCPPAT)。通过模拟得到所提出方法的第一类错误率和检验效能。最后，通过一个来自北美的类风湿关节炎基因数据对新提出方法进行实例验证。
　　(2)参考Wang等人提出的最大化似然比的检验方法，针对不同的X染色体失活模式（随机失活、偏倚失活和逃逸失活），将一个个体来自父亲和母亲的等位基因分别进行赋值。对女性的等位基因进行赋值时，用r∈[0，2]来衡量X染色体失活的偏倚程度，然后通过男性的等位基因赋值来体现X染色体是否失活。而这种赋值方式的效应量可以体现印记效应是否存在。这样可以将X染色体失活和遗传印记效应的信息同时合并到Logistic回归模型中。对不同的r值，得到不同的LR值，记为LRr。取LRr的最大值作为检验统计量，注意到其在无关联的原假设下的分布未知。所以，还采用了一种permutation的方法得到新提出方法的经验P值进行统计推断。相应的方法记为XCII方法(association test incorporating X chromosome inactivation and imprinting effects)。通过模拟得到所提出方法的第一类错误率和检验效能。
　　结果:
　　(1)在不同的模拟背景下，包括两种不同的样本量、两组不同的等位基因频率、三种不同的近交系数和五种不同的遗传印记效应模型，模拟结果表明所提出的方法可以很好地控制第一类错误率。并且，新方法的检验效能比现有的XPAT方法要高出很多。通过Monte Carlo重抽样方法对缺失基因型进行反推，XMCPPAT方法在存在缺失基因型个体时能控制第一类错误率。再者，由于女性中近交系数的改变对XPPAT和XMCPPAT方法的结果几乎不产生影响，因而新提出的XPPAT和XMCPPAT方法不受女性中Hardy-Weinberg平衡律（Hardy-Weinberg equilibrium，HWE）是否成立的影响。值得注意的是，对适用于缺失数据的XMCPPAT方法，采用了三种方式对等位基因频率进行处理，包括:真实的等位基因频率、由女性奠基者估算出的等位基因频率和由男性奠基者估算出的等位基因频率，将统计量分别记为XMCPPATt、XMCPPATf和XMCPPATm。另外，基于完整数据，用XPPATfull进行检验，作为所有方法的金标准。从结果上看，应用女性奠基者估算出的等位基因频率的XMCPPATf方法与基于完整数据的XPPATfull方法以及基于真实的等位基因频率的XMCPPATt方法检验效能很接近。这意味着，XMCPPATt和XMCPPATf方法可以捕获大部分缺失基因型个体的信息。然而，XMCPPATm不能很好地控制第一类错误率。
　　(2)模拟结果表明，XCII方法可以将第一类错误率稳定在设定的显著性水平附近，说明了XCII方法用以检验合并X染色体失活与印记效应信息的位点与疾病之间关联的有效性。在检验效能方面，当X染色体失活和遗传印记效应同时存在时，XCII方法的检验效能整体上比Wang的方法高。当仅有X染色体失活而印记效应不存在时，XCII方法的检验效能仍比Wang的方法高。因此，考虑X染色体上的关联分析研究时，更推荐使用第三章新提出的XCII方法。
　　结论:
　　(1)新提出的XPPAT和XMCPPAT检验统计量可以有效地检验基于家系数据的X染色体质量性状位点上的印记效应，XMCPPAT方法比现有的方法检验效能更高，更适用于对包含缺失基因型个体的数据进行印记效应的检验。
　　(2)考虑X染色体上的关联分析研究时，很多疾病同时跟X染色体失活和遗传印记效应相关联，相比只考虑X染色体失活的Wang的方法，XCII方法的表现更好。因此，推荐使用同时考虑X染色体失活和印记效应的XCII方法。

目录

摘要

第一章 前言

1．1 基本概念

1．2 关联分析

1．3 遗传印记效应

1．4 X染色体失活

1．5 论文组成部分

参考文献

第二章 基于家系数据的X染色体质量性状位点E遗传印记效应的检验方法

2．1 引言

2．2 方法

2．3 结果

2．4 附录

参考文献

第三章 同时合并X染色体失活和印记效应的关联分析方法

3．1 引言

3．2 方法

3．3 模拟结果

3．4 附录

参考文献

第四章 讨论与结论

4．2 同时合并X染色体失活和印记效应的关联分析方法

4．3 研究不足

第五章 研究展望

攻读学位期间成果

致谢

声明