Maintenance of hematopoietic stem cells by tyrosine-unphosphorylated STAT5 and JAK inhibition.

Adult haematopoietic stem cells (HSCs) are responsible for the lifelong production of blood and immune cells, a process regulated by extracellular cues including cytokines. Many cytokines signal through the conserved JAK/STAT pathway, in which tyrosine-phosphorylated STATs (pSTATs) function as transcription factors. STAT5 is a pivotal downstream mediator of several cytokines known to regulate haematopoiesis but its function in the HSC compartment remains poorly understood. Here, we show that STAT5-deficient HSCs exhibit an unusual phenotype: reduced multi-lineage repopulation and self-renewal, combined with reduced exit from quiescence and increased differentiation. This was driven not only by loss of canonical pSTAT5 signalling, but also by loss of distinct transcriptional functions mediated by STAT5 lacking canonical tyrosine phosphorylation (uSTAT5). Consistent with this concept, expression of an unphosphorylatable STAT5 mutant constrained wild-type HSC differentiation, promoted their maintenance and upregulated transcriptional programs associated with quiescence and stemness. The JAK1/2 inhibitor, ruxolitinib, which increased the uSTAT5:pSTAT5 ratio, had similar effects on murine HSC function: it constrained HSC differentiation and proliferation, promoted HSC maintenance and upregulated transcriptional programs associated with stemness. Ruxolitinib also enhanced serial replating of normal human HSPCs, CALR-mutant murine HSCs and HSPCs obtained from patients with myelofibrosis. Our results therefore reveal a previously unrecognized interplay between pSTAT5 and uSTAT5 in the control of HSC function and highlight JAK inhibition as a potential strategy for enhancing HSC function during ex vivo culture. Increased levels of uSTAT5 may also contribute to the failure of JAK inhibitors to eradicate myeloproliferative neoplasms.

Inflammation perturbs hematopoiesis by remodeling specific compartments of the bone marrow niche.

Hematopoietic stem and progenitor cells (HSPC) are regulated by interactions with stromal cells in the bone marrow (BM) cavity, which can be segregated into two spatially defined central marrow (CM) and endosteal (Endo) compartments. However, the importance of this spatial compartmentalization for BM responses to inflammation and neoplasia remains largely unknown. Here, we extensively validate a combination of scRNA-seq profiling and matching flow cytometry isolation that reproducibly identifies 7 key CM and Endo populations across mouse strains and accurately surveys both niche locations. We demonstrate that different perturbations exert specific effects on different compartments, with type I interferon responses causing CM mesenchymal stromal cells to adopt an inflammatory phenotype associated with overproduction of chemokines modulating local monocyte dynamics in the surrounding microenvironment. Our results provide a comprehensive method for molecular and functional stromal characterization and highlight the importance of altered stomal cell activity in regulating hematopoietic responses to inflammatory challenges.

Axin1 and Axin2 regulate the WNT-signaling landscape to promote distinct mesoderm programs.

How distinct mesodermal lineages - extraembryonic, lateral, intermediate, paraxial and axial - are specified from pluripotent epiblast during gastrulation is a longstanding open question. By investigating AXIN, a negative regulator of the WNT/ß-catenin pathway, we have uncovered new roles for WNT signaling in the determination of mesodermal fates. We undertook complementary approaches to dissect the role of WNT signaling that augmented a detailed analysis of Axin1;Axin2 mutant mouse embryos, including single-cell and single-embryo transcriptomics, with in vitro pluripotent Epiblast-Like Cell differentiation assays. This strategy allowed us to reveal two layers of regulation. First, WNT initiates differentiation of primitive streak cells into mesoderm progenitors, and thereafter, WNT amplifies and cooperates with BMP/pSMAD1/5/9 or NODAL/pSMAD2/3 to propel differentiating mesoderm progenitors into either posterior streak derivatives or anterior streak derivatives, respectively. We propose that Axin1 and Axin2 prevent aberrant differentiation of pluripotent epiblast cells into mesoderm by spatially and temporally regulating WNT signaling levels.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Adaptation to ex vivo culture reduces human hematopoietic stem cell activity independently of the cell cycle.

ABSTRACT: Loss of long-term hematopoietic stem cell (LT-HSC) function ex vivo hampers the success of clinical protocols that rely on culture. However, the kinetics and mechanisms through which this occurs remain incompletely characterized. In this study, through time-resolved single-cell RNA sequencing, matched in vivo functional analysis, and the use of a reversible in vitro system of early G1 arrest, we defined the sequence of transcriptional and functional events that occur during the first ex vivo division of human LT-HSCs. We demonstrated that the sharpest loss in LT-HSC repopulation capacity happens early on, between 6 and 24 hours of culture, before LT-HSCs commit to cell cycle progression. During this time window, LT-HSCs adapt to the culture environment, limit the global variability in gene expression, and transiently upregulate gene networks involved in signaling and stress responses. From 24 hours, LT-HSC progression past early G1 contributes to the establishment of differentiation programs in culture. However, contrary to the current assumptions, we demonstrated that the loss of HSC function ex vivo is independent of cell cycle progression. Finally, we showed that targeting LT-HSC adaptation to culture by inhibiting the early activation of JAK/STAT signaling improves HSC long-term repopulating function ex vivo. Collectively, our study demonstrated that controlling early LT-HSC adaptation to ex vivo culture, for example, via JAK inhibition, is critically important to improve HSC gene therapy and expansion protocols.

Autophagy counters inflammation-driven glycolytic impairment in aging hematopoietic stem cells.

Autophagy is central to the benefits of longevity signaling programs and to hematopoietic stem cell (HSC) response to nutrient stress. With age, a subset of HSCs increases autophagy flux and preserves regenerative capacity, but the signals triggering autophagy and maintaining the functionality of autophagy-activated old HSCs (oHSCs) remain unknown. Here, we demonstrate that autophagy is an adaptive cytoprotective response to chronic inflammation in the aging murine bone marrow (BM) niche. We find that inflammation impairs glucose uptake and suppresses glycolysis in oHSCs through Socs3-mediated inhibition of AKT/FoxO-dependent signaling, with inflammation-mediated autophagy engagement preserving functional quiescence by enabling metabolic adaptation to glycolytic impairment. Moreover, we show that transient autophagy induction via a short-term fasting/refeeding paradigm normalizes glycolytic flux and significantly boosts oHSC regenerative potential. Our results identify inflammation-driven glucose hypometabolism as a key driver of HSC dysfunction with age and establish autophagy as a targetable node to reset oHSC regenerative capacity.

Iron dysregulation and inflammatory stress erythropoiesis associates with long-term outcome of COVID-19.

Persistent symptoms following SARS-CoV-2 infection are increasingly reported, although the drivers of post-acute sequelae (PASC) of COVID-19 are unclear. Here we assessed 214 individuals infected with SARS-CoV-2, with varying disease severity, for one year from COVID-19 symptom onset to determine the early correlates of PASC. A multivariate signature detected beyond two weeks of disease, encompassing unresolving inflammation, anemia, low serum iron, altered iron-homeostasis gene expression and emerging stress erythropoiesis; differentiated those who reported PASC months later, irrespective of COVID-19 severity. A whole-blood heme-metabolism signature, enriched in hospitalized patients at month 1-3 post onset, coincided with pronounced iron-deficient reticulocytosis. Lymphopenia and low numbers of dendritic cells persisted in those with PASC, and single-cell analysis reported iron maldistribution, suggesting monocyte iron loading and increased iron demand in proliferating lymphocytes. Thus, defects in iron homeostasis, dysregulated erythropoiesis and immune dysfunction due to COVID-19 possibly contribute to inefficient oxygen transport, inflammatory disequilibrium and persisting symptomatology, and may be therapeutically tractable.

Cis inhibition of NOTCH1 through JAGGED1 sustains embryonic hematopoietic stem cell fate.

Hematopoietic stem cells (HSCs) develop from the hemogenic endothelium (HE) in the aorta- gonads-and mesonephros (AGM) region and reside within Intra-aortic hematopoietic clusters (IAHC) along with hematopoietic progenitors (HPC). The signalling mechanisms that distinguish HSCs from HPCs are unknown. Notch signaling is essential for arterial specification, IAHC formation and HSC activity, but current studies on how Notch segregates these different fates are inconsistent. We now demonstrate that Notch activity is highest in a subset of, GFI1 + , HSC-primed HE cells, and is gradually lost with HSC maturation. We uncover that the HSC phenotype is maintained due to increasing levels of NOTCH1 and JAG1 interactions on the surface of the same cell (cis) that renders the NOTCH1 receptor from being activated. Forced activation of the NOTCH1 receptor in IAHC activates a hematopoietic differentiation program. Our results indicate that NOTCH1-JAG1 cis-inhibition preserves the HSC phenotype in the hematopoietic clusters of the embryonic aorta.

A time- and single-cell-resolved model of murine bone marrow hematopoiesis.

The paradigmatic hematopoietic tree model is increasingly recognized to be limited, as it is based on heterogeneous populations largely defined by non-homeostatic assays testing cell fate potentials. Here, we combine persistent labeling with time-series single-cell RNA sequencing to build a real-time, quantitative model of in vivo tissue dynamics for murine bone marrow hematopoiesis. We couple cascading single-cell expression patterns with dynamic changes in differentiation and growth speeds. The resulting explicit linkage between molecular states and cellular behavior reveals widely varying self-renewal and differentiation properties across distinct lineages. Transplanted stem cells show strong acceleration of differentiation at specific stages of erythroid and neutrophil production, illustrating how the model can quantify the impact of perturbations. Our reconstruction of dynamic behavior from snapshot measurements is akin to how a kinetoscope allows sequential images to merge into a movie. We posit that this approach is generally applicable to understanding tissue-scale dynamics at high resolution.

DNA methylation restricts coordinated germline and neural fates in embryonic stem cell differentiation.

As embryonic stem cells (ESCs) transition from naive to primed pluripotency during early mammalian development, they acquire high DNA methylation levels. During this transition, the germline is specified and undergoes genome-wide DNA demethylation, while emergence of the three somatic germ layers is preceded by acquisition of somatic DNA methylation levels in the primed epiblast. DNA methylation is essential for embryogenesis, but the point at which it becomes critical during differentiation and whether all lineages equally depend on it is unclear. Here, using culture modeling of cellular transitions, we found that DNA methylation-free mouse ESCs with triple DNA methyltransferase knockout (TKO) progressed through the continuum of pluripotency states but demonstrated skewed differentiation abilities toward neural versus other somatic lineages. More saliently, TKO ESCs were fully competent for establishing primordial germ cell-like cells, even showing temporally extended and self-sustained capacity for the germline fate. By mapping chromatin states, we found that neural and germline lineages are linked by a similar enhancer dynamic upon exit from the naive state, defined by common sets of transcription factors, including methyl-sensitive ones, that fail to be decommissioned in the absence of DNA methylation. We propose that DNA methylation controls the temporality of a coordinated neural-germline axis of the preferred differentiation route during early development.

Tracking early mammalian organogenesis - prediction and validation of differentiation trajectories at whole organism scale.

Early organogenesis represents a key step in animal development, during which pluripotent cells diversify to initiate organ formation. Here, we sampled 300,000 single-cell transcriptomes from mouse embryos between E8.5 and E9.5 in 6-h intervals and combined this new dataset with our previous atlas (E6.5-E8.5) to produce a densely sampled timecourse of >400,000 cells from early gastrulation to organogenesis. Computational lineage reconstruction identified complex waves of blood and endothelial development, including a new programme for somite-derived endothelium. We also dissected the E7.5 primitive streak into four adjacent regions, performed scRNA-seq and predicted cell fates computationally. Finally, we defined developmental state/fate relationships by combining orthotopic grafting, microscopic analysis and scRNA-seq to transcriptionally determine cell fates of grafted primitive streak regions after 24 h of in vitro embryo culture. Experimentally determined fate outcomes were in good agreement with computationally predicted fates, demonstrating how classical grafting experiments can be revisited to establish high-resolution cell state/fate relationships. Such interdisciplinary approaches will benefit future studies in developmental biology and guide the in vitro production of cells for organ regeneration and repair.

Single-cell multi-omics analysis of COVID-19 patients with pre-existing autoimmune diseases shows aberrant immune responses to infection.

In COVID-19, hyperinflammatory and dysregulated immune responses contribute to severity. Patients with pre-existing autoimmune conditions can therefore be at increased risk of severe COVID-19 and/or associated sequelae, yet SARS-CoV-2 infection in this group has been little studied. Here, we performed single-cell analysis of peripheral blood mononuclear cells from patients with three major autoimmune diseases (rheumatoid arthritis, psoriasis, or multiple sclerosis) during SARS-CoV-2 infection. We observed compositional differences between the autoimmune disease groups coupled with altered patterns of gene expression, transcription factor activity, and cell-cell communication that substantially shape the immune response under SARS-CoV-2 infection. While enrichment of HLA-DRlow CD14+ monocytes was observed in all three autoimmune disease groups, type-I interferon signaling as well as inflammatory T cell and monocyte responses varied widely between the three groups of patients. Our results reveal disturbed immune responses to SARS-CoV-2 in patients with pre-existing autoimmunity, highlighting important considerations for disease treatment and follow-up.

Expression profiling of cerebrospinal fluid identifies dysregulated antiviral mechanisms in multiple sclerosis.

Despite the overwhelming evidence that multiple sclerosis is an autoimmune disease, relatively little is known about the precise nature of the immune dysregulation underlying the development of the disease. Reasoning that the CSF from patients might be enriched for cells relevant in pathogenesis, we have completed a high-resolution single-cell analysis of 96 732 CSF cells collected from 33 patients with multiple sclerosis (n = 48 675) and 48 patients with other neurological diseases (n = 48 057). Completing comprehensive cell type annotation, we identified a rare population of CD8+ T cells, characterized by the upregulation of inhibitory receptors, increased in patients with multiple sclerosis. Applying a Multi-Omics Factor Analysis to these single-cell data further revealed that activity in pathways responsible for controlling inflammatory and type 1 interferon responses are altered in multiple sclerosis in both T cells and myeloid cells. We also undertook a systematic search for expression quantitative trait loci in the CSF cells. Of particular interest were two expression quantitative trait loci in CD8+ T cells that were fine mapped to multiple sclerosis susceptibility variants in the viral control genes ZC3HAV1 (rs10271373) and IFITM2 (rs1059091). Further analysis suggests that these associations likely reflect genetic effects on RNA splicing and cell-type specific gene expression respectively. Collectively, our study suggests that alterations in viral control mechanisms might be important in the development of multiple sclerosis.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes.

Knowledge of locations and activities of cis -regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our V al i dated S ystematic I ntegrati on (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Preleukemic single-cell landscapes reveal mutation-specific mechanisms and gene programs predictive of AML patient outcomes.

Acute myeloid leukemia (AML) and myeloid neoplasms develop through acquisition of somatic mutations that confer mutation-specific fitness advantages to hematopoietic stem and progenitor cells. However, our understanding of mutational effects remains limited to the resolution attainable within immunophenotypically and clinically accessible bulk cell populations. To decipher heterogeneous cellular fitness to preleukemic mutational perturbations, we performed single-cell RNA sequencing of eight different mouse models with driver mutations of myeloid malignancies, generating 269,048 single-cell profiles. Our analysis infers mutation-driven perturbations in cell abundance, cellular lineage fate, cellular metabolism, and gene expression at the continuous resolution, pinpointing cell populations with transcriptional alterations associated with differentiation bias. We further develop an 11-gene scoring system (Stem11) on the basis of preleukemic transcriptional signatures that predicts AML patient outcomes. Our results demonstrate that a single-cell-resolution deep characterization of preleukemic biology has the potential to enhance our understanding of AML heterogeneity and inform more effective risk stratification strategies.

Autophagy counters inflammation-driven glycolytic impairment in aging hematopoietic stem cells.

Aging of the hematopoietic system promotes various blood, immune and systemic disorders and is largely driven by hematopoietic stem cell (HSC) dysfunction ( 1 ). Autophagy is central for the benefits associated with activation of longevity signaling programs ( 2 ), and for HSC function and response to nutrient stress ( 3,4 ). With age, a subset of HSCs increases autophagy flux and preserves some regenerative capacity, while the rest fail to engage autophagy and become metabolically overactivated and dysfunctional ( 4 ). However, the signals that promote autophagy in old HSCs and the mechanisms responsible for the increased regenerative potential of autophagy-activated old HSCs remain unknown. Here, we demonstrate that autophagy activation is an adaptive survival response to chronic inflammation in the aging bone marrow (BM) niche ( 5 ). We find that inflammation impairs glucose metabolism and suppresses glycolysis in aged HSCs through Socs3-mediated impairment of AKT/FoxO-dependent signaling. In this context, we show that inflammation-mediated autophagy engagement preserves functional quiescence by enabling metabolic adaptation to glycolytic impairment. Moreover, we demonstrate that transient autophagy induction via a short-term fasting/refeeding paradigm normalizes glucose uptake and glycolytic flux and significantly improves old HSC regenerative potential. Our results identify inflammation-driven glucose hypometabolism as a key driver of HSC dysfunction with age and establish autophagy as a targetable node to reset old HSC glycolytic and regenerative capacity. One-Sentence Summary: Autophagy compensates for chronic inflammation-induced metabolic deregulation in old HSCs, and its transient modulation can reset old HSC glycolytic and regenerative capacity.

Genome-wide transcription factor-binding maps reveal cell-specific changes in the regulatory architecture of human HSPCs.

Hematopoietic stem and progenitor cells (HSPCs) rely on a complex interplay among transcription factors (TFs) to regulate differentiation into mature blood cells. A heptad of TFs (FLI1, ERG, GATA2, RUNX1, TAL1, LYL1, LMO2) bind regulatory elements in bulk CD34+ HSPCs. However, whether specific heptad-TF combinations have distinct roles in regulating hematopoietic differentiation remains unknown. We mapped genome-wide chromatin contacts (HiC, H3K27ac, HiChIP), chromatin modifications (H3K4me3, H3K27ac, H3K27me3) and 10 TF binding profiles (heptad, PU.1, CTCF, STAG2) in HSPC subsets (stem/multipotent progenitors plus common myeloid, granulocyte macrophage, and megakaryocyte erythrocyte progenitors) and found TF occupancy and enhancer-promoter interactions varied significantly across cell types and were associated with cell-type-specific gene expression. Distinct regulatory elements were enriched with specific heptad-TF combinations, including stem-cell-specific elements with ERG, and myeloid- and erythroid-specific elements with combinations of FLI1, RUNX1, GATA2, TAL1, LYL1, and LMO2. Furthermore, heptad-occupied regions in HSPCs were subsequently bound by lineage-defining TFs, including PU.1 and GATA1, suggesting that heptad factors may prime regulatory elements for use in mature cell types. We also found that enhancers with cell-type-specific heptad occupancy shared a common grammar with respect to TF binding motifs, suggesting that combinatorial binding of TF complexes was at least partially regulated by features encoded in DNA sequence motifs. Taken together, this study comprehensively characterizes the gene regulatory landscape in rare subpopulations of human HSPCs. The accompanying data sets should serve as a valuable resource for understanding adult hematopoiesis and a framework for analyzing aberrant regulatory networks in leukemic cells.

Mapping interindividual dynamics of innate immune response at single-cell resolution.

Common genetic variants across individuals modulate the cellular response to pathogens and are implicated in diverse immune pathologies, yet how they dynamically alter the response upon infection is not well understood. Here, we triggered antiviral responses in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-sequencing. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), a statistical approach designed to identify nonlinear dynamic genetic effects across transcriptional trajectories of cells. This approach identified 1,275 expression quantitative trait loci (local false discovery rate 10%) that manifested during the responses, many of which were colocalized with susceptibility loci identified by genome-wide association studies of infectious and autoimmune diseases, including the OAS1 splicing quantitative trait locus in a COVID-19 susceptibility locus. In summary, our analytical approach provides a unique framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

An atlas of rabbit development as a model for single-cell comparative genomics.

Traditionally, the mouse has been the favoured vertebrate model for biomedical research, due to its experimental and genetic tractability. However, non-rodent embryological studies highlight that many aspects of early mouse development, such as its egg-cylinder gastrulation and method of implantation, diverge from other mammals, thus complicating inferences about human development. Like the human embryo, rabbits develop as a flat-bilaminar disc. Here we constructed a morphological and molecular atlas of rabbit development. We report transcriptional and chromatin accessibility profiles for over 180,000 single cells and high-resolution histology sections from embryos spanning gastrulation, implantation, amniogenesis and early organogenesis. Using a neighbourhood comparison pipeline, we compare the transcriptional landscape of rabbit and mouse at the scale of the entire organism. We characterize the gene regulatory programmes underlying trophoblast differentiation and identify signalling interactions involving the yolk sac mesothelium during haematopoiesis. We demonstrate how the combination of both rabbit and mouse atlases can be leveraged to extract new biological insights from sparse macaque and human data. The datasets and computational pipelines reported here set a framework for a broader cross-species approach to decipher early mammalian development, and are readily adaptable to deploy single-cell comparative genomics more broadly across biomedical research.