Generalized estimating equations for genome-wide association studies using longitudinal phenotype data

Many longitudinal cohort studies have both genome-wide measures of genetic variation and repeated measures of phenotypes and environmental exposures. Genome-wide association study analyses have typically used only cross-sectional data to evaluate quantitative phenotypes and binary traits. Incorporation of repeated measures may increase power to detect associations, but also requires specialized analysis methods. Here we discuss one such method – generalized estimating equations (GEE) – in the contexts of analysis of main effects of rare genetic variants and analysis of gene-environment interactions. We illustrate the potential for increased power using GEE analyses instead of cross-sectional analyses. We also address challenges that arise, such as the need for small-sample corrections when the minor allele frequency of a genetic variant and/or the prevalence of an environmental exposure is low. To illustrate methods for detection of gene-drug interactions on a genome-wide scale, using repeated measures data, we conduct single-study analyses and meta-analyses across studies in three large cohort studies participating in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium – the Atherosclerosis Risk in Communities (ARIC) study, the Cardiovascular Health Study (CHS), and the Rotterdam Study (RS)

Drug-gene interactions of antihypertensive medications and risk of incident cardiovascular disease: a pharmacogenomics study from the CHARGE consortium

Background Hypertension is a major risk factor for a spectrum of cardiovascular diseases (CVD), including myocardial infarction, sudden death, and stroke. In the US, over 65 million people have high blood pressure and a large proportion of these individuals are prescribed antihypertensive medications. Although large long-term clinical trials conducted in the last several decades have identified a number of effective antihypertensive treatments that reduce the risk of future clinical complications, responses to therapy and protection from cardiovascular events vary among individuals. Methods Using a genome-wide association study among 21,267 participants with pharmaceutically treated hypertension, we explored the hypothesis that genetic variants might influence or modify the effectiveness of common antihypertensive therapies on the risk of major cardiovascular outcomes. The classes of drug treatments included angiotensin-converting enzyme inhibitors, beta-blockers, calcium channel blockers, and diuretics. In the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, each study performed array-based genome-wide genotyping, imputed to HapMap Phase II reference panels, and used additive genetic models in proportional hazards or logistic regression models to evaluate drug-gene interactions for each of four therapeutic drug classes. We used meta-analysis to combine study-specific interaction estimates for approximately 2 million single nucleotide polymorphisms (SNPs) in a discovery analysis among 15,375 European Ancestry participants (3,527 CVD cases) with targeted follow-up in a case-only study of 1,751 European Ancestry GenHAT participants as well as among 4,141 African-Americans (1,267 CVD cases). Results Although drug-SNP interactions were biologically plausible, exposures and outcomes were well measured, and power was sufficient to detect modest interactions, we did not identify any statistically significant interactions from the four antihypertensive therapy meta-analyses (Pinteraction > 5.0×10−8). Similarly, findings were null for meta-analyses restricted to 66 SNPs with significant main effects on coronary artery disease or blood pressure from large published genome-wide association studies (Pinteraction ≥ 0.01). Our results suggest that there are no major pharmacogenetic influences of common SNPs on the relationship between blood pressure medications and the risk of incident CVD

Genome-wide analysis of mitochondrial DNA copy number reveals loci implicated in nucleotide metabolism, platelet activation, and megakaryocyte proliferation

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10(–15)) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10(–8)) and mtDNA replication (p = 1.2 × 10(–7)). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10(–4)). SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00439-021-02394-w

Genome-wide Association Study of Susceptibility to Particulate Matter–Associated QT Prolongation

BACKGROUND: Ambient particulate matter (PM) air pollution exposure has been associated with increases in QT interval duration (QT). However, innate susceptibility to PM-associated QT prolongation has not been characterized. OBJECTIVE: To characterize genetic susceptibility to PM-associated QT prolongation in a multi-racial/ethnic, genome-wide association study (GWAS). METHODS: Using repeated electrocardiograms (1986–2004), longitudinal data on PM<10 μm in diameter (PM10), and generalized estimating equations methods adapted for low-prevalence exposure, we estimated approximately 2.5×106 SNP×PM10 interactions among nine Women’s Health Initiative clinical trials and Atherosclerosis Risk in Communities Study subpopulations (n=22,158), then combined subpopulation-specific results in a fixed-effects, inverse variance-weighted meta-analysis. RESULTS: A common variant (rs1619661; coded allele: T) significantly modified the QT-PM10 association (p=2.11×10−8). At PM10 concentrations >90th percentile, QT increased 7 ms across the CC and TT genotypes: 397 (95% confidence interval: 396, 399) to 404 (403, 404) ms. However, QT changed minimally across rs1619661 genotypes at lower PM10 concentrations. The rs1619661 variant is on chromosome 10, 132 kilobase (kb) downstream from CXCL12, which encodes a chemokine, stromal cell-derived factor 1, that is expressed in cardiomyocytes and decreases calcium influx across the L-type Ca2+ channel. CONCLUSIONS: The findings suggest that biologically plausible genetic factors may alter susceptibility to PM10-associated QT prolongation in populations protected by the U.S. Environmental Protection Agency’s National Ambient Air Quality Standards. Independent replication and functional characterization are necessary to validate our findings. https://doi.org/10.1289/EHP34

Pharmacogenomics study of thiazide diuretics and QT interval in multi-ethnic populations : the cohorts for heart and aging research in genomic epidemiology

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common single-nucleotide polymorphisms (SNPs) modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, transethnic and cross-phenotype genome-wide analyses of European (66%), African American (15%) and Hispanic (19%) populations (N = 78 199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P <5 x 10(-8)), we found suggestive evidence (P <5 x 10(-6)) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (for example, NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide-SNP interactions.Peer reviewe

Obesity Partially Mediates the Diabetogenic Effect of Lowering LDL Cholesterol

OBJECTIVE LDL cholesterol (LDLc)-lowering drugs modestly increase body weight and type 2 diabetes risk, but the extent to which the diabetogenic effect of lowering LDLc is mediated through increased BMI is unknown. RESEARCH DESIGN AND METHODS We conducted summary-level univariable and multivariable Mendelian randomization (MR) analyses in 921,908 participants to investigate the effect of lowering LDLc on type 2 diabetes risk and the proportion of this effect mediated through BMI. We used data from 92,532 participants from 14 observational studies to replicate findings in individual-level MR analyses. RESULTS A 1-SD decrease in genetically predicted LDLc was associated with increased type 2 diabetes odds (odds ratio [OR] 1.12 [95% CI 1.01, 1.24]) and BMI (b 5 0.07 SD units [95% CI 0.02, 0.12]) in univariable MR analyses. The multivariable MR analysis showed evidence of an indirect effect of lowering LDLc on type 2 diabetes through BMI (OR 1.04 [95% CI 1.01, 1.08]) with a proportion mediated of 38% of the total effect (P 5 0.03). Total and indirect effect estimates were similar across a number of sensitivity analyses. Individual-level MR analyses confirmed the indirect effect of lowering LDLc on type 2 diabetes through BMI with an estimated proportion mediated of 8% (P 5 0.04). CONCLUSIONS These findings suggest that the diabetogenic effect attributed to lowering LDLc is partially mediated through increased BMI. Our results could help advance understanding of adipose tissue and lipids in type 2 diabetes pathophysiology and inform strategies to reduce diabetes risk among individuals taking LDLc-lowering medications

Ancestry-specific associations identified in genome-wide combined-phenotype study of red blood cell traits emphasize benefits of diversity in genomics

Background: Quantitative red blood cell (RBC) traits are highly polygenic clinically relevant traits, with approximately 500 reported GWAS loci. The majority of RBC trait GWAS have been performed in European- or East Asian-ancestry populations, despite evidence that rare or ancestry-specific variation contributes substantially to RBC trait heritability. Recently developed combined-phenotype methods which leverage genetic trait correlation to improve statistical power have not yet been applied to these traits. Here we leveraged correlation of seven quantitative RBC traits in performing a combined-phenotype analysis in a multi-ethnic study population. Results: We used the adaptive sum of powered scores (aSPU) test to assess combined-phenotype associations between ~ 21 million SNPs and seven RBC traits in a multi-ethnic population (maximum n = 67,885 participants; 24% African American, 30% Hispanic/Latino, and 43% European American; 76% female). Thirty-nine loci in our multi-ethnic population contained at least one significant association signal (p 5%) across all ancestral populations. Nineteen additional independent association signals were identified at seven known loci (HFE, KIT, HBS1L/MYB, CITED2/FILNC1, ABO, HBA1/2, and PLIN4/5). For example, the HBA1/2 locus contained 14 conditionally independent association signals, 11 of which were previously unreported and are specific to African and Amerindian ancestries. One variant in this region was common in all ancestries, but exhibited a narrower LD block in African Americans than European Americans or Hispanics/Latinos. GTEx eQTL analysis of all independent lead SNPs yielded 31 significant associations in relevant tissues, over half of which were not at the gene immediately proximal to the lead SNP. Conclusion: This work identified seven loci containing multiple independent association signals for RBC traits using a combined-phenotype approach, which may improve discovery in genetically correlated traits. Highly complex genetic architecture at the HBA1/2 locus was only revealed by the inclusion of African Americans and Hispanics/Latinos, underscoring the continued importance of expanding large GWAS to include ancestrally diverse populations. © 2020 The Author(s)

Pharmacogenomics of statin-related myopathy:Meta-analysis of rare variants from whole-exome sequencing

AIMS:Statin-related myopathy (SRM), which includes rhabdomyolysis, is an uncommon but important adverse drug reaction because the number of people prescribed statins world-wide is large. Previous association studies of common genetic variants have had limited success in identifying a genetic basis for this adverse drug reaction. We conducted a multi-site whole-exome sequencing study to investigate whether rare coding variants confer an increased risk of SRM. METHODS AND RESULTS:SRM 3-5 cases (N = 505) and statin treatment-tolerant controls (N = 2047) were recruited from multiple sites in North America and Europe. SRM 3-5 was defined as symptoms consistent with muscle injury and an elevated creatine phosphokinase level >4 times upper limit of normal without another likely cause of muscle injury. Whole-exome sequencing and variant calling was coordinated from two analysis centres, and results of single-variant and gene-based burden tests were meta-analysed. No genome-wide significant associations were identified. Given the large number of cases, we had 80% power to identify a variant with minor allele frequency of 0.01 that increases the risk of SRM 6-fold at genome-wide significance. CONCLUSIONS:In this large whole-exome sequencing study of severe statin-related muscle injury conducted to date, we did not find evidence that rare coding variants are responsible for this adverse drug reaction. Larger sample sizes would be required to identify rare variants with small effects, but it is unclear whether such findings would be clinically actionable