| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Hypertension. 2006;47:1147.)
© 2006 American Heart Association, Inc.
Original Articles |
From the Division of Population Genetics and Prevention (D.Gu, S.S., D.Ge, S.C., J.H.), Fu Wai Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China; National Human Genome Center at Beijing (D.Gu, S.S., B.Q.), Beijing, China; Georgia Prevention Institute (D.Ge), Department of Pediatrics, Medical College of Georgia, Augusta; and Institute of Biophysics (B.L., R.C.), Chinese Academy of Sciences, Beijing, China.
Correspondence to Dongfeng Gu, Division of Population Genetics and Prevention, Cardiovascular Institute and Fu Wai Hospital, No. 167 Beilishi Rd, Beijing, 100037, Peoples Republic of China. E-mail gudf{at}yahoo.com
| Abstract |
|---|
|
|
|---|
Key Words: hypertension, essential case-control studies genetics
| Introduction |
|---|
|
|
|---|
Although numerous genetic variations have shown association with hypertension, these associations are often not reproducible. For example, studies of the angiotensin-converting enzyme (ACE) insertion/deletion polymorphism, an extensively investigated variant, have a large number of both positive and negative reports.511 These inconsistent findings might be explained in part by the genetic and environmental heterogeneity among different ethnic groups.12,13 On the other hand, the failure to replicate some single-locus results might be because of an underlying genetic architecture in which genegene interactions are the norm rather than the exception.1416 That is, the effects of the variants under study might be masked by the effects of unstudied variant(s) that affect the phenotype, too.17 Therefore, tests for joint effects of multiple candidate variants may provide more information in the search for hypertension susceptibility genes.
With regard to the biological process of blood pressure regulation, we focused on 11 candidate genes associated with: (1) reninangiotensinaldosterone system (RAAS), including ACE, angiotensin II receptor type I (AGTR1), and aldosterone synthase (CYP11B2)15,18; (2) sympathetic nervous system, including
-1 adrenergic receptor 1A (ADRA1A), ß-2 adrenergic receptor (ADRB2), and tyrosine hydroxylase (TH)4,19,20; (3) lipoprotein metabolism, including lipoprotein lipase (LPL)21; (4) intracellular messengers, including G protein ß polypeptide 3 (GNB3) and epithelium nitric oxide synthase (NOS3) 22,23; and (5) sodium and electrolyte balance, including G proteincoupled receptor kinase 4 (GRK4) and protein kinase lysinedeficient 4 (WNK4)24,25 (as shown in Table 1).
|
In the present study, a large-scale evaluation of 33 candidate gene polymorphisms was undertaken in a Han Chinese casecontrol cohort. We tested for single locus by using typical logistic regression, as well as multilocus association using 2 statistical methods: classification and regression trees (CART)26 and multivariate adaptive regression splines (MARS).27 Both single locus and multilocus analyses were applied to test our hypothesis that these candidate genes under study may contribute to the etiology of hypertension independently and/or through complex interactions.
| Methods |
|---|
|
|
|---|
160 mm Hg and/or diastolic blood pressure
100 mm Hg and 490 age- and gender-matched unrelated healthy control subjects. Three blood pressure measurements were obtained according to a standard protocol recommended by the American Heart Association.29 Subjects with a clinical history of secondary hypertension, coronary heart disease, and diabetes were excluded from this study.
Selection of Candidate Genes and Polymorphisms
We selected 11 candidate genes from biochemical pathways that have been implicated in the development and progression of hypertension. We selected 33 single-nucleotide polymorphisms (SNPs) of these genes based on previous evidence of potential functionality, validated allele frequency, and sequence-proven allelic variation. Detailed information of 11 candidate genes and 33 SNPs is shown in Table 1.
DNA was extracted from leukocytes using a standard phenolchloroform method. All of the SNPs were genotyped using standard polymerase chain reaction/restriction fragment length polymorphism or direct sequencing methods.
Statistical Analyses
We sought evidence of association between each of the 33 SNPs, as well as their interactions and hypertension. Genomic control was used to examine the potential impact of population structure (
) in this study.30,31 To avoid the possible correlations among SNPs in each candidate gene, only 1 SNP was randomly selected from those genes with multiple markers. This process was repeated 1000 times, and the mean value of
was calculated. Our data showed that the mean effect was 1.12, indicating that no correction was necessary.
Single Locus Analyses
First, for descriptive purposes, crude allele and genotype frequencies for each SNP were calculated, and HardyWeinberg equilibrium was evaluated by using a goodness-of-fit test in controls. To avoid assumptions regarding modes of inheritance, all of the analyses were performed using additive, dominant, or recessive modes of each SNP. SNPs with a nominal P<0.1 were presented in the initial results. Second, we performed forward-stepwise multivariable logistic regression analyses using a nominal P value cut point of 0.01 to evaluate the independent effect of each SNP on hypertension. Multiple testing was adjusted using the Bonferroni correction.
The pattern of pairwise linkage disequilibrium (LD) between SNPs within each candidate gene was measured by D and r2, using the software GOLD.32 The haplo.score approach, as outlined by Schaid et al,33 was performed to assess the potential effects of haplotypes within each gene. This method models an individuals phenotype as a function of each inferred haplotype, weighted by their estimated probability, to test the global effects of haplotypes, as well as the individual effect.
Multilocus Analyses
The programs of CART and MARS (Salford Systems) were used to test for potential genegene interaction and thereby to identify specific locus combinations of interest for further investigation and replication.
CART26 can evaluate the relative significance of each predictor, and iteratively subdivides data to build a hierarchical classification model for an optimal combination of independent variables. The strength of CART is its ability to detect high-level interactions among the predictor variables. We ran the program with the following parameters: Gini index as a splitting criterion; a maximum tree depth of 4, indicating that 4-way interaction was allowed; and a minimum terminal node size of 50, ensuring that the null expected number of cases per terminal node would be
25.
MARS27 is a generalization of stepwise linear regression that is particularly suited for high-dimensional problems in which many independent variables might be modeled. MARS has the advantage that additive as well as interactive effects can be included in the models. In this study, 1-way (individual effect), 2-way, 3-way, and 4-way interaction models were considered. The maximum number of basis functions for each model was set to 20, 30, 40, and 60, respectively.
In both CART and MARS, the optimal models were selected in a similar 2-stage process. First, an overly large model was to fit the data by adding new nodes (CART) or basis functions (MARS). Second, the overfitted model was pruned back to a more optimal size. We used 10-fold cross-validation to evaluate overall model fit. The CART or MARS model was developed using randomly divided nine tenths of the data and then evaluated on the remaining one tenth of the subjects. To reduce the variability because of the random stratification into 10 strata, the 10-fold cross-validation process was repeated, and the results were averaged.
| Results |
|---|
|
|
|---|
The demographic and clinical characteristics of all of the individuals are shown in Table 2. The case group had significantly higher body mass index, systolic blood pressure, diastolic blood pressure, serum total cholesterol levels, triglyceride levels, and glucose levels than the controls, as well as lower high-density lipoprotein cholesterol levels. There were no significant differences between the cases and controls for smoking and drinking status.
|
Pairwise LD coefficients were displayed in Table I, available in an online supplement at http://hyper.ahajournals.org. Genotype frequencies for cases and controls, as well as odds ratios (OR) associated with additive, dominant, or recessive modes of inheritance for each of the 33 individual polymorphisms, were calculated. Table 3 presents the finding with a nominal univariable P<0.1 for association with hypertension. As shown, 11 polymorphisms were found in this initial analysis. In the following forward-stepwise multivariable logistic regression analyses, the ADRA1A*G2547C polymorphism (OR, 3.63; 95% CI, 1.43 to 9.23; P=0.0067), ADRB2*Q27E polymorphism (OR, 0.45; 95% CI, 0.298 to 0.68; P=0.0001), ADRB2*R16G polymorphism (OR, 1.46; 95% CI, 1.11 to 1.93; P=0.0078), GRK4*A486V polymorphism (OR, 2.05; 95% CI, 1.50 to 2.79; P=5.1x106), and TH* rs2070762 polymorphism (OR, 2.86; 95% CI, 1.99 to 4.10; P=1.17x108) were found to be independent predictors of hypertension. After the conservative Bonferroni correction (simultaneous adjustment for 33 comparisons), the dominant effects of ADRB2*Q27E and TH*rs2070762 and the recessive effect of GRK4*A486V remained significant associations with hypertension, indicating that these 3 SNPs are likely to represent independent effects.
|
Consistent with the single locus results, the haplotype analysis also found that only the ADRB2 gene, TH gene, and GRK4 gene showed significant effects for global haplotypes after Bonferroni correction (adjustment for 11 candidate genes). For instance, the haplotypes with the A486V C allele within the GRK4 gene had higher frequencies in cases than in controls, whereas the haplotypes with the A486V T allele had lower frequencies. Similar results were observed in the ADRB2 and TH genes (See Table II, available online).
The final CART model selected after pruning is shown in the Figure. The number and percentage of cases (case=1) and controls (case=0) are shown for each node. Five terminal nodes were fit. As shown in the Figure, the 3 most important predictors were the polymorphisms selected in the logistic model. The first split was according to the TH*rs2070762 genotype, indicating a dominant effect. Those with the heterozygous and homozygous mutant genotypes (C allele carriers) were further split according to the GRK4*A486V genotype, in a recessive effect. The next 2 split nodes were ADRB2*Q27E and ADRA1A*C2254G, in dominant and recessive effects, respectively. With those in terminal node 1 serving as the reference group, Table 4 presents naive ORs for terminal nodes 2, 3, 4, and 5. These results suggested that TH*rs2070762 exerted the greatest impact on hypertension risk, followed by GRK4*A486V and AGRB2*Q27E. The CART model also indicated some potential interactions among TH, GRK4, ADRB2, and ADRA1A genes. For example, the combination of TH*rs2070762 C allele and GRK4*A486V VV genotype significantly increased the risk of hypertension (Figure, terminal node 5). However, some decreased risk of hypertension was found among the subjects with TH*rs2070762 C allele, GRK4*A486V A allele, and AGRB2*Q27E E allele (Figure, terminal node 2).
|
|
The final 1-way model chosen using the MARS was completely consistent with the logistic model. The individual dominant effects of TH*rs2070762 and ADRB2*Q27E and recessive effect of GRK4*A486V were found to be independently associated with hypertension (Table 4). The MARS model with 2-way interactions included 2 additive interactions of TH*rs2070762(hm) and ADRB2*Q27E(w) and GRK4*A486V(m) and GNB3*A-350G(hm). Here, w denotes the homozygous wild genotype, h denotes heterozygous genotype, and m denotes homozygous mutant genotype. The naive ORs are shown in Table 4. When 3-way interactions were allowed, a different MARS model was selected. This model contained a dominant effect of ADRB2*Q27E(w); three 2-way interactions: GRK4*A486V(m)-GNB3*A-350G(hm), TH*rs2070762(hm)-GNB3*A-350G(hm), and CYP11B2*IC(hm)-GRK4*A486V(m); and a 3-way interaction of ADRB2*R16G(hm)-TH*rs2070762(hm)-GNB3*A-350G(hm) (Table 4). Although 4-way interactions were allowed, none were selected in the final model (data not shown).
We also used logistic regression to fit the models including all of the terms selected in the MARS models, as well as all of the lower-order interactions and individual effect terms. For 2-way interaction models, however, when the individual effect terms (TH*rs2070762, ADRB2*Q27E, GRK4*A486V, and GNB3*A-350G) were included in the model, the interactions of these terms were no more significant, indicating that these interactions depend on the individual effect terms. The similar results were observed in the 3-way interaction model, except for an interaction of CYP11B2 and GRK4, with marginal significance (P=0.09; data not shown).
To investigate whether there are interactions among those polymorphisms displaying no individual effects on hypertension, we further constructed the models excluding those 3 SNPs with individual effects (TH*rs2070762, ADRB2*Q27E, and GRK4*A486V). Only one 2-way interaction, CYP11B2*IC(hm)-AGTR1*A-1138T(w), was found to be significantly associated with hypertension. This interaction effect was supported by the crude data (Table 5). When the subjects were stratified by the CYP11B2*IC polymorphism, there was an effect of the AGTR1*A-1138T polymorphism only among those with the IC C allele. Compared with the AGTR1*A-1138T T allele carriers, the OR was 2.10 (95% CI, 1.26 to 3.51; P=0.002) for those individuals with AA homozygotes. The logistic regression analysis also indicated that this interaction was independent of the other 3 individual predictors (Table 6).
|
|
| Discussion |
|---|
|
|
|---|
To control for false-positive findings, several approaches were considered. First, we included only Northern Han Chinese subjects, who have been proven to have a significantly higher incidence and prevalence of hypertension and higher mean levels of blood pressure.28,34 Our genomic control analysis revealed no evidence of population stratification in this data. Second, all of the candidate genes selected have a substantial priori probability of involvement in hypertension. Finally, we used conservative Bonferroni correction and 10-fold cross-validation in single-locus and multilocus analyses, respectively, to control the false-positive findings potentially because of the multiple testing and overfitted model.
Although further functional studies are necessary to fully understand the underlying biological mechanisms of the observed genetic associations, our findings are biologically plausible. There are numerous studies demonstrating that essential hypertension is accompanied by sympathetic activation.4,19,20 It has been shown that noradrenaline plays an important role in the regulation of blood pressure, and the increased plasma catecholamine levels have been observed in the development of hypertension. As the rate-limiting enzyme in catecholamine biosynthesis, the TH gene is considered as a logic candidate for the etiology of hypertension.4,35 The present study found a SNP (rs2070762) located in intron 13 strongly associated with an increased risk of hypertension.
As an important component of sympathetic system, the ADRB2 gene has been implicated in the pathogenesis of hypertension, on the basis of both studies, suggesting altered ß2-mediated changes of cardiovascular functions36 and molecular genetic studies.37 Results from linkage and association studies support this gene participating in blood pressure regulation and development of hypertension.3739
With regard to GRK4, this gene has been implicated in essential hypertension, because it participates in the desensitization of the D1 receptor, which leads to sodium retention.24 The role of the GRK4 gene in the dopaminergic system in hypertension has recently been thoroughly reviewed.40,41 The current data indicated that GRK4*A486V was associated with hypertension in recessive mode, after Bonferroni correction, which was consistent with previous reports in white subjects.42
We previously reported linkage and association with hypertension and blood pressure on chromosome 8p22 in our hypertensive families.21 Two candidate genes for hypertension in this region, ADRA1A and LPL, were included in this study. No association with hypertension was observed for LPL gene polymorphisms. However, we found that a novel SNP, C2547G, located in the 3'-untranslated region of ADRA1A gene, was significantly associated with hypertension before Bonferroni correction, with the minor G allele more frequently observed in cases (2.3%) than in controls (0.8%). Because of the low frequency of the minor allele, this association should be interpreted with caution.
Our further exploratory multilocus analyses provided some suggestive joint effects of TH-ADRB2, TH-GNB3, GRK4-GNB3, and TH-ADRB2-GNB3. An interesting finding is the interaction of CYP11B2 and GRK4, with marginal significance (P=0.09) in logistic regression. A recent study constructed in Ghanaian subjects, including 13 polymorphisms of 8 genes, reported an interaction of GRK4 and ACE associated with hypertension.43 Among Japanese subjects, the best combination that was predictive of hypertension included GRK4, ACE, and CYP11B2; however, for low-renin hypertension in Japanese subjects, the best genetic model was also reported, which included only GRK4 and CYP11B2.41 Williams et al43 indicated that robustness for blood pressure might be maintained by genes in different but potentially compensatory pathways, vasoconstriction and sodium balance.
To investigate the potential interactions among those polymorphisms displaying no individual effects on hypertension, we further constructed the models excluding those 3 individual predictors (TH*rs2070762, ADRB2*Q27E, and GRK4*A486V). An interactive effect on hypertension was found between CYP11B2*IC and AGTR1*A-1138T, both from RAAS. This joint effect was supported by the crude data (Table 5) and logistic regression analysis (Table 6). The effects of combinations of RAAS gene polymorphisms on blood pressure and hypertension have also been investigated in previous studies, in which the combinations of polymorphisms were associated with the risk of hypertension, although no individual effect of each isolated genotype was detected.15,16 Recent studies indicated that the CYP11B2 gene and protein, locally presented in kidney, were regulated by low salt intake and angiotensin II type 1 receptor.44 These findings suggested an epistatic interaction in these 2 or more candidate genes in RAAS.
This study attempted to apply advanced statistical methods to address 2 common problems encountered in genetic dissection of complex traits. The first is the multiple testing, which is almost inevitable in large-scale gene mapping efforts. One approach used in single-locus analyses was conservative Bonferroni correction to control for the false-positive results. The other approach used in multilocus analyses was cross-validation to prune the overfitted model. The second is the high-dimensional problem caused by gene-gene interaction, which has been widely accepted as an important contributor to the complexity of mapping complex disease genes.17 In the present study, we used CART and MARS to fit the potential interactive models. Both methods offer advantages over traditional logistic regression in that they may discover interactions of genes displaying no strong individual effects45 and have been used in genetic association studies.4648
Perspectives
The present study provides further evidence that several functional polymorphisms within candidate genes act individually or together in the etiology of essential hypertension. Our data also demonstrate the use of nonparametric methods, such as CART and MARS, in detecting a genegene interaction effect for a complex disease. Further replications in larger independent samples are warranted. In addition, functional studies to prove the true existence of an interaction between the CYP11B2 and AGTR1 genes are also required in the future.
| Acknowledgments |
|---|
Received November 22, 2005; first decision December 19, 2005; accepted March 15, 2006.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. Zeng, V. A. M. Villar, G. M. Eisner, S. M. Williams, R. A. Felder, and P. A. Jose G Protein-Coupled Receptor Kinase 4: Role in Blood Pressure Regulation Hypertension, June 1, 2008; 51(6): 1449 - 1455. [Full Text] [PDF] |
||||
![]() |
J. A. Staessen, T. Kuznetsova, H. Zhang, M. Maillard, M. Bochud, S. Hasenkamp, J. Westerkamp, T. Richart, L. Thijs, X. Li, et al. Blood Pressure and Renal Sodium Handling in Relation to Genetic Variation in the DRD1 Promoter and GRK4 Hypertension, June 1, 2008; 51(6): 1643 - 1650. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Wang, I. Armando, L. D. Asico, C. Escano, X. Wang, Q. Lu, R. A. Felder, C. G. Schnackenberg, D. R. Sibley, G. M. Eisner, et al. The elevated blood pressure of human GRK4{gamma} A142V transgenic mice is not associated with increased ROS production Am J Physiol Heart Circ Physiol, May 1, 2007; 292(5): H2083 - H2092. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Ge, H. Zhu, Y. Huang, F. A. Treiber, G. A. Harshfield, H. Snieder, and Y. Dong Multilocus Analyses of Renin-Angiotensin-Aldosterone System Gene Variants on Blood Pressure at Rest and During Behavioral Stress in Young Normotensive Subjects Hypertension, January 1, 2007; 49(1): 107 - 112. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Hypertension Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2006 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |