β-Adrenergic Receptor Gene Polymorphisms and β-Blocker Treatment Outcomes in Hypertension

doi:10.1038/clpt.2008.139

Clin Pharmacol Ther. Author manuscript; available in PMC 2009 December 1.

Published in final edited form as:

Clin Pharmacol Ther. 2008 December; 84(6): 715–721.

Published online 2008 July 9. doi: 10.1038/clpt.2008.139

PMCID: PMC2675574

NIHMSID: NIHMS109967

β-Adrenergic Receptor Gene Polymorphisms and β-Blocker Treatment Outcomes in Hypertension

MA Pacanowski,¹ Y Gong,¹ RM Cooper-DeHoff,² NJ Schork,³ MD Shriver,⁴ TY Langaee,¹ CJ Pepine,² and JA Johnson^1,², For the INVEST Investigators

¹Department of Pharmacy Practice and Center for Pharmacogenomics, University of Florida College of Pharmacy, Gainesville, Florida, USA

²Division of Cardiovascular Medicine, Department of Medicine, University of Florida College of Medicine, Gainesville, Florida, USA

³Department of Molecular and Experimental Medicine, Scripps Research Institute, La Jolla, California, USA

⁴Department of Anthropology, Pennsylvania State University, University Park, Pennsylvania, USA.

Correspondence: JA Johnson (Email: johnson/at/cop.ufl.edu)

The publisher's final edited version of this article is available at Clin Pharmacol Ther

See other articles in PMC that cite the published article.

Abstract

Numerous studies have demonstrated that β₁- and β₂-adrenergic receptor gene (ADRB1 and ADRB2) variants influence cardiovascular risk and β-blocker responses in hypertension and heart failure. We evaluated the relationship between ADRB1 and ADRB2 haplotypes, cardiovascular risk (death, nonfatal myocardial infarction (MI ), and nonfatal stroke), and atenolol-based vs. verapamil sustained-release (SR )-based antihypertensive therapy in 5,895 coronary artery disease (CAD) patients. After an average of 2.8 years, death rates were higher in patients carrying the ADRB1 Ser49-Arg389 haplotype (hazard ratio (HR ) 3.66, 95% confidence interval (95% CI) 1.68–7.99). This mortality risk was significant in patients randomly assigned to verapamil SR (HR 8.58, 95% CI 2.06–35.8) but not atenolol (HR 2.31, 95% CI 0.82–6.55), suggesting a protective role for the β-blocker. ADRB2 haplotype associations were divergent within the treatment groups but did not remain significant after adjustment for multiple comparisons. ADRB1 haplotype variation is associated with mortality risk, and β-blockers may be preferred in subgroups of patients defined by ADRB1 or ADRB2 polymorphisms.

Cardiovascular disease is the leading cause of morbidity and mortality in the United States.¹ Evidence is accumulating that genetic polymorphisms may be predictive of cardiovascular risk.²^,³ Moreover, several studies have identified genetic factors that influence blood pressure and metabolic responses to β-blockers, thiazide diuretics, and renin–angiotensin system antagonists.⁴ Whether such pharmacogenetic differences translate to differences in the clinical outcome of antihypertensive therapy is less clear, particularly when patients receive multiple drugs that are titrated to a target blood pressure.⁵ A pharmacogenetic approach to treating hypertension could not only reduce the number and cost of medications but also reduce morbidity and mortality if the outcome of drug treatment differs by genotype.

Single-nucleotide polymorphisms (SNPs) in the genes encoding the adrenergic receptors have functional and physiological consequences.⁶ β₁-adrenergic receptors are important in regulating heart rate and contractility as well as renin release in the kidney. Nonsynonymous SNPs in the β₁-adrenergic receptor gene (ADRB1; Ser49Gly and Arg389Gly) result in differences in agonist-mediated downregulation and signaling activity.⁷ Similarly, β₂-adrenergic receptors play an important role in cardiac function, metabolism, and vascular tone. Nonsynonymous SNPs in the β₂-adrenergic receptor gene (ADRB2; Arg16Gly and Gln27Glu) affect expression, agonist responsiveness, and desensitization.⁸^,⁹ Previous epidemiologic studies have suggested that ADRB1 and ADRB2 SNPs correlate with a variety of intermediate cardiovascular phenotypes,¹⁰ cardiovascular risk,¹¹^,¹² and β-blocker responses in hypertension¹³^,¹⁴ and heart failure.¹⁵^,¹⁶

Taken together, genetically determined differences in β-adrenergic receptor function may influence cardiovascular risk as well as the degree of risk reduction conferred by drugs that inhibit adrenergic activity. Therefore, we investigated the influence of ADRB1 and ADRB2 haplotype variation on the incidence of death, nonfatal myocardial infarction (MI), and nonfatal stroke as well as the pharmacogenetics of β-blocker– and calcium channel blocker–based antihypertensive therapy in the INternational VErapamil SR/Trandolapril STudy— GENEtic Substudy (INVEST-GENES). Based on the literature, we hypothesized that patients with the ADRB1 Ser49-Arg389 haplotype or the ADRB2 haplotype containing the Arg16 and Gln27 alleles would be at relatively higher risk for cardiovascular events¹¹^,¹² and that atenolol would be beneficial as compared with sustained-release verapamil (verapamil SR) in patients with the ADRB1 Ser49-Arg389 haplotype.

RESULTS

Study population and baseline characteristics

INVEST-GENES consisted of an ethnically diverse, elderly population of patients with hypertension and documented coronary artery disease (CAD), including a large proportion of women and individuals with diabetes. Patients were randomly assigned to receive atenolol or verapamil SR. Trandolapril and hydrochlorothiazide were added as needed to control blood pressure. Demographic and clinical characteristics did not differ significantly by treatment strategy (Table 1). Haplotypes for at least one gene were inferred for a total of 5,895 patients, including 5,817 who were successfully genotyped at both ADRB1 loci (Ser49Gly and Arg389Gly) and 5,877 who were successfully genotyped at two of three ADRB2 loci (Gly16Arg, Gln27Glu, and 523C>A). Genotype concordance was >97% at all loci. All loci were in Hardy–Weinberg equilibrium except ADRB2 Gly16Arg in Hispanic patients (χ² = 6.12, P < 0.05). The minor deviation from Hardy–Weinberg equilibrium in Hispanics appeared to be due to admixture, as this was the most admixed population in our sample. In Hispanics with >80% European ancestry, there was no evidence of deviation from Hardy–Weinberg equilibrium (P = 0.47), which was selected as the cutoff because the black and white subjects in this sample have, on average, >80% ancestry from a single continental population. Low to moderate levels of linkage disequilibrium were noted between the SNPs in ADRB1 (r² = 0.04–0.19) and ADRB2 (r² = 0.07–0.45). The two variant loci in ADRB1 formed three common haplotypes, and the three variant loci in ADRB2 also formed three common haplotypes (Table 2). The genotype and haplotype distributions differed significantly by race/ ethnicity (Table 2).

Table 1

Baseline characteristics

Table 2

Minor allele and haplotype frequencies by race/ethnicity

The average follow-up period for the primary outcome was 2.8 ± 0.7 years. The treatment strategies did not differ in terms of the primary outcome (defined as the first occurrence of death, nonfatal MI, or nonfatal stroke; hazard ratio (HR) 0.81, 95% confidence interval (95% CI) 0.63–1.05, P = 0.12). The median doses of atenolol and verapamil SR at the end of the study were 100 and 240 mg/day, respectively. Hydrochlorothiazide was used in 74% of atenolol-treated patients and 61% of verapamil SR–treated patients, whereas trandolapril was used in 70% of atenolol-treated patients and 78% of verapamil SR–treated patients. Mean on-treatment systolic blood pressure values were not significantly different across treatment strategies and haplotype groups (data not shown). Mean on-treatment diastolic blood pressure was statistically lower in the verapamil SR strategy group (77.7 mm Hg vs. 78.2 mm Hg, P = 0.003), patients with the ADRB1 Ser49-Arg389 haplotype (0 copies, 78.4 mm Hg; 1 copy, 77.9 mm Hg; 2 copies, 77.6 mm Hg; P = 0.003), and patients with the ADRB2 Gly16-Glu27-523C haplotype (0 copies, 78.4 mm Hg; 1 copy, 77.8 mm Hg; 2 copies, 76.6 mm Hg; P < 0.0001). In all cases the differences were <2 mm Hg and therefore not likely to be clinically important. Mean on-treatment heart rate did not differ according to ADRB1 haplotype and was statistically, but not clinically, lower in patients with the ADRB2 Gly16-Glu27-523C haplotype (0 copies, 72.2 beats/min; 1 copy, 71.8 beats/min; 2 copies, 71.7 beats/min; P = 0.005).

Associations of ADRB1 with primary and secondary outcomes

The model based on the Ser49-Arg389 haplotype best characterized the risk for the primary outcome (log-rank P = 0.02), whereas the other two common haplotypes, Ser49-Gly389 and Gly49-Arg389, were not associated with the primary outcome (Supplementary Table S1 online). The haplotype risk was similar between patients with one or two copies of the Ser49-Arg389 (Figure 1), and these patients had a relatively higher risk for the primary outcome relative to non-carriers (Ser49-Arg389 carriers vs. noncarriers, HR 1.51, 95% CI 1.07–2.12, P = 0.02). This association was driven entirely by mortality (Ser49-Arg389 carriers vs. noncarriers, HR 3.66, 95% CI 1.68–7.99, P = 0.001), whereas no association was noted for nonfatal MI or nonfatal stroke (Figure 1). The cause of death was adjudicated, although the exact cause of noncardiovascular death was not specified. The Ser49-Arg389 haplotype was significantly associated with both cardiovascular and noncardiovascular mortality (data not shown). The increased risk of death among carriers of the Ser49-Arg389 haplotype was consistent across racial/ethnic groups and after adjustment for ancestry (Ser49-Arg389 carriers vs. noncarriers: whites, HR 2.32, 95% CI 0.84–6.44, P = 0.1; Hispanics, HR 2.95, 95% CI 0.89–9.78, P = 0.08; blacks, HR not calculated because all events occurred in carriers). Given the similarity in effect across racial/ethnic groups and the lack of apparent population stratification, all subsequent analyses are presented for the combined population to maintain power, controlling for race/ethnicity. In genotype-based analyses, the mortality association was detectable for Arg389Gly under a recessive model (P = 0.03); the data suggest that ADRB1 haplotype is more informative than the individual SNPs.

Figure 1

Associations of the ADRB1 Ser49-Arg389 haplotype with primary and secondary outcomes. Hazard ratios are based on reduced model adjusted for age, sex, and race/ethnicity. *Crude incidence per 1,000 patient years. MI, myocardial infarction.

ADRB2 associations with primary outcomes

None of the ADRB2 haplotype models revealed differential risk for the primary outcome in the overall population, with the lowest P value being 0.39 (Supplementary Table S1 online). Given the lack of statistical trends for a main effect for the primary outcome, secondary outcomes were not tested.

Pharmacogenetic associations with ADRB1 and ADRB2

Because genetic associations for ADRB1 with the primary outcome were driven by differences in death, pharmacogenetic analyses focused on this outcome. The increase in mortality risk among patients with one or two copies of the Ser49-Arg389 haplotype was significant in patients randomly assigned to verapamil SR (HR 8.58, 95% CI 2.06–35.8, P = 0.003) but not in patients assigned to atenolol (HR 2.31, 95% CI 0.82–6.55, P = 0.11). These data suggest that atenolol offsets the risk associated with the ADRB1 haplotype (P_interaction = 0.19; Figure 2). The point estimates were minimally affected when adjusted for clinical covariates and secondary drug use (data not shown). Adjustment for ancestry did not change the magnitude or direction of any of these associations (data not shown). These associations remained significant even when using the Bonferroni-corrected P value. Blood pressure and heart rate values were similar across haplotype and drug use strata (Figure 2). The median atenolol and verapamil doses did not differ according to ADRB1 haplotype, although Ser49-Arg389 carriers were less likely than noncarriers to receive a second-line drug in either strategy (atenolol + hydrochlorothiazide 73% vs. 78%, P = 0.007; verapamil SR + trandolapril 69% vs. 75%, P = 0.006).

Figure 2

All-cause mortality and mean on-treatment blood pressure by ADRB1 Ser49-Arg389 haplotype and atenolol/verapamil sustained-release (SR) therapy. HR, hazard ratio; 95% CI, 95% confidence interval; S49-R389, Ser49-Arg389 haplotype; AT, atenolol; VE, verapamil (more ...)

Pharmacogenetic analysis of common ADRB2 haplotypes revealed that the risk for the primary outcome differed significantly across the Gly16-Glu27-523C haplotype in verapamil SR–treated patients but not in atenolol-treated patients (P_interaction = 0.05; Figure 3). None of the other haplotype models was associated with the primary outcome within the treatment strategies. The haplotype-associated risks were driven largely by mortality (P_interaction = 0.11) and nonfatal MI (P_interaction = 0.06; Figure 3). On the basis of the Bonferroni-corrected P value, none of the subgroup associations would be defined as significant, despite the significant P values for the interaction term. Genotype-based analysis revealed that this was driven largely by the 523C>A genotype, but similar trends that were consistent with the haplotype association were noted for the Gly16Arg and Gln27Glu SNPs (data not shown), again supporting the view that the haplotype analysis is the more powerful one. The median atenolol and verapamil SR doses did not differ by ADRB2 haplotype; nor did the rates of hydrochlorothiazide or trandolapril use (data not shown).

Figure 3

Primary and secondary outcomes by ADRB2 Gly16-Glu27-523C haplotype and antihypertensive drug therapy. *Crude incidence per 1,000 patient years. †Includes only patients ever exposed to atenolol (94% in atenolol strategy). ‡Hazard ratios (more ...)

As an exploratory analysis, mortality was modeled for the randomized drugs on the basis of both ADRB1 and ADRB2 haplotype information. The analysis revealed that patients with at least one copy of the Ser49-Arg389 haplotype and zero copies of the Gly16-Glu27-523C haplotype (representing 42% of the study population) had better outcomes when treated with atenolol than with verapamil SR (HR 0.42, 95% CI 0.21–0.82, P = 0.01). Comparing this result to the HR of 0.64 when considering the ADRB1 gene alone suggests that a consideration of both genes may be even more informative for identifying those most likely to benefit from β-blocker therapy.

DISCUSSION

Amino-acid changing SNPs in the adrenergic receptors have previously been shown to be associated with cardiovascular and drug response phenotypes. In treated hypertensive patients with CAD, we identified a significant association between SNPs in ADRB1 and the incidence of all-cause death, showing that patients with haplotypes bearing the common alleles (Ser49 and Arg389) were at relatively higher risk. The risk associated with this haplotype was significantly reduced by β-blocker therapy, but not by calcium channel blocker therapy, a finding consistent with the understanding of the functional consequences of these genetic polymorphisms. Specifically, the Arg389 form of the receptor has been documented in several in vitro and ex vivo studies to be associated with increased coupling of the β₁-adrenergic receptor to G protein, leading to greater adenylyl cyclase activation.¹⁷^,¹⁸ The Ser49 form of the receptor has most consistently been associated with resistance to receptor downregulation.¹⁹^,²⁰ Therefore, the Ser49-Arg389 haplotype would be expected to be most responsive to activation by catecholamines, and consequently a greater response to β-blockade with this haplotype would also be expected.

Consistent with our primary hypothesis, patients with at least one copy of the haplotype containing the wild-type Ser49 and Arg389 alleles were at relatively higher risk for the primary outcome than patients with variant alleles. The main haplotype association was driven by more than a threefold difference in the rate of all-cause mortality.²¹ Our findings are in line with those of Iwai et al.,²² who reported an association between the Arg389 allele and MI, but other observational studies have not supported a major influence of ADRB1 variants on cardiovascular risk.¹²^,²³ The lack of consistency may be attributable to differences in patient populations studied (i.e., those with or without overt CAD)²³ or prevalent use of β-blockers.¹² Our data support the use of β-blockers as being a potential confounder for detecting the genetic association with outcomes, given that there were no differences in outcomes by genotype among the atenolol-treated patients. Overall, the results of the current investigation are in line with the known functionality of these variants and the widely recognized adverse consequences of chronic sympathetic activation.

The pharmacogenetic findings regarding the Ser49-Arg389 haplotype are consistent not only with the in vitro and ex vivo studies documenting the functional basis of these polymorphisms but also with the existing β-blocker pharmacogenetics literature. ²⁴^,²⁵ This study extends those findings to include outcomes of antihypertensive therapy. In our investigation, patients carrying the Ser49-Arg389 haplotype derived a significant survival benefit from the use of a β-blocker. We previously demonstrated that hypertensive patients who were Ser49-Arg389 homozygotes experienced a significantly greater blood pressure response to β-blockers than those with haplotypes containing a variant allele.¹³ With variations on this theme, four other studies have corroborated the finding that patients with the Arg389 allele or Ser49-Arg389 haplotype show a greater blood pressure response to β-blockers.⁸^,¹⁴^,²⁶^,²⁷ It follows that, in certain subpopulations, hypertension may have a strong adrenergic component that is particularly amenable to β-blocker therapy. Similarly, studies in heart failure patients have also suggested the Arg389 allele or Ser49-Arg389 haplotype is associated with the greatest improvement in ejection fraction after initiation of β-blocker therapy.¹⁵^,²⁸^,²⁹ Liggett et al. recently demonstrated that Arg389 homozygous heart failure patients derived a significant survival benefit from bucindolol as compared with placebo, whereas in Gly389 carriers the outcomes with bucindolol and placebo were similar.³⁰

Taken together, these studies suggest that the ADRB1 gene is an important pharmacogenetic target for β-blocker response. The literature suggests that polymorphisms in the gene influence intermediate-response phenotypes (e.g., blood pressure reduction and ejection fraction improvement) along with mortality outcomes in hypertension and heart failure. The specific findings from this study suggest that β-blockers may be the preferred antihypertensive therapy in hypertensive CAD patients who are Ser49-Arg389 carriers.

It should be noted that carriers of the major alleles (Ser49 and Arg389) were at relatively increased risk and benefited from β-blocker therapy, and, conversely, those with a variant on both chromosomes were at lower risk. This is in contrast to the apparent inheritance patterns of some of the other β-blocker pharmacogenetics studies for ADRB1 described above, in which a dominant model was typically assumed, with Ser49-Arg389 homozygotes exhibiting the greatest β-blocker response and carriers of at least one variant allele/haplotype having a lesser response. However, very few studies have been sufficiently powered to test for mode of inheritance, given that variant homozygotes, typically, have been minimally represented. In these cases a dominant model (whereby variant carriers were often collapsed to a single group) was typically pursued out of statistical necessity. Thus, although the mode of inheritance of this association seems to be in contrast to some of the β-blocker pharmacogenetics literature, it can also be concluded that the literature has not revealed a clear pattern of inheritance for this gene and the resulting phenotypes. This is the largest ADRB1 pharmacogenetic study to date and, from this perspective, may have been the best powered to assess for mode of inheritance. Additionally, with complex phenotypes, it is possible that different phenotypes will exhibit different inheritance patterns.¹⁵^,³⁰

We also identified a statistical interaction between atenolol and verapamil SR and ADRB2 haplotypes. However, this did not meet our threshold for significance after adjusting for multiple comparisons. This association was driven primarily by divergent risks for both death and nonfatal MI, and none of the associations appeared to be driven by differences in blood pressure. If the findings for ADRB2 are validated, knowledge of this haplotype may further enhance the ability to identify patients who might benefit from β-blocker therapy.

Limitations

As a cohort study nested in a randomized trial with adjudicated end points, INVEST-GENES has several advantages over population-based studies, and the results are generalizable to other CAD populations managed with contemporary interventions. However, the current investigation also has limitations that deserve consideration. First, despite the large size of the study population, the event rate was low, and this study may have been underpowered to evaluate gene–drug interactions, particularly for the individual end points and within certain subgroups of patients (i.e., individual end point by drug by gene by race). Second, the INVEST-GENES population is racially/ ethnically diverse. To control for potential confounding by population stratification, we considered analyses separately by race and by inclusion of ancestry informative marker data. These various analyses suggest that our findings are not confounded by population stratification. Third, the use of trandolapril and hydrochlorothiazide was, by design, different in the two treatment strategies and may have influenced the results. However, the findings were similar in the expanded model that adjusted for exposure to these drugs. Last, replication in independent cohorts is desirable. In the absence of another study with randomized drug therapy and comparably rigorous follow-up and phenotype definition, we must rely on the existing evidence from the laboratory and endophenotype studies.³¹

Conclusion

Identifying genetic markers for cardiovascular risk has the potential to improve cardiovascular risk stratification and identify those requiring more aggressive management of hypertension and related chronic diseases. Common SNPs in the genes encoding the β₁- and β₂-adrenergic receptors alter receptor activity and have physiological consequences. Consistent with the known functionality of the β₁-adrenergic receptor variants, we identified an association between ADRB1 haplotypes and the risk of death. More importantly, our data suggest that β-blockers offset this mortality risk, in keeping with observations that patients bearing the wild-type alleles are more responsive to β-blocker therapy in settings of blood pressure lowering, improvement in ejection fraction, and survival in heart failure. ADRB2 variants were similarly associated with treatment outcomes, but given the inconsistencies in the literature, these findings require independent replication. The pharmacogenetic evidence for β-blockers and adrenergic receptor genes is highly convincing, particularly for ADRB1, and our data suggest that a patient’s genotype could influence the antihypertensive drug choice independent of blood pressure responses.

METHODS

INVEST-GENES design and participants

INVEST was a prospective, randomized, open-label, blinded-end point (PROBE) trial designed to compare antihypertensive treatment outcomes in 22,576 patients. The INVEST-GENES cohort consisted of 5,979 patients from 184 sites in the United States and Puerto Rico who provided DNA samples and additional written informed consent for genetic studies. The details of the INVEST methods and main outcomes have been previously reported.³² Briefly, INVEST enrolled hypertensive patients over the age of 50 years with documented CAD. Patients were randomly assigned to a verapamil SR- or an atenolol-based treatment strategy. Trandolapril was recommended for all patients with heart failure, renal dysfunction, or diabetes. Hydrochlorothiazide and/or trandolapril were added as needed to achieve blood pressure targets defined in the sixth report of the Joint National Committee on Prevention, Detection, Evaluation, and Treatment of High Blood Pressure. Details on the addition of study drugs and dose titration may be found in the original INVEST publication. ³² Patients were followed every 6 weeks for the first 6 months and every 6 months thereafter until 2 years after the last patient was enrolled. Blood pressure control and cardiovascular outcomes were similar between the treatment strategies in the main trial.³²

Outcomes

The primary outcome was a composite of the first occurrence of all-cause mortality, nonfatal MI, or nonfatal stroke. Secondary outcomes included the individual components of the primary outcome. Events were adjudicated by an independent committee that was blinded to treatment strategy.

Genotyping

Buccal tissue samples were obtained by mouthwash, and genomic DNA was isolated using the Gentra Systems PureGene kit. Patients were genotyped for two variants in ADRB1 (Ser49Gly (145A>G), rs1801252; and Arg389Gly (1165C>G), rs1801253) and three variants in ADRB2 (Gly16Arg (46G>A), rs1042713; Gln27Glu (79A>G), rs1042714; and Arg175Arg (523C>A), rs1042718) using pyrosequencing (Biotage, Uppsala, Sweden) and TaqMan allelic discrimination (Applied Biosystems, Foster City, CA). Genotype accuracy was verified by genotyping 5–10% randomly selected duplicate samples for each SNP on the alternate platform. Ancestry informative markers (87 total) were genotyped using either allele-specific PCR with universal energy transfer labeled primers or competitive allele-specific PCR at Prevention Genetics (Marshfield, WI).

Statistics

Hardy–Weinberg equilibrium was tested for each racial/ethnic group using χ² analysis. Haplotypes were reconstructed separately for each racial/ethnic group from raw genotype data using PHASE software (version 2.1)³³ for patients who were successfully genotyped at both ADRB1 loci and at least two ADRB2 loci. Each haplotype was coded according to the number of copies (zero, one, or two). Linkage disequilibrium (r²) between the SNPs in each racial/ethnic group was estimated using Haploview.³⁴ Demographic, baseline clinical characteristics, and mean on-treatment (from 6 months to the last follow-up) blood pressure levels and heart rates were compared by haplotype to identify potentially confounding factors, using χ² tests for categorical data and t-tests, analysis of variance, or a nonparametric equivalent for continuous data.

Given that INVEST-GENES was a cohort of patients from a randomized clinical trial, that the candidate genes encode the drug targets, and that the hypotheses are pharmacogenetic in nature (in an effort to validate published associations), we screened for primary outcome associations both in the overall population and within the randomized treatment arms. The main effects of each haplotype were first evaluated in the intention-to-treat population for the primary outcome based on haplotype copy number (zero, one, or two) using Kaplan–Meier analysis with pooled log-rank tests. Cox proportional hazards regression was performed to estimate HR and 95% CI for each copy of the haplotype relative to zero copies of the haplotype. The regression model was initially adjusted for race/ethnicity, age, sex, and treatment strategy (reduced model). The following covariates were subsequently entered into the model using the stepwise procedure if P < 0.1 and retained if P < 0.05 (expanded model): history of heart failure, MI, diabetes, stroke or transient ischemic attack, renal insufficiency, dyslipidemia, left ventricular hypertrophy, peripheral vascular disease, stable angina, unstable angina, arrhythmia, cancer, or ever having smoked; body mass index; and baseline systolic and diastolic blood pressure values. The reduced model is presented unless the point estimates differ from the expanded model by more than 0.1. The threshold for significance in the screening analysis was set at P < 0.05 based on the a priori genetic and pharmacogenetic hypotheses that follow the existing clinical and experimental data in the literature for both genes and the studied polymorphisms. The proportionality of haplotype effects was evaluated by examining Schoenfeld residuals.

To reduce multiple comparisons, only significant haplotype associations with the primary outcome were followed by analysis of the secondary end points and pharmacogenetic relationships. To account for comparisons incurred in the stepwise analysis approach, the final significant P value for log-rank tests was defined as P ≤ 0.003, which was arrived at on the basis of a Bonferroni correction for 15 comparisons as follows: ADRB1: one gene test for primary outcome (one in overall population), one haplotype tested for three secondary outcomes, one haplotype tested for one secondary outcome in two treatment strategies (six tests for ADRB1 in total), and ADRB2: three gene tests for primary outcome (one in overall population and two within treatment strategies), one haplotype tested for three secondary outcomes in two treatment strategies (nine tests for ADRB2 in total). This correction might be viewed as being overly conservative, given that an assumption of the Bonferroni correction is independence of all the tests, and these tests were not all independent. Additionally, the most appropriate inheritance model (i.e., additive, dominant, or recessive) was selected on the basis of visual inspection of per-allele point estimates in an effort to maintain power in these analyses. Genotype-based analyses were considered exploratory.

Verapamil SR and atenolol were started at baseline in 100 and 94% of patients in the respective strategies. Consequently, pharmacogenetic analyses examined haplotype risks in the intention-to-treat population for patients in whom verapamil SR or atenolol was ever prescribed. Associations were evaluated by Kaplan–Meier analysis with pairwise log-rank tests and by Cox proportional hazard regression stratified by exposure or with haplotype–drug interaction terms. Given the differences in time to initiation and the overlapping use of trandolopril and hydrochlorothiazide, the expanded model for pharmacogenetic analyses also adjusted for average trandolapril and hydrochlorothiazide doses as time-varying covariates.

To control for the potential of population stratification in our racially and ethnically diverse population, we used a total panel of 87 ancestry informative markers, selected to show large allele-frequency differences across three parental populations (West Africans, Native Americans, and Europeans) from a large panel of more than 10,000 SNPs.³⁵ The maximum-likelihood criterion was then used to estimate each patient’s individual genomic ancestry proportions on these three axes, and these terms were included in the statistical models for patients with data for at least 50% of the ancestry informative markers. All analyses were performed with and without stratification by race/ethnicity because of reduced power in the racial/ethnic subgroup analysis. Data are presented for the combined population unless the results differed across racial/ ethnic strata.

All statistical analyses were performed using SAS version 9.1 (SAS Institute, Cary, NC).

Supplementary Material

Table 1: SUPPLEMENTARY MATERIAL

is linked to the online version of the paper at http://www.nature.com/cpt

Click here to view.^{(55K, doc)}

ACKNOWLEDGMENTS

We thank Lynda Stauffer and Kathy Eberst (posthumously) for processing and genotyping the samples, Robert Kolb for coordinating INVEST sites for INVEST-GENES participation, and all the INVEST patients who also agreed to participate in the genetic substudy. The phenotype and genotype data are available at http://www.PharmGKB.org; PS203971, PS205525, PS205530, PS203972, PS203974, PS205526, PS205527, PS205529, PS206277, PS207726, PS207727, PS207729, PS207730, PS207731, PS207732, PS207733, PS207734, and PS207735. This project was funded by National Institutes of Health grants HL074730, HL69758, and RR017568; a grant from Abbott Pharmaceuticals; and American Heart Association Postdoctoral Fellowship 0625619B. Clinical trials registration: http://clinicaltrials.gov, identifier NCT00133692.

Footnotes

Conflict of interest

Drs Johnson, Pepine, Cooper-DeHoff, and Langaee received grant funding from Abbott Laboratories. Drs Cooper-DeHoff and Pepine along with the University of Florida hold U.S. patent 5,991,731 related to INVEST. Dr Pepine has been a consultant for Abbott Laboratories. The other authors declared no conflict of interest.

References

1. Rosamond W, et al. Heart disease and stroke statistics—2007 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation. 2006;115:e69–e171. [PubMed]

2. Ginsburg GS, Donahue MP, Newby LK. Prospects for personalized cardiovascular medicine: the impact of genomics. J. Am. Coll. Cardiol. 2005;46:1615–1627. [PubMed]

3. Samani NJ, et al. Genomewide association analysis of coronary artery disease. N. Engl. J. Med. 2007;357:443–453. [PMC free article] [PubMed]

4. Johnson JA, Turner ST. Hypertension pharmacogenomics: current status and future directions. Curr. Opin. Mol. Ther. 2005;7:218–225. [PubMed]

5. Davis BR, et al. Antihypertensive therapy, the alpha-adducin polymorphism, and cardiovascular disease in high-risk hypertensive persons: the Genetics of Hypertension-Associated Treatment Study. Pharmacogenomics J. 2007;7:112–122. [PubMed]

6. Kirstein SL, Insel PA. Autonomic nervous system pharmacogenomics: a progress report. Pharmacol. Rev. 2004;56:31–52. [PubMed]

7. Sandilands A, Yeo G, Brown MJ, O’Shaughnessy KM. Functional responses of human beta1 adrenoceptors with defined haplotypes for the common 389R>G and 49S>G polymorphisms. Pharmacogenetics. 2004;14:343–349. [PubMed]

8. Dishy V, et al. The effect of common polymorphisms of the beta2-adrenergic receptor on agonist-mediated vascular desensitization. N. Engl. J. Med. 2001;345:1030–1035. [PubMed]

9. Drysdale CM, et al. Complex promoter and coding region beta 2-adrenergic receptor haplotypes alter receptor expression and predict in vivo responsiveness. Proc. Natl. Acad. Sci. USA. 2000;97:10483–10488. [PubMed]

10. Bengtsson K, et al. Polymorphism in the beta(1)-adrenergic receptor gene and hypertension. Circulation. 2001;104:187–190. [PubMed]

11. Heckbert SR, et al. Beta2-adrenergic receptor polymorphisms and risk of incident cardiovascular events in the elderly. Circulation. 2003;107:2021–2024. [PubMed]

12. Lanfear DE, Jones PG, Marsh S, Cresci S, McLeod HL, Spertus JA. Beta2-adrenergic receptor genotype and survival among patients receiving beta-blocker therapy after an acute coronary syndrome. JAMA. 2005;294:1526–1533. [PubMed]

13. Johnson JA, Zineh I, Puckett BJ, McGorray SP, Yarandi HN, Pauly DF. Beta 1-adrenergic receptor polymorphisms and antihypertensive response to metoprolol. Clin. Pharmacol. Ther. 2003;74:44–52. [PubMed]

14. Liu J, et al. Beta1-Adrenergic receptor polymorphisms influence the response to metoprolol monotherapy in patients with essential hypertension. Clin. Pharmacol. Ther. 2006;80:23–32. [PubMed]

15. Mialet Perez J, et al. Beta 1-adrenergic receptor polymorphisms confer differential function and predisposition to heart failure. Nat. Med. 2003;9:1300–1305. [PubMed]

16. Terra SG, et al. Beta1-adrenergic receptor polymorphisms and left ventricular remodeling changes in response to beta-blocker therapy. Pharmacogenet. Genomics. 2005;15:227–234. [PubMed]

17. Liggett SB. Beta2-adrenergic receptor polymorphisms and sudden cardiac death: a signal to follow. Circulation. 2006;113:1818–1820. [PubMed]

18. Mason DA, Moore JD, Green SA, Liggett SB. A gain-of-function polymorphism in a G-protein coupling domain of the human beta 1-adrenergic receptor. J. Biol. Chem. 1999;274:12670–12674. [PubMed]

19. Levin MC, Marullo S, Muntaner O, Andersson B, Magnusson Y. The myocardium-protective Gly-49 variant of the beta 1-adrenergic receptor exhibits constitutive activity and increased desensitization and downregulation. J. Biol. Chem. 2002;277:30429–30435. [PubMed]

20. Rathz DA, Brown KM, Kramer LA, Liggett SB. Amino acid 49 polymorphisms of the human beta1-adrenergic receptor affect agonistpromoted trafficking. J. Cardiovasc. Pharmacol. 2002;39:155–160. [PubMed]

21. A randomized trial of propranolol in patients with acute myocardial infarction. I. Mortality results. JAMA. 1982;247:1707–1714. [PubMed]

22. Iwai C, et al. Arg389Gly polymorphism of the human beta1-adrenergic receptor in patients with nonfatal acute myocardial infarction. Am. Heart. J. 2003;146:106–109. [PubMed]

23. White HL, et al. An evaluation of the beta-1 adrenergic receptor Arg389Gly polymorphism in individuals at risk of coronary events. A WOSCOPS substudy. Eur. Heart. J. 2002;23:1087–1092. [PubMed]

24. Pacanowski MA, Johnson JA. PharmGKB submission update: IX. ADRB1 gene summary. Pharmacol. Rev. 2007;59:2–4. [PMC free article] [PubMed]

25. Shin J, Johnson JA. Pharmacogenetics of beta-blockers. Pharmacotherapy. 2007;27:874–887. [PMC free article] [PubMed]

26. Liu J, et al. Gly389Arg polymorphism of beta1-adrenergic receptor is associated with the cardiovascular response to metoprolol. Clin. Pharmacol. Ther. 2003;74:372–379. [PubMed]

27. Karlsson J, et al. Beta1-adrenergic receptor gene polymorphisms and response to beta1-adrenergic receptor blockade in patients with essential hypertension. Clin. Cardiol. 2004;27:347–350. [PubMed]

28. Terra SG, et al. Association between beta-adrenergic receptor polymorphisms and their G-protein-coupled receptors with body mass index and obesity in women: a report from the NHLBI-sponsored WISE study. Int. J. Obes. (Lond) 2005;29:746–754. [PubMed]

29. Yu WP, Lou M, Deng B, Song HM, Wang HB. [Beta1-adrenergic receptor (Arg389Gly) polymorphism and response to bisoprolol in patients with chronic heart failure] Zhonghua Xin Xue Guan Bing Za Zhi. 2006;34:776–780. [PubMed]

30. Liggett SB, et al. A polymorphism within a conserved beta(1)-adrenergic receptor motif alters cardiac function and beta-blocker response in human heart failure. Proc. Natl. Acad. Sci. USA. 2006;103:11288–11293. [PubMed]

31. Chanock SJ, et al. Replicating genotype-phenotype associations. Nature. 2007;447:655–660. [PubMed]

32. Pepine CJ, et al. A calcium antagonist vs a non-calcium antagonist hypertension treatment strategy for patients with coronary artery disease. The International Verapamil-Trandolapril Study (INVEST): a randomized controlled trial. JAMA. 2003;290:2805–2816. [PubMed]

33. Pritchard JK, Stephens M, Rosenberg NA, Donnelly P. Association mapping in structured populations. Am. J. Hum. Genet. 2000;67:170–181. [PubMed]

34. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. [PubMed]

35. Shriver MD, et al. Large-scale SNP analysis reveals clustered and continuous patterns of human genetic variation. Hum. Genomics. 2005;2:81–89. [PubMed]