Abstract

It is well known from epidemiology that values of indices describing physiological state in a given age may influence human morbidity and mortality risks. Studies of connection between aging and life span suggest a possibility that dynamic properties of age trajectories of the physiological indices could also be important contributors to morbidity and mortality risks. In this paper we use data on longitudinal changes in body mass index, diastolic blood pressure, pulse pressure, pulse rate, blood glucose, hematocrit, and serum cholesterol in the Framingham Heart Study participants, to investigate this possibility in depth. We found that some of the variables describing individual dynamics of the age-associated changes in physiological indices influence human longevity and exceptional health more substantially than the variables describing physiological state. These newly identified variables are promising targets for prevention aiming to postpone onsets of common elderly diseases and increase longevity.

1. Introduction

Individual age trajectories of physiological indices result from complicated interplay among genetic and environmental (including behavioral) factors taking place during the aging process and so, they may differ substantially among individuals in cohort. Despite this fact the average age trajectories for the same index follow remarkable regularities. Figure 1 shows the average age trajectories of selected physiological indices evaluated from the data on the original cohort of the Framingham Heart Study (FHS).

One can see from this figure that some indices tend to change monotonically with age: the level of blood glucose (BG) increases almost monotonically; the pulse pressure (PP) increases from age 40 till age 85, then levels off and shows a tendency to decline only at later ages. The age trajectories of other indices are nonmonotonic: they tend to increase first and then decline. Physiological average body mass index (BMI) increases up to about age 70 and then declines, diastolic blood pressure (DBP) increases till age 55–60 and then declines, serum cholesterol (SCH) increases till age 50 in males, and till age 70 in females and then declines, pulse rate (PR) increases till age 55 in males and till age 45 in females and then declines, hematocrit (HC) declines after age 70 in both sexes. With small variations, these general patterns are similar in males and females.

The effects of these indices on mortality risk were studied in [13]. It was found that these effects are gender and age specific. The fact that the age dependence affects the shape of mortality risk function provided important insights into the mechanisms by which aging process affects the decline in stress resistance in individuals [46]. It was also found that dynamic properties of individual age trajectories of physiological indices may differ dramatically from one individual to the next.

Researchers continue the debates about determinants of the aging rate and about possible contribution of this rate to life span and healthy life span [717]. Since the rate of aging literally means the rate of changing with age, it would be reasonable to assume that individual differences in the aging rate are to be manifested in variability of dynamic properties of individual age trajectories of physiological indices. And if individual aging rate affects life span and healthy life span then one can expect that dynamic characteristics of such trajectories will affect morbidity and mortality risks.

A number of studies available in the literature support the view about the importance of using dynamic properties of individual age associated changes in physiological indices as the characteristics of aging process that predict morbidity and mortality risks, in addition to the use of the age-specific baseline measurements [1824].

In this paper we investigate the effects of selected parameters describing the dynamic properties of the age trajectories of seven physiological indices on consequent morbidity and mortality risks in participants of the FHS original cohort.

2. Data and Method

2.1. The Framingham Heart Study (FHS)

The FHS Original cohort was launched in 1948 (Exam 1), with 5,209 respondents (55% females) aged 28–62 years residing in Framingham, Massachusetts, who had not yet developed overt symptoms of cardiovascular disease, and continued to the present with biennial examinations (29 exams to date, data from exams 1–25 were used in this study) that include detailed medical history, physical exams, and laboratory tests.

2.2. Phenotypic Traits

Phenotypic traits collected in the FHS cohorts over 60 years and relevant to our analyses include life span, ages at onset of diseases (with the emphasis on cardiovascular diseases (CVD), cancer, and diabetes mellitus), as well as indices characterizing physiological state. The occurrence of diseases (CVD and cancer) and death has been followed through continuous surveillance of hospital admissions, death registries, clinical exams, and other sources, so that all the respective events are included in the study. We used data on onset of CVD, cancer (calculated from the followup data) and diabetes (defined as the age at the first exam when an individual has a value of BG exceeding 140 mg/dl and/or takes insulin and/or oral hypoglycemic agent) to define the age at onset of “unhealthy life” as the minimum of ages at onset of these three diseases. If an individual did not contract any of these diseases during the observation period than the individual was considered censored at the age of the last followup or death. Individuals who had any of the diseases before the first FHS exam were excluded from the analyses of “unhealthy life.” Data on physiological indices include random blood glucose (BG, exams 1–4, 6, 8–10, 13–23), body mass index (BMI, exams 1–25), diastolic blood pressure (DBP, exams 1–25), hematocrit (HC, exams 4–21), pulse pressure (PP, exams 1–25), pulse rate (PR, exams 1, 4–25), and serum cholesterol (SCH, exams 1–11, 13–15, 20, 22–25).

2.3. Definitions of “Dynamic” Risk Factors

We investigated dynamic properties of individual age trajectories of seven physiological indices mentioned above to select factors (referred to as “dynamic” risk factors) capable of affecting mortality risk and risk of onset of “unhealthy life.” BG was excluded from the list of indices for analyses of onset of “unhealthy life” because in the FHS data the onset of diabetes is specifically defined from the values of BG.

First, we evaluated the effect of the rate of changes in physiological indices at ages 40–60 on mortality risk and risk of onset of “unhealthy life” at ages 60+. For this purpose, we approximated the individual trajectories of those physiological indices that have a nearly linear dynamics (both for females and males) at ages 40–60 (BG, BMI, HC, and PP) by a linear function of the form , where is age and is the value of a physiological index at age . Individuals having less than 5 observations of respective index at ages 40–60 were excluded from the analyses. As a result, we have estimates of three risk factors for each individual and each index: an initial value of an index at age 40 (i.e., , referred to as “ ” in Tables 1 and 2 and the text below), the rate of change in the physiological index at ages 40–60 ( , “ ”), and the mean of absolute values of residuals, that is, deviations of observed values of an index from those approximated by a linear function at ages 40–60 (“ ”). The joint effect of these risk factors on mortality and incidence of “unhealthy life” was estimated (separately for each physiological index) by the Cox proportional hazards model with delayed entry (the left truncation time was defined as the maximum of the age at the first FHS exam and 60). Respectively, individuals with ages at death (onset of “unhealthy life”)/censoring below 60 were excluded from the analyses. Note that although the use of linear functions for describing individual aging-related changes is a rough approximation of monotonic changes, it captures important dynamic risk factor—the average rate of change of individual index at the age interval between 40 and 60 years.

Second, we evaluated the effect of dynamic characteristics of physiological indices with nonmonotonic age trajectories on mortality risk and risk of onset of “unhealthy life.” For this purpose, we approximated the age trajectories of such indices (BMI, DBP, HC, PR, and SCH) by two linear functions. The first one approximates the increase in the trajectory at the initial interval [ , ]: , where is age and is the value of a physiological index at age . The second one approximates the subsequent decline in the trajectory at the interval [ , ] after reaching the maximum value at age : . The intervals [ , ] for the fit were defined empirically for each index and sex as follows: for PR (females); for PR (males) and SCH (males); for BMI (males) and DBP (females and males); for SCH (females); for BMI (females) and HC (females and males). Note that the following restrictions on parameters were used in the estimation procedures: , , and to ensure the appropriate shape of the fit. Individuals having less than 6 observations of respective index at ages [ , ] and those having estimates of , at the boundary of allowable values (i.e., nearly zero) were excluded from the analyses. As a result, we have estimates of six risk factors for each individual and each index: an initial value of an index at age (i.e., , referred to as “ ” in Tables 3 and 4 and the text below), the rate of increase in the physiological index at ages [ , ] ( , “Left Slope”), the maximal value of the index approximated by two linear functions describing increase and decline in respective individual indices ( , “Max Index”), age at reaching the maximal value of the index ( , “Age Max”), the rate of decline in the index at ages [ , ] ( , “Right Slope,” see also Figure 2 for illustration), and the mean of absolute values of residuals, that is, deviations of observed values of an index from those approximated by two linear functions at ages [ , ] (“ ”). The joint effect of these risk factors on mortality and incidence of “unhealthy life” was estimated (separately for each physiological index) by the Cox proportional hazards model with delayed entry (the left truncation time was defined as the maximum of the age at the first FHS exam and ). Respectively, individuals with ages at death (onset of “unhealthy life”)/censoring below were excluded from the analyses. If is different for females and males for some index, then the maximum of the two values was used in the (sex-adjusted) model applied to that index. Note that all these calculations were performed for individual age trajectories of respective indices. As a result, each individual is now characterized by a vector of dynamic parameters.

We also evaluated the empirical (Kaplan-Meier) estimates of survival functions (and probabilities of staying free of the diseases defining the onset of “unhealthy life”) for individuals with different values of the dynamic risk factors based on the indices with nonmonotonic trajectories (separately for females and males). For each physiological index and each dynamic risk factor (“Age Max,” “Max Index,” “ ,” “Left Slope,” “Right Slope,” and “ ”), we calculated the values of the risk factor in all eligible individuals from the sample using the procedure described above. Then we evaluated the medians of such empirical distributions of risk factors, separately for females and males. These median values were used to define the sex-specific strata for estimation of survival curves. We assigned individuals of respective sex into one of two strata depending on whether the value of the risk factor for this individual is below (this stratum is denoted as “lower half” in Figures 36) or above (denoted as “upper half” in Figures 36) the (sex-specific) median value. In case of an odd number of individuals, the individual with the value of the risk factor equal to the median was assigned to the upper stratum. Then we calculated the Kaplan-Meier estimates of survival curves (conditional at the sex- and index-specific ages ) for individuals in these two strata. Note that individuals with ages at death (onset of “unhealthy life”)/censoring below were excluded from the analyses, as described above. Respective graphs are shown in Figures 36. For example, the median value of the right slopes calculated for BMI in females equals −0.103. Hence, individuals from the stratum denoted as “lower half” in the upper left graph of Figure 3 are females with values of the right slope of BMI smaller than −0.103. Individuals belonging to the stratum named “upper half” in the upper left graph of Figure 3 are females with values of the right slope of BMI larger than −0.103. The conditional survival curves for the two strata presented in this figure deal with individuals survived until age 75 years, which is the value of for BMI in females, as described above.

2.4. Statistical Analyses

Statistical analyses and graphic output were performed with SAS/STAT (SAS Institute Inc.) and MATLAB (MathWorks Inc.) software packages. values for the regression parameters in the tables were calculated using the Wald chi-square statistic with respect to a chi-square distribution with one degree of freedom using SAS/STAT PROC PHREG. The log-rank test was used to test the null hypotheses about the equality of the empirical survival curves in the strata. Respective values are shown in Figures 36 (SAS/STAT PROC LIFETEST was used for these purposes).

3. Results

3.1. Effect of Individual Dynamics of Physiological Indices at Ages 40–60 on Mortality Risk and Risk of Onset of “Unhealthy Life” at Ages 60+

As described in Section 2, we evaluated the effect of individual dynamics of physiological indices at ages 40–60 on mortality risk and risk of onset of “unhealthy life” at ages 60+ for those indices that have a nearly linear pattern of change at the age interval 40–60 for both females and males. Table 1 shows the estimates of the joint effect of these risk factors on mortality as evaluated by the Cox proportional hazards model. One can see from this table that the variability around the average linear trajectory (“ ”) and the average rate of change between ages 40 and 60 ( ”) are significant risk factors for mortality for all indices. The significance is highest ( ) for the slopes of HC and PP. The initial value of an index at age 40 (“ ”) is also a highly significant ( ) risk factor for mortality for HC and PP (i.e., higher values of respective index at age 40 correspond to higher risk of death compared to smaller values of this index), being nonsignificant for BG.

The effect of these dynamic characteristics on incidence of “unhealthy life” is similar (see Table 2). However, the variability is significant only for PP. Note that the effect of variable “Sex” on both mortality and risk of onset of “unhealthy life” is significant and that the risk for males is higher than those for females.

3.2. Effect of Dynamic Characteristics of Physiological Indices with Nonmonotonic Age Trajectories on Mortality Risk and Risk of Onset of “Unhealthy Life”

For indices with nonmonotonic age trajectories, we evaluated the maximum value of respective index, age at which this maximum is reached, the intercept, and the left and right slopes of the linear functions approximating the increase and decline of respective indices as described in Section 2. Tables 3 and 4 show the estimates of the joint effect of these risk factors on mortality and incidence of “unhealthy life” as evaluated by the Cox proportional hazards model. One can see from Table 3 that the effect of the rate of decline in the index after reaching the maximum (“Right Slope”) on mortality risk is significant for all indices (the highest significance, , is observed for BMI and DBP). In this case, the faster decline in the index after reaching the maximum corresponds to a significant increase in mortality risk (note that the values of “Right Slope” are negative by definition, see Section 2). On the contrary, the rate of increase in the index before reaching the maximum (“Left Slope”) and the initial value from which the increase has started (“ ”) are not significant risk factors for mortality for any index. The estimated maximal value of the index reached is also a significant risk factor for mortality in case of DBP, PR (both have ), and HC ( ). This means that the larger (sex-adjusted) maximal values of respective indices (reached at “Age Max”) correspond to a significant increase in mortality risk. The age at reaching the maximum is itself a significant risk factor for mortality for these indices. The negative values of respective estimates indicate that the later an individual reaches the maximum, the smaller mortality risk is. The variability in all indices except PR significantly affects mortality risks (the larger variability corresponds to higher mortality risks).

The effect of these dynamic characteristics on risk of onset of “unhealthy life” is less pronounced than that on mortality risks. Table 4 shows that the rate of decline after reaching the maximum (“Right Slope”), age at reaching the maximum (“Age Max”), and the variability become nonsignificant for all indices. The maximal value reached is significant only for DBP ( ) and PR ( ). The rate of increase of BMI before reaching the maximum becomes significant ( ), with the faster rate of increase corresponding to a higher risk of onset of “unhealthy life.” Note again that the effect of variable “Sex” on both mortality and risk of onset of “unhealthy life” is significant and that the risk for males is higher than that for females.

3.3. Effect of Dichotomized Dynamic Characteristics of Physiological Indices with Nonmonotonic Age Trajectories

We also evaluated the Kaplan-Meier estimates of survival functions for individuals with different values of the dynamic risk factors based on the indices with nonmonotonic trajectories dividing the entire sex-specific samples into strata representing individuals with the values of the index in the lower and upper halves of the empirical distribution of respective index (see Section 2).

Figure 3 shows the estimates of survival functions for females and males having the average rate of decline of different physiological indices after reaching the maximum (“Right Slope”) from the lower and upper halves of empirical distributions of this risk factor for respective indices. One can see from this figure that the lower absolute values of the slope (i.e., the lower rates of the postmaximum decline) in individuals from the upper half of the distribution are associated with better survival for all indices except SCH for females (nonsignificant results for PR for both sexes are not shown). The highest significance ( ) is observed for BMI in females and SCH in males.

Figure 4 illustrates similar estimates in case of “variability” of different physiological indices (the mean of absolute values of residuals, that is, deviations of observed values of an index from those approximated by two linear functions at respective age intervals, see Section 2). The higher “variability” of trajectories of BMI, DBP, and SCH at respective age intervals in individuals from the upper half of the distribution result in worse survival for both females and males (nonsignificant results for HC and PR are not shown). The highest significance ( ) is observed for DBP in both females and males.

Later ages at reaching the maximal value of DBP and PR in females from the upper half of the distribution are associated with better survival (Figure 5), however this was not observed for males. The higher estimated maximal values of these indices in individuals from the upper half of the distribution correspond to worse survival for both females and males. All other indices did not produce any significant results and are not shown in Figure 5.

Similar calculations for probabilities of staying free of the diseases defining the onset of “unhealthy life” revealed a more mosaic picture. The most consistent results were observed for DBP (Figure 6).

The higher initial values of DBP at age 45 and the higher estimated values of DBP reached in individuals from the upper halves of respective distributions are associated with worse chances of staying free of the “unhealthy life,” for both sexes. The lower rates of the postmaximum decline of DBP in females, but not males, from the upper half of the distribution correspond to better chances of staying free of the “unhealthy life” (Figure 6). In addition, earlier ages at reaching the maximum of SCH for females ( ), the smaller estimated maximal values of BMI ( ) and HC ( ) for females, and PR ( ) and SCH ( ) for males, and a smaller initial value of BMI at age 55 ( ) for females and a smaller initial value of SCH at age 40 ( ) for males, were related to better chances of staying free of the “unhealthy life.”

3.4. Sensitivity Analyses

We should note that the question about the effect of the quality of estimates is important given that at most 11 observations for the monotone indices or 15 observations for non-monotone indices were used (note that for non-monotone indices data from 30-year intervals [ , ], where , were used for calculating dynamic risk factors). These observations were used to estimate two parameters (those of the linear regression) for monotone indices and four parameters (age at reaching the maximal value of the index, intercept, and two slopes) for non-monotone indices. To partly reduce the effect of a poor fit due to a small number of longitudinal observations, we removed those individuals having less than 5 (less than 6 for non-monotone indices) observations from analyses. Clearly, the results could change had we used different numbers for the minimal allowable numbers of observations. To test how sensitive our results are to such changes, we performed sensitivity analyses with different minimal allowable numbers of observations: 4, 6, and 7 for monotone indices, and 5, 7, and 8 for non-monotone ones. The results showed that the effect of dynamic risk factors calculated from individual trajectories of monotone indices at ages 40–60 on mortality risk and risk of onset of “unhealthy life” at ages 60+ is stable across different studies. All risk factors for which the estimates of the regression parameters were significant ( ) in the original study, exhibited similar significant effects on mortality risk or risk of onset of “unhealthy life” in the other studies. Despite some variability in the values of the estimates across the studies, the “direction” of the effect (i.e., the sign of the estimate) was the same in all such cases. In the sensitivity analysis with the largest cutoff value, the values were somewhat larger in some cases which may be explained by a smaller sample size compared to the original study. The effects of dynamic risk factors calculated from individual trajectories of physiological indices with nonmonotonic patterns on mortality risk and risk of onset of “unhealthy life” showed similar stability across the studies as the monotone indices. In all cases (except for the right slope of PR which appeared either marginally significant or nonsignificant in the other studies), risk factors for which the estimates of the regression parameters were significant ( ) in the original study exhibited similar significant effects on mortality risk or risk of onset of “unhealthy life” in other studies. The same observation regarding the sensitivity analysis based on the largest cutoff value was true for non-monotone risk factors too.

4. Discussion

An increase in mortality rate with age is traditionally associated with progressing aging. This influence is mediated by the aging-associated changes in thousands of biological and physiological variables, some of which have been measured in aging studies. The fact that the age trajectories of some of such variables differ among individuals with short and long life spans and healthy life spans indicates that dynamic properties of respective indices affect the life history traits. Our analyses of the FHS data clearly demonstrate that the values of physiological indices at age 40 are significant contributors to both life span and healthy life span (as show the estimates of in Tables 1 and 2), suggesting that normalizing these variables around the age 40 is important for preventing age-associated morbidity and mortality later in life. Two dynamic parameters, and , also have significant effect on mortality risk (the former being more important predictor in case of healthy life span). These data suggest that keeping physiological indices stable over the years of life could be as important as their normalizing around the age 40. Thus, a slower change in an index with age is likely to indicate the slower aging and the lower morbidity and mortality risks.

Table 3 shows that dynamic properties of the indices that change nonmonotonically with age significantly contribute to mortality risks and further demonstrates the importance of maintaining stability of physiological state in aging humans: the lower rate of decline in an index after reaching the age at maximum means the more beneficial effect on all-cause mortality.

The fact that the effect of the studied dynamic characteristics on risks of “unhealthy life” onset (Table 4) is less pronounced than that on all-cause mortality risk may indicate that the dynamic characteristics reflect basic aging-related processes in body that result in increasing nonspecific vulnerability to death with age rather than in increasing vulnerability to a particular pathology.

The review of the literature (below) supports our findings with respect to importance of taking into account longitudinal changes in physiological indices when evaluating/predicting morbidity and mortality risks. One should note, however, that the impact and comparative contributions of dynamic parameters (left and right slopes, variability, intercept) on mortality risks were evaluated in our study for the first time. In our two recent publications we demonstrated that individuals who have different rates of aging related changes in BG levels also differ in longevity [2, 3].

The effects of aging associated changes in serum cholesterol on coronary and all-cause mortality were evaluated in Finnish Cohorts of the Seven Countries study [22]. Men with greatest declines in the cholesterol levels had increased cardiovascular and all-cause mortality compared with men with least change in the levels. In the Paris Prospective Study [23], it was shown that not only a low baseline total cholesterol level but also its decline over time was associated with a higher cancer mortality.

Similar to SCH, high blood pressure (BP) is a major risk factor for CVD. A study of two independent French male cohorts suggested that longitudinal changes in systolic and diastolic BP may be more accurate determinants of cardiovascular risks than baseline BP measures. In both cohorts, the group with a long-term increase in systolic and a decrease in diastolic BP (i.e., with increase in pulse pressure) had the highest relative risk of mortality from CVD compared to the group with no changes in either systolic or diastolic BP, independently of absolute values of BP or other risk factors [25]. Since this study included only males, it is important to note that changes in pulse pressure may in principle have different effects on mortality risks in males and females [26, 27].

The heart rate (HR) is one more index characterizing functioning of cardiovascular system. Prognostic importance of its baseline values as well as variability during 24-hour HR monitoring in patients with heart disease and in general population is recognized [2830]. Contrarily, the prognostic role of the long-term and age-related dynamics of HR is not sufficiently investigated and respective studies are limited. A recent study of the effects of HR at baseline, final HR, and HR change during followup, on survival of patients attending the Glasgow Blood Pressure Clinic revealed that the highest risk of all-cause mortality was in patients who had increased their HR by 5 bpm at the end of followup, as compared with those who had a consistently high (high-high) or low (low-low) HR. Authors concluded that change in HR during the followup is a better predictor of mortality risk in hypertensive patients than baseline or final HR [31].

The body mass index (BMI) is, probably, the most intensively studied index in connection with health and survival. Over recent decades, many studies addressed the effect of BMI dynamics on morbidity and mortality, especially the effect of losing body weight in overweight/obese individuals on risk factors for CVD and diabetes. It was shown that overweight adults who lost weight over 9 years had more favorable (lower) total and LDL cholesterol levels compared to normal-weight control, but less favorable BG levels [32]. In other studies weight loss was associated with excess mortality when compared with weight stability, even when controlled for confounding due to diseases known to cause both weight loss and increased mortality [33, 34].

It was also shown that the weight stability was associated with a lower mortality risk as compared with weight change (gain or loss) [35, 36]. Nilsson et al. [37], showed that in men with decreasing BMI during 16 years of followup the noncancer mortality was higher compared to BMI-stable men. Authors hypothesized that involuntary weight loss in otherwise healthy people could be a sign of premature aging, which in turn caused a nonspecific increase in mortality risk. In other studies, baseline weight and weight change had independent effects on total mortality, with both the associations exhibiting a U-shaped relation [38, 39].

Note that the seven physiological indices used in this paper do not exhaust the list of all possible physiological risk factors for mortality and morbidity. Therefore, the dynamic characteristics calculated from these seven indices cannot explain the entire variability in human life span and healthy life span. Other indices and risk factors can be explored on their association with mortality/morbidity risk if measurements of such indices are available in a longitudinal study for a substantially long-time period. See for example [40] where midlife risk factors were investigated for a cohort of Japanese American men with 40 years follow up.

5. Conclusion

In sum, our results indicate that the dynamic characteristics of age-related changes in physiological variables are important predictors of morbidity and mortality risks in the aging individuals. Previously published epidemiological findings are generally in concert with our results, which clearly indicates the need for further detailed studies of the dynamic parameters of aging related changes in human body with further application of these principles to the prevention strategies. We showed that the rate of changes in physiological state at the age interval between 40 and 60 years may serve as a good predictor of morbidity and mortality risks later in life. For nonmonotonically changing indices, the rates of decline after reaching the maximum, the maximal values, and the age at the maximum are important predictors of morbidity and mortality risks.

Senescence is likely to be the key player in physiological and biological changes observed in aging humans. The dynamic properties of these changes contain important information about the individual aging processes. This information, however, can be masked by the effects of compensatory adaptation and remodeling developing in response to the primary aging process. Studying mechanisms of such adaptation and its connection to morbidity and mortality risks is important for better understanding factors and mechanisms affecting long and healthy life.

Acknowledgments

The research reported in this paper was supported by the National Institute on Aging Grants R01AG027019 and R01AG028259. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute on Aging or the National Institutes of Health. The Framingham Heart Study (FHS) is conducted and supported by the National Heart, Lung and Blood Institute (NHLBI) in collaboration with the FHS Investigators. This paper was prepared using a limited access dataset obtained from the NHLBI and does not necessarily reflect the opinions or views of the FHS or the NHLBI.