Article Text

PDF

Genome-wide association study of perioperative myocardial infarction after coronary artery bypass surgery
  1. Miklos D Kertai1,
  2. Yi-Ju Li2,3,
  3. Yen-Wei Li2,
  4. Yunqi Ji2,
  5. John Alexander4,5,
  6. Mark F Newman1,5,
  7. Peter K Smith6,
  8. Diane Joseph5,
  9. Joseph P Mathew1,
  10. Mihai V Podgoreanu1,5
  11. for the Duke Perioperative Genetics and Safety Outcomes (PEGASUS) Investigative Team
  1. 1Division of Cardiothoracic Anesthesiology, Duke University, Durham, North Carolina, USA
  2. 2Department of Biostatistics and Bioinformatics, Duke University, Durham, North Carolina, USA
  3. 3Duke Molecular Physiology Institute; Duke University, Durham, North Carolina, USA
  4. 4Division of Cardiology, Duke University, Durham, North Carolina, USA
  5. 5Duke Clinical Research Institute; Duke University, Durham, North Carolina, USA
  6. 6Cardiac Surgery; Duke University, Durham, North Carolina, USA
  1. Correspondence to Dr Mihai V Podgoreanu; mihai.podgoreanu{at}duke.edu

Abstract

Objectives Identification of patient subpopulations susceptible to develop myocardial infarction (MI) or, conversely, those displaying either intrinsic cardioprotective phenotypes or highly responsive to protective interventions remain high-priority knowledge gaps. We sought to identify novel common genetic variants associated with perioperative MI in patients undergoing coronary artery bypass grafting using genome-wide association methodology.

Setting 107 secondary and tertiary cardiac surgery centres across the USA.

Participants We conducted a stage I genome-wide association study (GWAS) in 1433 ethnically diverse patients of both genders (112 cases/1321 controls) from the Genetics of Myocardial Adverse Outcomes and Graft Failure (GeneMAGIC) study, and a stage II analysis in an expanded population of 2055 patients (225 cases/1830 controls) combined from the GeneMAGIC and Duke Perioperative Genetics and Safety Outcomes (PEGASUS) studies. Patients undergoing primary non-emergent coronary bypass grafting were included.

Primary and secondary outcome measures The primary outcome variable was perioperative MI, defined as creatine kinase MB isoenzyme (CK-MB) values ≥10× upper limit of normal during the first postoperative day, and not attributable to preoperative MI. Secondary outcomes included postoperative CK-MB as a quantitative trait, or a dichotomised phenotype based on extreme quartiles of the CK-MB distribution.

Results Following quality control and adjustment for clinical covariates, we identified 521 single nucleotide polymorphisms in the stage I GWAS analysis. Among these, 8 common variants in 3 genes or intergenic regions met p<10−5 in stage II. A secondary analysis using CK-MB as a quantitative trait (minimum p=1.26×10−3 for rs609418), or a dichotomised phenotype based on extreme CK-MB values (minimum p=7.72×10−6 for rs4834703) supported these findings. Pathway analysis revealed that genes harbouring top-scoring variants cluster in pathways of biological relevance to extracellular matrix remodelling, endoplasmic reticulum-to-Golgi transport and inflammation.

Conclusions Using a two-stage GWAS and pathway analysis, we identified and prioritised several potential susceptibility loci for perioperative MI.

  • SURGERY
  • GENETICS

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/

Statistics from Altmetric.com

Strengths and limitations of this study

  • This is the first genome-wide association study of perioperative myocardial infarction, using prospective cohorts of cardiac surgical patients, standard definitions of the primary phenotype and full adjustment for non-genetic risk factors.

  • We conducted comprehensive and complementary single marker and pathway-based genome-wide association analyses.

  • The study is powered to detect relatively large effect sizes.

  • Rare genetic variant effects not analysed.

  • Predominantly Caucasian cohort, thus findings cannot be generalised to other populations.

Introduction

Despite advances in surgical techniques and pharmacological therapy, the incidence of myocardial infarction (MI) after coronary artery bypass grafting (CABG) remains as high as 19%, and is associated with increased mortality and long-term morbidity.1 Strategies to identify subpopulations of patients at risk for developing large myocardial infarcts on the one hand, or those displaying an intrinsic cardioprotective state on the other hand, remain high-priority knowledge gaps and could inform selection of more specific protective agents.2

The evidence for heritability of MI is striking, supported both by family studies and, recently, by a number of well powered and replicated genome-wide association studies (GWAS), which primarily implicate common genetic variants at the 9p21 locus in multiple racial groups.3–6 However, although family-based methods are not practical for studying perioperative MI (PMI), its genetic basis is strongly suggested by several observations, including wide variability in incidence and severity that is poorly explained by clinical and procedural risk factors, different racial susceptibility profiles and results from preclinical animal models. Indeed, extensive genetic variability has been found in biological pathways implicated in the pathophysiology of postoperative MI, such as the complex acute inflammatory response to cardiac surgery. Mounting evidence for heritability of a proinflammatory state suggests that individual genetic history may also significantly modulate the magnitude of postoperative inflammatory response after cardiac surgery.7 Yet only a few studies have identified allelic associations with altered susceptibility to myocardial ischaemia-reperfusion injury in cardiac surgical populations, all based on a candidate gene association approach.8–12 Thus, the overall influence of common genetic variation on the incidence of PMI remains poorly understood.

Recently, integrated testing of genes involved in the same biological pathway has emerged as an alternative strategy for evaluating the combined effects of multiple genetic variants with small effect size on a disease phenotype.13–15 Given the polygenic nature of disease susceptibility, this approach is increasingly being used to identify groups of gene variants with shared cellular function that are enriched for disease, while also improving the statistical power of GWAS. In this study, we adopted this strategy by first employing genome-wide association methodology to identify common genetic variants associated with PMI after CABG, followed by pathway-based analyses to uncover biological mechanisms of relevance to PMI.

Methods

The study design and reporting of the results follow the recommendations by ‘Strengthening the Reporting of Genetic Association Studies’ (STREGA).16 We performed a joint two-stage17 GWAS combined with a pathway analysis approach.

Patient populations

The stage I cohort (discovery cohort) comprised 1493 ethnically diverse subjects who underwent an isolated CABG with cardiopulmonary bypass for the first time and were enrolled between 2002 and 2003 in the Institutional Review Board (IRB) approved Genetics of Myocardial Outcomes and Graft Failure (GeneMAGIC) ancillary substudy of the multicentre Project of Ex-vivo Vein Graft Engineering via Transfection (PREVENT-IV).18 Of those, 1433 patients met eligibility requirements after applying quality control criteria and excluding patients with missing genotypes or phenotypic information.

In stage II, we tested the top genetic variants in an expanded data set in which an additional 622 patients of self-reported European ancestry were added to the discovery data set, leading to a total of 2055 patients. The additional patients underwent CABG with cardiopulmonary bypass between 1997 and 2006 as part of the Perioperative Genetics and Safety Outcomes Study, an IRB-approved longitudinal study at the Duke University Medical Center.11 ,12

Definition of PMI

PMI was defined according to the universal definition of MI19 as an elevation in the plasma level of creatine kinase MB isoenzyme (CK-MB) that was >10 times the upper limit of normal, as measured by a core laboratory within 24 h after surgery, and that was not attributable to an intervening clinical event or preoperative MI (adjudicated by the PREVENT-IV Clinical Events Committee).18

Genotyping and quality controls

Genomic DNA was isolated from whole blood or saliva using standard procedures. Genotyping in both cohorts was performed on the Illumina Human610-Quad BeadChip at the Duke Genomic Analysis Facility. Sample and genotype quality control of data flow included assessment of call rates, gender check, cryptic relatedness, SNP missingness and the Hardy-Weinberg equilibrium, as previously described.20 We used principal components derived from the EIGENSTRAT21 method to control for population stratification (see online supplementary methods).

Statistical analysis

Univariate regression analysis was performed to test differences in demographic, clinical and procedural characteristics between patients with and without postoperative MI; statistically significant covariates were subsequently used to adjust genetic association tests. Genome-wide association analyses performed in the stage I cohort used multivariable logistic regression models implemented in PLINK 1.07, assuming an additive genetic model, and including significant clinical covariates and the top 10 principal components to adjust for population stratification (see online supplementary table S1). Statistical significance for stage I analyses was a priori arbitrarily defined as a two-tailed p<0.001, to balance between the overly conservative Bonferroni correction and type II error, given that we had an a priori defined replication data set to obviate type I error. In stage II analyses, the same clinical covariate and principal component adjustments (re-estimated for the expanded data set) were applied, with statistical significance defined as a Bonferroni adjusted p value of 0.05/number of SNPs identified in stage I, to control the overall type I error rate. Haplotype association analysis was performed in the stage II cohort for genes tagged by the significant SNPs, adjusting for the same clinical covariates and principal components. Secondary analyses were performed for the top SNPs identified in the stage II cohort, using the continuous phenotype of CK-MB values and a dichotomised extreme CK-MB phenotype, which included participants within the first and fourth quartiles of the CK-MB distribution. For CK-MB, a linear regression model was applied with adjustment for the same set of covariates as in the primary analysis. For the extreme CK-MB subset, we reassigned trait status as ‘0’ for participants within the first CK-MB quartile, and ‘1’ for participants within the fourth CK-MB quartile from the stage II cohort. Logistic regression with the same set of covariates was performed. Finally, we employed pathway analysis to prioritise association results and provide biological interpretation, using functional ontology enrichment analysis tools implemented in MetaCore (GeneGO, St. Joseph, MI; see online supplementary methods).

Results

PMI was observed in 112 of 1433 patients (7.8%) and in 225 of 2055 patients (10.9%) in stage I and II cohorts, respectively (see online supplementary table S2). Aortic cross-clamp time, number of coronary artery grafts, and procedures in addition to CABG were significantly associated with PMI in stage I analyses. Stage II analyses showed that extracardiac arteriopathy and year of surgery also played significant roles. All the clinical variables identified above were included in multivariable logistic regression models to adjust the genetic association tests for potential confounders (see online supplementary table S2).

After applying quality control criteria (see online supplementary methods: A.1. Quality control of data flow), 534 390 markers were analysed for association with PMI in stage I. While none of the SNPs reached genome-wide significance (figure 1), 521 SNPs met the a priori defined discovery threshold of p<0.001 (minimum p=2.76×10−6, rs2044061 on chromosome 8) and were subsequently analysed in stage II (see online supplementary table S10). A quantile-quantile (Q-Q) plot of observed versus expected p values showed that the population substructure was well adjusted for (see online supplementary figure S2).

Figure 1

Manhattan plot of genome-wide association with perioperative myocardial infarction in stage I analysis. The x axis represents the genome in physical order (coloured by chromosome); the y axis showing –log10(p) for all single nucleotide polymorphisms (SNPs). None of the SNPs reached genome-wide significance (p<9.09×10−8), but 521 SNPs met the prespecified discovery threshold p<0.001 for inclusion in stage II analyses (minimum p=2.76×10−6, rs2044061 on chromosome 8).

In stage II analyses, eight of the 521 SNPs met the Bonferroni correction threshold (p<9.6×10−5 for 521 SNPs). The top 2 SNPs (rs10454444 and rs10913237), located in the pregnancy-associated plasma protein A2 gene (PAPPA2), were in high linkage disequilibrium, with p values for stage II of 2.43×10−6 (OR 0.46; 95% CI 0.33 to 0.63) and 2.47×10−6 (OR 0.46; 95% CI 0.33 to 0.63), respectively. In addition, one intronic SNP in histone deacetylase-4 (HDAC4), two in the SEC24 family, member D (SEC24D, a member of the cytoplasmic coat protein complex II, COPII), and two located in intergenic regions were associated with PMI (table 1).

Table 1

Top eight SNP associated with postoperative myocardial infarction

When CK-MB was tested as a quantitative trait, all eight SNPs remained significantly associated with plasma levels of CK-MB. However, rs4834703 in SEC24D (p=7.72×10−6) and rs10200850 in HDAC4 (p=8.26×10−5) showed the strongest association, followed by the 2 SNPs in PAPPA2 (p=0.0001; table 2). When the dichotomised extreme CK-MB phenotype was studied in the stage II cohort, only the 2 SNPs located in SEC24D (rs4834703 and rs6822035, p=0.002) and rs609418 in the intergenic region between RFXAP and SMAD9 (p=0.001) remained nominally significant (table 2). Overall, the rs4834703 SNP in SEC24D consistently showed strong association signals across postoperative MI, quantitative CK-MB and the extreme CK-MB traits.

Table 2

Association of the top eight SNP with CK-MB as a quantitative trait and as an extreme phenotype in the joint analysis data set

The haplotype structures of genic regions surrounding the significant SNPs in PAPPA2, HDAC4 and SEC24D are shown in online supplementary figure S3, and the results of haplotype analysis are shown in table 3 and online supplemental table S3. Haplotype analysis performed for the linkage disequilibrium blocks containing the significant markers showed that the A-A haplotype (rs6822035, rs10518325) in SEC24D had the most significant association with postoperative MI (p=5.54×10−7, OR 1.87; 95% CI 1.46 to 2.39, table 3).

Table 3

Estimated haplotype frequencies in the SEC24D gene and results of association tests with incidence of PMI in the stage II analysis cohort (n=2055)

In multivariate risk factor-adjusted stage I and II analyses, we found no evidence for association between genetic variation at the 9p21 locus and incident postoperative MI, including the rs10116277 SNP previously associated with MI and mortality in non-surgical20 and cardiac surgical8 ,9 cohorts (see online supplementary figure S4).

The top 10 enriched pathway maps for the stage I analysis SNPs are shown in online supplementary figure S5. The most significant were ‘Cell adhesion: extracellular matrix remodelling’ and ‘Cytoskeleton remodelling: TGF, WNT, and cytoskeletal remodelling’ (enrichment p value=1.9×10−8 and 7.3×10−6, respectively). The most significant canonical pathway maps in stage I and II comparative analyses were ‘Immune response: NFAT signalling and leucocyte interaction’ and ‘Cell adhesion: extracellular matrix remodelling’ (p=7.2×10−4 and 1.1×10−3, respectively; see online supplementary tables S4 and 5).

Discussion

We present the first report of a two-stage GWAS involving 225 PMI cases and 1830 controls. After accounting for clinical and procedural covariates, we identified eight significant SNPs mapped to three genes (PAPPA2, HDAC4, SEC24D) and two intergenic regions. The most significant association with PMI was exhibited by rs10454444 in PAPPA2 (p=2.43×10−6) in single-marker analyses, by SEC24D A-A (rs6822035, rs10518325, p=5.54×10−7) in haplotype analyses, and by SEC24D SNPs in secondary analyses when CK-MB was evaluated as a quantitative trait (rs4834703, p=2.43×10−6) or as an extreme CK-MB phenotype (rs4834703 and rs6822035, p=0.002).

Among these eight SNPs, rs4834703 in SEC24D showed the most consistently strong association with all three phenotypes evaluated in this study. The transport protein SEC24D is an integral component of cytoplasmic COPII transport machinery, a key player in vesicle trafficking of secretory proteins from the endoplasmic reticulum (ER) to the Golgi apparatus for delivery to downstream compartments. COPII is responsible for cargo sorting and vesicle morphogenesis, with roles in modulating ER exit, cell surface transport, lipid secretion and cholesterol biosynthesis, and function of G protein-coupled receptors. Conditions of ischaemia, oxidative injury or acute phase-response result in ER stress through accumulation of misfolded proteins, which leads to activation of the unfolded protein response (UPR) signalling pathway. If the protective mechanisms activated by the UPR are insufficient, cells die by apoptosis and autophagy. However, altered expression or function of SEC24 proteins could also explain ER trapping of misfolded proteins under conditions of ER stress,22 ,23 with SEC24D being the only isoform implicated in extracellular matrix secretion.24 Animal models of defects in SEC24D involve characteristic skeletal malformations, whereas its complete disruption results in early-embryonic lethality.25 However, there are no reports of human diseases associated with genetic variation in SEC24D.

PAPPA2 is a metalloproteinase that regulates local insulin-like growth factor (IGF) bioavailability by specifically cleaving IGF-binding protein 5 (IGFBP-5). In experimental models of myocardial ischaemia, IGFBP-5 inhibits both myocardial IGF-1, with implications for inflammation-linked angiogenesis and repair processes,26 and IGF-2, limiting its cardioprotective effects following ischaemia reperfusion.27 ,28 The related protein PAPPA1 cleaves IGFBP-4 and is activated and released from vulnerable atherosclerotic plaques. PAPPA1 has been extensively studied as a cardiovascular risk biomarker for the diagnosis and prognosis of acute coronary syndrome;29 ,30 however, the role of PAPPA2 in cardiovascular biology has not been previously reported. Consistent with our single-marker analysis findings, the GWAS pathway analysis identified significant contributions of variants in other IGF system component genes, namely IGF-2 and the IGF-1 receptor (see online supplementary figure S6 and table S5).

HDAC4, a member of the HDAC family, mediates changes in the chromatin structure by removing acetyl groups from the core histones, resulting in transcriptional repression. HDAC4 is highly expressed in the myocardium, where it plays an important role in the regulation of gene expression and is involved in myocardial cell cycle progression, differentiation and apoptosis.31 Experimental HDAC inhibition is associated with a profound reduction in ischaemia-induced myocardial cell death,31 by triggering preconditioning effects and promoting myocardial repair. Genetic variants in HDACs are important determinants of susceptibility to cardiovascular diseases,32 suggesting functional roles for HDAC4 in modulating perioperative myonecrosis, and their potential use in predicting individual patient responsiveness to HDAC inhibition.

One of the intergenic SNPs associated with PMI (rs609418) is located near the mother against the decapentaplegic homologue 9 (SMAD9) gene, part of the transforming growth factor β (TGF-β) signalling pathway, which is markedly activated in the infarcted myocardium. Members of the TGF-β superfamily transduce their signal from the membrane to the nucleus via a distinct combination of transmembrane receptors and downstream effectors—the SMAD proteins.33 TGF-β plays important and complex roles in regulating postinfarction inflammatory responses. In animal models, signalling through the SMAD transcription factors is associated with resolution of inflammation, repression of cytokine and chemokine gene synthesis, and protection against myocardial ischaemia-reperfusion injury.34 Consistent with this single-locus gene association result, the bone morphogenetic protein pathway, which also transduces its signals via a SMAD9-dependent cascade, was identified as one of the top-scoring pathways in our GWAS pathway analysis (p=5.2×10−3, see online supplementary table S7). Although no direct functional roles are currently attributed to this intergenic region, a query of the Regulome and Haploreg databases reports that rs609418 is located within active regulatory elements (GATA2 transcription factor binding site by ENCODE ChIP-seq and altering regulatory motifs in Gfi1 and Mef2 by the Position-Weight Matrix, see online supplementary figure S7).

Surprisingly, we have been unable to replicate previously reported associations between common genetic variants at the 9p21 locus and risk for PMI or mortality after CABG.8 ,9 Many different reasons may account for this, including inadequate sample size, variation in study design, differences in allele frequencies or variability in the definition of PMI phenotype.35 Indeed, our observation is consistent with other studies showing that, although genetic variants at the 9p21 locus are associated with incident coronary artery disease (CAD), they may not be associated with the actual risk of MI.36 This discrepancy between convincing associations of 9p21 with a greater burden of CAD but not with MI in the presence of underlying CAD has been further confirmed by nested case–control studies37 as well as meta-analyses.5 ,38 ,39 Taken together, in subjects with established CAD, such as those included in our study, any lack of association of the 9p21 locus with subsequent MI could have resulted from the presence of CAD in carriers and non-carriers of the risk variants at the 9p21 locus, which seems to primarily mediate an atherosclerotic phenotype.

The strengths of our study are (1) a relatively large population of cardiac surgery patients, (2) a prospective cohort design and (3) a combination of complementary single-marker and pathway-based genome-wide association analyses. This approach allowed us to identify genetic variants that carry only a small disease risk individually, but that jointly can contribute relatively large effects on PMI susceptibility. Furthermore, an application of results from pathway-based analysis may add structure to interpreting genomic data and allow exploration of cellular processes that functionally underpin the observed associations. Finally, by using this approach, replication of association findings at the gene and pathway level is much easier compared to replication at the individual SNP level.13 Of note, most of the genes identified through pathway enrichment analysis in this study encode targets for therapeutic drugs that have already been developed. Thus, by improving risk assessment and identifying allele-specific therapeutic responses, further investigation of loci and pathways prioritised in this study could yield actionable results for enhancing perioperative cardioprotection.

Several limitations are worth mentioning. Power calculations (see online supplementary results B.6.) show that, based on the current sample size and incidence of PMI, our study can detect a genotypic relative risk approximately 2 with 80% power (assuming a variant with a 10% minor allele frequency and a realistic linkage disequilibrium between the tested marker and the causal locus D′=0.8). Thus, although our study is the largest genetic association study of PMI conducted to date, it is powered to detect only common variants with relatively large effect sizes. Most published genetic association studies of PMI have reported larger effect sizes compared with ambulatory populations (OR range 1.79–3.97).8 ,10 ,11 The possibility of rare genetic variants that drive a pronounced clinical phenotype was not explored in this study, because only variants with minor allele frequencies >0.05 were assessed. Also, functional studies to further elucidate the potential biological effects of the SNPs identified were not feasible due to lack of plasma or tissue availability in the GeneMAGIC study. Finally, patients enrolled in our cohorts were predominantly Caucasian, and therefore our findings cannot be generalised to other ethnic groups.

In conclusion, we report the first GWAS in a cohort of patients at risk for MI after CABG surgery. On the basis of our integrated approach utilising primary (PMI) and secondary phenotypes (quantitative CK-MB, extreme CK-MB), single-marker analysis and pathway analysis, we identified several polymorphisms in the insulin growth factor system implicated in the regulation of extracellular matrix remodelling, as well as the ER-to-Golgi secretory pathway potentially involved in adaptive responses to ER stress. As other GWAS of PMI cohorts publish their results, we intend to collaborate on conducting a meta-analysis for this particular phenotype. While our GWAS results are intriguing, follow-up studies are needed to translate these initial findings into biological insights that could lead to predictive and therapeutic advances in perioperative care. For instance, in the regions of confirmed associations, causal variants will only occasionally be among those directly genotyped. Moreover, GWAS detect almost exclusively the effects of common SNPs, offering limited power to capture any rare and structural variants, such as insertions, deletions, inversions and translocations. Detailed sequencing may be necessary to further characterise the genetic variations in the PMI-associated regions to enhance the identification of causal variants. Furthermore, as in most cases, the functions of identified genes and their variants, as well as the mechanisms by which they may contribute to PMI pathophysiology, are largely unknown. Currently, very few existing groups can bridge from genetics to molecular biology or cell physiology and to disease or novel therapy; we propose to examine the mechanisms by which variation in the observed genes is involved in myocardial ischaemia-reperfusion injury by examining their functional effects in the human myocardium from patients undergoing cardiac surgery, as well in our previously described preclinical rodent and swine models of cardioplegic arrest and cardiopulmonary bypass.40 The use of such animal models would allow pharmacological studies targeting the identified pathways to explore the mechanism by which novel cardioprotective drugs would attenuate myocardial ischaemia-reperfusion injury, with the ultimate goal of developing personalised cardioprotective strategies.

References

View Abstract
  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors MDK contributed to the analysis and interpretation of data, as well as the drafting of the manuscript. Y-JL contributed to the analysis and interpretation of data and a critical revision of the manuscript. Y-JL and Y-WL contributed to the quality control and analysis of the data. JA, MFN, PKS, DJ and JPM contributed to the acquisition of data and a critical review of the manuscript. MVP contributed to the conception and design, analysis and interpretation of data, critical manuscript review and approved the final version of the submitted manuscript.

  • Funding This work was supported in part by the National Institutes of Health grants R01-HL075273 and R01-HL092071 (to MVP) and by the American Heart Association grants 0256342U and 9951185U (to JPM), and 0120492U (to MVP).

  • Competing interests None.

  • Patient consent Obtained.

  • Ethics approval The Duke University School of Medicine Institutional Review Board.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.