| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
(Hypertension. 2003;41:1027.)
© 2003 American Heart Association, Inc.
Scientific Contributions |
From the Department of Preventive Medicine and Epidemiology, Loyola University Stritch School of Medicine (X.Z., R.C., A.L., D.K.), Maywood, Ill; the McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine (Y.-P.C.C., D.Y., A.C.), Baltimore, Md; and the Division of Hypertension, University of Michigan School of Medicine (A.W.), Ann Arbor, Mich.
Correspondence to Dr Xiaofeng Zhu, Department of Preventive Medicine and Epidemiology, Loyola University Medical Center, 2160 S First Ave, Maywood, IL 60153. E-mail xzhu1{at}lumc.edu or Dr Aravinda Chakravarti, McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21287. E-mail aravinda@jhmi.edu
| Abstract |
|---|
|
|
|---|
Key Words: hypertension, genetic angiotensin-converting enzyme haplotypes angiotensin renin-angiotensin system case-control studies
| Introduction |
|---|
|
|
|---|
The inconsistencies observed in RAS gene-association studies demonstrate the challenges of dissecting complex multifactorial traits like hypertension by single-locus analysis. Inconsistencies might arise from the inadequate power of small sample sizes, population structure, varying effects of several disease-predisposing variants, gene-environment interactions, or poor study design. Each of these factors can hinder the detection of a modest contribution of an individual locus to a trait such as hypertension. Overall, the role of RAS gene variants in hypertension remains unclear, and a more thorough understanding of the patterns of linkage disequilibrium within RAS candidate genes is necessary. Haplotype-based analysis, which considers all of the variants segregating at these loci, can provide additional precision in studies of complex disease. However, accommodating all genetic and nongenetic information in the analytic process is still a formidable challenge.
| Methods |
|---|
|
|
|---|
|
Statistical Analyses
Descriptive statistics were obtained with SAS software (SAS Institute). Because antihypertensive medications are sometime prescribed to nonhypertensive individuals, hypertension was thus defined as treated individuals with an SBP>120 mm Hg or a DBP >70 mm Hg or persons not under treatment but with both SBP >140 mm Hg and DBP >90 mm Hg. Cases and controls were randomly selected from siblings within a family, with 1 sib sampled from each family. Because hypertensive individuals are, on average, older than normotensive individuals in this study, the oldest unaffected sib in a family was selected as the control. Hardy-Weinberg equilibrium (HWE) at each SNP was assessed by the
2 test with one degree of freedom26 in cases and controls separately. Single-locus tests of association between an SNP and hypertension were performed with a standard contingency
2 test for a case-control design. We estimated haplotype frequencies by the maximum-likelihood method from genotype data through the use of the expectation-maximization (E-M) algorithm under the assumption of HWE.2730 A likelihood ratio (LR) statistic was computed by testing equality of the allele or haplotype frequency for cases and controls treated separately versus combined. The null distribution of the LR was calculated by randomizing case and control status in our samples.31 We also compared individual haplotype frequencies between the cases and controls by
2 statistics from a series of simple 2x2 contingency tables by combining all other haplotypes and assessed the probability value with the permutation test.
Our data are drawn from nuclear families in which many parents have missing genotypes and multiple affected or unaffected siblings are genotyped. By selecting 1 sib per family in our case-control analysis, substantial power is sacrificed in this sampling scheme. Because case-control analysis is subject to population stratification, we used the transmission/disequilibrium test (TDT) for further analysis. The TDT statistic proposed by Clayton32 and used by the computer program TRANSMIT is therefore suitable. This method decomposes the full likelihood of the trio of parent and siblings ascertained by the affected offspring into 2 parts: one that depends on population stratification and one that does not. The test of no linkage or association is based on a partial score function that omits the terms most influenced by hidden population stratification. This approach handles missing parental genotype information by using the genotypes from unaffected sibs to infer the genotypes of missing parental genotypes, and it allows for multiple affected sibs. The significance level, indicated by P values, is assessed by a bootstrap procedure based on 10 000 replicates. For multiple polymorphisms, TRANSMIT also reconstructs haplotypes with an E-M algorithm and tests individual haplotype as well as the overall transmission distortion. Cervino and Hill33 used simulations to show that TRANSMIT is robust against population stratification and is more powerful than other TDT approaches, such as STDT34 and RCTDT.35 To further verify the robustness of TRANSMIT to population stratification in our samples, we tested linkage and association in 100 microsatellite markers with an average intermarker distance of 35 cM extracted from a previous genome scan. The number of significant tests can be considered a guide to the probability that TRANSMIT protects us against population stratification. We also performed TRANSMIT by pooling African Americans and European Americans together. To eliminate the possible inflation of a type I error due to the effect of pooling the 2 different samples, we calculated the P value in this way. We retained the pedigree structure and the individuals affected status and simulated the genotypes according to the actual allele frequencies in African American and European American samples separately. We then merged the 2 samples to calculate the test statistic. The P values for the pooled sample are calculated on the basis of 2000 replications.
Analysis with haplotypes of multiple loci can be more powerful than associations based on a single locus. Recent studies that examined the structure of the human genome have documented that linkage disequilibrium is nonrandom and that most of the genome lies in blocks where 3 to 5 haplotypes account for >80% of the segregating variation.3638 In our recent study of linkage disequilibrium of RAS genes, we identified similar patterns.25 Accordingly, we designed a test of the hypothesis that common variants found in the RAS genes contribute to susceptibility to hypertension and relied on the information provided by these haplotype blocks.
| Results |
|---|
|
|
|---|
|
Case-Control Analysis
We selected cases and controls as follows. When there were hypertensive sibs in a family, we randomly selected 1 sib per family as a case, and this family was then excluded from control selection. We then selected 1 sib as a control per family from the remaining families. Because hypertensives were generally older than normotensives, we selected the oldest qualified sib in a control family as a control. Table 2 presents the characteristics of our selected cases and controls. The distributions of body mass index, SBP, and DBP did not change as a result of selecting another sib as a case from families with >2 hypertensive sibs (data not shown). We then carried out HWE tests for all loci among cases and controls separately. After multiple corrections, only AGTR1 A44221G significantly deviated from HWE (P=0.0075 after correction) among the African American hypertensives. We next compared the genotype frequencies between cases and controls for each locus. This analysis indicated that AGTR1 C43732T and A44221G in African Americans and ACE A-239T and AGT C4072T in European American were associated with hypertension at a nominal level of significance (P<0.05).
|
Single-Locus TDT Analysis
Table 3 summarizes the results of testing for associations for the 25 SNPs among African Americans and European Americans separately and in the combined sample. The Obs# column represents the observed number of transmissions from both homozygous and heterozygous parents to the affected offspring, and the Exp# column represents the number of transmissions that would be expected under the null hypothesis of no association. Under the null hypothesis, the score test statistic follows a
2 distribution with one degree of freedom, but we report a simulated probability value based on 10 000 replications.
|
We observed that allele T at REN C4021T (P=0.0006) and G at A44221G (P=0.0006) for AGTR1 were transmitted to the affected offspring significantly more often than expected (Table 3) in the African Americans. The other loci with nominal transmission distortion were REN C-3212T, AGTR1 C43732T, and ACE A-239T in African Americans (P<0.05). In European Americans, ACE C8342T, A10593G, A11599G, G12292G, and C19329T were observed with nominal transmission distortions (P<0.05). Interestingly, some polymorphisms have different transmission distortion in the African Americans and European Americans. For example, the A allele of ACE A-239T is more often transmitted to affected offspring in the African American families but not among the European Americans. In the combined sample, REN C4021T, AGTR1 A44221G, and ACE C8342T, G12292G, and A15990G showed evidence of transmission distortion (P<0.05).
Haplotype Analysis
In a previous study,25 we determined the haplotype structure of the 4 RAS genes. To summarize those findings, SNPs that demonstrate significant linkage disequilibrium with each other are grouped to construct haplotypes, and these haplotype blocks, or linkage disequilibrium domains, are boxed and illustrated in the Figure. A single functional mutation existing in a candidate gene would fall within 1 of these blocks. Thus, haplotypes in the block could explain most of the variation attributable to functional alleles.
Case-Control Haplotype Analysis
Haplotype frequencies for the selected SNPs in a block were estimated with an E-M algorithm for cases and controls separately. Haplotype frequencies within a block were compared between cases and controls by the
2 square test. The overall haplotype distributions between cases and controls were also compared by the LR test. Significance levels (P values) were obtained by permuting cases and controls. Only the block consisting of C43732T and A44221G in AGTR1 was significantly associated with hypertension (P=0.005).
Haplotype TDT Analysis
We next applied the program TRANSMIT to reconstruct haplotypes and test the association between a haplotype block and hypertension from family data. Table 4 summarizes the results of testing linkage and association between a haplotype block and hypertension. We first report the results for African Americans. In the block consisting of C-4021T and C-3212T in REN, haplotype CC is transmitted to the affected offspring significantly less often than expected (P=0.0004), whereas TC and TT are transmitted more often (P<0.03). A global test also suggests that this block is associated with hypertension (P=0.0025). For the block consisting of C43732G and A44221G in AGTR1, haplotype CG is overtransmitted to the hypertensive offspring. The permutation test shows that this discrepancy is highly significant (P=0.0002), and a global test adjusting for the multiple haplotypes is confirmatory (P=0.0038). For ACE, the haplotype consisting of A11599G, A15990G, and A20060G represents the major haplotypes in this block, after excluding A-239T;25 these 3 SNPs also individually show a significant result in the TDT analysis. Considering the haplotypes consisting of these 3 SNPs, only haplotype AAA was transmitted more often than expected under the null hypothesis (P=0.006); the global test was not significant. When we added A-239T, however, the significance level was improved (P=0.0002), and a significant global test was also observed (P=0.047).
|
In our smaller sample of European Americans, similar TDT analysis was performed with the haplotypes. Significant transmission distortion was found in REN and ACE. In REN, haplotype CCTC was overtransmitted to hypertensives (P=0.01), whereas for ACE, haplotype GGG consisting of A10593G, A11599G, and A12292G is overtransmitted. We also observed that the P value for haplotype TCTC is very small, but the corresponding haplotype frequency is too small to consider this P value statistically significant. Global tests show that both genes are associated with hypertension. We did not find any significant association between haplotype and hypertension status in either AGT or AGTR1.
We then performed TDT analyses with TRANSMIT by using 100 microsatellite markers randomly selected from our previous genome-wide scan.24 Only 4 of the 100 markers were in transmission distortion (P<0.05) in African Americans and none in European Americans. Therefore, the type I error of the TDT analysis with TRANSMIT is reasonably well controlled when analyses are performed in African Americans and European Americans separately.
Finally, we performed the TDT analysis in the combined samples from both samples, but we obtained the probability values by simulating the genotypes separately in both samples and retaining the pedigree structures and phenotypes. The haplotype blocks in African Americans were shorter and consisted of a subset of the SNPs in European Americans; thus, we performed analyses based on the haplotype structure found in African Americans. Again, the results show significant evidence of linkage and association between hypertension and REN and AGTR1 (global P<0.05) and marginal evidence between hypertension and ACE (global P value=0.067).
| Discussion |
|---|
|
|
|---|
On the basis of this reasoning, we performed case-control and family-based TDT analyses between RAS genes and hypertension in African Americans and European Americans. In the single-locus analysis with unrelated individuals as controls, we identified only A44221G of AGTR1 as being associated with hypertension in African Americans, with the G allele increasing the risk of hypertension. The odds ratio associated with G was 2.41 (95% confidence interval, 1.264.63). In the single-locus family-based TDT analysis, A44221G in AGTR1 remained strongly significant, with the G allele being transmitted to hypertensive offspring more often than expected in African Americans. This result is consistent with the analysis in unrelated case-controls, but the controls in TDT are different from the unrelated controls. A44221G also significantly deviated from HWE in cases. A deviation from HWE among affected individuals implies disease heterogeneity and marker-disease linkage disequilibrium.42,43 Therefore, A44221G, located on exon 5 of the gene, could be an important polymorphism associated with hypertension. Although the A44221G SNP does not alter protein structure (proline
proline), it might alter the function of a nearby regulatory element or be in linkage equilibrium with another causative variant that is directly involved in hypertension susceptibility. We had very little power to detect linkage and association for this marker in European Americans. The T allele of C-4021T in REN was also transmitted more often than expected in African Americans but not in European Americans, and this might also have been due to the small sample size of European Americans. The combined European American and African American results demonstrate significant linkage and association between hypertension and REN C-4021T, AGTR1 A44221G, and ACE C8342T, A12292G, and A15990G.
We next performed haplotype analysis based on blocks previously defined by Zhu et al.25 Because of limited haplotype diversity in each haplotype block, we were able to restrict tests to the major haplotypes, which reduced the number of comparisons. We first compared the haplotype frequencies between cases and controls. The permutation test showed that the haplotype frequencies in cases were significantly different from controls only for a block consisting of C43732T and A44221G of AGTR1 in African Americans. However, a strong association between haplotype blocks and hypertension was found when haplotype TDT analysis with this block information was performed. In African Americans, REN, AGTR1, and ACE showed significant evidence of association with hypertension. REN and ACE genes also showed significant evidence in European Americans. No significant evidence for a role of the AGTR1 gene was found in European Americans, perhaps due to the small sample size and the rarity of this haplotype contributing to hypertension. Pooling African Americans and European Americans also resulted in significant associations between RAS genes and hypertension.
In the single-locus TDT analysis, we found different transmission distortions associated with the ACE A-239T polymorphisms in African Americans and European Americans. This inconsistent result can be explained in the haplotype analysis. From the pooled result in Table 4, it is apparent that haplotype AAAA of the ACE gene is present primarily among African Americans because its observed count in the pooled data is almost equal to that in African Americans (53.5 vs 55). Haplotype AAAA is overtransmitted to affected offspring. Therefore, the overtransmission of allele A of A-239T to affected offspring in African Americans is due to the overtransmission of haplotype AAAA. Haplotype AAGA is contributed by both African Americans and European Americans and is undertransmitted to affected offspring. Hence, the undertransmission of allele A of A-239T to affected offspring in European Americans is due to the undertransmission of haplotype AAGA. This observation further supports the contention that the different associations observed in the single-marker analyses in the 2 populations might be attributable to inheritance patterns of different haplotypes on which the marker resides.
The family-based TDT method is apparently more powerful than the unrelated case-control analysis in our data sets. First, family-based TDT uses all affected sibs, but only 1 affected sib was selected in the unrelated case-control analysis. The sample size of a family-based TDT is therefore larger than that of unrelated case-control analysis. Second, family-based TDT uses untransmitted alleles as controls; environmental factors are therefore better matched between cases and controls. In comparison, unrelated case-control analysis might be inadequately matched for factors such as age, body mass index, and other important factors associated with hypertension.
Perspectives
Our findings suggest that using haplotype blocks can reduce inconsistencies observed in single-marker analysis. When we compared the frequency of a single marker in the cases and controls, we did not have sufficient power to detect genetic associations, with the exception of the polymorphism A44221G in AGTR1. However, using haplotype blocks dramatically improved the power to detect a disease susceptibility region. With access to very large numbers of SNPs44,45 and improved understanding of the haplotype structure and linkage disequilibrium patterns of the genome,3638,46 haplotype-based association analysis, ranging from the whole genome to large sets of candidate genes, should greatly aid in the dissection and characterization of the genetic basis of common diseases.
| Acknowledgments |
|---|
Received December 16, 2002; first decision January 16, 2003; accepted March 14, 2003.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. Moore, P. Dicker, J. K. O'Brien, M. Stojanovic, R. M. Conroy, A. Treumann, E. T. O'Brien, D. Fitzgerald, D. Shields, and A. V. Stanton Renin Gene Polymorphisms and Haplotypes, Blood Pressure, and Responses to Renin-Angiotensin System Inhibition Hypertension, August 1, 2007; 50(2): 340 - 347. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Su, L. Lee, X. Li, J. Lv, Y. Hu, S. Zhan, W. Cao, L. Mei, Y.-M. Tang, D. Wang, et al. Association Between Angiotensinogen, Angiotensin II Receptor Genes, and Blood Pressure Response to an Angiotensin-Converting Enzyme Inhibitor Circulation, February 13, 2007; 115(6): 725 - 732. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. E. Dickson and C. D. Sigmund Genetic Basis of Hypertension: Revisiting Angiotensinogen Hypertension, July 1, 2006; 48(1): 14 - 20. [Full Text] [PDF] |
||||
![]() |
B. Baudin Polymorphism in angiotensin II receptor genes and hypertension Exp Physiol, May 1, 2005; 90(3): 277 - 282. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Dong, X. Wang, H. Zhu, F. A. Treiber, and H. Snieder Endothelin-1 Gene and Progression of Blood Pressure and Left Ventricular Mass: Longitudinal Findings in Youth Hypertension, December 1, 2004; 44(6): 884 - 890. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.-J. Wu, F.-T. Chiang, W. J. Chen, P.-H. Liu, K.-L. Hsu, J.-J. Hwang, L.-P. Lai, J.-L. Lin, C.-D. Tseng, and Y.-Z. Tseng Three single-nucleotide polymorphisms of the angiotensinogen gene and susceptibility to hypertension: single locus genotype vs. haplotype analysis Physiol Genomics, April 13, 2004; 17(2): 79 - 86. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Siani, P. Russo, F. Paolo Cappuccio, R. Iacone, A. Venezia, O. Russo, G. Barba, L. Iacoviello, and P. Strazzullo Combination of Renin-Angiotensin System Polymorphisms Is Associated With Altered Renal Sodium Handling and Hypertension Hypertension, March 1, 2004; 43(3): 598 - 602. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Hypertension Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2003 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |