Single-cell analyses of aging, inflammation and senescence.

Single-cell gene expression (transcriptomics) data are becoming robust and abundant, and are increasingly used to track organisms along their life-course. This allows investigation into how aging affects cellular transcriptomes, and how changes in transcriptomes may underlie aging, including chronic inflammation (inflammaging), immunosenescence and cellular senescence. We compiled and tabulated aging-related single-cell datasets published to date, collected and discussed relevant findings, and inspected some of these datasets ourselves. We specifically note insights that cannot (or not easily) be based on bulk data. For example, in some datasets, the fraction of cells expressing p16 (CDKN2A), one of the most prominent markers of cellular senescence, was reported to increase, in addition to its upregulated mean expression over all cells. Moreover, we found evidence for inflammatory processes in most datasets, some of these driven by specific cells of the immune system. Further, single-cell data are specifically useful to investigate whether transcriptional heterogeneity (also called noise or variability) increases with age, and many (but not all) studies in our review report an increase in such heterogeneity. Finally, we demonstrate some stability of marker gene expression patterns across closely similar studies and suggest that single-cell experiments may hold the key to provide detailed insights whenever interventions (countering aging, inflammation, senescence, disease, etc.) are affecting cells depending on cell type.

Introduction

Aging research offers exciting prospects. For example, in recent years, research on cellular senescence, one of the hallmarks of aging (Lopez-Otin et al., 2013), gained momentum with the finding that the elimination of senescent cells can improve health and survival in mice (Baar et al., 2017; Baker et al., 2016; Xu et al., 2018). Around the same time, single-cell analyses gained momentum (Svensson et al., 2018; Zappia et al., 2018). Now, both fields are increasingly coming together. As we will see, aging, inflammation and senescence have different effects on different cell-types, providing a direct justification for single-cell approaches. In turn, the processes underlying aging are often subtle, and will provide interesting challenges for the development of single-cell methods. And indeed, the development of molecular clocks, that is biomarker signatures predicting age, or aging, based on single-cell data requires an interdisciplinary approach. We start with a brief review of the development of the two fields as separate areas of research. In the main part, we will summarize, compare and contrast the results of recently published investigations of aging, aging-related chronic inflammation (inflammaging) and senescence (that is, cellular senescence and immunosenescence) based on single-cell RNA-sequencing (scRNA-seq). We will specifically focus on markers of cellular senescence. We do not consider other omics data, as these are still few and far between. We will describe multi-tissue studies and studies of whole blood, as well as of specific organ tissues, such as pancreas, brain, lung, skin and muscle. We will also describe some single-cell studies based on in-vitro material.

Investigating single-cell data with a focus on aging processes, or, more generally, along a time axis, all transcriptomics data are necessarily cross-sectional on the cell level: no single cell can be investigated twice. On the individual level, it is possible in principle to repeatedly take blood or tissue samples from the same individual. Longitudinal (multi-)omics studies have been done using bulk transcriptomics, e.g., (Schussler-Fiorenza Rose et al., 2019), but longitudinal in vivo studies at the single cell level, with relevance to aging processes, are still lacking; longitudinal samples of cancer biopsies were subjected to single-cell sequencing, however (Hamza et al., 2019; Maynard et al., 2019). Thus, the insights reported here regarding changes across time are cross-sectional, based on different human or animal donors for different time-points, who come with their inter-individual variability. This consideration is important regarding the signal in the data, referring to patterns of expression changes associated with aging processes; we would expect that this signal would increase if longitudinal data were investigated. Cross-sectionally at least, an increase in average heterogeneity of gene expression with age can be found in most of the single-cell data that are specifically investigated regarding this aspect; it is also known as transcriptional noise or variability and estimated in some single-cell studies as the correlation of gene expression patterns among cells and similar observations are long known from bulk data (Bahar et al., 2006; Işıldak et al., 2019; Martinez-Jimenez et al., 2017).

In some single-cell analyses, specific sets of marker genes are explicitly established to define cells, though in most cases, cell type definitions are done based on a priori published marker gene data. In that case, individual marker genes, or sets of these, can still be used to characterize the existing cells further, as well as new-found clusters of cells. Whether such further descriptions necessitate to define “sub-cell-types” or not may be debated; we just note that “cell types” are man-made constructs based on bona-fide boundaries or thresholds with an eye on utility and applicability, not scientific rigor. In our case, our default assumption is that aging processes do not imply the change of one cell type to another; rather we speak about a specific cell type in various states due to age- or aging-related change, which usually amounts to a deterioration of function. For example, we consider that cells of a certain cell type can be in “early” or “late” as well as “partial” or “full” states of, e.g., cellular senescence, depending on their expression of marker genes associated with cellular senescence.

Moreover, combining omics data with endpoint data, we can aim to identify biomarkers that predict these endpoints. Standard aging-related endpoints are related to mortality or morbidity, where the latter can be defined as reflecting disease and/or dysfunction. Biomarker candidates can then be some specific single-cell omics measurement, or sets (signatures) of such measurements, or sum scores and the like. A biomarker of aging then is expected to predict an endpoint better than chronological age (see also Fuellen et al. (2019)). Transferring these ideas to the single-cell world will require appropriate longitudinal single-cell datasets, including the clinical characterization of the donor animals or humans in terms of mortality or morbidity, and there will be exciting opportunities for biomarker discovery and validation once these data will become available (Rajewsky et al. (2020)).

Cellular senescence, inflammation and aging

While we review all aging-related single-cell datasets that we could identify, we will pay specific attention to chronic inflammation (inflammaging), and, most specifically, we will focus on the inflammaging-related hallmark of aging that is known as cellular senescence. Cellular senescence was first described almost 60 years ago when it was discovered that human diploid fibroblasts have a finite replicative potential in culture, after which the cells enter a state of irreversible replicative arrest (Hayflick and Moorhead, 1961). This phenomenon, later called replicative senescence, is linked to the constant loss of telomeric DNA that is associated with each cell division (Harley et al., 1990). Because of the 5′->3′ directionality of DNA-polymerases, the replication machinery is not able to copy the very ends of linear chromosomes. The ensuing loss of DNA, as time progresses, eventually triggers cell cycle checkpoints that preclude further cell divisions. This type of replicative senescence, also known as the Hayflick limit, exists in many species and correlates with maximum lifespan (Röhme, 1981). Consequently, it has given rise to the telomere theory of aging (Harley, 1991).

However, in recent years it became clear that telomere attrition is only one of several mechanisms that can induce cellular senescence. Various types of stress, such as reactive oxygen species, radiation, DNA damaging agents and the unfolded protein response are among such mechanisms (Campisi, 2013; Coppe et al., 2010). Thus, cellular senescence appears to be a special stress response state of the cell that is accompanied by broad morphological and biochemical changes (McHugh and Gil, 2018). Transcriptional changes include an up-regulation of tumor suppressor genes and a down-regulation of cell-cycle promoting genes (Chatsirisupachai et al., 2019; Wiley et al., 2017). Importantly, senescent cells frequently secrete an inflammatory mix of cytokines, growth factors and matrix metalloproteinases, which form the so-called senescence-associated secretory phenotype (SASP) (Coppe et al., 2010; de Keizer, 2017; Gorgoulis et al., 2019). This paracrine signaling has a range of negative effects involving tissue remodeling and also tumorigenesis, as well as chronic inflammation (inflammaging) (Coppe et al., 2010; Franceschi et al., 2000, 2018; Sun et al., 2018). In a nutshell, senescent cells are non-dividing but secretory, damaged, and metabolically altered cells, which frequently escape apoptosis. Molecular identification of senescent cells is not trivial, since the senescent state induced by different triggers in different tissues displays a certain degree of heterogeneity (Hernandez-Segura et al., 2017). However, key markers used for identification are a senescence-associated form of beta-(β)-galactosidase (Dimri et al., 1995) and p16 (CDKN2A), which is expressed by a majority of senescent cells (Campisi and d’Adda di Fagagna, 2007; Collins and Sedivy, 2003).

Clear support that senescent cells are not only correlated with aging and diseases, but instead are causally involved, as already suggested by Yanai and Fraifeld (2018), comes for instance from recent studies which transplanted senescent cells from old into young mice (Xu et al., 2018). This resulted in persistent functional impairment as well as spread of cellular senescence to host tissues. Another strong line of evidence comes from experiments that actually removed senescent cells from aged mice (Baar et al., 2017; Baker et al., 2016; Xu et al., 2018). In each case, an increase in lifespan and a delay of typical age related diseases has been observed. A good summary of the involvement of senescent cells in human diseases can be found in McHugh and Gil (2018), and a good yet brief account of the state of the art is van Deursen (2019). Cellular senescence is a frequently cited hallmark of aging (Lopez-Otin et al., 2013), and it is no surprise that it is implicated in the morbidity and co-morbidity of many (if not most or even all) age-associated diseases.

Single-cell analyses

Many fundamental questions in biology and medicine require the understanding of molecular processes at single-cell resolution. How does a zygote change and multiply to become a full-fledged organism made up of trillions of cells including many types of cells, tissues, and organs? Which cell within a tumor is the first ancestral cell that developed into a cancerous state and what molecular changes caused that cell to gain uncontrollable growth and proliferation? What kind of immune cells start accumulating at a damaged blood vessel and eventually cause the rupture of an arterial plaque that causes stroke or acute myocardial infarction? How close is a cell to reaching the Hayflick limit - is it still functional or is it in a senescent state? These and many other questions have been traditionally interrogated by studying single cells, but only few genes or few cells at a time, because, technologically it was not possible until recently to have an unbiased view of the whole genomes, transcriptomes, or epigenomes of single cells at an affordable manner at scale.

Whole transcriptome sequencing of a single cell was first achieved in 2009 (Tang et al., 2009). The single-cell sequencing technologies have since evolved quite rapidly from sequencing merely a manually picked individual cell (Tang et al., 2009), to massively multiplexed sequencing of hundreds of thousands of cells (Svensson et al., 2018). Even more recent protocols have achieved sequencing more than 2 million cells in one experiment at a cost of $0.01 per cell (Cao et al., 2019), which reflects the rapid decline in sequencing costs at single cell resolution.

Currently, a variety of competing single-cell RNA-seq protocols exist, each with its own advantages and trade-offs (Ding et al., 2019; Mereu et al., 2020; Svensson et al., 2017). Moreover, along with the sequencing protocols, single-cell sequencing data analysis methods have also seen a rapid rise. Multiple tools have been developed for almost every aspect of the data analysis. A resource (Zappia et al., 2018) that tracks the available single-cell sequencing data analysis tools from the scientific community, as of writing this manuscript, has tracked more than 630 tools; a new tool is added to the compendium almost every other day. Having many alternatives for each data analysis step is on one hand a luxury, on the other hand a challenge. Fortunately, there have been benchmarking studies carried out for each step of the workflow. Thus, one can make an informed decision about which tool to use depending on the requirements of the project at hand.

At least one benchmarking study exists for each step of the single-cell sequencing and data analysis workflow: sample preparation and sequencing (Mereu et al., 2020), alignment of raw reads (Baruzzo et al., 2017) or alignment-free quantification of read counts (Wu et al., 2018), count matrix pre-processing and normalization (Cole et al., 2019) -- this may include imputation for drop-outs (Zhang and Zhang, 2020), or batch correction (Tran et al., 2020) --, clustering (Duo et al., 2018), differential expression analysis to determine cluster-specific markers (Soneson and Robinson, 2018), cell type annotation (Abdelaal et al., 2019), dimension reduction (Sun et al., 2019), visualization, and downstream analyses, e.g., based on biological network data (Gawel et al., 2019; Mohammadi et al., 2019). As decisions made at each step of this workflow are likely to affect the downstream options and possible outcomes, combinatorial assessment of each component of this workflow has been assessed by various studies (Germain et al., 2020; Tian et al., 2019; Vieth et al., 2019). In some cases, a pre-defined set of tools and parameters can boost analysis speed and reproducibility (Wurmus et al., 2018). However, given the amount of tools and computational approaches and the fast development of the field, the analysis of single cell data is a moving target.

The common denominator of the benchmark studies could be summarized as that there is no one single combination of experimental protocols and computational tools that clearly dominate other options in all possible scenarios; the choices are context dependent, they need to be made according to the biological question at hand, and there must be multiple quality control checkpoints along the workflow. However, most of the single-cell sequencing data analysis challenges have been addressed, and the scientific community is converging onto a set of best practices for data analysis (Luecken and Theis, 2019). Thus, the last decade witnessed single-cell RNA-seq (scRNA-seq) becoming a mature, affordable, accessible technology suited to address the fundamental questions in biology. Finally, data from single-cell transcriptomics can be used to de-convolute bulk transcriptomics data (Dong et al., 2020), i.e., it is also a good companion to earlier technology and furthers the understanding of data already available.

Single-cell analyses of aging and senescence

Most publications on single-cell analyses of aging and senescence were published in the last two years. A first type of study is based on in-vivo data, investigating the effects of aging mostly in mice and humans, for a specific tissue and sometimes for a set of tissues (“multi-tissue”). In such studies, the data may be inspected specifically for markers of cellular senescence such as p16 (CDKN2A) or p21 (CDKN1A). A second type of study specifically investigates collections of senescent cells, in-vitro, identifying or inspecting specific biological mechanisms. We will first present the in-vivo studies.

Multi-tissue studies

At time of writing, two multi-tissue in-vivo mouse studies are available: the Tabula Muris Senis (Almanzar et al., 2020) and the Calico Study (Kimmel et al., 2019), see Table 1 for details.

Table 1

List of single-cell publications and data sets related to aging and senescence.

In-vivo datasets
Tabula Muris Senis (Almanzar et al., 2020)
Tissues:	bladder, bone marrow, brain (cerebellum, cortex, hippocampus and striatum), fat (brown, gonadal, mesenteric and subcutaneous), heart and aorta, kidney, large intestine, limb muscle and diaphragm, liver, lung, mammary gland, pancreas, skin, spleen, thymus, tongue and trachea
Animals/Cells:	∼530 000 cells from 30 male & female C57BL/6 J N mice: 6 age groups: 1 mo, 3 mo, 18 mo, 21 mo, 24 mo, 30 mo
Heterogeneity:	Not investigated
Exp. setup:	10x Genomics, also FACS-based for the 3 mo & 24 mo age groups; Nova-seq 6000
Remarks:	The number of p16-expressing cells doubled in old mice, and in the 10x Genomics data, p16 expression doubled as well.
Calico Study (Kimmel et al., 2019)
Tissues:	kidney, lung, spleen
Animals/Cells:	∼55 000 cells from 7 male C57Bl/6 mice: 7 mo, and 22/23 mo
Heterogeneity:	Not investigated
Exp. setup:	10X Genomics; HiSeq 4000
Remarks:	No change in p16 expression, nor of any specific senescence-related gene signature.
Immune Aging Atlas (Martinez-Jimenez et al., 2017)
Tissues:	blood (CD4+ T-cells)
Animals/Cells:	∼1 500 cells from 3 mo and 21 mo old male mice (C57/BL6 and CAST/EiJ strains; number of animals not stated)
Heterogeneity:	increasing cell-to-cell transcriptional variability, based on a Mann-Whitney-Wilcoxon test of distributional differences
Exp. setup:	Fluidigm protocol; Illumina HiSeq2500
Remarks:	Recording of immunological activation response, but no specific investigation of marker genes of aging or senescence.
Aging Hematopoietic Stem Cell Atlas (Grover et al., 2016)
Tissues:	blood: hematopoietic stem cells
Animals/Cells:	135 cells from male and female CD45.2 C57BL/6 mice: 2/3 mo and 20–25 mo
Heterogeneity:	Not investigated
Exp. setup:	Fluidigm C1 System, HiSeq 2000
Remarks:	HSC aging shows upregulation of platelet lineage gene numbers and expression
Aging T-cell Atlas (Elyahu et al., 2019)
Tissues:	blood (CD4+ T-cells)
Animals/Cells:	∼24 000 cells from 8 C57BL/6 mice (of unspecified gender), 2–3 mo and 22–24 mo
Heterogeneity:	Not investigated
Exp. setup:	10x Genomics GemCode Chromium, Illumina NextSeq 500
Remarks:	Investigation of subsets of cells associated with chronic inflammation and immunity decline, but no specific investigation of marker genes of aging or senescence.
Aging Artery Atlas (Zhang et al., 2020)
Tissues:	aorta and coronary arteries
Animals/Cells:	∼8 000 cells from 16 male & female cynomolgus monkeys, 4–6 y and 18–21 y
Heterogeneity:	increasing transcriptional noise, following the criterion of Enge et al. (2017)
Exp. setup:	10x Genomics, HiSeq4000
Remarks:	FOXO3A may be a key regulator of arterial aging.
Aging Pancreas Atlas (Enge et al., 2017)
Tissues:	pancreas
Animals/Cells:	2 544 cells from 8 human donors (both genders), 1 mo to 56 y old
Heterogeneity:	increasing transcriptional noise, following their criterion, but no correlation to mutational load
Exp. setup:	Illumina NextSeq, Smartseq2, also FACS
Remarks:	p16 was found upregulated in a higher fraction of cells.
Pancreatic Cancer Atlas (Elyada et al., 2019)
Tissues:	pancreatic cancer tissue + adjacent tissue
Animals/Cells:	21 200 cells from 6 patients, 64–87 y old (and 11 000 cells from 4 KPC model mice)
Heterogeneity:	Not investigated
Exp. setup:	10X Genomics Chromium; Illumina HiSeq 4000
Remarks:	Inflammatory CAFs were found, expressing interleukins and cytokines that are also characteristic for the SASP.
Aging/T1D Beta Cell Atlas (Thompson et al., 2019)
Tissues:	pancreas
Animals/Cells:	∼21 000 cells from 9 female NOD mice, 8 w, 14 w & 16 w
Heterogeneity:	Not investigated
Exp. setup:	10X Genomics Chromium; Illumina NovaSeq 6000
Remarks:	Islet cells become senescent and express SASP proteins (marked in particular by upregulation of p21 and Ccxl10).
Aging Brain Atlas (Ximerakis et al., 2019)
Tissues:	brain
Animals/Cells:	∼50 000 cells from 16 male C57BL/6 J mice, 2/3 mo and 21–23 mo
Heterogeneity:	no broadly increasing transcriptional variability, based on coefficient of variation of expression for selected sets of transcripts
Exp. setup:	10X Genomics Chromium, Illumina NextSeq500
Remarks:	Senescence pathway activation was found in a cluster of endothelial cells.
Aging Cerebromicrovascular Atlas (Kiss et al., 2020)
Tissues:	brain
Animals/Cells:	∼9000 cells from 6 male C57BL/6 mice, 3 mo and 28 mo
Heterogeneity:	Not investigated
Exp. setup:	10X Genomics Chromium; Illumina NovaSeq 6000
Remarks:	∼10 % of cells display cellular senescence, based on a precompiled list of markers.
Aging Microglia Atlas (Hammond et al., 2019)
Tissues:	brain
Animals/Cells:	∼76 000 microglia cells from male & female C57BL/6 J mice, embryonic, and at 4 d, 5 d, 30 d, 100 d and 540 d.
Heterogeneity:	increasing cellular diversity, potentially due to increasing transcriptional heterogeneity
Exp. setup:	10X Genomics Chromium, Illumina NextSeq500
Remarks:	Recording of the upregulation of genes linked to disease and inflammaging, but no specific investigation of marker genes of cellular senescence.
Aging Lung Atlas (Angelidis et al., 2019)
Tissues:	lung
Animals/Cells:	∼15 000 cells from 15 C57BL/6 mice, 3 mo and 24 mo
Heterogeneity:	increasing transcriptional noise, similar to Enge et al. (2017)
Exp. setup:	DropSeq; Illumina HiSeq4000
Remarks:	Increased cholesterol biosynthesis in type-2 pneumocytes and lipofibroblasts and altered relative frequency of airway epithelial cells were observed.
Lung Aging/Fibrosis Atlas (Reyfman et al., 2019)
Tissues:	lung
Animals/Cells:	∼76 000 cells from 16 humans (both genders), 8 lung tissue donors (21–63 y) and 8 patients (37–72 y) with pulmonary fibrosis
Heterogeneity:	Not investigated
Exp. setup:	Illumina HiSeq 4000; also bulk RNA-sequencing of flow-sorted cells
Remarks:	Senescent cells were identified as a rare cell population by a 1,311-gene-signature (but not by a 4-gene signature including p16).
Lung Fibrosis Mouse Model (Yao et al., 2019)
Tissues:	lung
Animals/Cells:	11 686 cells from 8 transgenic mice (both genders), 10–12 w
Heterogeneity:	Not investigated
Exp. setup	Illumina Novaseq 6000, FACS sorted cells, 10X Genomics Chromium
Remarks:	Alveolar type 2 (AT2) stem cells isolated from IPF lung tissue exhibit characteristic transcriptomic features of cellular senescence.
Aging Skin Atlas (Solé-Boldo et al., 2019)
Tissues:	skin
Animals/Cells:	∼15 000 cells from 5 male humans, 25–27 y and 53–70 y
Heterogeneity:	less clarity in GO annotations, potentially due to increasing transcriptional heterogeneity
Exp. setup:	10X Genomics Chromium, Illumina HiSeq 4000
Remarks:	An upregulation of “skin aging-associated secreted proteins” (SAASP) was observed.
Aging Muscle Stem Cell Atlas (Hernando-Herraez et al., 2019)
Tissues:	muscle
Animals/Cells:	377 cells from 4 male mice (transgenic Pax7-nGFP17 strain, on a (C57BL/6;DBA2 F1/JRj) genetic background), 1.5-2.1 mo and 23.3–27 mo
Heterogeneity:	increasing variability, based on the pairwise correlation of gene expression of cells, both within and between individual mice, correlated with methylation
Exp. setup:	FACS-sorting for muscle stem marker; bead transcription & amplification using Smartseq2
Remarks:	Quiescent cells were investigated, and accordingly, no change in expression of senescence markers Cdkn2a (p16) & Cdkn1b (p27) was observed.
In-vitrodatasets
HCA2 Fibroblast Atlas (Tang et al., 2019)
Tissues:	HCA2 fibroblasts, after replicative senescence after 38, 48 and 71 population doublings
Animals/Cells:	1 200 + 400 human cells
Heterogeneity:	Not investigated
Exp. setup:	Drop-seq; Library prep: Illumina Nextera XT
Remarks:	Investigation of senescence (beta-galactosidase, rate of EdU incorporation), role of telomeres.
Senescence Entry Cell Atlas (Zirkel et al., 2018)
Tissues:	HUVECs (mesodermal) IMR90 s (endodermal), and MSCs (multipotent)
Animals/Cells:	5 200 cells from 8 human donors & 2 cell lines
Heterogeneity:	Not investigated
Exp. setup:	10x Genomics Chromium; Illumina HiSeq 4000
Remarks:	Investigation of senescence (beta-galactosidase, reduced proliferation, 6 specific CpGs, HMGB2 and CTCF).
mTOR shRNA Atlas (Aarts et al., 2017)
Tissues:	primary human fibroblasts (IMR90) expressing reprogramming factors
Animals/Cells:	in-vitro shRNA screen
Heterogeneity:	Not investigated
Exp. setup:	Illumina Nextera XT; Illumina HiSeq 2500
Remarks:	Investigation of senescence (p16, p21, p19, SASP) & reprogramming and of the role of mTOR.

The tissues/animals/cells investigated are described as well as the experimental setup. Statements with respect to transcriptional heterogeneity, and observations regarding aging/inflammaging and cellular senescence are also provided. Abbreviations: d, day; w, week; mo, month; y, year; Exp, experimental.

The Tabula Muris Senis is an extension of the Tabula Muris study (Schaum et al., 2018) beyond the three-month (3 mo) time-point. It features 18 tissues at a maximum of 6 time points, from 1 month to 30 months, reporting gene expression for ∼530 000 cells. Of note, the fraction of cells expressing the senescence marker p16 (Cdkn2a) was found to double in older mice. A subsequent analysis of these data (Zhang et al., 2019) noted that many “global aging genes” (whose expression changes significantly with age in more than half of the cell types) are related to cellular senescence, citing App, Ctnnb1, Mapk1, Rac1, Arf1, and Junb in particular. On the pathway level, the (senescence-related) mTOR pathway is downregulated; this insight matches the overall trend of transcriptional downregulation, but it leads to some discussion in the paper regarding the plausibility of this finding, possibly as being adaptive. Also, Ubq is the only well-characterized gene among the top 20 global aging genes, suggesting that proteostasis plays a prominent role overall in mouse aging. We took the Tabula Muris Senis as the starting point to demonstrate that this data set is, like most others, a treasure trove that can (still) be mined in various ways, generating further insights of relevance to inflammaging and cellular senescence, see Box 1 .

We processed the Tabula Muris Senis data (Almanzar et al., 2020) (see Supplement) and investigated 10 senescence-associated genes that are expressed in more cells from old mice, as compared to young mice (see also Box 2). Curiously, the by far largest differences in the fraction of cells expressing these markers can be found in the liver, and among the cells of the liver, in the Kupffer cells (see Fig. 1 ). In particular, Interleukin 1b (IL1b), which is one of the cytokines involved in inflammaging (Franceschi et al., 2017) and a component of the senescence-associated secretory phenotype (SASP), can be found in almost 100 % of Kupffer cells in old mice, but only in ∼37 % of these cells in young mice. We find that IL1b is also upregulated in other liver cells, e.g., hepatic stellate cells, but only in 50 % of these. Kupffer cells are resident liver macrophages (Hunt et al., 2019). Aged Kupffer cells have a reduced capacity for phagocytosis and autophagy, and display an increased inflammatory phenotype contributing to inflammaging (Hunt et al., 2019). Further, IL1b was shown to be involved in inflammation and aging in the liver context, at least in case of alcoholic steato-hepatitis (Morsiani et al., 2019). As a small caveat, we acknowledge that merely looking into the differences in terms of fraction of cells expressing a marker between the young and old animals might not accurately reflect the actual amount of cellular senescence. For that purpose, one needs to additionally account for inter-individual and inter-generational cell type compositional differences, which is beyond the scope of this review.

Fig. 1

Fraction of cells expressing 10 aging-related marker genes in liver cells of the Tabula Muris Senis.

Alt-text: Box 1

The Calico Study (Kimmel et al., 2019) analyzed ∼55 000 single cells from lung, kidney and spleen of young (7 mo) and old (21 mo) mice. They find that age-related changes of gene expression are largely cell type specific and depend less on the type of tissue. However, they also found some 200 genes that are differentially expressed in more than five cell types. These genes are preferentially involved in transport of proteins to the endoplasmic reticulum as well as in inflammation. Inflammation pathways were found upregulated in immune as well as non-immune cells, in agreement with inflammaging. Interestingly, this study did not find an age related increase of the fraction of cells that express p16, nor of p16-expression on average.

Single-tissue studies

Further studies up to now focused on the blood (Elyahu et al., 2019; Grover et al., 2016; Martinez-Jimenez et al., 2017), the pancreas (Elyada et al., 2019; Enge et al., 2017; Thompson et al., 2019), the brain (Hammond et al., 2019; Kiss et al., 2020; Ximerakis et al., 2019) and the lung (Angelidis et al., 2019; Reyfman et al., 2019), see Table 1 and the next sections. Of note, comparing data of the Tabula Muris Senis and some of the single-tissue studies, we note high consistency among studies when looking at aging/senescence-related markers, see Box 2 .

For this review, we checked the consistency of marker gene expression across studies. In particular, we compared the Tabula Muris Senis lung cell data to the Aging Lung Atlas (Angelidis et al., 2019), and the Tabula Muris Senis T cell data (bone marrow, spleen, and thymus) to the Aging T-cell Atlas (Elyahu et al., 2019), for 10 senescence-associated genes that are expressed in more cells from old mice, as compared to young mice (see also Box 1). As can be seen in the Fig. 2 , the cell fractions, as well as the increases in cell fraction in old mice, are roughly the same in case of the T cells, and the increases are mostly comparable in case of the lung cells. In case of lung, the cell fractions featuring a marker are, however, much lower overall than in the Tabula Muris Senis, which can be explained by the differences in the depth of sequencing between the two studies (an average of 2200 genes per cell were detected in the Tabula Muris Senis data, while only 446 genes per cell were detected in the Aging Lung Atlas).

Fig. 2

Fraction of cells expressing 10 markers in lung cells and T cells of the Tabula Muris Senis, compared to the Aging Lung Atlas and the Aging T-cell Atlas.

Alt-text: Box 2

Blood & vascular system

Due to its easy accessibility, among the first “tissues” inspected in a single-cell fashion was blood.

In what we call an “Immune Aging Atlas”, (Martinez-Jimenez et al., 2017) studied naive CD4+ T cells in young (3 mo) and old (21 mo) mice by investigating the transcriptional pattern of ∼1500 cells before and after in-vitro immune stimulation. In old mice, most genes involved in this activation process responded to the stimulus, but with a reduced magnitude compared to young mice, reflecting immunosenescence. The authors also found that old mice showed a larger cell-to-cell transcriptional variability. Moreover, ∼10 % of the genes were differentially expressed between young and old, and these genes were often only expressed in a small subset of cells. This may lead to the observed variability, which the authors see as a hallmark of the aging process. The data of this paper was recently re-analyzed by (Hao et al., 2019), who proposed an “optimal gene filtering algorithm for scRNAseq data” (OGFSC). Using their algorithm, the originally reported reduction of the transcriptional response in old mice could no longer be established. This highlights the influence that the choice of a specific bioinformatics method can have. Thus, not only differences of the experimental conditions, but also differences in the analysis pipeline can affect the reproducibility of results.

An “Aging Hematopoietic Stem Cell (HSC) Atlas” (Grover et al., 2016) of young (2–3 month) and old (20–25 month) mice investigated age dependent transcriptional signatures, showing that the previously observed age associated decrease of lymphoid cells and erythrocytes is accompanied by the accumulation of a specific “old” HSC subgroup of cells showing an upregulation of platelet lineage gene expression. The respective HSCs predominantly generate platelets but recover their ability to increase lymphoid cell generation when the key platelet programming factor FOG-1 is compromised. The authors discovered an aging-related subclass of HSC that explains the apparent decline in engrafting ability, exploiting the single-cell resolution of their data.

The “Aging T-Cell Atlas” (Elyahu et al., 2019) investigated the transcriptomes of ∼24 000 CD4+ T cells in young and old mice and reports the drop of naive CD4+ cells typical for immunosenescence, but also a surprising accumulation of three CD4+ subtypes that occur only rarely in young animals, namely cytotoxic, exhausted and activated regulatory T cells, characterized in part by their expression of inflammatory cytokines (primarily IL-27, IFNb, and IL-6), possibly reflecting inflammaging in the old mice.

In what we refer to as an “Aging Artery Atlas”, Zhang et al. (2020) studied the vascular system of cynomolgus monkeys at the single-cell level to investigate differences between aortas and coronary arteries in young (∼6 year-old) and old (∼20 year-old) animals. Using almost 8000 endothelial cells, they identified on average 3100 genes per cell, and among these they found 515 downregulated vs. 259 upregulated genes that were differentially expressed in at least one type of vascular cell. Bioinformatics analysis indicated that FOXO3A is a key regulator of arterial aging and inactivation of this transcription factor in human endothelial cells in vitro reproduced many of the phenotypic changes observed in aged monkey arteries. Also in this study, an increase of transcriptional noise associated with age was found. Following Enge et al. (2017), transcriptional noise was defined by the ratio of biological to technical variation, based on correlating gene expression in individual cells with the average of the corresponding cell type; technical variation was in turn based on a spike-in approach.

Pancreas

The human pancreas is a focus of surprisingly many single-cell studies, and at the same time, cellular senescence is implicated in many aspects of development and therapy of pancreatic cancer (specifically PDAC, pancreatic ductal adenocarcinoma) and diabetes (both type 1 and type 2) (Elyada et al., 2019; Enge et al., 2017; Thompson et al., 2019).

An “Aging Pancreas Atlas” was thus developed (Enge et al., 2017) by performing single cell transcriptomics of ∼2500 human pancreas cells from donors spanning six decades of age. The authors found an increasing transcriptional noise (defined as described above, see the “Aging Artery Atlas”) and an increasing mutational load with age, but they concluded that the number of mutations per cell is insufficient to explain the observed transcriptional noise, though epigenetic changes may do so. The senescence marker p16 (CDKN2A) was investigated as a well-known senescence marker and was found in a higher fraction of cells in old compared to young donors, though its expression was not found to increase on average (over all cells). We may speculate that, even though the fraction of senescent cells increases in the aging pancreas, the cells surrounding these are not affected due to large amounts of extracellular matrix, which may block any paracrine effects. The Aging Pancreas Atlas data were further analyzed by (Chen et al., 2019), but without any consideration of age, with the aim of learning classification rules for six different pancreatic cell types. For this purpose they used various machine learning techniques such as Monte Carlo feature selection and support vector machines. Using only the best 7 features (biomarkers) they reached a classification accuracy of 0.966. It would be of interest to use the same or similar methods to learn and test classification rules for predicting donor age, compare their performance, and to understand which genes underlie such classification.

As a more specialized resource, the “Pancreatic Cancer Atlas” is a detailed study of various populations of cancer-associated fibroblasts (CAFs) implicated in PDAC (Elyada et al., 2019). Often, CAFs promote tumor growth by modifying collagen cross-linking, secreting diverse cytokines, chemokines and growth factors and recruiting immune-suppressive cells, while other studies indicated that CAFs actually restrict tumor growth (Ozdemir et al., 2014; Rhim et al., 2014). The single-cell studies by (Elyada et al., 2019), investigating human and mouse cells, help to come closer to an explanation, by showing that CAFs are actually a heterogeneous cell population consisting of different subtypes such as myofibroblastic CAFs (myCAF), inflammatory CAFs (iCAF) and antigen presenting CAFs (apCAF). Regarding inflammaging and cellular senescence, iCAFs are worth noting, since they express interleukins and cytokines that are also characteristic for the SASP.

Finally, the “Aging/T1D Beta Cell Atlas” (Thompson et al., 2019) describes how cellular senescence of islet cells is involved in the development of type 1 diabetes (T1D) in humans and in the non-obese diabetic (NOD) mouse model. For this purpose the authors performed scRNA-seq of ∼21 000 cells of 8, 14 and 16 week old NOD mice; specifically they collected CD45+ immune-cell-depleted but insulin-expressing cells from the pancreas islets. Remarkably, the authors find that removal of such cells via senolytic compounds (inhibiting BCL-2) was able to prevent the appearance of diabetes in these mice, pointing to a failure of the immune system to remove such senescent cells in case of T1D and to senolysis as a remedy. In the single-cell data, the senescence marker p21 (Cdkn1a) was found enriched in a cluster of cells in the diseased state at 14/16 weeks, whose overall transcriptional profile was “reminiscent of a senescence stress response”. In turn, the beta cell maturation marker Ucn3 was found depleted in that cluster.

Brain

Single cell studies are especially useful for studying the brain, since this complex organ contains a multitude of different cell types that are difficult to separate. However, brain omics of human tissue are almost always and necessarily post-mortem, which may pose problems caused by artifact patterns (Ferreira et al., 2018).

In the “Aging Brain Atlas”, brain tissue of young and old mice was studied by single-cell transcriptomics of ∼50 000 cells (Ximerakis et al., 2019). By differential gene expression analysis, the authors found that expression patterns do not change uniformly across cell types. The coefficient of variation of expression for selected sets of transcribed genes differed between young and old cells in many cell types, but the direction of change was not consistent among cell types, suggesting no broad increase of transcriptional variability. The latter observation makes this study the only one in this review that does not report a general increase in transcriptional heterogeneity. Another noteworthy result is that in most brain cell types the transcription of ribosomal proteins appears to be upregulated. In their study the authors also found that aging- and senescence-related pathways were upregulated, specifically in a cluster of cells identified as endothelial cells. A web interface is available, at https://portals.broadinstitute.org/single_cell/study/aging-mouse-brain.

Another study investigating the brain is the “Aging Cerebromicrovascular Atlas” (Kiss et al., 2020). It describes the age-induced dysregulation of microvascular endothelial function in the mouse brain, which contributes to diseases such as vascular cognitive impairment and Alzheimer’s disease, based on a new approach to identify senescent cerebromicrovascular endothelial cells. Extracted and digested brain tissues isolated from mice were separated by centrifugation and filtering, and senescent cerebromicrovascular endothelial cells were identified through clustering by specific expression patterns and by a precompiled list of senescence markers such as Cdkn2a and Bmi1. ∼10 % of the cerebromicrovascular endothelial cells experience cellular senescence at a biological age equivalent to a 75-year-old human; microglia are similarly affected.

Finally, the “Aging Microglia Atlas” (Hammond et al., 2019) describes microglia, the macrophages of the brain, that have been historically classified as resting, M1 proinflammatory, or M2 anti-inflammatory, based on in vitro stimulation methods. As the resting-M1-M2 classification system might be oversimplifying the real phenotypic heterogeneity of these cells in vivo, Hammond et al. (2019) carried out single-cell transcriptomics of mouse microglia ranging from embryogenesis to old age. They found that microglia subpopulations became less heterogeneous in adulthood, but heterogeneity increased again during aging. Some aged microglia expressed the Ccl4 chemokine, and based on the authors observations, we suggest that this cytokine contributes to inflammaging.

Lung

The lung has become a key tissue of recent single-cell investigations due to the Covid-19 pandemic. For example, we recently described (Fuellen et al., 2020) a single-cell dataset investigating the lung immune microenvironment, based on bronchoalveolar lavage fluid of Covid-19 patients (Liao et al., 2020). The arrival of longitudinal single-cell data helping to investigate the interconnections between Covid-19 disease susceptibility & sequelae, in particular inflammaging, is only a question of time.

Investigating lung aging in mice (“Aging Lung Atlas”), scRNA-seq was combined with in-depth proteomics by (Angelidis et al., 2019). They sequenced ∼15 000 cells of young (3 month) and old (24 month) animals and identified 30 different cell types. Their mass spectrometry based proteomics is based on bulk tissue, but it could still be related to their scRNA-seq data, finding a significant correlation between fold changes at the transcriptome and proteome level. However, there were notable exceptions such as collagen IV, which showed an age-related decline at the mRNA level, but an increase at the protein level. As in (Enge et al., 2017), increased transcriptional noise was noted. Enrichment analyses pointed to increased cholesterol biosynthesis in some cell types, and an altered relative frequency of airway epithelial cells was found as a hallmark of lung aging. Finally, chronic inflammation / inflammaging was observed, putatively related to the deregulation of lipid homeostasis and epithelial cellular senescence. The “Aging Lung Atlas” is available online and can be searched via a webtool (see https://theislab.github.io/LungAgingAtlas/).

A human “Lung Aging/Fibrosis Atlas” is presented by (Reyfman et al., 2019) based on the analysis of ∼76 000 single cells from eight healthy donors and eight patients with pulmonary fibrosis, of various ages. In their clustering analysis (by t-Distributed Stochastic Neighbor Embedding, t-SNE), the authors did not find a cluster of cells that specifically express senescence markers (they inspected p16, IL-6, Serpine1 and GLB), but a 1311-gene signature of replicative senescence scored highest in a population of alveolar type II stem cells that originated almost exclusively from the fibrotic lung tissue. Furthermore, this gene signature could be seen in a small fraction of practically all cell types present in the lung. This is consistent with the finding that many cell types can develop a senescent phenotype, but that this phenotype is very heterogeneous depending on the original cell type (Hernandez-Segura et al., 2017). Notably, a recent study using a transgenic mouse model investigated the role of alveolar type II stem cells for lung fibrosis further by artificially inducing cellular senescence in this cell type (Yao et al., 2019), initiating progressive lung fibrosis, which could be blocked by a senolytic treatment with dasatinib and quercetin. In this loss-of-function study, the authors performed scRNA-seq of ∼10 000 epithelial cells collected by fluorescence-activated cell sorting (FACS).

Skin

Studies that investigate the aging process using skin have to take special care, since “photo-aging” due to the exposure of humans to sunlight can strongly influence aging-related processes, acting as a major confounder. Thus, Solé-Boldo et al. (2019) collected cells from the sun protected inguinoiliac skin region of male human donors for their “Skin Aging Atlas”, comprising more than 15 000 cells. They found 14 known cell types in their sample and among these, they could identify 4 subtypes of fibroblasts based on clustering (t-SNE) plots and analysis of GO term enrichments, which they called secretory-papillary fibroblasts, secretory-reticular fibroblasts, mesenchymal fibroblasts and pro-inflammatory fibroblasts. Furthermore, they compared the samples from young (25 & 27 years) to old (53, 69, 70 years) donors to investigate age-dependent changes, revealing that fibroblasts feature a subtype-specific set of “skin-aging associated secreted proteins” (SAASP), as a set of proteins secreted by the fibroblasts from old human donors and already described previously (Waldera Lupa et al., 2015). The functional identity of cells appeared to become blurred, since fewer GO terms with less significant p-values were found in case of the old samples; higher transcriptional heterogeneity may be an explanation for this observation.

Muscle

In their “Muscle Aging Atlas”, Hernando-Herraez et al. (2019) used 377 single muscle stem cells from young (∼2 months) and old (∼25 months) mice to measure the transcriptome as well as DNA methylation, in the same cells at the same time. Muscle stem cells were specifically selected as the subject of the study, also motivated by their assumed “deep” quiescence; a transgene was used to sort out these cells based on the marker Pax7. Based on the pairwise correlation of gene expression of cells, both within and between individual mice, the authors found an increased variability at the transcriptome level comparing young versus old mice; some immunostaining was done and an increased heterogeneity on the protein level was also observed. Further, the correlation between the degree of methylation of promoter regions, and gene expression, at the single-cell level was negative, supporting the idea that age related changes in the methylome may be triggering the age related changes in the transcriptome. Assessing epigenetic age using adult muscle stem cells yielded surprising results, though. While the estimated age in the young mice was around 9 weeks, the age prediction for the muscle stem cells of the old mice was 10–11 weeks. Thus, epigenetic clocks based on bulk tissue may not confer the full picture, if some cell types (e.g. muscle stem cells) maintain a young epigenetic age, but dwindle in number in old individuals. Further studies involving single-cell methylation data from different (stem) cell populations are needed to deeper investigate single-cell-based epigenetic clocks and their ability to predict chronological or biological age.

In-vitro studies

Telomere erosion was investigated by the in-vitro study of senescent cells reported in (Tang et al., 2019). They cultured HCA2 fibroblasts until replicative senescence was reached and performed single cell transcriptomics on a total of 1200 cells from early, middle and late passage. They also irradiated early passage cells to induce SIPS (stress induced premature senescence) and analyzed 400 of those cells. Visualization via t-SNE clearly separated cells from different passages and especially late passage senescent cells appear to be more heterogeneous than irradiation-induced senescent cells. This is a further confirmation that, to a certain degree, cellular senescence is not only different in different cell types, but also depends on the stressor that caused the senescent state. Nevertheless, (Tang et al., 2019) identify a senescence core signature of up-regulated genes that include p21 (CDKN1A) and members of the SASP.

The entry into senescence was also investigated in another in-vitro study using umbilical vein endothelial cells (HUVECs), fetal lung fibroblasts (IMR90 s) and mesenchymal stromal cells (MSCs) (Zirkel et al., 2018). The levels of HMGB2, which is involved in the maintenance of the 3D chromatin structure, diminish in the nuclei of old and senescent tissues and the authors hypothesized that this decline triggers an entry into cellular senescence. They demonstrate the decline of HMGB2 before the appearance of established senescence markers, by scRNA-seq on 8300 proliferating and 5200 senescent HUVECS. Cellular senescence can also be caused by transcription-factor-induced reprogramming stress, (Banito et al., 2009), which can lower the efficiency of iPSC generation, and a combination of an shRNA library intervention approach with single-cell sequencing was presented by (Aarts et al., 2017), based on primary human fibroblasts (IMR90) expressing reprogramming factors. This allowed to identify those shRNAs that are most suitable to suppress this type of cellular senescence. The four top scoring target genes were CDKN1A, MTOR, MYOT and UBE2E1. While CDKN1A (p21), and MTOR are known to influence senescence, the other two, being part of sarcomeres and being an E2 ubiquitin-conjugating enzyme, were new candidates.

Overall, it is an open yet relevant question how faithfully the in-vitro models of cellular senescence reflect the all-important in-vivo situation, and single-cell analyses may enable a more fine-grained comparison. At least, the markers of cellular senescence first identified in-vitro are often, but not always, following the same or similar plausible patterns of age-related change in gene expression in-vivo, for bulk and single-cell data alike.

Conclusions

Overall, an increasing frequency of cells expressing aging/senescence markers, and/or an overall increasing expression of aging/senescence markers (albeit often at low absolute levels in the cells in which they are expressed), and an increase in transcriptional heterogeneity are the most frequent patterns revealed by the aging-related single-cell studies we reviewed. A potential corollary of the latter is that cellular identity tends to become less constrained, or even lost. This observation matches a number of recent insights, outside of single-cell transcriptomics. For example, single-cell investigations of global chromatin modification data revealed an increased epigenetic variation with aging (Cheung et al., 2018). Moreover, the loss of cell identity due to chromatin modifiers attending double-strand breaks was recently suggested to be one of the main drivers of aging (Yang et al., 2019). This proposal connects transcriptional heterogeneity to cellular senescence, which is often caused by double-strand breaks and other kinds of DNA damage.

There is a lack of interventional studies supported by single-cell transcriptomics. This topic will be fascinating to investigate, promising insights into the drug response for different cell types for the same individual at identical conditions. Another prospect lies in the investigation of transcriptomic profiles that are preserved when bringing these cells in culture. This ties observational studies to experimental follow-up investigations, given that a conserved “context-free” subset of characteristic gene expression changes may be a primary target for new drugs to modulate aging-associated (and other) phenotypes (Struckmann et al., 2020). And likely, such sets of co-regulated genes, established for specific cell types, will be recognized as motifs in their own right, with their own functional annotation, in analogy to what domain databases have contributed to our understanding of biological sequences.

Ten years back, de Magalhaes et al. (2010) reviewed the state of Next-Generation Sequencing in aging research in this journal. Technological progress has been breathtaking since then, and Next-Generation Sequencing has become possible at single-cell resolution, on a large scale. Aging affects every cell, and to some degree every cell is affected differently by aging. In turn, all cells together make up the aging metazoan organism, and their contribution to the aging phenotype is not uniform. It will be exciting to watch the progress that can be made by the magnifying glass provided by scRNA-seq, to aging research and beyond.

Funding

GF is supported by the (FKZ 01ZX1903A, 03V0396 and 1DS19049), the (Aging with elegans, grant agreement 633589), the (FU-583/5-1), Karls Erdbeerhof, Rövershagen, Germany, and the Interdisziplinäre Fakultät (Department AGIS) of Rostock University.

Declaration of Competing Interest

The authors declare that no competing interests exist.