key: cord-0016537-06qyxg29
authors: Kim, Sohee; Kim, Chuna
title: Transcriptomic Analysis of Cellular Senescence: One Step Closer to Senescence Atlas
date: 2021-03-31
journal: Mol Cells
DOI: 10.14348/molcells.2021.2239
sha: 11bbb4cf71359f17eb5d52b5e01eb07249c70027
doc_id: 16537
cord_uid: 06qyxg29

Senescent cells that gradually accumulate during aging are one of the leading causes of aging. While senolytics can improve aging in humans as well as mice by specifically eliminating senescent cells, the effect of the senolytics varies in different cell types, suggesting variations in senescence. Various factors can induce cellular senescence, and the rate of accumulation of senescent cells differ depending on the organ. In addition, since the heterogeneity is due to the spatiotemporal context of senescent cells, in vivo studies are needed to increase the understanding of senescent cells. Since current methods are often unable to distinguish senescent cells from other cells, efforts are being made to find markers commonly expressed in senescent cells using bulk RNA-sequencing. Moreover, single-cell RNA (scRNA) sequencing, which analyzes the transcripts of each cell, has been utilized to understand the in vivo characteristics of the rare senescent cells. Recently, transcriptomic cell atlases for each organ using this technology have been published in various species. Novel senescent cells that do not express previously established marker genes have been discovered in some organs. However, there is still insufficient information on senescent cells due to the limited throughput of the scRNA sequencing technology. Therefore, it is necessary to improve the throughput of the scRNA sequencing technology or develop a way to enrich the rare senescent cells. The in vivo senescent cell atlas that is established using rapidly developing single-cell technologies will contribute to the precise rejuvenation by specifically removing senescent cells in each tissue and individual.

Various types of damages, such as those in tissues, cells, and molecules, accumulate during aging and often lead to cellular senescence (Soares et al., 2014; Zhang et al., 2015) . Senescent cells are secretory cells that are still metabolically active, although their cell cycle is stably stopped (Dörr et al., 2013) . Senescent cells accumulate with age and have a significant impact on the health and longevity of an individual (Biran et al., 2017; Burd et al., 2013; Jeyapalan et al., 2007; Zhu et al., 2014) . Cellular senescence is induced by various stresses. For example, telomere shortening occurs due to repetitive cell division, leading to a kind of senescence called replicative senescence, or Hayflick's limit (Biran et al., 2017; Burd et al., 2013; Jeyapalan et al., 2007; Zhu et al., 2014) . Other stresses, including the activation of oncogenes such as Ras, the inhibition of tumor suppressor genes such as p53, the accumulation of DNA damage, chromatin disruption, and reactive oxygen species, can trigger senescence.

However, cellular senescence can also occur in specific physiological conditions, such as development and wound healing (Demaria et al., 2014; Storer et al., 2013) . Senescent cells affect the surrounding cells by secreting senescence-associated secretory phenotypes (SASPs), including cytokines, chemokines, growth factors, and matrix metalloproteases (Acosta et al., 2008; Coppé et al., 2008; Krtolica et al., 2001) . Depending on the physiological context, SASPs can be involved in diverse mechanisms, such as immune cell recruitment, inflammation, and extracellular matrix (ECM) remodeling (Borodkina et al., 2018; Coppé et al., 2008; Xu et al., 2015) . For example, SASPs can result in immune activation, growth arrest, and differentiation, processes that are essential for tissue restoration. SASPs can also result in cell growth, migration, and invasion, which have adverse effects on cancer. This context difference also exists in young and old tissues. In young tissues, the effects of SASPs are temporary because senescent cells can be eliminated by the immune cells recruited by pro-inflammatory SASPs (Demaria et al., 2014; Kang et al., 2011; Xue et al., 2007) ; thus, senescent cells facilitate wound healing, tumor suppression, and the restoration of tissue homeostasis (Campisi, 2001; Demaria et al., 2014; Storer et al., 2013) . Conversely, as we get older, the removal of senescent cells becomes slower, resulting in the gradual accumulation of senescent cells. In other words, the effect of SASPs persists, likely resulting in chronic inflammation, aging-related deterioration or disease, and tumorigenesis (Chen et al., 2020; Coppé et al., 2008; Krizhanovsky et al., 2008; Nathan and Ding, 2010) . It remains unclear whether the accumulation of senescent cells inhibits the immune surveillance system or if the aging of the immune surveillance system leads to the accumulation of senescent cells (Ray and Yung, 2018) .

While senescence has a beneficial aspect in facilitating development and repairing damage, it has another aspect that induces organ dysfunction and aging. Therefore, there have been various attempts to reverse aging by regulating the accumulation of senescent cells and SASPs (Paez-Ribes et al., 2019) . Among these efforts, senolytics attempt to improve aging by selectively eliminating senescent cells. For example, when p16-expressing senescent cells are specifically eliminated in aged mice (p16-3MR, INK/ATTAC, etc.), aging-related dysfunctions are improved, and healthspan is extended (Pajvani et al., 2005) . Most senolytic drugs induce the death of apoptosis-resistant senescent cells via a pro-apoptotic mechanism (Lozano-Torres et al., 2019) . Recently, senolytics have been extended to chimeric antigen receptor (CAR) T-cell therapy that precisely target surface molecules expressed on various senescent cells. CAR T-cell therapy could inhibit liver fibrosis and restored liver function in old mice (Amor et al., 2020) .

However, this strategy does not work effectively against all types of senescent cells since cellular senescence is induced by various mechanisms to have different transcription signatures (Hernandez-Segura et al., 2017) . In particular, p16, a marker for senescent cells, is not only universally expressed in various senescent cell types but is also expressed after the cell cycle is temporarily stopped (Marthandan et al., 2014) . In addition, the representative senolytic drugs, such as dasatinib, quercetin, and navitoclax, had different effects depending on the cell type (Zhu et al., 2015; . These data suggest that some degree of heterogeneity exists in senescent cells. Therefore, it would be essential to consider the heterogeneity of senescent cells to eliminate senescent cells more efficiently.

When irreparable DNA damage occurs in a cell, it dies through mechanisms such as apoptosis or necrosis or activates senescence (Childs et al., 2014) . The cell cycle of senescent cells is arrested through the DNA-damage-response dependent p53/p21 pathway or the p16 INK4a /RB pathway (Serrano et al., 1993) . The factors determining the choice between apoptosis and senescence remain elusive. Still, given that the pro-apoptotic factors PUMA and BIM are upregulated in senescent cells (Baar et al., 2017) , the senescent cell is likely committed to apoptosis but unable to enter apoptosis. Accordingly, in senescent cells, pro-inflammatory SASP components, such as interleukin (IL)-6, IL-8, and IL-1, suppress apoptosis via NF-κB by regulating the expression of anti-apoptotic B cell lymphoma proteins (Bcl) such as Bcl-2 and Bcl-Xl (Catz and Johnson, 2001; Gabellini et al., 2008; Liu et al., 2016; Lozano-Torres et al., 2019) . Additionally, the increase in the translational burden of the SASPs in senescent cells leads to the unfolded protein response, which is one of the major causes of senescence-specific morphological changes, including flattened shape, enlarged nuclear and nucleolus, and accumulated lysosomal senescence-associated β-galactosidase (SA-β-gal) (Druelle et al., 2016; Muñoz-Espín and Serrano, 2014) .

Various methods have been developed to detect senescent cells based on their characteristics. The SA-β-gal staining assay and quantitative analysis of biomarkers, such as p16, p21, and p53, are commonly used (Dimri et al., 1995; Krishnamurthy et al., 2004; Ressler et al., 2006; Sharpless and Sherr, 2015) . However, these methods also detect non-senescent cells. For example, both quiescent cells, which are induced by serum starvation or pH change in a culture environment, and confluent cells are also stained positive in the SA-β gal staining assay (Yang and Hu, 2005) . In addition, p16 is highly expressed in non-senescent macrophages and pRb-negative cells (Benassi et al., 1999; Hall et al., 2017; Hara et al., 1996) . Furthermore, while premature senescence, induced by the telomere deprotection, involves p53/p21 and p16 pathways in human cells, it only involves the p53/p21 pathway in mouse cells (Smogorzewska and de Lange, 2002) . In humans, during aging, p16 and p21 show tissue-specific expression patterns (Casella et al., 2019; Idda et al., 2020) . Therefore, the current detection method is not only ineffective in distinguishing between the senescent and non-senescent cells but also unable to illustrate the molecular mechanisms of senescence between species. Further understanding of senescent cells is required to develop methods that detect senescent cells specifically.

Several studies have been tried to identify common markers observed in all senescent cells through unbiased transcriptomic in vitro analysis ( Fig. 1, Tables 1 and 2). The differentially expressed genes (DEGs) analysis of various human fibroblasts induced by different stimuli revealed senescence transcriptomic signatures, which were commonly changed in all fibroblasts and did not overlap with the transcriptomic signatures of quiescent cells. The core senescence signature was defined by comparing the signatures between fibroblast and the three other cell types from DEG analyses. Interestingly, this signature did not include classical senescence marker genes, such as p16, p21, LMNB1, and members of SASP, but rather genes related to DNA damage checkpoints and mitosis, such as PLK3, CCND1, DYNLT3, and CHMP5 (Hernandez-Segura et al., 2017).

In addition, in another independently performed transcriptomic analysis, universally expressed genes were found in various cell lines, such as lung fibroblast, umbilical vein endothelial, and alveolar endothelial cells, using various cellular senescence induction stimuli, such as replicative, oncogene-induced, ionizing radiation-induced, or doxorubicin-induced senescence. According to this analysis, the expression level of senescence markers, such as p16 and p21, slightly or did not commonly change (Casella et al., 2019) . When comparing these in vitro senescence transcriptomes, only 7 genes were consistently changed regardless of cell type or senescence induction method (Fig. 1, Table 1 ). These analyses suggest that previously known senescence markers might be insufficient as common markers.

Interestingly, bulk RNA expression was found to be highly correlated to the cell type more than the method of inducing senescence. In other words, regardless of the senescence inducer, such as telomere shortening, ionizing radiation, doxorubicin treatment, and the epxression of the HRAS G12V oncogene, the endothelial cells (HUVECs, HAECs) and fibroblasts (WI-38, IMR-90) are grouped, respectively (Casella et al., 2019) . This observation suggests that cellular senescence is likely cell-type-specific. Furthermore, even within the same cell type, a transcriptomic feature related to the origin of cell types is more prominent than the senescence-related transcriptomic changes. These results suggest that the in vivo transcriptomic signature of cellular senescence, i.e., the occurrence of senescence in specific cell types in specific tissues, may be more intricate and diverse depending on the spatial diversity. Two main factors likely lead to the in vivo heterogeneity of senescent cells. First, during the RNA-sequencing analysis of gene expression changes in various mouse organs, most aging-related changes were grouped by organ (Zhou et al., 2020) . However, the different organs have varying epigenetic programs and microenvironments and undergo diverging aging processes (Ashapkin et al., 2017) . As a result, the senescence in an organ appears heterogeneous.

Second, senescent cells are not static; their transcriptomic and epigenomic states may change over time (Fig. 2) . By tracking transcriptomic changes over time, the senescence induced by ionizing radiation was divided into three stages: early senescence, after 4 days of senescence induction; intermediate senescence, after 10 days of senescence induction; and late senescence, after 20 days of senescence induction. The early, intermediate, and late stages were enriched for transcript related to DNA damage responses, signaling mediated by p38-MAPK, and cell cycle arrest, respectively. In addition, the SASP-encoding genes also significantly changed 

The composition of the senescent cell is different for each organ. Due to aging, heterogeneity between cells increases, and some cells become senescent. Although senescent cells are a rare population in tissue, SASP secreted by senescent cells not only affects neighboring cells but also travels through the blood to affect other tissues. The percentage of senescent cells is high in adipose tissue and low in the lung. SAT, subcutaneous adipose tissue.

over time (Hernandez-Segura et al., 2017) . Therefore, studies that consider the effects of temporal and spatial context are required to understand in vivo cellular senescence. Since senescent cells rarely exist in organs, it has been challenging to track and analyze their transcription signatures in vivo (Fig. 2) (Biran et al., 2017) . The SASPs secreted by senescent cells vary depending on the senescence induction stimuli and cell types (Basisty et al., 2020) . In vivo SASPs mediate senescent cell clearance by recruiting immune cells, maintaining tissue integrity, i.e., preventing excessive fibrosis by ECM degradation enzymes, and reprogramming adjacent cells through an NFκB-driven secretome like IL-6 (Chiche et al., 2017; Kang et al., 2011; Krizhanovsky et al., 2008; Mosteiro et al., 2016; Xue et al., 2007) . Also, the persistence of the effects of SASPs could create a chronic inflammatory microenvironment that is favorable to cancer formation and senescence induction (Chen et al., 2005; Coppé et al., 2006; Perez-Mancera et al., 2014) . Moreover, SASPs can travel through the blood and systemically drive aging-related dysfunctions (Fig. 2) (Baar et al., 2017; Childs et al., 2016; Wiley et al., 2019) . Therefore, the contents of SASPs are likely more diverse in vivo than in vitro, adding to the complexity in analyzing senescent cells in vivo. It will be difficult to fully explain the complexity of in vivo senescent cells with only in vitro studies.

Studies based on bulk RNA-sequencing have been used to track the longitudinal changes in transcriptomes in various mice tissues and expand the understanding of in vivo aging. For example, (Schaum et al., 2019) tracked aging-related transcriptomic changes in 17 murine organs. At 12 months (mid-life), the transcripts in specific organs, such as subcutaneous and gonadal adipose tissues, changed more rapidly than in other organs. The changes in adipose tissue during aging can lead to inflammation in various tissues (Palmer and Kirkland, 2016) . Indeed, genes related to the immune response continued to be activated throughout the organs from mid-life. The temporal clustering of the changes in genes showed that genes related to the extracellular matrix, mitochondrion, and protein folding were decreased with age, and the genes related to stress and immune responses were increased with age. However, these changes were not commonly detected in whole organs. In particular, no changes in the immune response genes were detected in bones, the pancreas, or the skin. No changes in genes related to protein folding were detected in white blood cells, the marrow, and the intestine. Since various transcriptomic changes, such as those in genes related to the extracellular matrix, protein folding, and immune response, were highly enriched in adipose tissues, adipose tissue is likely not only a source of inflammation but also an organ undergoing the most changes during aging.

It is unclear whether the transcriptomic changes in each organ during aging are due to a cell-intrinsic change or a change in cellular composition. Indeed, tracking changes in the cell populations, such as the accumulation of senescent cells, the infiltration of immune cells, and the reduction of stem cells, which can significantly influence changes in aging transcriptome, are challenging in bulk RNA-sequencing. Single-cell RNA-sequencing (scRNA-seq) technology, which can identify intra-organ and inter-organ heterogeneity by tracking changes in the transcriptome of each individual cell in the aging process, could be applied to aging research (Fig. 3) .

Unlike bulk RNA-seq, scRNA-seq needs to distinguish every single cell and molecule from a mixed sample. Most scRNAseq methods enrich the mRNA using poly(T) oligonucleotides with barcodes that can tag distinct cells and molecules. Specifically, the cell barcode identifies each cell to facilitate multiplexing. The unique molecule identifier (UMI) is a barcode that Fig. 3 . The single-cell RNA-sequencing of human tissue. scRNA-seq has been performed only in a limited number of human tissues. Due to the morphology and rarity of senescent cells, senescence can be underestimated by scRNA-seq. Since senescent cells may be lost occur during organ dissociation and library construction, it is necessary to enrich senescent cells or increase throughput to understand senescent cells thoroughly.

can identify each molecule, thereby correcting amplification biases. Meanwhile, template switching (TS) is used in the reverse transcription step to obtain full-length transcripts. When a cDNA molecule reaches the 5′-end of an mRNA molecule, a few nucleotides are added to the 3′-end of the cDNA by a specific reverse transcriptase. Then, a unique TS oligo binds to this site and uses it as a template, so that cDNA synthesis continues to the 5′-end of TS oligo. Afterward, the cDNA molecule is amplified by polymerase chain reaction or in vitro transcription and can be subject to the full-length or 3′-or 5′ -end mRNA sequencing.

Various microfluidic-and plate-based methods have been rapidly developed based on this process. For example, platebased methods, such as Smart-seq2, STRT-seq, and Quartzseq, have the advantage of higher cell capture rate and are applied in full-length or 3′-or 5′-end transcript sequencing. On the other hand, InDrops and Drop-seq, which are microfluidic system-based on single-cell encapsulation in droplets, can easily scale up the number of cells, thus increasing the possibility of detecting rare cells (Lafzi et al., 2018) . Since each method has strengths and weaknesses, the appropriate method should be selected while considering the cell type, costs, efficiency, and scalability.

Recently, the Tabula Muris Senis (Mouse Aging Cell Atlas), which analyzes the aging in mice with scRNA-seq, has been released (Almanzar et al., 2020) (Table 3 ). In this study, 23 organs from 19 male and 11 female mice have been analyzed at 6 time points. The study observed the appearance of aging-specific cell types in various tissues, such as the liver and bladder, and the immune infiltration of various tissues. In addition, the number of T and B cells in gonadal adipose tissues increased with age, and the B cell cluster highly expressing immunoglobulin J was found specifically in old mice. Therefore, scRNA-seq enables the monitoring of the changes in cell population and the discovery of new rare populations associated with aging.

However, scRNA-seq has difficulty detecting senescent cells, which account for only 2% of specific in vivo tissues (Biran et al., 2017) . In the Tabula Muris Senis dataset, the p16-expressing cell fraction at 24 months was doubled of that in 3 months, and the expression level of p16 also increased by approximately 2-fold. However, according to a murine scRNA-seq study of 3 organs at 2 time points con-ducted by Calico, the fraction of p16-expressing cells and the average p16 expression level did not increase during aging (Kimmel et al., 2019) (Table 3 ). This inconsistency suggests that p16-expressing senescent cells are quite rare and are likely to be lost during tissue dissociation and sequencing library construction. Alternatively, the accumulation pattern of p16-expressing senescent cells may differ between individuals.

Due to ethical issues, scRNA-seq analysis has been performed only in a limited number of human tissues (Table 4 and references therein). In human lung tissues obtained from donors and pulmonary fibrosis patients, rare cell populations, such as senescent cells, appeared when pulmonary fibrosis progressed with age (Reyfman et al., 2019) . Interestingly, senescence markers such as CDKN2A, IL6, SERPINE1, and GLB1, were not specifically expressed, but the transcriptomic signature of replicative senescence was enriched in senescent cells. These data suggest that in vivo, human senescent cells may have a different transcriptomic signature from in vitro senescent cells. Therefore, more in vivo longitudinal studies are needed to elucidate the transcriptomic signature of senescent cells and their effects on the surrounding tissues and organismal aging.

Senescent cells, which have heterogeneous transcriptomic signature, are accumulated with aging. Moreover, the senescent cells can affect the surrounding cells and the whole body by SASPs. We need to understand not only the features but also the ecosystem of the senescent cells. However, the throughput of scRNA-seq, which has only captured approximately 10,000 cells to date, needs to be improved to capture the rare population of senescent cells thoroughly (Fig. 3) . Also, a technology that enriches senescent cells needs to be developed. In addition, it will be necessary to establish an organ-specific, in vivo senescent cell marker gene other than p16, which may not be a suitable marker for every context. Furthermore, the degree of heterogeneity of senescent cells in vivo and the effect of SASPs on cells located in neighbor and other tissues remain elusive. Utilizing single-cell spatial transcriptomics will help elucidate the origin of the senescent cells and the range of cells are affected by SASP. in dermal fibroblast, promote senescence phenotypes including increased SA-β-gal-positive cells and increased inflammation.

*Specific ages are not defined.

( Fig. 2 ). In addition, utilizing single-cell spatial transcriptomics to study various tissues longitudinally can help expand the understanding of the systemic effects of senescent cells (Roy et al., 2020) . A detailed understanding of senescent cells in vivo will enable the discovery of new targets that can be used to specifically eliminate senescent cells and the establishment of new strategies to reverse aging.

Chemokine signaling via the CXCR2 receptor reinforces senescence

A single-cell transcriptomic atlas characterizes ageing tissues in the mouse

Senolytic CAR T cells reverse senescence-associated pathologies

Aging as an epigenetic phenomenon

Targeted apoptosis of senescent cells restores tissue homeostasis in response to chemotoxicity and aging

A proteomic atlas of senescence-associated secretomes for aging biomarker development

Alteration of pRb/p16/cdk4 regulation in human osteosarcoma

Quantitative identification of senescent cells in aging and disease

Social Life" of senescent cells: what is SASP and why study it?

Monitoring tumorigenesis and senescence in vivo with a p16 INK4a -luciferase model

Cellular senescence as a tumor-suppressor mechanism

Transcriptome signature of cellular senescence

Transcriptional regulation of bcl-2 by nuclear factor κB and its significance in prostate cancer

Epithelial cell senescence induces pulmonary fibrosis through Nanog-mediated fibroblast activation

Crucial role of p53-dependent cellular senescence in suppression of Pten-deficient tumorigenesis

Injury-induced senescence enables in vivo reprogramming in skeletal muscle

Senescence and apoptosis: dueling or complementary cell fates?

Senescent intimal foam cells are deleterious at all stages of atherosclerosis

Secretion of vascular endothelial growth factor by primary human fibroblasts at senescence

Senescence-associated secretory phenotypes reveal cell-nonautonomous functions of oncogenic RAS and the p53 tumor suppressor

An essential role for senescent cells in optimal wound healing through secretion of PDGF-AA

A reference single-cell transcriptomic atlas of human skeletal muscle tissue reveals bifurcated muscle stem cell populations

A biomarker that identifies senescent human cells in culture and in aging skin in vivo

Synthetic lethal metabolic targeting of cellular senescence in cancer therapy

ATF6α regulates morphological changes associated with senescence in human fibroblasts

Single-cell analysis of human pancreas reveals transcriptional signatures of aging and somatic mutation patterns

Involvement of nuclear factor-kappa B in bcl-xL-induced interleukin 8 expression in glioblastoma

p16 (Ink4a) and senescence-associated β-galactosidase can be induced in macrophages as part of a reversible response to physiological stimuli

Single-cell RNA sequencing of microglia throughout the mouse lifespan and in the injured brain reveals complex cell-state changes

Regulation of p16 CDKN2 expression and its implications for cell immortalization and senescence

Unmasking transcriptional heterogeneity in senescent cells

Survey of senescent cell markers with age in human tissues

Accumulation of senescent cells in mitotic tissue of aging primates

Senescence surveillance of pre-malignant hepatocytes limits liver cancer development

Murine single-cell RNA-seq reveals cellidentity-and tissue-specific trajectories of aging

Ink4a/Arf expression is a biomarker of aging

Senescence of activated stellate cells limits liver fibrosis

Senescent fibroblasts promote epithelial cell growth and tumorigenesis: a link between cancer and aging

Tutorial: guidelines for the experimental design of single-cell RNA sequencing studies

BH4 domain of Bcl-2 as a novel target for cancer therapy

The chemistry of senescence

Long-term quiescent fibroblast cells transit into senescence

Tissue damage and senescence provide critical signals for cellular reprogramming in vivo

Cellular senescence: from physiology to pathology

Nonresolving inflammation

Targeting senescent cells in translational medicine

Fat apoptosis through targeted activation of caspase 8: a new mouse model of inducible and reversible lipoatrophy

Aging and adipose tissue: potential interventions for diabetes and regenerative medicine

Inside and out: the activities of senescence in cancer

Immune senescence, epigenetics and autoimmunity

p16 INK4A is a robust in vivo biomarker of cellular aging in human skin

Single-cell transcriptomic analysis of human lung provides insights into the pathobiology of pulmonary fibrosis

A blueprint for characterizing senescence

Ageing hallmarks exhibit organ-specific temporal signatures

A new regulatory motif in cell-cycle control causing specific inhibition of cyclin D/CDK4

Forging a signature of in vivo senescence

Different telomere damage signaling pathways in human and mouse cells

Aging and DNA damage in humans: a meta-analysis study

Senescence is a developmental mechanism that contributes to embryonic growth and patterning

SILAC analysis reveals increased secretion of hemostasis-related factors by senescent cells

Targeting senescent cells enhances adipogenesis and metabolic function in old age

Senescence and tumour clearance is triggered by p53 restoration in murine liver carcinomas

The limitations and validities of senescence associated-β-galactosidase activity as an aging marker for human foreskin fibroblast Hs68 cells

A single-cell transcriptome atlas of the aging human and macaque retina

The four layers of aging

A human circulating immune cell landscape in aging and COVID-19

A landscape of murine long non-coding RNAs reveals the leading transcriptome alterations in adipose tissue during aging

Cellular senescence and the senescent secretory phenotype in age-related chronic diseases

Identification of a novel senolytic agent, navitoclax, targeting the Bcl-2 family of anti-apoptotic factors

The Achilles' heel of senescent cells: from transcriptome to senolytic drugs

A single-cell transcriptomic atlas of human skin aging

S.K. and C.K. wrote the manuscript.

The authors have no potential conflicts of interest to disclose.

https://orcid.org/0000-0002-1581-6729 Chuna Kim https://orcid.org/0000-0003-2337-1106