PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Circ Res. Author manuscript; available in PMC 2017 September 22.
Published in final edited form as:
PMCID: PMC5609485
NIHMSID: NIHMS732221

Identifying Novel Gene Variants in Coronary Artery Disease and Shared Genes with Several Cardiovascular Risk Factors

Abstract

Rationale

Coronary Artery Disease (CAD) is a critical determinant of morbidity and mortality. Previous studies have identified several cardiovascular disease (CVD) risk factors, which may partly arise from a shared genetic basis with CAD, and thus be useful for discovery of CAD genes.

Objective

We aimed to improve discovery of CAD genes, and inform the etiologic relationship between CAD and several CVD risk factors using a shared polygenic signal-informed statistical framework.

Methods and Results

Using genome-wide association studies (GWAS) summary statistics and shared polygenic pleiotropy-informed conditional and conjunctional false discovery rate (FDR) methodology, we systematically investigated genetic overlap between CAD and 8 traits related to CVD risk factors: low density lipoprotein (LDL) cholesterol, high density lipoprotein (HDL) cholesterol, triglycerides (TG), type 2 diabetes (T2D), C-reactive protein (CRP), body mass index (BMI), systolic blood pressure (SBP) and type 1 diabetes (T1D). We found significant enrichment of single nucleotide polymorphisms (SNPs) associated with CAD as a function of their association with LDL, HDL, TG, T2D, CRP, BMI, SBP and T1D. Applying the conditional FDR method to the enriched phenotypes, we identified 67 novel loci associated with CAD (overall conditional FDR < 0.01). Further, we identified 53 loci with significant effects in both CAD and at least one of LDL, HDL, TG, T2D, CRP, SBP and T1D.

Conclusions

The observed polygenic overlap between CAD and cardio-metabolic risk factors indicates an etiological relation that warrants further investigation. The new gene loci identified implicate novel genetic mechanisms related to CAD.

Keywords: Genome-wide association study, genetic pleiotropy, comorbid disorder, lipids, cardiovascular genomics, DNA polymorphisms, genetic epidemiology
Subject terms: Chronic Ischemic Heart Disease, Genetics, Association Studies, Hypertension

INTRODUCTION

Coronary artery disease (CAD) is a leading cause of death worldwide. The development of CAD is influenced by both genetic and environmental factors, as evident by its high heritability (40–50%), shown in twin and family studies1. Genome-wide association studies (GWAS) in CAD have identified a total of 46 genetic variants reaching genome-wide significance for CAD2. However, the identified genetic variants explain only a small proportion of estimated heritability2, i.e. only a small amount of the familial clustering of CAD. This apparent paradox is widely seen across GWAS for complex traits and is termed “the missing heritability problem”3, 4. However, recent discoveries suggest that existing GWAS can capture more of the heritability due to common variants if proper statistical tools are used57.

Hypertension8, obesity9, abdominal fat10, diabetes11, dyslipidemia1214, inflammation as reflected by high levels of C-reactive protein (CRP)15, are associated with CAD. Several studies have found overlapping pathophysiology16, but the underlying shared genetic factors and the extent of the polygenic overlap across these phenotypes are mainly unknown. We have developed an analytical framework for complex traits building on the polygenic overlap17 between two or more phenotypes6. This method has the potential to capture more of the polygenic effects in complex traits18, and has successfully been applied to psychiatric6, cardiovascular19, neurological diseases20 and cancer21. This ‘shared polygenic signal’ method could be particularly informative in CAD, a disease with known co-morbidities and overlapping pathophysiology with related cardiovascular and metabolic disorders2, 2225.

We used this approach to leverage the power of multiple large genomic studies to describe the extent of the polygenic overlap and identify overlapping SNPs between CAD and 8 associated traits and cardiovascular disease (CVD) risk factors where recent GWAS results are available: low density lipoprotein (LDL) cholesterol26, high density lipoprotein (HDL) cholesterol26, triglycerides (TG)26, type 2 diabetes (T2D)27, CRP28, body mass index (BMI)29, systolic blood pressure (SBP)30, 31 and type 1 diabetes (T1D)32. By combining data from these different GWAS, we hypothesized that the shared polygenic signal approach can improve discovery of CAD genes, and inform the etiologic relationship between CAD and CVD risk factors.

METHODS

Participant samples

We obtained summary statistics from large-scale genomic studies (p-values and risk allele when available) from public access websites or through collaboration with investigators. The summary statistics are based on the Metabochip33 for CAD2 (n=194,427 including 63,746 cases) and T2D27 (n=149,830), and standard GWAS for LDL26 (n=95,454), HDL26 (n=99,900), TG26 (n=96,568), BMI29 (n=123,865), SBP31 (n=203,056) and T1D32 (n=16,559), and CRP (n=66,185)28. Details on the inclusion criteria and phenotype characteristics of the different GWAS are described in the original publications.

There were some overlapping controls between CAD and T2D and also between CAD and T1D. In both instances this was mainly due to the inclusion of one or more sub-studies employing a shared control design (e.g. used by the Wellcome Trust Case Control Consortium and deCODE Genetics)34 (see Online Table I). There was also some sample overlap between CAD and LDL, HDL, TG, BMI and SBP (Online Table I). Note that even without raw data, an upper bound for the amount of sample overlap is obtainable from the original publications by comparing the sub-study definitions and samples sizes for CAD and each secondary trait (correlation of uncorrected test statistics due to sample overlap is given in Online Table I; see LeBlanc et al (in prep) for details).

The Women’s Genome Health Study (WGHS), initiated in 1992, is an ongoing prospective cohort including 23,294 initially healthy North American women of European ancestry with whole genome genotype data and follow-up formajor incident health events, including myocardial infarction (MI) and coronary heart disease (CHD; composed of MI, CHD death, and coronary revascularization) are recorded35. Over the approximately, 20 years of follow-up, there were 387 and 1007 cases respectively of incident MI and CHD among the 23,294 women.

The relevant institutional review boards or ethics committees approved the research protocol of the individual GWAS and all participants provided written informed consent.

Statistical analyses

We use Matlab (version R2013a) for all statistical analysis unless otherwise indicated. First, we looked for evidence of overlapping polygenic signal for CAD and each secondary trait. In the absence of an overlapping polygenic signal, the expectation is that the p-value distribution for CAD is independent from the p-values in the secondary trait. The dependency of the p-value distribution for CAD on each secondary trait can be visually explored using conditional quantile-quantile plots to evaluate genetic ‘enrichment’ in CAD as a function of a secondary phenotype. Quantile-quantile plots are a descriptive tool for visualizing the difference between an observed distribution and a theoretical distribution. With GWAS, quantiles of the observed (nominal) p-values, denoted by ‘p’, are plotted on the y-axis, with the quantiles of the theoretical null distribution (i.e. the uniform distribution), here denoted by ‘q’, on the x-axis. Conventionally, the -log10 transform is used to emphasize tail areas. If there is no deviation from the null distribution and thus no true genetic association present, a quantile-quantile plot falls on the 1:1 line. Leftward deflections of the observed distribution from the null line reflect increased tail probabilities in the distribution of the test statistics, and consequently an over-abundance of low p-values compared to that expected by chance, termed ‘enrichment’. Here, we constructed conditional quantile-quantile plots to investigate if enrichment in the primary phenotype (CAD) is related to significance in a given secondary phenotype, as visualized by a leftward deflection from the null line on the conditional quantile-quantile plot. A conditional quantile-quantile plot was separately constructed for CAD and each of the 8 secondary traits. To test for statistical significance associated with these conditional quantile-quantile plots, we used the Anderson-Darling test21. In brief, this is a statistical test of whether a given sample of data is drawn from a given probability distribution and allows us to determine if an observed leftward deflection is statistically significant (for additional details see 21). In this case, we used set of SNPs (GWAS p>0.1 in the secondary trait), i.e., SNPs that are signal depleted in the secondary trait, as the comparison set.

Second, once statistically significant enrichment was confirmed, we computed conditional False Discovery Rates (FDR), a statistical framework that leverages shared polygenic signal6, 18, to improve the discovery of SNPs for the primary trait of interest, CAD. The standard FDR is designed to control the expected proportion of incorrectly rejected null hypotheses, and is employed to correct for multiple comparisons. An extension of the standard FDR is the conditional FDR6, which in our application, is used to incorporate information from GWAS summary statistics of a second phenotype. The conditional FDR is defined as the probability of a SNP being null in the first phenotype given that the p-values in the first and second phenotype are as small as or smaller than the observed ones (see Supplemental Methods). Importantly, ranking SNPs according to conditional FDR re-orders SNPs compared to their raw CAD p-values, and this new ranking favors SNPs showing signal in both CAD and the given secondary trait. In contrast, the standard FDR does not re-rank the SNPs compared to their raw CAD p-values, but instead suggests a different significance cut-off compared to the Bonferroni correction.

In additional analysis, we computed the conjunctional FDR18 to detect loci showing strong evidence of association with both CAD and the given secondary trait. Low values in conditional FDR can be driven by association with both phenotypes or with the primary phenotype alone, whereas low values in conjunctional FDR are driven by association with both phenotypes.

The application and interpretation of FDR-based methodology is more challenging for post-GWAS specialized SNP panels such as the Metabochip33. The standard FDR is widely applied in GWAS where any given SNP is assumed to have the same prior probability of association as all other SNPs. The Metabochip (~200,000 SNPs) is designed to follow up SNPs of interest relating to metabolic and cardiovascular traits, including fine mapping around genome-wide significant SNPs. As such, the true positives (and the false positives) come in large dependent clumps. Large-scale dependence in the signal can lead to biased FDR36. To correct for this bias, we used an LD-pruned set of SNPs to estimate the conditional FDR distribution, which was then used for estimating the conditional FDR for the full SNP set (see the Supplemental Methods for details of this estimation procedure). To visualize the conditional and conjunctional FDR, we constructed Manhattan plots. Detailed information on conditional quantile-quantile plots, Manhattan plots, as well as conditional and conjunctional FDR can be found in earlier reports6, 18 and/or in the Supplement.

The conditional FDR assumes independent samples for CAD and each of the secondary traits. However, several of the participants were included in both a secondary trait GWAS and in the CAD study. Partially overlapping subjects between studies leads to dependencies between the test statistics for different traits for a given SNP under the null hypothesis37. We estimated the expected correlation of the cross-trait GWAS test statistics under the null hypothesis of no genetic associations using a similar method to the one described for GWAS meta-analysis37, 38 and corrected for the estimated correlation due to shared subjects using the Mahalanobis transformation (LeBlanc et al in prep). These corrected test statistics were used in all further analysis.

Stratified replication rate

As an internal validation of stratified enrichment, we performed a stratified replication rate analysis using methods described previously,18 where the contributing studies of the CARDIoGRAMplusC4D Consortium were repeatedly divided into independent discovery and validation sets. The purpose of this analysis is to show that an observed pattern of stratified enrichment is not due to spurious effects. In brief, we randomly selected half of the studies (24) for the discovery set, and used the remaining studies for replication, and repeated this procedure 200 times. For each SNP in the replication set and the discovery set, we computed a meta-analysis test statistic (Liptak’s method). For the discovery set, we calculated the associated two-tailed p-values, whereas for the replication samples they were converted to one-tailed p-values in order to preserve the direction of effect in the discovery sample. We then created a vector of -log10(p-value) cutoffs and binned SNPs according to their p-values in the discovery set SNPs. For each bin, we kept track of their respective p-values in the replication set. We can then calculate the replication rate for each bin as defined by the proportion of SNPs in that bin which has a replication p value < 0.05. We checked for stratified replication rates by plotting the replication rate curves for four strata based on significance in each secondary trait, using the same strata definitions as for the conditional quantile-quantile plots.

Independent validation

For all novel CAD SNPs identified in the conditional FDR analysis, we checked for nominal replication (p<0.05) in the WGHS. Since the WGHS data is collected prospectively, we conducted age-adjusted Cox regression over approximately 20 yrs of follow-up ending in 2013 for both MI and CHD.

Expression quantitative trait loci (eQTL) annotation

We tested whether the novel CAD SNPs discovered in the current study are associated with genotype-dependent gene expression in various tissue types. Such SNPs are known as eQTLs. To this end, we cross-referenced our novel findings from the conditional FDR analysis with three cis-eQTL databases: in whole blood39 (the most powerful eQTL database available), adipose tissue40 (relevant for metabolic disease) and lymphoblastoid cells (LCL)40. The whole blood eQTL data has been collected in a large collaborative effort n=5311 samples, the adipose and LCL eQTLs are from a sample size of approximately n=850. We considered a SNP to be an eQTL using an FDR q-value cutoff of 0.05. The FDR q-values were already available for whole blood, while for adipose tissue and LCL we downloaded the publically-available eQTL data and calculated q-values using the qvalue() package available from Bioconductor (version 2.14) in R (version 3.1.1).

Biological pathway analysis

To better understand the biological context of our results, we conducted an Ingenuity Pathway Analysis (IPA, QIAGEN Redwood City, www.qiagen.com/ingenuity) including all previously reported CAD genes and the nearest annotated gene for each novel SNP reported in our study. The available molecules and/or relationships in the IPA Knowledge Base for mammal (humans, mouse or rat) were considered. We set the confidence filter to relationships where the confidence is experimentally observed. We allowed a maximum size of 35 genes for generating networks and we allowed up to 25 networks in the overall analysis. IPA computes a score for each network according to the fit of that network to a set of focus gene and p-values are calculated using the right-tailed Fisher’s exact test.

RESULTS

We used a two-step analysis strategy. First we assessed overlapping polygenic enrichment for CAD and each of the other traits via conditional quantile-quantile plots, and applied the Anderson-Darling test to define which of the 8 secondary traits show significant polygenic overlap. This test requires the direction of the association and as this information was unavailable for SBP, we relied on a visual inspection of the conditional quantile-quantile plot for SBP. As illustrated in Online Table II, all testable traits showed significant enrichment after Bonferroni correction for 21 tests and SBP showed strong visual evidence for enrichment. Therefore all 8 secondary traits were retained for the second step of the analysis. Second we applied conditional and conjunctional FDR methods to identify new CAD risk loci and to identify overlapping loci between CAD and each of the 8 associated traits. Overall FDR thresholds of 0.01 and 0.05 were chosen for conditional and conjunctional FDR respectively. Conservatively adjusting for the 8 secondary traits being considered21 this translated to thresholds of 0.01/8 and 0.05/8 for conditional and conjunctional FDR.

Conditional quantile-quantile plots for CAD conditioned on nominal p-values of association with LDL, CRP, T1D and T2D showed significant enrichment across different levels of significance (Figure 1). Similar significant enrichment patterns were seen for HDL, TG, SBP and BMI (see Online Figure I). The increasing leftward shift with more strictly defined strata based on nominal p-values of associated phenotypes suggests a greater proportion of true associations for a given nominal CAD p-value. This is indicative of cross-trait polygenic enrichment. As illustrated in Figure 1, panel A: LDL, the proportion of SNPs in the −log10(pLDL) ≥ 3 category reaching a given significance level (e.g., −log10(pCAD) > 6) is much greater than for the all SNPs category, indicating a high level of enrichment (Figure 1).

Figure 1
Shared Polygenic Enrichment

Stratified replication rates were observed for all secondary traits with the exception of BMI (Online Figure II), indicating that the observed enrichment in the conditional quantile-quantile plots is also associated with increased replication rates. The observed pattern of stratified enrichment does not result from spurious effects, and replication rate is increased by conditioning on significance in each of the secondary traits, with the possible exception of BMI.

Conditional and conjunctional FDR were calculated for CAD paired with each of the 8 secondary phenotypes showing enrichment. The results of each analysis were filtered as follows. First, we filtered the lists of significant SNPs by their linkage disequilibrium patterns as observed in the 1000 Genomes41 dataset and report only the most significant result per annotated gene. We considered a SNP to be an independent finding if the linkage disequilibrium, defined using r2, was less than 0.2 with all other SNPs in the filtered list. Second, we further filtered the list of significant SNPs for novelty with respect to previously published CAD SNPs. We filtered out any previously reported genes and SNPs, including SNPs in linkage disequilibrium (r2>0.2) with those previously reported SNPs. Thus, the list of significant SNPs presented in Table 1 represent, to the best of our knowledge, independent novel SNPs for CAD. The corresponding conditional Manhattan plot is given in Figure 2.

Figure 2
Conditional FDR Manhattan plot’ of −log10 (FDR)* values for Coronary Artery Disease (CAD) alone (black), and −log10 (conditional FDR) for CAD given type 2 diabetes (T2D; CAD|T2D; navy blue), CAD given type 1 diabetes (T1D; ...
Table 1
Conditional FDR (<0.01), after controlling for multiple testing across phenotypes.

Over all 8 secondary traits, we identify 101 SNPs associated with CAD, 67 of which have not previously associated with CAD (previously reported SNPs not shown). Many of these new loci are located in regions with borderline significant association with CAD in previous studies42 as is evident by the CAD association p-value column given in Table 1. Of interest, several of the identified loci are found across the conditional analysis from several risk factors. These loci are not found using standard methods applying a genome-wide Bonferroni correction.

We looked to the WGHS for independent validation of these 67 new CAD SNPs and 12 of these show nominal replication for at least one endpoint (CHD or MI); see Online Table III.

Of the 67 novel CAD loci, 32 show genotype-dependent gene expression in whole blood regulating the expression of 57 unique genes and 42 of these 67 SNPs would not have been detected using the standard (unconditioned) FDR. We found evidence for 16 and 18 loci having an eQTL effect in adipose tissue and LCL respectively (Table 2). For six of these loci we observed an eQTL effect on the same gene in both whole blood and adipose tissue. Interestingly, 18 loci show an effect on the gene expression of more than one gene.

Table 2
Novel CAD SNPs that are also eQTLs in blood and/or adipose tissue and/or lymphoblastoid cells.

To further evaluate genetic overlap, we used the conjunctional FDR to identify SNPs with significant effects in both CAD and its associated risk factors. The conjunctional Manhattan plot for CAD is shown in Online Figure II. We identified 53 loci achieving conjunctional FDR<0.05, after adjustment for using multiple risk factors and pruning the results in the same manner as for the conditional FDR (Online Table IV; corresponding z-scores in Online Table V).

Follow-up Ingenuity Pathways Analysis (IPA) identified highly significantly associated “Top Canonical Pathways” relevant to CAD (e.g. LXR (liver X receptor)/RXR (retinoid X receptor) as well as FXR (Farnesoid X Receptor)/RXR Activation and Atherosclerosis Signaling); (Online Table V). Additionally, in “Top Diseases and Bio Function” CAD relevant diseases and functions are on top (Cardiovascular Disease and Lipid Metabolism) in the subgroups “Diseases and Disorders” and “Molecular and Cellular Functions”.

DISCUSSION

Combining data from large-scale genomic studies from different phenotypes in a conditional FDR framework, we show polygenic overlap between CAD and several CVD risk factor phenotypes and identify 67 novel CAD susceptibility loci. Further, conjunctional FDR analysis identified 53 novel loci associated with both CAD and the CVD risk factors LDL, HDL, TG, T1D, T2D, CRP and SBP. Importantly, we validated the conditional FDR approach by showing that replication rates in independent CAD sub-studies increase as a function of p-value in each secondary trait, with the possible exception of BMI. Further, we see nominal replication for 12/67 SNPs in the WGHS. Overall, these results suggest that a proportion of the clinically and epidemiologically observed association between these phenotypes can be explained by overlapping genetic loci (pleiotropy) and not simply shared environmental risk factors. Further, the findings provide further evidence that CAD is a highly polygenic disease.

Our findings of polygenic overlap provide novel insights into the relationship between CAD and major CVD risk factors. We demonstrate an interesting genetic dissociation among these risk factors and CAD, with strong enrichment for lipids, inflammation and metabolic disorders. The combination of dyslipidemia (i.e., high TG and LDL cholesterol and low HDL cholesterol), T2D, and high blood pressure forms the metabolic syndrome1214, 43, 44, and all of these factors (particularly LDL) showed strong genetic overlap with CAD. This is in agreement with a recent reports suggesting a common genetic basis for regulation of lipid and glucose homeostasis45, while previous studies did not show common genes for the different components of the metabolic syndrome46, but revealed strong lipid gene contribution. It is further supported by the pathway analysis that identified “Atherosclerosis Signaling” and “FXR/RXR Activation” among the three most relevant pathways. Genes activated by the FXR has been shown to influence vascular tension and regulate the unloading of cholesterol from foam cells47. Another important finding is the overlap between CAD and T2D. Based on conditional analysis of these two phenotypes, 21 novel loci were identified. This is in line with previous single gene studies suggesting a genetic link between T2D and CAD48.

The strong shared polygenic signal between LDL and CAD emphasizes the important role of LDL in CAD development, and support the notion that risk genes for atherosclerosis, such as LDL genes, are causal for CAD as recently suggested49. Finally, two of the phenotypes most strongly overlapping with CAD were CRP and T1D, two immune related phenotypes. CRP is regarded as a reliable marker of systemic inflammation and its role as a biomarker in CAD has been attributed to its ability to reflect up-stream inflammatory pathways. However, the finding in the present study suggests that the link between CRP and CAD may also reflect overlapping genetic loci. T1D is related to auto-immune mechanisms and its genetic overlap with CAD underlines the important role of the immune system in CAD, and could be due to a large number of overlapping genes between immune and lipid phenotypes50. In fact, the bidirectional interaction between inflammation and lipids is regarded as a phenotypic hallmark of atherosclerosis, and our findings suggest that this phenotype could reflect overlapping genes between these two interacting pathophysiological arms of atherogenesis. The pathway analysis revealed “LXR/RXR Activation”, as the top ranked canonical pathway. LXR/RXR are heterodimer nuclear receptors/transcription factors. LXR acts as a cholesterol sensor, and LXR pathway activation has been shown to stimulate lipogenesis and hypertriglyceridemia51. LXR/RXR can also modulate inflammatory responses to cholesterol exposure and could represent a regulator of the interaction between lipids and inflammation, being the most important pathway in the pathogenesis of CAD.

In the original CAD GWAS and follow-up Metabochip study, 46 loci were identified2. By combining the original CAD results with the CVD risk factor phenotypes GWAS, we identified 101 significant loci associated with CAD, of which 67 are novel, using the conditional FDR approach. Even though the original CAD study was quite large2, the increased power provided by additional GWAS of associated phenotypes together with the conditional FDR method more than doubled gene discovery. The novel SNPs discovered here contribute to explaining more of the missing heritability for CAD, but we cannot quantify how much more is explained since we are working at the summary-statistic level. These findings underline the cost-effectiveness of the current statistical methods and highlight several interesting genes in CAD pathology. IL1F10 (interleukin 1 family, member 10 (theta)) was identified in the pathway analysis of the CAD GWAS, it is known to bind IL1R and stimulate NF-kB pathway. VEGFA (vascular endothelial growth factor A) is well known in the CVD field, but to the best of our knowledge, this has never been shown in genetic studies. SLC18A1 (solute carrier family 18 (vesicular monoamine transporter), member 1) has been implicated in neuropsychiatric disorders, but not previously in CVD. SERPINH1 (serpin peptidase inhibitor, clade H (heat shock protein 47), member 1, (collagen binding protein 1)) is a heat shock protein, known to be involved in atherosclerosis. ILF3 (interleukin enhancer binding factor 3, 90kDa) is a matrix metallo proteinase, well studied in the CVD research field, and the findings of ILF10 and ILF3 underscore the role of the IL-1 cytokine family in CAD.

Although nearest-gene annotation can be informative, the vast majority of discovered SNPs are located outside coded DNA regions52. Therefore, annotating the identified genetic variants to the correct causal genes for the phenotype of interest often remains challenging52. One of the potential mechanisms whereby SNPs may affect phenotype variations is through altered gene expression. We successfully identified eQTL effects in whole blood, LCL and adipose tissue, suggesting these genes as potential causal candidates. Of interest, some of the genetic variants showed an effect on the gene expression of more than one gene. We speculate that the shared effect of the genetic variants on the phenotypes under study might be explained by the regulation of several different genes, but further studies would be necessary connect the genes with altered gene expression seen in Table 2 to the clinical phenotypes. Moreover, the majority of the genes regulated by the genetic variants were different from the nearest annotated gene. Given that the original whole blood has markedly different power and used different statistical eQTL definitions than the LCL and adipose tissue eQTL studies, a detailed cross-tissue comparison is not possible. Further studies are needed to determine the functional mechanisms involved in the novel CAD loci identified here.

There are certain limitations associated with the present results. Due to the overlap in some of the GWAS samples examined, we cannot completely exclude the contribution from environmental or behavioral factors. The shared participants between genomic studies could also affect the findings. However, we did adjust for overlapping subjects, and used strict FDR thresholds to account for the 8 secondary traits. Although clinical comorbidity and shared pathophysiology between these phenotypes poses a challenge for the interpretation of the basis of the shared polygenetic signals, their utility for increasing the power to detect new loci for CAD is not affected. The question remains if the identified shared genes are independent of other phenotypes (biological pleiotropy), or that the current findings are results of overlapping phenotypes (mediated by other phenotypes), as several of these risk factors can be co-occurring (mediated pleiotropy)53. However, it appears reasonable to interpret our findings as reflecting the existence of shared genetically determined pathophysiological processes across CAD and the associated phenotypes. In general, FDR methodology is a less conservative approach to multiple testing than Bonferroni correction. However, by using the conditional FDR, we are not simply relaxing the significance threshold, but are increasing power and incorporating useful information from a second trait into the analysis, allowing us to identify the SNPs more likely to replicate. We have not strictly replicated all of these findings in independent samples but we have shown that replication rates increase by conditioning on significance in the secondary traits, and have shown that 12 SNPs nominally replicate in the WGS. While the prospective design of the WGHS makes it suitable for validation of the candidate CAD associations, the numbers of incident events of MI and CAD were much smaller than in the discovery sample, which was composed of a preponderance of men compared with the all-female composition of the WGHS. However, in spite of much lower power and possibility of differences according to sex, the WGHS is the largest and most relevant independent dataset we were able to access and we found nominal association for novel CAD 12 loci.

In conclusion, we found substantial polygenic overlap between CAD and several related conditions, importantly LDL, T2D and CRP, providing more evidence for fundamental etiological relationship between these phenotypes that cannot be explained by lifestyle factors. The 67 novel CAD loci identified here provide new insight into genetic mechanisms of CAD and may form the basis for earlier diagnosis and new prevention and treatment strategies.

Novelty and Significance

What Is Known?

  • Previous work has identified 46 genetic risk variants associated with coronary artery disease (CAD)
  • Genetic data for traits with overlapping pathophysiology can be combined to improve power for identifying novel genetic risk variants.

What New Information Does This Article Contribute?

  • We identified 67 new genetic risk variants for CAD.
  • CAD and several cardio-metabolic traits share a large number of genetic risk factors.

Clinical and epidemiological evidence suggests a relationship between CAD and cardio-metabolic traits. In the presence of a shared polygenic signal (i.e., a large number of shared risk variants each with a small effect), traits with overlapping pathophysiology with CAD can be used in combination with novel statistical methodology to improve discovery of variants associated with CAD. Using large-scale genetic data from CAD and genetic data from hypertension, obesity, abdominal fat, diabetes, dyslipidemia, and inflammation (C-reactive protein), we found a polygenic overlap between CAD and each of these related traits. We identified 67 novel CAD risk variants and 53 risk variants jointly associated with CAD and at least one other related trait. These results highlight the importance of shared polygenic risk factors between coronary artery disease and cardiovascular risk factors. Our findings provide important insights into molecular mechanisms underlying coronary artery disease and have potential implications for prevention and treatment strategies.

Supplementary Material

306629R2 Online Data Supplement

Acknowledgments

The authors would like to thank the DIAGRAM Consortium, the CARDIoGRAMplusC4D Consortium, the Global Lipids Genetics Consortium, The International Consortium for Blood Pressure GWAS, the GEFOS Consortium, the Giant Consortium, the Type 1 Diabetes Consortium and the CHARGE Consortium Inflammation Working Group for the summary statistics GWAS data.

SOURCES OF FUNDING

The present project was supported by the Research Council of Norway (213837, 223273); the South East Norway Health Authority (2013–123); the Kristian Gerhard Jebsen Foundation; in addition to the National Institutes of Health grants R01AG031224, R01EB000790, and RC2DA29475, and the European Union project tegrat@senegVC; a grant from the Leducq Foundation, CADgenomics; and BMBF, Forschungskonsortien zur Systemmedizin, e:AtheroSysMED; DFG, SFB 1123. The WGHS is supported by HL043851, HL080467, and HL099355 from the National Heart, Lung, and Blood Institute and CA047988 from the National Cancer Institute with collaborative scientific support and funding for genotyping provided by Amgen. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Nonstandard Abbreviations and Acronyms

BMI
body mass index
CAD
coronary artery disease
CRP
C-reactive protein
CVD
cardiovascular disease
eQTL
expression quantitative trait locus
FDR
false discovery rate
GWAS
genome-wide association study
HDL
high density lipoprotein
LCL
lymphoblastoid cells
LDL
low density lipoprotein
SBP
systolic blood pressure
SNP
single nucleotide polymorphism
T1D
type 1 diabetes
T2D
type 2 diabetes
TG
triglycerides

Footnotes

This manuscript was sent to Elizabeth M. McNally, Consulting Editor, for review by expert referees, editorial decision, and final disposition.

DISCLOSURES

None reported.

References

1. Peden JF, Farrall M. Thirty-five common variants for coronary artery disease: The fruits of much collaborative labour. Hum Mol Genet. 2011;20:R198–205. [PMC free article] [PubMed]
2. CAD Consortium. Deloukas P, Kanoni S, Willenborg C, Farrall M, Assimes TL, Thompson JR, Ingelsson E, Saleheen D, Erdmann J, Goldstein BA, Stirrups K, Konig IR, et al. Large-scale association analysis identifies new risk loci for coronary artery disease. Nat Genet. 2013;45:25–33. [PMC free article] [PubMed]
3. Eichler EE, Flint J, Gibson G, Kong A, Leal SM, Moore JH, Nadeau JH. Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010;11:446–450. [PMC free article] [PubMed]
4. Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, McCarthy MI, Ramos EM, Cardon LR, Chakravarti A, Cho JH, Guttmacher AE, Kong A, et al. Finding the missing heritability of complex diseases. Nature. 2009;461:747–753. [PMC free article] [PubMed]
5. Li C, Yang C, Gelernter J, Zhao H. Improving genetic risk prediction by leveraging pleiotropy. Hum Genet. 2014;133:639–650. [PMC free article] [PubMed]
6. Andreassen OA, Thompson WK, Schork AJ, Ripke S, Mattingsdal M, Kelsoe JR, Kendler KS, O’Donovan MC, Rujescu D, Werge T, Sklar P, Roddey JC, Chen C-H, et al. Improved detection of common variants associated with schizophrenia and bipolar disorder using pleiotropy-informed conditional false discovery rate. PLoS Genet. 2013;9:e1003455. [PMC free article] [PubMed]
7. Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham JM, de Andrade M, Feenstra B, Feingold E, Hayes MG, Hill WG, Landi MT, Alonso A, et al. Genome partitioning of genetic variation for complex traits using common snps. Nat Genet. 2011;43:519–525. [PMC free article] [PubMed]
8. van den Hoogen PC, Feskens EJ, Nagelkerke NJ, Menotti A, Nissinen A, Kromhout D. The relation between blood pressure and mortality due to coronary heart disease among men in different parts of the world. Seven countries study research group. N Engl J Med. 2000;342:1–8. [PubMed]
9. The Global Burden of Metabolic Risk Factors for Chronic Diseases Collaboration. Metabolic mediators of the effects of body-mass index, overweight, and obesity on coronary heart disease and stroke: A pooled analysis of 97 prospective cohorts with 1.8 million participants. Lancet. 2014 [PMC free article] [PubMed]
10. Romero-Corral A, Montori VM, Somers VK, Korinek J, Thomas RJ, Allison TG, Mookadam F, Lopez-Jimenez F. Association of bodyweight with total mortality and with cardiovascular events in coronary artery disease: A systematic review of cohort studies. Lancet. 2006;368:666–678. [PubMed]
11. Group AS. Cushman WC, Evans GW, Byington RP, Goff DC, Jr, Grimm RH, Jr, Cutler JA, Simons-Morton DG, Basile JN, Corson MA, Probstfield JL, Katz L, Peterson KA, et al. Effects of intensive blood-pressure control in type 2 diabetes mellitus. N Engl J Med. 2010;362:1575–1585. [PMC free article] [PubMed]
12. Libby P, Theroux P. Pathophysiology of coronary artery disease. Circulation. 2005;111:3481–3488. [PubMed]
13. Messerli FH, Williams B, Ritz E. Essential hypertension. Lancet. 2007;370:591–603. [PubMed]
14. Eckel RH, Grundy SM, Zimmet PZ. The metabolic syndrome. Lancet. 2005;365:1415–1428. [PubMed]
15. Lagrand WK, Visser CA, Hermens WT, Niessen HW, Verheugt FW, Wolbink GJ, Hack CE. C-reactive protein as a cardiovascular risk factor: More than an epiphenomenon? Circulation. 1999;100:96–102. [PubMed]
16. Coffman TM. Under pressure: The search for the essential mechanisms of hypertension. Nat Med. 2011;17:1402–1409. [PubMed]
17. Wagner GP, Zhang J. The pleiotropic structure of the genotype-phenotype map: The evolvability of complex organisms. Nat Rev Genet. 2011;12:204–213. [PubMed]
18. Andreassen OA, Djurovic S, Thompson WK, Schork AJ, Kendler KS, O’Donovan MC, Rujescu D, Werge T, van de Bunt M, Morris AP, McCarthy MI, International Consortium for Blood Pressure Group, Diabetes Genetics Replicaton and Meta-analysis Consortium, Psychiatric Genomics Consortium Schizophrenia Working Group et al. Improved detection of common variants associated with schizophrenia by leveraging pleiotropy with cardiovascular-disease risk factors. Am J Hum Genet. 2013;92:197–209. [PubMed]
19. Andreassen OA, McEvoy LK, Thompson WK, Wang Y, Reppe S, Schork AJ, Zuber V, Barrett-Connor E, Gautvik K, Aukrust P, Karlsen TH, Djurovic S, Desikan RS, et al. Identifying common genetic variants in blood pressure due to polygenic pleiotropy with associated phenotypes. Hypertension. 2014;63:819–826. [PMC free article] [PubMed]
20. Andreassen OA, Harbo HF, Wang Y, Thompson WK, Schork AJ, Mattingsdal M, Zuber V, Bettella F, Ripke S, Kelsoe JR, Kendler KS, O’Donovan MC, Sklar P, et al. Genetic pleiotropy between multiple sclerosis and schizophrenia but not bipolar disorder: Differential involvement of immune-related gene loci. Mol Psychiatry. 2014 [PMC free article] [PubMed]
21. Andreassen OA, Zuber V, Thompson WK, Schork AJ, Bettella F, Djurovic S, Desikan RS, Mills IG, Dale AM. Shared common variants in prostate cancer and blood lipids. Int J Epidemiol. 2014 [PMC free article] [PubMed]
22. Wilson PW, D’Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB. Prediction of coronary heart disease using risk factor categories. Circulation. 1998;97:1837–1847. [PubMed]
23. Haffner SM, Lehto S, Ronnemaa T, Pyorala K, Laakso M. Mortality from coronary heart disease in subjects with type 2 diabetes and in nondiabetic subjects with and without prior myocardial infarction. N Engl J Med. 1998;339:229–234. [PubMed]
24. Malik S, Wong ND, Franklin SS, Kamath TV, L’Italien GJ, Pio JR, Williams GR. Impact of the metabolic syndrome on mortality from coronary heart disease, cardiovascular disease, and all causes in united states adults. Circulation. 2004;110:1245–1250. [PubMed]
25. Hansson GK. Inflammation, atherosclerosis, and coronary artery disease. N Eng J Med. 2005;352:1685–1695. [PubMed]
26. Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, Pirruccello JP, Ripatti S, Chasman DI, Willer CJ, Johansen CT, Fouchier SW, Isaacs A, et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010;466:707–713. [PMC free article] [PubMed]
27. Morris AP, Voight BF, Teslovich TM, Ferreira T, Segre AV, Steinthorsdottir V, Strawbridge RJ, Khan H, Grallert H, Mahajan A, Prokopenko I, Kang HM, Dina C, Esko T, et al. Large-scale association analysis provides insights into the genetic architecture and pathophysiology of type 2 diabetes. Nat Genet. 2012;44:981–990. [PMC free article] [PubMed]
28. Dehghan A, Dupuis J, Barbalic M, Bis JC, Eiriksdottir G, Lu C, Pellikka N, Wallaschofski H, Kettunen J, Henneman P, Baumert J, Strachan DP, Fuchsberger C, et al. Meta-analysis of genome-wide association studies in >80 000 subjects identifies multiple loci for c-reactive protein levels. Circulation. 2011;123:731–738. [PMC free article] [PubMed]
29. Speliotes EK, Willer CJ, Berndt SI, Monda KL, Thorleifsson G, Jackson AU, Allen HL, Lindgren CM, Luan J, Magi R, Randall JC, Vedantam S, Winkler TW, et al. Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet. 2010;42:937–948. [PMC free article] [PubMed]
30. Heid IM, Jackson AU, Randall JC, Winkler TW, Qi L, Steinthorsdottir V, Thorleifsson G, Zillikens MC, Speliotes EK, Magi R, Workalemahu T, White CC, Bouatia-Naji N, et al. Meta-analysis identifies 13 new loci associated with waist-hip ratio and reveals sexual dimorphism in the genetic basis of fat distribution. Nat Genet. 2010;42:949–960. [PMC free article] [PubMed]
31. Ehret GB, Munroe PB, Rice KM, Bochud M, Johnson AD, Chasman DI, Smith AV, Tobin MD, Verwoert GC, Hwang SJ, Pihur V, Vollenweider P, O’Reilly PF, et al. Genetic variants in novel pathways influence blood pressure and cardiovascular disease risk. Nature. 2011;478:103–109. [PMC free article] [PubMed]
32. Barrett JC, Clayton DG, Concannon P, Akolkar B, Cooper JD, Erlich HA, Julier C, Morahan G, Nerup J, Nierras C, Plagnol V, Pociot F, Schuilenburg H, et al. Genome-wide association study and meta-analysis find that over 40 loci affect risk of type 1 diabetes. Nat Genet. 2009;41:703–707. [PMC free article] [PubMed]
33. Voight BF, Kang HM, Ding J, Palmer CD, Sidore C, Chines PS, Burtt NP, Fuchsberger C, Li Y, Erdmann J, Frayling TM, Heid IM, Jackson AU, et al. The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits. PLoS Genet. 2012;8:e1002793. [PMC free article] [PubMed]
34. Witte JS. Genome-wide association studies and beyond. Annu Rev Public Health. 2010;31:9–20. 24 p following 20. [PMC free article] [PubMed]
35. Ridker PM, Chasman DI, Zee RY, Parker A, Rose L, Cook NR, Buring JE. Rationale, design, and methodology of the women’s genome health study: A genome-wide association study of more than 25,000 initially healthy american women. Clin Chem. 2008;54:249–255. [PubMed]
36. Chen L, Storey JD. Relaxed significance criteria for linkage analysis. Genetics. 2006;173:2371–2381. [PubMed]
37. Lin DY, Sullivan PF. Meta-analysis of genome-wide association studies with overlapping subjects. Am J Hum Genet. 2009;85:862–872. [PubMed]
38. Zaykin DV, Kozbur DO. P-value based analysis for shared controls design in genome-wide association studies. Genet Epidemiol. 2010;34:725–738. [PMC free article] [PubMed]
39. Westra HJ, Peters MJ, Esko T, Yaghootkar H, Schurmann C, Kettunen J, Christiansen MW, Fairfax BP, Schramm K, Powell JE, Zhernakova A, Zhernakova DV, Veldink JH, et al. Systematic identification of trans eqtls as putative drivers of known disease associations. Nat Genet. 2013;45:1238–1243. [PMC free article] [PubMed]
40. Grundberg E, Meduri E, Sandling JK, Hedman AK, Keildson S, Buil A, Busche S, Yuan W, Nisbet J, Sekowska M, Wilk A, Barrett A, Small KS, et al. Global analysis of DNA methylation variation in adipose tissue from twins reveals links to disease-associated variants in distal regulatory elements. Am J Hum Genet. 2013;93:876–890. [PubMed]
41. Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, Hurles ME, McVean GA. A map of human genome variation from population-scale sequencing. Nature. 2010;467:1061–1073. [PMC free article] [PubMed]
42. Ehret GB. Genome-wide association studies: Contribution of genomics to understanding blood pressure and essential hypertension. Curr Hypertens Rep. 2010;12:17–25. [PMC free article] [PubMed]
43. D’Agostino RB, Sr, Vasan RS, Pencina MJ, Wolf PA, Cobain M, Massaro JM, Kannel WB. General cardiovascular risk profile for use in primary care: The framingham heart study. Circulation. 2008;117:743–753. [PubMed]
44. Conroy RM, Pyorala K, Fitzgerald AP, Sans S, Menotti A, De Backer G, De Bacquer D, Ducimetiere P, Jousilahti P, Keil U, Njolstad I, Oganov RG, Thomsen T, et al. Estimation of ten-year risk of fatal cardiovascular disease in europe: The score project. Eur Heart J. 2003;24:987–1003. [PubMed]
45. Albert JS, Yerges-Armstrong LM, Horenstein RB, Pollin TI, Sreenivasan UT, Chai S, Blaner WS, Snitker S, O’Connell JR, Gong DW, Breyer RJ, 3rd, Ryan AS, McLenithan JC, Shuldiner AR, et al. Null mutation in hormone-sensitive lipase gene and risk of type 2 diabetes. N Eng J Med. 2014;370:2307–2315. [PMC free article] [PubMed]
46. Kristiansson K, Perola M, Tikkanen E, Kettunen J, Surakka I, Havulinna AS, Stancakova A, Barnes C, Widen E, Kajantie E, Eriksson JG, Viikari J, Kahonen M, Lehtimaki T, et al. Genome-wide screen for metabolic syndrome susceptibility loci reveals strong lipid gene contribution but no evidence for common genetic basis for clustering of metabolic syndrome traits. Circ Cardiovasc Genet. 2012;5:242–249. [PMC free article] [PubMed]
47. Porez G, Prawitt J, Gross B, Staels B. Bile acid receptors as targets for the treatment of dyslipidemia and cardiovascular disease. J Lipid Res. 2012;53:1723–1737. [PMC free article] [PubMed]
48. Qi L, Qi Q, Prudente S, Mendonca C, Andreozzi F, di Pietro N, Sturma M, Novelli V, Mannino GC, Formoso G, Gervino EV, Hauser TH, Muehlschlegel JD, et al. Association between a genetic variant related to glutamic acid metabolism and coronary heart disease in individuals with type 2 diabetes. JAMA. 2013;310:821–828. [PMC free article] [PubMed]
49. Do R, Willer CJ, Schmidt EM, Sengupta S, Gao C, Peloso GM, Gustafsson S, Kanoni S, Ganna A, Chen J, Buchkovich ML, Mora S, Beckmann JS, et al. Common variants associated with plasma triglycerides and risk for coronary artery disease. Nat Genet. 2013;45:1345–1352. [PMC free article] [PubMed]
50. Andreassen OA, Desikan RS, Wang Y, Thompson WK, Schork AJ, Zuber V, Doncheva NT, Ellinghaus E, Albrecht M, Mattingsdal M, Franke A, Lie BA, Mills I, et al. Abundant genetic overlap between blood lipids and immune-mediated diseases indicates shared molecular genetic mechanisms. PLoS One. 2015;10:e0123057. [PMC free article] [PubMed]
51. Parikh M, Patel K, Soni S, Gandhi T. Liver x receptor: A cardinal target for atherosclerosis and beyond. J Atheroscler Thromb. 2014;21:519–531. [PubMed]
52. Hindorff LA, Sethupathy P, Junkins HA, Ramos EM, Mehta JP, Collins FS, Manolio TA. Potential etiologic and functional implications of genome-wide association loci for human diseases and traits. Proc Natl Acad Sci U S A. 2009;106:9362–9367. [PubMed]
53. Solovieff N, Cotsapas C, Lee PH, Purcell SM, Smoller JW. Pleiotropy in complex traits: Challenges and strategies. Nat Rev Genet. 2013;14:483–495. [PMC free article] [PubMed]