Why an integrated view of gene expression studies on hematopoiesis in mouse aging is better than the sum of their parts.

Globally, the human population is aging, with an increased proportion of people in "old age" (over 60 years). This trend leads to a growing demand in aging research, stimulating studies in animal models such as mice, fish, and invertebrates. Recently, we published a research summary on the aging of hematopoietic stem cells (HSCs) in C57BL/6 mice based on 12 gene expression datasets. Here, I discuss in greater detail the added value of taking an integrated view, rather than considering each publication separately, to determine genes involved in aging. Considerable variation exists between lists of differentially expressed (DE) genes in HSCs, comparing young and old mice. This variation can result from factors such as inconsistent definitions of "young" and "old", technical variations and variations between laboratory mouse strains. We previously demonstrated that the variation between gene lists could be circumvented by forming a unified list of DE genes-the "aging list"-with citation indexes attached. The most frequently detected DE genes [approximately 200 most cited, which we named the "aging signature" (AS)] were highly consistent across publications. Gene Ontology classification of the AS list identified additional sources of variation between studies: one comes from the specifics of how the data are collected and analyzed; another comes from inconsistencies between how we define the gene categories. As discussed, overcoming these variations is the next challenge toward an integral approach to our systematic knowledge of the aging process.

Panhematopoietic RNA barcoding enables kinetic measurements of nucleate and anucleate lineages and the activation of myeloid clones following acute platelet depletion.

BACKGROUND: Platelets and erythrocytes constitute over 95% of all hematopoietic stem cell output. However, the clonal dynamics of HSC contribution to these lineages remains largely unexplored. RESULTS: We use lentiviral genetic labeling of mouse hematopoietic stem cells to quantify output from all lineages, nucleate, and anucleate, simultaneously linking these with stem and progenitor cell transcriptomic phenotypes using single-cell RNA-sequencing. We observe dynamic shifts of clonal behaviors through time in same-animal peripheral blood and demonstrate that acute platelet depletion shifts the output of multipotent hematopoietic stem cells to the exclusive production of platelets. Additionally, we observe the emergence of new myeloid-biased clones, which support short- and long-term production of blood cells. CONCLUSIONS: Our approach enables kinetic studies of multi-lineage output in the peripheral blood and transcriptional heterogeneity of individual hematopoietic stem cells. Our results give a unique insight into hematopoietic stem cell reactivation upon platelet depletion and of clonal dynamics in both steady state and under stress.

Taz protects hematopoietic stem cells from an aging-dependent decrease in PU.1 activity.

Specific functions of the immune system are essential to protect us from infections caused by pathogens such as viruses and bacteria. However, as we age, the immune system shows a functional decline that can be attributed in large part to age-associated defects in hematopoietic stem cells (HSCs)-the cells at the apex of the immune cell hierarchy. Here, we find that the Hippo pathway coactivator TAZ is potently induced in old HSCs and protects these cells from functional decline. We identify Clca3a1 as a TAZ-induced gene that allows us to trace TAZ activity in vivo. Using CLCA3A1 as a marker, we can isolate "young-like" HSCs from old mice. Mechanistically, Taz acts as coactivator of PU.1 and to some extent counteracts the gradual loss of PU.1 expression during HSC aging. Our work thus uncovers an essential role for Taz in a previously undescribed fail-safe mechanism in aging HSCs.

Measures of Clonal Hematopoiesis: Are We Missing Something?

Clonal Hematopoiesis (CH) is a common, age-related phenomenon of growing scientific interest, due to its association with hematologic malignancy, cardiovascular disease and decreased overall survival. CH is commonly attributed to the preferential outgrowth of a mutant hematopoietic stem cell (HSC) with enhanced fitness, resulting in clonal imbalance. In-depth understanding of the relation between HSC clonal dynamics, CH and hematologic malignancy requires integration of fundamental lineage tracing studies with clinical data. However, this is hampered by lack of a uniform definition of CH and by inconsistency in the analytical methods used for its quantification. Here, we propose a conceptual and analytical framework for the definition and measurement of CH. First, we transformed the conceptual definition of CH into the CH index, which provides a quantitative measure of clone numbers and sizes. Next, we generated a set of synthetic data, based on the beta-distribution, to simulate clonal populations with different degrees of imbalance. Using these clonal distributions and the CH index as a reference, we tested several established indices of clonal diversity and (in-)equality for their ability to detect and quantify CH. We found that the CH index was distinct from any of the other tested indices. Nonetheless, the diversity indices (Shannon, Simpson) more closely resembled the CH index than the inequality indices (Gini, Pielou). Notably, whereas the inequality indices mainly responded to changes in clone sizes, the CH index and the tested diversity indices also responded to changes in the number of clones in a sample. Accordingly, these simulations indicate that CH can result not only by skewing clonal abundancies, but also by variation in their overall numbers. Altogether, our model-based approach illustrates how a formalized definition and quantification of CH can provide insights into its pathogenesis. In the future, use of the CH index or Shannon index to quantify clonal diversity in fundamental as well as clinical clone-tracing studies will promote cross-disciplinary discussion and progress in the field.

A comprehensive transcriptome signature of murine hematopoietic stem cell aging.

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.

Niche derived netrin-1 regulates hematopoietic stem cell dormancy via its receptor neogenin-1.

Haematopoietic stem cells (HSCs) are characterized by their self-renewal potential associated to dormancy. Here we identify the cell surface receptor neogenin-1 as specifically expressed in dormant HSCs. Loss of neogenin-1 initially leads to increased HSC expansion but subsequently to loss of self-renewal and premature exhaustion in vivo. Its ligand netrin-1 induces Egr1 expression and maintains quiescence and function of cultured HSCs in a Neo1 dependent manner. Produced by arteriolar endothelial and periarteriolar stromal cells, conditional netrin-1 deletion in the bone marrow niche reduces HSC numbers, quiescence and self-renewal, while overexpression increases quiescence in vivo. Ageing associated bone marrow remodelling leads to the decline of netrin-1 expression in niches and a compensatory but reversible upregulation of neogenin-1 on HSCs. Our study suggests that niche produced netrin-1 preserves HSC quiescence and self-renewal via neogenin-1 function. Decline of netrin-1 production during ageing leads to the gradual decrease of Neo1 mediated HSC self-renewal.

Python for gene expression.

Genome biology shows substantial progress in its analytical and computational part in the last decades. Differential gene expression is one of many computationally intense areas; it is largely developed under R programming language. Here we explain possible reasons for such dominance of R in gene expression data. Next, we discuss the prospects for Python to become competitive in this area of research in coming years. We indicate that Python can be used already in a field of a single cell differential gene expression. We pinpoint still missing parts in Python and possibilities for improvement.

Persistent expression of microRNA-125a targets is required to induce murine hematopoietic stem cell repopulating activity.

MicroRNAs (miRs) are small noncoding RNAs that regulate gene expression posttranscriptionally by binding to the 3' untranslated regions of their target mRNAs. The evolutionarily conserved microRNA-125a (miR-125a) is highly expressed in both murine and human hematopoietic stem cells (HSCs), and previous studies have found that miR-125 strongly enhances self-renewal of HSCs and progenitors. In this study we explored whether temporary overexpression of miR-125a would be sufficient to permanently increase HSC self-renewal or, rather, whether persistent overexpression of miR-125a is required. We used three complementary in vivo approaches to reversibly enforce expression of miR-125a in murine HSCs. Additionally, we interrogated the underlying molecular mechanisms responsible for the functional changes that occur in HSCs on overexpression of miR-125a. Our data indicate that continuous expression of miR-125a is required to enhance HSC activity. Our molecular analysis confirms changes in pathways that explain the characteristics of miR-125a overexpressing HSCs. Moreover, it provides several novel putative miR-125a targets, but also highlights the complex molecular changes that collectively lead to enhanced HSC function.

Clonal Analysis of Patient-Derived Samples Using Cellular Barcodes.

Cellular barcoding is a relatively simple method that allows quantitative assessment of the clonal dynamics of normal, nonmalignant hematopoietic stem cells and of leukemia. Cellular barcodes are (semi-)random synthetic DNA sequences of a fixed length, which are used to uniquely mark and track cells over time. A successful barcoding experiment consists of several essential steps, including library production, transfection, transduction, barcode retrieval, and barcode data analysis. Key challenges are to obtain sufficient number of barcoded cells to conduct experiments and reliable barcode data analysis. This is especially relevant for experiments using primary leukemia cells (which are of limited availability and difficult to transduce), when studying low levels of chimerism, or when the barcoded cell population is sorted in different smaller subpopulations (e.g., lineage contribution of normal hematopoietic stem cells in murine xenografts). In these settings, retrieving accurate barcode data from low input material using standard PCR amplification techniques might be challenging and more sophisticated approaches are required. In this chapter we describe the procedures to transfect and transduce patient-derived leukemia cells, to retrieve barcoded data from both high and low input material, and to filter barcode data from sequencing noise prior to quantitative clonal analysis.

Detection of chemotherapy-resistant patient-derived acute lymphoblastic leukemia clones in murine xenografts using cellular barcodes.

Clonal heterogeneity fuels leukemia evolution, therapeutic resistance, and relapse. Upfront detection of therapy-resistant leukemia clones at diagnosis may allow adaptation of treatment and prevention of relapse, but this is hampered by a paucity of methods to identify and trace single leukemia-propagating cells and their clonal offspring. Here, we tested methods of cellular barcoding analysis, to trace the in vivo competitive dynamics of hundreds of patient-derived leukemia clones upon chemotherapy-mediated selective pressure. We transplanted Nod/Scid/Il2Rγ-/- (NSG) mice with barcoded patient-derived or SupB15 acute lymphoblastic leukemia (ALL) cells and assessed clonal responses to dexamethasone, methotrexate, and vincristine, longitudinally and across nine anatomic locations. We illustrate that chemotherapy reduces clonal diversity in a drug-dependent manner. At end-stage disease, methotrexate-treated patient-derived xenografts had significantly fewer clones compared with placebo-treated mice (100 ± 10 vs. 160 ± 15 clones, p = 0.0005), while clonal complexity in vincristine- and dexamethasone-treated xenografts was unaffected (115 ± 33 and 150 ± 7 clones, p = NS). Using tools developed to assess differential gene expression, we determined whether these clonal patterns resulted from random clonal drift or selection. We identified 5 clones that were reproducibly enriched in methotrexate-treated patient-derived xenografts, suggestive of pre-existent resistance. Finally, we found that chemotherapy-mediated selection resulted in a more asymmetric distribution of leukemia clones across anatomic sites. We found that cellular barcoding is a powerful method to trace the clonal dynamics of human patient-derived leukemia cells in response to chemotherapy. In the future, integration of cellular barcoding with single-cell sequencing technology may allow in-depth characterization of therapy-resistant leukemia clones and identify novel targets to prevent relapse.

C/EBPα and GATA-2 Mutations Induce Bilineage Acute Erythroid Leukemia through Transformation of a Neomorphic Neutrophil-Erythroid Progenitor.

Acute erythroid leukemia (AEL) commonly involves both myeloid and erythroid lineage transformation. However, the mutations that cause AEL and the cell(s) that sustain the bilineage leukemia phenotype remain unknown. We here show that combined biallelic Cebpa and Gata2 zinc finger-1 (ZnF1) mutations cooperatively induce bilineage AEL, and that the major leukemia-initiating cell (LIC) population has a neutrophil-monocyte progenitor (NMP) phenotype. In pre-leukemic NMPs Cebpa and Gata2 mutations synergize by increasing erythroid transcription factor (TF) expression and erythroid TF chromatin access, respectively, thereby installing ectopic erythroid potential. This erythroid-permissive chromatin conformation is retained in bilineage LICs. These results demonstrate that synergistic transcriptional and epigenetic reprogramming by leukemia-initiating mutations can generate neomorphic pre-leukemic progenitors, defining the lineage identity of the resulting leukemia.

Correction: Quantitative distribution of patient-derived leukemia clones in murine xenografts revealed by cellular barcodes.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Donor-to-Donor Heterogeneity in the Clonal Dynamics of Transplanted Human Cord Blood Stem Cells in Murine Xenografts.

Umbilical cord blood (UCB) provides an alternative source of hematopoietic stem cells (HSCs) for allogeneic transplantation. Administration of sufficient donor HSCs is critical to restore recipient hematopoiesis and to maintain long-term polyclonal blood formation. However, due to lack of unique markers, the frequency of HSCs among UCB CD34+ cells is the subject of ongoing debate, urging for reproducible strategies for their counting. Here, we used cellular barcoding to determine the frequency and clonal dynamics of human UCB HSCs and to determine how data analysis methods affect these parameters. We transplanted lentivirally barcoded CD34+ cells from 20 UCB donors into Nod/Scid/IL2Ry-/- (NSG) mice (n = 30). Twelve recipients (of 8 UCB donors) engrafted with >1% GFP+ cells, allowing for clonal analysis by multiplexed barcode deep sequencing. Using multiple definitions of clonal diversity and strategies for data filtering, we demonstrate that differences in data analysis can change clonal counts by several orders of magnitude and propose methods to improve their consistency. Using these methods, we show that the frequency of NSG-repopulating cells was low (median ∼1 HSC/104 CD34+ UCB cells) and could vary up to 10-fold between donors. Clonal patterns in blood became increasingly consistent over time, likely reflecting initial output of transient progenitors, followed by long-term HSCs with stable hierarchies. The majority of long-term clones displayed multilineage output, yet clones with lymphoid- or myeloid-biased output were also observed. Altogether, this study uncovers substantial interdonor and analysis-induced variability in the frequency of UCB CD34+ clones that contribute to post-transplant hematopoiesis. As clone tracing is increasingly relevant, we urge for universal and transparent methods to count HSC clones during normal aging and upon transplantation.

MiR-125a enhances self-renewal, lifespan, and migration of murine hematopoietic stem and progenitor cell clones.

Expansion of hematopoietic stem cells (HSCs) is a 'holy grail' of regenerative medicine, as successful stem cell transplantations depend on the number and quality of infused HSCs. Although many attempts have been pursued to either chemically or genetically increase HSC numbers, neither clonal analysis of these expanded cells nor their ability to support mature blood lineages has been demonstrated. Here we show that miR-125a, at the single cell level, can expand murine long-term repopulating HSCs. In addition, miR-125a increases clone longevity, clone size and clonal contribution to hematopoiesis. Unexpectedly, we found that miR-125a expanded HSCs clones were highly homogenously distributed across multiple anatomical sites. Interestingly, these miR-125a overexpressing cells had enhanced mobility and were more frequently detected in the spleen. Our study reveals a novel, cell-intrinsically controlled mechanism by which HSC migration is regulated.

CBX7 Induces Self-Renewal of Human Normal and Malignant Hematopoietic Stem and Progenitor Cells by Canonical and Non-canonical Interactions.

In this study, we demonstrate that, among all five CBX Polycomb proteins, only CBX7 possesses the ability to control self-renewal of human hematopoietic stem and progenitor cells (HSPCs). Xenotransplantation of CBX7-overexpressing HSPCs resulted in increased multi-lineage long-term engraftment and myelopoiesis. Gene expression and chromatin analyses revealed perturbations in genes involved in differentiation, DNA and chromatin maintenance, and cell cycle control. CBX7 is upregulated in acute myeloid leukemia (AML), and its genetic or pharmacological repression in AML cells inhibited proliferation and induced differentiation. Mass spectrometry analysis revealed several non-histone protein interactions between CBX7 and the H3K9 methyltransferases SETDB1, EHMT1, and EHMT2. These CBX7-binding proteins possess a trimethylated lysine peptide motif highly similar to the canonical CBX7 target H3K27me3. Depletion of SETDB1 in AML cells phenocopied repression of CBX7. We identify CBX7 as an important regulator of self-renewal and uncover non-canonical crosstalk between distinct pathways, revealing therapeutic opportunities for leukemia.

MicroRNA High Throughput Loss-of-Function Screening Reveals an Oncogenic Role for miR-21-5p in Hodgkin Lymphoma.

BACKGROUND/AIMS: Classical Hodgkin lymphoma (cHL) is among the most frequent lymphoma subtypes. The tumor cells originate from crippled germinal center (GC)-B cells that escaped from apoptosis. MicroRNAs (miRNAs) play important roles in B-cell maturation and aberrant expression of miRNAs contributes to the pathogenesis of cHL. Our aim was to identify oncogenic miRNAs relevant for growth of cHL using a high-throughput screening approach. METHODS: A lentiviral pool of 63 miRNA inhibition constructs was used to identify miRNAs essential to cell growth in three cHL cell lines in duplicate. As a negative control we also infected cHL cell lines with a lentiviral barcoded empty vector pool consisting of 222 constructs. The abundance of individual constructs was followed over time by a next generation sequencing approach. The effect on growth was confirmed using individual GFP competition assays and on apoptosis using Annexin-V staining. Our previously published Argonaute 2 (Ago2) immunoprecipitation (IP) data were used to identify target genes relevant for cell growth / apoptosis. Luciferase assays and western blotting were performed to confirm targeting by miRNAs. RESULTS: Four miRNA inhibition constructs, i.e. miR-449a-5p, miR-625-5p, let-7f-2-3p and miR-21-5p, showed a significant decrease in abundance in at least 4 of 6 infections. In contrast, none of the empty vector constructs showed a significant decrease in abundance in 3 or more of the 6 infections. The most abundantly expressed miRNA, i.e. miR-21-5p, showed significantly higher expression levels in cHL compared to GC-B cells. GFP competition assays confirmed the negative effect of miR-21-5p inhibition on HL cell growth. Annexin-V staining of cells infected with miR-21-5p inhibitor indicated a significant increase in apoptosis at day 7 and 9 after viral infection, consistent with the decrease in growth. Four miR-21-5p cell growth- and apoptosis-associated targets were AGO2-IP enriched in cHL cell lines and showed a significant decrease in expression in cHL cell lines in comparison to normal GC-B cells. For the two most abundantly expressed, i.e. BTG2 and PELI1, we confirmed targeting by miR-21-5p using luciferase assays and for PELI1 we also confirmed this at the protein level by western blotting. CONCLUSION: Using a miRNA loss-of-function high-throughput screen we identified four miRNAs with oncogenic effects in cHL and validated the results for the in cHL abundantly expressed miR-21-5p. MiR-21-5p is upregulated in cHL compared to GC-B cells and protects cHL cells from apoptosis possibly via targeting BTG2 and PELI1.

Left or right? Directions to stem cell engraftment.

In this issue of JEM, Wu et al. (https://doi.org/10.1084/jem.20171341) use genetic barcoding of macaque hematopoietic stem cells to demonstrate that, after transplantation, HSCs are very asymmetrically distributed and uncover a thymus-independent pathway for mature T cell production in the bone marrow.

Clonal selection and asymmetric distribution of human leukemia in murine xenografts revealed by cellular barcoding.

Genetic and phenotypic heterogeneity of human leukemia is thought to drive leukemia progression through a Darwinian process of selection and evolution of increasingly malignant clones. However, the lack of markers that uniquely identify individual leukemia clones precludes high-resolution tracing of their clonal dynamics. Here, we use cellular barcoding to analyze the clonal behavior of patient-derived leukemia-propagating cells (LPCs) in murine xenografts. Using a leukemic cell line and diagnostic bone marrow cells from 6 patients with B-progenitor cell acute lymphoblastic leukemia, we demonstrate that patient-derived xenografts were highly polyclonal, consisting of tens to hundreds of LPC clones. The number of clones was stable within xenografts but strongly reduced upon serial transplantation. In contrast to primary recipients, in which clonal composition was highly diverse, clonal composition in serial xenografts was highly similar between recipients of the same donor and reflected donor clonality, supporting a deterministic, clone-size-based model for clonal selection. Quantitative analysis of clonal abundance in several anatomic sites identified 2 types of anatomic asymmetry. First, clones were asymmetrically distributed between different bones. Second, clonal composition in the skeleton significantly differed from extramedullary sites, showing similar numbers but different clone sizes. Altogether, this study shows that cellular barcoding and xenotransplantation providea useful model to study the behavior of patient-derived LPC clones, which provides insights relevant for experimental studies on cancer stem cells and for clinical protocols for the diagnosis and treatment of leukemia.

Ectopic miR-125a Expression Induces Long-Term Repopulating Stem Cell Capacity in Mouse and Human Hematopoietic Progenitors.

Umbilical cord blood (CB) is a convenient and broadly used source of hematopoietic stem cells (HSCs) for allogeneic stem cell transplantation. However, limiting numbers of HSCs remain a major constraint for its clinical application. Although one feasible option would be to expand HSCs to improve therapeutic outcome, available protocols and the molecular mechanisms governing the self-renewal of HSCs are unclear. Here, we show that ectopic expression of a single microRNA (miRNA), miR-125a, in purified murine and human multipotent progenitors (MPPs) resulted in increased self-renewal and robust long-term multi-lineage repopulation in transplanted recipient mice. Using quantitative proteomics and western blot analysis, we identified a restricted set of miR-125a targets involved in conferring long-term repopulating capacity to MPPs in humans and mice. Our findings offer the innovative potential to use MPPs with enhanced self-renewal activity to augment limited sources of HSCs to improve clinical protocols.