- Research
- Open access
- Published:
Evaluation of T-cell aging-related immune phenotypes in the context of biological aging and multimorbidity in the Health and Retirement Study
Immunity & Ageing volume 19, Article number: 33 (2022)
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.
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).
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).
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).
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.
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.
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
Gruver AL, Hudson LL, Sempowski GD. Immunosenescence of ageing. J Pathol. 2007;211(2):144–56.
Pawelec G. Age and immunity: What is “immunosenescence”? Exp Gerontol. 2017;2018(105):4–9.
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.
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).
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.
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.
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.
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.
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).
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.
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.
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.
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
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.
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.
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.
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.
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.
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.
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).
Heeringa SG, Connor J. Technical Description of the Health and Retirement Study Sample Design. Ann Arbor: Institute for Social Research, University of Michigan; 1995.
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.
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.
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.
Klemera P, Doubal S. A new approach to the concept and computation of biological age. Mech Ageing Dev. 2006;127(3):240–8.
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).
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.
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.
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.
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.
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.
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.
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).
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.
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).
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).
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.
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.
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.
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.
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).
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).
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.
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.
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.
Klebanoff CA, Gattinoni L, Restifo NP. CD8+ T-cell memory in tumor immunology and immunotherapy. Immunol Rev. 2006;211:214–24.
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.
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.
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
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
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.
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.
About this article
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12979-022-00290-z