The predictive value of childhood blood pressure values for adult elevated blood pressure ()
1. INTRODUCTION
An increase in prevalence of HBP in children may lead to an increase in the prevalence of HBP in adults, assuming that BP in childhood and adolescence tracks with fidelity into adulthood. Rosen et al. reported a prevalence of HBP in 33% of men and 50% of women enrolled in a multi-ethnic study of the development of atherosclerosis in asymptomatic adults [1]. This high prevalence constitutes a public health threat because HBP engenders structural changes in the heart which frequently lead to congestive heart failure.
The variability of blood pressure makes it difficult to identify a hypertensive child at a single visit. Longitudinal modeling is an appealing way to build a predictive model for childhood blood pressure. Such an approach necessitates longitudinal blood pressure data that span the age range from childhood to adulthood so that the model’s predictions are based on empirical longitudinal data. We have previously applied autoregressive models to longitudinal pulmonary function data [2,3]. These models were subsequently enhanced by the inclusion of random effects and were applied to longitudinal blood pressure data from the Fels Longitudinal Study (FLS) with follow-up from 1930 to 1989 [4]. Beckett et al. reported that damped autoregressive models and the autoregressive moving average model (ARMA) (1, 1) both provided a good fit the longitudinal data from 2 to 60 years of age [4]. We applied longitudinal models to analyze serial childhood BP data for their ability to predict HBP at age 35 years in the FLS and reported that a child’s systolic blood pressure (SBP) measured on a single occasion and found to be above the 50th percentile of the CDC/NCHS 1999-2000 dataset predicts odds ratios of developing HBP as an adult ranging from 1.8 to 4.4 depending on sex and age [5].
The usual algorithm is to define children as hypertensive 1) if their systolic and/or diastolic BP is ≥95th percentile for age, gender and height on 3 consecutive visits [6]. Children are designated as pre-hypertensive if their systolic and/or diastolic BP is ≥90th percentile, or 2) if their SBP/DBP is above 120/80 and 3) they are not hypertensive. In this study we extend the analyses from birth to 80 years of age and the follow-up of the FLS cohorts through 2010, and we consider extensions of autoregressive models to account for measurement error and individual level time-to-time dependence. We include concurrent measures of BMI and height as timevarying covariates and include use of antihypertensive medication as an additional covariate in the models. The objectives are to:
1) Describe the patterns of change from childhood through late adulthood in SBP and DBP using longitudinal models;
2) Evaluate the effects of concurrent BMI, height, and antihypertensive medications on BP.
3) Employ fitted models to predict adult HBP or preHBP from childhood values in the same individuals;
4) Determine the risk of adult HBP for individual childhood values ≥ 0th percentile of the 4th Report relative to the risk of adult HBP for individual childhood values at the 50th percentile.
2. SUBJECTS AND METHODS
2.1. Study Sample
The FLS is the world’s largest and oldest longitudinal observational study of human growth and body composition. [7,8] Initiated in 1929 in Yellow Springs OH, the FLS has been housed at Wright State University in Dayton, OH since 1977. Of the total enrollment of 2079 subjects eight percent have been lost to follow-up and 16 percent have died. Data of the lost and deceased participants are used where appropriate. Childhood measurements made from birth through age 7 years include anthropometrics and blood pressure. These data are recorded during scheduled examinations at birth, 1, 3, 6, 9, and 12 months, then semi-annually to 18 years, and biennially thereafter. Reliability in the FLS is excellent with reliability coefficients for most variables well above 90%.
2.2. Measurement Protocols
The body size measurements were taken following recommendations in the Anthropometric Standardization Reference Manual [9]. Weight was measured to 0.1 kg using a SECA scale. Height was measured to 0.1 cm using a Holtain stadiometer. Body Mass Index (BMI = weight/height2; kg/m2) was calculated. All the measurements were taken twice, with a third measurement taken if the difference between the first two exceeds established tolerance (0.3 kg for weight and 0.5 cm for height), and the average values are used for analysis.
SBP and DBP were measured with a mercury sphygmomanometer semiannually from ages two through 18 years and every two years thereafter. SBP and fourth and fifth phase DBP were measured by trained observers with the participant seated. After 1974 a rigorously standardized protocol similar to that used in the Multiple Risk Factors Intervention Trial (MRFIT, 1974) was followed [10]. Arm circumferences were measured, and cuffs were used that were long enough to completely encircle the arm and wide enough to cover two-thirds of the length of the upper arm. Three measurements were taken at a single examination, and the average of the second and third readings was used for data analysis.
The data from participants using medications were included in our model since antihypertensive medications are designed to attenuate HBP. Data generated by participants with HBP who are not taking medication were analyzed with the data for participants who do not have elevated BP. We established an index indicating which participants have HBP but take no medication; have HBP and take appropriate medication; and have no elevated BP. This index was incorporated as a covariate in our analyses.
2.3. Statistical Analysis
The longitudinal analyses used data from multiple examinations over time for each individual. The analyses included periods of time ranging from 2 to 80 years. Individual serial lifetime SBP measurements were modeled as a quadratic function of age and BMI and a linear function of height. We sought to determine how closely individual BP measurements taken at two different ages were correlated.
If BP measurements on different occasions follow a multivariate normal distribution, then the mean and covariance structure are sufficient to describe the process of change in BP over time. In particular, these parameters can be used to address whether elevated BP measurements in childhood presage HBP in adulthood.
The serial changes over time for SBP were modeled using a covariance pattern model, which allows for a relaxing of the independence assumption on the error terms. Independent variables in this model included a subject’s age, BMI, height, and a dichotomous variable indicating whether a subject was on any blood pressure medications. The age (years) of the subject was each centered at age 40 (Age-40) and a quadratic function term was included to adjust for any non-linearity in the relationship between a subject’s age and their SBP measure. A similar strategy was used for BMI (kg/m2), but with BMI centered at 20 (BMI-20). A subject’s concurrent serial height (cm) was centered at 150 (Ht-150), and no quadratic term was included.
In a linear regression model, the error term is assumed to be independent across observations and normally distributed with a mean of zero and constant variance. In the covariance pattern model, the independence assumption is relaxed, and the error terms are allowed to be related within each individual. The vector of the error term still has zero expectation, but is allowed to have a variance-covariance matrix with a pre-specified structure. Three covariance structures were considered:
1) Compound Symmetry [CS]. In this model, all measurements at different times are assumed to have the same correlation. The covariance matrix has elements within each individual:, where 0 ≤ ρ ≤ 1.
2) Autoregressive (1) [AR (1)]. The correlation between measurements at examinations 1 year apart is ρ, 2 years apart is ρ2 and so on. Within individual covariance has the form: σ2∙ρ|i-j|.
3) Autoregressive-moving average (1, 1) [ARMA (1, 1)]. The correlation between measurements taken one year apart and the correlation between examinations greater than one year apart are parameterized separately. In the covariance matrices fit in this analysis the elements have the following form:
A similar but separate model was fit using DBP as the response. The means and covariances were estimated using the MIXED procedure from SAS. The analysis used measurements made at various ages as well as timevarying covariates, such as height and BMI and nontime-varying covariates, such as antihypertensive medication. It was assumed that missing values were missing at random and were estimated by a multiple imputation procedure. The maximum likelihood estimate (MLE) assumes that the mechanism that causes data to be incomplete is ignorable as defined by Rubin [11]. Comparisons between models were carried out by the Akaike Information Criterion (AIC) [12,13].
3. RESULTS
The demographic characteristics of the participants are shown in Table 1" target="_self"> Table 1. Almost one-half of the study subjects were less than 17 years at initial measurement and 40% were ≥40 years old when final measurements were taken and were followed for ≥20 years. Ninety-eight percent were white, 965 were males, and 1114 were females. One hundred and eighty-eight men and 145 women of the 2079 FLS subjects have been prescribed antihypertensive medication.
The mean values for SBP and DBP at each age were plotted from age 2 to 80 years as shown in Figure 1. The best fits were the first order autoregressive moving average criterion value based on the minimum AIC. The ARMA (1, 1) is an extension of autoregressive models that account for measurement error and thus perform better than the remaining two models: AR (1) and con-
Table 1. Demographic characteristics of blood pressure measurements taken in participants in the Fels Longitudinal Study (N = 2079), southwestern Ohio, 1929-2010.
stant correlation (CS) (Table 2).
This implies that the moving average variance structure, rather than the CS or the strict autoregressive (AR 1) model fits the data better. The ARMA (1, 1) covariance and the model parameters are shown in Table 3. In males and females SBP increases with age in a quadratic trend, and a significant quadratic trend exists for the relationship between BMI and SBP. Once an individual is taking antihypertensive medication, SBP decreases significantly. The results are shown for estimating the covariance parameters for the ARMA (1, 1) model, where ρ is the
Figure 1. Mean systolic blood pressure and diastolic blood pressure levels for 965 males and 1144 females in the FLS.
Table 2. Model fit criterion. Most desirable/lowest value in bold.
autoregressive parameter, and γ is the parameter associated with dependence of subsequent observations. All of the covariance parameters from the final model are shown in Table 3.
The ARMA (1, 1) model allows the identification of children at risk for elevated BP as adults. The probability of risk of a condition in adults, given certain values in the same individuals as children, can be estimated as shown in Beckett et al. [4]. The conditional probability of Y2 being larger than a criterion value (e.g., SBP > 130 mmHg) given Y1 of SBP of a certain value can be derived the fitted model, assuming that Y1 and Y2 follow a bivariate normal distribution. For example, Table 4 shows that if a 15-year-old male has a SBP of 115 mm Hg, the likelihood of his SBP exceeding 130 mmHg at age 35 years is 0.28. The corresponding risk for a female is much lower at 0.06.
In boys the probability of elevated BP at age 35 years is 0.46 with a SBP of 130 mmHg at age 15 years. The trends for the probabilities of elevated BP in adulthood for females were similar. For males, the probabilities increase from 0.18 to 0.33 with increased SBP at age 15 years from 105 mmHg to 120 mmHg, corresponding to a change in risk ratio from 2.3 to 4.3. Although the probabilities are much lower in girls, the corresponding increase in risk ratio for females ranges from 2.6 to 5.7.
Application of Findings
In practice, the conditional probability of the risk of adult values exceeding a cutpoint given a childhood value, , is useful in screening programs. From the selected fitted models, the mean and the covariance of repeated measurements of blood pressure are ascertained, and the prediction of adult blood pressure, SBP ≥ 130 mmHg or DBP ≥ 85 mmHg can be made from childhood measurements, assuming that the measurements are normally distributed.
The probabilities of elevated BP at age 35 years with reference percentiles at the 50th, 90th, 95th, and 99th from the Fourth Report on the Diagnosis, Evaluation, and Treatment of High Blood Pressure in Children and Adolescents are shown in Table 5: e.g., if a 9-year-old boy has a SBP at the 90th percentile, the likelihood of his SBP at 35 years exceeding 130 mmHg is 30% [6]. In men, the probabilities of elevated BP at age 35 are highly dependent on their childhood blood pressure. In females, the probability of having elevated BP at age 35 is low regardless of childhood values.
Table 6 presents risks and risk ratios for developing elevated DBP in adulthood given childhood BP measurements at the 50th, 90th, 95th and 99th percentiles. In general the risks of developing elevated DBP are less at all ages and percentiles than for SBP.
Figure 2 is shaded differentially to indicate age ranges
Table 3. Parameter estimates for systolic and diastolic blood pressure from ARMA [1, 1] models in participants in the Fels Longitudinal Study, southwestern Ohio, 1929-2010. Information for the SBP and DBP are in separate sections.
Table 4. Predictive value of systolic and diastolic blood pressure at age 15 for risk of elevated BP (SBP > 130; DBP > 85 mmHg) at age 35 in participants in the Fels Longitudinal Study. Southwestern Ohio, 1929-2010.
Figure 2. Selected percentiles for systolic blood pressure in males and females from the Fourth report on the Diagnosis, Evaluating, Treatment of High Blood Pressure in Children and Adolescents [6]. Segments of 50th, 90th, 95th, and 99th percentile lines are differentially shaded to indicate differences in the probability that systolic blood pressure at 35 y will be >130 mmHg [6].
Table 5. Predictive value of childhood systolic blood pressure percentiles at 50th, 90th, 95th, and 99th for age at the median height for risk of elevated BP at 35 years of age in males and females [6].
Table 6. Predictive value of childhood diastolic blood pressure percentiles at 50th, 90th, 95th, and 99th for age at the median height for risk of elevated BP at 35 years of age in males and females [6].
during which the probability of elevated BP (SBP > 130 mmHg) at 35 y was either <10%, 10% - 19.9%, 20% - 29.9%, 30% - 39.9%, 40% - 49.9%, and ≥50%. For example, it can be seen in Figure 2 that a male aged 12 y with a SBP of 123 mmHg is at the 95th percentile for Fourth Report data [6]. The corresponding SBP value is toward the upper end of the distribution for the general US male population. Figure 3 also show the likelihood of HBP at 35 y depending on the DBP percentile during childhood.
4. DISCUSSION
The increasing prevalence of HBP is a public health threat because HBP engenders structural changes in the heart which may lead to congestive heart failure (CHF). The extent of the problem was documented by Coca et al. who reported that 73% of a population-based sample of 946 adults with essential HBP met the echocardiographic criteria of de Simone et al for LVH [14,15].
An increase in prevalence of elevated BP in children may lead to an increase in the prevalence of HBP in adults, assuming that BP in childhood and adolescence tracks with fidelity into adulthood. To address this tracking question definitively, Chen & Wang performed a meta-analysis of 50 cohort studies of measurements of BP over time and reported the average BP tracking correlation coefficients were 0.38 for SBP and 0.28 for DBP [16]. A more recent meta-analysis of BP tracking was performed in 29 cohort studies by Toschke et al who reported that the overall adjusted average SBP tracking coefficient was 0.37 between ages 10 and 20 years in males and females [17]. Rosner and coworkers reported that elevated BP in childhood is the strongest predictor of HBP in adulthood, outweighing the influence of childhood height and weight [18,19].
Little is known about the ability of childhood blood pressure to predict adult elevated BP because of the paucity of serial data collected on the same individuals over several decades. The FLS data enabled us to link HBP in adults to BP in the same individuals as children. The ability to predict HBP in adults by measuring BP in children is of clinical and public health importance, especially in view of the increasing prevalence of HBP in children in this country [20].
One approach to analyzing long-term serial blood pressure data is to build a longitudinal model for blood pressure which spans the range from childhood to adulthood and to make predictions based on the model. We have previously considered autoregressive models for longitudinal data which were applied to longitudinal pulmonary function data [2,3]. These models were subsequently enhanced with the inclusion of random effects and applied to longitudinal BP data from the FLS with follow-up from 1930 to 1989 [4]. In the current study we extend follow-up of this cohort through 2010; consider extensions of autoregressive models to account for measurement error; and incorporate BMI, height and antihypertensive medications to improve the accuracy of prediction. The ARMA (1, 1) model that takes into account measurement error was selected to predict adult elevated BP from childhood blood pressure values in the same individuals.
The present study shows that elevated BP during childhood is an important risk factor for HBP at 35 y in males and in females. The predictability of adulthood elevated BP from childhood blood pressure values in the
Figure 3. Selected percentiles for diastolic blood pressure in males and females from the Fourth Report on the Diagnosis, Evaluating, Treatment of High Blood Pressure in Children and Adolescents [6]. Segments of 50th, 90th, 95th, and 99th percentile lines are differentially shaded to indicate differences in the probability that systolic blood pressure at 35 y will be >130 mmHg [6].
present study has limitations. The present study analyzed data from a large longitudinal study of non-Hispanic Caucasians in SW Ohio that were collected between 1929 and 2010. However, analysis of this extensive longitudinal data set may elucidate biologic relationships that apply to all races and ethnicities.
There is a potential cohort effect or secular trend that needs to be taken into account when using these results in clinical or public health assessments. A secular trend toward higher blood pressure in children has been reported by Muntner et al. [21].
It is recommended that the findings reported here be applied to non-Hispanic Caucasian children with blood pressures greater than the 90th percentile for age, sex and height. The present findings may also be applied to children of other US ethnic groups, but the probabilities may differ from those for non-Hispanic Caucasian children. It is probable that the accuracy of the prediction in the present study could be improved by incorporating parental data, dietary patterns, smoking behavior, and physical activity.
4.1. Perspectives
The present study takes advantage of the unique FLS data set with a rich 80-year history of serial collection of data on BP and direct measurements of body composition. The data sets are characterized by high reliability. The ability to predict adult elevated BP by measuring levels of blood pressure in childhood as early as two years of age is of clinical and public health importance, especially in males whose average risk at age 15 years for elevated BP at age 35 years is more than four time that of females with comparable childhood blood pressures. Such prediction may signal the need for corrective interventions that could be made during childhood. This approach is preferable to waiting until adulthood when the deleterious clustering of risk factors has been established, and when lifestyle changes are difficult to implement.
4.2. Innovation and Significance
In most longitudinal studies of BP the period of observation usually does not extend from childhood to adulthood. In the present study, the period of serial observations of BP, BMI and height extends from infancy to old age. The longitudinal model that best fits the empirical data includes error functions, covariates, and imputation. As a result, the AR (1,1) model provides the most accurate predictions extant of HBP in adulthood from given values in childhood. This work is significant because the prevalence of childhood obesity has tripled in the general population over the past three decades, and has quadrupled among several ethnic and racial minority groups [22]. The increased prevalence of childhood obesity has been accompanied by an increase in the prevalence of high blood pressure (HBP) in children and adolescents [21,23]. We have found that childhood values of SBP between the 50th and the 95th percentiles predict HBP in adulthood, Heretofore little attention has been paid to children with BP levels in this range, currently considered to be normal.
5. ACKNOWLEDGEMENTS
This study was supported by National Institutes of Health grants: R01HL40619, U01HL101064, R01HL113697, and HD 12252, Bethesda, Maryland.