An empirical comparison of meta-analysis and mega-analysis of individual participant data for identifying gene-environment interactions

pre-printMeta-analysis combining results from multiple studies is a standard practice in GWAS. For genetic main effects, meta-analysis has been shown to provide comparable results as mega-analysis that jointly analyzes the pooled data from the available studies. Gene-environment interaction (GEI) studies are an important component of genetic epidemiology research since they can explain a part of the missing heritability, elucidate the biological networks underlying disease risk, and identify individuals at high risk for disease. However, it is not known whether meta- and mega-analyses of interactions also yield comparable results. In this study, we investigate whether both approaches provide comparable results for identifying interaction effects using empirical data from 4 studies: the Framingham Heart Study, GENOA, HERITAGE and HyperGEN. We performed meta-analysis of cohort-specific results and mega-analysis by analyzing the pooled data from all 4 studies. We used the standard 1 degree of freedom (df) test of main effect only, the 1 df test of the interaction effect (in the presence of main effect), and the joint 2 df test of main and interaction effects. We found that the results from meta- and mega-analyses were highly consistent for all three tests. The correlation between -log (p) values from the two analyses was 0.89 for the 1 df main effect, 0.90 for the 1 df interaction test, and 0.91 for the joint 2 df test. Although mega-analysis provided slightly better results as expected, both yielded very similar results for the most promising SNPs. Moreover, mega-analysis is not always feasible especially in very large and diverse consortia since pooling of raw data may be limited by the terms of the informed consent. Our study illustrates that meta-analysis can be an effective approach also for identifying interactions in very large consortia without losing appreciable power

Gene-alcohol interactions identify several novel blood pressure loci including a promising locus near SLC16A9

Alcohol consumption is a known risk factor for hypertension, with recent candidate studies implicating gene-alcohol interactions in blood pressure (BP) regulation. We used 6,882 (predominantly) Caucasian participants aged 20 to 80 years from the Framingham SHARe (SNP Health Association Resource) to perform a genome-wide analysis of SNP-alcohol interactions on BP traits. We used a two-step approach in the ABEL suite to examine genetic interactions with three alcohol measures [ounces of alcohol consumed per week, drinks consumed per week, and the number of days drinking alcohol per week] on four BP traits [systolic (SBP), diastolic (DBP), mean arterial (MAP), and pulse (PP) pressure]. In the first step, we fit a linear mixed model of each BP trait onto age, sex, BMI, and antihypertensive medication while accounting for the phenotypic correlation among relatives. In the second step, we conducted 1 degree-of-freedom (df) score tests of the SNP main effect, alcohol main effect, and SNP-alcohol interaction using the maximum likelihood estimates of the parameters from the first step. We then calculated the joint 2 df score test of the SNP main effect and SNP-alcohol interaction using MixABEL. The effect of SNP rs10826334 (near SLC16A9) on SBP was significantly modulated by both the number of alcoholic drinks and the ounces of alcohol consumed per week (p-values of 1.27E-08 and 3.92E-08, respectively). Each copy of the G-allele decreased SBP by 3.79 mmHg in those consuming 14 drinks per week versus a 0.461 mmHg decrease in non-drinkers. Index SNPs in 20 other loci exhibited suggestive (p-value≤1E-06) associations with BP traits by the 1 df interaction test or joint 2df test, including 3 rare variants, one low-frequency variant, and SNPs near/in genes ESRRG, FAM179A, CRIPT-SOCS5, KAT2B,ADCY2, GLI3, ZNF716, SLIT1, PDE3A, KERA-LUM, RNF219-AS1, CLEC3A , FBX015, and IGSF5. SNP -alcohol interactions may enhance discovery of novel variants with large effects that can be targete

Comparison of two methods for analysis of gene-environment interactions in longitudinal family data: The Framingham heart study

Gene–environment interaction (GEI) analysis can potentially enhance gene discovery for common complex traits. However, genome-wide interaction analysis is computationally intensive. Moreover, analysis of longitudinal data in families is much more challenging due to the two sources of correlations arising from longitudinal measurements and family relationships. GWIS of longitudinal family data can be a computational bottleneck. Therefore, we compared two methods for analysis of longitudinal family data: a methodologically sound but computationally demanding method using the Kronecker model (KRC) and a computationally more forgiving method using the hierarchical linear model (HLM). The KRC model uses a Kronecker product of an unstructured matrix for correlations among repeated measures (longitudinal) and a compound symmetry matrix for correlations within families at a given visit. The HLM uses an autoregressive covariance matrix for correlations among repeated measures and a random intercept for familial correlations. We compared the two methods using the longitudinal Framingham heart study (FHS) SHARe data. Specifically, we evaluated SNP–alcohol (amount of alcohol consumption) interaction effects on high density lipoprotein cholesterol (HDLC). Keeping the prohibitive computational burden of KRC in mind, we limited the analysis to chromosome 16, where preliminary cross-sectional analysis yielded some interesting results. Our first important finding was that the HLM provided very comparable results but was remarkably faster than the KRC, making HLM the method of choice. Our second finding was that longitudinal analysis provided smaller P-values, thus leading to more significant results, than cross-sectional analysis. This was particularly pronounced in identifying GEIs. We conclude that longitudinal analysis of GEIs is more powerful and that the HLM method is an optimal method of choice as compared to the computationally (prohibitively) intensive KRC method

The role of SNP-loop diuretic interactions in hypertension across ethnic groups in HyperGEN

Blood pressure (BP) is significantly influenced by genetic factors; however, less than 3% of the BP variance has been accounted for by variants identified from genome-wide association studies (GWAS) of primarily European-descent cohorts. Other genetic influences, including gene-environment (GxE) interactions, may explain more of the unexplained variance in BP. African Americans (AA) have a higher prevalence and earlier age of onset of hypertension (HTN) as compared with European Americans (EA); responses to anti-hypertensive drugs vary across race groups. To examine potential interactions between the use of loop diuretics and HTN traits, we analyzed systolic (SBP) and diastolic (DBP) blood BP from 1,222 AA and 1,231 EA participants in the Hypertension Genetic Epidemiology Network (HyperGEN). Population-specific score tests were used to test associations of SBP and DBP, using a panel of genotyped and imputed single nucleotide polymorphisms (SNPs) for African Americans (2.9 million SNPs) and European Americans (2.3 million SNPs). Several promising loci were identified through gene-loop diuretic interactions, although no SNP reached genome-wide significance after adjustment for genomic inflation. In AA, SNPs in or near the genes NUDT12, CHL1, GRIA1, CACNB2, and PYHIN1 were identified for SBP, and SNPs near ID3 were identified for DBP. For EA, promising SNPs for SBP were identified in ESR1and for DBP in SPATS2L and EYA2. Among these SNPs, none were common across phenotypes or population groups. Biologic plausibility exists for many of the identified genes, suggesting that these are candidate genes for regulation of BP and/or anti-hypertensive drug response. The lack of genome-wide significance is understandable in this small study employing gene-drug interactions. These findings provide a set of prioritized SNPs/candidate genes for future studies in HTN. Studies in more diversified population samples may help identify previously missed variants

Multi-omics insights into the biological mechanisms underlying statistical gene-by-lifestyle interactions with smoking and alcohol consumption

Though both genetic and lifestyle factors are known to influence cardiometabolic outcomes, less attention has been given to whether lifestyle exposures can alter the association between a genetic variant and these outcomes. The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium\u27s Gene-Lifestyle Interactions Working Group has recently published investigations of genome-wide gene-environment interactions in large multi-ancestry meta-analyses with a focus on cigarette smoking and alcohol consumption as lifestyle factors and blood pressure and serum lipids as outcomes. Further description of the biological mechanisms underlying these statistical interactions would represent a significant advance in our understanding of gene-environment interactions, yet accessing and harmonizing individual-level genetic and \u27omics data is challenging. Here, we demonstrate the coordinated use of summary-level data for gene-lifestyle interaction associations on up to 600,000 individuals, differential methylation data, and gene expression data for the characterization and prioritization of loci for future follow-up analyses. Using this approach, we identify 48 genes for which there are multiple sources of functional support for the identified gene-lifestyle interaction. We also identified five genes for which differential expression was observed by the same lifestyle factor for which a gene-lifestyle interaction was found. For instance, in gene-lifestyle interaction analysis, the T allele of rs6490056