Statistical genetics of aging-related genomic and phenotypic change

为了帮助分析和理解受许多基因和环境因素影响的与衰老相关的复杂特征，我们提出了三种统计算法，用于全基因组基因分型和高通量测序研究。我们建议的新统计方法提供了分析其他类型数据的方法，例如来自测序的线粒体DNA(MtDNA)变体，或用于全基因组关联研究(Gwas)的X染色体上的变体，以及具有更复杂结构的数据(例如，多个相关性状)。为了测试这些算法，我们利用撒丁岛项目的特殊功能(参见年度报告AG000675-07)，该项目收集了&gt；300数量性状的纵向数据以及撒丁岛创始人群体的全基因组遗传数据。 要分析线粒体DNA变异及其对衰老相关性状的可能影响，为核DNA开发的基因型调用和分析程序是不够的，因为每个细胞都有100-10,000个mtDNA拷贝，这些拷贝可以在任何位置变化(异质性)，因此可以在不同拷贝中的任何位置拥有4个碱基中的每一个。我们已经开发了一种算法，旨在识别mtDNA的变异；它结合了每个序列读取中每个碱基的测序错误率，并且灵活地允许在所有个体的变异位置上有不同的等位基因部分。我们的方法进一步适用于环形线粒体基因组，这是与大多数映射算法假设的线性染色体的关键区别。我们正在从2000名撒丁岛计划参与者的全基因组测序中评估淋巴细胞mtDNA序列中的同质性和异质性。到目前为止的结果提供了关于线粒体DNA单倍群和撒丁岛同质性和异质性遗传的信息。正如预期的那样，母亲和她们的孩子基本上分享了所有的同质性，但异质性的比例较小。总体异质性随着年龄的增长而增加，但到目前为止估计的斜率很小，在20岁至80岁之间平均增加1个异质性，次要等位基因比例阈值为4%。 为了利用相关性状之间的相关性，从而增加遗传学研究的统计能力，我们正在开发一种方法来寻找对多个数量性状具有多效性的基因/变异。我们的方法将一组相关性状和一组SNPs(例如，由基因边界定义)投影到它们各自的正交主成分中，以便我们能够使用新的汇总统计量来联合测试性状和SNPs之间的关联。由于正交性，关联的重要性可以在零假设下使用模拟来有效地评估，而不是计算密集的排列。我们将我们的方法应用于撒丁岛项目的数据，其中我们首先关注了三个血脂特性，即高密度脂蛋白、低密度脂蛋白和甘油三酯，并使用RefSeq基因边界将SNPs分组为基因单元。为了证明我们的方法能够识别与多个血脂性状相关的基因，我们使用了发表在Teslovich等人的《自然》(2010)上的结果，其中95个血脂基因座被识别，21个基因座与多个血脂性状相关。我们能够证明我们的方法可以丰富这些多效性基因座：例如，当用我们的方法对95个基因座进行排序时，前20%的基因座包含了原始研究中40%的多效性基因座。我们的方法能够检测多个性状和多个遗传变异之间的联合关联。当一个基因对多个性状有适度的影响时，它将具有显著的优势，而标准的GWAS无法检测到这些性状。

To help to analyze and understand aging-related "complex" traits that are affected by many genes and environmental factors, we propose to develop three statistical algorithms for the analyses of genome-wide genotyping and high-throughput sequencing studies. Our proposed new statistical methods provide means to analyze additional types of data e.g., mitochondrial DNA (mtDNA) variants from sequencing, or variants on the X chromosome for genome-wide association studies (GWAS) and data with more complicated structures (e.g., multiple related traits). To test these algorithms, we take advantage of the special features of the SardiNIA project (see Annual Report AG000675-07), which has collected longitudinal data for >300 quantitative traits together with the whole-genome genetic data in the founder Sardinia population. To analyze mitochondrial DNA variation and its possible effects on aging-related traits, the genotype-calling and analytic programs developed for nuclear DNA are not adequate, because each cell has 100-10,000 mtDNA copies that can vary at any site (heteroplasmy), and can therefore have each of the 4 bases at any position in various copies. We have developed an algorithm that is targeted to identify variants in mtDNA; it incorporates the sequencing error rate of each base in each sequence read and is flexible to allow for different allele fractions at a variant site across all individuals. Our procedure is further adapted to the circular mitochondrial genome, a key difference from the linear chromosomes assumed by most mapping algorithms. We are assessing homoplasmies and heteroplasmies in mtDNA sequences of lymphocytes from whole-genome sequencing of 2,000 SardiNIA Project participants. The results to date provide information about mtDNA haplogroups and the inheritance of homo- and heteroplasmies in Sardinia. As expected, mothers and their children share essentially all homoplasmies but a lesser proportion of heteroplasmies. The overall heteroplasmy increases with age, but the slope is small in the estimates thus far, yielding an average increase of 1 heteroplasmy between ages 20 and 80 with the minor allele fraction threshold at 4%. To take advantage of correlations between related traits and hence to increase statistical power for genetic studies, we are developing a method to search for genes/variants that have pleiotropic effects on multiple quantitative traits. Our method projects a group of related traits and a set of SNPs (defined, for example, by gene boundaries) into their respective orthogonal principal components such that we are able to jointly test the association between traits and SNPs using a new summary statistic. Because of the orthogonality, the significance of association can be efficiently evaluated using simulations under the null hypothesis instead of more computationally intensive permutations. We apply our method to the SardiNIA project data, where we have first focused on three lipid traits HDL, LDL, and Triglycerides, and use RefSeq gene boundaries to group SNPs into gene units. To show that our method is able to identify genes that are associated with more than one blood lipid trait, we use results published in Teslovich et al., Nature (2010), in which 95 loci for blood lipids were identified and 21 loci were associated with multiple lipid traits. We are able to show that our method can enrich those pleiotropic loci: for example, when the 95 loci are ranked by our method, the top 20% loci include 40% of the pleiotropic loci in the original study. Our method is able to detect joint associations between multiple traits and multiple genetic variants. It will have significant advantages when a gene has moderate effects on multiple traits that a standard GWAS is unable to detect.