Cameron Fraser
Brady M. Owen- 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:
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. 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. 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. 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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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, GreeceReviewed by:
Lisa Wagar, University of California, Irvine, United StatesLuca 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