Statistical genetics of aging-related genomic and phenotypic change

为了帮助分析和理解受许多基因和环境因素影响的与衰老相关的“复杂”性状，我们沿着开发用于全基因组基因分型分析和高通量测序研究的统计算法的路径。我们提出的新的计算工具提供了分析其他类型数据的方法，例如，识别线粒体DNA （mtDNA）变异，并从全基因组序列中有效地估计mtDNA拷贝数。对于算法的实验测试，我们正在利用撒丁岛项目的特殊优势（见年度报告AG000675）来帮助组装撒丁岛创始人人群的线粒体序列数据和多种表型数据。与此同时，我们正在与NIA和其他合作者进行全基因组关联研究（GWAS）和流行病学分析，以研究撒丁岛队列的一系列年龄相关特征。例如，在过去的一年里，我们参与了与虚弱相关的特征（步行速度、握力和骨密度）和听力能力作为年龄和性别的函数的流行病学分析。

To help to analyze and understand aging-related "complex" traits that are affected by many genes and environmental factors, we have followed the path of developing statistical algorithms for the analyses of genome-wide genotyping and high-throughput sequencing studies. Our proposed new computational tools provide means to analyze additional types of data e.g., to identify mitochondrial DNA (mtDNA) variants and to estimate mtDNA copy number efficiently from whole-genome sequences. For experimental tests of the algorithms, we are capitalizing on the special advantages of the SardiNIA project (see Annual Report AG000675) to help in the assembly of mitochondrial sequence data and multiple phenotypic data in the founder Sardinian population. At the same time, we are carrying out genome-wide association studies (GWAS) and epidemiological analyses with NIA and other collaborators for a series of age-related traits in the Sardinian cohort. In the past year, for example, this has involved us in epidemiological analyses for frailty-related traits (walking speed, grip strength and bone density) and hearing capacity as a function of age and sex. In order to conduct analyses on large-scale consortium data to study mtDNA variation and copy number, we have developed two computational programs, providing a general solution for the analysis of mtDNA dynamics based on whole-genome sequencing studies. One program (mitoCaller) is designed specifically to identify mtDNA variants; the other (mitoCalc) infers mtDNA copy number in a cell directly from genome sequences. Applying the programs to leukocyte sequences of 2,000 SardiNIA participants, we have shown that heteroplasmies (mtDNA variants with more than one allele at a site) increase with age, and that copy number is relatively highly heritable and is correlated with metabolic traits, particularly central fat levels. In more recent work, we have increased the speed ofmitoCalc 100-fold (fastMitoCalc). The new program is being applied to white cells of 50,000 deep sequenced individuals (TOPMed program, NHLBI), for GWAS on copy number. We have also initiated analyses of possible effects of mtDNA copy number or variants on traits and diseases. In another study, we have created a program that uses machine learning methods to measure physiological rates of aging of individuals. We assess the extent to which an individual's physiological age could be determined as a composite score inferred from a broad range of biochemical and physiological data. Data were collected in the SardiNIA population study. We use machine learning strategies on data for 6,000 Sardinian participants, who ranged in age from 12 to 81. The best predictive models are determined from multiple combinations of dimensionality reduction, classification, and regression algorithms. They reach very strong correlations (R > 0.9) between predicted and actual ages, and show relative stability in successive visits of the same individuals (R>0.5). We then define an Effective Rate of Aging (ERA) for each participant, a trait measured as the ratio of an individual's predicted age to his/her chronological age. The inference that individuals have a characteristic rate of aging is supported by findings that in the SardiNIA cohort, the inferred values of ERA shows genetic heritability of 40%. This has been sufficient to initiate genome-wide association studies that identify genetic variants influencing the rate of aging. In an ongoing collaborative effort, we study a special structural feature of DNA, G-quadruplex (G4) structures as potential DNA roadblocks that perturb mitochondrial replication machinery. Computational analysis of mtDNA variants associated with predicted G4-forming structures in two large population cohorts are assessing G4 mtDNA as a factor contributing to genetic variability and disease risk. At the same time, biochemical studies are being performed in a reconstituted system to assess mitochondrial replisome progression through difficult-to-replicate sequences characterized by G4 structure. The study of possible mtDNA replication blockades complements analysis of the unique SardiNIA and InChianti populations to identify G4-associated mutations that affect age-related traits.