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.

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.

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.

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.

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.

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.

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.

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.

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.

Proteomic analysis of circulating immune cells identifies cellular phenotypes associated with COVID-19 severity.

Certain serum proteins, including C-reactive protein (CRP) and D-dimer, have prognostic value in patients with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Nonetheless, these factors are non-specific, providing limited mechanistic insight into the peripheral blood mononuclear cell (PBMC) populations that drive the pathogenesis of severe COVID-19. To identify cellular phenotypes associated with disease, we performed a comprehensive, unbiased analysis of total and plasma-membrane PBMC proteomes from 40 unvaccinated individuals with SARS-CoV-2, spanning the whole disease spectrum. Combined with RNA sequencing (RNA-seq) and flow cytometry from the same donors, we define a comprehensive multi-omic profile for each severity level, revealing that immune-cell dysregulation progresses with increasing disease. The cell-surface proteins CEACAMs1, 6, and 8, CD177, CD63, and CD89 are strongly associated with severe COVID-19, corresponding to the emergence of atypical CD3+CD4+CEACAM1/6/8+CD177+CD63+CD89+ and CD16+CEACAM1/6/8+ mononuclear cells. Utilization of these markers may facilitate real-time patient assessment by flow cytometry and identify immune populations that could be targeted to ameliorate immunopathology.

Genotoxic aldehyde stress prematurely ages hematopoietic stem cells in a p53-driven manner.

Aged hematopoietic stem cells (HSCs) display diminished self-renewal and a myeloid differentiation bias. However, the drivers and mechanisms that underpin this fundamental switch are not understood. HSCs produce genotoxic formaldehyde that requires protection by the detoxification enzymes ALDH2 and ADH5 and the Fanconi anemia (FA) DNA repair pathway. We find that the HSCs in young Aldh2-/-Fancd2-/- mice harbor a transcriptomic signature equivalent to aged wild-type HSCs, along with increased epigenetic age, telomere attrition, and myeloid-biased differentiation quantified by single HSC transplantation. In addition, the p53 response is vigorously activated in Aldh2-/-Fancd2-/- HSCs, while p53 deletion rescued this aged HSC phenotype. To further define the origins of the myeloid differentiation bias, we use a GFP genetic reporter to find a striking enrichment of Vwf+ myeloid and megakaryocyte-lineage-biased HSCs. These results indicate that metabolism-derived formaldehyde-DNA damage stimulates the p53 response in HSCs to drive accelerated aging.

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.

Physioxia improves the selectivity of hematopoietic stem cell expansion cultures.

Hematopoietic stem cells (HSCs) are a rare type of hematopoietic cell that can entirely reconstitute the blood and immune system after transplantation. Allogeneic HSC transplantation (HSCT) is used clinically as a curative therapy for a range of hematolymphoid diseases; however, it remains a high-risk therapy because of its potential side effects, including poor graft function and graft-versus-host disease (GVHD). Ex vivo HSC expansion has been suggested as an approach to improve hematopoietic reconstitution in low-cell dose grafts. Here, we demonstrate that the selectivity of polyvinyl alcohol (PVA)-based mouse HSC cultures can be improved using physioxic culture conditions. Single-cell transcriptomic analysis helped confirm the inhibition of lineage-committed progenitor cells in physioxic cultures. Long-term physioxic expansion also afforded culture-based ex vivo HSC selection from whole bone marrow, spleen, and embryonic tissues. Furthermore, we provide evidence that HSC-selective ex vivo cultures deplete GVHD-causing T cells and that this approach can be combined with genotoxic-free antibody-based conditioning HSCT approaches. Our results offer a simple approach to improve PVA-based HSC cultures and the underlying molecular phenotype, and highlight the potential translational implications of selective HSC expansion systems for allogeneic HSCT.

Stromal niche inflammation mediated by IL-1 signalling is a targetable driver of haematopoietic ageing.

Haematopoietic ageing is marked by a loss of regenerative capacity and skewed differentiation from haematopoietic stem cells (HSCs), leading to impaired blood production. Signals from the bone marrow niche tailor blood production, but the contribution of the old niche to haematopoietic ageing remains unclear. Here we characterize the inflammatory milieu that drives both niche and haematopoietic remodelling. We find decreased numbers and functionality of osteoprogenitors at the endosteum and expansion of central marrow LepR+ mesenchymal stromal cells associated with deterioration of the sinusoidal vasculature. Together, they create a degraded and inflamed old bone marrow niche. Niche inflammation in turn drives the chronic activation of emergency myelopoiesis pathways in old HSCs and multipotent progenitors, which promotes myeloid differentiation and hinders haematopoietic regeneration. Moreover, we show how production of interleukin-1ß (IL-1ß) by the damaged endosteum acts in trans to drive the proinflammatory nature of the central marrow, with damaging consequences for the old blood system. Notably, niche deterioration, HSC dysfunction and defective regeneration can all be ameliorated by blocking IL-1 signalling. Our results demonstrate that targeting IL-1 as a key mediator of niche inflammation is a tractable strategy to improve blood production during ageing.

Relaxin/insulin-like family peptide receptor 4 (Rxfp4) expressing hypothalamic neurons modulate food intake and preference in mice.

OBJECTIVE: Insulin-like peptide 5 (INSL5) signalling, through its cognate receptor relaxin/insulin-like family peptide receptor 4 (RXFP4), has been reported to be orexigenic, and the high fat diet (HFD) preference observed in wildtype mice is altered in Rxfp4 knock-out mice. In this study, we used a new Rxfp4-Cre mouse model to investigate the mechanisms underlying these observations. METHODS: We generated transgenic Rxfp4-Cre mice and investigated central expression of Rxfp4 by RT-qPCR, RNAscope and intraparenchymal infusion of INSL5. Rxfp4-expressing cells were chemogenetically manipulated in global Cre-reporter mice using designer receptors exclusively activated by designer drugs (DREADDs) or after stereotactic injection of a Cre-dependent AAV-DIO-Dq-DREADD targeting a population located in the ventromedial hypothalamus (RXFP4VMH). Food intake and feeding motivation were assessed in the presence and absence of a DREADD agonist. Rxfp4-expressing cells in the hypothalamus were characterised by single-cell RNA-sequencing (scRNAseq) and the connectivity of RXFP4VMH cells was investigated using viral tracing. RESULTS: Rxfp4-Cre mice displayed Cre-reporter expression in the hypothalamus. Active expression of Rxfp4 in the adult mouse brain was confirmed by RT-qPCR and RNAscope. Functional receptor expression was supported by cyclic AMP-responses to INSL5 application in ex vivo brain slices and increased HFD and highly palatable liquid meal (HPM), but not chow, intake after intra-VMH INSL5 infusion. scRNAseq of hypothalamic RXFP4 neurons defined a cluster expressing VMH markers, alongside known appetite-modulating neuropeptide receptors (Mc4r, Cckar and Nmur2). Viral tracing demonstrated RXFP4VMH neural projections to nuclei implicated in hedonic feeding behaviour. Whole body chemogenetic inhibition (Di-DREADD) of Rxfp4-expressing cells, mimicking physiological INSL5-RXFP4 Gi-signalling, increased intake of the HFD and HPM, but not chow, whilst activation (Dq-DREADD), either at whole body level or specifically within the VMH, reduced HFD and HPM intake and motivation to work for the HPM. CONCLUSION: These findings identify RXFP4VMH neurons as regulators of food intake and preference, and hypothalamic RXFP4 signalling as a target for feeding behaviour manipulation.