Skip to main content

BRIEF RESEARCH REPORT article

Front. Aging, 30 May 2024
Sec. Genetics, Genomics and Epigenomics of Aging

Naïve T-cell decline is a significant contributor to expression changes in ageing blood

  • Systematic Medicine, Melbourne, VIC, Australia

No clear consensus has emerged from the literature on the gene expression changes that occur in human whole blood with age. In this study we compared whole blood ageing genes from the published literature with data on gene specificity for leukocyte subtypes. Surprisingly we found that highly ranked ageing genes were predominantly expressed by naïve T cells, with limited expression from more common cell types. Highly ranked ageing genes were also more likely to have decreased expression with age. Taken together, it is plausible that much of the observed gene expression changes in whole blood is reflecting the decline in abundance of naïve T cells known to occur with age, rather than changes in transcription rates in common cell types. Correct attribution of the gene expression changes that occur with age is essential for understanding the underlying mechanisms.

1 Introduction

Ageing is a deleterious process that is inevitable for any multicellular organism that survives for enough time. In humans it is associated with an increased burden of disease (Lopez-Otin et al., 2013) and is a considerable challenge for health systems as lifespans increase (Prince et al., 2015).

Several studies have investigated differences in gene expression between young and aged individuals to elucidate the mechanisms of ageing. Whole blood is commonly used for human gene expression studies as it is one of the easiest sample types to obtain and the RNA profile can be rapidly stabilized (Asare et al., 2008). However, there is no clear consensus on the genes or pathways that are differentially expressed with age in whole blood (Supplementary Table S1).

Some studies have suggested that changes in the relative proportion of blood cell types has a considerable influence on the transcriptional changes observed in blood with age. Nakamura et al. (2012) reported that of the 16 age related genes identified, most were strongly associated with lymphocyte lineages. Limitations of this study include the small number of ageing genes identified, and the use of in vitro cell culture data to determine gene specificity. Pellegrino-Coppola et al. (2021) proposed a regression model to correct for changes in cell composition and found this reduced the number of differentially expressed genes. Jonkman et al. (2022) identified a decrease in expression of genes associated with naïve T cells and an increase in expression of genes associated with activated T cells.

In this study we sought to characterize the relative contribution of different leukocyte subtypes to the expression of age associated genes in whole blood. Age associated genes were identified from a review of published literature, and gene specificity was determined using ex vivo expression data.

2 Methods

2.1 Selecting age associated genes

Age associated genes used in this analysis were derived from Peters et al. (2015) (Peters). Several published studies were considered (Supplementary Table S1) however, only Peters meet the following eligibility criteria:

Studies collected peripheral whole blood from humans, and stabilized the samples using Tempus, PAXGene or equivalent technology. The rationale for this was to eliminate potential confounding factors that could change gene expression signatures after sample collection. Eligible studies had at least 500 unique donors that spanned an age range of at least 25–70 years old and were reasonably representative of the general population. To be included, studies also had to quantify at least half of all protein coding genes using either gene expression arrays or RNA-Seq. Group assignment for differential expression needed to be based on chronological age, with a list of differentially expressed genes publicly available or obtainable from the authors on request. For cases where the same data was analyzed in multiple studies, only the larger study was included.

Peters determined differentially expressed genes by conducting a meta-analysis of 13 independent cohort studies (Supplementary Table S2), containing a total of 14,983 unique donors. For each of the cohorts, differentially expressed genes were determined using linear regression analysis, with potentially confounding variables modelled as random variables. For 6 of the 13 studies, counts of granulocytes, lymphocytes and monocytes were modelled as random variables. Despite this we believe the Peters differentially expressed genes are suitable for this analysis as only three broad leukocyte categories were modelled, and only for a minority of cohorts within the meta-analysis. Additionally, up to 10 random variables were included in each model, reducing the degree to which any one covariate could influence the model (Zhang, 2014).

2.2 Gene sets

Peters identified 1,497 age associated genes, and gave each gene a ranking reflecting the strength of association with donor age. Because of the high power achieved from the large sample size, many of the lower ranked genes were associated with small effect sizes and may be less biologically relevant. Due to this consideration, we focused on two gene sets for our analysis.

The first gene set included the 20 most highly ranked ageing genes (Supplementary Table S3), which corresponded with a Z-score of at least half the highest ranked gene. The second gene set included all genes reported by Peters to be differentially expressed with age. 1,459 out of 1,497 (97.5%) could be mapped to proteins in the Human Blood Atlas, and were used in this analysis.

2.3 Attributing gene expression to leukocyte subtypes

The specificity of genes to leukocyte subtypes was determined using data from the Human Blood Atlas (Uhlen et al., 2019). The Human Blood Atlas is an open-access database containing genome wide single cell expression data for protein coding genes for 18 leukocyte subtypes. The markers used by the Human Blood Atlas to define the cell types are summarized in Supplementary Table S4.

The Human Blood Atlas used the same blood samples, sample preparation protocol and sequencing pipeline for all 18 leukocyte subtypes. In addition, a normalization strategy was employed with the specific objective of facilitating comparison of expression values between cell types.

Using the normalized expression data from the Human Blood Atlas, the proportion of expression (x) of each protein coding gene attributable to each cell type was estimated using the following formula:

xi,j=nTPMi,jΣnTPMi

Where (i) is the protein coding gene, (j) is the leukocyte subtype and (nTPM) is the normalized transcripts per million reported by the Human Blood Atlas.

2.4 Statistical analysis

The median proportion of expression attributable to a cell type for a gene set was compared with the median for all protein coding genes using the one-sample Wilcoxon test (α = 0.05). When performing multiple comparisons, p-values were adjusted for multiple testing using the Benjamini & Hochberg approach (Benjamini and Hochberg, 1995).

The proportion of differentially expressed genes reported to decrease or increase with age was compared with a theoretical distribution of 0.5 using the two tailed binomial test.

All data processing and statistical analysis was performed using R (v4.3.1) (R Core Team, 2023) with the packages tidyverse (v2.0.0) (Wickham et al., 2019), biomaRt (v2.58.0) (Durinck et al., 2009), cumstats (v1.0) (Erdely and Castillo, 2017) and ggpubr (v0.6.0) (Kassambara, 2003). Code used for this analysis can be accessed at https://github.com/systematicmedicine/Naive-cell-publication/.

3 Results

3.1 Highly ranked ageing genes predominantly expressed by naïve T cells

Highly ranked ageing genes in whole blood were predominantly expressed by CD4+ and CD8+ naïve T cells (Figure 1). This is surprising as naïve T cells are a relatively rare sub population within whole blood. 87% of the expression of CD248 (Peters rank 1) and 91% of the expression of LRNN3 (Peters rank 2) was attributable to naïve T cells. For the top 20 ranked genes, the median expression attributable to naïve T cells was 35%, 3.8 fold higher than would be expected from a random selection of 20 genes (p = 8.1 × 10−12). The median proportion of gene expression attributable to naïve T-cells reduced as lower ranked genes were included (Figure 2). However, this remained higher for the entire set of 1,459 ageing genes, compared to the median for all protein coding genes (p = 2.9 × 10−24).

Figure 1
www.frontiersin.org

Figure 1. Percentage of gene expression attributable to each of the 18 leukocyte subtypes, for the 20 highest ranked ageing genes reported in Peters. Genes listed in rank order (highest at top). Naïve CD4+ and CD8+ T-cells account for considerably more expression of highly ranked ageing genes than would be expected by chance. Statistical significance was assessed with one-sample Wilcoxon tests. The stars indicate statistical significance: ***p ≤ 0.001, **p ≤ 0.01, *p ≤ 0.05. Dashed red line corresponds to median for all protein coding genes. Direction of expression change refers to direction of gene expression change reported by Peters (negative denotes decreased expression with age).

Figure 2
www.frontiersin.org

Figure 2. Median percentage of gene expression attributable to naïve T cells (CD4+ and CD8+) for ageing genes reported by Peters. Statistical significance was assessed for 1,459 ageing genes with a one-sample Wilcoxon test (p = 2.9 × 10−24). Red dashed line corresponds with median for all protein coding genes.

3.2 Highly ranked ageing genes negatively associated with age

Of the top 20 ranked ageing genes reported by Peters, 17 (85%) had decreased expression with age. Compared with a theoretical distribution of 50% (half of differentially expressed genes increase with age, half decrease) this is unlikely to occur by chance (p = 2.6 × 10−3). For all 1,497 ageing genes reported by Peters, 60% had decreased expression with age which is also unlikely to occur by chance (p = 1.7 × 10−14). The abundance of naïve T cells is widely accepted to decline with age (Britanova et al., 2014; Li et al., 2019; den Braber et al., 2012; van der Geest et al., 2015; Nasi et al., 2006), potentially explaining this observation.

3.3 Highly ranked ageing genes are expressed less than expected in common leukocyte subtypes

Several common leukocyte subtypes were weakly associated with the expression of highly ranked ageing genes, especially myeloid lineages (Figure 1). Expression attributable to basophils, monocytes (classical and intermediate), eosinophils, neutrophils and myeloid dendritic cells were all significantly lower than expected for both gene sets (Table 1). If changes in transcription rates of these common cell types was a major contributor to age related expression changes, we would expect genes they express to feature more prominently in the highly ranked ageing genes.

Table 1
www.frontiersin.org

Table 1. Median proportion of gene expression attributable to each cell type, for ageing gene sets and all protein coding genes.

3.4 Ageing genes that increase in expression with age associated with several T-cell lineages

The majority of ageing genes used in this study decrease in expression with age. When restricting the analysis to the subset of ageing genes that increase expression with age, a significant association was found with several T-cell lineages. For the 20 highest ranked genes that increase expression with age (Supplementary Table S5), expression attributable to GdT-cells, CD8 T cells (memory and naïve) and natural killer cells were significantly higher than expected (Supplementary Figure S1).

4 Discussion

A common assumption is that differential gene expression is primarily driven by changes in cellular transcription. While this is often the case, in heterogenous tissues such as whole blood it can also be driven by changes in the relative proportion of cell types.

This study found that the genes with the strongest association with age were predominately expressed by naïve T-cells, and that most of these age associated genes decreased in expression with age. Given that naïve T-cells are known to decline in abundance with age (Britanova et al., 2014; Li et al., 2019; den Braber et al., 2012; van der Geest et al., 2015; Nasi et al., 2006), we propose that the largest gene expression changes seen in ageing blood may reflect the reduction in naïve-T cells rather than a change in transcription profiles of common cell types.

This has important implications for fundamental research, as incorrect attribution of observed gene expression changes could lead to invalid conclusions being drawn about the underlying mechanisms of ageing. This potential confounding factor of cell composition may also apply to other tissue types and be more difficult to identify if gene expression data from subpopulations is unavailable or less robust. Other age correlated measures, such as DNA methylation, chromatin accessibility and protein abundance may also be confounded by age related changes in tissue composition. For example, naïve T cells have a lower epigenetic age than other blood cell types (Jonkman et al., 2022; Tomusiak et al., 2023).

The findings of this study also have important implications for translational research, especially transcriptomic age prediction models (clocks). Transcriptomic clocks trained on bulk whole blood gene expression data, with no regard for changes in cell composition, may be little more than predictors of naïve T cell decline. Such models would be poor surrogates of biological age and an unsuitable tool for drug discovery. For translational applications, it may be best to use measures that have a clear mechanistic link to the phenotype being targeted (Vincent et al., 2015), rather than black box predictors.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.proteinatlas.org/about/download. Code used for this analysis can be accessed at https://github.com/systematicmedicine/Naive-cell-publication/.

Ethics statement

Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.

Author contributions

CF: Conceptualization, Data curation, Formal Analysis, Methodology, Writing–original draft, Writing–review and editing. BO: Conceptualization, Writing–review and editing, Methodology.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was funded by Systematic Medicine Pty Ltd. The authors had complete autonomy over study design, analysis, interpretation of data, the writing of this article and decision to publish.

Acknowledgments

We would like to acknowledge Scott Needham and James Phie for their feedback on the draft manuscript.

Conflict of interest

Authors CF and BO were employed by Systematic Medicine.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fragi.2024.1389789/full#supplementary-material

References

Asare, A. L., Kolchinsky, S. A., Gao, Z., Wang, R., Raddassi, K., Bourcier, K., et al. (2008). Differential gene expression profiles are dependent upon method of peripheral blood collection and RNA isolation. BMC Genomics 9, 474. doi:10.1186/1471-2164-9-474

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57 (1), 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Britanova, O. V., Putintseva, E. V., Shugay, M., Merzlyak, E. M., Turchaninova, M. A., Staroverov, D. B., et al. (2014). Age-related decrease in TCR repertoire diversity measured with deep and normalized sequence profiling. J. Immunol. 192 (6), 2689–2698. doi:10.4049/jimmunol.1302064

PubMed Abstract | CrossRef Full Text | Google Scholar

den Braber, I., Mugwagwa, T., Vrisekoop, N., Westera, L., Mogling, R., de Boer, A. B., et al. (2012). Maintenance of peripheral naive T cells is sustained by thymus output in mice but not humans. Immunity 36 (2), 288–297. doi:10.1016/j.immuni.2012.02.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Durinck, S., Spellman, P. T., Birney, E., and Huber, W. (2009). Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat. Protoc. 4 (8), 1184–1191. doi:10.1038/nprot.2009.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Erdely, A., and Castillo, I. (2017) Cumulative descriptive statistics.

Google Scholar

Jonkman, T. H., Dekkers, K. F., Slieker, R. C., Grant, C. D., Ikram, M. A., van Greevenbroek, M. M. J., et al. (2022). Functional genomics analysis identifies T and NK cell activation as a driver of epigenetic clock progression. Genome Biol. 23 (1), 24. doi:10.1186/s13059-021-02585-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kassambara, A. (2003) Ggpubr R package: ggplot2-based publication ready plots.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopez-Otin, C., Blasco, M. A., Partridge, L., Serrano, M., and Kroemer, G. (2013). The hallmarks of aging. Cell 153 (6), 1194–1217. doi:10.1016/j.cell.2013.05.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakamura, S., Kawai, K., Takeshita, Y., Honda, M., Takamura, T., Kaneko, S., et al. (2012). Identification of blood biomarkers of aging by transcript profiling of whole blood. Biochem. Biophys. Res. Commun. 418 (2), 313–318. doi:10.1016/j.bbrc.2012.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Nasi, M., Troiano, L., Lugli, E., Pinti, M., Ferraresi, R., Monterastelli, E., et al. (2006). Thymic output and functionality of the IL-7/IL-7 receptor system in centenarians: implications for the neolymphogenesis at the limit of human life. Aging Cell 5 (2), 167–175. doi:10.1111/j.1474-9726.2006.00204.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pellegrino-Coppola, D., Claringbould, A., Stutvoet, M., Consortium, B., Boomsma, D. I., Ikram, M. A., et al. (2021). Correction for both common and rare cell types in blood is important to identify genes that correlate with age. BMC Genomics 22 (1), 184. doi:10.1186/s12864-020-07344-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Peters, M. J., Joehanes, R., Pilling, L. C., Schurmann, C., Conneely, K. N., Powell, J., et al. (2015). The transcriptional landscape of age in human peripheral blood. Nat. Commun. 6, 8570. doi:10.1038/ncomms9570

PubMed Abstract | CrossRef Full Text | Google Scholar

Prince, M. J., Wu, F., Guo, Y., Gutierrez Robledo, L. M., O'Donnell, M., Sullivan, R., et al. (2015). The burden of disease in older people and implications for health policy and practice. Lancet 385 (9967), 549–562. doi:10.1016/S0140-6736(14)61347-7

PubMed Abstract | CrossRef Full Text | Google Scholar

R Core Team (2023) A language and environment for statistical computing.

Google Scholar

Tomusiak, A., Floro, A., Tiwari, R., Riley, R., Matsui, H., Andrews, N., et al. (2023). Development of a novel epigenetic clock resistant to changes in immune cell composition. bioRxiv. doi:10.1101/2023.03.01.530561

CrossRef Full Text | Google Scholar

Uhlen, M., Karlsson, M. J., Zhong, W., Tebani, A., Pou, C., Mikes, J., et al. (2019). A genome-wide transcriptomic analysis of protein-coding genes in human blood cells. Science. 366 (6472), eaax9198. doi:10.1126/science.aax9198

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Geest, K. S., Abdulahad, W. H., Teteloshvili, N., Tete, S. M., Peters, J. H., Horst, G., et al. (2015). Low-affinity TCR engagement drives IL-2-dependent post-thymic maintenance of naive CD4+ T cells in aged humans. Aging Cell 14 (5), 744–753. doi:10.1111/acel.12353

PubMed Abstract | CrossRef Full Text | Google Scholar

Vincent, F., Loria, P., Pregel, M., Stanton, R., Kitching, L., Nocka, K., et al. (2015). Developing predictive assays: the phenotypic screening "rule of 3. Sci. Transl. Med. 7 (293), 293ps15. doi:10.1126/scitranslmed.aab1201

PubMed Abstract | CrossRef Full Text | Google Scholar

Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L. D. A., François, R., et al. (2019). Welcome to the tidyverse. J. Open Source Softw. 4 (43), 1686. doi:10.21105/joss.01686

CrossRef Full Text | Google Scholar

Zhang, Z. (2014). Too much covariates in a multivariable model may cause the problem of overfitting. J. Thorac. Dis. 6 (9), E196–E197. doi:10.3978/j.issn.2072-1439.2014.08.33

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: ageing, gene expression, transcriptomics, ageing clocks, naïve T cell, leukocyte, whole blood, confounding variables

Citation: Fraser C and Owen BM (2024) Naïve T-cell decline is a significant contributor to expression changes in ageing blood. Front. Aging 5:1389789. doi: 10.3389/fragi.2024.1389789

Received: 23 February 2024; Accepted: 13 May 2024;
Published: 30 May 2024.

Edited by:

George A. Garinis, University of Crete, Greece

Reviewed by:

Lisa Wagar, University of California, Irvine, United States
Luca Pangrazzi, University of Innsbruck, Austria
Shouneng Peng, Icahn School of Medicine at Mount Sinai, United States
Ye-Eun Yoo, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea

Copyright © 2024 Fraser and Owen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Cameron Fraser, cameron.fraser@systematicmedicine.com, info@systematicmedicine.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.