Discovering Genetic Interactions in Large-Scale Association Studies by Stage-wise Likelihood Ratio Tests

Despite the success of genome-wide association studies in medical genetics, the underlying genetics of many complex diseases remains enigmatic. One plausible reason for this could be the failure to account for the presence of genetic interactions in current analyses. Exhaustive investigations of interactions are typically infeasible because the vast number of possible interactions impose hard statistical and computational challenges. There is, therefore, a need for computationally efficient methods that build on models appropriately capturing interaction. We introduce a new methodology where we augment the interaction hypothesis with a set of simpler hypotheses that are tested, in order of their complexity, against a saturated alternative hypothesis representing interaction. This sequential testing provides an efficient way to reduce the number of non-interacting variant pairs before the final interaction test. We devise two different methods, one that relies on a priori estimated numbers of marginally associated variants to correct for multiple tests, and a second that does this adaptively. We show that our methodology in general has an improved statistical power in comparison to seven other methods, and, using the idea of closed testing, that it controls the family-wise error rate. We apply our methodology to genetic data from the PRO-CARDIS coronary artery disease case/control cohort and discover three distinct interactions. While analyses on simulated data suggest that the statistical power may suffice for an exhaustive search of all variant pairs in ideal cases, we explore strategies for a priori selecting subsets of variant pairs to test. Our new methodology facilitates identification of new disease-relevant interactions from existing and future genome-wide association data, which may involve genes with previously unknown association to the disease. Moreover, it enables construction of interaction networks that provide a systems biology view of complex diseases, serving as a basis for more comprehensive understanding of disease pathophysiology and its clinical consequences.</p

Genetic loci on chromosome 5 are associated with circulating levels of interleukin-5 and eosinophil count in a European population with high risk for cardiovascular disease

IL-5 is a Th2 cytokine which activates eosinophils and is suggested to have an atheroprotective role. Genetic variants in the IL5 locus have been associated with increased risk of CAD and ischemic stroke. In this study we aimed to identify genetic variants associated with IL-5 concentrations and apply a Mendelian randomisation approach to assess IL-5 levels for causal effect on intima-media thickness in a European population at high risk of coronary artery disease. We analysed SNPs within robustly associated candidate loci for immune, inflammatory, metabolic and cardiovascular traits. We identified 2 genetic loci for IL-5 levels (chromosome 5, rs56183820, BETA = 0.11, P = 6.73E−5 and chromosome 14, rs4902762, BETA = 0.12, P = 5.76E−6) and one for eosinophil count (rs72797327, BETA = −0.10, P = 1.41E−6). Both chromosome 5 loci were in the vicinity of the IL5 gene, however the association with IL-5 levels failed to replicate in a meta-analysis of 2 independent cohorts (rs56183820, BETA = 0.04, P = 0.2763, I2 = 24, I2 − P = 0.2516). No significant associations were observed between SNPs associated with IL-5 levels or eosinophil count and IMT measures. Expression quantitative trait analyses indicate effects of the IL-5 and eosinophil-associated SNPs on RAD50 mRNA expression levels (rs12652920 (r2 = 0.93 with rs56183820) BETA = −0.10, P = 8.64E−6 and rs11739623 (r2 = 0.96 with rs72797327) BETA = −0.23, P = 1.74E−29, respectively). Our data do not support a role for IL-5 levels and eosinophil count in intima-media thickness, however SNPs associated with IL-5 and eosinophils might influence stability of the atherosclerotic plaque via modulation of RAD50 levels

Causal relevance of blood lipid fractions in the development of carotid atherosclerosis: Mendelian randomization analysis.

BACKGROUND: Carotid intima-media thickness (CIMT), a subclinical measure of atherosclerosis, is associated with risk of coronary heart disease events. Statins reduce progression of CIMT and coronary heart disease risk in proportion to the reduction in low-density lipoprotein cholesterol. However, interventions targeting triglycerides (TGs) or high-density lipoprotein cholesterol (HDL-C) have produced inconsistent effects on CIMT and coronary heart disease risk, making it uncertain whether such agents are ineffective for coronary heart disease prevention or whether CIMT is an inadequate marker of HDL-C or TG-mediated effects. We aimed to determine the causal association among the 3 major blood lipid fractions and common CIMT using mendelian randomization analysis. METHODS AND RESULTS: Genetic scores specific for low-density lipoprotein cholesterol, HDL-C, and TGs were derived based on single nucleotide polymorphisms from a gene-centric array in ≈5000 individuals (Cardiochip scores) and from a genome-wide association meta-analysis in >100 000 individuals (Global Lipids Genetic Consortium scores). These were used as instruments in a mendelian randomization analysis in 2 prospective cohort studies. A genetically predicted 1 mmol/L higher low-density lipoprotein cholesterol concentration was associated with a higher common CIMT by 0.03 mm (95% confidence interval, 0.01-0.04) and 0.04 mm (95% confidence interval, 0.02-0.06) based on the Cardiochip and Global Lipids Genetic Consortium scores, respectively. HDL-C and TGs were not causally associated with CIMT. CONCLUSIONS: Our findings confirm a causal relationship between low-density lipoprotein cholesterol and CIMT but not with HDL-C and TGs. At present, the suitability of CIMT as a surrogate marker in trials of cardiovascular therapies targeting HDL-C and TGs is questionable and requires further study

Genome-wide trans-ancestry meta-analysis provides insight into the genetic architecture of type 2 diabetes susceptibility.

To further understanding of the genetic basis of type 2 diabetes (T2D) susceptibility, we aggregated published meta-analyses of genome-wide association studies (GWAS), including 26,488 cases and 83,964 controls of European, east Asian, south Asian and Mexican and Mexican American ancestry. We observed a significant excess in the directional consistency of T2D risk alleles across ancestry groups, even at SNPs demonstrating only weak evidence of association. By following up the strongest signals of association from the trans-ethnic meta-analysis in an additional 21,491 cases and 55,647 controls of European ancestry, we identified seven new T2D susceptibility loci. Furthermore, we observed considerable improvements in the fine-mapping resolution of common variant association signals at several T2D susceptibility loci. These observations highlight the benefits of trans-ethnic GWAS for the discovery and characterization of complex trait loci and emphasize an exciting opportunity to extend insight into the genetic architecture and pathogenesis of human diseases across populations of diverse ancestry

Sex-Specific Effects of Adiponectin on Carotid Intima-Media Thickness and Incident Cardiovascular Disease

Background-Plasma adiponectin levels have previously been inversely associated with carotid intima-media thickness (IMT), a marker of subclinical atherosclerosis. In this study, we used a sex-stratified Mendelian randomization approach to investigate whether adiponectin has a causal protective influence on IMT. Methods and Results-Baseline plasma adiponectin concentrationwas tested for association with baseline IMT, IMT progression over 30 months, and occurrence of cardiovascular events within 3 years in 3430 participants (women, n=1777; men, n=1653) with high cardiovascular risk but no prevalent disease. Plasma adiponectin levels were inversely associated with baseline mean bifurcation IMT after adjustment for established risk factors (beta=-0.018, Pless than0.001) in men but not in women (beta=-0.006, P=0.185; P for interaction=0.061). Adiponectin levels were inversely associated with progression of mean common carotid IMT in men (beta=-0.0022, P=0.047), whereas no association was seen in women (0.0007, P=0.475; P for interaction=0.018). Moreover, we observed that adiponectin levels were inversely associated with coronary events in women (hazard ratio 0.57, 95% CI 0.37 to 0.87) but not in men (hazard ratio 0.82,95% CI0.54 to 1.25). Agenescore of adiponectin-raisingalleles in6loci, reported recently inalarge multi-ethnic metaanalysis, was inversely associated with baseline mean bifurcation IMT in men (beta=-0.0008, P=0.004) but not in women (beta=-0.0003, P=0.522; P for interaction=0.007). Conclusions-This report provides some evidence for adiponectin protecting against atherosclerosis, with effects being confined to men; however, compared with established cardiovascular risk factors, the effect of plasma adiponectin was modest. Further investigation involving mechanistic studies is warranted.Funding Agencies|European Commission [QLG1-CT-2002-00896]; Swedish Heart-Lung Foundation; Swedish Research Council [8691, 0593]; Knut and Alice Wallenberg Foundation; Foundation for Strategic Research; Stockholm County Council [592229]; Karolinska Institutet; Stockholm County Council; European Union Framework Programme 7 for the Innovative Medicine Initiative [IMI/115006]; Academy of Finland [110413]; British Heart Foundation [RG2008/08, RG2008/014]; Italian Ministry of Health (Ricerca Corrente); Uppsala University; Uppsala University Hospital; Swedish Research Council for Infrastructures; Swedish Heart-Lung Foundation [20120600, 20130399]; Tore Nilsson foundation; Gamla Tjanarinnor foundation; Thurings foundation; Stiftelsen for Gamla Tjanarinnor; Ake Wiberg foundation; Tore Nilssons foundation; Magnus Bergvall Foundation; Foundation for Old Servants; Ministry of Education and Culture in Finland; Vasterbotten County Council; Swedish Heart and Lung Foundation; National Excellence Program [TAMOP 4.2.4.A/1-11-1-2012-0001]; European Union; European Social Fund; UK Medical Research Council [K013351]; Economic and Social Research Council; Academy of Finland; University College London Genetics Institute</p

Identification of a novel proinsulin-associated SNP and demonstration that proinsulin is unlikely to be a causal factor in subclinical vascular remodelling using Mendelian randomisation

Background and aims Increased proinsulin relative to insulin levels have been associated with subclinical atherosclerosis (measured by carotid intima-media thickness (cIMT)) and are predictive of future cardiovascular disease (CVD), independently of established risk factors. The mechanisms linking proinsulin to atherosclerosis and CVD are unclear. A genome-wide meta-analysis has identified nine loci associated with circulating proinsulin levels. Using proinsulin-associated SNPs, we set out to use a Mendelian randomisation approach to test the hypothesis that proinsulin plays a causal role in subclinical vascular remodelling. Methods We studied the high CVD-risk IMPROVE cohort (n = 3345), which has detailed biochemical phenotyping and repeated, state-of-the-art, high-resolution carotid ultrasound examinations. Genotyping was performed using Illumina Cardio-Metabo and Immuno arrays, which include reported proinsulin-associated loci. Participants with type 2 diabetes (n = 904) were omitted from the analysis. Linear regression was used to identify proinsulin-associated genetic variants. Results We identified a proinsulin locus on chromosome 15 (rs8029765) and replicated it in data from 20,003 additional individuals. An 11-SNP score, including the previously identified and the chromosome 15 proinsulin-associated loci, was significantly and negatively associated with baseline IMTmean and IMTmax (the primary cIMT phenotypes) but not with progression measures. However, MR-Eggers refuted any significant effect of the proinsulin-associated 11-SNP score, and a non-pleiotropic SNP score of three variants (including rs8029765) demonstrated no effect on baseline or progression cIMT measures. Conclusions We identified a novel proinsulin-associated locus and demonstrated that whilst proinsulin levels are associated with cIMT measures, proinsulin per se is unlikely to have a causative effect on cIMT

Genetic fine mapping and genomic annotation defines causal mechanisms at type 2 diabetes susceptibility loci.

We performed fine mapping of 39 established type 2 diabetes (T2D) loci in 27,206 cases and 57,574 controls of European ancestry. We identified 49 distinct association signals at these loci, including five mapping in or near KCNQ1. 'Credible sets' of the variants most likely to drive each distinct signal mapped predominantly to noncoding sequence, implying that association with T2D is mediated through gene regulation. Credible set variants were enriched for overlap with FOXA2 chromatin immunoprecipitation binding sites in human islet and liver cells, including at MTNR1B, where fine mapping implicated rs10830963 as driving T2D association. We confirmed that the T2D risk allele for this SNP increases FOXA2-bound enhancer activity in islet- and liver-derived cells. We observed allele-specific differences in NEUROD1 binding in islet-derived cells, consistent with evidence that the T2D risk allele increases islet MTNR1B expression. Our study demonstrates how integration of genetic and genomic information can define molecular mechanisms through which variants underlying association signals exert their effects on disease

The genetics of blood pressure regulation and its target organs from association studies in 342,415 individuals

To dissect the genetic architecture of blood pressure and assess effects on target-organ damage, we analyzed 128,272 SNPs from targeted and genome-wide arrays in 201,529 individuals of European ancestry and genotypes from an additional 140,886 individuals were used for validation. We identified 66 blood pressure loci, of which 17 were novel and 15 harbored multiple distinct association signals. The 66 index SNPs were enriched for cis-regulatory elements, particularly in vascular endothelial cells, consistent with a primary role in blood pressure control through modulation of vascular tone across multiple tissues. The 66 index SNPs combined in a risk score showed comparable effects in 64,421 individuals of non-European descent. The 66-SNP blood pressure risk score was significantly associated with target-organ damage in multiple tissues, with minor effects in the kidney. Our findings expand current knowledge of blood pressure pathways and highlight tissues beyond the classic renal system in blood pressure regulation

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.Peer reviewe