• A comprehensive transcriptome analysis of aged murine HSC identifies an aging signature.

  • The HSC aging signature is highly enriched for cell membrane–related transcripts and identifies age-associated heterogeneity.

We surveyed 16 published and unpublished data sets to determine whether a consistent pattern of transcriptional deregulation in aging murine hematopoietic stem cells (HSC) exists. Despite substantial heterogeneity between individual studies, we uncovered a core and robust HSC aging signature. We detected increased transcriptional activation in aged HSCs, further confirmed by chromatin accessibility analysis. Unexpectedly, using 2 independent computational approaches, we established that deregulated aging genes consist largely of membrane-associated transcripts, including many cell surface molecules previously not associated with HSC biology. We show that Selp (P-selectin), the most consistent deregulated gene, is not merely a marker for aged HSCs but is associated with HSC functional decline. Additionally, single-cell transcriptomics analysis revealed increased heterogeneity of the aged HSC pool. We identify the presence of transcriptionally “young-like” HSCs in aged bone marrow. We share our results as an online resource and demonstrate its utility by confirming that exposure to sympathomimetics or deletion of Dnmt3a/b molecularly resembles HSC rejuvenation or aging, respectively.

Aging is a natural and time-dependent process, leading to the functional decline of many cell types and tissues, including those of the hematopoietic system. Aged hematopoietic stem cells (HSCs) have less regenerative potential, and their overall contribution to formation of mature blood populations, at the clonal level, is restricted and lineage skewed when compared with their young counterparts, suggesting stem cell exhaustion, one of the hallmarks of aging.1-9 Transplantation of aged HSCs into young recipients does not improve HSC function,3,5  indicating that cell-intrinsic aspects play a pivotal role in HSC aging. The differential expression of genes is a prominent readout of multiple cell-intrinsic factors, and various studies have compared the transcriptomes of young and aged HSCs.10-24 These studies have suggested a plethora of distinct aging-associated mechanisms, including aberrant regulation of genes involved in lymphoid and platelet differentiation14,22 ; elevated inflammation and stress responses12,17 ; chromatin remodeling/DNA methylation12 ,23-25; changes in cell cycle and maintenance of quiescence15,16,19 ; and age-associated replicative stress in HSCs.13  It remains unclear whether all suggested factors equally contribute to HSC aging and why different studies report different mechanisms.

Although a substantial number of HSC aging transcriptome studies has been published, no cross-validation of all transcriptomic HSC aging data has yet been carried out. In fact, it is unknown whether a common HSC aging transcriptome exists, and if so, which genes are involved. In order to address this question, we selected 15 studies in which transcriptomes of young and aged murine HSCs were compared and added to these our own recently published transcriptome analysis.26 

Despite uncovering considerable heterogeneity among different studies, we were able to identify a robust and biologically relevant HSC transcriptomic aging signature (AS). This signature is mainly composed of genes that encode for membrane-associated proteins, many with unknown functions in HSC biology. A large fraction of AS genes are upregulated in aged HSCs, indicating transcriptional activation in this population. We also show that aging increases transcriptional heterogeneity among single HSCs and provide evidence that transcriptionally “young” HSCs exist in aged mice. Furthermore, we demonstrate that the top age-associated gene, Selp (P-selectin), contributes to functional decline of HSCs. Finally, we provide an open Web-based resource to the community and use independent data sets to demonstrate how our approach can be used to validate novel biological paradigms.

Mice

All experiments were approved by the Central Commission for Animal Testing and the Animal Ethical Committee. Young (4-5 months) and aged (20-24 months) C57BL/6, C57BL/6.SJL (10-14 weeks), and C57BL/6J-kitW-41J/kitW-41J (W41) mice (originally from European division of Harlan) were obtained from the Central Animal Facility at the University Medical Center Groningen. TetOP-H2B-GFP mice (provided by Rob Coppes, University Medical Center Groningen [UMCG], Groningen, The Netherlands). H2B-GFP expression was induced by providing mice (6 weeks) with food containing doxycycline for 6 weeks. Mice were then fed with normal food for 14 weeks. Mice were kept in group cages and housed under conventional temperature- and day cycle–controlled conditions.

Flow cytometry

Bone marrow was isolated from the tibia, femur, pelvis, sternum, and spine by crushing, and red blood cells were lysed with erylysis buffer. All samples were analyzed on a BD FACSCanto II (BD Biosciences) and sorted on MoFlo Astrios or XDP cell sorters (Beckman Coulter). Detailed information of the antibodies used for isolation and identifications can be found in the supplemental Data (available on the Blood Web site).

Retroviral constructs and transduction

Selp was amplified from C57BL/6J long-term HSC (LT-HSC) complementary DNA and then cloned into pSF912–IRES–GFP vector (originally as pSF91 from C. Baum, Hannover Medical School, Hannover, Germany, further modified in our laboratory). Virus was produced by transfecting Platinum-E retroviral packaging cells (Cell Biolabs) according to the manufacturer’s protocol. Viral supernatant was harvested 48 hours after transfection. LT-HSCs were transduced as previously described.6  Cells were harvested 20 to 22 hours after transduction for transplantation.

BM transplantation

For SELPlow and SELPhigh transplantation, 1000 SELPlow and SELPhigh LT-HSCs were isolated from an old CD45.2+ C57BL/6 mouse and transplanted together with 2 million CD45.1+ W41 bone marrow into lethally irradiated (9 Gy) CD45.1+ C57BL/6.SJL recipients, respectively. For Selp overexpression transplantation, empty vector (SF91–IRES–GFP empty) and SELP OE CD45.2+ cells harvested after retroviral transduction were transplanted together with 2 million CD45.1+ W41 bone marrow into lethally irradiated CD45.1+ C57BL/6.SJL recipients, respectively.

Assay for transposase-accessible chromatin using sequencing (ATAC-seq)

A total of 2.5 × 104 LT-HSCs were sorted from young and aged mice in biological replicates. Library preparation was performed according to Buenrostro et al.27  In brief, cells were lysed and tagmented for 1 hour at 37°C. Samples were then polymerase chain reaction amplified using provided adapter primers for 11 or 12 cycles. Following, samples were size selected for fragments of 150 to 1000 bp. Samples were pooled and sequenced to obtain 30 to 50 × 106 reads using a NextSeq 500. Data are deposited under accession number GSE166674.

Immunofluorescence

Images were acquired on a Leica Sp8 confocal microscope and quantified on Fiji ImageJ.

RNA content

A total of 4000 to 6000 LT-HSCs were seeded onto spots on an adhesion immunofluorescent slide (VWR). For RNA staining, cells were fixed with 100% methanol for 10 minutes and washed with phosphate-buffered saline. Cells were subsequently stained with 1:10 000 SYTO RNA Select Green Fluorescent Stain (Thermo Fisher Scientific) for 20 minutes at room temperature. After washing, coverslips were mounted with ProLong Diamond Antifade Mountant with 4′,6-diamidino-2-phenylindole (DAPI).

RNA polymerase II

For RNA polymerase II staining, cells were fixed and permeabilized with Fixation/Permeabilization Solution Buffer (BD Bioscience) for 20 minutes on ice. After washing, cells were blocked first with Endogenous Biotin-Blocking Kit (Thermo Fisher Scientific) followed by 4% bovine serum albumin for 30 minutes. Cells were stained with 1:100 biotin mouse monoclonal RNA polymerase II antibody (Novus Bio) at 4°C overnight. Cells were then washed with 0.1% Triton-X-100 and stained with 1:1000 streptavidin Alexa Fluor 488 secondary antibody. After washing, coverslips were mounted with ProLong Diamond Antifade Mountant with DAPI.

Further information

More detailed information about the methods used can be found in the supplemental Data.

Additional resources

Additional resources include the following: AS Web tool site (https://eriba.umcg.nl/agingsignature) and GitHub with all resources (https://github.com/LeonidBystrykh/data-for-manuscript).

Systematic collection of publications reporting gene expression changes in aged HSCs

We took advantage of PubMed’s API E-utility to select all informative studies on HSC aging in C57BL/6 mice (Figure 1A; supplemental Data). This search resulted in 16 published studies that met all criteria and included our own recently published data set26  (Figure 1B; supplemental Table 1). The various studies used different sequencing platforms, ranging from microarray hybridization, bulk RNA sequencing to single-cell RNA sequencing (scRNA-seq) (Figure 1B; supplemental Table 1) and included mainly 2 age groups (Figure 1B; supplemental Figure 1A). We employed 2 alternative approaches to perform the analysis. First, we compared the self-reported, preanalyzed lists of DEGs from the selected manuscripts and referred to this as “meta-analysis.” In parallel, we reanalyzed all available raw data in the most uniform way; these are referred to as “reanalysis” (Figure 1A; supplemental Table 1).

Figure 1.

Systematic analysis reveals heterogeneity in transcriptomic studies of young and aged murine HSCs. (A) Selection of young and aged HSC transcriptomic studies. Representation of decision tree used to acquire and analyze the different publications used for the systematic analysis. (B) Overview of the age groups used in each of the selected studies. The age group of mice used in each study (y-axis) is represented in months (x-axis). Color codes refer to the transcriptome platform used. Squares next to each study indicate in which kind of analysis each study was included. (C) Heterogeneity across different transcriptomic studies. Number of differentially expressed (DE) genes (DEGs) reported by each study. Left panel shows the number of DEGs as reported by the authors in their respective manuscripts. Right panel shows the number of DEGs identified by us using a uniform reanalysis. Red bars indicate the number of genes that are upregulated upon aging, while blue bars indicate downregulated genes. In the reanalysis panel, “h-Chambers” indicates half Chambers, meaning that just half of the DEGs were taken into account. (D) Ranking of DEGs across different studies. Horizontal bars represent the number of DEGs (x-axis, the actual number of DEG genes are shown in each bar), ranked by the number of studies in which they were found (occurrence in studies [consistency], y-axis). Orange and green bars indicate meta-analysis and reanalysis, respectively. (E) The list of genes gets more consistent with the increase number of studies added. Simulation of different list assembly using combinations of subsets of 12 studies (gradient color from blue to red). As the number of studies increased, the coefficient of variation (y-axis) decreases. Dots represent the different values and mean consistency for individual genes (see supplemental Data for further details). (F) Predicted resolution to discover consistently reported genes. The model demonstrates that >90% of all HSC aging genes have been identified (y-axis) in the 2 analytic approaches. Purple and green dashed lines indicate meta-analysis and reanalysis, respectively. (G) Extent of overlap between all independent studies and the HSC AS (centered red circle, AS). The size of the circles represents the number of DEGs identified in each study, the number of overlapping DEGs with the AS is represented next to each study identification, and the distance from the center represents the degree of similarity of each study to the HSC AS. Studies that show a higher degree of similarity are closer to the AS.

Figure 1.

Systematic analysis reveals heterogeneity in transcriptomic studies of young and aged murine HSCs. (A) Selection of young and aged HSC transcriptomic studies. Representation of decision tree used to acquire and analyze the different publications used for the systematic analysis. (B) Overview of the age groups used in each of the selected studies. The age group of mice used in each study (y-axis) is represented in months (x-axis). Color codes refer to the transcriptome platform used. Squares next to each study indicate in which kind of analysis each study was included. (C) Heterogeneity across different transcriptomic studies. Number of differentially expressed (DE) genes (DEGs) reported by each study. Left panel shows the number of DEGs as reported by the authors in their respective manuscripts. Right panel shows the number of DEGs identified by us using a uniform reanalysis. Red bars indicate the number of genes that are upregulated upon aging, while blue bars indicate downregulated genes. In the reanalysis panel, “h-Chambers” indicates half Chambers, meaning that just half of the DEGs were taken into account. (D) Ranking of DEGs across different studies. Horizontal bars represent the number of DEGs (x-axis, the actual number of DEG genes are shown in each bar), ranked by the number of studies in which they were found (occurrence in studies [consistency], y-axis). Orange and green bars indicate meta-analysis and reanalysis, respectively. (E) The list of genes gets more consistent with the increase number of studies added. Simulation of different list assembly using combinations of subsets of 12 studies (gradient color from blue to red). As the number of studies increased, the coefficient of variation (y-axis) decreases. Dots represent the different values and mean consistency for individual genes (see supplemental Data for further details). (F) Predicted resolution to discover consistently reported genes. The model demonstrates that >90% of all HSC aging genes have been identified (y-axis) in the 2 analytic approaches. Purple and green dashed lines indicate meta-analysis and reanalysis, respectively. (G) Extent of overlap between all independent studies and the HSC AS (centered red circle, AS). The size of the circles represents the number of DEGs identified in each study, the number of overlapping DEGs with the AS is represented next to each study identification, and the distance from the center represents the degree of similarity of each study to the HSC AS. Studies that show a higher degree of similarity are closer to the AS.

Close modal
Table 1.

Top 10 AS genes with experimental validated function in HSC aging

Gene symbolConsistency
(out of 12)
Average log2FCFunctionPhenotypeReferences
Itgb3 1.25 ± 0.48 Cell adhesion; signal transduction Itgb3high identifies myeloid-biased HSCs 17,39 
Vwf 1.37 ± 0.66 Cell adhesion; signal transduction Vwfhigh identifies myeloid-biased HSCs 14,40,41 
Neo1 1.84 ± 1.18 Cell adhesion; signal transduction Neo1high identifies more dormant HSCs 26,42 
Alcam 1.39 ± 0.62 Cell adhesion Alcam+ identifies myeloid-biased HSCs and with reduced potential 43 
CD86 −1.39 ± 0.36 Signal transduction Required for lymphoid differentiation 44 
Gadd45g 1.54 ± 0.99 Stress response Gadd45g/ shows increased self-renewal and deficient megakaryocytic differentiation 45 
Satb1 −1.21 ± 0.50 Epigenetic effector Required for lymphoid differentiation 46 
Sdpr 1.88 ± 1.61 Transcription regulation Involved in PRC2 expression regulation 47 
Slamf1 0.92 ± 0.16 Cell adhesion; signal transduction Slamf1high identifies myeloid-biased HSCs with reduced potential 48,49 
Dnmt3b −1.09 ± 0.62 Epigenetic regulator Dnmt3a/b−/− shows differentiation blockage/increased self-renewal 29,50,51 
Gene symbolConsistency
(out of 12)
Average log2FCFunctionPhenotypeReferences
Itgb3 1.25 ± 0.48 Cell adhesion; signal transduction Itgb3high identifies myeloid-biased HSCs 17,39 
Vwf 1.37 ± 0.66 Cell adhesion; signal transduction Vwfhigh identifies myeloid-biased HSCs 14,40,41 
Neo1 1.84 ± 1.18 Cell adhesion; signal transduction Neo1high identifies more dormant HSCs 26,42 
Alcam 1.39 ± 0.62 Cell adhesion Alcam+ identifies myeloid-biased HSCs and with reduced potential 43 
CD86 −1.39 ± 0.36 Signal transduction Required for lymphoid differentiation 44 
Gadd45g 1.54 ± 0.99 Stress response Gadd45g/ shows increased self-renewal and deficient megakaryocytic differentiation 45 
Satb1 −1.21 ± 0.50 Epigenetic effector Required for lymphoid differentiation 46 
Sdpr 1.88 ± 1.61 Transcription regulation Involved in PRC2 expression regulation 47 
Slamf1 0.92 ± 0.16 Cell adhesion; signal transduction Slamf1high identifies myeloid-biased HSCs with reduced potential 48,49 
Dnmt3b −1.09 ± 0.62 Epigenetic regulator Dnmt3a/b−/− shows differentiation blockage/increased self-renewal 29,50,51 

Individual studies identify highly variable numbers of DEGs

Strikingly, whereas some studies reported thousands of genes to be differentially expressed, others only reported dozens (Figure 1C). When DEGs from the meta-analysis and reanalysis approaches were compared, both demonstrated a significant correlation between fold-change (FC) expression values, indicating that reanalysis of the raw data did not deviate drastically from the originally reported data sets (supplemental Figure 1C). The observed study-to-study variations might result from differences in mouse gender, group size, exact age, or HSC immunophenotyping; intrinsic differences between expression platforms; or unknown environmental effects, including differences in mouse facilities or breeding companies from which the mice originated (supplemental Table 1). Principal-component analysis performed in the reanalyzed sets demonstrated that studies clustered according to their sequencing platform (supplemental Figure 1B), although degree of their separation was minor. Overall, this analysis demonstrates that the variability between studies partially has a technical origin, but the discrepancy between studies is not fully explained.

Most reported age-associated HSC genes have low reproducibility

Collectively, these 16 studies reported >6000 genes being differentially expressed between young and aged HSCs. Surprisingly, the vast majority of reported genes demonstrated poor reproducibility, as almost 80% of all putative age-associated genes were reported in a single study only (supplemental Figure 1D). Moreover, pairwise analysis of all different studies showed that some studies had little to no overlapping genes with any other study (supplemental Figure 1E-F).

Identification of an HSC aging signature (AS)

Despite limited reproducibility between individual data sets, we wondered whether it would be possible to identify consistently DEGs across all studies. To this end, we ranked all reported genes by their occurrence in different studies (consistency) (Figure 1D). This revealed a limited number of highly consistent genes that were frequently documented to be differentially expressed upon HSC aging. The single most consistently deregulated gene was Selp. We provide a complete gene list in supplemental Table 2. We restricted further analysis to genes that were consistently reported in ≥4 studies, resulting in a list of 142 genes for the meta-analysis and 220 genes for the reanalyzed data sets. We refer to this as the “aging signature” (supplemental Figure 1G; supplemental Table 2).

The AS is robust and provides higher resolution than individual studies

In order to assess the minimal number of publications required to generate a robust list of consistent aging genes, we modeled all possible combinations of most consistent genes from the reanalyzed data, starting from 1 data set at a time to all combinations of 11 out of 12 data sets (Figure 1E). This analysis generated multiple consistencies for the same genes across different permutations. When the analysis was performed with a small number of data sets, we observed a high degree of variation, which started to diminish when the number of data sets increased. When >6 data sets were used, the number and consistency of genes stabilized (Figure 1E; supplemental Data). Complementarily, we also modeled the predicted resolution (the number of aging genes found) for the same range of studies (Figure 1F). The model predicted that the number of studies that we used for our analyses (11 for meta-analysis and 12 for reanalysis) have high resolution; with the current number of included studies, >90% of all the possible age-associated genes have been found. This also implies that the gene collection in the HSC AS remains stable and that the addition of novel HSC transcriptomes is unlikely to significantly alter the AS.

Figure 2.

The HSC AS is highly enriched for membrane-associated genes. (A) Classification of the proteins-coding AS genes by cellular compartment. The localization of each protein was determined using UniProt classification. Colors of the bars represent upregulation (red) or repression (blue) upon aging of AS genes. Numbers next to the bars indicate the number of genes found in each category. (B) Network representation of the aging genes based on the cellular compartment where the encoded proteins are predicted to be localized. The size of each circle represents the number of genes found in each category. Connections between each cellular compartment reflect the number of genes which are shared per cellular compartment. (C-D) GO analysis of biological process for AS genes. Horizontal bar represent individual GO terms (y-axis) and their log-transformed false discovery rate (FDR) values (x-axis). Numbers in each bar demonstrate the number of genes found per GO term. Red (C) and blue (D) bars represent GO terms found to be significantly enriched in upregulated and downregulated AS genes, respectively. GO terms are selected based on an FDR <0.05 cutoff.

Figure 2.

The HSC AS is highly enriched for membrane-associated genes. (A) Classification of the proteins-coding AS genes by cellular compartment. The localization of each protein was determined using UniProt classification. Colors of the bars represent upregulation (red) or repression (blue) upon aging of AS genes. Numbers next to the bars indicate the number of genes found in each category. (B) Network representation of the aging genes based on the cellular compartment where the encoded proteins are predicted to be localized. The size of each circle represents the number of genes found in each category. Connections between each cellular compartment reflect the number of genes which are shared per cellular compartment. (C-D) GO analysis of biological process for AS genes. Horizontal bar represent individual GO terms (y-axis) and their log-transformed false discovery rate (FDR) values (x-axis). Numbers in each bar demonstrate the number of genes found per GO term. Red (C) and blue (D) bars represent GO terms found to be significantly enriched in upregulated and downregulated AS genes, respectively. GO terms are selected based on an FDR <0.05 cutoff.

Close modal

We evaluated the degree of similarity between the AS and the individual studies from which it is derived (Figure 1G). All studies demonstrated approximately the same similarity to the AS, with no clear outliers. No single study encompassed all aging genes, but collectively, they contribute to the robustness of the AS. We conclude that the AS can be used as a reference for DEGs upon HSC aging.

The AS reveals high enrichment for membrane-related processes in HSC transcriptional activation

We restricted subsequent analyses to the AS derived from the reanalyzed data and observed several remarkable features. First, approximately half of the aging genes encode for membrane-associated proteins (Figure 2A-B). These included many cell adhesion molecules or channels (Selp, Alcam, Cd34, Dsg2, Itgb3, and Clca3a1) and cytokine receptors (Csf2rb, Ebi3, Flt3, Ghr, Il1rapl2, and Osmr), many of which have not been previously associated with HSC biology.

Gene Ontology (GO) analysis identified a total of 46 enriched GO function terms (supplemental Table 3). As expected, the most significant GO terms revolve around membrane-associated functions (supplemental Figure 2A). Among upregulated processes, we observed significant enrichment of different GO terms associated with cell adhesion and activation and the organization of cell junction and projection (Figure 2C). Within the downregulated processes, we identified significant enrichment for DNA replication, cytokine signaling, tissue homeostasis, and cell proliferation (Figure 2D). Pathway analysis also demonstrated analogous pathways (supplemental Figure 2B).

We next collected all GO terms associated with every gene present in the AS (supplemental Figure 2C), resulting in >2000 GO terms. The majority of these terms were very specific and contained <20 AS genes (supplemental Figure 2C, lower panel). We then searched in the AS GO term database for key terms that represent previously described HSC aging mechanisms, including “inflammation,”12  “replicative stress,”13  “cell cycle,”15,16  and “epigenetic drift,”12,23,28  and indeed found these. However, the large majority of genes (149) were associated with “membrane” (supplemental Figure 2D), denoting that significant changes occur in the HSC membrane upon aging.

Secondly, we observed that a disproportional number (181 genes, ∼80%) from the AS were upregulated in aged HSCs (Figure 3A). We assessed whether this would be reflected in global transcriptional activation. Comparing young and aged LT-HSCs (LSK CD48 CD150+) (supplemental Figure 3A), we observed an approximately threefold higher overall RNA yield per cell in aged HSCs (supplemental Figure 3B). We confirmed this by intracellular staining of young and aged LT-HSCs with the RNA-specific dye SYTO (supplemental Figure 3C). In order to investigate if increased RNA abundance in aged LT-HSCs was directly linked to increased transcription, we performed intracellular staining for RNA polymerase II in young and aged HSCs and indeed observed increased RNA polymerase II activity in aged cells (Figure 3B).

Figure 3.

Skewing of HSC AS genes reflects transcriptional activation in aged HSCs. (A) Skewing of HSC AS genes. Volcano plot depicting expression levels of AS genes; 181 upregulated and 39 downregulated genes are represented. Dots represent the average of either adjusted P values (y-axis) and FCs (x-axis) per gene in different studies. (B) Protein levels of RNA polymerase II in individual young and aged LT-HSCs. Right panel shows the mean fluorescence intensity (MFI) quantification of young (n = 120 cells) and aged HSCs (n = 139 cells). Scale bars, 5 μm; mean ± standard deviation is shown. ****P < .0001. (C) Aged HSCs demonstrate more age-specific chromatin accessible sites than young HSCs. Volcano plot depicting differentially accessible sites (peaks) between young and aged HSCs measured by ATAC-seq. Dot colors demonstrate peaks either in genomic regions overlapping with coding (green) or noncoding (purple) annotations. Peaks with negative FC values are significantly more accessible in young HSCs, whereas peaks with positive FC values are more significantly accessible in aged HSCs. Gene symbols annotate peaks which overlap with gene bodies from some AS genes. (D) Transcriptional program is associated with ATAC-seq of young and aged HSCs. GO analysis of functional process for overlapping aged-specific ATAC-seq accessible sites with genes from the aging list. Horizontal bar represent individual GO terms (y-axis) and their FDR values (x-axis). Numbers in each bar demonstrate the number of genes found per GO term. GO terms are selected based on a FDR <0.05 cutoff.

Figure 3.

Skewing of HSC AS genes reflects transcriptional activation in aged HSCs. (A) Skewing of HSC AS genes. Volcano plot depicting expression levels of AS genes; 181 upregulated and 39 downregulated genes are represented. Dots represent the average of either adjusted P values (y-axis) and FCs (x-axis) per gene in different studies. (B) Protein levels of RNA polymerase II in individual young and aged LT-HSCs. Right panel shows the mean fluorescence intensity (MFI) quantification of young (n = 120 cells) and aged HSCs (n = 139 cells). Scale bars, 5 μm; mean ± standard deviation is shown. ****P < .0001. (C) Aged HSCs demonstrate more age-specific chromatin accessible sites than young HSCs. Volcano plot depicting differentially accessible sites (peaks) between young and aged HSCs measured by ATAC-seq. Dot colors demonstrate peaks either in genomic regions overlapping with coding (green) or noncoding (purple) annotations. Peaks with negative FC values are significantly more accessible in young HSCs, whereas peaks with positive FC values are more significantly accessible in aged HSCs. Gene symbols annotate peaks which overlap with gene bodies from some AS genes. (D) Transcriptional program is associated with ATAC-seq of young and aged HSCs. GO analysis of functional process for overlapping aged-specific ATAC-seq accessible sites with genes from the aging list. Horizontal bar represent individual GO terms (y-axis) and their FDR values (x-axis). Numbers in each bar demonstrate the number of genes found per GO term. GO terms are selected based on a FDR <0.05 cutoff.

Close modal

Finally, we evaluated the number of chromatin accessible sites by ATAC-seq. In agreement with the data above, aged HSCs showed an increased number of overall accessible sites compared with young HSCs (supplemental Figure 3D). In addition, differential analysis of accessible sites demonstrated strong skewing toward more open sites in aged HSCs (Figure 3C). We then speculated that aged-specific accessible sites could contribute to deregulated expression. To investigate this, we overlapped the age-restricted accessible sites with the genes in the aging list, which would reflect genes more accessible in aged HSCs and therefore upregulated, and performed a GO analysis (Figure 3D). This revealed high enrichment for transcription-related processes, further reinforcing our hypothesis that aged HSCs are transcriptionally activated.

Taken together, the AS suggests that aged HSCs undergo adaptation at cellular and functional levels through membrane-related processes and that aged HSCs demonstrate higher transcriptional activation.

The AS as a reference for age-associated studies

We aimed to test whether the AS can be used as a reference to (re)analyze existing and future HSC gene expression profiles. To this end, we developed 2 independent approaches where custom transcriptome sets can be compared with the AS.

First, we designed 2 scoring methods, the discrete gene set enrichment (dGSE) algorithm and a rank score (RS), both of which are able to cope with semiordered (discrete) lists of genes (supplemental Data) and can be used for testing bulk RNA-sequencing sets. To benchmark both methods, we individually tested the 12 HSC aging transcriptome studies that comprise the AS (Figure 4A). As expected, we observed that the different data sets had a wide range of enrichment scores when compared with the AS.

Figure 4.

Employing the AS as a bona fide transcriptomic reference. (A) Using RS and dGSE analyses to compare existing data sets to the AS. Left panel shows a graphical representation of the dGSE enrichment analyses from our own generated HSC transcriptome data compared with the AS. Right panel shows plots of RS (upper) and dGSE (lower) enrichment analyses of all individual studies that comprise the AS. Numbers refer to the studies shown on the right. (B) Searching for an AS using RS and dGSE analyses in independent data sets. Upper left panel shows a lack of enrichment for aging genes in an embryonic stem cell signature from Wong et al.30  Upper right panel shows enrichment for the top 250 genes (sorted by FDR value) of Kim et al (K.K., A.M.-M and B.v.E, manuscript submitted December 2020). Lower right panel shows significant enrichment for an HSC AS among genes downregulated in “rejuvenated” HSCs from Maryanovich et al.18  Lower right panel shows significant enrichment for a HSC signature among differentially expressed genes in HSCs isolated from young DNMT3A/B knockout mice from Challen et al.29  The right-most panels show combined RS and dGSE scores for each test set. (C) Single-cell AS scores. Separation of young (circles) and aged (triangles) single-cell transcriptomics and their respective AS score. From Grover scRNA-seq data (48 young and 68 aged cells), individual cells were separated by principal component 1 (PC1) (x-axis) and PC2 (y-axis) and color-coded according to their individual AS score. Left column displays only young cells; middle column shows just aged cells and right column shows both populations. Cells are colored score from low (blue) to high (yellow) AS scores. On the far-right panel, violin and boxplots demonstrate the overall AS score (y-axis) distribution of the young (gray) and aged (red) cells.

Figure 4.

Employing the AS as a bona fide transcriptomic reference. (A) Using RS and dGSE analyses to compare existing data sets to the AS. Left panel shows a graphical representation of the dGSE enrichment analyses from our own generated HSC transcriptome data compared with the AS. Right panel shows plots of RS (upper) and dGSE (lower) enrichment analyses of all individual studies that comprise the AS. Numbers refer to the studies shown on the right. (B) Searching for an AS using RS and dGSE analyses in independent data sets. Upper left panel shows a lack of enrichment for aging genes in an embryonic stem cell signature from Wong et al.30  Upper right panel shows enrichment for the top 250 genes (sorted by FDR value) of Kim et al (K.K., A.M.-M and B.v.E, manuscript submitted December 2020). Lower right panel shows significant enrichment for an HSC AS among genes downregulated in “rejuvenated” HSCs from Maryanovich et al.18  Lower right panel shows significant enrichment for a HSC signature among differentially expressed genes in HSCs isolated from young DNMT3A/B knockout mice from Challen et al.29  The right-most panels show combined RS and dGSE scores for each test set. (C) Single-cell AS scores. Separation of young (circles) and aged (triangles) single-cell transcriptomics and their respective AS score. From Grover scRNA-seq data (48 young and 68 aged cells), individual cells were separated by principal component 1 (PC1) (x-axis) and PC2 (y-axis) and color-coded according to their individual AS score. Left column displays only young cells; middle column shows just aged cells and right column shows both populations. Cells are colored score from low (blue) to high (yellow) AS scores. On the far-right panel, violin and boxplots demonstrate the overall AS score (y-axis) distribution of the young (gray) and aged (red) cells.

Close modal

We next used transcriptomic data sets from 3 unrelated studies. The first independent data set compared young and physiologically aged HSCs (K.K., A.M.-M. and B.v.E., manuscript submitted December 2020). The second involved HSCs that were “rejuvenated” by a sympathomimetic drug acting on adrenoreceptor β3.18  The third study used expression data from HSCs that were derived from Dmnt3a/b-deficient mice that show accelerated aging phenotypes.29  As a negative control, we included a transcriptomic signature of embryonic stem cells.30  As expected, the embryonic stem cell signature demonstrated low scores for aging genes. In contrast, the independent set of young and physiologically aged HSCs demonstrated enrichment when compared with the AS. When we compared the list of genes that were downregulated in “rejuvenated” old HSCs to the AS, we also found strong enrichment. The transcriptome of young HSCs isolated from Dmnt3a/b knockout mice also scored high for HSC aging genes. These data suggest that both β3 sympathomimetic signaling and loss of Dnmt3a/b molecularly phenocopy normal HSC aging, additionally confirming the value of our method (Figure 4B).

This approach allowed us to investigate the relationship of our HSC AS to a previously established multiorgan interspecies transcriptomic signature.31  Expectedly, we observed a mild enrichment between both signatures, denoting that some more general aging-associated genes are also present in HSC aging (supplemental Figure 4A). Genes present in both signatures are related to inflammatory response. We also investigated the involvement of transforming growth factor β (TGF-β) signaling, which has been implicated in HSC aging,21,23,32  in the AS. We compared different molecular signatures of TGF-β signaling33,34  to our AS, but none of the sets demonstrated high scores. Furthermore, only 5 TGF-β target genes from all sets (Bmpr1a, Id2, Jun, Itgb3, and Mef2c) are present in the AS. Thus, TGF-β signaling, from a transcriptional perspective, is not robustly involved in HSC aging (supplemental Figure 4A).

Interestingly, single-cell transplantations have revealed that the overall engraftment potential of aged HSCs is substantially lower than that of young HSCs.2  However, some aged HSCs display engraftment performance similar to young cells, suggesting that the aged HSC pool is functionally heterogenous. Thus, we wondered whether the AS can be used to also identify intrinsic variation between cells from the same population. To identify such “young-like” cells using single-cell transcriptomic studies, we assigned an AS score to individual cells (Figure 4C; supplemental Figure 4B). This revealed that separation of young and aged cells is in good agreement with their assigned AS scores, where aged cells score higher and young cells lower, also predicting their age. Interestingly, some aged cells had a low AS score and were close to the young cells cluster, suggesting that these cells are, at least transcriptionally, young-like.

Taken together, the HSC AS demonstrated to be a bode fine transcriptional reference to physiological aging and thus a valuable resource for the scientific community. It allowed us to demonstrate that the aged HSC pool is heterogeneous and that some aged cells are transcriptionally similar to young HSCs. We provide both (bulk and single-cell) analyses as a resource to the community at https://eriba.umcg.nl/agingsignature.

Machine learning algorithms independently identify the best age-associated gene predictors

Employing machine learning strategies in conjunction with scRNA-seq data have allowed for prediction and detection of different hematopoietic populations.35,36  We envisioned that by applying this approach to the age-related scRNA-seq data sets would allow for the identification of genes which best divide young and old HSCs. This approach would also serve as an independent and alternative approach to DEG analysis. We trained and tested 2 different machine learning algorithms, ADAboost and Random Forest (RNDforest), for identification of the best gene predictors. Interestingly, both algorithms determined that Nupr1, Selp, and Sult1a1 were the best predictors (supplemental Figure 5A). We tried to restrict the aging-related list to the best 20 age-associated genes. With those 20 genes, both algorithms performed well on identifying young and aged HSCs (Figure 5A-B; supplemental Figure 5B-C) when compared with the full list of AS genes. Attempts to further narrow down the age predictors from 20 to 5 genes, however, showed considerable deterioration in prediction capacity compared with larger lists of age predictors (supplemental Figure 5D). Out of 20 gene predictors, the majority are also membrane related (Figure 5C).

Figure 5.

Machine learning identifies best 20 gene predictors for separating young and aged HSCs in scRNA-seq. (A) Machine learning is able to predict with high accuracy young and aged HSCs. Top panel depicts a schematic representation of machine learning application into transcriptomic data from different scRNA-seq sets. Lower panel represents the output of one of the algorithm used (ADAboost) in individual cells from different sets. Young (circles) and aged (triangles) single cells were separated by PC1 (x-axis) and PC2 (y-axis) and color-coded according to the match between the measured aged of the cells and the predicted aged measured by the algorithm (gray for young cells, red for aged cells, and orange for mismatched cells). (B) Machine learning scores varies depending on which training set is used. Heatmap from ADAboost training depicting different overall scores for different scRNA-seq sets used. The overall score is color-coded from blue (lower scores) to yellow (higher scores). Sets on the x-axis (training sets) were used to train the algorithm, and the following sets on the y-axis (test sets) were scored according to training. (C) The best machine learning gene predictors also have high enrichment for membrane-associated proteins. Horizontal stacked plot comparing the list of the 20 best predictors extracted from machine learning algorithms (best predictors) and AS genes. The percentage of genes in each cellular location is represented as a percentage (x-axis) and divided by category (different colors). The 20 best predictor gene symbols are represented below. ER, endoplasmic reticulum.

Figure 5.

Machine learning identifies best 20 gene predictors for separating young and aged HSCs in scRNA-seq. (A) Machine learning is able to predict with high accuracy young and aged HSCs. Top panel depicts a schematic representation of machine learning application into transcriptomic data from different scRNA-seq sets. Lower panel represents the output of one of the algorithm used (ADAboost) in individual cells from different sets. Young (circles) and aged (triangles) single cells were separated by PC1 (x-axis) and PC2 (y-axis) and color-coded according to the match between the measured aged of the cells and the predicted aged measured by the algorithm (gray for young cells, red for aged cells, and orange for mismatched cells). (B) Machine learning scores varies depending on which training set is used. Heatmap from ADAboost training depicting different overall scores for different scRNA-seq sets used. The overall score is color-coded from blue (lower scores) to yellow (higher scores). Sets on the x-axis (training sets) were used to train the algorithm, and the following sets on the y-axis (test sets) were scored according to training. (C) The best machine learning gene predictors also have high enrichment for membrane-associated proteins. Horizontal stacked plot comparing the list of the 20 best predictors extracted from machine learning algorithms (best predictors) and AS genes. The percentage of genes in each cellular location is represented as a percentage (x-axis) and divided by category (different colors). The 20 best predictor gene symbols are represented below. ER, endoplasmic reticulum.

Close modal

The top aging gene, Selp, affects HSC functioning

The most consistent gene from the AS, Selp (supplemental Table 2), was also among the best machine learning predictors (Figure 5). Selp, which encodes for the membrane adhesion molecule P-selectin, was confirmed to be differently expressed also at the protein levels in young and aged hematopoietic progenitor subsets by flow cytometry (Figure 6A; supplemental Figure 6A). To investigate whether Selp affects HSC functioning, we competitively transplanted 1000 SELPlow or 1000 SELPhigh aged LT-HSCs (supplemental Figure 6B) into lethally irradiated young recipients. Whereas overall engraftment levels were similar for SELPlow and SELPhigh LT-HSCs, we found that SELPhigh LT-HSCs produced more myeloid cells, suggesting that SELPhigh HSCs are myeloid biased (Figure 6B; supplemental Figure 6C).

Figure 6.

Selp is the top AS candidate and contributes to HSC functional decline. (A) SELP protein expression level across different hematopoietic stem and progenitors (HSPCs) in young (gray) and aged (red) mice. Panel shows FCs of aged (n = 3 mice) samples normalized to their young (n = 3 mice) counterparts. Dots show biological replicates, while the bars represent the mean values for the given population. Flow cytometry plot for SELP expression in young and aged LT-HSCs population in the 3 biological replicates. (B) Selphigh LT-HSCs display a myeloid bias. The upper panel shows the experimental approach for competitive transplantation. SELPlow (blue) and SELPhigh (orange) LT-HSCs were purified from aged mice and competitively transplanted (n = 6-7 mice for each group). Overall engraftment (white blood cells [WBC]) and T cells levels were similar comparing both groups, whereas myeloid (My) and T-cell levels indicate that transplanted LT-HSC SELPhigh have myeloid bias. (C) SELP overexpression (OE) recipients show impairment of erythroid differentiation. The upper panel shows the experimental approach for competitive transplantation of SELP overexpression and empty vector (EV) HSCs. Lower panel displays peripheral blood (PB) output of the GFP+ (transduced cells) of both groups (n = 3 per group). Within the GFP+ population, the Ter-119+ population is significantly reduced. For all panels, mean ± standard deviation is shown. ns, nonsignificant; ****P < .0001; n indicates biological replicates.

Figure 6.

Selp is the top AS candidate and contributes to HSC functional decline. (A) SELP protein expression level across different hematopoietic stem and progenitors (HSPCs) in young (gray) and aged (red) mice. Panel shows FCs of aged (n = 3 mice) samples normalized to their young (n = 3 mice) counterparts. Dots show biological replicates, while the bars represent the mean values for the given population. Flow cytometry plot for SELP expression in young and aged LT-HSCs population in the 3 biological replicates. (B) Selphigh LT-HSCs display a myeloid bias. The upper panel shows the experimental approach for competitive transplantation. SELPlow (blue) and SELPhigh (orange) LT-HSCs were purified from aged mice and competitively transplanted (n = 6-7 mice for each group). Overall engraftment (white blood cells [WBC]) and T cells levels were similar comparing both groups, whereas myeloid (My) and T-cell levels indicate that transplanted LT-HSC SELPhigh have myeloid bias. (C) SELP overexpression (OE) recipients show impairment of erythroid differentiation. The upper panel shows the experimental approach for competitive transplantation of SELP overexpression and empty vector (EV) HSCs. Lower panel displays peripheral blood (PB) output of the GFP+ (transduced cells) of both groups (n = 3 per group). Within the GFP+ population, the Ter-119+ population is significantly reduced. For all panels, mean ± standard deviation is shown. ns, nonsignificant; ****P < .0001; n indicates biological replicates.

Close modal

To determine the role of Selp in the functioning of LT-HSCs, we overexpressed Selp in young LT-HSCs and competitively transplanted these in young recipients. Whereas Selp overexpression did not affect differentiation to leukocytes and platelets, erythroid production was almost completely abolished (Figure 6C). We also employed a Tet-On H2B-GFP reporter mouse to determine the dormancy or activation status of SELPlow and SELPhigh LT-HSCs. GFP fluorescence was significantly higher in SELPlow LT-HSCs, indicating that myeloid-biased SELPhigh LT-HSCs were more actively cycling (supplemental Figure 6D).

Thus, higher levels of SELP in LT-HSCs functionally contribute to myeloid bias, erythroid blockage, and HSC activation. Collectively, this demonstrates that SELP is not only a marker of aged HSCs but has the capacity to contribute to the functional impairment of HSCs.

In this study, we demonstrate that individual studies have limited resolution to identify genes altered in HSC aging, but collectively, they provide an unprecedented profile of transcriptomic alterations in HSC aging (Figure 1D). From this extensive repertoire of studies, we were able to identify a comprehensive, robust, and stable transcriptomic signature of HSC aging, which we validated by 3 independent data sets.

The HSC AS is highly enriched for membrane-associated transcripts (Figure 2A). Several of these have been previously associated with aspects of the HSC aging process and thus document the rigor of our approach (Table 1). However, the role of most other genes in HSC (aging) biology remains unknown. The strong involvement of cell-surface molecules does provide insight into the long-standing discussion whether HSC aging results from cell-intrinsic or cell-extrinsic molecular alterations. This enrichment of cell membrane–associated transcripts suggests that physiologically aged HSCs communicate differently with their immediate environment compared with their young counterparts.

Our study also provides novel insight into the transcriptional activation state of HSCs (Figure 3). In line with recent reports that loss of heterochromatin leads to age-associated phenotypes in HSCs,37  RNA-specific quantification, RNA polymerase II activity, and chromatin accessibility analyses all indicate that aged HSCs are, at a population level, more transcriptional active than young HSCs.

The AS can be employed as a bona fide transcriptional reference for physiological HSC aging for both bulk and single-cell transcriptomic studies (Figure 4). We provide protocols and a user-friendly and interactive website (https://agingsignature.webhosting.rug.nl/) to allow users to execute custom queries and compare new data sets against the AS that we present here.

Our approaches also allowed powerful analyses and uncovered novel perspectives into HSC aging from a single-cell perspective. By applying our AS score to single-cell transcriptomic studies, we determined that HSCs isolated from aged mice contain a fraction of cells which behave, at least transcriptionally, like young HSCs (Figure 4C). Further studies are necessary to understand if such cells indeed are functionally superior and also why such cells do not seem to “age.” We also successfully employed machine leaning to further identify best age-associated gene predictors (Figure 5).

Both computational approaches demonstrated that Selp expression is strongly associated with HSC aging. Previous data have shown that Selp-deficient HSCs display a competitive advantage over wild-type HSCs.38  Here, we document that high levels of SELP expression induced myeloid-biased, erythroid blockage and cell activation, all hallmarks of HSC aging8,9  (Figure 6). Together with Selp, which is a well-known inflammatory marker, Nupr1 is one of the top-listed predictors and belongs to a group of inflammatory and stress-responding transcription factors not earlier detected in relation to the aging phenotype. Although we were able to reduce the list of age-related genes, it is clear that HSC aging cannot be reduced to a single gene. Instead, our machine learning analysis suggests that multiple age-related genes are involved in the HSC aging process. The hypothesis of a single gene explaining most of the aging phenotype sounds unlikely, and it is more likely that multiple cell surface markers are required to prospectively separate young and aged HSCs.

Nevertheless, our data and those of others demonstrate that (at least some of) the aging genes that we identify not only are markers of aged HSCs but also are involved in age-dependent functional decline. The identification of the HSC AS does not resolve whether altered expression is the cause or consequence of loss of HSC functioning during aging. Future functional studies are required to assess whether perturbation of AS genes can improve functioning of aged stem cells. We believe that this resource is valuable in assisting the aging research community to verify involvement of genes in physiological HSC aging.

The authors thank T. Bijma, G. Mesander, and J. Teunis from University Medical Center Groningen (UMCG) Flowcytometry Unit facilities for their assistance on cell sorting; Klaas Sjollema from UMCG Microscopy and Imaging Center for assistance with confocal microscopy; and B. Bakker, M. Gerritsen, E. Verovskaya, and all members of the Ageing Biology and Stem Cells Laboratory for discussions.

This work was supported by the Netherlands Organization for Scientific Research/Mouse Clinic for Cancer and Aging, the Landsteiner Foundation for Blood Transfusion Research (LSBR1703), a China Student Council Fellowship (D.Y.), Marriage, a EU FP7 Marie Curie Initial Training Network (contract 316964), and ARCH (a European Union’s Horizon 2020 Research and Innovation Program) under Marie Skłodowska-Curie grant agreement 813091. B.v.E. was supported by grants from the Deutsche Forschungsgemeinschaft (EY 120/1-1), Else Kröner-Fresenius Foundation (2016_A58), and the German Cancer Aid (Deutsche Krebshilfe, 70113138). The Fritz Lipmann Institute is a member of the Leibniz Association and is financially supported by the Federal Government of Germany and the State of Thuringia.

Contribution: A.F.S., D.Y., S.L., G.d.H., and L.B. conceptualized the study; A.F.S., D.Y., and L.B. carried out the methodology; A.F.S., D.Y., S.S., A.A., K.K., A.M.-M., B.v.E., and L.B. performed the investigation; A.F.S., D.Y., S.S., K.K., A.M.-M., B.v.E., and L.B. performed the formal analysis; A.F.S., E.Z., and L.B. curated the data; E.Z. and A.A. collected the data; A.F.S. wrote the original manuscript draft; A.F.S., D.Y., G.d.H., and L.B. reviewed and edited the manuscript; and G.d.H. and L.B. supervised the study.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Correspondence: Gerald de Haan, Sanquin Research and Landsteiner Laboratory, Academic Medical Center, University of Amsterdam, Amsterdam, The Netherlands; e-mail: g.dehaan@sanquin.nl.

ATAC-seq data were deposited under the accession number GSE166674.

The online version of this article contains a data supplement.

There is a Blood Commentary on this article in this issue.

The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked "advertisement" in accordance with 18 USC section 1734.

1.
de Haan
G
,
Nijhof
W
,
Van Zant
G
.
Mouse strain-dependent changes in frequency and proliferation of hematopoietic stem cells during aging: correlation between lifespan and cycling activity
.
Blood.
1997
;
89
(
5
):
1543
-
1550
.
2.
Dykstra
B
,
Olthof
S
,
Schreuder
J
,
Ritsema
M
,
de Haan
G
.
Clonal analysis reveals multiple functional defects of aged murine hematopoietic stem cells
.
J Exp Med.
2011
;
208
(
13
):
2691
-
2703
.
3.
Morrison
SJ
,
Wandycz
AM
,
Akashi
K
,
Globerson
A
,
Weissman
IL
.
The aging of hematopoietic stem cells
.
Nat Med.
1996
;
2
(
9
):
1011
-
1016
.
4.
Rossi
DJ
,
Jamieson
CHM
,
Weissman
IL
.
Stems cells and the pathways to aging and cancer
.
Cell.
2008
;
132
(
4
):
681
-
696
.
5.
Sudo
K
,
Ema
H
,
Morita
Y
,
Nakauchi
H
.
Age-associated characteristics of murine hematopoietic stem cells
.
J Exp Med.
2000
;
192
(
9
):
1273
-
1280
.
6.
Verovskaya
E
,
Broekhuis
MJC
,
Zwart
E
, et al
.
Asymmetry in skeletal distribution of mouse hematopoietic stem cell clones and their equilibration by mobilizing cytokines
.
J Exp Med.
2014
;
211
(
3
):
487
-
497
.
7.
Verovskaya
E
,
Broekhuis
MJC
,
Zwart
E
, et al
.
Heterogeneity of young and aged murine hematopoietic stem cells revealed by quantitative clonal analysis using cellular barcoding
.
Blood.
2013
;
122
(
4
):
523
-
532
.
8.
López-Otín
C
,
Blasco
MA
,
Partridge
L
,
Serrano
M
,
Kroemer
G
.
The hallmarks of aging
.
Cell.
2013
;
153
(
6
):
1194
-
1217
.
9.
Singh
PP
,
Demmitt
BA
,
Nath
RD
,
Brunet
A
.
The genetics of aging: a vertebrate perspective
.
Cell.
2019
;
177
(
1
):
200
-
220
.
10.
Beerman
I
,
Rossi
DJ
.
Epigenetic control of stem cell potential during homeostasis, aging, and disease
.
Cell Stem Cell.
2015
;
16
(
6
):
613
-
625
.
11.
Bersenev
A
,
Rozenova
K
,
Balcerek
J
,
Jiang
J
,
Wu
C
,
Tong
W
.
Lnk deficiency partially mitigates hematopoietic stem cell aging
.
Aging Cell.
2012
;
11
(
6
):
949
-
959
.
12.
Chambers
SM
,
Shaw
CA
,
Gatza
C
,
Fisk
CJ
,
Donehower
LA
,
Goodell
MA
.
Aging hematopoietic stem cells decline in function and exhibit epigenetic dysregulation
.
PLoS Biol.
2007
;
5
(
8
):
e201
.
13.
Flach
J
,
Bakker
ST
,
Mohrin
M
, et al
.
Replication stress is a potent driver of functional decline in ageing haematopoietic stem cells
.
Nature.
2014
;
512
(
7513
):
198
-
202
.
14.
Grover
A
,
Sanjuan-Pla
A
,
Thongjuea
S
, et al
.
Single-cell RNA sequencing reveals molecular and functional platelet bias of aged haematopoietic stem cells
.
Nat Commun.
2016
;
7
(
1
):
11075
.
15.
Kirschner
K
,
Chandra
T
,
Kiselev
V
, et al
.
Proliferation drives aging-related functional decline in a subpopulation of the hematopoietic stem cell compartment
.
Cell Rep.
2017
;
19
(
8
):
1503
-
1511
.
16.
Kowalczyk
MS
,
Tirosh
I
,
Heckl
D
, et al
.
Single-cell RNA-seq reveals changes in cell cycle and differentiation programs upon aging of hematopoietic stem cells
.
Genome Res.
2015
;
25
(
12
):
1860
-
1872
.
17.
Mann
M
,
Mehta
A
,
de Boer
CG
, et al
.
Heterogeneous responses of hematopoietic stem cells to inflammatory stimuli are altered with age
.
Cell Rep.
2018
;
25
(
11
):
2992
-
3005
.
18.
Maryanovich
M
,
Zahalka
AH
,
Pierce
H
, et al
.
Adrenergic nerve degeneration in bone marrow drives aging of the hematopoietic stem cell niche [published correction appears in Nat Med. 2019;25:701]
.
Nat Med.
2018
;
24
(
6
):
782
-
791
.
19.
Noda
S
,
Ichikawa
H
,
Miyoshi
H
.
Hematopoietic stem cell aging is associated with functional decline and delayed cell cycle progression
.
Biochem Biophys Res Commun.
2009
;
383
(
2
):
210
-
215
.
20.
Norddahl
GL
,
Pronk
CJ
,
Wahlestedt
M
, et al
.
Accumulating mitochondrial DNA mutations drive premature hematopoietic aging phenotypes distinct from physiological stem cell aging
.
Cell Stem Cell.
2011
;
8
(
5
):
499
-
510
.
21.
Quéré
R
,
Saint-Paul
L
,
Carmignac
V
, et al
.
Tif1γ regulates the TGF-β1 receptor and promotes physiological aging of hematopoietic stem cells
.
Proc Natl Acad Sci USA.
2014
;
111
(
29
):
10592
-
10597
.
22.
Rossi
DJ
,
Bryder
D
,
Zahn
JM
, et al
.
Cell intrinsic alterations underlie hematopoietic stem cell aging
.
Proc Natl Acad Sci USA.
2005
;
102
(
26
):
9194
-
9199
.
23.
Sun
D
,
Luo
M
,
Jeong
M
, et al
.
Epigenomic profiling of young and aged HSCs reveals concerted changes during aging that reinforce self-renewal
.
Cell Stem Cell.
2014
;
14
(
5
):
673
-
688
.
24.
Wahlestedt
M
,
Norddahl
GL
,
Sten
G
, et al
.
An epigenetic component of hematopoietic stem cell aging amenable to reprogramming into a young state
.
Blood.
2013
;
121
(
21
):
4257
-
4264
.
25.
Beerman
I
,
Bock
C
,
Garrison
BS
, et al
.
Proliferation-dependent alterations of the DNA methylation landscape underlie hematopoietic stem cell aging
.
Cell Stem Cell.
2013
;
12
(
4
):
413
-
425
.
26.
Renders
S
,
Svendsen
AF
,
Panten
J
, et al
.
Niche derived netrin-1 regulates hematopoietic stem cell dormancy via its receptor neogenin-1
.
Nat Commun.
2021
;
12
(
1
):
608
.
27.
Buenrostro
JD
,
Giresi
PG
,
Zaba
LC
,
Chang
HY
,
Greenleaf
WJ
.
Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position
.
Nat Methods.
2013
;
10
(
12
):
1213
-
1218
.
28.
Hidalgo
I
,
Herrera-Merchan
A
,
Ligos
JM
, et al
.
Ezh1 is required for hematopoietic stem cell maintenance and prevents senescence-like cell cycle arrest
.
Cell Stem Cell.
2012
;
11
(
5
):
649
-
662
.
29.
Challen
GA
,
Sun
D
,
Mayle
A
, et al
.
Dnmt3a and Dnmt3b have overlapping and distinct functions in hematopoietic stem cells
.
Cell Stem Cell.
2014
;
15
(
3
):
350
-
364
.
30.
Wong
DJ
,
Liu
H
,
Ridky
TW
,
Cassarino
D
,
Segal
E
,
Chang
HY
.
Module map of stem cell genes guides creation of epithelial cancer stem cells
.
Cell Stem Cell.
2008
;
2
(
4
):
333
-
344
.
31.
de Magalhães
JP
,
Curado
J
,
Church
GM
.
Meta-analysis of age-related gene expression profiles identifies common signatures of aging
.
Bioinformatics.
2009
;
25
(
7
):
875
-
881
.
32.
Blank
U
,
Karlsson
S
.
TGF-β signaling in the control of hematopoietic stem cells
.
Blood.
2015
;
125
(
23
):
3542
-
3550
.
33.
Subramanian
A
,
Tamayo
P
,
Mootha
VK
, et al
.
Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles
.
Proc Natl Acad Sci USA.
2005
;
102
(
43
):
15545
-
15550
.
34.
Liberzon
A
,
Birger
C
,
Thorvaldsdóttir
H
,
Ghandi
M
,
Mesirov
JP
,
Tamayo
P
.
The Molecular Signatures Database (MSigDB) hallmark gene set collection
.
Cell Syst.
2015
;
1
(
6
):
417
-
425
.
35.
Hamey
FK
,
Göttgens
B
.
Machine learning predicts putative hematopoietic stem cells within large single-cell transcriptomics data sets
.
Exp Hematol.
2019
;
78
:
11
-
20
.
36.
Fidanza
A
,
Kingdom
U
,
Stumpf
P
, et al
.
Single-cell analyses and machine learning define hematopoietic progenitor and HSC-like cells derived from human PSCs
.
Blood.
2020
;
136
(
25
):
2893
-
2904
.
37.
Keenan
CR
,
Iannarella
N
,
Naselli
G
, et al
.
Extreme disruption of heterochromatin is required for accelerated hematopoietic aging
.
Blood.
2020
;
135
(
23
):
2049
-
2058
.
38.
Sullivan
C
,
Chen
Y
,
Shan
Y
, et al
.
Functional ramifications for the loss of P-selectin expression on hematopoietic and leukemic stem cells
.
PLoS One.
2011
;
6
(
10
):
e26246
.
39.
Rapisarda
V
,
Borghesan
M
,
Miguela
V
, et al
.
Integrin beta 3 regulates cellular senescence by activating the TGF-β pathway
.
Cell Rep.
2017
;
18
(
10
):
2480
-
2493
.
40.
Sanjuan-Pla
A
,
Macaulay
IC
,
Jensen
CT
, et al
.
Platelet-biased stem cells reside at the apex of the haematopoietic stem-cell hierarchy
.
Nature.
2013
;
502
(
7470
):
232
-
236
.
41.
Pinho
S
,
Marchand
T
,
Yang
E
,
Wei
Q
,
Nerlov
C
,
Frenette
PS
.
Lineage-biased hematopoietic stem cells are regulated by distinct niches
.
Dev Cell.
2018
;
44
(
5
):
634
-
641.e4
.
42.
Gulati
GS
,
Zukowska
M
,
Noh
JJ
, et al
.
Neogenin-1 distinguishes between myeloid-biased and balanced Hoxb5+ mouse long-term hematopoietic stem cells
.
Proc Natl Acad Sci USA.
2019
;
116
(
50
):
25115
-
25125
.
43.
Jeannet
R
,
Cai
Q
,
Liu
H
,
Vu
H
,
Kuo
Y-H
.
Alcam regulates long-term hematopoietic stem cell engraftment and self-renewal
.
Stem Cells.
2013
;
31
(
3
):
560
-
571
.
44.
Shimazu
T
,
Iida
R
,
Zhang
Q
, et al
.
CD86 is expressed on murine hematopoietic stem cells and denotes lymphopoietic potential
.
Blood.
2012
;
119
(
21
):
4889
-
4897
.
45.
Thalheimer
FB
,
Wingert
S
,
De Giacomo
P
, et al
.
Cytokine-regulated GADD45G induces differentiation and lineage selection in hematopoietic stem cells
.
Stem Cell Reports.
2014
;
3
(
1
):
34
-
43
.
46.
Satoh
Y
,
Yokota
T
,
Sudo
T
, et al
.
The Satb1 protein directs hematopoietic stem cell differentiation toward lymphoid lineages
.
Immunity.
2013
;
38
(
6
):
1105
-
1115
.
47.
Taiwo
O
,
Wilson
GA
,
Emmett
W
, et al
.
DNA methylation analysis of murine hematopoietic side population cells during aging
.
Epigenetics.
2013
;
8
(
10
):
1114
-
1122
.
48.
Oguro
H
,
Ding
L
,
Morrison
SJ
.
SLAM family markers resolve functionally distinct subpopulations of hematopoietic stem cells and multipotent progenitors
.
Cell Stem Cell.
2013
;
13
(
1
):
102
-
116
.
49.
Beerman
I
,
Bhattacharya
D
,
Zandi
S
, et al
.
Functionally distinct hematopoietic stem cells modulate hematopoietic lineage potential during aging by a mechanism of clonal expansion
.
Proc Natl Acad Sci USA.
2010
;
107
(
12
):
5465
-
5470
.
50.
Tadokoro
Y
,
Ema
H
,
Okano
M
,
Li
E
,
Nakauchi
H
.
De novo DNA methyltransferase is essential for self-renewal, but not for differentiation, in hematopoietic stem cells
.
J Exp Med.
2007
;
204
(
4
):
715
-
722
.
51.
Challen
GA
,
Mayle
A
,
Sun
D
, et al
.
Dnmt3b has few specific functions in adult hematopoietic stem cells but shows abnormal activity in the absence of Dnmt3a
.
Blood.
2013
;
122
(
21
):
734
.

Author notes

*

G.d.H. and L.B. contributed equally to this study.

Supplemental data

Sign in via your Institution