(Hypertension. 2003;41:207.)
© 2003 American Heart Association, Inc.
Scientific Contributions |
From the Centre for Genetic Epidemiology (J.S.C., J.L.H.) and the Department of Physiology (S.B.H.), The University of Melbourne, Parkville, Victoria, Australia.
Correspondence to Prof Stephen Harrap, Department of Physiology, The University of Melbourne, Parkville, Victoria 3010, Australia. E-mail s.harrap{at}unimelb.edu.au
| Abstract |
|---|
|
|
|---|
Key Words: antihypertensive therapy blood pressure genetics human epidemiology
| Introduction |
|---|
|
|
|---|
Among genome scans, some use measured pressure,3 but more commonly, treated individuals are excluded.47 Underlying hypertension is sometimes presumed and hypertensive values substituted for treated pressures.7,8 Observed blood pressure ranking is lost in these approaches, and adjustments do not necessarily reflect possible treatment effects. One study used a nonparametric algorithm to adjust BP while taking ranking into account.9 More recently, we have used a method by which we arbitrarily added 10 mm Hg and 5 mm Hg, respectively, to the measured systolic BP (SBP) and diastolic BP (DBP) pressures for medicated individuals in an attempt to reflect the anticipated effect of treatment.10 In this study, we use the data from the Victorian Family Heart Study (VFHS) to examine the effects on estimates of variance components when different approaches are used to adjust BP in treated subjects.
| Methods |
|---|
|
|
|---|
The 767 nuclear families included in this analysis comprised 2912 individuals (1431 male, 1481 female) with 805 sibling pairs (237 sister-sister, 200 brother-brother, 368 brother-sister). The mean age was 55.2 years (SD, 6.4) for fathers, 52.5 years (SD, 5.8) for mothers, and 24.0 years (SD, 3.7) for offspring. The sample included 150 offspring twin pairs, of whom 66 were monozygotic (MZ) and 84 dizygotic (DZ). A total of 244 subjects (8.4% of total) reported receiving BP-lowering medications. These comprised 16.7% (n=128) of fathers, 14.6% (n=112) of mothers, and 0.3% (n=4) of offspring.
Phenotype Adjustment
Mean SBP and DBP were initially adjusted for age by gender within each generation by using linear regression methods. The residuals from these analyses comprise what we term the original data. For untreated subjects, their original data were used in all variance component analysis. For treated subjects, we used either their original data or other approaches to adjustment, namely: (a) Exclusion: Treated subjects were excluded from analyses by coding their phenotypes as missing values. (b) Fixed substitution: Fixed values equivalent to the 95th percentile of the distribution for the relevant age, gender, and generation were substituted for treated subjects. (c) Random substitution: Randomly generated SBP values from a normal distribution of mean 150 mm Hg (SD, 5), constrained to boundaries between 140 and 160 mm Hg, and DBP values from a normal distribution of mean 95 mm Hg (SD, 2.5), constrained to boundaries of 90 mm Hg and 100 mm Hg, were substituted for treated subjects.7 (d) Fixed addition: Based on the known average treatment effects,11 fixed increments of 10 mm Hg SBP and 5 mm Hg DBP were added to treated pressures. After addition, SBP and DBP were adjusted for age by gender within each generation. (e) Stepped addition: To account for the number of drugs, stepped increments of 8/4 mm Hg, 14/10 mm Hg, and 20/16 mm Hg were added to measured SBP/DBP of treated subjects taking 1, 2, and 3 drug classes, respectively. These stepped increments were chosen to achieve an average increase of 10/5 mm Hg in treated subjects. After addition, SBP and DBP were adjusted for age by gender within each generation.
Statistical Methods
Variance component analyses were conducted under the assumption of a multivariate normal distribution within a family with the use of FISHER software.12,13 Standard errors and 95% CIs were calculated by using large sample normal approximations.14 The total variance was decomposed into additive genetic (
a2), shared environmental (
se2), and individual-specific (
e2) variance components. The details of these methods have been published elsewhere.10 Descriptive statistics of phenotypes were calculated with the use of S.A.G.E. software (S.A.G.E. 4.2, Department of Epidemiology and Biostatistics, Case Western Reserve University).
| Results |
|---|
|
|
|---|
|
|
Variance Component Analyses
For both SBP and DBP, the individual-specific component of variance (
e2) is relatively unchanged across the different methods of adjustment for treatment (Table 2). Compared with the original data, excluding treated subjects, resulted in a decrease in total variance, mostly caused by a decrease in the additive genetic component (
a2), both for SBP (by 12.5 mm Hg2) and to a lesser extent for DBP (by 2.7 mm Hg2). It also resulted in a decrease in the shared environmental component (
se2) for SBP (by 8.6 mm Hg2). The total familial component (
a2+
se2) decreased by 21.1 mm Hg2 for SBP.
|
In contrast, compared with the original data, fixed substitution resulted in an increase in total variance caused by increases in both
a2, by 11.5 mm Hg2 and 2.8 mm Hg,2 and in
se2, by 10.8 mm Hg2 and 2.4 mm Hg2 for SBP and DBP, respectively. Similar comments apply to random substitution, fixed addition, and stepped addition.
As a result of fixed substitution, the average increase of SBP was 17.2 mm Hg (SD, 17.6), but changes in individual SBP ranged from a reduction of 53.5 mm Hg to an increase of 64.5 mm Hg. For DBP, the range was -25.0 mm Hg to +35.5 mm Hg (average, +9.6 mm Hg; SD, 9.7). With random substitution, changes in SBP were from -56.7 mm Hg to +59.2 mm Hg (average, +14.4 mm Hg; SD, 17.3) and for DBP were from -22.5 mm Hg to +39.0 mm Hg (average, +10.9 mm Hg; SD, 10.1).
| Discussion |
|---|
|
|
|---|
If exclusion of treated subjects should be avoided, what should be done? The use of measured pressures from treated subjects is not ideal; they represent a biased distortion in a quantitative analysis because treatment lowers BP and is usually applied to those with the highest values. Indeed, treatable subjects might be expected to provide genetic and environmental clues to high BP, discoveries that depend on unadulterated relations between BP phenotypes and familial factors, including genotypes. A potential solution is to adjust measured pressures in treated subjects so that they better reflect the inherent untreated levels.
One approach is to assume that treated patients have underlying hypertension and substitute values typical of hypertensives for measured pressures. When we applied methods of either fixed or random substitution that others have used, there was an increase in mean BP and increases in both genetic and shared environmental components of variance. This indicates that treatment effects were obscuring the familial contributions to BP. The corollary is that treatment counteracts mechanisms controlled by these same genetic and shared environmental factors.
Although the substitution method increased the variance attributable to familial factors, the individual adjustments to BP were of concern. Because the mean measured pressures in treated subjects fell below 150/95 mm Hg, the average increases in pressures were relatively large (SBP, 17.3 mm Hg; DBP, 10.9 mm Hg) and exceeded the anticipated BP-lowering effects of an average of 1.5 BP-lowering drugs per treated individual.11 This questions the assumption that all treated subjects had underlying hypertension. These classes of drugs are also commonly used for other conditions such as myocardial ischemia, heart failure, fluid retention, and tremor. Furthermore, a proportion of treated patients might have had white-coat hypertension with otherwise normal pressures.15 Second and more importantly, BP ranking was scrambled with substitution, and some changes were extreme (up to 3 SD) and many of the treated pressures were replaced with lower values. This resulted from a compression of wide range of measured pressures in treated subjects (from 91.5 to 209.5 mm Hg for SBP and from 59.5 to 117.0 mm Hg for DBP) to a 95th percentile or a constrained distribution of
150/95 mm Hg. Indeed, for individual BP, the substitution approach can obscure the relations between phenotype and genotype and greatly weaken genetic linkage and linkage disequilibrium analyses that depend on precise individual pressure measures.16
Our addition method makes no assumptions regarding the underlying reasons for treatment and is based on reasonable expectations of treatment effects.11 Furthermore, the method avoids major or misdirected adjustments to BP and retains the ranking of treated subjects. Genetic and shared environmental variance components were maximized and individual-specific components minimized when we applied the addition methods, especially the stepped addition method that took into account number of drug classes an individual received.
These findings also have implications for the power of linkage and linkage disequilibrium studies, including genome scans. The likely polygenic nature of BP determination implies that the contribution of individual QTLs might be small. Therefore, strategies to maximize the genetic signal and the power of linkage and linkage disequilibrium mapping are important. The greater the heritability (the proportion of total variance attributable to the additive genetic component), the smaller the sample size required. A change in heritability from 0.4 to 0.5 might halve the number of sibling pairs required.17
Perspectives
The elimination of subjects treated with BP-lowering drugs from quantitative family analyses excludes an informative subgroup and diminishes sample size. Estimates of genetic and shared environmental effects become unreliable, and statistical power suffers. A more logical approach is one in which appropriate pressure increments are added to treated pressure levels, especially if commensurate with drug treatment. This approach not only retains subjects in the analyses but augments the familial genetic and shared environmental signals without increasing the noise from individual-specific sources of variation.
| Acknowledgments |
|---|
Received August 28, 2002; first decision September 20, 2002; accepted October 22, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
I. Hajjar, D. T. Lackland, L. A. Cupples, and L. A. Lipsitz Association Between Concurrent and Remote Blood Pressure and Disability in Older Adults Hypertension, December 1, 2007; 50(6): 1026 - 1032. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. J. Busst, K. J. Scurrah, J. A. Ellis, and S. B. Harrap Selective Genotyping Reveals Association Between the Epithelial Sodium Channel {gamma}-Subunit and Systolic Blood Pressure Hypertension, October 1, 2007; 50(4): 672 - 678. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Zintzaras, G. Kitsios, D. Kent, N. J. Camp, L. Atwood, P. N. Hopkins, and S. C. Hunt Genome-Wide Scans Meta-Analysis for Pulse Pressure Hypertension, September 1, 2007; 50(3): 557 - 564. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Rao, L. Zhang, J. Wessel, K. Zhang, G. Wen, B. P. Kennedy, B. K. Rana, M. Das, J. L. Rodriguez-Flores, D. W. Smith, et al. Tyrosine Hydroxylase, the Rate-Limiting Enzyme in Catecholamine Biosynthesis: Discovery of Common Human Genetic Variants Governing Transcription, Autonomic Activity, and Blood Pressure In Vivo Circulation, August 28, 2007; 116(9): 993 - 1006. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Dauchet, E. Kesse-Guyot, S. Czernichow, S. Bertrais, C. Estaquio, S. Peneau, A.-C. Vergnaud, S. Chat-Yung, K. Castetbon, V. Deschamps, et al. Dietary patterns and blood pressure change over 5-y follow-up in the SU.VI.MAX cohort Am. J. Clinical Nutrition, June 1, 2007; 85(6): 1650 - 1656. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Saavedra The Challenge of Genetic Studies in Hypertension Circ. Res., May 25, 2007; 100(10): 1389 - 1393. [Full Text] [PDF] |
||||
![]() |
F. Rao, G. Wen, J. R. Gayen, M. Das, S. M. Vaingankar, B. K. Rana, M. Mahata, B. P. Kennedy, R. M. Salem, M. Stridsberg, et al. Catecholamine Release-Inhibitory Peptide Catestatin (Chromogranin A352-372): Naturally Occurring Amino Acid Variant Gly364Ser Causes Profound Changes in Human Autonomic Activity and Alters Risk for Hypertension Circulation, May 1, 2007; 115(17): 2271 - 2281. [Abstract] [Full Text] [PDF] |
||||
![]() |
J.-J. Hottenga, J. B. Whitfield, D. Posthuma, G. Willemsen, E. J.C. de Geus, N. G. Martin, and D. I. Boomsma Genome-Wide Scan for Blood Pressure in Australian and Dutch Subjects Suggests Linkage at 5P, 14Q, and 17P Hypertension, April 1, 2007; 49(4): 832 - 838. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. K. Rana, P. A. Insel, S. H. Payne, K. Abel, E. Beutler, M. G. Ziegler, N. J. Schork, and D. T. O'Connor Population-Based Sample Reveals Gene-Gender Interactions in Blood Pressure in White Americans Hypertension, January 1, 2007; 49(1): 96 - 106. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Iwai, K. Kajimoto, Y. Kokubo, and H. Tomoike Extensive Genetic Analysis of 10 Candidate Genes for Hypertension in Japanese Hypertension, November 1, 2006; 48(5): 901 - 907. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Franceschini, J. W. MacCluer, H. H.H. Goring, S. A. Cole, K. M. Rose, L. Almasy, V. Diego, S. Laston, E. T. Lee, B. V. Howard, et al. A Quantitative Trait Loci-Specific Gene-by-Sex Interaction on Systolic Blood Pressure Among American Indians: The Strong Heart Family Study Hypertension, August 1, 2006; 48(2): 266 - 270. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Bianchi Genetic variations of tubular sodium reabsorption leading to "primary" hypertension: from gene polymorphism to clinical symptoms Am J Physiol Regulatory Integrative Comp Physiol, December 1, 2005; 289(6): R1536 - R1549. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. R Burton, K. J Scurrah, M. D Tobin, and L. J Palmer Covariance components models for longitudinal family data Int. J. Epidemiol., October 1, 2005; 34(5): 1063 - 1077. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Bochud, P. Bovet, R. C. Elston, F. Paccaud, C. Falconnet, M. Maillard, C. Shamlaye, and M. Burnier High Heritability of Ambulatory Blood Pressure in Families of East African Descent Hypertension, March 1, 2005; 45(3): 445 - 450. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Ho, A. Pinto, S. D. Hall, D. A. Flockhart, L. Li, T. C. Skaar, P. Cadman, D. T. O'Connor, U. Wagner, N. S. Fineberg, et al. Association Between the CYP3A5 Genotype and Blood Pressure Hypertension, February 1, 2005; 45(2): 294 - 298. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Kupper, G. Willemsen, H. Riese, D. Posthuma, D. I. Boomsma, and E. J.C. de Geus Heritability of Daytime Ambulatory Blood Pressure in an Extended Twin Design Hypertension, January 1, 2005; 45(1): 80 - 85. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. de Lange, T. D. Spector, and T. Andrew Genome-Wide Scan for Blood Pressure Suggests Linkage to Chromosome 11, and Replication of Loci on 16, 17, and 22 Hypertension, December 1, 2004; 44(6): 872 - 877. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. B. Harrap, J. S. Cui, Z. Y. H. Wong, and J. L. Hopper Familial and Genomic Analyses of Postural Changes in Systolic and Diastolic Blood Pressure Hypertension, March 1, 2004; 43(3): 586 - 591. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. J. Camp, P. N. Hopkins, S. J. Hasstedt, H. Coon, A. Malhotra, R. M. Cawthon, and S. C. Hunt Genome-Wide Multipoint Parametric Linkage Analysis of Pulse Pressure in Large, Extended Utah Pedigrees Hypertension, September 1, 2003; 42(3): 322 - 328. [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. |