Xiuqing Guo

Actionable exomic incidental findings in 6503 participants: challenges of variant classification

Recommendations for laboratories to report incidental findings from genomic tests have stimulated interest in such results. In order to investigate the criteria and processes for assigning the pathogenicity of specific variants and to estimate the frequency of such incidental findings in patients of European and African ancestry, we classified potentially actionable pathogenic single-nucleotide variants (SNVs) in all 4300 European- and 2203 African-ancestry participants sequenced by the NHLBI Exome Sequencing Project (ESP). We considered 112 gene-disease pairs selected by an expert panel as associated with medically actionable genetic disorders that may be undiagnosed in adults. The resulting classifications were compared to classifications from other clinical and research genetic testing laboratories, as well as with in silico pathogenicity scores. Among European-ancestry participants, 30 of 4300 (0.7%) had a pathogenic SNV and six (0.1%) had a disruptive variant that was expected to be pathogenic, whereas 52 (1.2%) had likely pathogenic SNVs. For African-ancestry participants, six of 2203 (0.3%) had a pathogenic SNV and six (0.3%) had an expected pathogenic disruptive variant, whereas 13 (0.6%) had likely pathogenic SNVs. Genomic Evolutionary Rate Profiling mammalian conservation score and the Combined Annotation Dependent Depletion summary score of conservation, substitution, regulation, and other evidence were compared across pathogenicity assignments and appear to have utility in variant classification. This work provides a refined estimate of the burden of adult onset, medically actionable incidental findings expected from exome sequencing, highlights challenges in variant classification, and demonstrates the need for a better curated variant interpretation knowledge base.

Set-based tests for genetic association in longitudinal studies

Genetic association studies with longitudinal markers of chronic diseases (e.g., blood pressure, body mass index) provide a valuable opportunity to explore how genetic variants affect traits over time by utilizing the full trajectory of longitudinal outcomes. Since these traits are likely influenced by the joint effect of multiple variants in a gene, a joint analysis of these variants considering linkage disequilibrium (LD) may help to explain additional phenotypic variation. In this article, we propose a longitudinal genetic random field model (LGRF), to test the association between a phenotype measured repeatedly during the course of an observational study and a set of genetic variants. Generalized score type tests are developed, which we show are robust to misspecification of within-subject correlation, a feature that is desirable for longitudinal analysis. In addition, a joint test incorporating gene-time interaction is further proposed. Computational advancement is made for scalable implementation of the proposed methods in large-scale genome-wide association studies (GWAS). The proposed methods are evaluated through extensive simulation studies and illustrated using data from the Multi-Ethnic Study of Atherosclerosis (MESA). Our simulation results indicate substantial gain in power using LGRF when compared with two commonly used existing alternatives: (i) single marker tests using longitudinal outcome and (ii) existing gene-based tests using the average value of repeated measurements as the outcome.

ExomeChip-Wide Analysis of 95 626 Individuals Identifies 10 Novel Loci Associated With QT and JT Intervals

Background: QT interval, measured through a standard ECG, captures the time it takes for the cardiac ventricles to depolarize and repolarize. JT interval is the component of the QT interval that reflects ventricular repolarization alone. Prolonged QT interval has been linked to higher risk of sudden cardiac arrest. Methods and Results: We performed an ExomeChip-wide analysis for both QT and JT intervals, including 209 449 variants, both common and rare, in 17 341 genes from the Illumina Infinium HumanExome BeadChip. We identified 10 loci that modulate QT and JT interval duration that have not been previously reported in the literature using single-variant statistical models in a meta-analysis of 95 626 individuals from 23 cohorts (comprised 83 884 European ancestry individuals, 9610 blacks, 1382 Hispanics, and 750 Asians). This brings the total number of ventricular repolarization associated loci to 45. In addition, our approach of using coding variants has highlighted the role of 17 specific genes for involvement in ventricular repolarization, 7 of which are in novel loci. Conclusions: Our analyses show a role for myocyte internal structure and interconnections in modulating QT interval duration, adding to previous known roles of potassium, sodium, and calcium ion regulation, as well as autonomic control. We anticipate that these discoveries will open new paths to the goal of making novel remedies for the prevention of lethal ventricular arrhythmias and sudden cardiac arrest.

DNA methylation age is associated with an altered hemostatic profile in a multi-ethnic meta-analysis

Many hemostatic factors are associated with age and age-related diseases; however, much remains unknown about the biological mechanisms linking aging and hemostatic factors. DNA methylation is a novel means by which to assess epigenetic aging, which is a measure of age and the aging processes as determined by altered epigenetic states. We used a meta-analysis approach to examine the association between measures of epigenetic aging and hemostatic factors, as well as a clotting time measure. For fibrinogen, we performed European and African ancestry-specific meta-analyses which were then combined via a random effects meta-analysis. For all other measures we could not estimate ancestry-specific effects and used a single fixed effects meta-analysis. We found that 1-year higher extrinsic epigenetic age as compared with chronological age was associated with higher fibrinogen (0.004 g/L/y; 95% confidence interval, 0.001-0.007; P = .01) and plasminogen activator inhibitor 1 (PAI-1; 0.13 U/mL/y; 95% confidence interval, 0.07-0.20; P = 6.6 × 10-5) concentrations, as well as lower activated partial thromboplastin time, a measure of clotting time. We replicated PAI-1 associations using an independent cohort. To further elucidate potential functional mechanisms, we associated epigenetic aging with expression levels of the PAI-1 protein encoding gene (SERPINE1) and the 3 fibrinogen subunit-encoding genes (FGA, FGG, and FGB) in both peripheral blood and aorta intima-media samples. We observed associations between accelerated epigenetic aging and transcription of FGG in both tissues. Collectively, our results indicate that accelerated epigenetic aging is associated with a procoagulation hemostatic profile, and that epigenetic aging may regulate hemostasis in part via gene transcription.

Genetically determined plasma lipid levels and risk of diabetic retinopathy: A mendelian randomization study

Results from observational studies examining dyslipidemia as a risk factor for diabetic retinopathy (DR) have been inconsistent. We evaluated the causal relationship between plasma lipids and DR using a Mendelian randomization approach. We pooled genome-wide association studies summary statistics from 18 studies for two DR phenotypes: any DR (N = 2,969 case and 4,096 control subjects) and severe DR (N = 1,277 case and 3,980 control subjects). Previously identified lipid-associated single nucleotide polymorphisms served as instrumental variables. Meta-analysis to combine the Mendelian randomization estimates from different cohorts was conducted. There was no statistically significant change in odds ratios of having any DR or severe DR for any of the lipid fractions in the primary analysis that used single nucleotide polymorphisms that did not have a pleiotropic effect on another lipid fraction. Similarly, there was no significant association in the Caucasian and Chinese subgroup analyses. This study did not show evidence of a causal role of the four lipid fractions on DR. However, the study had limited power to detect odds ratios less than 1.23 per SD in genetically induced increase in plasma lipid levels, thus we cannot exclude that causal relationships with more modest effect sizes exist.

Evaluation of the causal effect of fibrinogen on incident coronary heart disease via Mendelian randomization

Background: Fibrinogen is an essential hemostatic factor and cardiovascular disease risk factor. Early attempts at evaluating the causal effect of fibrinogen on coronary heart disease (CHD) and myocardial infraction (MI) using Mendelian randomization (MR) used single variant approaches, and did not take advantage of recent genome-wide association studies (GWAS) or multi-variant, pleiotropy robust MR methodologies. Methods and Findings: We evaluated evidence for a causal effect of fibrinogen on both CHD and MI using MR. We used both an allele score approach and pleiotropy robust MR models. The allele score was composed of 38 fibrinogen-associated variants from recent GWAS. Initial analyses using the allele score incorporated data from 11 European-ancestry prospective cohorts to examine incidence CHD and MI. We also applied 2 sample MR methods with data from a prevalent CHD and MI GWAS. Results are given in terms of the hazard ratio (HR) or odds ratio (OR), depending on the study design, and associated 95% confidence interval (CI). In single variant analyses no causal effect of fibrinogen on CHD or MI was observed. In multi-variant analyses using incidence CHD cases and the allele score approach, the estimated causal effect (HR) of a 1 g/L higher fibrinogen concentration was 1.62 (CI = 1.12, 2.36) when using incident cases and the allele score approach. In 2 sample MR analyses that accounted for pleiotropy, the causal estimate (odds ratio) was reduced to 1.18 (CI = 0.98, 1.42) and 1.09 (CI = 0.89, 1.33) in the 2 most precise (smallest CI) models, out of 4 models evaluated. In the 2 sample MR analyses for MI, there was only very weak evidence of a causal effect in only 1 out of 4 models. Conclusions: A small causal effect of fibrinogen on CHD is observed using multi-variant MR approaches which account for pleiotropy, but not single variant MR approaches. Taken together, results indicate that even with large sample sizes and multi-variant approaches, MR analyses still cannot exclude the null when estimating the causal effect of fibrinogen on CHD, but that any potential causal effect is likely to be much smaller than observed in epidemiological studies.Background Fibrinogen is an essential hemostatic factor and cardiovascular disease risk factor. Early attempts at evaluating the causal effect of fibrinogen on coronary heart disease (CHD) and myocardial infraction (MI) using Mendelian randomization (MR) used single variant approaches, and did not take advantage of recent genome-wide association studies (GWAS) or multi-variant, pleiotropy robust MR methodologies. Methods and Findings We evaluated evidence for a causal effect of fibrinogen on both CHD and MI using MR. We used both an allele score approach and pleiotropy robust MR models. The allele score was composed of 38 fibrinogen-associated variants from recent GWAS. Initial analyses using the allele score incorporated data from 11 European-ancestry prospective cohorts to examine incidence CHD and MI. We also applied 2 sample MR methods with data from a prevalent CHD and MI GWAS. Results are given in terms of the hazard ratio (HR) or odds ratio (OR), depending on the study design, and associated 95% confidence interval (CI). In single variant analyses no causal effect of fibrinogen on CHD or MI was observed. In multi-variant analyses using incidence CHD cases and the allele score approach, the estimated causal effect (HR) of a 1 g/L higher fibrinogen concentration was 1.62 (CI = 1.12, 2.36) when using incident cases and the allele score approach. In 2 sample MR analyses that accounted for pleiotropy, the causal estimate (OR) was reduced to 1.18 (CI = 0.98, 1.42) and 1.09 (CI = 0.89, 1.33) in the 2 most precise (smallest CI) models, out of 4 models evaluated. In the 2 sample MR analyses for MI, there was only very weak evidence of a causal effect in only 1 out of 4 models. Conclusions A small causal effect of fibrinogen on CHD is observed using multi-variant MR approaches which account for pleiotropy, but not single variant MR approaches. Taken together, results indicate that even with large sample sizes and multi-variant approaches MR analyses still cannot exclude the null when estimating the causal effect of fibrinogen on CHD, but that any potential causal effect is likely to be much smaller than observed in epidemiological studies. Author Summary Initial Mendelian Randomization (MR) analyses of the causal effect of fibrinogen on coronary heart disease (CHD) utilized single variants and did not take advantage of modern, multivariant approaches. This manuscript provides an important update to these initial analyses by incorporating larger sample sizes and employing multiple, modern multi-variant MR approaches to account for pleiotropy. We used incident cases to perform a MR study of the causal effect of fibrinogen on incident CHD and the nested outcome of myocardial infarction (MI) using an allele score approach. Then using data from a case-control genome-wide association study for CHD and MI we performed two sample MR analyses with multiple, pleiotropy robust approaches. Overall, the results indicated that associations between fibrinogen and CHD in observational studies are likely upwardly biased from any underlying causal effect. Single variant MR approaches show little evidence of a causal effect of fibrinogen on CHD or MI. Multi-variant MR analyses of fibrinogen on CHD indicate there may be a small positive effect, however this result needs to be interpreted carefully as the 95% confidence intervals were still consistent with a null effect. Multi-variant MR approaches did not suggest evidence of even a small causal effect of fibrinogen on MI.

Adiposity-Independent Effects of Aging on Insulin Sensitivity and Clearance in Humans and Mice

Aims/hypothesis Aging is associated with impaired insulin sensitivity and increased prevalence of type 2 diabetes. However, it remains unclear whether aging-related insulin resistance is due to age per se, or increased adiposity associated with advanced age. In the present study, we investigate the impact of aging on insulin sensitivity independent of changes in body composition. Methods Cohorts of C57BL/6J male mice at 4-8 months of age (‘young’) and 18-27 mo (‘aged’) exhibiting similar body composition were characterized with static (plasma glucose and insulin levels) and dynamic (glucose and insulin tolerance tests) measures of glucose metabolism on chow and high-fat diets. Insulin sensitivity was assessed by hyperinsulinemic-euglycemic clamp analysis. The relationship between aging and insulin resistance in humans was investigated in 1,250 non-diabetic Mexican-American individuals who underwent hyperinsulinemic-euglycemic clamps. Results In mice with similar body composition, age had no detrimental effect on plasma glucose and insulin levels. However, aged mice demonstrated mildly, but reproducibly, improved glucose tolerance on both chow and high-fat diets due to increased glucose-stimulated insulin secretion. Moreover, hyperinsulinemic-euglycemic clamps revealed impaired insulin sensitivity and reduced insulin clearance in aged mice on both diets. Consistent with results in the mouse, age remained an independent determinant of insulin resistance after adjustment for body composition in Mexican-Americn males. Advanced age was also associated with diminished insulin clearance, but this effect was dependent on increased BMI. Conclusions/interpretation This study demonstrates for the first time that aging per se impairs insulin sensitivity independent of adiposity in mice and humans. These results raise the possibility that the pathogenetic mechanisms of age-related and obesity-associated insulin resistance are distinct. Abbreviations BAI body adiposity index GEE generalized estimating equations HF high-fat IQR interquartile range MCRI metabolic clearance rate of insulin T2D type 2 diabetes

Rare-variant association tests in longitudinal studies, with an application to the Multi-Ethnic Study of Atherosclerosis (MESA)

Over the past few years, an increasing number of studies have identified rare variants that contribute to trait heritability. Due to the extreme rarity of some individual variants, gene‐based association tests have been proposed to aggregate the genetic variants within a gene, pathway, or specific genomic region as opposed to a one‐at‐a‐time single variant analysis. In addition, in longitudinal studies, statistical power to detect disease susceptibility rare variants can be improved through jointly testing repeatedly measured outcomes, which better describes the temporal development of the trait of interest. However, usual sandwich/model‐based inference for sequencing studies with longitudinal outcomes and rare variants can produce deflated/inflated type I error rate without further corrections. In this paper, we develop a group of tests for rare‐variant association based on outcomes with repeated measures. We propose new perturbation methods such that the type I error rate of the new tests is not only robust to misspecification of within‐subject correlation, but also significantly improved for variants with extreme rarity in a study with small or moderate sample size. Through extensive simulation studies, we illustrate that substantially higher power can be achieved by utilizing longitudinal outcomes and our proposed finite sample adjustment. We illustrate our methods using data from the Multi‐Ethnic Study of Atherosclerosis for exploring association of repeated measures of blood pressure with rare and common variants based on exome sequencing data on 6,361 individuals.

An Empirical Comparison of Joint and Stratified Frameworks for Studying G × E Interactions: Systolic Blood Pressure and Smoking in the CHARGE Gene-Lifestyle Interactions Working Group: Joint vs. Stratified Framework for Interaction

Studying gene‐environment (G × E) interactions is important, as they extend our knowledge of the genetic architecture of complex traits and may help to identify novel variants not detected via analysis of main effects alone. The main statistical framework for studying G × E interactions uses a single regression model that includes both the genetic main and G × E interaction effects (the “joint” framework). The alternative “stratified” framework combines results from genetic main‐effect analyses carried out separately within the exposed and unexposed groups. Although there have been several investigations using theory and simulation, an empirical comparison of the two frameworks is lacking. Here, we compare the two frameworks using results from genome‐wide association studies of systolic blood pressure for 3.2 million low frequency and 6.5 million common variants across 20 cohorts of European ancestry, comprising 79,731 individuals. Our cohorts have sample sizes ranging from 456 to 22,983 and include both family‐based and population‐based samples. In cohort‐specific analyses, the two frameworks provided similar inference for population‐based cohorts. The agreement was reduced for family‐based cohorts. In meta‐analyses, agreement between the two frameworks was less than that observed in cohort‐specific analyses, despite the increased sample size. In meta‐analyses, agreement depended on (1) the minor allele frequency, (2) inclusion of family‐based cohorts in meta‐analysis, and (3) filtering scheme. The stratified framework appears to approximate the joint framework well only for common variants in population‐based cohorts. We conclude that the joint framework is the preferred approach and should be used to control false positives when dealing with low‐frequency variants and/or family‐based cohorts.

General Framework for Meta‐Analysis of Haplotype Association Tests

For complex traits, most associated single nucleotide variants (SNV) discovered to date have a small effect, and detection of association is only possible with large sample sizes. Because of patient confidentiality concerns, it is often not possible to pool genetic data from multiple cohorts, and meta‐analysis has emerged as the method of choice to combine results from multiple studies. Many meta‐analysis methods are available for single SNV analyses. As new approaches allow the capture of low frequency and rare genetic variation, it is of interest to jointly consider multiple variants to improve power. However, for the analysis of haplotypes formed by multiple SNVs, meta‐analysis remains a challenge, because different haplotypes may be observed across studies. We propose a two‐stage meta‐analysis approach to combine haplotype analysis results. In the first stage, each cohort estimate haplotype effect sizes in a regression framework, accounting for relatedness among observations if appropriate. For the second stage, we use a multivariate generalized least square meta‐analysis approach to combine haplotype effect estimates from multiple cohorts. Haplotype‐specific association tests and a global test of independence between haplotypes and traits are obtained within our framework. We demonstrate through simulation studies that we control the type‐I error rate, and our approach is more powerful than inverse variance weighted meta‐analysis of single SNV analysis when haplotype effects are present. We replicate a published haplotype association between fasting glucose‐associated locus (G6PC2) and fasting glucose in seven studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and we provide more precise haplotype effect estimates.