Next Article in Journal
RNA Splicing in Cancer and Targeted Therapies
Next Article in Special Issue
Pharmacogenomics of Dementia: Personalizing the Treatment of Cognitive and Neuropsychiatric Symptoms
Previous Article in Journal
Three Blind Moles: Molecular Evolutionary Insights on the Tempo and Mode of Convergent Eye Degeneration in Notoryctes typhlops (Southern Marsupial Mole) and Two Chrysochlorids (Golden Moles)
Previous Article in Special Issue
An Exploratory Application of Multilayer Networks and Pathway Analysis in Pharmacogenomics
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 11
10.3390/genes14112019
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction

by
Jasmine A. Luzum
1,2,*,
Alessandra M. Campos-Staffico
1,
Jia Li
3,
Ruicong She
3,
Hongsheng Gui
2,
Edward L. Peterson
3,
Bin Liu
3,
Hani N. Sabbah
4,
Mark P. Donahue
5,
William E. Kraus
5,
L. Keoki Williams
2 and
David E. Lanfear
2,4
1
Department of Clinical Pharmacy, University of Michigan College of Pharmacy, Ann Arbor, MI 48109, USA
2
Center for Individualized and Genomic Medicine Research (CIGMA), Henry Ford Health System, Detroit, MI 48202, USA
3
Department of Public Health Sciences, Henry Ford Health System, Detroit, MI 48202, USA
4
Heart and Vascular Institute, Henry Ford Health System, Detroit, MI 48202, USA
5
School of Medicine, Duke University, Durham, NC 27710, USA
*
Author to whom correspondence should be addressed.
Genes 2023, 14(11), 2019; https://doi.org/10.3390/genes14112019
Submission received: 6 July 2023 / Revised: 20 October 2023 / Accepted: 22 October 2023 / Published: 28 October 2023
(This article belongs to the Special Issue Pharmacogenomics: Challenges and Future)
Browse Figures
Review Reports Versions Notes

Abstract

:
In patients with heart failure with reduced ejection fraction (HFrEF), individual responses to beta-blockers vary. Candidate gene pharmacogenetic studies yielded significant but inconsistent results, and they may have missed important associations. Our objective was to use an unbiased genome-wide association study (GWAS) to identify loci influencing beta-blocker survival benefit in HFrEF patients. Genetic variant × beta-blocker exposure interactions were tested in Cox proportional hazards models for all-cause mortality stratified by self-identified race. The models were adjusted for clinical risk factors and propensity scores. A prospective HFrEF registry (469 black and 459 white patients) was used for discovery, and linkage disequilibrium (LD) clumped variants with a beta-blocker interaction of p < 5 × 10−5, were tested for Bonferroni-corrected validation in a multicenter HFrEF clinical trial (288 black and 579 white patients). A total of 229 and 18 variants in black and white HFrEF patients, respectively, had interactions with beta-blocker exposure at p < 5 × 10−5 upon discovery. After LD-clumping, 100 variants and 4 variants in the black and white patients, respectively, remained for validation but none reached statistical significance. In conclusion, genetic variants of potential interest were identified in a discovery-based GWAS of beta-blocker survival benefit in HFrEF patients, but none were validated in an independent dataset. Larger cohorts or alternative approaches, such as polygenic scores, are needed.

1. Introduction

Heart failure (HF) is a global pandemic that affects at least 26 million individuals [1]. Despite improvements in HF therapies over the past few decades, the one-year mortality rate is still approximately 30% globally [1]. Beta-blockers are one of the few pharmacologic therapies that can significantly reduce morbidity and mortality in patients with heart failure with reduced ejection fraction (HFrEF) [2]. Thus, beta-blockers are a cornerstone for guideline-recommended HFrEF therapy [2]. The landmark beta-blocker clinical trials for HFrEF demonstrated significantly improved outcomes, on average, in a large patient sample [3,4,5]; however, a substantial portion of HFrEF patients did not receive any benefit. Only ~22% of HFrEF patients demonstrated a marked and sustained improvement in their left ventricular ejection fraction (LVEF) with beta-blocker therapy [6,7]. In a randomized, double-blind clinical trial of HFrEF patients, the mean ± sd of the change in LVEF, after approximately 12 months of treatment, was +10.9% ± 11.0% and +7.2% ± 7.7% for carvedilol and metoprolol, respectively [8]. With the global prevalence of HF continuing to rise, the need to identify HFrEF patients that are most likely to benefit from this potentially life-saving therapy is of paramount importance.
Unfortunately, clinical characteristics do not predict which HFrEF patients will benefit from beta-blocker therapy. The clinical subgroup analyses of beta-blocker clinical trials showed similar benefit regardless of a patient’s age, sex, comorbidities, or HFrEF etiology [3,4,5]. More recent studies have shown that the benefits of beta-blockers are similar, regardless of a patient’s race [9]. Therefore, research efforts have focused on identifying novel factors that may predict beta-blocker benefits in HFrEF patients, such as genetics [10,11]. Pharmacogenetic studies of beta-blockers in HFrEF have primarily used the candidate gene approach (i.e., only a few genes are selected for analysis in the genome, based on a priori hypotheses that genes will play a role). Beta-blockers inhibit the neurohormonal drivers of HF pathophysiology via the sympathetic adrenergic and renin–angiotensin–aldosterone systems [12,13,14,15]. Accordingly, pharmacogenetic studies of beta-blockers in HFrEF patients have primarily focused on candidate genetic variants from these systems [10,11]. However, the results from candidate gene association studies have been inconsistent [10,11], which negates clinical utility, and suggests that additional and/or alternative genetic variants may be involved. The major limitation of candidate gene studies is that they can miss thousands of other genes in the genome, which may play a role. However, a genome-wide association study (GWAS) of beta-blocker survival benefit in patients with HFrEF has not previously been conducted. Therefore, the objective of this study was to use a GWAS to discover and validate novel genetic variant(s) that significantly interact with beta-blocker survival benefit in patients with HFrEF.

2. Materials and Methods

2.1. Discovery Dataset

The discovery dataset came from a prospective HF pharmacogenomic registry (HFGR) designed to discover novel ways to predict HF prognoses and responses to HF therapies [9]. The registry began enrolling patients in October 2007, within the Henry Ford Health System in Detroit, MI, USA, and enrolment was completed in March 2015. The inclusion criteria were as follows: age ≥ 18 years; a diagnosis of HF using the Framingham Heart Study criteria [16]; and at least one documented LVEF < 50% measured via echocardiography, nuclear stress tests, or radionuclide blood pool imaging. The LVEF cutoff was chosen to reflect patients with systolic HF because the study was designed and started prior to more recent reclassifications, which suggest that HFrEF should be defined as a LVEF ≤ 40% [2]. All patients were covered by the affiliated health insurance of the Henry Ford Health System, the Health Alliance Plan, which allowed access to claims data. Patients were excluded if they required dialysis or supplementary oxygen. Detailed clinical information (e.g., demographics, physical examination results, past medical history, laboratory values, functional status, and medication use) and blood samples were collected at the time of enrolment. Patient deaths were identified using the Social Security Administration Death Master File, the Michigan State Division of Vital Records, and administrative data from the Henry Ford Health System through 30 April 2017. The study received approval from the Institutional Review Board of the Henry Ford Health System. All patients provided written informed consent prior to study participation.

2.2. Validation Dataset

Data from the “Heart Failure: A Controlled Trial Investigating Outcomes of Exercise Training” (HF-ACTION) trial [17] were used for validation. Briefly, HF-ACTION was a multicenter, randomized, controlled trial of aerobic exercise training in medically stable outpatients with HFrEF. The trial included patients with LVEF ≤ 35% and New York Heart Association (NYHA) class II to IV symptoms despite optimal medical therapy for at least 6 weeks. Patients were randomized from April 2003 to February 2007 within the United States, Canada, and France. Exclusion criteria included major comorbidities or limitations that could interfere with exercise training; recent (≤6 weeks) major cardiovascular events or planned (≤6 months) procedures; performance of regular exercise training; or use of devices that limited the ability to achieve target heart rates. The intervention consisted of 36 supervised exercise sessions followed by home-based training versus usual care alone. The primary outcome was a composite of all-cause mortality or hospitalization, and the pre-specified secondary outcomes included all-cause mortality, cardiovascular mortality or cardiovascular hospitalization, and cardiovascular mortality or heart failure hospitalization. The protocol was reviewed and approved by the appropriate institutional review board or ethics committee for each participating center and by the coordinating center’s institutional review board. The exercise intervention did not significantly affect the primary or secondary outcomes in the protocol-specified primary analysis.

2.3. Beta-Blocker Exposure

In the HFGR, pharmacy claims data were used to calculate beta-blocker exposure as described previously [18]. Briefly, beta-blocker doses were standardized into daily dose equivalents using doses targeted in clinical HFrEF trials or using maximum recommended daily doses if HFrEF trial data were unavailable (e.g., atenolol). The following daily doses were used for each agent: atenolol, 100 mg; bisoprolol, 10 mg; carvedilol, 50 mg; labetalol, 600 mg; and metoprolol, 200 mg (all formulations). Beta-blocker exposure was calculated over a 6-month period by multiplying the standardized daily dose equivalents by the total quantity of medication dispensed during the period, and then dividing by the total number of days during the period. By utilizing pharmacy claims data, our beta-blocker exposure values also partially account for patient adherence over an extended 6-month period. We have previously shown that this method more strongly correlates with relevant HFrEF outcomes (e.g., heart rate, hospitalization, mortality) than single time point calculations (e.g., discharge medication status) for beta-blocker exposure [19]. In HF-ACTION, beta-blocker doses were standardized into daily dose equivalents the same as in the HFGR, but only the dose at baseline (and not time-varying) was used in the analysis.

2.4. Genotyping

The same genotyping methods were used for both the HFGR and HF-ACTION datasets. Blood samples were processed immediately and stored at −70 °C. Genotyping was performed by the University of Michigan genotyping core lab (Ann Arbor, MI, USA) using the Axiom Biobank Genotyping Array (Affymetrix®; Santa Clara, CA, USA) [9]. This array includes the following ~600 K genetic variants: 1) ~300 K genome-wide variants with minor allele frequencies >1%; (2) ~250 K low frequency (<1%) coding variants from global exome sequencing projects; and 3) an additional ~50 K variants to improve African ancestry coverage (Yoruba in Ibadan, Nigeria [YRI] booster). Imputation was performed using the Michigan Imputation Server with Minimac3 [20] and the 1000 Genomes Project as the reference panel. The following filters were used for genotyping quality control: call rate < 95%, minor allele frequency (MAF) < 0.01, Hardy–Weinberg Equilibrium (HWE) p-value < 1.0−8, imputation score r2 > 0.8, duplicate samples identified by identity-by-state distances, and monomorphic variants.

2.5. Statistical Analysis

Baseline characteristics were described by mean ± standard deviation for continuous variables, and counts and percentages for categorical variables. Baseline characteristics were compared by the Student’s t-test, chi-square test, or Fisher’s exact test as appropriate. Proportion of YRI ancestry was estimated using ADMIXTURE software v 1.3.0. The GWAS was performed in the discovery dataset using Cox proportional hazards regression models with the primary outcome of all-cause mortality stratified by self-identified black and white race. The models consisted of each individual genetic variant, beta-blocker exposure, and the multiplicative interaction term between each genetic variant and beta-blocker exposure (genetic variant * beta-blocker exposure). The genetic variants were tested using additive genetic models. Beta-blocker exposure was included as a time-updating variable, and the models were adjusted for the following covariates: the Meta-Analysis Global Group in Chronic Heart Failure (MAGGIC) risk score (excluding beta-blocker as an input variable) [21], beta-blocker propensity score, and the first two principal components. Treatment with angiotensin inhibitors, including angiotensin-converting enzyme inhibitors and angiotensin receptor blockers, were accounted for in the MAGGIC and propensity scores. Principal components (PC) analysis was performed using EIGENSOFT [22,23]. Propensity scores for beta-blocker exposure were calculated as previously described [9,24]. Genomic control was further used to correct for test statistic inflation due to population stratification [25]. The GWAS was performed using the survival package in R, and the other statistical analyses were performed using SAS v9.4. Variants of interest were defined as having a covariate-adjusted interaction p-value < 5 × 10−5 in the HFGR. Those variants were clumped by linkage disequilibrium ([LD] r2 > 0.8) using PLINK and then tested for validation in similar Cox proportional hazards regression models in the HF-ACTION dataset. Variants were considered validated if the Bonferroni-corrected p-value was statistically significant in HF-ACTION (i.e., p = 0.05/# of LD-clumped variants).

3. Results

3.1. Baseline Characteristics

Baseline characteristics are displayed in Table 1. In both the discovery and validation datasets, there were fewer white than black women. In addition, white patients were older, had a higher frequency of ischemic etiology and atrial fibrillation, and had higher MAGGIC risk scores [26] compared with black patients.
In the discovery dataset only, white patients had greater LVEF, while a higher frequency of diabetes, higher systolic blood pressure and heart rate, and higher serum creatinine were observed in black patients. In the validation dataset only, white patients had a higher frequency of COPD and higher levels of NT pro-BNP, while black patients had a higher frequency of stroke, body mass index, and beta-blocker exposure. There were no significant differences in duration of follow-up and frequency of deaths between white and black patients in both datasets. The mean proportion of YRI ancestry in the self-reported black and white patients in HFGR was 83.5% and 1.0%, respectively. The mean proportion of YRI ancestry in self-reported black and white patients in HF-ACTION was 88.9% and 1.2%, respectively.

3.2. Discovery GWAS

The quantile–quantile (QQ) plots are displayed in Figure 1. A total of 229 variants in the black patients (Figure 1A) had interactions with beta-blocker exposure at p < 5 × 10−5 and 18 variants in the white patients (Figure 1B). The complete lists of variants with p < 5 × 10−5 and their MAF are provided in the Supplementary Material. Manhattan plots are presented in Figure 2 (black patients) and Figure 3 (white patients). In the black patients, one variant (rs117032090 on chromosome 7) had p = 1.3 × 10−11.

3.3. Validation of GWAS Results

After LD clumping, 100 of the 229 variants with p < 5 × 10−5 in the black patients and 4 of the 18 variants with p < 5 × 10−5 in the white patients remained for validation testing in HF-ACTION. None of the variants met the Bonferroni-corrected level of statistical significance in HF-ACTION, in neither the black patients (p < 0.0005) nor white patients (p < 0.0125). The interaction p-value for the highly significant variant rs117032090 in the HFGR was 0.9999 in HF-ACTION. One variant was suggestive for an association in the black patients: rs16844448 had an interaction p-value = 8.1 × 10−6 in HFGR and p = 0.005 in HF-ACTION. The intronic rs16844448 variant is in LRP1B, the gene encoding low-density lipoprotein (LDL) receptor-related protein 1B. Table 2 shows the results from the Cox proportional hazards regression model in black patients for rs16844448.

4. Discussion

To the best of our knowledge, this is the first GWAS designed to discover and validate novel genetic variants related to beta-blocker survival benefit in black and white patients with HFrEF. According to our findings, 229 and 18 variants with covariate-adjusted interaction p-values < 5 × 10−5 were identified in black and white patients, respectively. After LD clumping, one-hundred variants in black patients and four variants in white patients with p-values < 5 × 10−5 were moved forward for validation testing. None of the variants met the Bonferroni-corrected level of statistical significance in black (p < 0.0005) nor white patients (p < 0.0125) in the validation testing. However, one variant (rs16844448) was found to be suggestive of association with beta-blocker survival benefit in black patients with HFrEF.
This suggestive variant (rs16844448) is in LRP1B, which encodes the low-density lipoprotein (LDL) receptor-related protein-1B (LRP1B). HaploReg v4.1 shows that this intronic variant (rs16844448) alters five regulatory motifs for CEBPG, Foxa_disc5, Nkx2_11, Nkx2-4, PLZF, and RXRA_known2 [27], but its functionality has not yet otherwise been described. However, the LRP1B gene is located on 2q22.1, a genomic region previously implicated in premature coronary atherosclerotic disease [28,29] and heart failure [30]. Also, this suggestive variant (rs16844448) is in strong linkage disequilibrium with two other intronic variants (rs67721025 and rs70988415) in black individuals [27].
Although multiple genetic variants were identified and selected for validation, none of them were ultimately identified as a statistically significant genetic variant related to beta-blocker survival benefit in patients with HFrEF. We speculate the following reasons that could potentially have contributed to the failed validation of these genetic variants. First, the clinical characteristics in the discovery (HFGR) and validation (HF-ACTION) datasets differed. Patients enrolled in HF-ACTION had more severe HFrEF than those in HFGR, as patients in HF-ACTION had lower LVEF and higher MAGGIC risk scores. Also, exposure to beta-blockers differed between the datasets. Patients in HF-ACTION were exposed to a higher daily dose of beta-blockers than those in HFGR, not to mention that patients in the HF-ACTION clinical trial most likely had greater adherence to drug treatment than patients in the observational HFGR study. Even though we aimed to identify robust and generalizable genetic variants for different populations of patients with HFrEF, we cannot disregard those differences in disease severity and drug exposure as factors that can intrinsically modify the clinical benefit of beta-blockers in HFrEF patients.
Despite using methods such as LD clumping and p-value thresholding, some of the top hits identified in the discovery GWAS were considered false positive after attempted validation. These false positive variants could have somehow correlated with mortality by chance due to the high number of genetic variants tested in the GWAS. On the other hand, it is possible that there are genetic variants with smaller effect sizes that significantly interact with the survival benefit of beta-blockers in patients with HFrEF, but they were not detected in the GWAS because our study may have been underpowered. There are a few factors that could have affected the power of our study. First, the relatively small sample sizes in both the discovery and validation datasets can be explained by the challenge of enrolling patients with HFrEF, which is considered the deadliest subtype of heart failure [31], and with whole-genome data available. Second, although beta-blockers are guideline-recommended drugs with undeniable benefits for HFrEF [2], prescription filling rates for these drugs in patients with HFrEF are still considered low in the U.S. [32]. These two reasons combined contribute greatly to the fact that studies of the survival benefit of beta-blockers in HFrEF, such as ours, have small sample sizes. Moreover, stratifying the patients by black and white race further reduced the sample size, and testing for drug–gene interaction rather than the main effect further contributed to a reduction in statistical power [33].
Although this small sample size limitation exists in our study, it is noteworthy that pharmacogenomic GWASs with even smaller samples were successful in identifying and validating novel genetic variants involved in drug response. This is possible because the effect sizes of pharmacogenetic variants can be extremely large compared to the genetic variants affecting common, complex diseases. In contrast to common, complex diseases, drug responses have had less time to be affected by selective evolutionary pressures, as most drugs have only been invented in the past 100 years. As an example, the study that identified the HLA-B*5701 genotype as the determinant of flucloxacillin-induced liver injury included only 333 individuals (of which 51 were cases) and reproduced its findings in a cohort with only 23 cases [34]. The odds ratio for the HLA-B*5701 variant in this small sample was 80.6 (p = 9.0 × 10−19). Considering that the polygenic score we developed was validated in four independent datasets [35,36], this corroborates the notion that the survival benefit of beta-blockers is not only polygenic but is also formed by several small-effect genetic variants.
The interest in genomic studies that may explain the different responses to beta-blockers is not recent. In 2015, the first pharmacogenomic GWAS study failed to identify significant genetic variants that could reduce the risk of cardiovascular outcomes in patients with hypertension treated with beta-blockers [37]. However, subsequent GWASs have identified genetic variants associated with uncontrolled [38] and systolic blood pressure response [39] in patients with hypertension treated with beta-blockers. In other GWAS studies, two novel genetic variants were identified and associated with changes in heart rate in response to beta-blockers [40] as well as a genetic variant that has been associated with lower odds for new-onset diabetes in beta-blocker users [41]. To our knowledge, this is the first GWAS aimed at identifying genetic variants that predict the survival benefit from beta-blockers in HFrEF patients. Our findings could differ from these other studies mainly because the underlying mechanisms driving beta-blocker benefit in heart failure are different than in hypertension, as they do not involve lowering blood pressure [42].
Some limitations must be taken into consideration when reading this study. Low statistical power may have limited the discovery of genetic variants by the GWAS, especially for those with smaller effects. It is worth mentioning that the large clinical trials that demonstrated survival benefit of beta-blockers in heart failure date back to 20 years ago—when DNA samples from participants were not routinely collected—hindering pharmacogenomic analyses. In addition, the discovery dataset used in the GWAS was derived from a single center, and the beta-blocker exposure assessed in the validation dataset was related only to the dose at baseline—and was not time-varying. Data distinguishing potential etiologies of HFrEF other than ischemic vs. non-ischemic were not available.

5. Conclusions

Many variants of interest were identified in this discovery GWAS of beta-blocker survival benefit in patients with HFrEF, but none of the variants were validated in an independent dataset.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes14112019/s1, the complete list of variants with p < 5 × 10−5 and their annotations.

Author Contributions

Conceptualization, D.E.L.; data curation, B.L.; formal analysis, J.L., R.S., H.G. and E.L.P.; funding acquisition, J.A.L., A.M.C.-S., H.N.S., L.K.W. and D.E.L.; methodology, J.A.L., J.L., R.S., H.G., E.L.P., B.L., H.N.S., M.P.D., W.E.K., L.K.W. and D.E.L.; supervision, J.L., L.K.W. and D.E.L.; validation, J.A.L., J.L., R.S., E.L.P., M.P.D., W.E.K. and D.E.L.; writing—original draft, J.A.L. and A.M.C.-S.; writing—review and editing, J.A.L., A.M.C.-S., J.L., R.S., H.G., E.L.P., B.L., H.N.S., M.P.D., W.E.K., L.K.W. and D.E.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was primarily supported by grant # R01HL103871 to Dr. Lanfear from the National Heart, Lung, and Blood Institute of the NIH, with other support for Dr. Lanfear—R01HL132154 and P50MD017351; Dr. Luzum—K08HL146990 and L30HL110279; Dr. Sabbah—P01HL074237 and R01HL132154; Dr. Williams—R01HL118267 and R01HL141845; and Dr. Campos-Staffico—F32HL162231 and the Michigan Institute for Clinical & Health Research (U074551). Dr. Williams is also supported by the National Institute of Allergy and Infectious Diseases (R01AI079139) and the National Institute of Diabetes and Digestive and Kidney Diseases (R01DK064695 and R01DK113003). Dr. Luzum was also supported by a Futures Grant from the American College of Clinical Pharmacy.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the Henry Ford Health System (#4562).

Informed Consent Statement

All participants provided written voluntary informed consent.

Data Availability Statement

The data from HFGR are publicly available via the database of Genotypes and Phenotypes (https://www-ncbi-nlm-nih-gov.proxy.lib.umich.edu/gap (accessed on 23 October 2023). The remaining data that support the findings of this study may be made available upon reasonable request from qualified researchers trained in human subject confidentiality protocols by contacting D.E.L. ([email protected]) at Henry Ford Hospital.

Conflicts of Interest

Dr. Lanfear is a consultant for Amgen, Janssen, Ortho Diagnostics, and Duke Clinical Research Institute (Novartis); has participated in clinical trials from Amgen, Bayer, and Janssen; and has submitted a provisional patent request for the β-blocker polygenic score. Dr. Luzum is a consultant for Ariel Precision Medicine. The other authors report no conflict.

References

  1. Savarese, G.; Lund, L.H. Global Public Health Burden of Heart Failure. Card. Fail. Rev. 2017, 3, 7–11. [Google Scholar] [CrossRef]
  2. Maddox, T.M.; Januzzi, J.L.; Allen, L.A.; Breathett, K.; Butler, J.; Davis, L.L.; Fonarow, G.C.; Ibrahim, N.E.; Lindenfeld, J.; Masoudi, F.A.; et al. 2021 Update to the 2017 ACC Expert Consensus Decision Pathway for Optimization of Heart Failure Treatment: Answers to 10 Pivotal Issues About Heart Failure With Reduced Ejection Fraction. J. Am. Coll. Cardiol. 2021, 77, 772–810. [Google Scholar] [CrossRef] [PubMed]
  3. MERIT-HF Investigators. Effect of Metoprolol CR/XL in Chronic Heart Failure: Metoprolol CR/XL Randomised Intervention Trial in Congestive Heart Failure (MERIT-HF). Lancet 1999, 353, 2001–2007. [Google Scholar] [CrossRef]
  4. CIBIS Investigators and Committees. A Randomized Trial of Beta-Blockade in Heart Failure. The Cardiac Insufficiency Bisoprolol Study (CIBIS). Circulation 1994, 90, 1765–1773. [Google Scholar] [CrossRef] [PubMed]
  5. Packer, M.; Bristow, M.R.; Cohn, J.N.; Colucci, W.S.; Fowler, M.B.; Gilbert, E.M.; Shusterman, N.H. The Effect of Carvedilol on Morbidity and Mortality in Patients with Chronic Heart Failure. U.S. Carvedilol Heart Failure Study Group. N. Engl. J. Med. 1996, 334, 1349–1355. [Google Scholar] [CrossRef] [PubMed]
  6. Metra, M.; Nodari, S.; Parrinello, G.; Giubbini, R.; Manca, C.; Dei Cas, L. Marked Improvement in Left Ventricular Ejection Fraction during Long-Term Beta-Blockade in Patients with Chronic Heart Failure: Clinical Correlates and Prognostic Significance. Am. Heart J. 2003, 145, 292–299. [Google Scholar] [CrossRef]
  7. Lupón, J.; Díez-López, C.; de Antonio, M.; Domingo, M.; Zamora, E.; Moliner, P.; González, B.; Santesmases, J.; Troya, M.I.; Bayés-Genís, A. Recovered Heart Failure with Reduced Ejection Fraction and Outcomes: A Prospective Study. Eur. J. Heart Fail. 2017, 19, 1615–1623. [Google Scholar] [CrossRef]
  8. Metra, M.; Giubbini, R.; Nodari, S.; Boldi, E.; Modena, M.G.; Dei Cas, L. Differential Effects of Beta-Blockers in Patients with Heart Failure: A Prospective, Randomized, Double-Blind Comparison of the Long-Term Effects of Metoprolol versus Carvedilol. Circulation 2000, 102, 546–551. [Google Scholar] [CrossRef]
  9. Luzum, J.A.; Peterson, E.; Li, J.; She, R.; Gui, H.; Liu, B.; Spertus, J.A.; Pinto, Y.M.; Williams, L.K.; Sabbah, H.N.; et al. Race and Beta-Blocker Survival Benefit in Patients With Heart Failure: An Investigation of Self-Reported Race and Proportion of African Genetic Ancestry. J. Am. Heart Assoc. 2018, 7, e007956. [Google Scholar] [CrossRef]
  10. Talameh, J.A.; Lanfear, D. Pharmacogenetics in Chronic Heart Failure: New Developments and Current Challenges. Curr. Heart Fail. Rep. 2012, 9, 23–32. [Google Scholar] [CrossRef]
  11. Talameh, J.A.; McLeod, H.L.; Adams, K.F.; Patterson, J.H. Genetic Tailoring of Pharmacotherapy in Heart Failure: Optimize the Old, While We Wait for Something New. J. Card. Fail. 2012, 18, 338–349. [Google Scholar] [CrossRef]
  12. Bristow, M.R. Beta-Adrenergic Receptor Blockade in Chronic Heart Failure. Circulation 2000, 101, 558–569. [Google Scholar] [CrossRef] [PubMed]
  13. Holmer, S.R.; Hengstenberg, C.; Mayer, B.; Engel, S.; Löwel, H.; Riegger, G.A.; Schunkert, H. Marked Suppression of Renin Levels by Beta-Receptor Blocker in Patients Treated with Standard Heart Failure Therapy: A Potential Mechanism of Benefit from Beta-Blockade. J. Intern. Med. 2001, 249, 167–172. [Google Scholar] [CrossRef]
  14. Packer, M. The Neurohormonal Hypothesis: A Theory to Explain the Mechanism of Disease Progression in Heart Failure. J. Am. Coll. Cardiol. 1992, 20, 248–254. [Google Scholar] [CrossRef]
  15. Sackner-Bernstein, J.D.; Mancini, D.M. Rationale for Treatment of Patients with Chronic Heart Failure with Adrenergic Blockade. JAMA 1995, 274, 1462–1467. [Google Scholar] [CrossRef]
  16. McKee, P.A.; Castelli, W.P.; McNamara, P.M.; Kannel, W.B. The Natural History of Congestive Heart Failure: The Framingham Study. N. Engl. J. Med. 1971, 285, 1441–1446. [Google Scholar] [CrossRef] [PubMed]
  17. O’Connor, C.M.; Whellan, D.J.; Lee, K.L.; Keteyian, S.J.; Cooper, L.S.; Ellis, S.J.; Leifer, E.S.; Kraus, W.E.; Kitzman, D.W.; Blumenthal, J.A.; et al. Efficacy and Safety of Exercise Training in Patients with Chronic Heart Failure: HF-ACTION Randomized Controlled Trial. JAMA 2009, 301, 1439–1450. [Google Scholar] [CrossRef] [PubMed]
  18. Lanfear, D.E.; Hrobowski, T.N.; Peterson, E.L.; Wells, K.E.; Swadia, T.V.; Spertus, J.A.; Williams, L.K. Association of β-Blocker Exposure with Outcomes in Heart Failure Differs between African American and White Patients. Circ. Heart Fail. 2012, 5, 202–208. [Google Scholar] [CrossRef]
  19. Lanfear, D.; Peterson, E.; Wells, K.; Williams, L.K. Discharge Medication Status Compares Poorly with Claims-Based Outpatient Medication Exposure Estimates. Circ. Cardiovasc. Qual. Outcomes 2011, 4, AP234. [Google Scholar] [CrossRef]
  20. Das, S.; Forer, L.; Schönherr, S.; Sidore, C.; Locke, A.E.; Kwong, A.; Vrieze, S.I.; Chew, E.Y.; Levy, S.; McGue, M.; et al. Next-Generation Genotype Imputation Service and Methods. Nat. Genet. 2016, 48, 1284–1287. [Google Scholar] [CrossRef]
  21. Pocock, S.J.; Ariti, C.A.; McMurray, J.J.V.; Maggioni, A.; Køber, L.; Squire, I.B.; Swedberg, K.; Dobson, J.; Poppe, K.K.; Whalley, G.A.; et al. Predicting Survival in Heart Failure: A Risk Score Based on 39 372 Patients from 30 Studies. Eur. Heart J. 2013, 34, 1404–1413. [Google Scholar] [CrossRef]
  22. Price, A.L.; Patterson, N.J.; Plenge, R.M.; Weinblatt, M.E.; Shadick, N.A.; Reich, D. Principal Components Analysis Corrects for Stratification in Genome-Wide Association Studies. Nat. Genet. 2006, 38, 904–909. [Google Scholar] [CrossRef] [PubMed]
  23. Patterson, N.; Price, A.L.; Reich, D. Population Structure and Eigenanalysis. PLoS Genet. 2006, 2, e190. [Google Scholar] [CrossRef] [PubMed]
  24. D’Agostino, R.B. Propensity Score Methods for Bias Reduction in the Comparison of a Treatment to a Non-Randomized Control Group. Stat. Med. 1998, 17, 2265–2281. [Google Scholar] [CrossRef]
  25. Devlin, B.; Roeder, K.; Wasserman, L. Genomic Control for Association Studies: A Semiparametric Test to Detect Excess-Haplotype Sharing. Biostatistics 2000, 1, 369–387. [Google Scholar] [CrossRef]
  26. Rich, J.D.; Burns, J.; Freed, B.H.; Maurer, M.S.; Burkhoff, D.; Shah, S.J. Meta-Analysis Global Group in Chronic (MAGGIC) Heart Failure Risk Score: Validation of a Simple Tool for the Prediction of Morbidity and Mortality in Heart Failure With Preserved Ejection Fraction. J. Am. Heart Assoc. 2018, 7, e009594. [Google Scholar] [CrossRef] [PubMed]
  27. HaploReg v4.1. Available online: https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php (accessed on 3 March 2022).
  28. Pajukanta, P.; Cargill, M.; Viitanen, L.; Nuotio, I.; Kareinen, A.; Perola, M.; Terwilliger, J.D.; Kempas, E.; Daly, M.; Lilja, H.; et al. Two Loci on Chromosomes 2 and X for Premature Coronary Heart Disease Identified in Early- and Late-Settlement Populations of Finland. Am. J. Hum. Genet. 2000, 67, 1481–1493. [Google Scholar] [CrossRef] [PubMed]
  29. Divers, J.; Palmer, N.D.; Langefeld, C.D.; Brown, W.M.; Lu, L.; Hicks, P.J.; Smith, S.C.; Xu, J.; Terry, J.G.; Register, T.C.; et al. Genome-Wide Association Study of Coronary Artery Calcified Atherosclerotic Plaque in African Americans with Type 2 Diabetes. BMC Genet. 2017, 18, 105. [Google Scholar] [CrossRef]
  30. Kulminski, A.M.; He, L.; Culminskaya, I.; Loika, Y.; Kernogitski, Y.; Arbeev, K.G.; Loiko, E.; Arbeeva, L.; Bagley, O.; Duan, M.; et al. Pleiotropic Associations of Allelic Variants in a 2q22 Region with Risks of Major Human Diseases and Mortality. PLoS Genet. 2016, 12, e1006314. [Google Scholar] [CrossRef]
  31. Lam, C.S.P.; Gamble, G.D.; Ling, L.H.; Sim, D.; Leong, K.T.G.; Yeo, P.S.D.; Ong, H.Y.; Jaufeerally, F.; Ng, T.P.; Cameron, V.A.; et al. Mortality Associated with Heart Failure with Preserved vs. Reduced Ejection Fraction in a Prospective International Multi-Ethnic Cohort Study. Eur. Heart J. 2018, 39, 1770–1780. [Google Scholar] [CrossRef]
  32. Loop, M.S.; van Dyke, M.K.; Chen, L.; Safford, M.M.; Kilgore, M.L.; Brown, T.M.; Durant, R.W.; Levitan, E.B. Low Utilization of Beta-Blockers Among Medicare Beneficiaries Hospitalized for Heart Failure With Reduced Ejection Fraction. J. Card. Fail. 2019, 25, 343–351. [Google Scholar] [CrossRef]
  33. Aiken, L.S.; West, S.G. Multiple Regression: Testing and Interpreting Interactions; Sage: Newbury Park, CA, USA, 1991. [Google Scholar]
  34. Daly, A.K.; Donaldson, P.T.; Bhatnagar, P.; Shen, Y.; Pe’er, I.; Floratos, A.; Daly, M.J.; Goldstein, D.B.; John, S.; Nelson, M.R.; et al. HLA-B*5701 Genotype Is a Major Determinant of Drug-Induced Liver Injury Due to Flucloxacillin. Nat. Genet. 2009, 41, 816–819. [Google Scholar] [CrossRef]
  35. Lanfear, D.E.; Luzum, J.A.; She, R.; Gui, H.; Donahue, M.P.; O’Connor, C.M.; Adams, K.F.; Sanders-van Wijk, S.; Zeld, N.; Maeder, M.T.; et al. Polygenic Score for Beta-Blocker Survival Benefit in European Ancestry Patients with Reduced Ejection Fraction Heart Failure. Circ. Heart Fail. 2020, 13, e007012. [Google Scholar] [CrossRef] [PubMed]
  36. Lanfear, D.E.; Luzum, J.A.; She, R.; Li, J.; Sabbah, H.N.; Zeld, N.; Liu, B.; Peterson, E.; Williams, L.K. Validation of a Polygenic Score for Beta-Blocker Survival Benefit in Patients With Heart Failure Using the United Kingdom Biobank. Circ. Genom. Precis. Med. 2023, 16, e003835. [Google Scholar] [CrossRef]
  37. Bis, J.C.; Sitlani, C.; Irvin, R.; Avery, C.L.; Smith, A.V.; Sun, F.; Evans, D.S.; Musani, S.K.; Li, X.; Trompet, S.; et al. Drug-Gene Interactions of Antihypertensive Medications and Risk of Incident Cardiovascular Disease: A Pharmacogenomics Study from the CHARGE Consortium. PLoS ONE 2015, 10, e0140496. [Google Scholar] [CrossRef]
  38. Magvanjav, O.; Gong, Y.; McDonough, C.W.; Chapman, A.B.; Turner, S.T.; Gums, J.G.; Bailey, K.R.; Boerwinkle, E.; Beitelshees, A.L.; Tanaka, T.; et al. Genetic Variants Associated With Uncontrolled Blood Pressure on Thiazide Diuretic/β-Blocker Combination Therapy in the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) and INVEST (International Verapamil-SR Trandolapril Study) Trials. J. Am. Heart Assoc. 2017, 6, e006522. [Google Scholar] [CrossRef]
  39. Singh, S.; Warren, H.R.; Hiltunen, T.P.; McDonough, C.W.; El Rouby, N.; Salvi, E.; Wang, Z.; Garofalidou, T.; Fyhrquist, F.; Kontula, K.K.; et al. Genome-Wide Meta-Analysis of Blood Pressure Response to Β1-Blockers: Results From ICAPS (International Consortium of Antihypertensive Pharmacogenomics Studies). J. Am. Heart Assoc. 2019, 8, e013115. [Google Scholar] [CrossRef] [PubMed]
  40. Shahin, M.H.; Conrado, D.J.; Gonzalez, D.; Gong, Y.; Lobmeyer, M.T.; Beitelshees, A.L.; Boerwinkle, E.; Gums, J.G.; Chapman, A.; Turner, S.T.; et al. Genome-Wide Association Approach Identified Novel Genetic Predictors of Heart Rate Response to β-Blockers. J. Am. Heart Assoc. 2018, 7, e006463. [Google Scholar] [CrossRef]
  41. Chang, S.-W.; McDonough, C.W.; Gong, Y.; Johnson, T.A.; Tsunoda, T.; Gamazon, E.R.; Perera, M.A.; Takahashi, A.; Tanaka, T.; Kubo, M.; et al. Genome-Wide Association Study Identifies Pharmacogenomic Loci Linked with Specific Antihypertensive Drug Treatment and New-Onset Diabetes. Pharmacogenomics J. 2018, 18, 106–112. [Google Scholar] [CrossRef]
  42. Bhatt, A.S.; DeVore, A.D.; DeWald, T.A.; Swedberg, K.; Mentz, R.J. Achieving a Maximally Tolerated β-Blocker Dose in Heart Failure Patients: Is There Room for Improvement? J. Am. Coll. Cardiol. 2017, 69, 2542–2550. [Google Scholar] [CrossRef] [PubMed]
Genes 14 02019 g001
Figure 1. QQ plots in self-identified black patients (A) and white patients (B). The solid red lines are the p-values expected by chance alone. The black points are the observed p-values for each locus tested in the GWAS.
Figure 1. QQ plots in self-identified black patients (A) and white patients (B). The solid red lines are the p-values expected by chance alone. The black points are the observed p-values for each locus tested in the GWAS.
Genes 14 02019 g001
Genes 14 02019 g002
Figure 2. Manhattan plot in black patients. The horizontal red line indicates p-value = 1 × 10−8. The horizontal blue line indicates p = 5 × 10−5. Variants of interest were defined as p < 5 × 10−5. The red circle indicates p = 1.3 × 10−11 for rs117032090 on chromosome 7.
Figure 2. Manhattan plot in black patients. The horizontal red line indicates p-value = 1 × 10−8. The horizontal blue line indicates p = 5 × 10−5. Variants of interest were defined as p < 5 × 10−5. The red circle indicates p = 1.3 × 10−11 for rs117032090 on chromosome 7.
Genes 14 02019 g002
Genes 14 02019 g003
Figure 3. Manhattan plot in white patients. The horizontal red line indicates p-value = 1 × 10−8. The horizontal blue line indicates p = 5 × 10−5. Variants of interest were defined as p < 5 × 10−5.
Figure 3. Manhattan plot in white patients. The horizontal red line indicates p-value = 1 × 10−8. The horizontal blue line indicates p = 5 × 10−5. Variants of interest were defined as p < 5 × 10−5.
Genes 14 02019 g003
Table 1. Baseline characteristics of the discovery (HFGR) and validation (HF-ACTION) datasets stratified by self-identified race.
Table 1. Baseline characteristics of the discovery (HFGR) and validation (HF-ACTION) datasets stratified by self-identified race.
CharacteristicsDiscovery Dataset: HFGRValidation Dataset: HF-ACTION
Black (n = 469)White (n = 459)1 pBlack (n = 288)White (n = 579)1 p
Female, n (%)193 (41%)141 (31%)0.001131 (45%)137 (22%)<0.001
Age, (years)65 ± 1271 ± 10<0.00154 ± 1361 ± 12<0.001
LVEF, (%)33 ± 1136 ± 10<0.00125 ± 825 ± 80.827
Ischemic etiology, n (%)153 (33%)254 (55%)<0.00189 (31%)357 (62%)<0.001
COPD, n (%)105 (22%)111 (24%)0.56924 (8.2%)85 (14%)0.019
Atrial fibrillation or flutter, n (%)89 (19%)169 (37%)<0.00146 (16%)153 (25%)0.002
Stroke, n (%)51 (11%)36 (7.8%)0.14147 (16%)48 (7.8%)<0.001
Diabetes, n (%)219 (47%)172 (37%)0.00599 (34%)185 (30%)0.288
Body mass index, (kg/m2)31 ± 731 ± 70.14033 ± 830 ± 6<0.001
SBP, (mmHg)132 ± 23126 ± 22<0.001115 ± 18115 ± 180.859
HR, (bpm)72 ± 1370 ± 130.01171 ± 1271 ± 110.688
NT pro-BNP, (pg/mL)2962 ± 32603170 ± 30680.3211266 ± 20091831 ± 24840.002
Serum creatinine, (mg/dL)1.39 ± 1.081.18 ± 0.61<0.0011.31 ± 0.661.33 ± 0.760.675
2 MAGGIC risk score17.6 ± 7.418.6 ± 7.00.02818.4 ± 5.720.7 ± 6.6<0.001
Beta-blocker exposure, (mg/day)0.22 ± 0.320.22 ± 0.320.9110.76 ± 0.470.60 ± 0.42<0.001
Length of follow-up, (days)1105 ± 6951118 ± 7070.780984 ± 374989 ± 3780.853
Deaths, n (%)112 (24%)116 (25%)0.67743 (15%)86 (15%)>0.999
1 p-values are for the comparison between self-identified black and white patients within each dataset. Bolded p-values indicate p < 0.05. 2 MAGGIC risk score was calculated without beta-blockers as an input variable. COPD—chronic obstructive pulmonary disease; HF-ACTION—Heart Failure: A Controlled Trial Investigating Outcomes of Exercise Training trial [17]; HFGR—HF pharmacogenomic registry; HR—heart rate; LVEF—left ventricular ejection fraction; MAGGIC—Meta-Analysis Global Group in Chronic Heart Failure risk score [21]; NT pro-BNP—N-terminal pro b-type natriuretic peptide; SBP—systolic blood pressure.
Table 2. Cox proportional hazards regression model results for the rs16844448 × beta-blocker exposure interaction term in black patients in HFGR and HF-ACTION.
Table 2. Cox proportional hazards regression model results for the rs16844448 × beta-blocker exposure interaction term in black patients in HFGR and HF-ACTION.
HFGR (n = 469)HF-ACTION (n = 288)
VariableCoeff.HR95% CIp-ValueCoeff.HR95% CIp-Value
rs16844448 × beta-blocker exposure interaction4.373.715.4–353.58.1 × 10−64.055.13.4–865.80.005
MAGGIC (excluding beta-blocker)0.11.11–1.2<0.0010.11.11–1.2<0.001
Beta-blocker propensity score0.051.00.9–1.20.6000.041.00.8–1.40.800
CI—confidence interval; Coeff.—coefficient; HF-ACTION—Heart Failure: A Controlled Trial Investigating Outcomes of Exercise Training [17]; HFGR—Heart failure pharmacogenomic registry [9]; HR—hazard ratio; MAGGIC—Meta-Analysis Global Group in Chronic Heart Failure risk score [21].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luzum, J.A.; Campos-Staffico, A.M.; Li, J.; She, R.; Gui, H.; Peterson, E.L.; Liu, B.; Sabbah, H.N.; Donahue, M.P.; Kraus, W.E.; et al. Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction. Genes 2023, 14, 2019. https://doi.org/10.3390/genes14112019

AMA Style

Luzum JA, Campos-Staffico AM, Li J, She R, Gui H, Peterson EL, Liu B, Sabbah HN, Donahue MP, Kraus WE, et al. Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction. Genes. 2023; 14(11):2019. https://doi.org/10.3390/genes14112019

Chicago/Turabian Style

Luzum, Jasmine A., Alessandra M. Campos-Staffico, Jia Li, Ruicong She, Hongsheng Gui, Edward L. Peterson, Bin Liu, Hani N. Sabbah, Mark P. Donahue, William E. Kraus, and et al. 2023. "Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction" Genes 14, no. 11: 2019. https://doi.org/10.3390/genes14112019

APA Style

Luzum, J. A., Campos-Staffico, A. M., Li, J., She, R., Gui, H., Peterson, E. L., Liu, B., Sabbah, H. N., Donahue, M. P., Kraus, W. E., Williams, L. K., & Lanfear, D. E. (2023). Genome-Wide Association Study of Beta-Blocker Survival Benefit in Black and White Patients with Heart Failure with Reduced Ejection Fraction. Genes, 14(11), 2019. https://doi.org/10.3390/genes14112019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop