Abstract
Background
Inter-relationships between multimorbidity and geriatric syndromes are poorly understood. This study assesses heterogeneity in joint trajectories of somatic disease, functional status, cognitive performance, and depressive symptomatology.
Methods
We analyzed 16 years of longitudinal data from the Health and Retirement Study (HRS, 1998-2016) for n = 11,565 older adults (≥65 years) in the United States. Group-based mixture modeling identified latent clusters of older adults following similar joint trajectories across domains.
Results
We identified four distinct multidimensional trajectory groups: (1) Minimal Impairment with Low Multimorbidity (32.7% of the sample; mean = 0.60 conditions at age 65, 2.1 conditions at age 90) had limited deterioration; (2) Minimal Impairment with High Multimorbidity (32.9%; mean = 2.3 conditions at age 65, 4.0 at age 90) had minimal deterioration; (3) Multidomain Impairment with Intermediate Multimorbidity (19.9%; mean = 1.3 conditions at age 65, 2.7 at age 90) had moderate depressive symptomatology and functional impariments with worsening cognitive performance; (4) Multidomain Impairment with High Multimorbidity (14.1%; mean = 3.3 conditions at age 65; 4.7 at age 90) had substantial functional limitation and high depressive symptomatology with worsening cognitive performance. Black and Hispanic race/ethnicity, lower wealth, lower education, male sex, and smoking history were significantly associated with membership in the two Multidomain Impairment classes.
Conclusions
There is substantial heterogeneity in combined trajectories of interrelated health domains in late life. Membership in the two most impaired classes was more likely for minoritized older adults.
Introduction
For older adults, having two or more chronic diseases (multimorbidity) is associated with poor health outcomes and high costs to the health care system which exceed the risks attributable to any individual constituent disease.1,2 Various studies document these relationships whether conditions are clinically-concordant (e.g., hypertension, heart disease), or discordant (e.g., arthritis, cancer).3,4 The development of multimorbidity is associated with additional clinical phenotypes marked by syndromes that involve impairment across multiple body systems. 5 Still, inter-relationships between changes in multimorbidity and important geriatric syndromes—such as disability, frailty, and cognitive impairment—are not well-characterized. It is critical to better specify these relationships to understand aging processes, incorporate what matters to older adults,6,7 and engineer programs that preserve independence.7,8
The medical complexity and workload involved in managing various health maintenance tasks and coordinating care between multiple providers can be substantial and overwhelming for patients,9–11 especially during transitions to serious illness. 12 This complexity has particularly significant effects among older adults with limited financial means, low levels of education, and inadequate access to timely health care services.9,13 Understanding the unique trajectories of aging that encompass multiple, person-centered health domains is therefore crucial to developing targeted interventions and guiding investments in programs that promote healthy aging.
Research has begun to focus on the synergistic relationships with other facets of wellbeing common in geriatric practice, such as mobility, mentation, and multicomplexity (elements of the “5 Ms Framework”).7,14,15 Older adults’ ability to manage chronic diseases, maintain physical and mental functioning, and address social-emotional needs is critical to attaining high quality of life. Nevertheless, comparatively little work has focused on longitudinal changes across multiple health domains simultaneously to examine aging processes from a comprehensive, whole-person perspective. 16 Long-term assessment may clarify the relationship between domains important to sustaining independence. In addition, variability across multiple health domains may be strongly related to divergent experiences and differential access to socioeconomic resources over the life course. Such cumulative inequality is evident in findings that minoritized older adults are living extended periods of their lifespan in sicker, more disabled states.17–20 To develop and deploy interventions to remediate these observed disparities, it is important to assess multiple intersecting domains of health for vulnerable groups.
This study aims to identify distinct multidimensional trajectories of aging across four health domains that together, convey a holistic view of wellbeing in late life: multimorbidity burden, functional status, cognitive performance, and depressive symptomatology. Because we anticipate that these processes are heterogeneous for population groups due to cumulative inequality occurring across the lifespan, we hypothesize that minoritized older adults and those with low socioeconomic status are more likely to follow aging trajectories characterized by substantial impairment across these health domains.
Methods
Data Source
The Health and Retirement Study (HRS) is a nationally-representative prospective cohort study of community-dwelling adults age 51 years and older in the United States. 21 The HRS is a biennial longitudinal study of over 43,000 people 51 years of age and older that began in 1992 and designed to follow middle-aged and older Americans in pre- and post-retirement age until death. Therefore the HRS provides rich population-based data that capture economic and health transitions of adults in mid and late life (https://hrs.isr.umich.edu/about). Because of cross-wave inconsistency in the assessment of key variables prior to 1998, we utilized all available waves from 1998 through 2016. This study was approved by our university’s Institutional Review Board (IRB ID# STUDY00017034, STUDY00019414).
Study Sample
All community-dwelling respondents aged 65-95 years who participated in the HRS survey during the study period (1998-2016) were eligible for inclusion. Similar to other studies, we limited our sample to respondents who provided a minimum of three observations to allow for stable estimation of non-linear trajectories across multiple outcomes.22,23 A total of 22,272 HRS respondents interviewed between 1998 and 2016 were age-eligible and living in the community. Of those, 15,030 respondents provided the required minimum of three waves of outcome data. We then excluded 382 respondents with proxy interviews at all interview waves and 2,712 respondents with clinically-inconsistent chronic disease self-reports that we were unable to adjudicate using a previously developed multistep adjudication method. 24 Additionally, we excluded 230 non-Hispanic respondents who had missing race/ethnicity data or reported “other” as their racial/ethnic category because of heterogeneity within this racial/ethnic identity group, and 141 respondents with missing data on other baseline sociodemographic characteristics. Finally, we dropped proxy responses at the wave-level because measures of depressive symptoms were not assessed during proxy interviews. The final sample consisted of 11,565 respondents who contributed a total of 66,952 repeated observations between 1998 and 2016, with an average of 5.79 (SD=2.20) observations per respondent.
Outcome Measures
Multimorbidity was operationalized as a count of self-reported somatic diseases (range: 0–7) available and assessed at each wave: heart disease (including myocardial infarction, coronary heart disease, angina, congestive heart failure, other heart problems), hypertension, stroke (excluding transient ischemic attack [TIA]), diabetes, arthritis, lung disease (including chronic bronchitis, emphysema, excluding asthma), and cancer (including any malignant tumors excluding skin cancer). They were assessed at baseline with, “Has a doctor ever told you that you have…”, and at follow-up waves with, “Since we last talked with you, has a doctor told you that you have…”. Exclusion decisions (i.e, TIA, asthma, skin cancer) were made by the HRS.
Functional status was measured at each survey wave assessing self-reported limitations in performing six activities of daily living (ADLs) (walking across a room, dressing, bathing, eating, getting in or out of bed, and using the toilet), and five instrumental activities of daily living (IADLs) (using a telephone, managing money, taking medications, shopping for groceries, and preparing hot meals). Each item was coded as 0 if the respondent reported no difficulty and 1 if the reported having any difficulty in performing the activity. Items were summed to yield a combined score (range 0-11), with higher scores indicating greater functional limitation. 25
Cognitive performance was measured at each survey wave using the 27-point HRS cognitive scale26–28 consisting of: 1) immediate and delayed 10-noun free recall to assess memory; 2) serial sevens subtraction to assess working memory; and 3) counting backwards to assess speed of mental processing. Lower scores represent worse cognitive performance, with a score of ≤11 indicative of cognitive impairment. 28
Depressive symptoms were measured at each wave using the 8-item Centers for Epidemiologic Research Depression (CES-D) scale,29–31 which consists of eight dichotomous items assessing the presence of specific symptoms for much of the prior week. The total number of endorsed symptoms are summed to scores ranging from 0-8. Higher scores indicate greater depressive symptoms. Previous research has shown that a score of four or more is a valid indicator of clinically-significant depressive symptomatology. 32
Predictors
We compared distributions of sociodemographic and health-related factors between identified multidimensional trajectory groups for: mutually-exclusive race/ethnicity categories for non-Hispanic White (White), non-Hispanic Black (Black), and Hispanic; female sex; educational attainment (less than high school education, high school graduate, some college or greater than college education); nativity (foreign-born); household wealth quartile (derived from baseline net worth in US dollars); baseline coupled status (married/living with a partner); baseline smoking status (current, previous, or never smoker), and baseline body mass index (BMI) category (underweight, healthy weight, overweight, obese). BMI was calculated according to the established formula (BMI=weight [pounds] x 703 / height^2 [inches]). The BMI categories were defined as underweight (BMI < 18.5), healthy weight (BMI=18.5 to < 25.0), overweight (BMI=25 to < 30.0), and obese (BMI ≥ 30). 33 Years of age were used as the time metric in the trajectory models, therefore, age was not included as a potential predictor of trajectory group membership.
Statistical Analysis
We used group-based trajectory modeling (GBTM) to identify subgroups that followed distinct patterns of joint change in trajectories of multimorbidity, functional status, cognitive performance, and depressive symptoms. Briefly, GBTM is a semi-parametric, finite mixture modeling approach designed to identify clusters of individuals who follow similar trajectories over time regarding some outcome of interest. 34 In the current study, we employ a recent generalization of GBTM, group-based multidimensional trajectory modeling, which identifies clusters of individuals who share common trajectories across multiple longitudinal indicators. 35 Model parameters are estimated via maximum likelihood and are robust to missing data under the missing at random (MAR) assumption.
Based on the empirical distributions of the outcome measures, we selected the most appropriate link function for each outcome from those currently implemented in GBTM software. 34 Specifically, we used a censored normal model specification for multimorbidity and cognitive performance, and a zero-inflated Poisson model for functional status and depressive symptoms. We modeled the trajectory of each outcome as a quadratic function of age. Because the model search process is computationally challenging when examining trajectories across multiple outcomes, we used a two-step approach to identify the best-fitting multidimensional trajectory model. 35 First, we fit a series of separate GBTM models for each outcome, successively increasing the number of trajectory groups from one to six, then we identified the optimal number of groups for each outcome using the following criteria: 1) reduction in Bayesian information criterion (BIC); 2) average posterior probability of group membership > 0.7 for all groups; 3) odds of correct classification > 5.0 for all groups; 4) group size >5% of total sample. 35 Next, we fit multidimensional trajectory models across the range of optimal group numbers identified in the first step and selected the final model by comparing the fit indices and qualitatively evaluating whether each increase in group number resulted in the identification of a trajectory group that was substantively distinct from those previously identified.
After selecting the optimal multidimensional trajectory model, we assigned each respondent to the trajectory group for which they had the highest posterior probability of membership (modal assignment). We calculated descriptive statistics stratified by trajectory group and conducted multinomial logistic regression to examine the association of baseline factors with assigned trajectory group membership. As our primary factor of interest, we examined the association of race/ethnicity with group membership in an unadjusted model. We then adjusted this model for the predictors described above and examined interactions between race/ethnicity and sex, education, and net worth. To assess for potential bias in either coefficient estimates or standard errors introduced by the probabilistic nature of trajectory group assignment, we conducted sensitivity analyses comparing the results of the multinomial regression model with identically specified models where: 1) group membership was determined by modal assignment and respondents were weighted using the posterior probability of assigned group membership, and 2) group membership was determined by proportional assignment, where respondents were assigned to each trajectory group in proportion to their posterior probability of group assignment. Analyses were conducted in Stata 17, with multidimensional trajectory models fit using the ‘traj’ add-on package. 36
Results
Sample characteristics
Baseline descriptive characteristics for participants overall and by multidimensional trajectory groups, Health & Retirement Study 1998-2016.
Abbreviations: SD= standard deviation; IQR= interquartile range.
Multidimensional trajectories
Fit statistics for the three-, four-, and five-group multidimensional trajectory model, Health and Retirement Study 1998-2016
The four groups of distinct joint trajectories of multimorbidity accumulation, functional status, cognitive performance, and depressive symptoms are displayed in Figure 1 and the estimated trajectory parameters (intercepts, linear terms, and quadratic terms) are provided in Table 3. Minimal Impairment with Low Multimorbidity (Group 1), representing 32.7% (n=3771) of the sample, had, on average, low rates of chronic disease at age 65 and little accumulation over time; no ADL/IADL limitation at age 65 and low levels of limitation in late life; high cognitive performance at age 65 that, while declining, does not cross the threshold for cognitive impairment by age 90; and almost no depressive symptomatology throughout the observed age span. Multidimensional trajectories of functional limitations, number of chronic conditions, cognitive performance, and depressive symptoms, Health and Retirement Study 1998-2016. Group-based multitrajectory model parameter estimates, Health and Retirement Study 1998-2016 Abbreviations: ADL = activities of daily living; IADL = instrumental activities of daily living; CES-D = Center for Epidemiological Studies Depression scale. *p < 0.05, **p < 0.01, ***p < 0.001.
Minimal Impairment with High Multimorbidity (Group 2; n=3878; 32.9%), is distinguished from the minimal impairment group by a higher multimorbidity level from the outset and throughout the observed age span. This group reaches the multimorbidity threshold by age 65, on average, and has greater accumulation of chronic conditions with increasing age. However, despite higher multimorbidity burden, members of this group developed few ADL/IADL limitations, maintained cognitive performance, and had few depressive symptoms across the observed age span.
Older adults in Multidomain Impairment with Intermediate Multimorbidity (Group 3; n=2291; 20.2%) had higher average multimorbidity burden than Group 1, Minimal Impairment with Low Multimorbidity—but lower average multimorbidity burden than Group 2, Minimal Impairment with High Multimorbidity. This group also displayed poor cognitive performance at age 65, rapid cognitive decline with advancing age, and high depressive symptomatology across the observed age span. Notably, while this group had, on average, good ADL/IADL function in early aging (65-75 years), the number of functional limitations sharply increases after age 75.
Finally, 14.2% (n=1625) of the sample comprised Multidomain Impairment with High Multimorbidity (Group 4). Respondents in this group had, on average, high levels of multimorbidity—more than three of seven chronic conditions at age 65 and the highest number of chronic conditions by age 90. Similar to the Multidomain Impairment with Intermediate Multimorbidity, members of this group demonstrated cognitive performance that was, on average, lower at age 65 than in the two low impairment groups, and exhibited rapid declination with advancing age. Compared with all other groups, this group experienced the greatest ADL/IADL functional limitation across the entire observed age span. Members of this group had, on average, two existing ADL/IADL limitations at age 65 and six by age 90. Additionally, this group had, on average, the highest number of depressive symptoms at age 65 and across the observed age span.
Multidimensional trajectory group membership characteristics
Sociodemographic and health-related characteristics for each multidimensional trajectory group are also presented in Table 1. There was a disproportionate number of Black and Hispanic respondents in both the Multidomain Impairment with Intermediate Multimorbidity and the Multidomain Impairment with High Multimorbidity groups. Overall, respondents from minoritized racial/ethnic backgrounds, low educational attainment, low wealth, current or prior smokers, and obese respondents were overrepresented in the two Multidomain Impairment trajectory groups (Groups 3 and 4). Tables A1-2 present outcome and chronic disease means and medians at selected ages—65, 75, and 90 years of age—for each of the trajectory groups. Mean outcome levels were comparatively worse at each reported age across outcomes for the Multidomain Impairment with High Multimorbidity group, especially when compared with the Minimal Impairment with Low Multimorbidity group (Table A1). Equivalently aged participants in the Multidomain Impairment with High Multimorbidity group had much higher prevalence of each of the chronic diseases compared with the other three groups (Table A2). The Minimal Impairment with High Multimorbidity group had higher prevalence across all diseases and particularly for heart disease, hypertension, diabetes, lung disease, and cancer compared with the Minimal Impairment with Low Multimorbidity and the Multidomain Impairment with Intermediate Multimorbidity groups.
Multidimensional trajectory group membership predictors
In unadjusted and adjusted multinomial logistic regression models, the healthiest and most functional group, Minimal Impairment with Low Multimorbidity (Group 1), served as the reference group. Consistent with the high posterior probabilities observed for assigned class membership, sensitivity analyses examining the impact of uncertainty in class assignment found negligible differences in the magnitude and statistical significance of coefficients between weighted and unweighted models, and, therefore, we report the results from unweighted models using modal class assignment. Compared with White respondents, Black respondents had greater odds of being in Multidomain Impairment with Intermediate Multimorbidity (Group 3) (OR=3.59, 95% CI=3.06,4.21) or Multidomain Impairment with High Multimorbidity (Group 4) (OR=4.30, 95% CI=3.63,5.09) than being in the reference group in the unadjusted model. Similarly, compared with White respondents, Hispanic respondents had greater odds of being in Multidomain Impairment with Intermediate Multimorbidity (Group 3) (OR=3.0, 95% CI=2.48,3.64) or Multidomain Impairment with High Multimorbidity (Group 4) (OR=3.29, 95% CI=2.68,4.04) in the unadjusted model.
Multinomial logistic regression of multidimensional trajectory group membership, Health and Retirement Study 1998-2016
Notes. Reference trajectory group is Minimal Impairment with Low Multimorbidity group. Referent wealth quartile 4 represents highest wealth quartile.
AbbreviationsOR = odds ratio; CI = confidence interval; NH = non-Hispanic.
* p < 0.05, ** p < 0.01, *** p < 0.001.
After adjusting for covariates, female respondents were more likely than male respondents to be in Multidomain Impairment with High Multimorbidity (Group 4) (OR=1.33, 95% CI=1.15,1.53) rather than in Minimal Impairment with Low Multimorbidity (Group 1), although there were no significant associations between sex and membership in other trajectory groups. Being coupled was associated with greater odds of being in Minimal Impairment with High Multimorbidity (Group 2) (OR=1.32, 95% CI=1.18,1.48) or Multidomain Impairment with High Multimorbidity (Group 4) (OR=1.19, 95% CI=1.029,1.38) relative to the reference.
Educational attainment and household wealth were both significant factors of group membership after controlling for other variables in the model. Relative to high school graduates, respondents with greater than a high school education were significantly less likely to be in Multidomain Impairment with Intermediate Multimorbidity (Group 3) (OR=0.62, 95% CI=0.53,0.73) or Multidomain Impairment with High Multimorbidity (Group 4) (OR=0.71, 95% CI=0.59,0.87). Conversely, respondents who did not finish high school were less likely to be in Minimal Impairment with High Multimorbidity (Group 2) (OR=0.83, 95% CI=0.72,0.96) but more likely to be in Multidomain Impairment with Intermediate Multimorbidity (Group 3) (OR=2.25, 95% CI=1.96,2.58) or Multidomain Impairment with High Multimorbidity (Group 4) (OR=2.12, 95% CI=1.82,2.48). There was a significant wealth gradient in group membership, with lower wealth associated with significantly greater odds of membership in either of the two Multidomain Impairment groups (Groups 3 and 4). Compared with the wealthiest respondents, persons in the lowest wealth quartile had 134% greater odds of being in Multidomain Impairment with Intermediate Multimorbidity (Group 3) (OR=2.34, 95% CI=1.94,2.82) and 494% greater odds of being in Multidomain Impairment with High Multimorbidity (Group 4) (OR=5.94, 95% CI=4.72,7.47).
Smoking history and increased BMI were both significantly associated with greater likelihood of being in the Multidomain Impairment groups (Groups 3 and 4). Compared with never smokers, both previous and current smokers had greater odds of being in Multidomain Impairment with Intermediate Multimorbidity (Group 3) and Multidomain Impairment with High Multimorbidity (Group 4), with current smokers having 96% greater odds and 202% greater odds, respectively. Of all the characteristics included in the fully adjusted model, obesity had the strongest association with membership in Multidomain Impairment with High Multimorbidity (Group 4), with obese respondents having markedly greater odds (OR=8.37, 95% CI=6.98,10.01) than healthy weight respondents of being in this group.
Discussion
This study broadens the understanding of the relationships between several important and interrelated health domains in late life. Our study identified four distinct latent classes of multidimensional trajectories of chronic disease, functional status, cognitive performance, and emotional wellbeing which, when examined in combination, reflected varying courses of multimorbidity burden and impairment. Older adults in the Minimal Impairment with Low Multimorbidity group were consistent with a healthy aging trajectory along the multiple domains of health examined in this study. Respondents in this group had high levels of function across multiple domains into late life. Similarly, older adults in the Minimal Impairment with High Multimorbidity fared well on many domains of health as they aged despite experiencing high multimorbidity burden.
Notably, while the Minimal Impairment with High Multimorbidity group had high and rising levels of multimorbidity, members appeared resilient to these high disease burdens because they exhibited preserved functional and cognitive status and experienced low depressive symptomatology well into old age. Relative to the referent minimally impaired group, membership in this group was associated with sociodemographic characteristics such as having lower wealth, and living in a coupled partnership, as well as lifestyle factors such as smoking and being overweight or obese. It is interesting that these correlates hint at mechanisms for high levels of multimorbidity—having fewer material resources, smoking, and obesity have been linked to the development of several chronic diseases 37 —but also suggest potential compensatory mechanisms—such as being in a coupled partnership—that may be associated with observed resilience on other health domains. While members may experience high disease burdens, the social support gained from these partnerships may provide pathways to functional, cognitive, and emotional health maintenance in the face of managing multiple chronic conditions. Future work should examine coupled partnerships more fully, including assessing the role of relationship quality and disruptions to partnerships on multidimensional aspects of health.
In contrast to the two Minimal Impairment groups, the remaining two latent classes had demonstrably worse trajectories of aging across multiple health domains. The Multidomain Impairment with Intermediate Multimorbidity group was characterized by poor cognitive performance throughout the ages observed, and a rapid uptick in the difficulty to perform daily life functions. Collectively, these multiple trajectories reinforce that impaired ADL/IADL function in later life may be associated with cognitive decline for this group. Finally, the Multidomain Impairment with High Multimorbidity group fared comparatively worse on every domain examined from age 65 onward. Members of this group demonstrated high levels of chronic somatic and mental health morbidity and poor function—cognitive as well as daily life function—throughout older age. The high levels of impairment across all four outcomes throughout the observed age span suggest that exposures earlier in the lifecourse (early life, early adulthood, and midlife) should be examined to further understand critical periods in these multidimensional courses of health and wellbeing into old age. This group was also more likely to have obese BMI at baseline and had the youngest average age at death. Members of this group may benefit from earlier life interventions to prevent obesity and chronic conditions associated with obesity and earlier mortality.
In multinomial logistic regression models comparing the four identified groups, we found that socially disadvantaged and minoritized populations were more likely to be members of the Multidomain Impairment groups. Black and Hispanic older adults were significantly more likely than White older adults to be members of either of the two Multidomain Impairment multidimensional trajectory groups that experienced high levels and rates of impairment and disease accumulation. Similarly, older adults with low socioeconomic status were significantly more likely to follow trajectories of substantial impairment across multiple health domains into old age. These vulnerable population groups are at higher risk of multidomain decline and as such, should be the focus of clinical interventions and public health programs aimed at modulating these declines.
Our study reinforces the importance of identifying changes across functional, physiological, and emotional health in older age.13,38–43 Wickrama and colleagues (2013) 42 assessed multidimensional health trajectories across chronic disease and functional domains, and found three latent classes characterized by maintained, persistently high, and deteriorating levels across these domains. In accordance with our study findings, socioeconomically disadvantaged and racial/ethnic minoritized groups were more likely to belong to adverse trajectory classes, reflecting significant health inequalities in late life. Also largely in alignment with our findings, Xu et al. (2015) 43 assessed cumulative changes across multiple physiological systems with the evaluation of multidimensional trajectories across physiological and functional health domains. They identified four distinct patterns of joint changes in functional, emotional, and cognitive performance and found stark educational differences in later old age, suggestive of social stratification in aging. Our findings extend beyond these studies by evaluating trajectories across a long follow-up period, and evaluating the intertwined courses of multimorbidity and common geriatric syndromes (disability and cognition) for older adults.
The study has several strengths. First, it contributes to the rapidly evolving multimorbidity literature by evaluating concurrent changes in key health domains. By identifying distinct classes of older adults who experience changes in their somatic multimorbidity burden alongside critical aspects to “what matters most” to older adults, our findings highlight heterogeneity across multiple markers of health in late life that need to be addressed. Second, the HRS is a large, ongoing, prospective health interview survey representative of the U.S. population of middle-aged and older adults. By virtue of the long follow-up time, multiple repeated observations for a large sample of Americans permits the study of longitudinal changes of these important health domains. Third, the prospective design and oversampling of minoritized older adults in HRS enables the examination of group membership characteristics for these groups.
Several limitations should also be noted. First, the data are self-reported. While this is a feature for many of the measures employed in these analyses (i.e., assessed cognitive performance, and reported depressive symptoms), chronic disease self-reports may be subject to misreporting. However, multiple studies have shown adequate concordance between patient self-reports and objectively-ascertained diagnoses.44–46 Second, we were limited by the number of chronic diseases assessed in the HRS. These diagnoses cover many common aging-related conditions, but there may be important distinctions with examination of more specific diagnoses. Future studies should also evaluate multidimensional trajectories that include both more specific diagnoses (e.g., ischemic heart disease) and a wider range of pertinent chronic diseases consistent with multimorbidity measurement frameworks.47,48 On a related matter, the biennial interview schedule may not sufficiently capture important transitions and end-of-life terminal declines. Additional work should evaluate functional changes in data sources that conduct more frequent assessments. Third, there is substantial heterogeneity in the Hispanic ethnicity category in country of origin or ancestry and which is associated with substantial divergent health status between Hispanic subgroups. 49 Fourth, there were insufficient numbers of other racial and ethnic groups to provide adequate power for modeling. Assessment of multidimensional trajectories among understudied racial/ethnic groups and within large racial/ethnic classifications should be examined in data sources that permit these evaluations. Furthermore, the exclusion of respondents who contributed fewer than three waves of data to facilitate the modeling of non-linear trajectories may have resulted in a sample biased toward more healthy respondents, and the presence of differential attrition rates across identified trajectory groups may have resulted in biased estimates of trajectory group size. 50 Finally, these analyses are not intended to assess causal effects. Rather, the multinomial regression analyses identify between-group differences in key sociodemographic and health-related factors. There may be important mediational relationships between included model covariates such as BMI and education or income that should be the focus of future work.
Our study has important implications. It suggests that high multimorbidity burden can be present with heterogeneous impacts on other aspects of aging in American older adults. While some distinct classes of older adults may experience high multimorbidity, there are two divergent pathways characterized by either resilience or substantial impairment throughout old age in relation to these high burdens. Our findings highlight the importance of targeting and remediating poor progression of multimorbidity, depressive symptoms, cognitive, and functional health domains for minoritized and low socioeconomic older adults in the US. This involves addressing upstream social factors and structural inequities that precede downstream morbidity, functional, cognitive, and affective health outcomes. These results would logically support early interventions to prevent the development and worsening of chronic disease burden, depressive symptoms, as well as cognitive and functional impairment, and involve not only healthcare interventions, but also those targeting key social determinants of health, for example, improving the quality of early-life education particularly in lower-income communities, and facilitating engagement with community-based services and supports into old age. Finally, our study points toward the need for future research to identify factors that enable resilience to high multimorbidity, particularly modifiable factors amenable to therapeutic intervention. There is a need to prioritize healthy aging initiatives that improve early detection, prevention, management, and support for physiological, emotional, cognitive, and functional health needs, particularly among minoritized and low socioeconomic groups.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Institute on Aging at the National Institutes of Health (Grant Number: R01AG055681 and RF1AG05845 to ARQ). Additional funding: National Institute on Aging/National Institutes of Health (R01AG047891 to HGA who also contributed from the Yale Claude D. Pepper Older Americans Independence Center P30AG021342 and P30AG066508; the University of Michigan Claude D. Pepper Older Americans Independence Center AG024824 and the Michigan Institute for Clinical and Health Research UL1TR000433 support to AB; and the Oregon Alzheimer’s Disease Research Center P30AG008017, P30AG066518 and P30AG024978 to JK). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Appendix
Mean and median values for outcomes by multidimensional trajectory group at selected ages, Health and Retirement Study 1998-2016
Minimal Impairment with Low Multimorbidity
Minimal Impairment with High Multimorbidity
Multidomain Impairment with Intermediate Multimorbidity
Multidomain Impairment with High Multimorbidity
Mean (SD)
Median (IQR)
Mean (SD)
Median (IQR)
Mean (SD)
Median (IQR)
Mean (SD)
Median (IQR)
Chronic Conditions
Age 65
0.6(0.55)
1(1)
2.3(0.86)
2(1)
1.33(0.70)
1(1)
3.34(1.19)
3(1)
Age 75
1.26(0.78)
1(1)
3.14(0.90)
3(1)
2.02(0.80)
2(1)
3.95(1.08)
4(2)
Age 90
2.1(0.98)
2(2)
4.01(0.77)
4(0)
2.72(0.87)
3(1)
4.65(0.98)
5(1)
ADL/IADL Limitations
Age 65
0.01(0.11)
0(0)
0.05(0.24)
0(0)
0.21(0.57)
0(0)
1.52(2.12)
1(2)
Age 75
0.05(0.27)
0(0)
0.14(0.47)
0(0)
0.7(1.24)
0(1)
2.68(2.55)
2(3)
Age 90
0.75(1.30)
0(1)
1.51(1.92)
1(2)
4.34(2.35)
4(3.5)
5.96(2.69)
5(3)
HRS Cognitive Scale
Age 65
18.08(3.33)
18(4)
17.14(3.44)
17(5)
13.78(4.15)
14(6)
13.63(4.46)
14(7)
Age 75
16.06(3.65)
16(5)
15.44(3.53)
16(5)
11.78(4.37)
12(6)
11.45(4.66)
12(7)
Age 90
11.64(3.92)
12(6)
11.73(3.86)
12(5)
7.67(4.40)
8(6)
7.35(3.23)
8(4)
CES-D
Age 65
0.41(0.85)
0(1)
0.71(1.10)
0(1)
1.58(1.59)
1(2)
2.49(1.90)
2(3)
Age 75
0.55(0.92)
0(1)
0.88(1.32)
0(1)
1.89(1.77)
1(3)
2.85(1.96)
3(3)
Age 90
1.18(1.42)
1(2)
1.4(1.46)
1(2)
2.57(1.74)
2(3)
2.98(1.92)
3(3)
Chronic disease prevalence by multitrajectory group and age, Health and Retirement Study 1998-2016
Minimal Impairment with Low Multimorbidity
Minimal Impairment with High Multimorbidity
Multidomain Impairment with Intermediate Multimorbidity
Multidomain Impairment with High Multimorbidity
Total
Heart disease
n=yes/total n (%)
n=yes/total n (%)
n=yes/total n (%)
n=yes/total n (%)
n=yes/total n (%)
Age 65-66
63/1696 (3.71)
704/2495 (28.22)
70/1117 (6.27)
546/1083 (50.42)
1383/6391 (21.64)
Age 75-76
277/2255 (12.28)
1161/2299 (50.50)
312/1424 (21.91)
570/862 (66.13)
2320/6840 (33.92)
Age 85-86
371/1393 (26.63)
565/752 (75.13)
266/600 (44.33)
175/202 (86.63)
1377/2947 (46.73)
Hypertension
Age 65-66
361/1696 (21.29)
1864/2495 (74.71)
558/1117 (49.96)
923/1083 (85.23)
3706/6391 (57.99)
Age 75-76
990/2255 (43.90)
2000/2299 (86.99)
919/1424 (64.54)
798/862 (92.58)
4707/6840 (68.82)
Age 85-86
810/1393 (58.15)
694/752 (92.29)
454/600 (75.67)
196/202 (97.03)
2154/2947 (73.09)
Stroke
Age 65-66
3/1696 (0.18)
97/2495 (3.89)
17/1117 (1.52)
198/1083 (18.28)
315/6391 (4.93)
Age 75-76
41/2255 (1.82)
257/2299 (11.18)
94/1424 (6.60)
266/862 (30.86)
658/6840 (9.62)
Age 85-86
77/1393 (5.53)
163/752 (21.68)
114/600 (19.00)
92/202 (45.54)
446/2947 (15.13)
Diabetes
Age 65-66
30/1696 (1.77)
582/2495 (23.33)
103/1117 (9.22)
503/1083 (46.45)
1218/6391 (19.06)
Age 75-76
145/2255 (6.43)
810/2299 (35.23)
254/1424 (17.84)
496/862 (57.54)
1705/6840 (24.93)
Age 85-86
106/1393 (7.61)
297/752 (39.49)
98/600 (16.33)
101/202 (50.00)
602/2947 (20.43)
Lung Disease
Age 65-66
14/1696 (0.83)
259/2495 (10.38)
41/1117 (3.67)
307/1083 (28.35)
621/6391 (9.72)
Age 75-76
83/2255 (3.68)
387/2299 (16.83)
105/1424 (7.37)
294/862 (34.11)
869/6840 (12.70)
Age 85-86
73/1393 (5.24)
183/752 (24.34)
58/600 (9.67)
58/202 (28.71)
372/2947 (12.62)
Arthritis
Age 65-66
517/1696 (30.48)
1863/2495 (74.67)
696/1117 (62.31)
983/1083 (90.77)
4059/6391 (63.51)
Age 75-76
1079/2255 (47.85)
1964/2299 (85.43)
1091/1424 (76.62)
820/862 (95.13)
4954/6840 (72.43)
Age 85-86
888/1393 (63.75)
707/752 (94.02)
516/600 (86.00)
201/202 (99.50)
2312/2947 (78.45)
Cancer
Age 65-66
57/1696 (3.36)
443/2495 (17.76)
55/1117 (4.92)
223/1083 (20.59)
778/6391 (12.17)
Age 75-76
256/2255 (11.35)
689/2299 (29.97)
169/1424 (11.87)
250/862 (29.00)
1364/6840 (19.94)
Age 85-86
241/1393 (17.30)
332/752 (44.15)
123/600 (20.50)
79/202 (39.11)
775/2947 (26.30)
