Profiling an integrated network of cellular senescence and immune resilience measures in natural aging: a prospective multi-cohort study
========================================================================================================================================

* Natalia Mitin
* Amy Entwistle
* Anne Knecht
* Susan L. Strum
* Allison Ross
* Kirsten Nyrop
* Hyman B. Muss
* Denis Tsygankov
* Joseph M. Raffaele

## Abstract

**Background** Biological aging begins decades before the onset of age-related clinical conditions and is mediated by both cellular senescence and declining adaptive immune function. These processes are functionally related with the rate of senescent cell accumulation dependent upon a balance between induction and immune clearance. We previously showed that biomarkers in these domains can identify patients at-risk of surgery-related adverse events. Here, we describe evidence of clinical relevance in early aging and metabolic phenotypes in a general adult population.

**Methods** We enrolled a total of 482 participants (ages 25-90) into two prospective, cross-sectional healthy aging cohorts. Expression of biomarkers of adaptive immune function and cellular senescence (SapereX) was measured in CD3+ T cells isolated from peripheral blood.

**Findings** We established a network of biomarkers of adaptive immune function that correlate with cellular senescence and associate with early aging phenotypes. SapereX immune components associated with a decrease in CD4+ T cells, an increase in cytotoxic CD8+ T cells, and a loss of CD8+ naïve T cells (Pearson correlation 0.3-0.6). These components also associated with a metric of immune resilience, an ability to withstand antigen challenge and inflammation. In contrast, SapereX components were only weakly associated with GlycanAge (Pearson correlation 0.03-0.15) and commonly used DNA methylation clocks (Pearson correlation 0-0.25). Finally, SapereX biomarkers, in particular p16, were associated with chronic inflammation and metabolic dysregulation.

**Interpretation** Measurement of SapereX biomarkers may capture essential elements of the relationship between cellular senescence and dysregulated adaptive immune function and may provide a benchmark for clinically relevant health decisions.

## Introduction

Aging is a complex process that begins decades before the onset of common age-related conditions such as frailty and overt chronic disease1–3. As such, many aspects of aging remain poorly understood. However, the immune system is a central mediator of organismal aging with well-described early changes beginning in the second or third decade. Therefore, biomarkers of immune system aging are needed to identify early vulnerabilities and allow targeted interventions to both repair immune system function and improve an individual’s aging trajectory.

Immune aging is comprised of immunosenescence, characterized by changes in the proportion and function of numerous immune cell types, as well as replicative or cellular senescence. Cellular senescence is an umbrella term describing a cellular response to a wide variety of stimuli including genomic instability, telomere attrition, epigenetic changes, and impaired proteostasis4. Cells in this state are characterized by stable growth arrest, resistance to apoptosis, and a complex, pro-inflammatory secretome. Cellular senescence promotes aging, and the removal of senescent cells reverses aging pathology5–7. Seminal work by Yousefzadeh and Niedernhofer showed that induction of cellular senescence solely in immune cells can induce cellular senescence in diverse cell and tissue types throughout the body, impairing their function8.

Cells that have undergone cellular senescence accumulate with age. The rate of accumulation depends on both the rate of induction by the stressors described above, as well as the capacity of the immune system, particularly cytotoxic CD8+ T cells, to surveil and clear these altered cells9,10. In biologically young, healthy organisms, cellular senescence is maintained in an adaptive balance. As the organism ages or the immune system is impaired due to other challenges, cellular senescence accumulates and causes dysfunction. While multiple lines of evidence from independent labs suggest a role for the adaptive immune system in the maintenance of cellular senescence, the relationship is complex and is just beginning to be understood.

Previous studies showed little association between cellular senescence and traditional measures of immunosenescence1,11. We took a different approach, interrogating a curated list of biomarkers known to reflect adaptive immune function and compared them to the expression of p16, the most uniquely associated biomarker of cellular senescence, as well as longitudinal changes in p16 (ref12, Justice et al, in print; Sapere Bio, unpublished). Biomarkers that correlated with p16 individually were retained for further analysis as potential biomarkers of immunosurveillance of senescent cells. In the first evidence of clinical relevance, the resulting network of biomarkers was used to identify patients at risk for adverse events associated with cardiac surgery12. Here, we continue to interrogate and characterize this set of biomarkers to establish a broader relationship between cellular senescence and adaptive immune system function. We now describe evidence of clinical relevance in early aging and metabolic phenotypes in a cohort of healthy adults.

## Materials and Methods

### Study design and participants

Two cohorts were enrolled. Community-dwelling adults were enrolled into a natural aging registry study (community cohort). Patients seeking physician-guided health interventions for longevity were enrolled into a long-term registry study (physician-guided cohort). Adults between the ages of 25-85 years old living in North Carolina were recruited into the community cohort. Exclusion criteria were diagnosis of an autoimmune disorder, history or current or chemotherapy, immunotherapy, or radiation therapy, history of solid organ or bone marrow transplants, dialysis, or an active infection (acute or chronic) for which antibiotics and/or antivirals were prescribed within the last 14 days. To ensure representation across age ranges among the 250 participants, we enrolled 30 participants in the age group of 25-34 years, 45 in 35-44, 50 in 45-55, 50 in 55-64, 45 in 65-74, and 30 in 75-85. The clinical study was registered with clinicaltrials.gov ([NCT05123859](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT05123859&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom)). The IRB of the University of North Carolina at Chapel Hill gave ethical approval for this study. Adult patients (25 years and older) seeking care at the Raffaele Medical Clinic were enrolled into a physician-guided cohort. There were no exclusion criteria for medical reasons. An independent IRB (Advarra) gave ethical approval for this study.

### Sample collection and biomarker measurements

Peripheral blood samples were collected at a lab facility (Any Lab Test Now, Durham, NC) or at the Raffaele Medical clinic using Sapere Bio sample collection kits that stabilize whole blood. 7.5ml of blood was used to isolate T cells as previously reported13. T cell pellets were stored at -80°C until further analysis. Gene expression of p14, p16, LAG3, CD28, and CD244 was analyzed by real-time qPCR and normalized to housekeeping genes. Positive and negative controls were included in each run, and Cts over 37 were considered below the limit of detection. Expression of each SapereX gene is reported as log2 of arbitrary units, as is standard for qPCR reporting.

### Statistical analyses

Descriptive statistics were reported as mean (standard deviation) for continuous variables and as frequency (%) for categorical variables. For all analyses, missing data for any reason was not imputed. If a measure was missing, the subject was excluded from summary statistics and analyses. Pairwise associations were quantified using the Pearson correlation coefficient and p-value (two-tailed). One-way ANOVA analyses were used where three or more groups were present. Two-tailed p-values of less than 0.05 were considered statistically significant. Statistical analyses were performed in SAS version 9.4 and JMP 12.2.0 (SAS Institute, Cary, NC).

## Results

### Study cohorts

In this study, a total of 482 participants were enrolled in one of two registries – 250 community-dwelling adults (>25 years) were enrolled in a community cohort and 232 patients seeking physician-guided health interventions for longevity were enrolled in a physician-guided cohort. As shown in Table 1, demographics and clinical characteristics were similar across cohorts with some exceptions. The community cohort had more female participants, higher BMI, and a higher rate of diabetes. The physician-guided cohort had a higher rate of hyperlipidemia and pre-diabetes as well as small statistically (but not clinically) significant differences in other measures such as blood pressure, heart rate, hemoglobin, and hematocrit. Notably, the physician-guided cohort had a higher neutrophil-to-lymphocyte ratio, a marker of systemic chronic inflammation14–16.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/T1)

Table 1. Participant characteristics

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/T2)

Table 2. Association between SapereX biomarkers and early pathology.

### Relationship between biomarkers of cellular senescence and immune function

We have previously identified a network of biomarkers that link cellular senescence, immune system function, predisposition to surgery-related adverse events12. This network included gene expression of p16, p14, LAG3, CD28, and CD244 measured in peripheral blood T lymphocytes. We refer to this network of biomarkers as SapereX. p16 is the biomarker most uniquely associated with cellular senescence and has been used as a marker of established senescent cells in studies linking cellular senescence, functional decline, and aging-related diseases17,18. A decrease in CD2819 and an increase in LAG3 and CD244 have all been implicated in the loss of immune function, particularly in CD8+ T cells. p14, a transcript related to p1620, is a regulator of cell cycle arrest and apoptosis through p53/p21. The pairwise association between biomarkers within SapereX network is shown in Figure 1. Overall, there was a large degree of association between all the components, with CD28 negatively correlating with the other markers. Interestingly, while p16 and p14 were associated with chronological age, other markers were weakly associated and had highly heterogeneous levels of expression in individuals within every age group. Therefore, there is a significant association between levels of a biomarker of cellular senescence (p16) and biomarkers of adaptive immune function.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/28/2023.08.25.23294589/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/F1)

Figure 1. 
Relationship between biomarkers within SapereX network. Scatterplot correlation matrix of expression levels (log2) of p16, p14, LAG3, CD28, CD244 as well as chronological age. The color of each correlation circle represents the correlation between each pair of variables on a scale from red (+1) to blue (-1). The size of each circle represents the significance test between the variables. A larger circle indicates a more significant relationship, and the Pearson correlation coefficient is shown as a number and a line of linear fit on the corresponding scatterplot. The histograms (diagonal across the matrix) show the distribution of each biomarker in the entire cohort.

### SapereX biomarker network and age-associated immune decline

We also examined the relationship between SapereX and conventional biomarkers of immune function. Traditionally, a decline in immune system function is associated with a decrease in lymphocytes with the corresponding increase in neutrophils, a decrease in CD4+ and an increase in CD8+ cytotoxic T cells, and more specifically an increase of CD8+CD28-T cells that have impaired proliferation and a decrease in CD8+CD95+ naïve T cells. CD4+, CD8+ T cell subset count data via flow cytometry was available in 168 participants in a physician-guided cohort and is shown in Figure 2. As expected, LAG3, CD244, and interestingly, p14 were associated with a decrease in CD4+ and an increase in CD8+ T cells, a decrease in CD8+CD95+ T cells and an increase in CD8+CD28-T cells. CD28 gene expression as measured within SapereX (presence of CD28) was inversely associated with all of the above. Expression levels of p16 were also positively associated with biomarkers of immune decline (a decrease in CD8+CD95+ T cells and an increase in CD8+CD28-T cells) (Figure 2) and also a decline in lymphocytes and an increase in neutrophils.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/28/2023.08.25.23294589/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/F2)

Figure 2. 
Relationship between SapereX biomarker network and established molecular metrics of age-associated immune decline. Percent of lymphocytes or neutrophils (from CBC labs) as well as CD4+, CD8+, CD8+CD28- and CD8+CD95+ T cells (measured by flow cytometry) were compared to expression levels of SapereX biomarkers. Scatterplot correlation matrix of expression levels (log2) of p16, p14, LAG3, CD28, CD244 as well as chronological age. The color of each correlation circle represents the correlation between each pair of variables on a scale from red (+1) to blue (-1). The size of each circle represents the significance test between the variables. A larger circle indicates a more significant relationship, and the Pearson correlation coefficient is shown as a number and a line of linear fit on the corresponding scatterplot. The histograms (diagonal across the matrix) show the distribution of each biomarker in the entire cohort.

### SapereX and other molecular aging biomarkers

Next, we examined the relationship between SapereX and other biomarkers of aging such as telomere length, intrinsic and extrinsic methylation clocks and GlycanAge. GlycanAge was particularly interesting given its association with immune aging and chronic inflammation. Telomere length was determined as the difference between granulocyte and lymphocyte telomere length measurements to account for individual differences in telomere length (ref). Intrinsic epigenetic aging clock (Fig.3 “Intrinsic”) represents changes in DNA methylation that best fit chronological age. Extrinsic epigenetic clocks (“Extrinsic”, Pace of Aging (PoA) and PhenoAge) are postulated to provide additional information on the aging process. Consistent with our previously published data, there was a weak association between p16 and epigenetic clocks, except for PoA (ref). p16 expression was also weakly positively associated with telomere shortening. There was no association between p16 and GlycanAge. p14 expression showed a similar pattern to p16. There was even less association between the remaining components of SapereX network and other aging biomarkers, except for CD244 and telomere shortening (r=0.34). While, GlycanAge showed no association with any component of the SapereX network, it was associated with the extrinsic epigenetic age clock (TruDiagnostic) and the PhenoAge biological age calculation. Overall, our data suggest that SapereX biomarkers capture different aspects of the aging process than epigenetic clocks, GlycanAge, and to some extent, telomere shortening.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/28/2023.08.25.23294589/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/F3)

Figure 3. 
Relationship between SapereX biomarker network and molecular biomarkers of aging. Telomere length, GlycanAge or epigenetic clocks (Pace of aging, intrinsic, extrinsic, PhenoAge) were compared to expression levels of SapereX biomarkers. Scatterplot correlation matrix of expression levels (log2) of p16, p14, LAG3, CD28, CD244 as well as chronological age. The color of each correlation circle represents the correlation between each pair of variables on a scale from red (+1) to blue (-1). The size of each circle represents the significance test between the variables. A larger circle indicates a more significant relationship, and the Pearson correlation coefficient is shown as a number and a line of linear fit on the corresponding scatterplot. The histograms (diagonal across the matrix) show the distribution of each biomarker in the entire cohort.

### SapereX is associated with immune resilience

A recent study of ∼ 48,500 participants developed a quantitative measure of immune resilience linked to longevity and defined as the ability of the immune system to withstand antigen challenges and regulate inflammation21. That study weighted CD4+ and CD8+ absolute T cell counts to derive progressively worsening immune grades, with grade I defined as a primordial state that is optimal for immune health. We utilized the same algorithm to categorize immune resilience grades for 168 participants with known CD4+ and CD8+ T cell counts (same subset as Figure 2). Then, we examined the distribution of SapereX biomarkers across grades of immune resilience. LAG3, p14 and CD244 expression progressively increased with the loss of immune resilience (Figure 4). CD28 expression progressively decreased. On the other hand, p16 and chronological age were not independently associated with changes in immune resilience. In support of that idea, PCA analysis (Figure 5) grouped p16 and chronological age in a separate cluster from the rest of the SapereX network.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/28/2023.08.25.23294589/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/F4)

Figure 4. 
A subset of the SapereX biomarker network is associated with immune resilience. Distribution of SapereX biomarkers, chronological age, or GlycanAge across categories of immune resilience. Immune grades (IHG) were determined based on CD4+ and CD8+ T cell counts and represent grade I, IIa, IIb, and IIc21. Participants in Grade I have the most resilience towards antigen-specific challenges and aging in general. p value for one-way ANOVA analysis among groups is shown on each plot.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/12/28/2023.08.25.23294589/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2023/12/28/2023.08.25.23294589/F5)

Figure 5. 
Principal components analysis of SapereX biomarkers and chronological age identified two clusters. P16 and age in cluster #1 and biomarkers associated with adaptive immune system decline and p14 in cluster #2.

Surprisingly, GlycanAge and an extrinsic methylation clock of immune aging (TruDiagnostic) showed no association with immune resilience (Figure 4).

### SapereX is associated with early clinical indicators of decline

We have established that the SapereX network is associated with a number of age-related changes in the adaptive immune system including a new measure of immune resilience. As previously mentioned, a decline in immune system function is linked to the earliest stages of organismal decline, well before the presentation of chronic diseases. We have examined the association between SapereX and clinical markers of physiologic decline. Given the lab results available in these cohorts, we chose biomarkers of metabolic health (insulin resistance, fasting glucose, HbA1C), cardiovascular health (total cholesterol), and neutrophil-to-lymphocyte ratio (NLR), a marker of low-grade systemic inflammation. p16, and not other SapereX components was associated with metabolic health (HbA1C and insulin resistance) and total cholesterol. p16, p14, LAG3, and to a lesser extent CD28 were associated with NLR. None of the SapereX biomarkers were associated with fasting glucose. Association of p16 with insulin resistance is consistent with data in mice where depletion of senescent cells improved insulin sensitivity22.

## Discussion

Aging research has historically focused on the very old and on diseases of old age. Recently, attention has shifted to the idea that biological aging begins early and, as such, investigations in subjects across a wide age range and outside the context of reactive disease care are foundational to geroscience. Here, we describe a network of biomarkers that capture both cellular senescence and T cell function and have shown their association with early indicators of physiologic decline. Given a number of recent studies demonstrating that T cells, specifically CD8+ T cells, are critical for immune surveillance of senescent cells (ref), our work provides an early description of markers that can define the relationship between these domains in humans.

Biomarkers in the SapereX network were selected based on previous studies (ref12, Justice et al, in print; Sapere Bio, unpublished). Selection criteria included an independent association with p16 and each marker’s capacity to capture distinct as well as potentially overlapping biological functions of the adaptive immune system. Multiple markers of presumably similar function (e.g., T cell exhaustion) were screened for association with p16 as well as the ability to change with longitudinal changes in p16. LAG3 and CD244 were both described as markers of T cell exhaustion and LAG3 is the newest target of immuno-oncology therapies.

CD244 appears to have a distinct role in T cell function, accumulating with age in CD8+ T cells where it is associated with increased proliferation and apoptosis, and terminal differentiation. In addition, CD244 has been shown to be an inhibitor of autophagy by directly interacting with and suppressing activity of the Beclin-1/Vps34 complex23. Autophagy is essential for the function of T cells including effector CD8+ T cell survival and memory formation and its activity declines with age24,25. While there is a possibility that the immune biomarkers we chose do not fully capture the ability of T cells to surveil senescent cells, our finding that immune biomarkers (LAG3, CD244, CD28) as well as p14 are strongly associated with immune resilience, ie, the ability of the organism to withstand immune challenge and resist inflammation, reinforce the importance of the selected markers. Our data are consistent with a study by Burd et al in healthy adults where p16 expression in T cells was associated with CD28 and CD244 and not other markers of T cell function such as KLRG1, CD57, CD45RA, and the PD-1 pathway11.

While SapereX immune markers were independently associated with p16, their relationship is likely to be highly interconnected. For example, in patients with cardiovascular disease, a second-degree interaction between p16 and CD244 was required to predict surgery-associated adverse events12. Our study design does not allow us to infer directionality in the relationship between cellular senescence and immune function. While dysfunction in the immune system can lead to accumulation of senescent cells, there is also evidence that accumulation of senescent cells can damage the immune system by inducing T cell exhaustion and reducing T cell proliferation26,27. Interestingly, of all the SapereX biomarkers, p14 is most strongly associated with p16, yet clusters with markers of T cell function in the PCA analysis and is associated with immune resilience. It is possible that p14 modulates both domains of the SapereX network-cellular senescence and immune function. p14 may affect cellular senescence through regulation of apoptosis via the p53/p21 axis28. Consistent with its putative role in cellular senescence, continuous p14 expression was required for cellular commitment to senescence in the epidermis29. Likewise, p14 may regulate immune function. Studies have shown that induction of p14 led to mitochondrial dysfunction and T cell exhaustion in tumor-infiltrating lymphocytes30.

Biomarkers in the SapereX network were mostly discordant with other markers of biological aging, including methylation clocks, telomere shortening, and GlycanAge. While global epigenetic changes are associated with the development and maintenance of the senescent phenotype31–33, the most commonly used methylation clocks do not correlate with cellular senescence (Sedrak 2023, in print, and Fig 3) or SapereX immune biomarkers. Telomere shortening was modestly associated with CD244 (r=0.34), possibly due to the link between CD244 and apoptosis. Not surprisingly, telomere shortening was also associated with p16 (r=0.28), but not other markers. Interestingly, neither GlycanAge nor the intrinsic epigenetic clock (TruDiagnostic), putative biomarkers of immune function and chronic inflammation, was associated with immune resilience (Figure 3). Likewise, neither GlycanAge nor the intrinsic epigenetic clock was associated with components of the SapereX network.

With increasing attention on interventions that may extend human healthspan and ongoing clinical trials to alter organismal cellular senescence and other aging pathways, it is increasingly important to understand the dynamics of these processes in humans. A current hypothesis in the field is that a decline in immune function may precede accumulation of senescent cells (ref). In some patients, addressing immune dysfunction may provide an indirect way of targeting cellular senescence and improving aging trajectory. Alternatively, understanding the interplay between the immune system and cellular senescence may improve patient selection and avoid harm in clinical trials. For example, senolytic therapy may be inappropriate for patients with already low cellular senescence and normal immune biomarkers, while patients with low cellular senescence but dysregulated immune function may benefit from interventions targeting the immune system. And vice versa, patients with high cellular senescence but normal immune biomarkers may be better suited to a senolytic approach absent an immune intervention.

Our study has some limitations including a relatively small sample size, particularly for clinical variables collected for a subset of participants, e.g insulin resistance. Another limitation is the cross-sectional design. Longitudinal data are needed to validate the presence of a temporal relationship between biomarkers of immune function and cellular senescence and that work is ongoing. In addition, CD11b+ dendritic cells appear to be essential mediators of surveillance of senescent cells by CD8+ T cells34. In this study, we have not captured measures of dendritic cell activation or function. However, we have previously demonstrated that plasma levels of suPAR, a protein secreted mostly by the immune cells35, contributes to the predictive ability of SapereX to identify clinical vulnerabilities mediated by uncontrolled inflammation in the setting of cardiac surgery.

There is a time-honored adage in medicine - ‘what gets measured, gets managed’. Aging, in its foundational stages, is rarely managed in part because clinically relevant molecular measures are lacking36. This despite growing recognition that progressive and dramatic loss of physiological resilience precedes frailty and chronic disease by many decades. The ability to measure and map changes in immune aging through cellular senescence and T cell function creates the potential for improved safety, efficacy, and personalization of diverse interventions aimed at improving health span, longevity, and functional performance throughout life.

## DATA AVAILABILITY

Due to the informed consent and data privacy policies, the clinical data are not publicly available, i.e., accessible for anyone, for any purpose without a review by the Central IRB on a project-by-project basis.

## FUNDING

This work was funded in part by grant from the NIH/NIA R21 AG070356 to DT. The funding sources had no involvement in the study design or in the collection, analysis and interpretation of data. The funding sources also had no involvement in the writing of this report or in the decision to submit the report for publication.

## DISCLOSURES

NM, AE, AK, SLS are employees of Sapere Bio, hold equity stake and inventors on patents

## Footnotes

*   Figure files have been reuploaded to correct image quality and sizing issues

*   Received August 25, 2023.
*   Revision received December 28, 2023.
*   Accepted December 28, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## REFERENCES

1.  Liu Y, Sanoff HK, Cho H, et al. Expression of p16(INK4a) in peripheral blood T-cells is a biomarker of human aging. Aging Cell 2009; 8: 439–48.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1474-9726.2009.00489.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19485966&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268213300009&link_type=ISI) 

2.  Elliott ML, Caspi A, Houts RM, et al. Disparities in the pace of biological aging among midlife adults of the same chronological age have implications for future frailty risk and policy. Nat Aging 2021; 1: 295–308.
    
    

3.  Belsky DW, Caspi A, Houts R, et al. Quantification of biological aging in young adults. Proceedings of the National Academy of Sciences of the United States of America 2015; 112: E4104–10.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEyLzMwL0U0MTA0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTIvMjgvMjAyMy4wOC4yNS4yMzI5NDU4OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

4.  He S, Sharpless NE. Senescence in Health and Disease. Cell 2017; 169: 1000–11.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2017.05.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28575665&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 

5.  Baker DJ, Wijshake T, Tchkonia T, et al. Clearance of p16Ink4a-positive senescent cells delays ageing-associated disorders. Nature 2011; 479: 232–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10600&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22048312&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298030800043&link_type=ISI) 

6.  Baker DJ, Childs BG, Durik M, et al. Naturally occurring p16INK4A -positive cells shorten healthy lifespan. Nature 2016; 530: 184–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature16932&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26840489&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 

7.  Xu M, Pirtskhalava T, Farr JN, et al. Senolytics improve physical function and increase lifespan in old age. Nature Medicine 2018; 24: 1246–56.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-018-0092-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29988130&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 

8.  Yousefzadeh MJ, Flores RR, Zhu Y, et al. An aged immune system drives senescence and ageing of solid organs. Nature 2021; 594: 100–5.
    
    

9.  Marin I, Serrano M, Pietrocola F. Recent insights into the crosstalk between senescent cells and CD8 T lymphocytes. Aging 2023; 9.
    
    

10. Kang T-W, Yevsa T, Woller N, et al. Senescence surveillance of pre-malignant hepatocytes limits liver cancer development. Nature 2011; 479: 547–51.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10599&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22080947&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000297285600055&link_type=ISI) 

11. Burd CE, Peng J, Laskowski BF, et al. Association of Epigenetic Age and p16INK4a With Markers of T-Cell Composition in a Healthy Cohort. The Journals of Gerontology: Series A 2020; **XX**: 1–5.
    
    

12. Entwistle A, Walker S, Knecht A, Strum SL, Mitin N. Signature of pre-operative biomarkers of cellular senescence to predict risk of cardiac and kidney adverse events after cardiac surgery. medRxiv 2023; published online April 3. DOI:10.1101/2023.04.03.23288081.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMy4wNC4wMy4yMzI4ODA4MXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMTIvMjgvMjAyMy4wOC4yNS4yMzI5NDU4OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

13. Mitin N, Nyrop KA, Strum SL, et al. A biomarker of aging, p16, predicts peripheral neuropathy in women receiving adjuvant taxanes for breast cancer. npj Breast Cancer 2022; 8: 103.
    
    

14. Zahorec R. Neutrophil-to-lymphocyte ratio, past, present and future perspectives. BLL 2021; 122: 474–88.
    
    

15. Nishijima TF, Deal AM, Lund JL, Nyrop KA, Muss HB, Sanoff HK. Inflammatory markers and overall survival in older adults with cancer. Journal of Geriatric Oncology 2019; 10: 279–84.
    
    

16. Fest J, Ruiter R, Ikram MA, Voortman T, Van Eijck CHJ, Stricker BH. Reference values for white blood-cell-based inflammatory markers in the Rotterdam Study: a population-based prospective cohort study. Sci Rep 2018; 8: 10566.
    
    

17. Liu Y, Sanoff HK, Cho H, et al. Expression of p16 INK4a in peripheral blood T-cells is a biomarker of human aging. Aging Cell 2009; 8: 439–48.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1474-9726.2009.00489.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19485966&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268213300009&link_type=ISI) 

18. Kim WY, Sharpless NE. The regulation of INK4/ARF in cancer and aging. Cell 2006; 127: 265–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2006.10.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17055429&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241827100015&link_type=ISI) 

19. Chen X, Liu Q, Xiang AP. CD8+CD28-T cells: not only age-related cells but a subset of regulatory T cells. Cell Mol Immunol 2018; 15: 734–6.
    
    

20. Ozenne P, Eymin B, Brambilla E, Gazzeri S. The ARF tumor suppressor: Structure, functions and status in cancer. International Journal of Cancer 2010; 127: 2239–47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.25511&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20549699&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283563900001&link_type=ISI) 

21. Ahuja SK, Manoharan MS, Lee GC, et al. Immune resilience despite inflammatory stress promotes longevity and favorable health outcomes including resistance to infection. Nat Commun 2023; 14: 3286.
    
    

22. Xu M, Palmer AK, Ding H, et al. Targeting senescent cells enhances adipogenesis and metabolic function in old age. eLife 2015; 4: e12997.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.12997&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26687007&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 

23. Chaudhary A, Leite M, Kulasekara BR, et al. Human Diversity in a Cell Surface Receptor that Inhibits Autophagy. Current Biology 2016; 26: 1791–801.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cub.2016.05.003&link_type=DOI) 

24. Xu X, Araki K, Li S, et al. Autophagy is essential for effector CD8+ T cell survival and memory formation. Nat Immunol 2014; 15: 1152–61.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ni.3025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25362489&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 

25. Macian F. Autophagy in T Cell Function and Aging. Front Cell Dev Biol 2019; 7: 213.
    
    

26. Janelle V, Neault M, Lebel MÈ, et al. p16INK4a Regulates Cellular Senescence in PD-1-Expressing Human T Cells. Frontiers in Immunology 2021; 12: 1–17.
    
    

27. Palacio L, Goyer M, Maggiorani D, et al. Restored immune cell functions upon clearance of senescence in the irradiated splenic environment. Aging Cell 2019; 18. DOI:10.1111/acel.12971.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/acel.12971&link_type=DOI) 

28. Yosef R, Pilpel N, Tokarsky-Amiel R, et al. Directed elimination of senescent cells by inhibition of BCL-W and BCL-XL. Nat Commun 2016; 7: 11190.
    
    

29. Tokarsky-Amiel R, Azazmeh N, Helman A, et al. Dynamics of Senescent Cell Formation and Retention Revealed by p14ARF Induction in the Epidermis. Cancer Research 2013; 73: 2829–39.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiY2FucmVzIjtzOjU6InJlc2lkIjtzOjk6IjczLzkvMjgyOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzEyLzI4LzIwMjMuMDguMjUuMjMyOTQ1ODkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

30. Koss B, Shields BD, Taylor EM, et al. Epigenetic Control of *Cdkn2a.Arf* Protects Tumor-Infiltrating Lymphocytes from Metabolic Exhaustion. Cancer Research 2020; 80: 4707–19.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1158/1538-7445.AM2020-4707&link_type=DOI) 

31. Cruickshanks HA, McBryan T, Nelson DM, et al. Senescent cells harbour features of the cancer epigenome. Nat Cell Biol 2013; 15: 1495–506.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncb2879&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24270890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F12%2F28%2F2023.08.25.23294589.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000327944200014&link_type=ISI) 

32. Lowe R, Overhoff MG, Ramagopalan SV, et al. The senescent methylome and its relationship with cancer, ageing and germline genetic variation in humans. Genome Biol 2015; 16: 194.
    
    

33. Xie W, Kagiampakis I, Pan L, et al. DNA Methylation Patterns Separate Senescence from Transformation Potential and Indicate Cancer Risk. Cancer Cell 2018; 33: 309–321.e5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ccell.2018.01.008&link_type=DOI) 

34. Marin I, Boix O, Garcia-Garijo A, et al. Cellular Senescence Is Immunogenic and Promotes Antitumor Immunity. Cancer Discovery 2023; 13: 410–31.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1158/2159-8290.CD-22-0523&link_type=DOI) 

35. Hindy G, Tyrrell DJ, Vasbinder A, et al. Increased soluble urokinase plasminogen activator levels modulate monocyte function to promote atherosclerosis. Journal of Clinical Investigation 2022; 132: e158788.
    
    

36. Ferrucci L, Gonzalez-Freire M, Fabbri E, et al. Measuring biological aging in humans: A quest. Aging Cell 2019; : 1–21.