Somatic epimutations enable single-cell lineage tracing in native hematopoiesis across the murine and human lifespan.

Current approaches to lineage tracing of stem cell clones require genetic engineering or rely on sparse somatic DNA variants, which are difficult to capture at single-cell resolution. Here, we show that targeted single-cell measurements of DNA methylation at single-CpG resolution deliver joint information about cellular differentiation state and clonal identities. We develop EPI-clone, a droplet-based method for transgene-free lineage tracing, and apply it to study hematopoiesis, capturing hundreds of clonal trajectories across almost 100,000 single-cells. Using ground-truth genetic barcodes, we demonstrate that EPI-clone accurately identifies clonal lineages throughout hematopoietic differentiation. Applied to unperturbed hematopoiesis, we describe an overall decline of clonal complexity during murine ageing and the expansion of rare low-output stem cell clones. In aged human donors, we identified expanded hematopoietic clones with and without genetic lesions, and various degrees of clonal complexity. Taken together, EPI-clone enables accurate and transgene-free single-cell lineage tracing at scale.

Integrative single-cell expression and functional studies unravels a sensitization to cytarabine-based chemotherapy through HIF pathway inhibition in AML leukemia stem cells.

Relapse remains a major challenge in the clinical management of acute myeloid leukemia (AML) and is driven by rare therapy-resistant leukemia stem cells (LSCs) that reside in specific bone marrow niches. Hypoxia signaling maintains cells in a quiescent and metabolically relaxed state, desensitizing them to chemotherapy. This suggests the hypothesis that hypoxia contributes to the chemoresistance of AML-LSCs and may represent a therapeutic target to sensitize AML-LSCs to chemotherapy. Here, we identify HIFhigh and HIFlow specific AML subgroups (inv(16)/t(8;21) and MLLr, respectively) and provide a comprehensive single-cell expression atlas of 119,000 AML cells and AML-LSCs in paired diagnostic-relapse samples from these molecular subgroups. The HIF/hypoxia pathway signature is attenuated in AML-LSCs compared with more differentiated AML cells but is more expressed than in healthy hematopoietic cells. Importantly, chemical inhibition of HIF cooperates with standard-of-care chemotherapy to impair AML growth and to substantially eliminate AML-LSCs in vitro and in vivo. These findings support the HIF pathway in the stem cell-driven drug resistance of AML and unravel avenues for combinatorial targeted and chemotherapy-based approaches to specifically eliminate AML-LSCs.

Biological basis of extensive pleiotropy between blood traits and cancer risk.

BACKGROUND: The immune system has a central role in preventing carcinogenesis. Alteration of systemic immune cell levels may increase cancer risk. However, the extent to which common genetic variation influences blood traits and cancer risk remains largely undetermined. Here, we identify pleiotropic variants and predict their underlying molecular and cellular alterations. METHODS: Multivariate Cox regression was used to evaluate associations between blood traits and cancer diagnosis in cases in the UK Biobank. Shared genetic variants were identified from the summary statistics of the genome-wide association studies of 27 blood traits and 27 cancer types and subtypes, applying the conditional/conjunctional false-discovery rate approach. Analysis of genomic positions, expression quantitative trait loci, enhancers, regulatory marks, functionally defined gene sets, and bulk- and single-cell expression profiles predicted the biological impact of pleiotropic variants. Plasma small RNAs were sequenced to assess association with cancer diagnosis. RESULTS: The study identified 4093 common genetic variants, involving 1248 gene loci, that contributed to blood-cancer pleiotropism. Genomic hotspots of pleiotropism include chromosomal regions 5p15-TERT and 6p21-HLA. Genes whose products are involved in regulating telomere length are found to be enriched in pleiotropic variants. Pleiotropic gene candidates are frequently linked to transcriptional programs that regulate hematopoiesis and define progenitor cell states of immune system development. Perturbation of the myeloid lineage is indicated by pleiotropic associations with defined master regulators and cell alterations. Eosinophil count is inversely associated with cancer risk. A high frequency of pleiotropic associations is also centered on the regulation of small noncoding Y-RNAs. Predicted pleiotropic Y-RNAs show specific regulatory marks and are overabundant in the normal tissue and blood of cancer patients. Analysis of plasma small RNAs in women who developed breast cancer indicates there is an overabundance of Y-RNA preceding neoplasm diagnosis. CONCLUSIONS: This study reveals extensive pleiotropism between blood traits and cancer risk. Pleiotropism is linked to factors and processes involved in hematopoietic development and immune system function, including components of the major histocompatibility complexes, and regulators of telomere length and myeloid lineage. Deregulation of Y-RNAs is also associated with pleiotropism. Overexpression of these elements might indicate increased cancer risk.

Clonally resolved single-cell multi-omics identifies routes of cellular differentiation in acute myeloid leukemia.

Inter-patient variability and the similarity of healthy and leukemic stem cells (LSCs) have impeded the characterization of LSCs in acute myeloid leukemia (AML) and their differentiation landscape. Here, we introduce CloneTracer, a novel method that adds clonal resolution to single-cell RNA-seq datasets. Applied to samples from 19 AML patients, CloneTracer revealed routes of leukemic differentiation. Although residual healthy and preleukemic cells dominated the dormant stem cell compartment, active LSCs resembled their healthy counterpart and retained erythroid capacity. By contrast, downstream myeloid progenitors constituted a highly aberrant, disease-defining compartment: their gene expression and differentiation state affected both the chemotherapy response and leukemia's ability to differentiate into transcriptomically normal monocytes. Finally, we demonstrated the potential of CloneTracer to identify surface markers misregulated specifically in leukemic cells. Taken together, CloneTracer reveals a differentiation landscape that mimics its healthy counterpart and may determine biology and therapy response in AML.

Disruption of stem cell niche-confined R-spondin 3 expression leads to impaired hematopoiesis.

Self-renewal and differentiation of stem and progenitor cells are tightly regulated to ensure tissue homeostasis. This regulation is enabled both remotely by systemic circulating cues, such as cytokines and hormones, and locally by various niche-confined factors. R-spondin 3 (RSPO3) is one of the most potent enhancers of Wnt signaling, and its expression is usually restricted to the stem cell niche where it provides localized enhancement of Wnt signaling to regulate stem cell expansion and differentiation. Disruption of this niche-confined expression can disturb proper tissue organization and lead to cancers. Here, we investigate the consequences of disrupting the niche-restricted expression of RSPO3 in various tissues, including the hematopoietic system. We show that normal Rspo3 expression is confined to the perivascular niche in the bone marrow. Induction of increased systemic levels of circulating RSPO3 outside of the niche results in prominent loss of early B-cell progenitors and anemia but surprisingly has no effect on hematopoietic stem cells. Using molecular, pharmacologic, and genetic approaches, we show that these RSPO3-induced hematopoietic phenotypes are Wnt and RSPO3 dependent and mediated through noncanonical Wnt signaling. Our study highlights a distinct role for a Wnt/RSPO3 signaling axis in the regulation of hematopoiesis, as well as possible challenges related to therapeutic use of RSPOs for regenerative medicine.

Cellular hierarchies predict drug response in acute myeloid leukemia.

In a recent Nature Medicine study, Zeng and colleagues integrate both genomic and stem cell models of acute myeloid leukemia (AML) by deconvoluting cellular hierarchies of more than 1,000 AML samples. This work introduces a framework capable of predicting drug responses to targeted therapies in future clinical trials.

Transcriptomic diversity in human medullary thymic epithelial cells.

The induction of central T cell tolerance in the thymus depends on the presentation of peripheral self-epitopes by medullary thymic epithelial cells (mTECs). This promiscuous gene expression (pGE) drives mTEC transcriptomic diversity, with non-canonical transcript initiation, alternative splicing, and expression of endogenous retroelements (EREs) representing important but incompletely understood contributors. Here we map the expression of genome-wide transcripts in immature and mature human mTECs using high-throughput 5' cap and RNA sequencing. Both mTEC populations show high splicing entropy, potentially driven by the expression of peripheral splicing factors. During mTEC maturation, rates of global transcript mis-initiation increase and EREs enriched in long terminal repeat retrotransposons are up-regulated, the latter often found in proximity to differentially expressed genes. As a resource, we provide an interactive public interface for exploring mTEC transcriptomic diversity. Our findings therefore help construct a map of transcriptomic diversity in the healthy human thymus and may ultimately facilitate the identification of those epitopes which contribute to autoimmunity and immune recognition of tumor antigens.

Single-cell proteo-genomic reference maps of the hematopoietic system enable the purification and massive profiling of precisely defined cell states.

Single-cell genomics technology has transformed our understanding of complex cellular systems. However, excessive cost and a lack of strategies for the purification of newly identified cell types impede their functional characterization and large-scale profiling. Here, we have generated high-content single-cell proteo-genomic reference maps of human blood and bone marrow that quantitatively link the expression of up to 197 surface markers to cellular identities and biological processes across all main hematopoietic cell types in healthy aging and leukemia. These reference maps enable the automatic design of cost-effective high-throughput cytometry schemes that outperform state-of-the-art approaches, accurately reflect complex topologies of cellular systems and permit the purification of precisely defined cell states. The systematic integration of cytometry and proteo-genomic data enables the functional capacities of precisely mapped cell states to be measured at the single-cell level. Our study serves as an accessible resource and paves the way for a data-driven era in cytometry.

Decline in IGF1 in the bone marrow microenvironment initiates hematopoietic stem cell aging.

Decline in hematopoietic stem cell (HSC) function with age underlies limited health span of our blood and immune systems. In order to preserve health into older age, it is necessary to understand the nature and timing of initiating events that cause HSC aging. By performing a cross-sectional study in mice, we discover that hallmarks of aging in HSCs and hematopoiesis begin to accumulate by middle age and that the bone marrow (BM) microenvironment at middle age induces and is indispensable for hematopoietic aging. Using unbiased approaches, we find that decreased levels of the longevity-associated molecule IGF1 in the local middle-aged BM microenvironment are a factor causing HSC aging. Direct stimulation of middle-aged HSCs with IGF1 rescues molecular and functional hallmarks of aging, including restored mitochondrial activity. Thus, although decline in IGF1 supports longevity, our work indicates that this also compromises HSC function and limits hematopoietic health span.

Identification of leukemic and pre-leukemic stem cells by clonal tracking from single-cell transcriptomics.

Cancer stem cells drive disease progression and relapse in many types of cancer. Despite this, a thorough characterization of these cells remains elusive and with it the ability to eradicate cancer at its source. In acute myeloid leukemia (AML), leukemic stem cells (LSCs) underlie mortality but are difficult to isolate due to their low abundance and high similarity to healthy hematopoietic stem cells (HSCs). Here, we demonstrate that LSCs, HSCs, and pre-leukemic stem cells can be identified and molecularly profiled by combining single-cell transcriptomics with lineage tracing using both nuclear and mitochondrial somatic variants. While mutational status discriminates between healthy and cancerous cells, gene expression distinguishes stem cells and progenitor cell populations. Our approach enables the identification of LSC-specific gene expression programs and the characterization of differentiation blocks induced by leukemic mutations. Taken together, we demonstrate the power of single-cell multi-omic approaches in characterizing cancer stem cells.

Combined single-cell and spatial transcriptomics reveal the molecular, cellular and spatial bone marrow niche organization.

The bone marrow constitutes the primary site for life-long blood production and skeletal regeneration. However, its cellular and spatial organization remains controversial. Here, we combine single-cell and spatially resolved transcriptomics to systematically map the molecular, cellular and spatial composition of distinct bone marrow niches. This allowed us to transcriptionally profile all major bone-marrow-resident cell types, determine their localization and clarify sources of pro-haematopoietic factors. Our data demonstrate that Cxcl12-abundant-reticular (CAR) cell subsets (Adipo-CAR and Osteo-CAR) differentially localize to sinusoidal and arteriolar surfaces, act locally as 'professional cytokine-secreting cells' and thereby establish peri-vascular micro-niches. Importantly, the three-dimensional bone-marrow organization can be accurately inferred from single-cell transcriptome data using the RNA-Magnet algorithm described here. Together, our study reveals the cellular and spatial organization of bone marrow niches and offers a systematic approach to dissect the complex organization of whole organs.

Human haematopoietic stem cell lineage commitment is a continuous process.

Blood formation is believed to occur through stepwise progression of haematopoietic stem cells (HSCs) following a tree-like hierarchy of oligo-, bi- and unipotent progenitors. However, this model is based on the analysis of predefined flow-sorted cell populations. Here we integrated flow cytometric, transcriptomic and functional data at single-cell resolution to quantitatively map early differentiation of human HSCs towards lineage commitment. During homeostasis, individual HSCs gradually acquire lineage biases along multiple directions without passing through discrete hierarchically organized progenitor populations. Instead, unilineage-restricted cells emerge directly from a 'continuum of low-primed undifferentiated haematopoietic stem and progenitor cells' (CLOUD-HSPCs). Distinct gene expression modules operate in a combinatorial manner to control stemness, early lineage priming and the subsequent progression into all major branches of haematopoiesis. These data reveal a continuous landscape of human steady-state haematopoiesis downstream of HSCs and provide a basis for the understanding of haematopoietic malignancies.

Inflammation-Induced Emergency Megakaryopoiesis Driven by Hematopoietic Stem Cell-like Megakaryocyte Progenitors.

Infections are associated with extensive platelet consumption, representing a high risk for health. However, the mechanism coordinating the rapid regeneration of the platelet pool during such stress conditions remains unclear. Here, we report that the phenotypic hematopoietic stem cell (HSC) compartment contains stem-like megakaryocyte-committed progenitors (SL-MkPs), a cell population that shares many features with multipotent HSCs and serves as a lineage-restricted emergency pool for inflammatory insults. During homeostasis, SL-MkPs are maintained in a primed but quiescent state, thus contributing little to steady-state megakaryopoiesis. Even though lineage-specific megakaryocyte transcripts are expressed, protein synthesis is suppressed. In response to acute inflammation, SL-MkPs become activated, resulting in megakaryocyte protein production from pre-existing transcripts and a maturation of SL-MkPs and other megakaryocyte progenitors. This results in an efficient replenishment of platelets that are lost during inflammatory insult. Thus, our study reveals an emergency machinery that counteracts life-threatening platelet depletions during acute inflammation.