Skip to main content

Evaluation of T-cell aging-related immune phenotypes in the context of biological aging and multimorbidity in the Health and Retirement Study

Abstract

Background

Cellular changes in adaptive immune system accompany the process of aging and contribute to an aging-related immune phenotype (ARIP) characterized by decrease in naïve T-cells (TN) and increase in memory T-cells (TM). A population-representative marker of ARIP and its associations with biological aging and age-related chronic conditions have not been studied previously.

Methods

We developed two ARIP indicators based on well understood age-related changes in T cell distribution: TN/(TCM (Central Memory) + TEM (Effector Memory) + TEFF (Effector)) (referred as TN/TM) in CD4 + and CD8 + T-cells. We compared them with existing ARIP measures including CD4/CD8 ratio and CD8 + TN cells by evaluating associations with chronological age and the Klemera Doubal measure of biological age (measured in years) using linear regression, multimorbidity using multinomial logistic regression and two-year mortality using logistic regression.

Results

CD8 + TN and CD8 + TN/TM had the strongest inverse association with chronological age (beta estimates: -3.41 and -3.61 respectively; p-value < 0.0001) after adjustment for sex, race/ethnicity and CMV status. CD4 + TN/TM and CD4 + TN had the strongest inverse association with biological age (β = -0.23; p = 0.003 and β = -0.24; p = 0.004 respectively) after adjustment for age, sex, race/ethnicity and CMV serostatus. CD4/CD8 ratio was not associated with chronological age or biological age. CD4 + TN/TM and CD4 + TN was inversely associated with multimorbidity. For CD4 + TN/TM, people with 2 chronic conditions had an odds ratio of for 0.74 (95%CI: 0.63–0.86 p = 0.0003) compared to those without any chronic conditions while those with 3 chronic conditions had an odds ratio of 0.75 (95% CI: 0.63–0.90; p = 0.003) after adjustment for age, sex, race/ethnicity, CMV serostatus, smoking, and BMI. The results for the CD4 + TN subset were very similar to the associations seen with the CD4 + TN/TM. CD4 + TN/TM and CD4 + TN were both associated with two-year mortality (OR = 0.80 (95% CI: 0.67–0.95; p = 0.01) and 0.81 (0.70–0.94; p = 0.01), respectively).

Conclusion

CD4 + TN/TM and CD4 + TN had a stronger association with biological age, age-related morbidity and mortality compared to other ARIP measures. Future longitudinal studies are needed to evaluate the utility of the CD4 + subsets in predicting the risk of aging-related outcomes.

Introduction

Aging is a complex process accompanied by changes in the immune system which might result in the acceleration of biological decline and incidence of chronic diseases [1]. Immune system aging reduces the ability of the immune system to mount effective immune responses against new infections, vaccinations, and antigens [2, 3]. Previous studies have consistently reported that aging-related changes in T cells [4,5,6] contribute to an aging-related immune phenotype (ARIP) that is characterized by a reduced T-cell repertoire, reduction of naïve T cells (TN), and accumulation of memory (TM) and effector T cells [7, 8]. Though the determinants of ARIP remain incompletely understood, age, sex, and exposure to cytomegalovirus (CMV) are major determinants of T cell subsets [9,10,11,12]. In our previous study [13], we evaluated the associations of individual T-cell subsets measured in the Health and Retirement Study with chronological age. The study showed that CD4 + and CD8 + TN decreased with age and CD4 + TN cells were substantially higher among CMV seronegative individuals and women. We observed that CD8 + TEFF, CD4 + TEFF, and CD4 + TEM cells were strongly influenced by CMV seropositivity, and CD8 + TEM increased with age. CD8 + and CD4 + TCM decreased with age although CD4 + TCM was higher among women compared to men and CD8 + TCM was higher among men. These results show that the inter-individual variability in T cell subsets is determined by several factors such as age, sex and CMV serostatus and individual T cell subsets are differentially influenced by the various factors. Optimal T-cell immunity involves harmonized action between multiple cell populations and a population level metric combining the cell populations will be biologically meaningful in understanding the aging immune system. In this regard, CD4/CD8 ratio is a well-established ARIP measure where the prevalence of an inverted CD4/CD8 ratio increases with age [14, 15]. However, we have shown previously that the prevalence of the inverted CD4/CD8 ratio increases with age only among CMV seropositive individuals and not among CMV seronegative individuals suggesting that this ratio is not a universal marker of ARIP as previously thought [13]. Thus, other measures of ARIP, beyond the inverted CD4/CD8 ratio, that are independent of CMV serostatus need to be identified. To this end, one study identified a cellular composite measure (IMM-AGE) that was associated with mortality, and another study developed an inflammatory aging clock (iAGE) utilizing soluble systemic chronic inflammation markers or cytokines which was associated with multimorbidity [16, 17]. However, the interpretation of cytokine levels is complicated by the fact that that circulating cytokine levels are influenced by secretion from multiple cell types and the IMM-Age measure that utilized a combination of flow cytometry and mass cytometry estimated differences in both innate and adaptive immune systems. However, since age-related and CMV-related changes are seen predominantly in the T-cell subsets, a composite measure focused on T cells and relatively easy to implement in large epidemiological studies is needed.

To address this gap, we used T-cell immunophenotyping data that measured 11 T cell subsets in a large, population-based sample and used well understood aging-related changes in T cell distribution to develop aging-related immune phenotype (ARIP) indicators that would be robust across sexes, race/ethnic groups and CMV serostatus. We developed ARIP measures based on a priori knowledge of (i) the role of chronological age in the immune system and (ii) the role of underlying biological mechanisms in the immune system such as CMV which increases susceptibility to age-related chronic conditions [18,19,20].We also evaluated the correlation between individual T-cell subsets and association of individual T-cell subsets with chronological age, biological age and multimorbidity to develop composite ARIP measures. We created two candidate measures: TN/ (TCM + TEM + TEFF) (referred as TN/ TM) in CD4 + and CD8 + T cells and benchmarked them against CD4/CD8 ratio, CD8 + TN and CD4 + TN by evaluating associations of all these measures with chronological age, biological age, multimorbidity (that includes several aging-related diseases such as diabetes, cancer, lung disease, stroke, and heart disease), and mortality.

Methods

Study Population

The Health and Retirement Study (HRS), supported by National Institute of Aging (NIA), is a nationally representative longitudinal survey of adults in the United States over the age of 50. The study began in 1992 with participants 51–61 years of age and now uses a steady-state design to replenish the sample every 6 years with younger cohorts [21, 22] with surveys/interviews being conducted biennially. The HRS sample employs a multi-stage area probability design with oversampling of African-American and Hispanic households at about twice the rate of Whites [21]. Whole blood was collected from HRS participants in the 2016 wave as part of the Venous Blood Study (VBS) and immunophenotyping was performed on samples from 9932 participants. After removing participants with missing age, sex, race/ethnicity, CMV seroprevalence, Venous Blood Study survey weights, and T cell data, 8603 participants were included in the study. After weighting the sample using Venous Blood Study survey weights, it is representative of a U.S older national population.

Immunophenotyping

Peripheral Blood Mononuclear cells (PBMCs) were isolated from whole blood and used for the measurement of immune cells [23]. All flow cytometry measurements were performed on an LSRII flow cytometer or a Fortessa X20 instrument (BD Biosciences, San Diego, CA). Immunophenotyping data was analyzed using OpenCyto and FlowAnnotator as described previously [24]. T-cells and ten T-cell subsets were evaluated in HRS. Panel 1 consisted of antibody cocktails targeting T cells and T cell subsets (Supplementary Table 1). T cell subsets were represented as percentages of their parent population. The detailed methods and cell subset definitions have been mentioned previously [13].

Ascertainment of multimorbidity and mortality

Prevalence of type II diabetes, stroke, lung disease, heart disease, and cancer were used to generate a 2016 multimorbidity outcome. Type II diabetes was defined as fasting glucose value >  = 126 mg/dL or use of diabetes medications or insulin. Stroke, lung disease, heart disease and cancer were obtained from self-report. Multimorbidity was defined as a categorical variable with four categories (no prevalent chronic conditions (reference), a single chronic condition, two chronic conditions, and three or more chronic conditions). Incident multimorbidity in 2018 and 2020 was ascertained by defining a baseline study population with participants who had a multimorbidity score of 0 or 1 in 2016. Among the participants in this baseline population, incident multimorbidity was defined as a multimorbidity score of 2 or higher in 2018 and 2020 using self-report of Type II diabetes, stroke, lung disease, heart disease, and cancer. Mortality assessed at 2018 was used for this analysis. This measurement is based on HRS tracking efforts when there is no interview in the current wave, or an exit interview was not obtained in the prior wave. The reported vital status is based on reports from someone capable of reporting a death.

Ascertainment of biological age and self-rated health

Biological age was developed using the Klemera Doubal [25] method as a combination of biomarkers which represents the decline in aging-related physiological functioning and susceptibility to disease in old age. It was previously validated in the NHANES cohort as a more reliable predictor of mortality than chronological age and other biological age algorithms [26]. The biomarkers used for calculating the biological age were: Cardiovascular function (Systolic Blood pressure), Metabolic markers (Total Cholesterol and Fasting Glucose), Inflammation (CMV and C-Reactive Protein), Kidney function (Serum Creatinine and Blood Urea Nitrogen), Liver function (Alkaline phosphatase and Albumin), Lung function (Peak flow). Biological age was calculated using the R package BioAge and the in-built function kdm_calc [27]. Self-rated health was measured in the Health and Retirement Study by asking the question “Would you say your health is excellent, very good, good, fair, or poor?” We combined the “excellent”, “very good” and “good” categories to represent a good score for health and the “fair” and “poor” categories to represent a bad score for health.

Ascertainment of participant characteristics and biological age

Chronological age (in years), sex (female/male), race/ethnicity (Hispanic, Non-Hispanic Whites, Non-Hispanic Black, and Non-Hispanic Other) were obtained from the HRS demographic data [22]. Smoking was self-reported and categorized as never smokers, former smokers, and current smokers. Height was measured in inches and weight was measured in pounds. Height and weight for all the participants were obtained for half the sample in 2014 and the other half in 2016. BMI was calculated by using the formula (Weight in lbs/ (Height in inches)2) *703. CMV seroprevalence was measured using the ratio of total IgG to anti-CMV IgG in serum using the Roche e411 immunoassay analyzer (Roche Diagnostics Corporation). The results are reported as non-reactive (< 0.5 COI), borderline (0.5 to < 1.0 COI) or reactive >  = 1.0 COI) where COI is cut-off interval [28]. This was used as a binary variable by combining the borderline and negative groups.

Statistical analysis

We evaluated two candidate ARIP measures (TN/TM CD4 + and CD8 + T cells) to compare with existing CD4/CD8 ratio and the individual cell type CD8 + TN as we have shown previously that it had a strong inverse association with chronological age [13]. We studied associations of the candidate ARIP and existing ARIP measures with chronological age, biological age, and multimorbidity. CD4 + TN/TM, CD8 + TN/TM and CD4/CD8 ratio were log-transformed due to skewed distributions.

Pearson correlation coefficients were calculated between individual T-cell subsets. To study the association of ARIP measures with chronological age, and biological age, survey linear regression models were used after adjustment for chronological age, sex, race/ethnicity, and CMV status. A biological age acceleration measure was also calculated by subtracting chronological age from biological age. The participants were categorized into “higher” category if the age acceleration measure was higher than 0, and the participants were grouped into “lower” category if the age acceleration measure was lower than 0. Using the “lower” group as a reference, logistic regression analysis was performed between ARIP measures and biological age after adjustment for age, sex, race/ethnicity and CMV status. To study the association of ARIP measures with prevalent multimorbidity in 2016, multinomial survey logistic regression models were used. We studied the association between ARIP measures and the individual components of multimorbidity (type II diabetes, hypertension, lung disease, heart disease, and stroke) by using a series of survey logistic regression models after adjustment for age, sex, race/ethnicity, CMV status, smoking status, and BMI. We used survey cox regression models to study the association of ARIP measures in incident multimorbidity in 2018 and 2020. We used survey logistic regression models to study the association of the ARIP measures with self-rated health and mortality in 2018. The survey models were adjusted for survey design parameters including strata and cluster details for sampling error and participant sample weights from the Venous Blood Study to account for sample design. In sensitivity analyses, three additional characterizations of ARIP were explored and results are presented in the supplementary material: (i) CD4 + TN/CD8 + TN; (ii) CD4 + TN + TCM; (iii) CD8 TN + TCM.

We used SAS version 9.4 (SAS Institute, Inc., Cary, NC) and R Statistical Analysis software version 4.0.0 for all the analyses.

Results

Participant characteristics

Among the 8603 participants included, the average age of participants was 68.65 years. Fifty-four percent of the participants were women; 10.02% of participants were Non-Hispanic Black and 9.01% were Hispanic; 11.87% of the participants were current smokers and the average BMI was 29.69 kg/m2. Sixty-four percent of the participants were seropositive for cytomegalovirus. In 2016, 48.13% of participants were without any chronic conditions, 32.82% participants with one chronic condition, 14.06% participants with two chronic conditions, and 4.98% participants with three or more chronic conditions. Heart disease was the most prevalent (24.83%) among the diseases used to define multimorbidity (Table 1). Six percent of the participants had incident multimorbidity in 2018 and 2020, and four percent of the participants were not alive in 2018 (Table 1). As expected, CD4 + TN was negatively correlated with CD4 + TCM, CD4 + TEFF and CD4 + TEM (R = -0.72, -0.3 and -0.22, respectively) (Supplementary Fig. 1). CD8 + TN was negatively correlated with CD8 + TEFF and CD8 + TEM (R = -0.67 and -0.26, respectively) (Supplementary Fig. 1). Based on the negative correlation measures of TN with TCM, TEFF, and TEM, the composite ARIP measures CD4 + and CD8 + TN /(TCM + TEM + TEFF) (referred as TN/ TM) were developed.

Table 1 Descriptive statistics of Health and Retirement Study participant characteristics measured in 2016 among 8603 participants with immunophenotyping data

Association of ARIP measures with chronological age

Overall, the CD8 + subsets were more strongly associated with chronological age than the CD4 + subsets. Among the individual cell subsets, CD4 + and CD8 + TN were negatively associated with chronological age (β = -0.32; p = 0.05 and β = -3.41; p < 0.0001). CD4 + and CD8 + TEFF (β = 0.66; p < 0.0001 and β = 2.58; p < 0.0001), CD4 + and CD8 + TEM (β = 1.81; p < 0.0001 and β = 1.7; p < 0.0001) and CD8 + TCM (β = 0.39; p = 0.04) were positively associated with chronological age (Supplementary Fig. 2). Among the ARIP measures, univariate analysis showed that CD4/CD8, CD4 + TN/TM, and CD4 + TN had the weakest association with chronological age whereas CD8 + TN/TM and CD8 + TN had the strongest association with chronological age (Supplementary Fig. 3). A one standard deviation (SD) increase in CD8 + TN/TM (p < 0.001) was associated with a 3.61 years lower chronological age (Fig. 1) after adjusting for sex, race/ethnicity, and CMV status. CD4/CD8 and CD4 + TN /TM were not associated with chronological age (p = 0.13 and p = 0.84, respectively). The other ARIP measures also had an inverse association with chronological age (Supplementary table 2).

Fig. 1
figure 1

Association of ARIP measures with chronological age. Age is used as the dependent variable in survey linear regression models. The beta estimates are estimated per one standard deviation unit of the ARIP measures. The models are adjusted for sex, race/ethnicity, and CMV status. The solid black line along 0 indicates no association

Association of ARIP measures with biological age

Among the individual T-cell subsets, CD4 + TN was negatively associated with biological age (0.24 years decrease in biological age for one standard deviation (SD) increase in CD4 + TN (p = 0.001) (Fig. 2), one SD increase in CD4 + TCM and CD4 + TEFF was associated with 0.18 years (p = 0.01) and 0.13 years (p = 0.02) higher biological age respectively (Supplementary Fig. 4). Among the CD8 + T cell subsets, only the CD8 + TCM subset was positively associated with biological age (0.12 years higher biological age for one SD increase in CD8 + TCM (p = 0.03)) (Supplementary Fig. 4). Univariate analysis showed that CD4/CD8, CD4 + TN/TM, and CD4 + TN had the weakest association with biological age whereas CD8 + TN/TM and CD8 + TN had the strongest association with biological age (Supplementary Fig. 5). However, after adjustment for covariates, CD4 + TN/TM had a strong negative association with biological age with 0.23 years lower biological age for one standard deviation increase in CD4 + TN/TM (p-value = 0.003) (Fig. 2). CD4/CD8 ratio, CD8 + TN/TM, and CD8 + TN were not associated with biological age. CD4 + TN/TM,CD4 + TN and CD4 + TCM remained associated with biological age acceleration (older category). One SD increase in CD4 + TN/TM was associated with 8% lower odds of having higher biological age as compared to chronological age (OR: 0.92 (95% CI: 0.84 – 1.00); p-value = 0.05), one SD increase in CD4 + TN was associated with 11% lower odds of having higher biological age (OR: 0.89 (95% CI: 0.83 – 0.95; p = 0.001) and one SD increase in CD4 + TCM was associated with 9% higher odds of having higher biological age (OR: 1.09 (95% CI: 1.02–1.17; p = 0.01). CD8 + TN/TM was not associated with biological age acceleration (OR: 1.04; (95% CI: 0.98 – 1.11; p-value = 0.19) (Supplementary table 3). Additional adjustment for CRP did not substantially change the associations between ARIP measures and biological age (data not shown).

Fig. 2
figure 2

Association of ARIP measures with biological age. Biological age is used as the dependent variable in survey linear regression models. The beta estimates are estimated per one standard deviation unit of the ARIP measures. The models are adjusted for chronological age, sex, race/ethnicity, and CMV status. The solid black line along 0 indicates no association

Association of ARIP measures with multimorbidity

CD4 + TN had an odds ratio of 0.85 (95% CI: 0.79–0.92; p < 0.0001) when comparing those with a single chronic condition (multimorbidity score = 1) vs. those without any chronic condition (multimorbidity score = 0), while the corresponding odds ratio for those with two and 3 + chronic conditions was 0.72 (95% CI: 0.65 – 0.79; p < 0.0001) and 0.74 (95% CI: 0.64–0.87; p = 0.0004) respectively (Fig. 3). Among the other individual T-cell subsets, CD4 + TCM was consistently positively associated across different categories of multimorbidity whereas the other T-cell subsets did not have a consistent association across multimorbidity categories (Supplementary Fig. 6). CD4 + TN/TM had an odds ratio of 0.89 (95% CI: 0.82–0.96; p = 0.005) when comparing those with a single chronic condition vs. those without any chronic condition, while the corresponding odds ratio was 0.74 (95% CI: 0.63 – 0.86; p = 0.0003) and 0.75 (95% CI: 0.63–0.90; p = 0.003) for those with two and 3 + chronic conditions, respectively (Fig. 3). CD4/CD8, CD8 + TN/TM and CD8 + TN were not consistently associated across the different categories of multimorbidity (Fig. 3). Further stratifying the highest level of the multimorbidity score into 3 and 4 + did not change the observed associations (data not shown). Among the chronic diseases included in the multimorbidity score, ARIP measures had the strongest association with prevalent cancer after multivariate adjustment. Lung disease, diabetes, stroke, and heart disease were not individually associated with the ARIP measures after multivariate adjustment (Fig. 4, Supplementary Table 4). CD4 + TCM was associated with incident multimorbidity over four years follow-up (OR: 1.09 (95% CI: 1.01 – 1.18; p-value = 0.03). CD4 + TN/TM and CD4 + TN also had association with incident multimorbidity in the consistent direction (OR: 0.91 (95% CI: 0.82- 1.02); p-value = 0.09 and OR = 0.91 (95% CI: 0.82 – 1.01); p-value = 0.09) though the associations were not statistically significant. (Supplementary Table 4).

Fig. 3
figure 3

Odds ratios and 95% CI of association of multimorbidity levels with ARIP markers per one SD unit increase in ARIP marker. Adjusted for age, sex, race/ethnicity, CMV status, smoking status, and BMI. The solid black line along OR of 1 indicates no association

Fig. 4
figure 4

Odds ratios and 95% CI of association of individual components of multimorbidity ARIP markers per one SD unit increase in ARIP marker. Adjusted for age, sex, race/ethnicity, CMV status, smoking status, and BMI. The solid black line along OR of 1 indicates no association

Association of ARIP measures with self-rated health and two-year mortality

All the ARIP measures were associated with self-reported health (bad score for health vs good score for health) with the strongest association in CD4 + TN/TM, CD4 + TN and CD4 + TCM with OR of 0.84 (95% CI: 0.78 – 0.91), 0.83 (95% CI: 0.77 – 0.89) and 1.14 (95% CI: 1.04–1.25; p = 0.006), respectively (Supplementary table 5). The other individual T cell subsets were not associated with self-rated health (data not shown). CD4 + TN/TM and CD4 + TN had the strongest association with mortality with an OR of 0.80 (95% CI: 0.67 – 0.95; p = 0.01) and 0.81 (95% CI: 0.70–0.94; p = 0.01), respectively. The other ARIP measures were not associated with mortality in 2018 (Fig. 5). Among the individual T-cell subsets, only CD8 + TEFF (OR = 1.27 (95% CI: 1.09–1.49; p = 0.003) was associated with mortality in 2018.

Fig. 5
figure 5

Odds ratios and 95% CI of association of mortality with ARIP markers per one SD unit increase in ARIP marker. Adjusted for age, sex, race/ethnicity, CMV status, smoking status, and BMI. The solid black line along OR of 1 indicates no association

Discussion

This study explored the associations of ARIP measures with chronological age, biological age, and multimorbidity outcomes. CD4 + TN/TM and CD4 + TN had the strongest association with biological age and mortality and was inversely associated with increasing levels of multimorbidity. Though CD8 + TN/TM was not associated with biological age nor mortality, CD8 + TN/TM had associations with individual conditions including heart disease, diabetes and cancer. More examination is needed to better understand the relationship between CD8 + TN/TM and the individual chronic conditions. CD4/CD8 ratio was not associated with chronological age or biological age and did not have a consistent association with multimorbidity (Table 2). The associations with increasing levels of multimorbidity, mortality and biological age indicates the CD4 + TN/TM and CD4 + TN may be the biomarkers to identify individuals at higher risk for in age acceleration and its associated morbidities and higher mortality.

Table 2 Summarizing the associations of the candidate and existing ARIP measures with chronological age, biological age, and multimorbidity

CD4 + TN/TM and CD4 + TN had associations of similar strengths with biological age, multimorbidity and mortality. CD4 + TN and CD4 + TCM were strongly negatively correlated. While CD4 + TN was inversely associated with biological age and multimorbidity, CD4 + TCM was positively associated with biological age and multimorbidity. These observed associations indicate antagonistic roles of CD4 + TN and CD4 + TCM in biological aging and the pathogenesis of age-related chronic conditions. Hence, a representation of ARIP as a combination of CD4 + TN/TM may be a better representation of the overall associations between T cell immunity and health outcomes when compared to describing these associations with individual T-cell subsets. However, CD4 + TN also shows similar associations with biological age, multimorbidity and mortality as CD4 + TN/TM.

Among the two new measures developed in this study, CD8 + TN/TM was associated with chronological age but not associated with biological age and higher chronic conditions. CD4 + TN/TM was not associated with chronological age but had a strong inverse association with biological age and with a higher number of prevalent co-occurring chronic conditions. The findings with CD4 + TN/TM are consistent with previous studies where CD4 + T cell subsets were associated with chronic diseases. In the Multi-Ethnic Study of Atherosclerosis (MESA), naïve and memory CD4 + T cells were cross-sectionally associated with type II diabetes and subclinical atherosclerosis [29, 30]. CD4 + T cells produce interferon- γ (IFN- γ) which contributes to inflammation, glucose intolerance, and insulin resistance in diet-induced obesity mice [31, 32]. Another study performed among type II diabetes patients found a significant reduction in the naïve pool in both CD4 + and CD8 + populations and a significant rise in TEM and TEFF populations of the CD4 + subset compared to age‐matched controls [33]. Decreased CD4 + TN/TM ratio has also been observed in 76 non-small cell lung cancer (NSCLC) patients compared to 28 age and sex-matched healthy volunteers [34]. CD4 + CD28null cells which have advanced effector functions have been found at an increased frequency in ischemic stroke patients which was also associated with stroke severity [35]. These studies are consistent with our findings that the distribution of naïve and memory CD4 + T cells may be important in determining age-related outcomes. In contrast, the CD8 + TN cells were inversely associated with chronological age after adjustment for sex, race/ethnicity, and CMV serostatus but not associated with biological age and a higher number of chronic conditions. Previous studies have also used CD4/CD8 ratio as a measure of immune aging. An increased prevalence of the inverted CD4/CD8 ratio was associated with short-term mortality in the OCTO immune longitudinal study, possibly due to confounding by CMV serostatus [36]. CD4/CD8 ratio was not associated with biological age and chronological age after adjustment for race/ethnicity, sex, and CMV status and was not consistently associated with a higher prevalence of multiple chronic diseases.

Although increasing chronological age is an important component of senescence, it does not directly measure the accelerated decline in health among individuals. Previous studies have shown that biological age predicts morbidity and mortality independent of chronological age in individuals [26, 37,38,39]. This is the first study to demonstrate an association between biological age and CD4 + TN/TM, CD4 + TN and CD4 + TCM. In the Rotterdam study, longitudinal associations of biological age with all-cause morbidity, stroke, cancer, and diabetes mellitus, suggests that biological age may predetermine disease occurrence [39]. The strong association between multiple CD4 + subsets including CD4 + TN/TM and biological age suggests that CD4 + TN/TM might be an important risk factor for future disease occurrence and this needs to be confirmed in future studies.

Heterogeneity in strength of the association between the different components of multimorbidity and the ARIP measures were observed but the direction of the association was consistent among the diseases examined. CD4 + TN/TM and CD4 + TN were lower among participants who had one chronic condition (multimorbidity score = 1) compared to those without any chronic conditions (multimorbidity score = 0). Since cancer was the third most frequent chronic condition (after heart disease and diabetes), this may reflect the strong association between cancer and CD4 + T cells. The distribution of chronic conditions when participants had two chronic conditions (multimorbidity score = 2) is similar to the overall prevalence of chronic conditions. Since the CD4 + TN/TM and CD4 + TN CD4 subsets were associated with the three most common chronic conditions (heart disease, diabetes and cancer), this could explain the significant association between these subsests and multimorbidity score of two. However, when participants had three or more chronic conditions (multimorbidity score = 3), lung disease had a higher prevalence compared to when participants had two chronic conditions (multimorbidity score = 2). Though none of the ARIP measures were significantly associated with lung disease, CD4 + TN/TM and CD4 + TN had the strongest association (OR = 0.88 and 0.90 respectively) with lung disease. This may explain why only CD4 + TN/TM and CD4 + TN was associated with a multimorbidity score of 3. Among the disease conditions which constituted multimorbidity, cancer had the strongest association with the ARIP measures. Various theories have been proposed describing the involvement of naïve and memory CD4 + T cells in tumor immunity. The tumor microenvironment conditions promote the differentiation of memory and naïve CD4 + T cells into CD4 + T regulatory cells which suppress antitumor immunity [40,41,42]. The concept of immunosurveillance in cancer indicates that the immune system can recognize and proactively remove precursors of cancer and naïve T cells may be involved in this process as they perform immunosurveillance roles to undergo proliferation in response to homeostatic signals [43, 44]. In general, TCM cells have been found to produce higher levels of cytokines which have increased efficiency in an antitumor response [45, 46]. We also observed that the ARIP measures were not significantly associated with incident multimorbidity after adjustment for age, sex, race/ethnicity, CMV status, smoking status, and BMI though the direction of association was consistent with what was observed in the cross-sectional analysis. This could be due to the limited number of incident multimorbidity events in the short follow-up period and the association between the ARIP measures and incident multimorbidity needs to be further evaluated in future studies. CD4 + TN/TM and CD4 + TN were associated with two-year mortality as well strengthening their role in determining biological aging.

Major strengths of this study include the standardized and rigorous immunophenotyping methods used for measuring the cell populations, representative sample including Hispanic, black, and white individuals, immunophenotyping data available for a large sample size (about 9938 individuals), and concurrent measurement of cytomegalovirus seroprevalence. The novelty of this study is the creation of composite ARIP measures as a marker for an aspect of aging and characterizing the measures using chronological age, biological age, and age-related health outcomes. A limitation of this study is the absence of longitudinal immunophenotype data. Measurement of longitudinal changes in the immune cell subsets would help in studying the change in immune cells with age within the same participant and would help us to better understand the temporal relationship between immune aging and disease. The majority of the multimorbidity measure is based on the self-report of the participant and hence could be influenced by factors such as recall, social desirability, and lack of diagnosis. A previous study found that the sensitivity of self-report of hypertension was 88.9% among blacks, 82.8% for whites, and 84.0% for Hispanic ethnicity, and specificity was 92.8% for the whites and 86% for the blacks [47] informing that the measurements are robust. Another study performed to validate self-reported cancer measures with Medicare claims data had a 73.2% sensitivity and 96.2% specificity [48].

In conclusion, the results suggest that composite measures of age-related immune phenotypes are more meaningful in determining biological aging and multimorbidity. In this study, we showed that CD4 + TN/TM and CD4 + TN had strong associations with biological age, multimorbidity and mortality and maybe a more consistent measure of immune aging compared to specific cell types such as CD8 + naïve T cells and the widely used CD4/CD8. CD4 + TN/TM and CD4 + TN can be used in future studies with longitudinal measurements of the immune cells to better understand the role of immune aging in age-related morbidity and mortality.

Availability of data and materials

The datasets generated and/or analyzed during the current study are available in the HRS repository at https://hrsdata.isr.umich.edu/data-products/public-survey-data and https://hrsdata.isr.umich.edu/data-products/sensitive-health.

References

  1. Gruver AL, Hudson LL, Sempowski GD. Immunosenescence of ageing. J Pathol. 2007;211(2):144–56.

    Article  CAS  Google Scholar 

  2. Pawelec G. Age and immunity: What is “immunosenescence”? Exp Gerontol. 2017;2018(105):4–9.

    Google Scholar 

  3. Oh SJ, Lee JK, Shin OS. Aging and the immune system: The impact of immunosenescence on viral infection, immunity and vaccine immunogenicity. Immune Netw. 2019;19(6):1–18.

    Article  CAS  Google Scholar 

  4. Weng N-P. Aging of the immune system: how much can the adaptive immune system adapt? Immunity. 2006;24(5):495–9 (https://pubmed.ncbi.nlm.nih.gov/16713964).

    Article  CAS  Google Scholar 

  5. Hakim FT, Gress RE. Immunosenescence: deficits in adaptive immunity in the elderly. Tissue Antigens. 2007;70(3):179–89. https://doi.org/10.1111/j.1399-0039.2007.00891.x.

    Article  CAS  PubMed  Google Scholar 

  6. Crooke SN, Ovsyannikova IG, Poland GA, Kennedy RB. Immunosenescence and human vaccine immune responses. Immun Ageing. 2019;16(1):25. https://doi.org/10.1186/s12979-019-0164-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Aiello A, Farzaneh F, Candore G, Caruso C, Davinelli S, Gambino CM, et al. Immunosenescence and its hallmarks: How to oppose aging strategically? A review of potential options for therapeutic intervention. Front Immunol. 2019;10(SEP):1–19.

    Google Scholar 

  8. Bischof J, Gärtner F, Zeiser K, Kunz R, Schreiner C, Hoffer E, et al. Immune Cells and Immunosenescence. Folia Biol (Praha). 2019;65(2):53–63.

    CAS  Google Scholar 

  9. Saule P, Trauet J, Dutriez V, Lekeux V, Dessaint J-P, Labalette M. Accumulation of memory T cells from childhood to old age: Central and effector memory cells in CD4+ versus effector memory and terminally differentiated memory cells in CD8+ compartment. Mech Ageing Dev. 2006;127(3):274–81 (https://www.sciencedirect.com/science/article/pii/S0047637405002757).

    Article  CAS  Google Scholar 

  10. Li M, Yao D, Zeng X, Kasakovski D, Zhang Y, Chen S, et al. Age related human T cell subset evolution and senescence. Immun Ageing. 2019;16(1):24. https://doi.org/10.1186/s12979-019-0165-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016;16(10):626–38. https://doi.org/10.1038/nri.2016.90.

    Article  CAS  PubMed  Google Scholar 

  12. Sylwester AW, Mitchell BL, Edgar JB, Taormina C, Pelte C, Ruchti F, et al. Broadly targeted human cytomegalovirus-specific CD4+ and CD8+ T cells dominate the memory compartments of exposed subjects. J Exp Med. 2005;202(5):673–85. https://doi.org/10.1084/jem.20050882.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Thyagarajan B, Faul J, Vivek S, Kim JK, Nikolich-Žugich J, Weir D, et al. Age-Related Differences in T-Cell Subsets in a Nationally Representative Sample of People Older Than Age 55: Findings From the Health and Retirement Study. J Gerontol Ser A 2021;glab300. https://doi.org/10.1093/gerona/glab300

  14. Strindhall J, Skog M, Ernerudh J, Bengner M, Löfgren S, Matussek A, et al. The inverted CD4/CD8 ratio and associated parameters in 66-year-old individuals: the Swedish HEXA immune study. Age (Omaha). 2013;35(3):985–91. https://doi.org/10.1007/s11357-012-9400-3.

    Article  CAS  Google Scholar 

  15. Vasson M-P, Farges M-C, Goncalves-Mendes N, Talvas J, Ribalta J, Winklhofer-Roob B, et al. Does aging affect the immune status? A comparative analysis in 300 healthy volunteers from France, Austria and Spain. Immun Ageing. 2013;10(1):38. https://doi.org/10.1186/1742-4933-10-38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Alpert A, Pickman Y, Leipold M, Rosenberg-Hasson Y, Ji X, Gaujoux R, et al. A clinically meaningful metric of immune age derived from high-dimensional longitudinal monitoring. Nat Med. 2019;25(3):487–95. https://doi.org/10.1038/s41591-019-0381-y.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Sayed N, Huang Y, Nguyen K, Krejciova-Rajaniemi Z, Grawe AP, Gao T, et al. An inflammatory aging clock (iAge) based on deep learning tracks multimorbidity, immunosenescence, frailty and cardiovascular aging. Nat Aging. 2021;1(7):598–615. https://doi.org/10.1038/s43587-021-00082-y.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Spyridopoulos I, Martin-Ruiz C, Hilkens C, Yadegarfar ME, Isaacs J, Jagger C, et al. CMV seropositivity and T-cell senescence predict increased cardiovascular mortality in octogenarians: results from the Newcastle 85+ study. Aging Cell. 2016;15(2):389–92. https://doi.org/10.1111/acel.12430.

    Article  CAS  PubMed  Google Scholar 

  19. Maussang D, Verzijl D, van Walsum M, Leurs R, Holl J, Pleskoff O, et al. Human cytomegalovirus-encoded chemokine receptor US28 promotes tumorigenesis. Proc Natl Acad Sci U S A. 2006;103(35):13068–73.

    Article  CAS  Google Scholar 

  20. Chen S, de Craen AJM, Raz Y, Derhovanessian E, Vossen ACTM, Westendorp RGJ, et al. Cytomegalovirus seropositivity is associated with glucose regulation in the oldest old. Results from the Leiden 85-plus Study. Immun Ageing. 2012;9(1):18 (https://pubmed.ncbi.nlm.nih.gov/22929089).

    Article  Google Scholar 

  21. Heeringa SG, Connor J. Technical Description of the Health and Retirement Study Sample Design. Ann Arbor: Institute for Social Research, University of Michigan; 1995.

    Book  Google Scholar 

  22. Sonnega A, Faul JD, Ofstedal MB, Langa KM, Phillips JWR, Weir DR. Cohort Profile: the Health and Retirement Study (HRS). Int J Epidemiol. 2014;43(2):576–85. https://doi.org/10.1093/ije/dyu067.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Barcelo H, Faul J, Crimmins E, Thyagarajan B. A Practical Cryopreservation and Staining Protocol for Immunophenotyping in Population Studies. Curr Protoc Cytom. 2018;84(1):e35.

    PubMed  PubMed Central  Google Scholar 

  24. Hunter-Schlichting D, Lane J, Cole B, Flaten Z, Barcelo H, Ramasubramanian R, et al. Validation of a hybrid approach to standardize immunophenotyping analysis in large population studies: The Health and Retirement Study. Sci Rep. 2020;10(1):8759. https://doi.org/10.1038/s41598-020-65016-x.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Klemera P, Doubal S. A new approach to the concept and computation of biological age. Mech Ageing Dev. 2006;127(3):240–8.

    Article  Google Scholar 

  26. Levine ME. Modeling the rate of senescence: can estimated biological age predict mortality more accurately than chronological age? J Gerontol A Biol Sci Med Sci. 2013;68(6):667–74 (https://pubmed.ncbi.nlm.nih.gov/23213031. 2012/12/03).

    Article  Google Scholar 

  27. Kwon D, Belsky DW. A toolkit for quantification of biological age from blood chemistry and organ function test data: BioAge. GeroScience. 2021;43(6):2795–808. https://doi.org/10.1007/s11357-021-00480-5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Crimmins EM, Faul JD, Thyagarajan B, Weir DR. Venous Blood Collection and Assay Protocol in the 2016 Health and Retirement Study. Ann Arbor: Survey Research Center, Institute for Social Research, University of Michigan; 2017.

    Google Scholar 

  29. Olson NC, Doyle MF, Jenny NS, Huber SA, Psaty BM, Kronmal RA, et al. Decreased Naive and Increased Memory CD4+ T Cells Are Associated with Subclinical Atherosclerosis: The Multi-Ethnic Study of Atherosclerosis. PLoS One. 2013;8(8):e71498. https://doi.org/10.1371/journal.pone.0071498.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Olson NC, Doyle MF, de Boer IH, Huber SA, Jenny NS, Kronmal RA, et al. Associations of Circulating Lymphocyte Subpopulations with Type 2 Diabetes: Cross-Sectional Results from the Multi-Ethnic Study of Atherosclerosis (MESA). PLoS One. 2015;10(10):e0139962. https://doi.org/10.1371/journal.pone.0139962.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Winer S, Chan Y, Paltser G, Truong D, Tsui H, Bahrami J, et al. Normalization of obesity-associated insulin resistance through immunotherapy. Nat Med. 2009;15(8):921–9. https://doi.org/10.1038/nm.2001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Rocha VZ, Folco EJ, Sukhova G, Shimizu K, Gotsman I, Vernon AH, et al. Interferon-gamma, a Th1 cytokine, regulates fat inflammation: a role for adaptive immunity in obesity. Circ Res. 2008;103(5):467–76.

    Article  CAS  Google Scholar 

  33. Lau EYM, Carroll EC, Callender LA, Hood GA, Berryman V, Pattrick M, et al. Type 2 diabetes is associated with the accumulation of senescent T cells. Clin Exp Immunol. 2019;197(2):205–13 (https://pubmed.ncbi.nlm.nih.gov/31251396. 2019/07/08).

    Article  CAS  Google Scholar 

  34. Yang P, Ma J, Yang X, Li W. Peripheral CD4+ naïve/memory ratio is an independent predictor of survival in non-small cell lung cancer. Oncotarget. 2017;8(48):83650–9.

    Article  Google Scholar 

  35. Tuttolomondo A, Pecoraro R, Casuccio A, Di Raimondo D, Buttà C, Clemente G, et al. Peripheral Frequency of CD4+ CD28− Cells in Acute Ischemic Stroke: Relationship With Stroke Subtype and Severity Markers. Medicine (Baltimore). 2015;94(20):e813 (https://journals.lww.com/md-journal/Fulltext/2015/05040/Peripheral_Frequency_of_CD4__CD28__Cells_in_Acute.19.aspx).

    Article  CAS  Google Scholar 

  36. Wikby A, Maxson P, Olsson J, Johansson B, Ferguson FG. Changes in CD8 and CD4 lymphocyte subsets, T cell proliferation responses and non-survival in the very old: the Swedish longitudinal OCTO-immune study. Mech Ageing Dev. 1998;102(2):187–98 (https://www.sciencedirect.com/science/article/pii/S0047637497001516).

    Article  CAS  Google Scholar 

  37. Liu Z, Kuo P-L, Horvath S, Crimmins E, Ferrucci L, Levine M. A new aging measure captures morbidity and mortality risk across diverse subpopulations from NHANES IV: A cohort study. PLOS Med. 2019;15(12):e1002718. https://doi.org/10.1371/journal.pmed.1002718.

    Article  Google Scholar 

  38. Bahour N, Cortez B, Pan H, Shah H, Doria A, Aguayo-Mazzucato C. Diabetes mellitus correlates with increased biological age as indicated by clinical biomarkers. GeroScience. 2022;44(1):415–27. https://doi.org/10.1007/s11357-021-00469-0.

    Article  CAS  PubMed  Google Scholar 

  39. Waziry R, Gras L, Sedaghat S, Tiemeier H, Weverling GJ, Ghanbari M, et al. Quantification of biological age as a determinant of age-related diseases in the Rotterdam Study: a structural equation modeling approach. Eur J Epidemiol. 2019;34(8):793–9.

    Article  CAS  Google Scholar 

  40. Tay RE, Richardson EK, Toh HC. Revisiting the role of CD4+ T cells in cancer immunotherapy—new insights into old paradigms. Cancer Gene Ther. 2021;28(1):5–17. https://doi.org/10.1038/s41417-020-0183-x.

    Article  CAS  PubMed  Google Scholar 

  41. Valzasina B, Piconese S, Guiducci C, Colombo MP. Tumor-Induced Expansion of Regulatory T Cells by Conversion of CD4<sup>+</sup>CD25<sup>−</sup> Lymphocytes Is Thymus and Proliferation Independent. Cancer Res. 2006;66(8):4488–95 (http://cancerres.aacrjournals.org/content/66/8/4488.abstract).

    Article  CAS  Google Scholar 

  42. Liu VC, Wong LY, Jang T, Shah AH, Park I, Yang X, et al. Tumor Evasion of the Immune System by Converting CD4<sup>+</sup>CD25<sup>−</sup> T Cells into CD4<sup>+</sup>CD25<sup>+</sup> T Regulatory Cells: Role of Tumor-Derived TGF-β. J Immunol. 2007;178(5):2883–92 (http://www.jimmunol.org/content/178/5/2883.abstract).

    Article  CAS  Google Scholar 

  43. Larbi A, Fulop T. From, “truly naïve” to “exhausted senescent” T cells: When markers predict functionality. Cytom Part A. 2014;85(1):25–35. https://doi.org/10.1002/cyto.a.22351.

    Article  CAS  Google Scholar 

  44. Zitvogel L, Tesniere A, Kroemer G. Cancer despite immunosurveillance: immunoselection and immunosubversion. Nat Rev Immunol. 2006;6(10):715–27. https://doi.org/10.1038/nri1936.

    Article  CAS  PubMed  Google Scholar 

  45. Sallusto F, Lenig D, Förster R, Lipp M, Lanzavecchia A. Two subsets of memory T lymphocytes with distinct homing potentials and effector functions. Nature [Internet]. 1999;401(6754):708–12. https://doi.org/10.1038/44385.

    Article  CAS  PubMed  Google Scholar 

  46. Klebanoff CA, Gattinoni L, Restifo NP. CD8+ T-cell memory in tumor immunology and immunotherapy. Immunol Rev. 2006;211:214–24.

    Article  CAS  Google Scholar 

  47. White K, Avendaño M, Capistrant BD, Robin Moon J, Liu SY, Maria GM. Self-reported and measured hypertension among older US- and foreign-born adults. J Immigr Minor Heal. 2012;14(4):721–6.

    Article  Google Scholar 

  48. Mullins MA, Kler JS, Eastman MR, Kabeto M, Wallner LP, Kobayashi LC. Validation of Self-reported Cancer Diagnoses Using Medicare Diagnostic Claims in the US Health and Retirement Study, 2000–2016. Cancer Epidemiol Biomarkers Prev. 2022;31(1):287-292. https://doi.org/10.1158/1055-9965.EPI-21-0835.

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Institute on Aging (NIA; R01 AG AG060110). The Health and Retirement Study is supported by NIA (U01 AG009740).

Author information

Authors and Affiliations

Authors

Contributions

R.R. worked on the data analysis and drafted the manuscript; B.T., R.R., J.F. and E.M.C. developed the hypothesis, and provided critical feedback to the manuscript. S.V., H.C.S.M, E.K. and J.N.-Z. provided critical feedback to the analysis strategy and edited the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Bharat Thyagarajan.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Boards at University of Michigan and the University of Minnesota (IRB#: 1601M82462).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

T cell subset definitions measured in the Health and Retirement Study. Figure S1. Pearson correlation heatmap between individual T-cell subsets measured in the Health and Retirement Study. Figure S2. Association of individual T-cell subsets with chronological age. Age is used as the dependent variable in survey linear regression models. The beta estimates are estimated per one standard deviation unit of individual T-cell subsets. The models are adjusted for sex, race/ethnicity, and CMV status. The solid black line along 0 indicates no association. Figure S3. Scatterplot of ARIP measures with chronological age with a LOWESS curve. Figure S4. Association of individual T-cell subsets with biological age. Biological age is used as the dependent variable in survey linear regression models. The beta estimates are estimated per one standard deviation unit of individual T-cell subsets. The models are adjusted for chronological sex, race/ethnicity, and CMV status. The solid black line along 0 indicates no association. Figure S5. Scatterplot of ARIP measures with biological age with a LOWESS curve. Table S2. Beta estimates of the association between the additional ARIP measures with chronological age and biological age. Table S3. Beta estimates of the association between ARIP measures and biological age acceleration measure after adjustment for age, sex, race/ethnicity, and CMV status. Figure S6. Odds ratios and 95% CI of association of multimorbidity levels with individual T-cell subsets per one SD unit increase in ARIP marker. Adjusted for age, sex, race/ethnicity, CMV status, smoking status, and BMI. The solid black line along OR of 1 indicates no association. Table S4. Estimates of association of immunosenescence measures with age-related outcomes. Table S5. Associations between self-rated health and ARIP measures after adjustment for age, sex, race/ethnicity, CMV status, smoking status, and BMI.

Additional file 2.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ramasubramanian, R., Meier, H.C.S., Vivek, S. et al. Evaluation of T-cell aging-related immune phenotypes in the context of biological aging and multimorbidity in the Health and Retirement Study. Immun Ageing 19, 33 (2022). https://doi.org/10.1186/s12979-022-00290-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12979-022-00290-z

Keywords