CellRank for directed single-cell fate mapping

Single-cell dissection of human hematopoietic reconstitution after allogeneic hematopoietic stem cell transplantation

Hematopoietic stem cell transplantation is an effective regenerative therapy for many malignant, inherited, or autoimmune diseases. However, our understanding of reconstituted hematopoiesis in transplant patients remains limited. Here, we uncover the reconstitution dynamics of human allogeneic hematopoietic stem and progenitor cells (HSPCs) at single-cell resolution after transplantation. Transplanted HSPCs underwent rapid and measurable changes during the first 30 days after transplantation, characterized by a strong proliferative response on the first day. Transcriptomic analysis of HSPCs enabled us to observe that immunoregulatory neutrophil progenitors expressing high levels of the S100A gene family were enriched in granulocyte colony-stimulating factor-mobilized peripheral blood stem cells. Transplant recipients who developed acute graft-versus-host disease (aGVHD) infused fewer S100A^high immunoregulatory neutrophil progenitors, immunophenotyped as Lin^-CD34⁺CD66b⁺CD177⁺, than those who did not develop aGVHD. Therefore, our study provides insights into the regenerative process of transplanted HSPCs in human patients and identifies a potential criterion for identifying patients at high risk for developing aGVHD early after transplant.

The pre-existing T cell landscape determines the response to bispecific T cell engagers in multiple myeloma patients

Bispecific T cell engagers (TCEs) have shown promise in the treatment of various cancers, but the immunological mechanism and molecular determinants of primary and acquired resistance to TCEs remain poorly understood. Here, we identify conserved behaviors of bone marrow-residing T cells in multiple myeloma patients undergoing BCMAxCD3 TCE therapy. We show that the immune repertoire reacts to TCE therapy with cell state-dependent clonal expansion and find evidence supporting the coupling of tumor recognition via major histocompatibility complex class I (MHC class I), exhaustion, and clinical response. We find the abundance of exhausted-like CD8⁺ T cell clones to be associated with clinical response failure, and we describe loss of target epitope and MHC class I as tumor-intrinsic adaptations to TCEs. These findings advance our understanding of the in vivo mechanism of TCE treatment in humans and provide the rationale for predictive immune-monitoring and conditioning of the immune repertoire to guide future immunotherapy in hematological malignancies.

The molecular consequences of androgen activity in the human breast

Estrogen and progesterone have been extensively studied in the mammary gland, but the molecular effects of androgen remain largely unexplored. Transgender men are recorded as female at birth but identify as male and may undergo gender-affirming androgen therapy to align their physical characteristics and gender identity. Here we perform single-cell-resolution transcriptome, chromatin, and spatial profiling of breast tissues from transgender men following androgen therapy. We find canonical androgen receptor gene targets are upregulated in cells expressing the androgen receptor and that paracrine signaling likely drives sex-relevant androgenic effects in other cell types. We also observe involution of the epithelium and a spatial reconfiguration of immune, fibroblast, and vascular cells, and identify a gene regulatory network associated with androgen-induced fat loss. This work elucidates the molecular consequences of androgen activity in the human breast at single-cell resolution.

Adipose tissue at single-cell resolution

Adipose tissue exhibits remarkable plasticity with capacity to change in size and cellular composition under physiological and pathophysiological conditions. The emergence of single-cell transcriptomics has rapidly transformed our understanding of the diverse array of cell types and cell states residing in adipose tissues and has provided insight into how transcriptional changes in individual cell types contribute to tissue plasticity. Here, we present a comprehensive overview of the cellular atlas of adipose tissues focusing on the biological insight gained from single-cell and single-nuclei transcriptomics of murine and human adipose tissues. We also offer our perspective on the exciting opportunities for mapping cellular transitions and crosstalk, which have been made possible by single-cell technologies.

Active eosinophils regulate host defence and immune responses in colitis

In the past decade, single-cell transcriptomics has helped to uncover new cell types and states and led to the construction of a cellular compendium of health and disease. Despite this progress, some difficult-to-sequence cells remain absent from tissue atlases. Eosinophils-elusive granulocytes that are implicated in a plethora of human pathologies1-5-are among these uncharted cell types. The heterogeneity of eosinophils and the gene programs that underpin their pleiotropic functions remain poorly understood. Here we provide a comprehensive single-cell transcriptomic profiling of mouse eosinophils. We identify an active and a basal population of intestinal eosinophils, which differ in their transcriptome, surface proteome and spatial localization. By means of a genome-wide CRISPR inhibition screen and functional assays, we reveal a mechanism by which interleukin-33 (IL-33) and interferon-γ (IFNγ) induce the accumulation of active eosinophils in the inflamed colon. Active eosinophils are endowed with bactericidal and T cell regulatory activity, and express the co-stimulatory molecules CD80 and PD-L1. Notably, active eosinophils are enriched in the lamina propria of a small cohort of patients with inflammatory bowel disease, and are closely associated with CD4+ T cells. Our findings provide insights into the biology of eosinophils and highlight the crucial contribution of this cell type to intestinal homeostasis, immune regulation and host defence. Furthermore, we lay a framework for the characterization of eosinophils in human gastrointestinal diseases.

The early neutrophil-committed progenitors aberrantly differentiate into immunoregulatory monocytes during emergency myelopoiesis

Inflammatory stimuli cause a state of emergency myelopoiesis leading to neutrophil-like monocyte expansion. However, their function, the committed precursors, or growth factors remain elusive. In this study we find that Ym1⁺Ly6C^hi monocytes, an immunoregulatory entity of neutrophil-like monocytes, arise from progenitors of neutrophil 1 (proNeu1). Granulocyte-colony stimulating factor (G-CSF) favors the production of neutrophil-like monocytes through previously unknown CD81⁺CX3CR1^lo monocyte precursors. GFI1 promotes the differentiation of proNeu2 from proNeu1 at the cost of producing neutrophil-like monocytes. The human counterpart of neutrophil-like monocytes that also expands in response to G-CSF is found in CD14⁺CD16^- monocyte fraction. The human neutrophil-like monocytes are discriminated from CD14⁺CD16^- classical monocytes by CXCR1 expression and the capacity to suppress T cell proliferation. Collectively, our findings suggest that the aberrant expansion of neutrophil-like monocytes under inflammatory conditions is a process conserved between mouse and human, which may be beneficial for the resolution of inflammation.

The pseudokinase Trib1 regulates the transition of exhausted T cells to a KLR + CD8 + effector state and its deletion improves checkpoint blockade

T cell exhaustion (T EX ) impairs the ability of T cells to clear chronic infection or cancer. While exhausted T cells are hypofunctional, some exhausted T cells retain effector gene signatures, a feature that is associated with expression of KLRs (killer lectin-like receptors). Although KLR + T cells may improve control of chronic antigen, the signaling molecules regulating this population are poorly understood. Using scRNA-seq, flow cytometry, RNA velocity, and scTCR-seq, we demonstrate that deleting the pseudokinase Trib1 shifts T EX towards CX3CR1 + intermediates (T INT ) with robust enrichment of KLR + CD8 + T cells (T KLR ) via clonal T cell expansion. These changes are associated with globally increased KLR gene expression throughout the exhaustion program. Further, Trib1 loss augments anti-PD-L1 blockade to improve viral clearance by expanding the T KLR population. Together, these data identify Trib1 as an important regulator of T cell exhaustion whose targeting enhances the KLR + effector state and improves the response to checkpoint inhibitor therapy.

Adult hippocampal neurogenesis in Alzheimer's disease: A roadmap to clinical relevance

Adult hippocampal neurogenesis (AHN) drops sharply during early stages of Alzheimer's disease (AD), via unknown mechanisms, and correlates with cognitive status in AD patients. Understanding AHN regulation in AD could provide a framework for innovative pharmacological interventions. We here combine molecular, behavioral, and clinical data and critically discuss the multicellular complexity of the AHN niche in relation to AD pathophysiology. We further present a roadmap toward a better understanding of the role of AHN in AD by probing the promises and caveats of the latest technological advancements in the field and addressing the conceptual and methodological challenges ahead.

Tbx2 and Tbx3 regulate cell fate progression of the otic vesicle for inner ear development

The morphogenesis of the otic vesicle (OV) to form inner ear organs serves as an excellent model system to understand cell fate acquisition on a single cell level. Tbx2 and Tbx3 (Tbx2/3) encode closely related T-box transcription factors that are expressed widely in the mammalian OV. Inactivation of both genes in the OV (Tbx2/3cKO) results in failed morphogenesis into inner ear organs. To understand the basis of these defects, single cell RNA-sequencing (scRNA-seq) was performed on the OV lineage, in controls versus Tbx2/3cKO embryos. We identified a multipotent population termed otic progenitors in controls that are marked by expression of the known otic placode markers Eya1, Sox2, and Sox3 as well as new markers Fgf18, Cxcl12, and Pou3f3. The otic progenitor population was increased three-fold in Tbx2/3cKO embryos, concomitant with dysregulation of genes in these cells as well as reduced progression to more differentiated states of prosensory and nonsensory cells. An ectopic neural population of cells was detected in the posterior OV of Tbx2/3cKO embryos but had reduced maturation to delaminated neural cells. As all three cell fates were affected in Tbx2/3cKO embryos, we suggest that Tbx2/3 promotes progression of multipotent otic progenitors to more differentiated cell types in the OV.

Reconstruction and deconstruction of human somitogenesis in vitro

The vertebrate body displays a segmental organization that is most conspicuous in the periodic organization of the vertebral column and peripheral nerves. This metameric organization is first implemented when somites, which contain the precursors of skeletal muscles and vertebrae, are rhythmically generated from the presomitic mesoderm. Somites then become subdivided into anterior and posterior compartments that are essential for vertebral formation and segmental patterning of the peripheral nervous system^1-4. How this key somitic subdivision is established remains poorly understood. Here we introduce three-dimensional culture systems of human pluripotent stem cells called somitoids and segmentoids, which recapitulate the formation of somite-like structures with anteroposterior identity. We identify a key function of the segmentation clock in converting temporal rhythmicity into the spatial regularity of anterior and posterior somitic compartments. We show that an initial 'salt and pepper' expression of the segmentation gene MESP2 in the newly formed segment is transformed into compartments of anterior and posterior identity through an active cell-sorting mechanism. Our research demonstrates that the major patterning modules that are involved in somitogenesis, including the clock and wavefront, anteroposterior polarity patterning and somite epithelialization, can be dissociated and operate independently in our in vitro systems. Together, we define a framework for the symmetry-breaking process that initiates somite polarity patterning. Our work provides a platform for decoding general principles of somitogenesis and advancing knowledge of human development.

Single-cell multi-omics and lineage tracing to dissect cell fate decision-making

The concept of cell fate relates to the future identity of a cell, and its daughters, which is obtained via cell differentiation and division. Understanding, predicting, and manipulating cell fate has been a long-sought goal of developmental and regenerative biology. Recent insights obtained from single-cell genomic and integrative lineage-tracing approaches have further aided to identify molecular features predictive of cell fate. In this perspective, we discuss these approaches with a focus on theoretical concepts and future directions of the field to dissect molecular mechanisms underlying cell fate.

Decoding anterior-posterior axis emergence among mouse, monkey, and human embryos

Anterior-posterior axis formation regulated by the distal visceral endoderm (DVE) and anterior visceral endoderm (AVE) is essential for peri-implantation embryogenesis. However, the principles of the origin and specialization of DVE and AVE remain elusive. Here, with single-cell transcriptome analysis and pseudotime prediction, we show that DVE and AVE independently originate from the specialized primary endoderm in mouse blastocysts. Along distinct developmental paths, these two lineages, respectively, undergo four representative states with stage-specific transcriptional patterns around implantation. Further comparative analysis shows that AVE, but not DVE, is detected in human and non-human primate embryos, defining differences in polarity formation across species. Moreover, stem cell-assembled human blastoids lack DVE or AVE precursors, implying that additional induction of stem cells with DVE/AVE potential could promote the current embryo-like models and their post-implantation growth. Our work provides insight into understanding of embryonic polarity formation and early mammalian development.

Dynamic visualization of high-dimensional data

Dimensionality reduction (DR) is commonly used to project high-dimensional data into lower dimensions for visualization, which could then generate new insights and hypotheses. However, DR algorithms introduce distortions in the visualization and cannot faithfully represent all relations in the data. Thus, there is a need for methods to assess the reliability of DR visualizations. Here we present DynamicViz, a framework for generating dynamic visualizations that capture the sensitivity of DR visualizations to perturbations in the data resulting from bootstrap sampling. DynamicViz can be applied to all commonly used DR methods. We show the utility of dynamic visualizations in diagnosing common interpretative pitfalls of static visualizations and extending existing single-cell analyses. We introduce the variance score to quantify the dynamic variability of observations in these visualizations. The variance score characterizes natural variability in the data and can be used to optimize DR algorithm implementations.

scFates: a scalable python package for advanced pseudotime and bifurcation analysis from single-cell data

SUMMARY

scFates provides an extensive toolset for the analysis of dynamic trajectories comprising tree learning, feature association testing, branch differential expression and with a focus on cell biasing and fate splits at the level of bifurcations. It is meant to be fully integrated into the scanpy ecosystem for seamless analysis of trajectories from single-cell data of various modalities (e.g. RNA and ATAC).

AVAILABILITY AND IMPLEMENTATION

scFates is released as open-source software under the BSD 3-Clause 'New' License and is available from the Python Package Index at https://pypi.org/project/scFates/. The source code is available on GitHub at https://github.com/LouisFaure/scFates/. Code reproduction and tutorials on published datasets are available on GitHub at https://github.com/LouisFaure/scFates_notebooks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A topographic atlas defines developmental origins of cell heterogeneity in the human embryonic lung

The lung contains numerous specialized cell types with distinct roles in tissue function and integrity. To clarify the origins and mechanisms generating cell heterogeneity, we created a comprehensive topographic atlas of early human lung development. Here we report 83 cell states and several spatially resolved developmental trajectories and predict cell interactions within defined tissue niches. We integrated single-cell RNA sequencing and spatially resolved transcriptomics into a web-based, open platform for interactive exploration. We show distinct gene expression programmes, accompanying sequential events of cell differentiation and maturation of the secretory and neuroendocrine cell types in proximal epithelium. We define the origin of airway fibroblasts associated with airway smooth muscle in bronchovascular bundles and describe a trajectory of Schwann cell progenitors to intrinsic parasympathetic neurons controlling bronchoconstriction. Our atlas provides a rich resource for further research and a reference for defining deviations from homeostatic and repair mechanisms leading to pulmonary diseases.

Single-cell and long-read sequencing to enhance modelling of splicing and cell-fate determination

Single-cell sequencing technologies have revolutionised the life sciences and biomedical research. Single-cell sequencing provides high-resolution data on cell heterogeneity, allowing high-fidelity cell type identification, and lineage tracking. Computational algorithms and mathematical models have been developed to make sense of the data, compensate for errors and simulate the biological processes, which has led to breakthroughs in our understanding of cell differentiation, cell-fate determination and tissue cell composition. The development of long-read (a.k.a. third-generation) sequencing technologies has produced powerful tools for investigating alternative splicing, isoform expression (at the RNA level), genome assembly and the detection of complex structural variants (at the DNA level). In this review, we provide an overview of the recent advancements in single-cell and long-read sequencing technologies, with a particular focus on the computational algorithms that help in correcting, analysing, and interpreting the resulting data. Additionally, we review some mathematical models that use single-cell and long-read sequencing data to study cell-fate determination and alternative splicing, respectively. Moreover, we highlight the emerging opportunities in modelling cell-fate determination that result from the combination of single-cell and long-read sequencing technologies.

Transcriptional networks of transient cell states during human prefrontal cortex development

The human brain is divided into various anatomical regions that control and coordinate unique functions. The prefrontal cortex (PFC) is a large brain region that comprises a range of neuronal and non-neuronal cell types, sharing extensive interconnections with subcortical areas, and plays a critical role in cognition and memory. A timely appearance of distinct cell types through embryonic development is crucial for an anatomically perfect and functional brain. Direct tracing of cell fate development in the human brain is not possible, but single-cell transcriptome sequencing (scRNA-seq) datasets provide the opportunity to dissect cellular heterogeneity and its molecular regulators. Here, using scRNA-seq data of human PFC from fetal stages, we elucidate distinct transient cell states during PFC development and their underlying gene regulatory circuitry. We further identified that distinct intermediate cell states consist of specific gene regulatory modules essential to reach terminal fate using discrete developmental paths. Moreover, using in silico gene knock-out and over-expression analysis, we validated crucial gene regulatory components during the lineage specification of oligodendrocyte progenitor cells. Our study illustrates unique intermediate states and specific gene interaction networks that warrant further investigation for their functional contribution to typical brain development and discusses how this knowledge can be harvested for therapeutic intervention in challenging neurodevelopmental disorders.

Simulation-based inference of differentiation trajectories from RNA velocity fields

We report Cytopath, a method for trajectory inference that takes advantage of transcriptional activity information from the RNA velocity of single cells to perform trajectory inference. Cytopath performs this task by defining a Markov chain model, simulating an ensemble of possible differentiation trajectories, and constructing a consensus trajectory. We show that Cytopath can recapitulate the topological and molecular characteristics of the differentiation process under study. In our analysis, we include differentiation trajectories with varying bifurcated, circular, convergent, and mixed topologies studied in single-snapshot as well as time-series single-cell RNA sequencing experiments. We demonstrate the capability to reconstruct differentiation trajectories, assess the association of RNA velocity-based pseudotime with actually elapsed process time, and identify drawbacks in current state-of-the art trajectory inference approaches.

Transcriptional advantage influence odorant receptor gene choice

Odorant receptors (ORs) obey mutual exclusivity and monoallelic mode of expression. Efforts are ongoing to decipher the molecular mechanism that drives the ‘one-neuron-one-receptor’ rule of olfaction. Recently, single-cell profiling of olfactory sensory neurons (OSNs) revealed the expression of multiple ORs in the immature neurons, suggesting that the OR gene choice mechanism is much more complex than previously described by the silence-all-and-activate-one model. These results also led to the genesis of two possible mechanistic models i.e. winner-takes-all and stochastic selection. We developed Reverse Cell Tracking (RCT), a novel computational framework that facilitates OR-guided cellular backtracking by leveraging Uniform Manifold Approximation and Projection embeddings from RNA Velocity Workflow. RCT-based trajectory backtracking, coupled with statistical analysis, revealed the OR gene choice bias for the transcriptionally advanced (highest expressed) OR during neuronal differentiation. Interestingly, the observed selection bias was uniform for all ORs across different spatial zones or their relative expression within the olfactory organ. We validated these findings on independent datasets and further confirmed that the OR gene selection may be regulated by Upf3b. Lastly, our RNA dynamics-based tracking of the differentiation cascade revealed a transition cell state that harbors mixed molecular identities of immature and mature OSNs, and their relative abundance is regulated by Upf3b.

Single-cell microglial transcriptomics during demyelination defines a microglial state required for lytic carcass clearance

BACKGROUND

Microglia regulate the response to injury and disease in the brain and spinal cord. In white matter diseases microglia may cause demyelination. However, how microglia respond and regulate demyelination is not fully understood.

METHODS

To understand how microglia respond during demyelination, we fed mice cuprizone-a potent demyelinating agent-and assessed the dynamics of genetically fate-mapped microglia. We then used single-cell RNA sequencing to identify and track the microglial subpopulations that arise during demyelination. To understand how microglia contribute to the clearance of dead oligodendrocytes, we ablated microglia starting at the peak of cuprizone-induced cell death and used the viability dye acridine orange to monitor apoptotic and lytic cell morphologies after microglial ablation. Lastly, we treated serum-free primary microglial cultures to model distinct aspects of cuprizone-induced demyelination and assessed the response.

RESULTS

The cuprizone diet generated a robust microglial response by week 4 of the diet. Single-cell RNA sequencing at this time point revealed the presence of several cuprizone-associated microglia (CAM) clusters. These clusters expressed a transcriptomic signature indicative of cytokine regulation and reactive oxygen species production with altered lysosomal and metabolic changes consistent with ongoing phagocytosis. Using acridine orange to monitor apoptotic and lytic cell death after microglial ablation, we found that microglia preferentially phagocytose lytic carcasses. In culture, microglia exposed to lytic carcasses partially recapitulated the CAM state, suggesting that phagocytosis contributes to this distinct microglial state during cuprizone demyelination.

CONCLUSIONS

Microglia serve multiple roles during demyelination, yet their transcriptomic state resembles other neurodegenerative conditions. The phagocytosis of cellular debris is likely a universal cause for a common neurodegenerative microglial state.

Single-cell transcriptome analyses reveal microglia types associated with proliferative retinopathy

Pathological angiogenesis is a major cause of irreversible blindness in individuals of all age groups with proliferative retinopathy (PR). Mononuclear phagocytes (MPs) within neovascular areas contribute to aberrant retinal angiogenesis. Due to their cellular heterogeneity, defining the roles of MP subsets in PR onset and progression has been challenging. Here, we aimed to investigate the heterogeneity of microglia associated with neovascularization and to characterize the transcriptional profiles and metabolic pathways of proangiogenic microglia in a mouse model of oxygen-induced PR (OIR). Using transcriptional single-cell sorting, we comprehensively mapped all microglia populations in retinas of room air (RA) and OIR mice. We have unveiled several unique types of PR-associated microglia (PRAM) and identified markers, signaling pathways, and regulons associated with these cells. Among these microglia subpopulations, we found a highly proliferative microglia subset with high self-renewal capacity and a hypermetabolic microglia subset that expresses high levels of activating microglia markers, glycolytic enzymes, and proangiogenic Igf1. IHC staining shows that these PRAM were spatially located within or around neovascular tufts. These unique types of microglia have the potential to promote retinal angiogenesis, which may have important implications for future treatment of PR and other pathological ocular angiogenesis-related diseases.

Fibroblast inflammatory priming determines regenerative versus fibrotic skin repair in reindeer

Adult mammalian skin wounds heal by forming fibrotic scars. We report that full-thickness injuries of reindeer antler skin (velvet) regenerate, whereas back skin forms fibrotic scar. Single-cell multi-omics reveal that uninjured velvet fibroblasts resemble human fetal fibroblasts, whereas back skin fibroblasts express inflammatory mediators mimicking pro-fibrotic adult human and rodent fibroblasts. Consequently, injury elicits site-specific immune responses: back skin fibroblasts amplify myeloid infiltration and maturation during repair, whereas velvet fibroblasts adopt an immunosuppressive phenotype that restricts leukocyte recruitment and hastens immune resolution. Ectopic transplantation of velvet to scar-forming back skin is initially regenerative, but progressively transitions to a fibrotic phenotype akin to the scarless fetal-to-scar-forming transition reported in humans. Skin regeneration is diminished by intensifying, or enhanced by neutralizing, these pathologic fibroblast-immune interactions. Reindeer represent a powerful comparative model for interrogating divergent wound healing outcomes, and our results nominate decoupling of fibroblast-immune interactions as a promising approach to mitigate scar.

Comprehensive single-cell transcriptional profiling defines shared and unique epithelial injury responses during kidney fibrosis

The underlying cellular events driving kidney fibrogenesis and metabolic dysfunction are incompletely understood. Here, we employed single-cell combinatorial indexing RNA sequencing to analyze 24 mouse kidneys from two fibrosis models. We profiled 309,666 cells in one experiment, representing 50 cell types/states encompassing epithelial, endothelial, immune, and stromal populations. Single-cell analysis identified diverse injury states of the proximal tubule, including two distinct early-phase populations with dysregulated lipid and amino acid metabolism, respectively. Lipid metabolism was defective in the chronic phase but was transiently activated in the very early stages of ischemia-induced injury, where we discovered increased lipid deposition and increased fatty acid β-oxidation. Perilipin 2 was identified as a surface marker of intracellular lipid droplets, and its knockdown in vitro disrupted cell energy state maintenance during lipid accumulation. Surveying epithelial cells across nephron segments identified shared and unique injury responses. Stromal cells exhibited high heterogeneity and contributed to fibrogenesis by epithelial-stromal crosstalk.

High precision, high throughput generation of droplets containing single cells

The Poisson limit is a major problem for the isolation of single cells in different single-cell technologies and applications. In droplet-based single-cell assays, a scheme that is increasingly popular, the intrinsic randomness during single-cell encapsulation in droplets requires most of the created droplets to be empty, which has a profound impact on the efficiency and throughput of such techniques, and on the predictability of the combinatory droplet assays. Here we present a simple passive microfluidic system overcoming this limitation with unprecedented efficacy, allowing the generation of single-cell droplets for a wide range of operating conditions, with extremely high throughput (more than 22 000 single-cell loaded droplets per minute) and with an extremely low fault ratio (doublets or empty droplets), applicable to any cells and deformable particles. This versatile technique will shift the paradigm of single-cell encapsulation and will impact single-cell sequencing, rare cell isolation, multicellular/bead studies in immunology or cancer biology, etc.

Dental niche cells directly contribute to tooth reconstitution and morphogenesis

Mammalian teeth develop from the inductive epithelial-mesenchymal interaction, an important mechanism shared by many organs. The cellular basis for such interaction remains elusive. Here, we generate a dual-fluorescence model to track and analyze dental cells from embryonic to postnatal stages, in which Pitx2⁺ epithelium and Msx1⁺ mesenchyme are sufficient for tooth reconstitution. Single-cell RNA sequencing and spatial mapping further revealed critical cellular dynamics during molar development, where tooth germs are organized by Msx1⁺Sdc1⁺ dental papilla and surrounding dental niche. Surprisingly, niche cells are more efficient in tooth reconstitution and can directly regenerate papilla cells through interaction with dental epithelium. Finally, from the dental niche, we identify a group of previously unappreciated migratory Msx1⁺ Sox9⁺ cells as the potential cell origin for dental papilla. Our results indicate that the dental niche cells directly contribute to tooth organogenesis and provide critical insights into the essential cell composition for tooth engineering.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

Biologic insights from single-cell studies of psoriasis and psoriatic arthritis

INTRODUCTION

Psoriasis (PSO) and psoriatic arthritis (PSA) represent a large burden of global inflammatory disease, but sustained treatment response and early diagnosis remain challenging. Both conditions arise from complex immune cell dysregulation. Single-cell techniques, including single-cell RNA sequencing (scRNA-seq), have revolutionized our understanding of pathogenesis by illuminating heterogeneous cell populations and their interactions.

AREAS COVERED

We discuss the transcriptional profiles and cellular interactions unique to PSO/PSA affecting T cells, myeloid cells, keratinocytes, innate lymphoid cells, and stromal cells. We also review advances, limitations, and future challenges associated with single-cell studies.

EXPERT OPINION

Following analyses of 22 single-cell studies, several themes emerged. A small subpopulation of cells can have a large impact on disease pathogenesis. Multiple cell types identified via scRNA-seq play supporting roles in PSO pathogenesis, contrary to the traditional paradigm focusing on IL-23/IL-17 signaling among dendritic cells and T cells. Immune cell states are dynamic, with psoriatic subpopulations aberrantly re-activating and differentiating into inflammatory phenotypes depending on surrounding signaling cues. Comparison of circulating immune cells with resident skin/joint cells has uncovered specific T cell clonotypes associated with the disease. Finally, machine learning models demonstrate great promise in identifying biomarkers to diagnose clinically ambiguous rashes and PSA at earlier stages.

Liver tumour immune microenvironment subtypes and neutrophil heterogeneity

The heterogeneity of the tumour immune microenvironment (TIME), organized by various immune and stromal cells, is a major contributing factor of tumour metastasis, relapse and drug resistance^1-3, but how different TIME subtypes are connected to the clinical relevance in liver cancer remains unclear. Here we performed single-cell RNA-sequencing (scRNA-seq) analysis of 189 samples collected from 124 patients and 8 mice with liver cancer. With more than 1 million cells analysed, we stratified patients into five TIME subtypes, including immune activation, immune suppression mediated by myeloid or stromal cells, immune exclusion and immune residence phenotypes. Different TIME subtypes were spatially organized and associated with chemokine networks and genomic features. Notably, tumour-associated neutrophil (TAN) populations enriched in the myeloid-cell-enriched subtype were associated with an unfavourable prognosis. Through in vitro induction of TANs and ex vivo analyses of patient TANs, we showed that CCL4⁺ TANs can recruit macrophages and that PD-L1⁺ TANs can suppress T cell cytotoxicity. Furthermore, scRNA-seq analysis of mouse neutrophil subsets revealed that they are largely conserved with those of humans. In vivo neutrophil depletion in mouse models attenuated tumour progression, confirming the pro-tumour phenotypes of TANs. With this detailed cellular heterogeneity landscape of liver cancer, our study illustrates diverse TIME subtypes, highlights immunosuppressive functions of TANs and sheds light on potential immunotherapies targeting TANs.

The fate of early perichondrial cells in developing bones

In endochondral bone development, bone-forming osteoblasts and bone marrow stromal cells have dual origins in the fetal cartilage and its surrounding perichondrium. However, how early perichondrial cells distinctively contribute to developing bones remain unidentified. Here we show using in vivo cell-lineage analyses that Dlx5+ fetal perichondrial cells marked by Dlx5-creER do not generate cartilage but sustainably contribute to cortical bone and marrow stromal compartments in a manner complementary to fetal chondrocyte derivatives under the regulation of Hedgehog signaling. Postnatally, Dlx5+ fetal perichondrial cell derivatives preferentially populate the diaphyseal marrow stroma with a dormant adipocyte-biased state and are refractory to parathyroid hormone-induced bone anabolism. Therefore, early perichondrial cells of the fetal cartilage are destined to become an adipogenic subset of stromal cells in postnatal diaphyseal bone marrow, supporting the theory that the adult bone marrow stromal compartments are developmentally prescribed within the two distinct cells-of-origins of the fetal bone anlage.

Single-cell transcriptome analysis reveals cellular heterogeneity in mouse intra- and extra articular ligaments

Ligaments are collagenous connective tissues that connect bones. Injury of knee ligaments, namely anterior cruciate ligament (ACL) and medial collateral ligament (MCL), is common in athletes. Both ligaments have important functions, but distinct regeneration capacities. The capacity for recovery after injury also diminishes with age. However, cellular heterogeneity in the ligaments remains unclear. Here, we profiled the transcriptional signatures of ACL and MCL cells in mice using single-cell RNA sequencing. These ligaments comprise three fibroblast types expressing Col22a1, Col12a1, or Col14a1, but have distinct localizations in the tissue. We found substantial heterogeneity in Col12a1- and Col14a1-positive cells between ACL and MCL. Gene Ontology analysis revealed that angiogenesis- and collagen regulation-related genes were specifically enriched in MCL cells. Furthermore, we identified age-related changes in cell composition and gene expression in the ligaments. This study delineates cellular heterogeneity in ligaments, serving as a foundation for identifying potential therapeutic targets for ligament injuries.

Cell heterogeneity in the mouse anterior cruciate ligament (ACL) and medial collateral ligament (MCL) is demonstrated using single-cell analysis with three types of fibroblasts identified, expressing Col14a1, Col12a1, or Col22a1.

UniTVelo: temporally unified RNA velocity reinforces single-cell trajectory inference

The recent breakthrough of single-cell RNA velocity methods brings attractive promises to reveal directed trajectory on cell differentiation, states transition and response to perturbations. However, the existing RNA velocity methods are often found to return erroneous results, partly due to model violation or lack of temporal regularization. Here, we present UniTVelo, a statistical framework of RNA velocity that models the dynamics of spliced and unspliced RNAs via flexible transcription activities. Uniquely, it also supports the inference of a unified latent time across the transcriptome. With ten datasets, we demonstrate that UniTVelo returns the expected trajectory in different biological systems, including hematopoietic differentiation and those even with weak kinetics or complex branches.

Deciphering molecular and cellular ex vivo responses to bispecific antibodies PD1-TIM3 and PD1-LAG3 in human tumors

BACKGROUND

Next-generation cancer immunotherapies are designed to broaden the therapeutic repertoire by targeting new immune checkpoints including lymphocyte-activation gene 3 (LAG-3) and T cell immunoglobulin and mucin-domain containing-3 (TIM-3). Yet, the molecular and cellular mechanisms by which either receptor functions to mediate its inhibitory effects are still poorly understood. Similarly, little is known on the differential effects of dual, compared with single, checkpoint inhibition.

METHODS

We here performed in-depth characterization, including multicolor flow cytometry, single cell RNA sequencing and multiplex supernatant analysis, using tumor single cell suspensions from patients with cancer treated ex vivo with novel bispecific antibodies targeting programmed cell death protein 1 (PD-1) and TIM-3 (PD1-TIM3), PD-1 and LAG-3 (PD1-LAG3), or with anti-PD-1.

RESULTS

We identified patient samples which were responsive to PD1-TIM3, PD1-LAG3 or anti-PD-1 using an in vitro approach, validated by the analysis of 659 soluble proteins and enrichment for an anti-PD-1 responder signature. We found increased abundance of an activated (HLA-DR⁺CD25⁺GranzymeB⁺) CD8⁺ T cell subset and of proliferating CD8⁺ T cells, in response to bispecific antibody or anti-PD-1 treatment. Bispecific antibodies, but not anti-PD-1, significantly increased the abundance of a proliferating natural killer cell subset, which exhibited enrichment for a tissue-residency signature. Key phenotypic and transcriptional changes occurred in a PD-1⁺CXCL13⁺CD4⁺ T cell subset, in response to all treatments, including increased interleukin-17 secretion and signaling toward plasma cells. Interestingly, LAG-3 protein upregulation was detected as a unique pharmacodynamic effect mediated by PD1-LAG3, but not by PD1-TIM3 or anti-PD-1.

CONCLUSIONS

Our in vitro system reliably assessed responses to bispecific antibodies co-targeting PD-1 together with LAG-3 or TIM-3 using patients' tumor infiltrating immune cells and revealed transcriptional and phenotypic imprinting by bispecific antibody formats currently tested in early clinical trials.

CXCL12 defines lung endothelial heterogeneity and promotes distal vascular growth

There is a growing amount of data uncovering the cellular diversity of the pulmonary circulation and mechanisms governing vascular repair after injury. However, the molecular and cellular mechanisms contributing to the morphogenesis and growth of the pulmonary vasculature during embryonic development are less clear. Importantly, deficits in vascular development lead to significant pediatric lung diseases, indicating a need to uncover fetal programs promoting vascular growth. To address this, we used a transgenic mouse reporter for expression of Cxcl12, an arterial endothelial hallmark gene, and performed single-cell RNA sequencing on isolated Cxcl12-DsRed+ endothelium to assess cellular heterogeneity within pulmonary endothelium. Combining cell annotation with gene ontology and histological analysis allowed us to segregate the developing artery endothelium into functionally and spatially distinct subpopulations. Expression of Cxcl12 is highest in the distal arterial endothelial subpopulation, a compartment enriched in genes for vascular development. Accordingly, disruption of CXCL12 signaling led to, not only abnormal branching, but also distal vascular hypoplasia. These data provide evidence for arterial endothelial functional heterogeneity and reveal conserved signaling mechanisms essential for pulmonary vascular development.

The tempo and mode of gene regulatory programs during bacterial infection

Innate immune recognition of bacterial pathogens is a key determinant of the ensuing systemic response, and host or pathogen heterogeneity in this early interaction can impact the course of infection. To gain insight into host response heterogeneity, we investigate macrophage inflammatory dynamics using primary human macrophages infected with Group B Streptococcus. Transcriptomic analysis reveals discrete cellular states within responding macrophages, one of which consists of four sub-states, reflecting inflammatory activation. Infection with six additional bacterial species-Staphylococcus aureus, Listeria monocytogenes, Enterococcus faecalis, Yersinia pseudotuberculosis, Shigella flexneri, and Salmonella enterica-recapitulates these states, though at different frequencies. We show that modulating the duration of infection and the presence of a toxin impacts inflammatory trajectory dynamics. We provide evidence for this trajectory in infected macrophages in an in vivo model of Staphylococcus aureus infection. Our cell-state analysis defines a framework for understanding inflammatory activation dynamics in response to bacterial infection.

Functional Heterogeneity of Bone Marrow Mesenchymal Stem Cell Subpopulations in Physiology and Pathology

Bone marrow mesenchymal stem cells (BMSCs) are multi-potent cell populations and are capable of maintaining bone and body homeostasis. The stemness and potential therapeutic effect of BMSCs have been explored extensively in recent years. However, diverse cell surface antigens and complex gene expression of BMSCs have indicated that BMSCs represent heterogeneous populations, and the natural characteristics of BMSCs make it difficult to identify the specific subpopulations in pathological processes which are often obscured by bulk analysis of the total BMSCs. Meanwhile, the therapeutic effect of total BMSCs is often less effective partly due to their heterogeneity. Therefore, understanding the functional heterogeneity of the BMSC subpopulations under different physiological and pathological conditions could have major ramifications for global health. Here, we summarize the recent progress of functional heterogeneity of BMSC subpopulations in physiology and pathology. Targeting tissue-resident single BMSC subpopulation offers a potentially innovative therapeutic strategy and improves BMSC effectiveness in clinical application.

Archetype tasks link intratumoral heterogeneity to plasticity and cancer hallmarks in small cell lung cancer

Small cell lung cancer (SCLC) tumors comprise heterogeneous mixtures of cell states, categorized into neuroendocrine (NE) and non-neuroendocrine (non-NE) transcriptional subtypes. NE to non-NE state transitions, fueled by plasticity, likely underlie adaptability to treatment and dismal survival rates. Here, we apply an archetypal analysis to model plasticity by recasting SCLC phenotypic heterogeneity through multi-task evolutionary theory. Cell line and tumor transcriptomics data fit well in a five-dimensional convex polytope whose vertices optimize tasks reminiscent of pulmonary NE cells, the SCLC normal counterparts. These tasks, supported by knowledge and experimental data, include proliferation, slithering, metabolism, secretion, and injury repair, reflecting cancer hallmarks. SCLC subtypes, either at the population or single-cell level, can be positioned in archetypal space by bulk or single-cell transcriptomics, respectively, and characterized as task specialists or multi-task generalists by the distance from archetype vertex signatures. In the archetype space, modeling single-cell plasticity as a Markovian process along an underlying state manifold indicates that task trade-offs, in response to microenvironmental perturbations or treatment, may drive cell plasticity. Stifling phenotypic transitions and plasticity may provide new targets for much-needed translational advances in SCLC. A record of this paper's Transparent Peer Review process is included in the supplemental information.

NKG2A and HLA-E define an alternative immune checkpoint axis in bladder cancer

Programmed cell death protein 1 (PD-1)/programmed death-ligand 1 (PD-L1)-blockade immunotherapies have limited efficacy in the treatment of bladder cancer. Here, we show that NKG2A associates with improved survival and responsiveness to PD-L1 blockade immunotherapy in bladder tumors that have high abundance of CD8+ T cells. In bladder tumors, NKG2A is acquired on CD8+ T cells later than PD-1 as well as other well-established immune checkpoints. NKG2A+ PD-1+ CD8+ T cells diverge from classically defined exhausted T cells through their ability to react to human leukocyte antigen (HLA) class I-deficient tumors using T cell receptor (TCR)-independent innate-like mechanisms. HLA-ABC expression by bladder tumors is progressively diminished as disease progresses, framing the importance of targeting TCR-independent anti-tumor functions. Notably, NKG2A+ CD8+ T cells are inhibited when HLA-E is expressed by tumors and partly restored upon NKG2A blockade in an HLA-E-dependent manner. Overall, our study provides a framework for subsequent clinical trials combining NKG2A blockade with other T cell-targeted immunotherapies, where tumors express higher levels of HLA-E.

Integrating temporal single-cell gene expression modalities for trajectory inference and disease prediction

BACKGROUND

Current methods for analyzing single-cell datasets have relied primarily on static gene expression measurements to characterize the molecular state of individual cells. However, capturing temporal changes in cell state is crucial for the interpretation of dynamic phenotypes such as the cell cycle, development, or disease progression. RNA velocity infers the direction and speed of transcriptional changes in individual cells, yet it is unclear how these temporal gene expression modalities may be leveraged for predictive modeling of cellular dynamics.

RESULTS

Here, we present the first task-oriented benchmarking study that investigates integration of temporal sequencing modalities for dynamic cell state prediction. We benchmark ten integration approaches on ten datasets spanning different biological contexts, sequencing technologies, and species. We find that integrated data more accurately infers biological trajectories and achieves increased performance on classifying cells according to perturbation and disease states. Furthermore, we show that simple concatenation of spliced and unspliced molecules performs consistently well on classification tasks and can be used over more memory intensive and computationally expensive methods.

CONCLUSIONS

This work illustrates how integrated temporal gene expression modalities may be leveraged for predicting cellular trajectories and sample-associated perturbation and disease phenotypes. Additionally, this study provides users with practical recommendations for task-specific integration of single-cell gene expression modalities.

Single-cell analyses of axolotl telencephalon organization, neurogenesis, and regeneration

Salamanders are tetrapod models to study brain organization and regeneration; however, the identity and evolutionary conservation of brain cell types are largely unknown. We delineated the cell populations in the axolotl telencephalon during homeostasis and regeneration using single-cell genomic profiling. We identified glutamatergic neurons with similarities to amniote neurons of hippocampus, dorsal and lateral cortex, and conserved γ-aminobutyric acid-releasing (GABAergic) neuron classes. We inferred transcriptional dynamics and gene regulatory relationships of postembryonic, region-specific neurogenesis and unraveled conserved differentiation signatures. After brain injury, ependymoglia activate an injury-specific state before reestablishing lost neuron populations and axonal connections. Together, our analyses yield insights into the organization, evolution, and regeneration of a tetrapod nervous system.

Single-cell transcriptomics identifies pathogenic T-helper 17.1 cells and pro-inflammatory monocytes in immune checkpoint inhibitor-related pneumonitis

Background

Immune checkpoint inhibitor (ICI)-related pneumonitis is the most frequent fatal immune-related adverse event associated with programmed cell death protein-1/programmed death ligand-1 blockade. The pathophysiology however remains largely unknown, owing to limited and contradictory findings in existing literature pointing at either T-helper 1 or T-helper 17-mediated autoimmunity. In this study, we aimed to gain novel insights into the mechanisms of ICI-related pneumonitis, thereby identifying potential therapeutic targets.

Methods

In this prospective observational study, single-cell RNA and T-cell receptor sequencing was performed on bronchoalveolar lavage fluid of 11 patients with ICI-related pneumonitis and 6 demographically-matched patients with cancer without ICI-related pneumonitis. Single-cell transcriptomic immunophenotyping and cell fate mapping coupled to T-cell receptor repertoire analyses were performed.

Results

We observed enrichment of both CD4+ and CD8+ T cells in ICI-pneumonitis bronchoalveolar lavage fluid. The CD4+ T-cell compartment showed an increase of pathogenic T-helper 17.1 cells, characterized by high co-expression of TBX21 (encoding T-bet) and RORC (ROR-γ), IFN-G (IFN-γ), IL-17A, CSF2 (GM-CSF), and cytotoxicity genes. Type 1 regulatory T cells and naïve-like CD4+ T cells were also enriched. Within the CD8+ T-cell compartment, mainly effector memory T cells were increased. Correspondingly, myeloid cells in ICI-pneumonitis bronchoalveolar lavage fluid were relatively depleted of anti-inflammatory resident alveolar macrophages while pro-inflammatory ‘M1-like’ monocytes (expressing TNF, IL-1B, IL-6, IL-23A, and GM-CSF receptor CSF2RA, CSF2RB) were enriched compared with control samples. Importantly, a feedforward loop, in which GM-CSF production by pathogenic T-helper 17.1 cells promotes tissue inflammation and IL-23 production by pro-inflammatory monocytes and vice versa, has been well characterized in multiple autoimmune disorders but has never been identified in ICI-related pneumonitis.

Conclusions

Using single-cell transcriptomics, we identified accumulation of pathogenic T-helper 17.1 cells in ICI-pneumonitis bronchoalveolar lavage fluid—a phenotype explaining previous divergent findings on T-helper 1 versus T-helper 17 involvement in ICI-pneumonitis—, putatively engaging in detrimental crosstalk with pro-inflammatory ‘M1-like’ monocytes. This finding yields several novel potential therapeutic targets for the treatment of ICI-pneumonitis. Most notably repurposing anti-IL-23 merits further research as a potential efficacious and safe treatment for ICI-pneumonitis.

Control of cell state transitions

Understanding cell state transitions and purposefully controlling them is a longstanding challenge in biology. Here we present cell state transition assessment and regulation (cSTAR), an approach for mapping cell states, modelling transitions between them and predicting targeted interventions to convert cell fate decisions. cSTAR uses omics data as input, classifies cell states, and develops a workflow that transforms the input data into mechanistic models that identify a core signalling network, which controls cell fate transitions by influencing whole-cell networks. By integrating signalling and phenotypic data, cSTAR models how cells manoeuvre in Waddington's landscape1 and make decisions about which cell fate to adopt. Notably, cSTAR devises interventions to control the movement of cells in Waddington's landscape. Testing cSTAR in a cellular model of differentiation and proliferation shows a high correlation between quantitative predictions and experimental data. Applying cSTAR to different types of perturbation and omics datasets, including single-cell data, demonstrates its flexibility and scalability and provides new biological insights. The ability of cSTAR to identify targeted perturbations that interconvert cell fates will enable designer approaches for manipulating cellular development pathways and mechanistically underpinned therapeutic interventions.

Probing cell identity hierarchies by fate titration and collision during direct reprogramming

Despite the therapeutic promise of direct reprogramming, basic principles concerning fate erasure and the mechanisms to resolve cell identity conflicts remain unclear. To tackle these fundamental questions, we established a single‐cell protocol for the simultaneous analysis of multiple cell fate conversion events based on combinatorial and traceable reprogramming factor expression: Collide‐seq. Collide‐seq revealed the lack of a common mechanism through which fibroblast‐specific gene expression loss is initiated. Moreover, we found that the transcriptome of converting cells abruptly changes when a critical level of each reprogramming factor is attained, with higher or lower levels not contributing to major changes. By simultaneously inducing multiple competing reprogramming factors, we also found a deterministic system, in which titration of fates against each other yields dominant or colliding fates. By investigating one collision in detail, we show that reprogramming factors can disturb cell identity programs independent of their ability to bind their target genes. Taken together, Collide‐seq has shed light on several fundamental principles of fate conversion that may aid in improving current reprogramming paradigms.

A bifurcation concept for B-lymphoid/plasmacytoid dendritic cells with largely fluctuating transcriptome dynamics

Hematopoiesis was considered a hierarchical stepwise process but was revised to a continuous process following single-cell RNA sequencing. However, the uncertainty or fluctuation of single-cell transcriptome dynamics during differentiation was not considered, and the dendritic cell (DC) pathway in the lymphoid context remains unclear. Here, we identify human B-plasmacytoid DC (pDC) bifurcation as large fluctuating transcriptome dynamics in the putative B/NK progenitor region by dry and wet methods. By converting splicing kinetics into diffusion dynamics in a deep generative model, our original computational methodology reveals strong fluctuation at B/pDC bifurcation in IL-7Rα⁺ regions, and LFA-1 fluctuates positively in the pDC direction at the bifurcation. These expectancies are validated by the presence of B/pDC progenitors in the IL-7Rα⁺ fraction and preferential expression of LFA-1 in pDC-biased progenitors with a niche-like culture system. We provide a model of fluctuation-based differentiation, which reconciles continuous and discrete models and is applicable to other developmental systems.

Club cells employ regeneration mechanisms during lung tumorigenesis

The high plasticity of lung epithelial cells, has for many years, confounded the correct identification of the cell-of-origin of lung adenocarcinoma (LUAD), one of the deadliest malignancies worldwide. Here, we employ lineage-tracing mouse models to investigate the cell of origin of Eml4-Alk LUAD, and show that Club and Alveolar type 2 (AT2) cells give rise to tumours. We focus on Club cell originated tumours and find that Club cells experience an epigenetic switch by which they lose their lineage fidelity and gain an AT2-like phenotype after oncogenic transformation. Single-cell transcriptomic analyses identified two trajectories of Club cell evolution which are similar to the ones used during lung regeneration, suggesting that lung epithelial cells leverage on their plasticity and intrinsic regeneration mechanisms to give rise to a tumour. Together, this study highlights the role of Club cells in LUAD initiation, identifies the mechanism of Club cell lineage infidelity, confirms the presence of these features in human tumours, and unveils key mechanisms conferring LUAD heterogeneity.

Heme-stress activated NRF2 skews fate trajectories of bone marrow cells from dendritic cells towards red pulp-like macrophages in hemolytic anemia

Heme is an erythrocyte-derived toxin that drives disease progression in hemolytic anemias, such as sickle cell disease. During hemolysis, specialized bone marrow-derived macrophages with a high heme-metabolism capacity orchestrate disease adaptation by removing damaged erythrocytes and heme-protein complexes from the blood and supporting iron recycling for erythropoiesis. Since chronic heme-stress is noxious for macrophages, erythrophagocytes in the spleen are continuously replenished from bone marrow-derived progenitors. Here, we hypothesized that adaptation to heme stress progressively shifts differentiation trajectories of bone marrow progenitors to expand the capacity of heme-handling monocyte-derived macrophages at the expense of the homeostatic generation of dendritic cells, which emerge from shared myeloid precursors. This heme-induced redirection of differentiation trajectories may contribute to hemolysis-induced secondary immunodeficiency. We performed single-cell RNA-sequencing with directional RNA velocity analysis of GM-CSF-supplemented mouse bone marrow cultures to assess myeloid differentiation under heme stress. We found that heme-activated NRF2 signaling shifted the differentiation of bone marrow cells towards antioxidant, iron-recycling macrophages, suppressing the generation of dendritic cells in heme-exposed bone marrow cultures. Heme eliminated the capacity of GM-CSF-supplemented bone marrow cultures to activate antigen-specific CD4 T cells. The generation of functionally competent dendritic cells was restored by NRF2 loss. The heme-induced phenotype of macrophage expansion with concurrent dendritic cell depletion was reproduced in hemolytic mice with sickle cell disease and spherocytosis and associated with reduced dendritic cell functions in the spleen. Our data provide a novel mechanistic underpinning of hemolytic stress as a driver of hyposplenism-related secondary immunodeficiency.

Single-nuclei isoform RNA sequencing unlocks barcoded exon connectivity in frozen brain tissue

Single-nuclei RNA sequencing characterizes cell types at the gene level. However, compared to single-cell approaches, many single-nuclei cDNAs are purely intronic, lack barcodes and hinder the study of isoforms. Here we present single-nuclei isoform RNA sequencing (SnISOr-Seq). Using microfluidics, PCR-based artifact removal, target enrichment and long-read sequencing, SnISOr-Seq increased barcoded, exon-spanning long reads 7.5-fold compared to naive long-read single-nuclei sequencing. We applied SnISOr-Seq to adult human frontal cortex and found that exons associated with autism exhibit coordinated and highly cell-type-specific inclusion. We found two distinct combination patterns: those distinguishing neural cell types, enriched in TSS-exon, exon-polyadenylation-site and non-adjacent exon pairs, and those with multiple configurations within one cell type, enriched in adjacent exon pairs. Finally, we observed that human-specific exons are almost as tightly coordinated as conserved exons, implying that coordination can be rapidly established during evolution. SnISOr-Seq enables cell-type-specific long-read isoform analysis in human brain and in any frozen or hard-to-dissociate sample.

Single-cell analysis of human basal cell carcinoma reveals novel regulators of tumor growth and the tumor microenvironment

How basal cell carcinoma (BCC) interacts with its tumor microenvironment to promote growth is unclear. We use singe-cell RNA sequencing to define the human BCC ecosystem and discriminate between normal and malignant epithelial cells. We identify spatial biomarkers of tumors and their surrounding stroma that reinforce the heterogeneity of each tissue type. Combining pseudotime, RNA velocity-PAGA, cellular entropy, and regulon analysis in stromal cells reveals a cancer-specific rewiring of fibroblasts, where STAT1, TGF-β, and inflammatory signals induce a noncanonical WNT5A program that maintains the stromal inflammatory state. Cell-cell communication modeling suggests that tumors respond to the sudden burst of fibroblast-specific inflammatory signaling pathways by producing heat shock proteins, whose expression we validated in situ. Last, dose-dependent treatment with an HSP70 inhibitor suppresses in vitro vismodegib-resistant BCC cell growth, Hedgehog signaling, and in vivo tumor growth in a BCC mouse model, validating HSP70's essential role in tumor growth and reinforcing the critical nature of tumor microenvironment cross-talk in BCC progression.

Mapping the developing human immune system across organs

Single-cell genomics studies have decoded the immune-cell composition of several human prenatal organs but were limited in understanding the developing immune system as a distributed network across tissues. We profiled nine prenatal tissues combining single-cell RNA sequencing, antigen-receptor sequencing, and spatial transcriptomics to reconstruct the developing human immune system. This revealed the late acquisition of immune effector functions by myeloid and lymphoid cell subsets and the maturation of monocytes and T cells prior to peripheral tissue seeding. Moreover, we uncovered system-wide blood and immune cell development beyond primary hematopoietic organs, characterized human prenatal B1 cells, and shed light on the origin of unconventional T cells. Our atlas provides both valuable data resources and biological insights that will facilitate cell engineering, regenerative medicine, and disease understanding.

Developmental landscape of human forebrain at a single-cell level identifies early waves of oligodendrogenesis

Oligodendrogenesis in the human central nervous system has been observed mainly at the second trimester of gestation, a much later developmental stage compared to oligodendrogenesis in mice. Here, we characterize the transcriptomic neural diversity in the human forebrain at post-conception weeks (PCW) 8-10. Using single-cell RNA sequencing, we find evidence of the emergence of a first wave of oligodendrocyte lineage cells as early as PCW 8, which we also confirm at the epigenomic level through the use of single-cell ATAC-seq. Using regulatory network inference, we predict key transcriptional events leading to the specification of oligodendrocyte precursor cells (OPCs). Moreover, by profiling the spatial expression of 50 key genes through the use of in situ sequencing (ISS), we identify regions in the human ventral fetal forebrain where oligodendrogenesis first occurs. Our results indicate evolutionary conservation of the first wave of oligodendrogenesis between mice and humans and describe regulatory mechanisms involved in human OPC specification.