Abstract
Background
Pulmonary arterial pressure (PAP) varies considerably in heart failure (HF) despite similar degrees of left ventricular (LV) dysfunction. Bradykinin alters vascular tone and common variations in the kinin B2 receptor (BDKRB2) gene exists. We hypothesized that genetic variation in this receptor would influence PAP in HF.
Methods
131 HF patients (>1yr history systolic HF), without COPD, not currently smoking, BMI < 40, without atrial fibrillation completed the study which included a blood draw for genotyping and neurohormones (ACE, A-II, Bradykinin, ANP, BNP, and catecholamines), an echocardiogram for cardiac function and systolic PAP (PAPsys).
Results
Mean LVEF was 29% ∓ 12%, NYHA class 2 ∓ 1, age 56 ∓ 12 yr, BMI 28 ∓ 5 kg/m2. Forty-six patients (35%) were homozygous for the +9 allele, 58 (44%) were heterozygous (+9/-9) and 27 (21%) were homozygous for the -9 allele of the BDKRB2. PAPsys averaged 42 ∓ 13, 38 ∓ 12, and 35 ∓ 11 mmHg for +9/+9, +9/-9 and -9/-9, respectively (p = 0.03). There was a trend towards gene effect for plasma ACE with the highest values in +9/+9 and lowest in -9/-9 patients (9.5 ∓ 10.7, 7.1 ∓ 8.7, and 5.4 ∓ 6.4 U/L, respectively, p = 0.06). There were no differences in plasma bradykinin or A-II, LVEF, or NYHA across genotypes.
Conclusion
These data suggest the +9/+9 polymorphism of the BDKRB2 receptor influences pulmonary vascular tone in stable HF.
Introduction
Patients with heart failure (HF) often develop pulmonary venous hypertension (PH) with an associated reactive component causing elevations in pulmonary vascular resistance (PVR). Pulmonary hypertension, secondary to HF, is a common result of systolic or diastolic dysfunction which results in a hemodynamic shift to the pulmonary circulation, elevated pulmonary venous pressure, and subsequent elevation of pulmonary artery pressure (PAP) and right ventricular pressure.1,2 In addition to this classical mechanism of elevated PAP in HF, there also appears to be a reactive component linked to a number of mitogenic and vasoactive mediators. The imbalance of mitogenic (endothelin-1, interleukin-1, vascular endothelial growth factor, etc.) and vasoactive (endothelin-1, thromboxane A, serotonin, etc.) mediators result in a predominantly vasoconstrictive atmosphere and can result in structural remodeling of the vascular endothelium and underlying smooth muscle. 3 Interestingly however, the degree of PH in HF is highly variable for a given degree of left ventricular dysfunction and disease severity, suggesting the possibility that genetic variation may influence the susceptibility to PH in HF.
As part of the renin-angiotensin-aldosterone system, bradykinin (BK) plays an important role in the cardiovascular system by influencing blood pressure and cell proliferation. 4 Bradykinin is a potent endogenous vasodilator nonapeptide (formed of nine amino acid residues), released from plasma globulins called kininogens. In humans, the biological action of bradykinin is mediated through the activation of two principle G-protein-coupled kinin receptor subtypes, B1 and B2. 5 The vascular B1 receptor is normally expressed very weakly but is markedly upregulated in the presence of inflammation, cardiovascular disease, and angiotensin converting enzyme (ACE) inhibition. 6 8 The endothelial cell associated bradykinin B2 receptor subtype (BDKRB2) is constitutively expressed in most tissues and is considered a much stronger mediator of vasodilation through increased production and release of nitric oxide at the endothelial level, however, it is most highly expressed in the pulmonary vasculature. 9 11 In humans the BDKRB2 gene has been mapped to chromosome 14q32. The gene is more than 25 kb in size and consists of three exons. The presence of a 9 bp deletion (-9) in the gene encoding the BDKRB2 is associated with expression of higher concentrations of receptor mRNA, suggesting its strong functional relevance. 12
Despite the growing understanding that bradykinin demonstrates strong vasodilatory properties in the systemic circulation, little is known regarding the influence of bradykinin, bradykinin receptor regulation, or bradykinin receptor genotype variants on the regulation of pulmonary vascular tonality. 13 15 Due to the relationships between bradykinin and systemic pressure regulation and the relative paucity of information on the genetic interactions in the pulmonary circulation, the purpose of this study was to test the hypothesis that HF patients homozygous for the +9 polymorphism of the BDKRB2 would exhibit elevated levels of systolic PAP.
Methods
Population characteristics
Patients were recruited prospectively from the Mayo Clinic heart failure service and the Cardiovascular Health Clinic over the period of 2000 to 2004. Inclusion criteria included all of the following: patients with a history of ischemic or dilated cardiomyopathy, stable HF symptoms (>3 months), duration of HF symptoms >1 year, left ventricular ejection fraction (EF) ≤50%, body mass index (BMI) <40 kg/m2, and non-smokers with a smoking history <15 pack-years. Patients were treated with standard optimized medications for HF at the time of the study. All participants gave written informed consent after being provided a description of study requirements. The study protocol was approved by the Mayo Clinic Institutional Review Board; all procedures of this investigation conformed to the principles outlined in the Declaration of Helsinki and followed Institutional and Health Insurance Portability and Accountability Act guidelines.
Echocardiographic evaluation
Resting supine Doppler and 2D echocardiographic measurements were performed at rest in the supine position according to the recommendations of the American Society of Echocardiography. 16 Images were acquired with an echocardiographic ultrasound system (Vivid 7, General Electric) and stored for off-line analysis (EchoPac software, version 4.0, General Electric). Maximal early velocity (E), maximal late velocity (A), early to late velocity ratio (E/A), and deceleration time of the transmitral inflow during diastole were obtained from a 2-dimensional apical window. We measured PAPsys as the sum of the highest velocity derived pressure gradient of the tricuspid regurgitant jet recorded from either the apical or parasternal positions and the estimated pressure of the right atrium evaluated from the inferior vena cava diameter and collapse using the simplified Bernoulli equation. 17 20
Left atrial (LA) dimension, left ventricular (LV) mass, LV internal dimension during systole and diastole (LVIDs and LVIDd, respectively), LV posterior wall thickness and left atrial end-diastolic dimension were assessed by 2D apical 2- and 4-chamber views. Estimates of mass were calculated using the formula of Troy and colleagues. 21 Left ventricular mass index was calculated as the quotient of LV mass and body surface area. Stroke volume (SV), cardiac output (CO), and LV mass were calculated according to the methods of Lewis et al. 22 The modified Simpson's rule was used to calculate EF. 23
Measurement of plasma neurohormones
Blood was drawn into Vacutainer tubes with and without EDTA and chilled until centrifuged at 4 °C at 2500 rpm for 10 minutes. Serum and plasma aliquots (1 mL) were transferred into 12 × 75 mm polystyrene tubes and frozen at -80 °C until assayed. Plasma angiotensin converting enzyme (ACE) activity was determined by spectrophotometric methods using reagents from Sigma (St. Louis, MO). Dopamine, epinephrine, norepinephrine, angiotensin II (A-II), and atrial natriuretic peptide (ANP) were assessed according to standard methods developed in the Mayo Clinic immunochemical core laboratory. Plasma brady-kinin was measured by radioimmunoassay (Phoenix Pharmaceuticals, Belmont, CA). Brain natriuretic peptide was measured by specific immunoradiometric assay using two monoclonal antibodies against human BNP, one recognizing a carboxy-terminal sequence and the other the ring structure of human BNP.
BDKRB2 genotyping
Buffy coat was obtained from whole blood collected in EDTA treated tubes and extracted using the Gentra Systems Puregene® DNA Isolation kit. Polymerase chain reaction (PCR) was enhanced for high G/C content with Q solution (Qiagen, Valencia, CA) and conducted using the following primer sequence; (forward) 5′-GCC CTT GAA AGA TGA GCT G-3’ and (reverse) 5′-AAC TCC CCA CGA CCA CAG-3’. 24 Briefly, 6.25 microliters Qiagen HotStartTaq 2X Master Mix in a final volume of 12.5 microliters including 15 ng of DNA, 1.5 mM magnesium chloride, 1.25 U HotStarTaq DNA polymerase, 200 μM each dNTP, and 5 pmoles of each primer was combined. After an initial 15 minutes at 95 °C to activate the HotStarTaq, the DNA was amplified by 35 cycles of 30 seconds at 94 °C, 30 seconds at 60 °C, and 1 minute at 72 °C. These cycles were followed by 10 minutes at 72 °C for final extension of PCR products. The PCR products were purified by ultrafiltration and sequenced in the Mayo Molecular Core Facility by a dye terminator sequencing reaction. The -9 allele produced a 266 bp PCR product whereas the +9 allele produced a 275 bp product.
Statistical analysis
Statistical analysis and graphic presentation were accomplished using SAS (v8.2, SAS Institute, Cary, NC) to compare differences across genotypes. Univariate linear regression models and Mantel Haenszel chi-squared tests for linear trends were constructed to investigate the associations of patient demographics, plasma neurohormones, and echocardiographic variables with bradykinin genotype considered as a 3-level ordinate variable. Multiple variable linear regression models were also constructed; adjusting covariates included age, gender, body weight, height, and exercise history. Hardy-Weinberg equilibrium was tested by chi square analysis. Statistical significance was set at an alpha level of 0.05 for all analyses. Data are presented as mean ∓ standard deviation (SD) unless otherwise indicated.
Results
Participant characteristics
Participant clinical characteristics and medication use at the time of the study, according to genotype, are reported in Table 1. The percentage of women across groups was not different (homozygous +9 = 37%, heterozygous = 34.5%, homozygous -9 = 33.3%). There were no trends across groups for age, height, weight, or BMI. There also were no trends across groups for smoking or exercise history. There was a significant trend in the percentage of the groups using beta-blocker therapy with 88.9% in the homozygous +9, 79.0% in the heterozygous and 55.6% in the homozygous -9 groups. Further, there was a trend towards significant difference for the percentage of the groups using diuretic therapy with 75.6% in the homozygous +9, 63.2% in the heterozygous and 55.6% in the homozygous -9 groups.
Participant Characteristics and Patient Medications According to Bradykinin B2 Receptor Gene Polymorphisms.
Plasma markers of neurohormonal activation and renin angiotensin-aldosterone system (RAAS) activity
Table 2 shows the plasma neurohormone levels related to the RAAS activity. Because of the rightly skewed nature of these neurohormones, the log transformation of the absolute data was analyzed. Although there was not a statistically significant trend across the groups for angiotensin converting enzyme (ACE) inhibitor use, there was a trend towards reduced plasma levels of ACE for the homozygous -9 group as compared to the heterozygous or homozygous +9 groups (Table 2). There were no other significant differences across the groups for plasma levels of bradykinin, angiotensin II, atrial natriuretic peptide (ANP), brain natriuretic peptide (BNP), dopamine, epinephrine, or norepinephrine.
Differences in Plasma Neurohormonal Markers and Markers of RAAS Activity According to Bradykinin B2 Receptor Gene Polymorphisms.
Echocardiographic evaluation
The results of the functional and structural echo-cardiographic evaluation are shown in Table 3. There were no significant differences across the groups for EF, CO, or SV. There were also no significant trends across the groups for diastolic function measures including E and A wave velocities, E/A ratio, or deceleration time. There was a significant trend in PAPsys across the three genotype groups, being highest in the homozygous +9 group, intermediate in the heterozygous group, and lowest in the homozygous -9 group. The were no clinical or statistically significant trends across groups for LV mass or LV mass index, LV posterior wall thickness, LV dimension during systole and diastole or left atrial dimension.
Differences in Functional and Structural Echocardiographic Measures According to Bradykinin B2 Receptor Gene Polymorphisms.
To assess the independence of the genotypes on PAPsys, a multiple linear regression model was constructed, adjusting for age, gender, height, weight, and exercise history. After adjustment, there was a significant association between the level of PAPsys and the homozygous +9 genotype (p = 0.02) (Fig. 1). After adjusting for the same five covariates, an increased plasma ACE level demonstrated a trend toward significance with the homozygous +9 patients (Fig. 2). Further, after additionally adjusting for ACE inhibitor use, this association became significant (p = 0.04).

Relationship between bradykinin receptor B2 (BDKRB2) gene polymorphism and systolic pulmonary artery pressure (PAPsys) in heart failure (HF) patients. Patients who are homozygous +9 for the BDKRB2 demonstrate significantly higher PAPsys compared to those who are heterozygous or homozygous -9.

Relationship between bradykinin receptor B2 (BDKRB2) gene polymorphism and plasma angiotensin converting enzyme (ACE) level in heart failure (HF) patients. Although patients who are homozygous for the -9 BDKRB2 polymorphism demonstrate slightly lower plasma ACE levels, there was no significant trend across the groups.
Discussion
Primary findings
This study evaluated the influence of the bradykinin B2 receptor genotype (BDKRB2) on pulmonary vascular pressure in patients with stable HF. Our results suggest patients homozygous for the +9 polymorphism demonstrate elevated PAPsys despite no significant differences in functional or structural cardiac parameters or heart failure severity. Further, the relationship between PAPsys and genotype polymorphism (p = 0.03) was strengthened after adjusting for potential confounding variables including age, height, weight, gender, and exercise history (p = 0.016). Although ACE levels were low due to chronic pharmacologic suppression, there remained evidence for a relationship between the BDKRB2 genotypes and plasma angiotensin converting enzyme levels. Further, the relationship between ACE levels and genotype was enhanced after adjusting for ACE inhibitor usage.
Renin angiotensin-aldosterone system and regulation of vascular tone
Patients with primary systolic HF have a reduced cardiac reserve which is associated with neurohumoral activation of the renin-angiotensin-aldosterone system and peripheral vasoconstriction. Through blocking the enzyme responsible for conversion of angiotensin I to angiotensin II, ACE inhibitor therapy limits the generation of angiotensin II and thereby reduces the associated vasoconstriction and sodium and water retention. Therapy with ACE inhibition has known benefits in patients with heart failure including improvements in functional capacity, morbidity, and mortality. 25 28 The administration of ACE inhibitor therapy in these patients causes systemic vasodilation, attributable to the loss of angiotensin II-mediated vasoconstriction. 29 Recently, however, it has been suggested that the vasodilatory properties of bradykinin also contribute to the ACE inhibitor-induced vasodilation in both the radial and coronary arteries.13,30 This attenuation of vascular endothelial dysfunction by ACE inhibition appears to result from bradykinin potentiation. 13
Bradykinin is a potent endothelium-dependent vasodilator that has a brief duration of action with a plasma half-life of approximately 15-30 seconds due to its rapid degradation by ACE.9,31 As such, ACE efficiently breaks down approximately 95% of bradykinin in a single passage through the pulmonary circulation. 32 Exogenous administration of bradykinin induces vasodilation of epicardial coronary and resistance arteries in humans mediated in part by nitric oxide and endothelium-derived hyperpolarizing factor. 30 35 The local and systemic vascular effects of exogenous bradykinin administration can be enhanced by ACE inhibition.32,36,37 In hypertensive patients and sodium-deplete volunteers, systemic infusion of a selective bradykinin B2 receptor antagonist attenuates the hypotensive effects of captopril similar to an angiotensin II type 1 receptor blocker. 38 This suggests that some of the short-term hypotensive effects of ACE inhibition are mediated via augmentation of endogenous bradykinin through the BDKRB2.
Bradykinin and vasculature function
Bradykinin demonstrates strong vasodilatory properties in the systemic circulation. Interestingly, little is known regarding the influence of bradykinin, bradykinin receptor regulation, or bradykinin receptor genotype variants on the regulation of pulmonary vascular tonality. In the systemic circulation, bradykinin plays an integral role in vascular regulation. Hornig and colleagues suggest that bradykinin contributes to flow mediated dilation of the brachial artery. These authors have shown that the independent infusion of the BDKRB2 antagonist icatibant results in a 33% reduction of flow-mediated dilation despite no significant effect on baseline arterial diameter or blood flow velocity. After infusion of the ACE inhibitor quinaprilat, flow-mediated dilation increased significantly by 46%. The coinfusion of icatibant and quinaprilat produced a reduction in flow-mediated dilation similar to icatibant alone suggesting bradykinin is fundamental in regulating vascular tone via the BDKRB2. 13
Antagonism of the BDKRB2 also attenuates the short-term blood pressure lowering effects of ACE inhibition. Squire and colleagues have shown, in healthy men on a normo-sodium diet, that icatibant elicits a significant increase in mean arterial pressure. This would suggest a role for bradykinin in basal control of blood pressure. Further, when icatibant was co-administered with the ACE inhibitor perindoprilat the normal reduction in mean arterial pressure was attenuated despite no effect on neurohormonal responses to perindoprilat. 15 Similarly, vasodilation in response to acetylcholine and bradykinin is blunted in resistance arteries of the forearm in patients with coronary artery disease with or without left ventricular dysfunction. Administration of quinaprilat enhanced the vasodilator response to bradykinin and suggests that ACE inhibition improves vasodilation through maximization of the bioavailability of bradykinin. 14
In the pulmonary circulation, Taraseviciene-Stewart and colleagues examined the influence of a long-acting BDKRB2 agonist on pulmonary artery pressure in rats with severe pulmonary hypertension. These authors demonstrated that the administration of the stable BDKRB2 agonist B9972 reduced pulmonary artery pressure, right ventricular hypertrophy, and decreased number of apoptosis derived plexiform lesions in the obliterated vessels. 11
Despite the understanding of bradykinins role in peripheral vascular tonality, little work has been done with regards to the specific polymorphisms of the BDKRB2 and vascular tonality either in the systemic or pulmonary circulations. In a prospective study, Dhamrait and colleagues demonstrated a significant contribution of the BDKRB2 +9/+9 polymorphism to the coronary artery disease risk associated with systemic hypertension. 39 In a more recent study, Pretorius et al. suggest that participants with the +9/+9 BDKRB2 polymorphism have higher systolic blood pressure and vascular resistance which may contribute to the increased left ventricular mass associated with hypertension. 40 These data suggest a role for the BDKRB2 +9/+9 polymorphism in changes in pulmonary pressures and the development of associated pulmonary hypertension. The results of the present study suggest that patients homozygous for the +9 polymorphism have higher PAPsys despite a higher proportion of these patients being treated with beta-blocker therapy. In this context two potential mechanisms for this difference are apparent. First, due to the genetic disposition of these patients, certain pharmacologic therapies may be less effective and require the addition of beta-blockers to stem the influence of co-morbidities. Pharmacogenetic interactions such as this have been suggested, however; the clarity of these interactions remains elusive and thus further research in this area is clearly warranted. 41 43 Second, the patients homozygous for the -9 polymorphism demonstrated less beta-blocker therapy. This course of pharmacologic disease management in these patients may be the result of less symptomology or greater effectiveness of other pharmacologic strategies. As our understanding of pharmacogenomics widens, these interactions may become clearer. 41 43
Limitations
A potential limitation to this study is the use of echo-derived measures of pulmonary artery pressure. Although this methodology has been well validated by a number of investigators, it does not allow for the evaluation of pulmonary vascular resistance thus limiting the robustness of our analysis.17–19,44 As this study group consisted of a relatively homogeneous patient group, with subjects excluded with morbid obesity, extensive smoking history, heart rhythm abnormalities (e.g. afib), and pacer dependency, follow-up studies would be valuable on a more heterogenous HF cohort.
Summary and conclusion
The results of the present study further the findings of others in a number of ways. First, despite maximization of pharmacologically mediated therapeutic regimens which include ACE inhibitor therapy, patients with heart failure continue to develop pulmonary arterial hypertension via a reactive component, most likely mediated through pulmonary arterial endothelial dysfunction. Second, independent of traditional measures cardiac structure and function or circulating plasma levels of neurohumoral markers indicative of activation of the renin-angiotensin-aldosterone system; patients who are homozygous +9 for the BDKRB2 polymorphism demonstrate significantly elevated levels of pulmonary artery systolic pressure. These results highlight the strength of the genetic contribution to the reactive component of secondary pulmonary hypertension in HF and provide more defined pathways to the development of pulmonary hypertension which may ultimately lead to more personalized treatment strategies in this population.
Disclosure
The authors report no conflicts of interest.
Footnotes
Acknowledgements
The authors of this study would like to Kathy O'Malley and Minelle Hulsebus for help with data acquisition and management. This work was supported in part by National Institutes of Health grants HL71478 and HL07111 and American Heart Association grant 0725715Z.
