(Hypertension. 2000;35:574.)
© 2000 American Heart Association, Inc.
Scientific Contributions |
From the Twin Research and Genetic Epidemiology Unit, St Thomas Hospital, London, UK (H.S., U.P., T.D.S.); Department of Cardiology, St Vincents Hospital, Sydney, Australia (C.S.H., R.P.K.); and Gemini Holdings plc, Cambridge, UK (P.J.K.).
| Abstract |
|---|
|
|
|---|
Key Words: genetics twins augmentation index blood pressure aging hypertrophy
| Introduction |
|---|
|
|
|---|
AI is defined as the ratio of augmented systolic pressure (due to the late systolic peak in the central pressure waveform) to pulse pressure8 and represents a measure of a combination of factors related to large arterial function. It is increased in association with arterial stiffness, caused by the degeneration of the arterial wall as a result of repetitive cyclic stress, which is a major determinant of increasing systolic pressure with age and in patients with hypertension.13 Stiffening of the main central arteries also leads to increased velocity of the pulse wave as it travels down the arteries, which causes wave reflections from peripheral sites to return earlier, increasing the dominance of the late systolic peak in the aortic pressure wave. This elevated central systolic pressure increases left ventricular load and myocardial oxygen demand, serves as a stimulus to hypertrophy, and hinders left ventricular ejection, predisposing to left ventricular failure.14 A second major influence on the AI is the intensity of wave reflection. This may explain the more prominent secondary (augmented) pressure peak in shorter subjects and in women compared with men.15 By increasing arterial stiffness, age also increases the intensity of wave reflection as a result of a decrease in the attenuation of reflected waves by the relatively stiffer arterial tree. The importance of AI is supported by studies that show this measure to be an independent determinant of left ventricular mass,9 16 which is considered a stronger predictor of cardiovascular events than BP.17
To date very little is known about the importance of genetic factors in the development of large arterial function. While both systolic (SBP) and diastolic BP (DBP) have been shown to be heritable to a variable extent,18 19 20 21 22 there is no similar information on arterial function itself. Applanation tonometry provides a means whereby information on arterial function, quantified by the AI, can be widely and accurately collected. A feature of AI, however, is its strong dependence on age,8 23 height,15 16 24 heart rate,15 25 and BP.15 26 It is therefore possible that any genetic influence on AI is, at least partly, due to genes in common with mean arterial pressure (MAP), heart rate, and height.
The aim of this study, therefore, was to quantify the genetic and environmental sources of individual differences in AI and to determine the extent to which they are mediated by or dependent on the effects of age, height, heart rate, and MAP. A classical twin design was used to answer these questions. AI and MAP were based on the central aortic pressure wave derived from the radial waveform as measured by applanation tonometry8 27 in a large sample of unselected female twins across a broad age range. Both univariate and multivariate quantitative genetic modeling techniques were used to analyze the data.
| Methods |
|---|
|
|
|---|
Measures
Subjects were interviewed and studied by trained research
nurses. Height was measured to the nearest 0.5 cm with the use of a
wall-mounted stadiometer. Weight (light clothing only) was measured to
the nearest 0.1 kg with digital scales. Brachial BP was measured with
an automated cuff sphygmomanometer (OMRON HEM713C) with the subject in
the seated position under standardized conditions. The average of 2 BP
readings was recorded.
The radial arterial pressure waveform was recorded with the subject supine with the use of applanation tonometry30 and commercially available acquisition and analysis software (SphygmoCor, PWV Medical). Applanation tonometry records fluctuations in the pressure wave by a high-fidelity transducer incorporated into the tip of a pencil-shaped probe. Radial artery recordings were taken because this artery is ideal for tonometry since the artery is supported by the bone of the radius behind, allowing accurate applanation (flattening) by the probe. From this recording, the central arterial pressure was derived with the use of a generalized transfer function, which has previously been shown to (1) give an accurate estimate of the central arterial pressure waveform and its characteristics,31 (2) be constant across individuals, and (3) be constant across conditions as diverse as exercise-related systolic hypertension and the hypotension associated with cardiac failure.32 33
The radial pressure waveform was calibrated to the measured brachial sphygmomanometric pressure, and MAP was derived by integration. The central DBP and MAP were set as equal to the brachial DBP and MAP in the generalized transfer function since it has been shown that these change little in transit down the circulation.34 35
The AI was based on the central aortic pressure wave derived from the radial waveform, as discussed above. As outlined previously, the central pressure waveform is characterized by an early and a late systolic peak in the majority of individuals. The AI is defined as the difference between the early and late pressure peaks divided by the pulse pressure (Figure 1).8 11 In our study, intraobserver and interobserver reproducibility for AI, expressed as mean intraclass correlations, were 0.82 and 0.84, respectively.
|
Analytical Approach
The aims of our analyses were 2-fold. First, to estimate
the influence of genetic factors on AI and a range of BP-related
variables, we applied univariate model-fitting
techniques. Second, to determine the extent to which the genetic
influence on AI is dependent on age, height, heart rate, and MAP, we
used multivariate modeling.
Univariate Twin Analysis
Details of univariate model fitting to twin data
have been described elsewhere.36 37 38 In short, the
technique is based on the comparison of the variance-covariance
matrices in MZ and DZ twin pairs and allows separation of the observed
phenotypic variance into additive (A) or dominant (D) genetic
components and shared (C) or unique (E) environmental components. The
latter also contains measurement error. Dividing each of these
components by the total variance yields the different standardized
components of variance, for example, heritability
(h2), which can be defined as the ratio of
additive genetic variance to total phenotypic variance. By
incorporating age into the model, the influence of age on the
phenotype can also be quantified.39 Extension of
this univariate (1 variable) model to a model including
multiple variables (multivariate) additionally
allows exploration of whether the origin of the covariance
between height, heart rate, MAP, and AI is genetically and/or
environmentally determined.
Multivariate Twin Analysis
A triangular or Cholesky decomposition36 40
including height, heart rate, MAP, AI, and age was used in the
multivariate model fitting. This model enabled us to
quantify which part of the (genetic or environmental) variance
components was specific to AI and which part was due to the influence
of height, heart rate, MAP, or age.
Figure 2 shows the triangular decomposition of the genetic and environmental factors for the 4 phenotypes (and age) included in the analysis. The number of latent factors equals the number of variables. Only 1 twin of a pair is shown. The first factor contributes to all 4 variables, the second factor influences the subsequent 3 variables, and so on. The last factor only influences AI. The observed phenotypes are shown in squares, and latent factors are shown in circles. Factor loadings of observed variables on the different latent factors are represented by the arrows. Factor loadings for the last factor (A4, C4, E4) reflect influences specific to AI and can be represented by hAI, cAI, and eAI respectively. The heritability of AI corrected for the influence of height, MAP, and age can be calculated as follows: h2cor=h2AI/(h2AI+c2AI+e2AI).
|
The genetic correlation (rg) between 2
traits gives an indication of the amount of overlap between (sets of)
genes influencing those traits. rg is
calculated as the (additive) genetic covariance
(COVA) between 2 traits divided by the square
root of the product of the total genetic variance components of
each of the traits. The genetic correlation between, for example, MAP
and AI therefore equals
rg=COVA(MAP, AI)/
(VAMAP*VAAI).
Shared and unique environmental correlations are calculated in a similar fashion.
Model Fitting
A series of submodels nested within the
full-parameter Cholesky model were fitted to the
multivariate variance-covariance matrices. The
significance of variance components A, C, D, and age was assessed by
testing the deterioration in model fit after each component was dropped
from the full model, leading to a model in which the pattern of
variances and covariances is explained by as few
parameters as possible. Submodels were compared with the
full model by hierarchical
2 tests. The
difference in
2 values between submodel and
full model is itself approximately distributed as
2, with degrees of freedom (df)
equal to the difference in df of submodel and full model.
Model selection was also guided by Akaikes information criterion
(AIC=
2-2df). The model with the
lowest AIC reflects the best balance between goodness of fit and
parsimony.
Statistical Software
Data handling and preliminary analyses were done with
STATA.41 Univariate and
multivariate quantitative genetic modeling were
performed with Mx.42
| Results |
|---|
|
|
|---|
General characteristics of the twins who were included in the final analyses are shown in Table 1. Characteristics were very similar in the 2 zygosity groups. Results after exclusion of all pairs of which at least 1 twin was on either antihypertensive or hypoglycemic medication (Table 1) were virtually identical to those of the entire group. The results presented here are therefore based on the entire group.
|
Univariate Analysis
For all measures, intraclass correlations in MZ twin pairs were
larger than those in DZ twin pairs, indicating substantial genetic
influences on all traits (Table 2). This
pattern was confirmed by the univariate model fitting.
Parameter estimates of the best fitting
univariate models and their 95% CI are shown in Table 3. The heritability of AI was estimated
at 37%, whereas age explained 19% of the variance. In the best model
for AI, C could be dropped without a deterioration in fit
(
2[1]=0.0; P=1.0), ie, the
contribution of shared environment to AI was not significant.
|
|
For all BP variables (including heart rate), a full ACE model was chosen as the best model, on the basis of the lowest AIC. Although this method of model selection led to a consistent picture for all BP traits with ACE models providing the most realistic explanation of the data, for several BP traits, h2, c2, or both were not always strictly significant, ie, the CI included zero (see Table 3). Heritabilities ranged from 13% to 25% in the BP variables, whereas shared environment and age explained 12% to 18% and 8% to 20% of the variance, respectively. No influence of age on heart rate was found. Height showed the expected large genetic component (60%) and a somewhat smaller influence of shared environment (25%).
Multivariate Analysis
Values of phenotypic correlations (n=1538) were as expected: AI
was positively correlated with age (r=0.44) and MAP
(r=0.33) and negatively with heart rate
(r=-0.28). Taller women had a lower AI
(r=-0.21), most likely due to later wave reflection from
more distant peripheral reflecting
sites.15 16 24
Results of the multivariate model fitting, including
height, heart rate, MAP, AI, and age, are shown in Table 4. Both hierarchical
2 tests and AIC values indicated that the ACE
model including age showed the best fit to the data. Neither C nor A
could be dropped from the full ACE model because both the AE
(
2[10]=37.0; P<0.001) and the CE
model (
2[10]=123.2; P<0.001)
showed a significant deterioration in fit. There was no indication of a
dominant genetic influence because the ACE model showed a better fit to
the data than the ADE model. The influence of age could not be dropped
from the ACE model without a significant reduction in fit
(
2[4]=318.2; P<0.001).
|
Parameter estimates of the best-fitting multivariate model are not shown because they were virtually identical to the univariate parameter estimates. Moderate genetic correlations between AI and MAP (rg=0.38), AI, and heart rate (rg=-0.38) and between AI and height (rg=-0.29) were found. Unique environmental correlations for the same pairs of variables were smaller, with values of 0.15, -0.32, and -0.02. These results indicate that the phenotypic correlations between especially AI and MAP and AI and height are predominantly due to shared genes influencing these traits.
Sources of variance in AI are shown in Table 5. The last column lists the amount of variance in AI explained by age, height, heart rate, MAP, and AI specific factors. Each of these sources of variance is additionally decomposed in genetic, shared, and unique environmental factors. In the best-fitting multivariate model, 36.5% of the variance in AI was determined by genetic factors, of which 5.6% was due to genes in common with MAP, 4.6% to genes in common with heart rate, and 3.1% to genes in common with height. Approximately 44% was due to unique environmental variation, of which 2.6% was due to environmental influences in common with MAP and 4.6% in common with heart rate. The contribution of shared environment to AI was negligible. It is clear from Table 5 that most of the total phenotypic variance in AI (60.2%) is explained by genetic and unique environmental factors specific to AI, with height (3.1%), heart rate (9.5%), and MAP (8.4%) contributing only relatively small parts of the variance. The remaining variation in AI (18.8%) could be attributed to the effect of age. The heritability of AI adjusted for age, height, heart rate, and MAP equalled the following: h2cor=0.232/(0.232+0.0+0.37)=0.38.
|
| Discussion |
|---|
|
|
|---|
The heritability estimates for central SBP and pulse pressure were both relatively small but similar to the estimate for peripheral SBP. Interestingly, age had a larger effect on central than on peripheral SBP, which may reflect the larger effect of aging on systolic augmentation.23 It may also be related to the greater amplification of the pulse between central and peripheral arteries in the more compliant arterial system of youth. As the arteries stiffen with age, amplification decreases, and central pressure more closely approximates peripheral pressure.45
Intraclass correlations for the BP-related variables were as high as one would expect for traits showing a considerable amount of intraindividual variation. Intraclass correlations for peripheral SBP and DBP were very similar to those reported in earlier studies, although heritability estimates in our study were somewhat lower.18 19 20 21 22 43 44 One reason for this difference may be that twin studies have little power to discriminate between additive genes or shared environment as sources of familial resemblance unless a multivariate approach is used.22 46 Earlier twin studies on BP were substantially smaller than our study and therefore lacked the power to detect a shared environmental component, which may have led to inflated estimates of heritability. Even with a twin study of our size, we had difficulty in making the distinction between additive genes and shared environment in univariate analyses of the BP traits. Another reason for the lower heritability estimates in our study may be the incorporation of age as a variance component in the model. We decided to incorporate age because, if unaccounted for, it can spuriously introduce a common environmental effect.39 We would have observed higher heritabilities if we had either ignored the effect of age or had performed the model-fitting analysis using age-adjusted BP values.
Several potential sources of bias may have had an influence on our results. Medication use is unlikely to have had an effect on the model-fitting results since these remained unchanged after exclusion of twins using antihypertensive and hypoglycemic medication. The mean values of general characteristics of MZ and DZ twins in this study were very similar. Only weight and the percentage of postmenopausal twins showed some slight differences between zygosities. Neither weight nor menopausal status was significantly associated with AI in a multiple regression after the inclusion of age (data not shown). Therefore, these minor differences between MZ and DZ twins are unlikely to have biased the results for AI.
The reported results are likely to be representative of the general population because basic characteristics of the twins were similar to a population-based sample of women participating in the Chingford cohort study in London.47 Furthermore, mortality of adult twins has been shown to be similar to that of the general population.48 However, since they were limited to women, the results of our study cannot be generalized to men. Hayward and Kelly15 observed that from an early age, AI was higher in women than in men. This is most likely related to more prominent wave reflection due to shorter stature in women, rather than due to gender differences in arterial stiffness, since women have been shown to have a more compliant arterial tree, at least until the menopause.49 Future twin studies aimed at systematically investigating gender differences in genetic influences on AI should include opposite-sex twin pairs in addition to same-sex male twins.50
Our estimate of heritability for AI provides a yardstick against which the contribution to AI of either candidate genes or novel loci can be evaluated. Since AI is thought to be a composite of arterial stiffness (determined by the structural components of the arterial wall) and wave reflection from peripheral sites, any gene influencing these pathways may be considered a candidate for AI. A fruitful strategy to detect novel loci may be to measure, in addition to AI, 1 or more related traits, for example, pulse wave velocity in a large number of unselected sib pairs or DZ twins for use in a genome screen. AI and pulse wave velocity may be influenced by a common set of genes, which would offer prospects for gene finding since it has been shown that multivariate genetic modeling increases the power to locate pleiotropic quantitative trait loci.51
In summary, AI has a significant genetic component, which is largely independent of the influence of BP, heart rate, height, and age. Finding genes for AI could help to unravel pathophysiological mechanisms causing left ventricular hypertrophy and may lead to improvements in prevention, diagnosis, and treatment of at-risk populations.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received May 24, 1999; first decision June 15, 1999; accepted September 15, 1999.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S. Munir, A. Guilcher, T. Kamalesh, B. Clapp, S. Redwood, M. Marber, and P. Chowienczyk Peripheral Augmentation Index Defines the Relationship Between Central and Peripheral Pulse Pressure Hypertension, January 1, 2008; 51(1): 112 - 118. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Kupper, J. Denollet, E. J. C. de Geus, D. I. Boomsma, and G. Willemsen Heritability of Type-D Personality Psychosom Med, September 1, 2007; 69(7): 675 - 681. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. D. O'Leary, R. L. Hughson, J. K. Shoemaker, D. K. Greaves, D. E. Watenpaugh, B. R. Macias, and A. R. Hargens Heterogeneity of responses to orthostatic stress in homozygous twins J Appl Physiol, January 1, 2007; 102(1): 249 - 254. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. F. Mitchell, A. L. DeStefano, M. G. Larson, E. J. Benjamin, M.-H. Chen, R. S. Vasan, J. A. Vita, and D. Levy Heritability and a Genome-Wide Linkage Scan for Arterial Stiffness, Wave Reflection, and Mean Arterial Pressure: The Framingham Heart Study Circulation, July 12, 2005; 112(2): 194 - 199. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. M. McEniery, Yasmin, S. Wallace, K. Maki-Petaja, B. McDonnell, J. E. Sharman, C. Retallick, S. S. Franklin, M. J. Brown, R. C. Lloyd, et al. Increased Stroke Volume and Aortic Stiffness Contribute to Isolated Systolic Hypertension in Young Adults Hypertension, July 1, 2005; 46(1): 221 - 226. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Laurent, P. Boutouyrie, and P. Lacolley Structural and Genetic Bases of Arterial Stiffness Hypertension, June 1, 2005; 45(6): 1050 - 1055. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Hernelahti, E. Levalahti, R. L. Simonen, J. Kaprio, U. M. Kujala, A. L. T. Uusitalo-Koskinen, M. C. Battie, and T. Videman Relative roles of heredity and physical activity in adolescence and adulthood on blood pressure J Appl Physiol, September 1, 2004; 97(3): 1046 - 1052. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Greenfield, K. Samaras, L. V. Campbell, A. B. Jenkins, P. J. Kelly, T. D. Spector, and C. S. Hayward Physical activity reduces genetic susceptibility to increased central systolic pressure augmentation: a study of female twins J. Am. Coll. Cardiol., July 16, 2003; 42(2): 264 - 270. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Oliver and D. J. Webb Noninvasive Assessment of Arterial Stiffness and Risk of Atherosclerotic Events Arterioscler. Thromb. Vasc. Biol., April 1, 2003; 23(4): 554 - 566. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. E. North, J. W. MacCluer, R. B. Devereux, B. V. Howard, T. K. Welty, L. G. Best, E. T. Lee, R. R. Fabsitz, and M. J. Roman Heritability of Carotid Artery Structure and Function: The Strong Heart Family Study Arterioscler. Thromb. Vasc. Biol., October 1, 2002; 22(10): 1698 - 1703. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Gardin, D. Brunner, P. J. Schreiner, X. Xie, C. L. Reid, K. Ruth, D. E. Bild, and S. S. Gidding Demographics and correlates of five-year change in echocardiographic left ventricular mass in young black and white adult men and women: the Coronary Artery Risk Development in Young Adults (CARDIA) Study J. Am. Coll. Cardiol., August 7, 2002; 40(3): 529 - 535. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. de Simone, F. Pasanisi, and F. Contaldo Link of Nonhemodynamic Factors to Hemodynamic Determinants of Left Ventricular Hypertrophy Hypertension, July 1, 2001; 38(1): 13 - 18. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Drayna, A. Manichaikul, M. d. Lange, H. Snieder, and T. Spector Genetic Correlates of Musical Pitch Recognition in Humans Science, March 9, 2001; 291(5510): 1969 - 1972. [Abstract] [Full Text] |
||||
![]() |
F. C. Luft Twins in Cardiovascular Genetic Research Hypertension, February 1, 2001; 37(2): 350 - 356. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||