Longitudinal analysis reveals that delayed bystander CD8+ T cell activation and early immune pathology distinguish severe COVID-19 from mild disease.

The kinetics of the immune changes in COVID-19 across severity groups have not been rigorously assessed. Using immunophenotyping, RNA sequencing, and serum cytokine analysis, we analyzed serial samples from 207 SARS-CoV2-infected individuals with a range of disease severities over 12 weeks from symptom onset. An early robust bystander CD8+ T cell immune response, without systemic inflammation, characterized asymptomatic or mild disease. Hospitalized individuals had delayed bystander responses and systemic inflammation that was already evident near symptom onset, indicating that immunopathology may be inevitable in some individuals. Viral load did not correlate with this early pathological response but did correlate with subsequent disease severity. Immune recovery is complex, with profound persistent cellular abnormalities in severe disease correlating with altered inflammatory responses, with signatures associated with increased oxidative phosphorylation replacing those driven by cytokines tumor necrosis factor (TNF) and interleukin (IL)-6. These late immunometabolic and immune defects may have clinical implications.

Local and systemic responses to SARS-CoV-2 infection in children and adults.

It is not fully understood why COVID-19 is typically milder in children1-3. Here, to examine the differences between children and adults in their response to SARS-CoV-2 infection, we analysed paediatric and adult patients with COVID-19 as well as healthy control individuals (total n = 93) using single-cell multi-omic profiling of matched nasal, tracheal, bronchial and blood samples. In the airways of healthy paediatric individuals, we observed cells that were already in an interferon-activated state, which after SARS-CoV-2 infection was further induced especially in airway immune cells. We postulate that higher paediatric innate interferon responses restrict viral replication and disease progression. The systemic response in children was characterized by increases in naive lymphocytes and a depletion of natural killer cells, whereas, in adults, cytotoxic T cells and interferon-stimulated subpopulations were significantly increased. We provide evidence that dendritic cells initiate interferon signalling in early infection, and identify epithelial cell states associated with COVID-19 and age. Our matching nasal and blood data show a strong interferon response in the airways with the induction of systemic interferon-stimulated populations, which were substantially reduced in paediatric patients. Together, we provide several mechanisms that explain the milder clinical syndrome observed in children.

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.

A single-cell molecular map of mouse gastrulation and early organogenesis.

Across the animal kingdom, gastrulation represents a key developmental event during which embryonic pluripotent cells diversify into lineage-specific precursors that will generate the adult organism. Here we report the transcriptional profiles of 116,312 single cells from mouse embryos collected at nine sequential time points ranging from 6.5 to 8.5 days post-fertilization. We construct a molecular map of cellular differentiation from pluripotency towards all major embryonic lineages, and explore the complex events involved in the convergence of visceral and primitive streak-derived endoderm. Furthermore, we use single-cell profiling to show that Tal1-/- chimeric embryos display defects in early mesoderm diversification, and we thus demonstrate how combining temporal and transcriptional information can illuminate gene function. Together, this comprehensive delineation of mammalian cell differentiation trajectories in vivo represents a baseline for understanding the effects of gene mutations during development, as well as a roadmap for the optimization of in vitro differentiation protocols for regenerative medicine.

Notch ligand Dll4 impairs cell recruitment to aortic clusters and limits blood stem cell generation.

Hematopoietic stem cells (HSCs) develop from the hemogenic endothelium in cluster structures that protrude into the embryonic aortic lumen. Although much is known about the molecular characteristics of the developing hematopoietic cells, we lack a complete understanding of their origin and the three-dimensional organization of the niche. Here, we use advanced live imaging techniques of organotypic slice cultures, clonal analysis, and mathematical modeling to show the two-step process of intra-aortic hematopoietic cluster (IACH) formation. First, a hemogenic progenitor buds up from the endothelium and undergoes division forming the monoclonal core of the IAHC. Next, surrounding hemogenic cells are recruited into the IAHC, increasing their size and heterogeneity. We identified the Notch ligand Dll4 as a negative regulator of the recruitment phase of IAHC. Blocking of Dll4 promotes the entrance of new hemogenic Gfi1+ cells into the IAHC and increases the number of cells that acquire HSC activity. Mathematical modeling based on our data provides estimation of the cluster lifetime and the average recruitment time of hemogenic cells to the cluster under physiologic and Dll4-inhibited conditions.

Myeloid progenitor cluster formation drives emergency and leukaemic myelopoiesis.

Although many aspects of blood production are well understood, the spatial organization of myeloid differentiation in the bone marrow remains unknown. Here we use imaging to track granulocyte/macrophage progenitor (GMP) behaviour in mice during emergency and leukaemic myelopoiesis. In the steady state, we find individual GMPs scattered throughout the bone marrow. During regeneration, we observe expanding GMP patches forming defined GMP clusters, which, in turn, locally differentiate into granulocytes. The timed release of important bone marrow niche signals (SCF, IL-1ß, G-CSF, TGFß and CXCL4) and activation of an inducible Irf8 and ß-catenin progenitor self-renewal network control the transient formation of regenerating GMP clusters. In leukaemia, we show that GMP clusters are constantly produced owing to persistent activation of the self-renewal network and a lack of termination cytokines that normally restore haematopoietic stem-cell quiescence. Our results uncover a previously unrecognized dynamic behaviour of GMPs in situ, which tunes emergency myelopoiesis and is hijacked in leukaemia.

CHD7 and Runx1 interaction provides a braking mechanism for hematopoietic differentiation.

Hematopoietic stem and progenitor cell (HSPC) formation and lineage differentiation involve gene expression programs orchestrated by transcription factors and epigenetic regulators. Genetic disruption of the chromatin remodeler chromodomain-helicase-DNA-binding protein 7 (CHD7) expanded phenotypic HSPCs, erythroid, and myeloid lineages in zebrafish and mouse embryos. CHD7 acts to suppress hematopoietic differentiation. Binding motifs for RUNX and other hematopoietic transcription factors are enriched at sites occupied by CHD7, and decreased RUNX1 occupancy correlated with loss of CHD7 localization. CHD7 physically interacts with RUNX1 and suppresses RUNX1-induced expansion of HSPCs during development through modulation of RUNX1 activity. Consequently, the RUNX1:CHD7 axis provides proper timing and function of HSPCs as they emerge during hematopoietic development or mature in adults, representing a distinct and evolutionarily conserved control mechanism to ensure accurate hematopoietic lineage differentiation.

Assessing the reliability of spike-in normalization for analyses of single-cell RNA sequencing data.

By profiling the transcriptomes of individual cells, single-cell RNA sequencing provides unparalleled resolution to study cellular heterogeneity. However, this comes at the cost of high technical noise, including cell-specific biases in capture efficiency and library generation. One strategy for removing these biases is to add a constant amount of spike-in RNA to each cell and to scale the observed expression values so that the coverage of spike-in transcripts is constant across cells. This approach has previously been criticized as its accuracy depends on the precise addition of spike-in RNA to each sample. Here, we perform mixture experiments using two different sets of spike-in RNA to quantify the variance in the amount of spike-in RNA added to each well in a plate-based protocol. We also obtain an upper bound on the variance due to differences in behavior between the two spike-in sets. We demonstrate that both factors are small contributors to the total technical variance and have only minor effects on downstream analyses, such as detection of highly variable genes and clustering. Our results suggest that scaling normalization using spike-in transcripts is reliable enough for routine use in single-cell RNA sequencing data analyses.

Enhanced hemato-endothelial specification during human embryonic differentiation through developmental cooperation between AF4-MLL and MLL-AF4 fusions.

The t(4;11)(q21;q23) translocation is associated with high-risk infant pro-B-cell acute lymphoblastic leukemia and arises prenatally during embryonic/fetal hematopoiesis. The developmental/pathogenic contribution of the t(4;11)-resulting MLL-AF4 (MA4) and AF4-MLL (A4M) fusions remains unclear; MA4 is always expressed in patients with t(4;11)+ B-cell acute lymphoblastic leukemia, but the reciprocal fusion A4M is expressed in only half of the patients. Because prenatal leukemogenesis manifests as impaired early hematopoietic differentiation, we took advantage of well-established human embryonic stem cell-based hematopoietic differentiation models to study whether the A4M fusion cooperates with MA4 during early human hematopoietic development. Co-expression of A4M and MA4 strongly promoted the emergence of hemato-endothelial precursors, both endothelial- and hemogenic-primed. Double fusion-expressing hemato-endothelial precursors specified into significantly higher numbers of both hematopoietic and endothelial-committed cells, irrespective of the differentiation protocol used and without hijacking survival/proliferation. Functional analysis of differentially expressed genes and differentially enriched H3K79me3 genomic regions by RNA-sequencing and H3K79me3 chromatin immunoprecipitation-sequencing, respectively, confirmed a hematopoietic/endothelial cell differentiation signature in double fusion-expressing hemato-endothelial precursors. Importantly, chromatin immunoprecipitation-sequencing analysis revealed a significant enrichment of H3K79 methylated regions specifically associated with HOX-A cluster genes in double fusion-expressing differentiating hematopoietic cells. Overall, these results establish a functional and molecular cooperation between MA4 and A4M fusions during human hematopoietic development.

Key regulators control distinct transcriptional programmes in blood progenitor and mast cells.

Despite major advances in the generation of genome-wide binding maps, the mechanisms by which transcription factors (TFs) regulate cell type identity have remained largely obscure. Through comparative analysis of 10 key haematopoietic TFs in both mast cells and blood progenitors, we demonstrate that the largely cell type-specific binding profiles are not opportunistic, but instead contribute to cell type-specific transcriptional control, because (i) mathematical modelling of differential binding of shared TFs can explain differential gene expression, (ii) consensus binding sites are important for cell type-specific binding and (iii) knock-down of blood stem cell regulators in mast cells reveals mast cell-specific genes as direct targets. Finally, we show that the known mast cell regulators Mitf and c-fos likely contribute to the global reorganisation of TF binding profiles. Taken together therefore, our study elucidates how key regulatory TFs contribute to transcriptional programmes in several distinct mammalian cell types.

Integrated genome-scale analysis of the transcriptional regulatory landscape in a blood stem/progenitor cell model.

Comprehensive study of transcriptional control processes will be required to enhance our understanding of both normal and malignant hematopoiesis. Modern sequencing technologies have revolutionized our ability to generate genome-scale expression and histone modification profiles, transcription factor (TF)-binding maps, and also comprehensive chromatin-looping information. Many of these technologies, however, require large numbers of cells, and therefore cannot be applied to rare hematopoietic stem/progenitor cell (HSPC) populations. The stem cell factor-dependent multipotent progenitor cell line HPC-7 represents a well-recognized cell line model for HSPCs. Here we report genome-wide maps for 17 TFs, 3 histone modifications, DNase I hypersensitive sites, and high-resolution promoter-enhancer interactomes in HPC-7 cells. Integrated analysis of these complementary data sets revealed TF occupancy patterns of genomic regions involved in promoter-anchored loops. Moreover, preferential associations between pairs of TFs bound at either ends of chromatin loops led to the identification of 4 previously unrecognized protein-protein interactions between key blood stem cell regulators. All HPC-7 data sets are freely available both through standard repositories and a user-friendly Web interface. Together with previously generated genome-wide data sets, this study integrates HPC-7 data into a genomic resource on par with ENCODE tier 1 cell lines and, importantly, is the only current model with comprehensive genome-scale data that is relevant to HSPC biology.

The human GFI136N variant induces epigenetic changes at the Hoxa9 locus and accelerates K-RAS driven myeloproliferative disorder in mice.

The coding single nucleotide polymorphism GFI136N in the human gene growth factor independence 1 (GFI1) is present in 3%-7% of whites and increases the risk for acute myeloid leukemia (AML) by 60%. We show here that GFI136N, in contrast to GFI136S, lacks the ability to bind to the Gfi1 target gene that encodes the leukemia-associated transcription factor Hoxa9 and fails to initiate histone modifications that regulate HoxA9 expression. Consistent with this, AML patients heterozygous for the GFI136N variant show increased HOXA9 expression compared with normal controls. Using ChipSeq, we demonstrate that GFI136N specific epigenetic changes are also present in other genes involved in the development of AML. Moreover, granulomonocytic progenitors, a bone marrow subset from which AML can arise in humans and mice, show a proliferative expansion in the presence of the GFI136N variant. In addition, granulomonocytic progenitors carrying the GFI136N variant allele have altered gene expression patterns and differ in their ability to grow after transplantation. Finally, GFI136N can accelerate a K-RAS driven fatal myeloproliferative disease in mice. Our data suggest that the presence of a GFI136N variant allele induces a preleukemic state in myeloid precursors by deregulating the expression of Hoxa9 and other AML-related genes.

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.

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.

The ATR inhibitor ceralasertib potentiates cancer checkpoint immunotherapy by regulating the tumor microenvironment.

The Ataxia telangiectasia and Rad3-related (ATR) inhibitor ceralasertib in combination with the PD-L1 antibody durvalumab demonstrated encouraging clinical benefit in melanoma and lung cancer patients who progressed on immunotherapy. Here we show that modelling of intermittent ceralasertib treatment in mouse tumor models reveals CD8+ T-cell dependent antitumor activity, which is separate from the effects on tumor cells. Ceralasertib suppresses proliferating CD8+ T-cells on treatment which is rapidly reversed off-treatment. Ceralasertib causes up-regulation of type I interferon (IFNI) pathway in cancer patients and in tumor-bearing mice. IFNI is experimentally found to be a major mediator of antitumor activity of ceralasertib in combination with PD-L1 antibody. Improvement of T-cell function after ceralasertib treatment is linked to changes in myeloid cells in the tumor microenvironment. IFNI also promotes anti-proliferative effects of ceralasertib on tumor cells. Here, we report that broad immunomodulatory changes following intermittent ATR inhibition underpins the clinical therapeutic benefit and indicates its wider impact on antitumor immunity.

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.

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.

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.