Trajectory-based differential expression analysis for single-cell sequencing data

Cyclin O controls entry into the cell-cycle variant required for multiciliated cell differentiation

Multiciliated cells (MCCs) ensure fluid circulation in various organs. Their differentiation is marked by the amplification of cilia-nucleating centrioles, driven by a genuine cell-cycle variant, which is characterized by wave-like expression of canonical and non-canonical cyclins such as Cyclin O (CCNO). Patients with CCNO mutations exhibit a subtype of primary ciliary dyskinesia called reduced generation of motile cilia (RGMC). Here, we show that Ccno is activated at the crossroads of the onset of MCC differentiation, the entry into the MCC cell-cycle variant, and the activation of the centriole biogenesis program. Its absence blocks the G1/S-like transition of the cell-cycle variant, interrupts the centriologenesis transcription program, and compromises the production of centrioles and cilia in mouse brain and human respiratory MCCs. Altogether, our study identifies CCNO as a core regulator of entry into the MCC cell-cycle variant and the interruption of this variant as one etiology of RGMC.

Cyclin switch tailors a cell cycle variant to orchestrate multiciliogenesis

Meiosis, endoreplication, and asynthetic fissions are variations of the canonical cell cycle where either replication or mitotic divisions are muted. Here, we identify a cell cycle variantconserved across organs and mammals, where both replication and mitosis are muted, and that orchestrates the differentiation of post-mitotic progenitors into multiciliated cells (MCCs). MCC progenitors reactivate most of the cell cycle transcriptional program but replace the temporal expression of cyclins E2 and A2 with non-canonical cyclins O and A1. In addition, the primary APC/C inhibitor Emi1 is silenced. Re-expressing cyclins E2 and A2 and/or Emi1 can induce partial replication or mitosis. This shows that a cell can co-opt the cell cycle genetic program and regulate only certain elements to qualitatively and quantitatively divert CDK activity toward differentiation rather than division. We propose this cell cycle variant to exploit the existence of a cytoplasmic-or centriolar-CDK threshold lower than the S-phase threshold.

TIME-CoExpress: Temporal Trajectory Modeling of Dynamic Gene Co-expression Patterns Using Single-Cell Transcriptomics Data

The rapid advancements of single-cell RNA sequencing (scRNAseq) technology provide high-resolution views of transcriptomic activity within a single cell. Most routine analyses of scRNAseq data focus on individual genes; however, the one-gene-at-a-time analysis is likely to miss meaningful genetic interactions. Gene co-expression analysis addresses this issue by identifying coordinated gene expression changes in response to cellular conditions, such as developmental or temporal trajectory. Identifying differential co-expression gene combinations along the cell temporal trajectory using scRNAseq data can provide deeper insight into the biological processes. Existing approaches for gene co-expression analysis assume a restrictive linear change of gene co-expression. In this paper, we propose a copula-based approach with proper data-driven smoothing functions to model non-linear gene co-expression changes along cellular temporal trajectories. Our proposed approach provides flexibility to incorporate characteristics such as over-dispersion and zero-inflation rate observed in scRNAseq data into the modeling framework. We conducted a series of simulation analyses to evaluate the performance of the proposed algorithm. We demonstrate the implementation of the proposed algorithm using a scRNAseq dataset and identify differential co-expression gene pairs along cell temporal trajectory in pituitary embryonic development comparing Nxn - / - mutated versus wild-type mice.

DUX4-stimulated genes define the antiviral response to herpesviruses in human trophoblasts

The placenta serves as a barrier against the vertical transmission of viruses, due to the protective functions of fetal-derived trophoblasts. Although some antiviral programs of trophoblasts are well-documented, our understanding of how trophoblasts respond to teratogenic viruses remains incomplete. To address this, we profiled the transcriptional responses of human trophoblast organoids to seven teratogenic viruses. We discovered that herpesviruses including HSV-1, HSV-2, and HCMV did not trigger an interferon (IFN) response but instead induced the expression of DUX4 and downstream target genes, termed DUX4-stimulated genes (DSGs). This program was uniquely, highly expressed in trophoblasts. Through single-cell RNA sequencing, we defined the trophoblast response to DUX4 and demonstrated that this program defines cells with low viral transcripts following HSV-1 infection. Similar to the IFN-stimulated genes, we observed that many DSGs with diverse predicted functions exhibited anti-herpesvirus activity. These findings establish DUX4 and DSGs as a critical antiviral defense mechanism in trophoblasts.

Fibroblastic reticular cells generate protective intratumoral T cell environments in lung cancer

Stringent control of T cell activity in the tumor microenvironment is essential for the generation of protective antitumor immunity. However, the identity, differentiation, and functions of the cells that create critical fibroblastic niches promoting tumor-infiltrating T cells remain elusive. Here, we show that CCL19-expressing fibroblastic reticular cells (FRCs) generate interconnected T cell environments (TEs) in human non-small cell lung cancer, including tertiary lymphoid structures and T cell tracks. Analysis of the FRC-T cell interactome in TEs indicated molecular networks regulating niche-specific differentiation of CCL19-expressing fibroblasts and T cell activation pathways. Single-cell transcriptomics and cell fate-mapping analyses in mice confirmed that FRCs in TEs originate from mural and adventitial progenitors. Ablation of intratumoral FRC precursors decreased antitumor T cell activity, resulting in reduced tumor control during coronavirus vector-based immunotherapy. In summary, specialized FRC niches in the tumor microenvironment govern the quality and extent of antitumor T cell immunity.

Molecular and cellular dynamics of the developing human neocortex

The development of the human neocortex is highly dynamic, involving complex cellular trajectories controlled by gene regulation1. Here we collected paired single-nucleus chromatin accessibility and transcriptome data from 38 human neocortical samples encompassing both the prefrontal cortex and the primary visual cortex. These samples span five main developmental stages, ranging from the first trimester to adolescence. In parallel, we performed spatial transcriptomic analysis on a subset of the samples to illustrate spatial organization and intercellular communication. This atlas enables us to catalogue cell-type-specific, age-specific and area-specific gene regulatory networks underlying neural differentiation. Moreover, combining single-cell profiling, progenitor purification and lineage-tracing experiments, we have untangled the complex lineage relationships among progenitor subtypes during the neurogenesis-to-gliogenesis transition. We identified a tripotential intermediate progenitor subtype-tripotential intermediate progenitor cells (Tri-IPCs)-that is responsible for the local production of GABAergic neurons, oligodendrocyte precursor cells and astrocytes. Notably, most glioblastoma cells resemble Tri-IPCs at the transcriptomic level, suggesting that cancer cells hijack developmental processes to enhance growth and heterogeneity. Furthermore, by integrating our atlas data with large-scale genome-wide association study data, we created a disease-risk map highlighting enriched risk associated with autism spectrum disorder in second-trimester intratelencephalic neurons. Our study sheds light on the molecular and cellular dynamics of the developing human neocortex.

Hepatic stellate cell single cell atlas reveals a highly similar activation process across liver disease aetiologies

Background & Aims

The progression of chronic liver disease (CLD) is characterized by excessive extracellular matrix deposition, disrupting hepatic architecture and function. Upon liver injury, hepatic stellate cells (HSCs) differentiate towards myofibroblasts and become inflammatory, proliferative and fibrogenic. To date, it is still unclear whether HSC activation is driven by similar mechanisms in different aetiologies.

Methods

HSCs from multiple publicly available single-cell RNA-sequencing datasets were annotated and merged into a single-cell HSC activation atlas. Spheroid co-cultures of primary mouse hepatocytes/HSCs (n = 5) and ELISAs on patient plasma samples (n = 80) were performed to validate the mechanistic insight obtained from the HSC atlas.

Results

We established an HSC activation atlas in which HSCs are clearly divided into three distinct transcriptomic profiles: quiescent HSCs, initiatory HSCs and myofibroblasts. These transcriptomic profiles are present in each of the investigated mouse liver injury models as well as in human CLDs, indicating that HSC activation is a conserved process. This activation process is driven by a core set of transcription factors independent of liver injury or species. Furthermore, we reveal novel ligands associated with activation of HSCs in multiple liver injury models and validate the profibrotic effect of parathyroid hormone. Finally, we identify COLEC10 as a conserved marker for quiescent HSCs and a biomarker of liver fibrosis in patients with different CLDs (p <0.0001).

Conclusions

We reveal unexpected similarities in the regulatory mechanisms of HSCs across diverse liver injury settings and species. The HSC activation atlas has the potential to provide novel insights into liver fibrosis and steer novel treatment options.

Impact and implications

This study establishes a single-cell atlas of hepatic stellate cells across various liver injuries, highlighting a conserved activation process between different injuries and across species. The discovery of novel activating ligands and the biomarker COLEC10 in human plasma could be used to enhance diagnostic and therapeutic strategies. Additionally, the conserved activation process supports the use of any mouse model for mechanistic studies and testing of new anti-fibrotic compounds, streamlining preclinical research efforts.

Gene-level alignment of single-cell trajectories

Single-cell data analysis can infer dynamic changes in cell populations, for example across time, space or in response to perturbation, thus deriving pseudotime trajectories. Current approaches comparing trajectories often use dynamic programming but are limited by assumptions such as the existence of a definitive match. Here we describe Genes2Genes, a Bayesian information-theoretic dynamic programming framework for aligning single-cell trajectories. It is able to capture sequential matches and mismatches of individual genes between a reference and query trajectory, highlighting distinct clusters of alignment patterns. Across both real world and simulated datasets, it accurately inferred alignments and demonstrated its utility in disease cell-state trajectory analysis. In a proof-of-concept application, Genes2Genes revealed that T cells differentiated in vitro match an immature in vivo state while lacking expression of genes associated with TNF signaling. This demonstrates that precise trajectory alignment can pinpoint divergence from the in vivo system, thus guiding the optimization of in vitro culture conditions.

Trajectory inference from single-cell genomics data with a process time model

Single-cell transcriptomics experiments provide gene expression snapshots of heterogeneous cell populations across cell states. These snapshots have been used to infer trajectories and dynamic information even without intensive, time-series data by ordering cells according to gene expression similarity. However, while single-cell snapshots sometimes offer valuable insights into dynamic processes, current methods for ordering cells are limited by descriptive notions of "pseudotime" that lack intrinsic physical meaning. Instead of pseudotime, we propose inference of "process time" via a principled modeling approach to formulating trajectories and inferring latent variables corresponding to timing of cells subject to a biophysical process. Our implementation of this approach, called Chronocell, provides a biophysical formulation of trajectories built on cell state transitions. The Chronocell model is identifiable, making parameter inference meaningful. Furthermore, Chronocell can interpolate between trajectory inference, when cell states lie on a continuum, and clustering, when cells cluster into discrete states. By using a variety of datasets ranging from cluster-like to continuous, we show that Chronocell enables us to assess the suitability of datasets and reveals distinct cellular distributions along process time that are consistent with biological process times. We also compare our parameter estimates of degradation rates to those derived from metabolic labeling datasets, thereby showcasing the biophysical utility of Chronocell. Nevertheless, based on performance characterization on simulations, we find that process time inference can be challenging, highlighting the importance of dataset quality and careful model assessment.

Single-cell atlas of Leishmania development in sandflies reveals the heterogeneity of transmitted parasites and their role in infection

Sandfly vectors transmit Leishmania through egestion of parasites into the host skin. The transmissible dose is shaped by Leishmania development in the sandfly gut, described as a sequential differentiation of promastigote morphotypes. Apart from isolated mammal-infective metacyclic promastigotes, little is known about the transcriptional programs and molecular markers for other stages coinhabiting the midgut in mature infections and cotransmitted by the sandfly bite. Here, we elucidate the single-cell transcriptomic complexity of Leishmania major colonizing its natural vector Phlebotomus duboscqi at early and late infection, providing markers for three transmissible stages. Contrary to prevailing models, our analyses indicate a nonlinear developmental progression, with bifurcation into either replicating early metacyclics or attached and detached haptomonads. We demonstrate that haptomonads constitute a key component of the transmitted inoculum and, along with nondividing late metacyclics, are infectious to and exacerbate the pathology in mice. Our single-cell analysis and validated markers will facilitate further studies on the Leishmania life cycle and its implications for vector-to-host transmission dynamics.

Single-cell analysis reveals the loss of FABP4-positive proliferating valvular endothelial cells relates to functional mitral regurgitation

BACKGROUND

Functional mitral regurgitation (MR) is a common form of mitral valve dysfunction that often persists even after surgical intervention, requiring reoperation in some cases. To advance our understanding of the pathogenesis of functional MR, it is crucial to characterize the cellular composition of the mitral valve leaflet and identify molecular changes in each cell subtype within the mitral valves of MR patients. Therefore, we aimed to comprehensively examine the cellular and molecular components of mitral valves in patients with MR.

METHODS

We conducted a single-cell RNA sequencing (scRNA-seq) analysis of mitral valve leaflets extracted from six patients who underwent heart transplantation. The cohort comprised three individuals with moderate-to-severe functional MR (MR group) and three non-diseased controls (NC group). Bioinformatics was applied to identify cell types, delineate cell functions, and explore cellular developmental trajectories and interactions. Key findings from the scRNA-seq analysis were validated using pathological staining to visualize key markers in the mitral valve leaflets. Additionally, in vitro experiments with human primary valvular endothelial cells were conducted to further support our results.

RESULTS

Our study revealed that valve interstitial cells are critical for adaptive valve remodelling, as they secrete extracellular matrix proteins and promote fibrosis. We discovered an abnormal decrease in a subpopulation of FABP4 (fatty acid binding protein 4)-positive proliferating valvular endothelial cells. The trajectory analysis identifies this subcluster as the origin of VECs. Immunohistochemistry on the expanded cohort showed a reduction of FABP4-positive VECs in patients with functional MR. Intervention experiments with primary cells indicated that FABP4 promotes proliferation and migration in mitral valve VECs and enhances TGFβ-induced differentiation.

CONCLUSIONS

Our study presented a comprehensive assessment of the mitral valve cellular landscape of patients with MR and sheds light on the molecular changes occurring in human mitral valves during functional MR. We found a notable reduction in the proliferating endothelial cell subpopulation of valve leaflets, and FABP4 was identified as one of their markers. Therefore, FABP4 positive VECs served as proliferating endothelial cells relates to functional mitral regurgitation. These VECs exhibited high proliferative and differentiative properties. Their reduction was associated with the occurrence of functional MR.

IRE1α silences dsRNA to prevent taxane-induced pyroptosis in triple-negative breast cancer

Chemotherapy is often combined with immune checkpoint inhibitor (ICIs) to enhance immunotherapy responses. Despite the approval of chemo-immunotherapy in multiple human cancers, many immunologically cold tumors remain unresponsive. The mechanisms determining the immunogenicity of chemotherapy are elusive. Here, we identify the ER stress sensor IRE1α as a critical checkpoint that restricts the immunostimulatory effects of taxane chemotherapy and prevents the innate immune recognition of immunologically cold triple-negative breast cancer (TNBC). IRE1α RNase silences taxane-induced double-stranded RNA (dsRNA) through regulated IRE1-dependent decay (RIDD) to prevent NLRP3 inflammasome-dependent pyroptosis. Inhibition of IRE1α in Trp53-/- TNBC allows taxane to induce extensive dsRNAs that are sensed by ZBP1, which in turn activates NLRP3-GSDMD-mediated pyroptosis. Consequently, IRE1α RNase inhibitor plus taxane converts PD-L1-negative, ICI-unresponsive TNBC tumors into PD-L1high immunogenic tumors that are hyper-sensitive to ICI. We reveal IRE1α as a cancer cell defense mechanism that prevents taxane-induced danger signal accumulation and pyroptotic cell death.

Molecular, cellular, and developmental organization of the mouse vomeronasal organ at single cell resolution

We have generated single cell transcriptomic atlases of vomeronasal organs (VNO) from juvenile and adult mice. Combined with spatial molecular imaging, we uncover a distinct, previously unidentified class of cells that express the vomeronasal receptors (VRs) and a population of canonical olfactory sensory neurons in the VNO. High-resolution trajectory and cluster analyses reveal the lineage relationship, spatial distribution of cell types, and a putative cascade of molecular events that specify the V1r, V2r, and OR lineages from a common stem cell population. The expression of vomeronasal and olfactory receptors follow power law distributions, but there is high variability in average expression levels between individual receptor and cell types. Substantial co-expression is found between receptors across clades, from different classes, and between olfactory and VRs, with nearly half from pairs located on the same chromosome. Interestingly, the expression of V2r, but not V1r, genes is associated with various transcription factors, suggesting distinct mechanisms of receptor choice associated with the two cell types. We identify association between transcription factors, surface axon guidance molecules, and individual VRs, thereby uncovering a molecular code that guides the specification of the vomeronasal circuitry. Our study provides a wealth of data on the development and organization of the accessory olfactory system at both cellular and molecular levels to enable a deeper understanding of vomeronasal system function.

MAIT cell plasticity enables functional adaptation that drives antibacterial immune protection

Mucosal-associated invariant T (MAIT) cells are known for their rapid effector functions and antibacterial immune protection. Here, we define the plasticity of interferon-γ (IFN-γ)-producing MAIT1 and interleukin-17A (IL-17A)-producing MAIT17 cell subsets in vivo. Whereas T-bet+ MAIT1 cells remained stable in all experimental settings, after adoptive transfer or acute Legionella or Francisella infection, RORγt+ MAIT17 cells could undergo phenotypic and functional conversion into both RORγt+T-bet+ MAIT1/17 and RORγt-T-bet+ MAIT1 cells. This plasticity ensured that MAIT17 cells played a dominant role in generating antibacterial MAIT1 responses in mucosal tissues. Single-cell transcriptomics revealed that MAIT17-derived MAIT1 cells were distinct from canonical MAIT1 cells yet could migrate out of mucosal tissues to contribute to the global MAIT1 pool in subsequent systemic infections. Human IL-17A-secreting MAIT cells also showed similar functional plasticity. Our findings have broad implications for understanding the role of MAIT cells in combatting infections and their potential utility in MAIT cell-targeted vaccines.

Multiomics Research: Principles and Challenges in Integrated Analysis

Multiomics research is a transformative approach in the biological sciences that integrates data from genomics, transcriptomics, proteomics, metabolomics, and other omics technologies to provide a comprehensive understanding of biological systems. This review elucidates the fundamental principles of multiomics, emphasizing the necessity of data integration to uncover the complex interactions and regulatory mechanisms underlying various biological processes. We explore the latest advances in computational methodologies, including deep learning, graph neural networks (GNNs), and generative adversarial networks (GANs), which facilitate the effective synthesis and interpretation of multiomics data. Additionally, this review addresses the critical challenges in this field, such as data heterogeneity, scalability, and the need for robust, interpretable models. We highlight the potential of large language models to enhance multiomics analysis through automated feature extraction, natural language generation, and knowledge integration. Despite the important promise of multiomics, the review acknowledges the substantial computational resources required and the complexity of model tuning, underscoring the need for ongoing innovation and collaboration in the field. This comprehensive analysis aims to guide researchers in navigating the principles and challenges of multiomics research to foster advances in integrative biological analysis.

Integration of spatial multiplexed protein imaging and transcriptomics in the human kidney tracks the regenerative potential timeline of proximal tubules

The organizational principles of nephronal segments are based on longstanding anatomical and physiological attributes that are closely linked to the homeostatic functions of the kidney. Novel molecular approaches have recently uncovered layers of deeper signatures and states in tubular cells that arise at various timepoints on the spectrum between health and disease. For example, a dedifferentiated state of proximal tubular cells with mesenchymal stemness markers is frequently seen after injury. The persistence of such a state is associated with failed repair. Here, we introduce a novel analytical pipeline applied to highly multiplexed spatial protein imaging to characterize proximal tubular subpopulations and neighborhoods in reference and disease human kidney tissue. The results were validated and extended through integration with spatial and single cell transcriptomics. We demonstrate that, in reference tissue, a large proportion of S1 and S2 proximal tubular epithelial cells express THY1, a mesenchymal stromal and stem cell marker that regulates differentiation. Kidney disease is associated with loss of THY1 and transition towards expression of PROM1, another stem cell marker shown recently to be linked to failed repair. We demonstrate that the trajectory of proximal tubular cells to THY1 expression is clearly distinct from that of PROM1, and that a state with PROM1 expression is associated with niches of inflammation. Our data support a model in which the interplay between THY1 and PROM1 expression in proximal tubules associates with their regenerative potential and marks the timeline of disease progression.

Chromatin remodelling in damaged intestinal crypts orchestrates redundant TGFβ and Hippo signalling to drive regeneration

Cell state dynamics underlying successful tissue regeneration are undercharacterized. In the intestine, damage prompts epithelial reprogramming into revival stem cells (revSCs) that reconstitute Lgr5+ intestinal stem cells (ISCs). Here single-nuclear multi-omics of mouse crypts regenerating from irradiation shows revSC chromatin accessibility overlaps with ISCs and differentiated lineages. While revSC genes themselves are accessible throughout homeostatic epithelia, damage-induced remodelling of chromatin in the crypt converges on Hippo and the transforming growth factor-beta (TGFβ) signalling pathway, which we show is transiently activated and directly induces functional revSCs. Combinatorial gene expression analysis further suggests multiple sources of revSCs, and we demonstrate TGFβ can reprogramme enterocytes, goblet and paneth cells into revSCs and show individual revSCs form organoids. Despite this, loss of TGFβ signalling yields mild regenerative defects, whereas interference in both Hippo and TGFβ leads to profound defects and death. Intestinal regeneration is thus poised for activation by a compensatory system of crypt-localized, transient morphogen cues that support epithelial reprogramming and robust intestinal repair.

Conserved B cell signaling, activation, and differentiation in porcine jejunal and ileal Peyer's patches despite distinct immune landscapes

Peyer's patches (PPs) are B cell-rich sites of intestinal immune induction, yet PP-associated B cell signaling, activation, and differentiation are poorly defined. Single-cell and spatial transcriptomics were completed to study B cells from porcine jejunum and ileum containing PPs. Intestinal locations had distinct immune landscapes, including more follicular B cells in ileum and increased MHC-II-encoding gene expression in jejunal B cells. Despite distinct landscapes, conserved B cell dynamics were detected across intestinal locations, including B cell signaling to CD4+ macrophages that are putative phagocytic, cytotoxic, effector cells and deduced routes of B cell activation/differentiation, including resting B cells migrating into follicles to replicate/divide or differentiate into antibody-secreting cells residing in intestinal crypts. A six-biomarker panel recapitulated transcriptomics findings of B cell phenotypes, frequencies, and spatial locations via ex vivo and in situ staining. Findings convey conserved B cell dynamics across intestinal locations containing PPs, despite location-specific immune environments. Results establish a benchmark of B cell dynamics for understanding intestinal immune induction important to promoting gut/overall health.

Ionizing radiation inhibits zebrafish embryo hatching through induction of tissue inhibitors of metalloproteinases (TIMPs) expression

Ionizing radiation (IR) has garnered growing attention because of its biological effects on aquatic organisms and humans. Here, we identify the most impacted organs and uncover the molecular mechanisms causing the changes in the context of vertebrate development using single-cell RNA sequencing. Alterations in cellular composition and biological functions were explored using transcriptomic profiling of zebrafish embryos exposed to 5 Gy. Single-cell RNA sequencing analyses unveiled notable shifts in the proportions of brain/central nervous system and hatching gland clusters. Although IR exposure led to increased expression of hatching enzymes, a significant but mild delay in hatching was observed following 5 Gy IR exposure. Gene Ontology analysis showed an increased expression of tissue inhibitors of metalloproteinases (TIMPs), known as matrix metalloproteinase inhibitors, which was confirmed via whole-mount in situ hybridization. Correlation analysis linked TIMPs to transcription factors cebpb and cebpd, which were significantly correlated post-IR exposure. Although no morphological changes were observed in some organs, including the brain, the study reveals substantial alterations in developing vertebrates. Notably, despite increased hatching enzymes, elevated TIMPs in the hatching gland suggest a regulatory mechanism impacting hatching activity. This research contributes to comprehending the ecological repercussions of IR exposure, emphasizing the importance of safety measures for aquatic ecosystems and overall environmental health.

A multimodal zebrafish developmental atlas reveals the state-transition dynamics of late-vertebrate pluripotent axial progenitors

Elucidating organismal developmental processes requires a comprehensive understanding of cellular lineages in the spatial, temporal, and molecular domains. In this study, we introduce Zebrahub, a dynamic atlas of zebrafish embryonic development that integrates single-cell sequencing time course data with lineage reconstructions facilitated by light-sheet microscopy. This atlas offers high-resolution and in-depth molecular insights into zebrafish development, achieved through the sequencing of individual embryos across ten developmental stages, complemented by reconstructions of cellular trajectories. Zebrahub also incorporates an interactive tool to navigate the complex cellular flows and lineages derived from light-sheet microscopy data, enabling in silico fate-mapping experiments. To demonstrate the versatility of our multimodal resource, we utilize Zebrahub to provide fresh insights into the pluripotency of neuro-mesodermal progenitors (NMPs) and the origins of a joint kidney-hemangioblast progenitor population.

Alveolar regeneration by airway secretory-cell-derived p63+ progenitors

Lung injury activates epithelial stem or progenitor cells for alveolar repair and regeneration. Unraveling the origin and fate of injury-induced progenitors is crucial for elucidating lung repair mechanisms. Here, we report that p63-expressing progenitors emerge upon bleomycin-induced mouse lung injury. Single-cell RNA sequencing and clonal analysis reveal that these p63+ progenitors proliferate rapidly and differentiate into alveolar type 1 and type 2 cells through different trajectories. Dual recombinase-mediated sequential genetic-lineage tracing demonstrates that p63+ progenitors originate from airway secretory cells and subsequently generate alveolar cells. Functionally, p63 activation is essential for efficient alveolar regeneration from secretory cells post injury. Our study identifies secretory-cell-derived p63+ progenitors as contributors to alveolar repair, suggesting a potential therapeutic avenue for lung regeneration following injury.

GraphPCA: a fast and interpretable dimension reduction algorithm for spatial transcriptomics data

The rapid advancement of spatial transcriptomics technologies has revolutionized our understanding of cell heterogeneity and intricate spatial structures within tissues and organs. However, the high dimensionality and noise in spatial transcriptomic data present significant challenges for downstream data analyses. Here, we develop GraphPCA, an interpretable and quasi-linear dimension reduction algorithm that leverages the strengths of graphical regularization and principal component analysis. Comprehensive evaluations on simulated and multi-resolution spatial transcriptomic datasets generated from various platforms demonstrate the capacity of GraphPCA to enhance downstream analysis tasks including spatial domain detection, denoising, and trajectory inference compared to other state-of-the-art methods.

Effects of intergenerational transmission of small intestinal bacteria cultured from stunted Bangladeshi children with enteropathy

Environmental enteric dysfunction (EED), a small intestinal disorder found at a high prevalence in stunted children, is associated with gut mucosal barrier disruption and decreased absorptive capacity due to reduced intact small intestinal villi1-4. To test the hypothesis that intergenerational transmission of a perturbed small intestinal microbiota contributes to undernutrition by inducing EED5, we characterized two consortia of bacterial strains cultured from duodenal aspirates from stunted Bangladeshi children with EED - one of which induced local and systemic inflammation in gnotobiotic female mice. Offspring of dams that received this inflammatory consortium exhibited immunologic changes along their gut that phenocopied features of EED in children. Single nucleus plus bulk RNA-sequencing revealed alterations in inter-cellular signaling pathways related to intestinal epithelial cell renewal, barrier integrity and immune function while analyses of cerebral cortex disclosed alterations in glial- and endothelial-neuronal signaling pathways that regulate neural growth/axonal guidance, angiogenesis and inflammation. Analysis of ultrasonic vocalization calls in gnotobiotic P5-P9 pups indicated increased arousal and perturbed neurodevelopment in the offspring of dams harboring the inflammation-inducing consortium. Cohousing experiments and follow-up screening of candidate disease-promoting bacterial isolates identified a strain typically found in the oral microbiota (Campylobacter concisus) as a contributor to enteropathy. Given that fetal growth was also impaired in the dams with the consortium that induced enteropathy, this preclinical model allows the effects of the human small intestinal microbiota on both pre- and postnatal development to be ascertained, setting the stage for identification of small intestinal microbiota-targeted therapeutics for (intergenerational) undernutrition.

Astrocyte allocation during brain development is controlled by Tcf4-mediated fate restriction

Astrocytes in the brain exhibit regional heterogeneity contributing to regional circuits involved in higher-order brain functions, yet the mechanisms controlling their distribution remain unclear. Here, we show that the precise allocation of astrocytes to specific brain regions during development is achieved through transcription factor 4 (Tcf4)-mediated fate restriction based on their embryonic origin. Loss of Tcf4 in ventral telencephalic neural progenitor cells alters the fate of oligodendrocyte precursor cells to transient intermediate astrocyte precursor cells, resulting in mislocalized astrocytes in the dorsal neocortex. These ectopic astrocytes engage with neocortical neurons and acquire features reminiscent of dorsal neocortical astrocytes. Furthermore, Tcf4 functions as a suppressor of astrocyte fate during the differentiation of oligodendrocyte precursor cells derived from the ventral telencephalon, thereby restricting the fate to the oligodendrocyte lineage in the dorsal neocortex. Together, our findings highlight a previously unappreciated role for Tcf4 in regulating astrocyte allocation, offering additional insights into the mechanisms underlying neurodevelopmental disorders linked to Tcf4 mutations.

A polarized FGF8 source specifies frontotemporal signatures in spatially oriented cell populations of cortical assembloids

Organoids generating major cortical cell types in distinct compartments are used to study cortical development, evolution and disorders. However, the lack of morphogen gradients imparting cortical positional information and topography in current systems hinders the investigation of complex phenotypes. Here, we engineer human cortical assembloids by fusing an organizer-like structure expressing fibroblast growth factor 8 (FGF8) with an elongated organoid to enable the controlled modulation of FGF8 signaling along the longitudinal organoid axis. These polarized cortical assembloids mount a position-dependent transcriptional program that in part matches the in vivo rostrocaudal gene expression patterns and that is lost upon mutation in the FGFR3 gene associated with temporal lobe malformations and intellectual disability. By producing spatially oriented cell populations with signatures related to frontal and temporal area identity within individual assembloids, this model recapitulates in part the early transcriptional divergence embedded in the protomap and enables the study of cortical area-relevant alterations underlying human disorders.

Detecting significant expression patterns in single-cell and spatial transcriptomics with a flexible computational approach

Gene expression data holds the potential to shed light on multiple biological processes at once. However, data analysis methods for single cell sequencing mostly focus on finding cell clusters or the principal progression line of the data. Data analysis for spatial transcriptomics mostly addresses clustering and finding spatially variable genes. Existing data analysis methods are effective in finding the main data features, but they might miss less pronounced, albeit significant, processes, possibly involving a subset of the samples. In this work we present SPIRAL: Significant Process InfeRence ALgorithm. SPIRAL is based on Gaussian statistics to detect all statistically significant biological processes in single cell, bulk and spatial transcriptomics data. The algorithm outputs a list of structures, each defined by a set of genes working simultaneously in a specific population of cells. SPIRAL is unique in its flexibility: the structures are constructed by selecting subsets of genes and cells based on statistically significant and consistent differential expression. Every gene and every cell may be part of one structure, more or none. SPIRAL also provides several visual representations of structures and pathway enrichment information. We validated the statistical soundness of SPIRAL on synthetic datasets and applied it to single cell, spatial and bulk RNA-sequencing datasets. SPIRAL is available at https://spiral.technion.ac.il/ .

Age, sex, and cell type-resolved hypothalamic gene expression across the pubertal transition in mice

BACKGROUND

The hypothalamus plays a central role in regulating puberty. However, our knowledge of the postnatal gene regulatory networks that control the pubertal transition in males and females is incomplete. Here, we investigate the age-, sex- and cell-type-specific gene regulation in the hypothalamus across the pubertal transition.

METHODS

We used RNA-seq to profile hypothalamic gene expression in male and female mice at five time points spanning the onset of puberty (postnatal days (PD) 12, 22, 27, 32, and 37). By combining this data with hypothalamic single nuclei RNA-seq data from pre- and postpubertal mice, we assigned gene expression changes to their most likely cell types of origin. In our colony, pubertal onset occurs earlier in male mice, allowing us to focus on genes whose expression is dynamic across ages and offset between sexes, and to explore the bases of sex effects.

RESULTS

Our age-by-sex pattern of expression enriched for biological pathways involved hormone production, neuronal activation, and glial maturation. Additionally, we inferred a robust expansion of oligodendrocytes precursor cells into mature oligodendrocytes spanning the prepubertal (PD12) to peri-pubertal (PD27) timepoints. Using spatial transcriptomic data from postpubertal mice, we observed the lateral hypothalamic area and zona incerta were the most oligodendrocyte-rich regions and that these cells expressed genes known to be involved in pubertal regulation.

CONCLUSION

Together, by incorporating multiple biological timepoints and using sex as a variable, we identified gene and cell-type changes that may participate in orchestrating the pubertal transition and provided a resource for future studies of postnatal hypothalamic gene regulation.

Comparative analysis of rhesus macaque and human placental organoids highlights evolutionary differences in placentation

Throughout evolution, the placenta has diversified in both structure and cellular composition while maintaining its fundamental function. Trophoblasts are fetal-derived cells responsible for nourishing and protecting the developing fetus and are a universal component of all placentas. While primate placentas exhibit many shared morphological features, species-specific differences in gene expression remain largely unexplored, primarily due to the lack of suitable in vitro models. To address this gap, we developed organoids from rhesus macaque placentas, including trophoblast and maternal-derived decidua types, and compared them with organoids derived from human placentas. We performed comparative single-cell RNA sequencing to delineate shared and distinct transcriptional signatures between rhesus macaque and human trophoblasts. We further defined the differentiation trajectories leading to the multinucleated syncytiotrophoblast, highlighting both shared and unique transcriptional signatures between rhesus macaque and human placentas. This work establishes novel in vitro models of the non-human primate placenta and characterizes distinct trophoblast gene profiles in rhesus macaques and humans, providing insights into interspecies variations in placental biology.

Integrative analyses of bulk and single-cell RNA-seq reveals the correlation between SPP1+ macrophages and resistance to neoadjuvant chemoimmunotherapy in esophageal squamous cell carcinoma

Neoadjuvant chemoimmunotherapy (NACI) has significant implications for the treatment of esophageal cancer. However, its clinical efficacy varies considerably among patients, necessitating further investigation into the underlying mechanisms. The rapid advancement of single-cell RNA sequencing (scRNA-seq) technology facilitates the analysis of patient heterogeneity at the cellular level, particularly regarding treatment outcomes. In this study, we first analyzed scRNA-seq data of esophageal squamous cell carcinoma (ESCC) following NACI, obtained from the Gene Expression Omnibus (GEO) database. After performing dimensionality reduction, clustering, and annotation on the scRNA-seq data, we employed CellChat to investigate differences in cell-cell communication among samples from distinct efficacy groups. The results indicated that macrophages in the non-responder exhibited stronger cell communication intensity compared to those in responders, with SPP1 and GALECTIN signals showing the most significant differences between the two groups. This finding underscores the crucial role of macrophages in the efficacy of NACI. Subsequently, reclustering of macrophages revealed that Mac-SPP1 may be primarily responsible for treatment resistance, while Mac-C1QC appears to promote T cell activation. Finally, we conducted transcriptome sequencing on ESCC tissues obtained from 32 patients who underwent surgery following NACI. Utilizing CIBERSORT, CIBERSORTx, and WGCNA, we analyzed the heterogeneity of tumor microenvironment among different efficacy groups and validated the correlation between SPP1+ macrophages and resistance to NACI in ESCC using publicly available transcriptome sequencing datasets. These findings suggest that SPP1+ macrophages may represent a key factor contributing to resistance against NACI in ESCC.

Unraveling the ecological landscape of mast cells in esophageal cancer through single-cell RNA sequencing

Background

Esophageal cancer (EC) is a major health issue, ranking seventh in incidence and sixth in mortality worldwide. Despite advancements in multidisciplinary treatment approaches, the 5-year survival rate for EC remains low at 21%. Challenges in EC treatment arise from late-stage diagnosis, high malignancy, and poor prognosis. Understanding the tumor microenvironment is critical, as it includes various cellular and extracellular components that influence tumor behavior and treatment response. Mast cells (MCs), as tissue-resident immune cells, play dual roles in tumor dynamics. High-throughput single-cell RNA sequencing offers a powerful tool for analyzing tumor heterogeneity and immune interactions, although its application in EC is limited.

Methods

In this study, we investigated the immune microenvironment of EC using single-cell RNA sequencing and established a comprehensive immune profile. We also performed analysis of upstream transcription factors and downstream pathway enrichment to further comprehensively decipher MCs in EC. Besides, we performed knockdown experiments to explore the role of epidermal growth factor receptor (EGFR) signaling pathway in MCs-tumor cell interactions, highlighting its potential as a prognostic marker. Finally, we constructed a prognostic model for EC, which provided valuable suggestions for the diagnosis and prognosis of EC.

Results

Our analysis identified 11 major cell types, of which MCs were particularly present in pericarcinoma tissues. Further grouping of the 5,001 MCs identified 8 distinct subtypes, including SRSF7-highly expressed MCs, which showed strong tumor preference and potential tumor-promoting properties. Moreover, we identified the key signaling receptor EGFR and validated it by in vitro knockdown experiments, demonstrating its cancer-promoting effects. In addition, we established an independent prognostic indicator, SRSF7+ MCs risk score (SMRS), which showed a correlation between high SMRS group and poor prognosis.

Conclusion

These findings illuminate the complex interactions within the tumor microenvironment of EC and suggest that targeting specific MCs subtypes, particularly via the EGFR signaling pathway, may present novel therapeutic strategies. This study establishes a comprehensive immune map of EC, offering insights for improved treatment approaches.

Tracking changes in functionality and morphology of repopulated microglia in young and old mice

BACKGROUND

Microglia (MG) are myeloid cells of the central nervous system that support homeostasis and instigate neuroinflammation in pathologies. Single-cell RNA sequencing (scRNA-seq) revealed the functional heterogeneity of MG in mouse brains. Microglia are self-renewing cells and inhibition of colony-stimulating factor 1 receptor (CSF1R) signaling depletes microglia which rapidly repopulate. The functions of repopulated microglia are poorly known.

METHODS

We combined scRNA-seq, bulk RNA-seq, immunofluorescence, and confocal imaging to study the functionalities and morphology of repopulated microglia.

RESULTS

A CSRF1R inhibitor (BLZ-945) depleted microglia within 21 days and a number of microglia was fully restored within 7 days, as confirmed by TMEM119 staining and flow cytometry. ScRNA-seq and computational analyses demonstrate that repopulated microglia originated from preexisting progenitors and reconstituted functional clusters but upregulated inflammatory genes. Percentages of proliferating, immature microglia displaying inflammatory gene expression increased in aging mice. Morphometric analysis of MG cell body and branching revealed a distinct morphology of repopulated MG, particularly in brains of old mice. We demonstrate that with aging some repopulated MG fail to reach the homeostatic phenotype. These differences may contribute to the deterioration of MG protective functions with age.

Single-cell analysis reveals hypoxia-induced immunosuppressive microenvironment in intrahepatic cholangiocarcinoma

The role of hypoxia in the tumor microenvironment of intrahepatic cholangiocarcinoma (iCCA) remains unclear. Here, we generated a comprehensive atlas of the entire tumor microenvironment and delineated the multifaceted cell-cell interactions to decipher hypoxia-induced pro-tumor immune suppression. We discovered hypoxia is significantly associated with iCCA progression via the activation of HIF1A expression. Moreover, hypoxia-dependent PPARγ-mediated fatty acid oxidation in APOE+ TAMs promoted M2 macrophage polarization by activating the HIF1A-PPARG-CD36 axis. These polarized APOE+ TAMs recruited Treg cell infiltration via the CCL3-CCR5 pair to form an immunosuppressive microenvironment. APOE+ TAMs tended to co-localize spatially with Treg cells in the malignant tissue based on spatial transcriptome data and immunofluorescence analysis results. We identified tumor-reactive CXCL13+ CD8-PreTex with specific high expression of ENTPD1 and ITGAE, which acted as precursors of CD8-Tex and had higher cytotoxicity, lower exhaustion, and more vigorous proliferation. Consequently, CXCL13+ CD8-PreTex functioned as a positive regulator of antitumor immunity by expressing the pro-inflammatory cytokines IFNG and TNF, associated with a better survival outcome. Our study reveals the mechanisms involved in hypoxia-induced immunosuppression and suggests that targeting precursor-exhausted CXCL13+CD8+ T cells might provide a pratical immunotherapeutic approach.

Inducible pluripotent stem cells to study human mast cell trajectories

Mast cells (MCs) are derived from CD34+ hematopoietic progenitors, consist of different subtypes, and are involved in several inflammatory conditions. However, our understanding of human MC developmental trajectories and subtypes has been limited by a scarcity of suitable cellular model systems. Herein, we developed an in vitro model of human MC differentiation from induced pluripotent stem cells (iPSC) to study human MC differentiation trajectories. Flow cytometry characterization of hemopoietic cells derived from the myeloid cells-forming complex (MCFC) revealed an initial increase in Lin- CD34+ hematopoietic progenitors within Weeks 1-3, followed by an increase in CD34- CD45RA- SSClow and SSChigh hematopoietic cells. The Lin- CD34+ hematopoietic progenitors consisted of SSClow CD45RA- CD123± c-Kit+ FcεRI+ populations that were β7-integrinhigh CD203c+ and β7-integrinhigh CD203c- cells consistent with CMPFcεRI+ cells. Flow cytometry and cytologic analyses of the CD34- Lin- (SSClow) population revealed hypogranular cell populations, predominantly characterized by CD45RA- CD123± c-Kit+ FcεRI- β7-integrinlow and CD45RA- CD123± c-Kit- FcεRI+ β7-integrinMid cells. Analyses of hypergranular SSChigh cells identified Lin- CD34- CD45RA- c-Kit+ FcεRI- and Lin- CD34- CD45RA- c-Kit+ FcεRI+ cells. scRNA-seq analysis of the cells harvested at week 4 of the MCFC culture revealed the presence of monocyte and granulocyte progenitors (n = 547 cells, 26.7 %), Erythrocyte / unknown (n = 85, 4.1 %), neutrophils / myelocytes (n = 211 cells, 10.2 %), mast cell progenitor 1 (n = 599, 29.1 %), mast cell progenitor 2 (n = 152, 7.4 %), committed mast cell precursor (n = 113, 5.5 %), and MCs (n = 353, 17.1 %). In silico analyses of the MC precursor and mature MC populations revealed transcriptionally distinct MC precursor subtype and mature MC states (CMA1+ and CMA1- subtypes). Culturing MC precursor populations in MC maturation media (mast cell media II) led to homogenous mature MC populations as evidenced by high expression of high-affinity IgE receptor, metachromatic granules, presence of MC granule proteins (Tryptase and Chymase) and activation following substance P stimulation and FcεRI crosslinking. This human iPSC-based approach generates MC precursors and phenotypically mature and functional MC populations. This system will be a useful model to generate human MC populations and broaden our understanding of MC biology and transcriptional regulation of MC differentiation trajectories.

Wnt-deficient and hypoxic environment orchestrates squamous reprogramming of human pancreatic ductal adenocarcinoma

Human pancreatic cancer is characterized by the molecular diversity encompassing native duct-like and squamous cell-like identities, but mechanisms underlying squamous transdifferentiation have remained elusive. To comprehensively capture the molecular diversity of human pancreatic cancer, we here profiled 65 patient-derived pancreatic cancer organoid lines, including six adenosquamous carcinoma lines. H3K27me3-mediated erasure of the ductal lineage specifiers and hijacking of the TP63-driven squamous-cell programme drove squamous-cell commitment, providing survival benefit in a Wnt-deficient environment and hypoxic conditions. Gene engineering of normal pancreatic duct organoids revealed that GATA6 loss and a Wnt-deficient environment, in concert with genetic or hypoxia-mediated inactivation of KDM6A, facilitate squamous reprogramming, which in turn enhances environmental fitness. EZH2 inhibition counterbalanced the epigenetic bias and curbed the growth of adenosquamous cancer organoids. Our results demonstrate how an adversarial microenvironment dictates the molecular and histological evolution of human pancreatic cancer and provide insights into the principles and significance of lineage conversion in human cancer.

Trajectory-centric framework TrajAtlas reveals multi-scale differentiation heterogeneity among cells, genes, and gene modules in osteogenesis

Osteoblasts, the key cells responsible for bone formation and the maintenance of skeletal integrity, originate from a diverse array of progenitor cells. However, the mechanisms underlying osteoblast differentiation from these multiple osteoprogenitors remain poorly understood. To address this knowledge gap, we developed a comprehensive framework to investigate osteoblast differentiation at multiple scales, encompassing cells, genes, and gene modules. We constructed a reference atlas focused on differentiation, which incorporates various osteoprogenitors and provides a seven-level cellular taxonomy. To reconstruct the differentiation process, we developed a model that identifies the transcription factors and pathways involved in differentiation from different osteoprogenitors. Acknowledging that covariates such as age and tissue type can influence differentiation, we created an algorithm to detect differentially expressed genes throughout the differentiation process. Additionally, we implemented methods to identify conserved pseudotemporal gene modules across multiple samples. Overall, our framework systematically addresses the heterogeneity observed during osteoblast differentiation from diverse sources, offering novel insights into the complexities of bone formation and serving as a valuable resource for understanding osteogenesis.

A single-cell transcriptomic dataset of pluripotent stem cell-derived astrocytes via NFIB/SOX9 overexpression

Astrocytes, the predominant glial cells in the central nervous system, play essential roles in maintaining brain function. Reprogramming induced pluripotent stem cells (iPSCs) to become astrocytes through overexpression of the transcription factors, NFIB and SOX9, is a rapid and efficient approach for studying human neurological diseases and identifying therapeutic targets. However, the precise differentiation path and molecular signatures of induced astrocytes remain incompletely understood. Accordingly, we performed single-cell RNA sequencing analysis on 64,736 cells to establish a comprehensive atlas of NFIB/SOX9-directed astrocyte differentiation from human iPSCs. Our dataset provides detailed information about the path of astrocyte differentiation, highlighting the stepwise molecular changes that occur throughout the differentiation process. This dataset serves as a valuable reference for dissecting uncharacterized transcriptomic features of NFIB/SOX9-induced astrocytes and investigating lineage progression during astrocyte differentiation. Moreover, these findings pave the way for future studies on neurological diseases using the NFIB/SOX9-induced astrocyte model.

Joint trajectory inference for single-cell genomics using deep learning with a mixture prior

Trajectory inference methods are essential for analyzing the developmental paths of cells in single-cell sequencing datasets. It provides insights into cellular differentiation, transitions, and lineage hierarchies, helping unravel the dynamic processes underlying development and disease progression. However, many existing tools lack a coherent statistical model and reliable uncertainty quantification, limiting their utility and robustness. In this paper, we introduce VITAE (Variational Inference for Trajectory by AutoEncoder), a statistical approach that integrates a latent hierarchical mixture model with variational autoencoders to infer trajectories. The statistical hierarchical model enhances the interpretability of our framework, while the posterior approximations generated by our variational autoencoder ensure computational efficiency and provide uncertainty quantification of cell projections along trajectories. Specifically, VITAE enables simultaneous trajectory inference and data integration, improving the accuracy of learning a joint trajectory structure in the presence of biological and technical heterogeneity across datasets. We show that VITAE outperforms other state-of-the-art trajectory inference methods on both real and synthetic data under various trajectory topologies. Furthermore, we apply VITAE to jointly analyze three distinct single-cell RNA sequencing datasets of the mouse neocortex, unveiling comprehensive developmental lineages of projection neurons. VITAE effectively reduces batch effects within and across datasets and uncovers finer structures that might be overlooked in individual datasets. Additionally, we showcase VITAE's efficacy in integrative analyses of multiomic datasets with continuous cell population structures.

Effector memory-type regulatory T cells display phenotypic and functional instability

Regulatory T cells (Treg cells) hold promise for sustainable therapy of immune disorders. Recent advancements in chimeric antigen receptor development and genome editing aim to enhance the specificity and function of Treg cells. However, impurities and functional instability pose challenges for the development of safe gene-edited Treg cell products. Here, we examined different Treg cell subsets regarding their fate, epigenomic stability, transcriptomes, T cell receptor repertoires, and function ex vivo and after manufacturing. Each Treg cell subset displayed distinct features, including lineage stability, epigenomics, surface markers, T cell receptor diversity, and transcriptomics. Earlier-differentiated memory Treg cell populations, including a hitherto unidentified naïve-like memory Treg cell subset, outperformed late-differentiated effector memory-like Treg cells in regulatory function, proliferative capacity, and epigenomic stability. High yields of stable, functional Treg cell products could be achieved by depleting the small effector memory-like Treg cell subset before manufacturing. Considering Treg cell subset composition appears critical to maintain lineage stability in the final cell product.

Molecular, Cellular, and Developmental Organization of the Mouse Vomeronasal organ at Single Cell Resolution

We have generated single cell transcriptomic atlases of vomeronasal organs (VNO) from juvenile and adult mice. Combined with spatial molecular imaging, we uncover a distinct, previously unidentified class of cells that express the vomeronasal receptors and a population of canonical olfactory sensory neurons in the VNO. High resolution trajectory and cluster analyses reveal the lineage relationship, spatial distribution of cell types, and a putative cascade of molecular events that specify the V1r, V2r, and OR lineages from a common stem cell population. The expression of vomeronasal and olfactory receptors follow power law distributions, but there is high variability in average expression levels between individual receptor and cell types. Substantial co-expression is found between receptors across clades, from different classes, and between olfactory and vomeronasal receptors, with nearly half from pairs located on the same chromosome. Interestingly, the expression of V2r, but not V1r, genes is associated with various transcription factors, suggesting distinct mechanisms of receptor choice associated with the two cell types. We identify association between transcription factors, surface axon guidance molecules, and individual VRs, thereby uncovering a molecular code that guides the specification of the vomeronasal circuitry. Our study provides a wealth of data on the development and organization of the accessory olfactory system at both cellular and molecular levels to enable a deeper understanding of vomeronasal system function.

LMNA-Related Dilated Cardiomyopathy: Single-Cell Transcriptomics during Patient-Derived iPSC Differentiation Support Cell Type and Lineage-Specific Dysregulation of Gene Expression and Development for Cardiomyocytes and Epicardium-Derived Cells with Lamin A/C Haploinsufficiency

LMNA-related dilated cardiomyopathy (DCM) is an autosomal-dominant genetic condition with cardiomyocyte and conduction system dysfunction often resulting in heart failure or sudden death. The condition is caused by mutation in the Lamin A/C (LMNA) gene encoding Type-A nuclear lamin proteins involved in nuclear integrity, epigenetic regulation of gene expression, and differentiation. The molecular mechanisms of the disease are not completely understood, and there are no definitive treatments to reverse progression or prevent mortality. We investigated possible mechanisms of LMNA-related DCM using induced pluripotent stem cells derived from a family with a heterozygous LMNA c.357-2A>G splice-site mutation. We differentiated one LMNA-mutant iPSC line derived from an affected female (Patient) and two non-mutant iPSC lines derived from her unaffected sister (Control) and conducted single-cell RNA sequencing for 12 samples (four from Patients and eight from Controls) across seven time points: Day 0, 2, 4, 9, 16, 19, and 30. Our bioinformatics workflow identified 125,554 cells in raw data and 110,521 (88%) high-quality cells in sequentially processed data. Unsupervised clustering, cell annotation, and trajectory inference found complex heterogeneity: ten main cell types; many possible subtypes; and lineage bifurcation for cardiac progenitors to cardiomyocytes (CMs) and epicardium-derived cells (EPDCs). Data integration and comparative analyses of Patient and Control cells found cell type and lineage-specific differentially expressed genes (DEGs) with enrichment, supporting pathway dysregulation. Top DEGs and enriched pathways included 10 ZNF genes and RNA polymerase II transcription in pluripotent cells (PP); BMP4 and TGF Beta/BMP signaling, sarcomere gene subsets and cardiogenesis, CDH2 and EMT in CMs; LMNA and epigenetic regulation, as well as DDIT4 and mTORC1 signaling in EPDCs. Top DEGs also included XIST and other X-linked genes, six imprinted genes (SNRPN, PWAR6, NDN, PEG10, MEG3, MEG8), and enriched gene sets related to metabolism, proliferation, and homeostasis. We confirmed Lamin A/C haploinsufficiency by allelic expression and Western blot. Our complex Patient-derived iPSC model for Lamin A/C haploinsufficiency in PP, CM, and EPDC provided support for dysregulation of genes and pathways, many previously associated with Lamin A/C defects, such as epigenetic gene expression, signaling, and differentiation. Our findings support disruption of epigenomic developmental programs, as proposed in other LMNA disease models. We recognized other factors influencing epigenetics and differentiation; thus, our approach needs improvement to further investigate this mechanism in an iPSC-derived model.

TrajectoryGeometry suggests cell fate decisions can involve branches rather than bifurcations

Differentiation of multipotential progenitor cells is a key process in the development of any multi-cellular organism and often continues throughout its life. It is often assumed that a bi-potential progenitor develops along a (relatively) straight trajectory until it reaches a decision point where the trajectory bifurcates. At this point one of two directions is chosen, each direction representing the unfolding of a new transcriptional programme. However, we have lacked quantitative means for testing this model. Accordingly, we have developed the R package TrajectoryGeometry. Applying this to published data we find several examples where, rather than bifurcate, developmental pathways branch. That is, the bipotential progenitor develops along a relatively straight trajectory leading to one of its potential fates. A second relatively straight trajectory branches off from this towards the other potential fate. In this sense only cells that branch off to follow the second trajectory make a 'decision'. Our methods give precise descriptions of the genes and cellular pathways involved in these trajectories. We speculate that branching may be the more common behaviour and may have advantages from a control-theoretic viewpoint.

Single-cell multi-omics reveal stage of differentiation and trajectory-dependent immunity-related gene expression patterns in human erythroid cells

The role of Erythroid cells in immune regulation and immunosuppression is one of the emerging topics in modern immunology that still requires further clarification as Erythroid cells from different tissues and different species express different immunoregulatory molecules. In this study, we performed a thorough investigation of human bone marrow Erythroid cells from adult healthy donors and adult acute lymphoblastic leukemia patients using the state-of-the-art single-cell targeted proteomics and transcriptomics via BD Rhapsody and cancer-related gene copy number variation analysis via NanoString Sprint Profiler. We found that human bone marrow Erythroid cells express the ARG1, LGALS1, LGALS3, LGALS9, and C10orf54 (VISTA) immunosuppressive genes, CXCL5, CXCL8, and VEGFA cytokine genes, as well as the genes involved in antimicrobial immunity and MHC Class II antigen presentation. We also found that ARG1 gene expression was restricted to the single erythroid cell cluster that we termed ARG1-positive Orthochromatic erythroblasts and that late Erythroid cells lose S100A9 and gain MZB1 gene expression in case of acute lymphoblastic leukemia. These findings show that steady-state erythropoiesis bone marrow Erythroid cells express myeloid signature genes even without any transdifferentiating stimulus like cancer.

Inferring histology-associated gene expression gradients in spatial transcriptomic studies

Spatially resolved transcriptomics has revolutionized RNA studies by aligning RNA abundance with tissue structure, enabling direct comparisons between histology and gene expression. Traditional approaches to identifying signature genes often involve preliminary data grouping, which can overlook subtle expression patterns in complex tissues. We present Spatial Gradient Screening, an algorithm which facilitates the supervised detection of histology-associated gene expression patterns without prior data grouping. Utilizing spatial transcriptomic data along with single-cell deconvolution from injured mouse cortex, and TCR-seq data from brain tumors, we compare our methodology to standard differential gene expression analysis. Our findings illustrate both the advantages and limitations of cluster-free detection of gene expression, offering more profound insights into the spatial architecture of transcriptomes. The algorithm is embedded in SPATA2, an open-source framework written in R, which provides a comprehensive set of tools for investigating gene expression within tissue.

The Impact of Low Protein Diet on the Molecular and Cellular Development of the Fetal Kidney

Background

Low nephron number has a direct impact on the development of hypertension and chronic kidney disease later in life. While intrauterine growth restriction caused by maternal low protein diet (LPD) is thought to be a significant cause of reduced nephron endowment in impoverished communities, its influence on the cellular and molecular processes which drive nephron formation are poorly understood.

Methods

We conducted a comprehensive characterization of the impact of LPD on kidney development using tomographic and confocal imaging to quantify changes in branching morphogenesis and the cellular and morphological features of nephrogenic niches across development. These analyses were paired with single-cell RNA sequencing to dissect the transcriptional changes that LPD imposes during renal development to affect nephron number.

Results

Single cell analysis at E14.5 and P0 revealed differences in the expression of genes and pathways involved in metabolism, cell cycle, epigenetic regulators and reciprocal inductive signals in most cell types analyzed, yielding imbalances and shifts in cellular energy production and cellular trajectories. In the nephron progenitor cells, LPD impeded cellular commitment and differentiation towards pre-tubular and renal vesicle structures. Confocal microscopy revealed a reduction in the number of pre-tubular aggregates and proliferation in nephron progenitor cells. We also found changes in branching morphogenesis, with a reduction in cell proliferation in the ureteric tips as well as reduced tip and tip parent lengths by optical projection tomography which causes patterning defects.

Conclusions

This unique profiling demonstrates how a fetal programming defect leads to low nephron endowment which is intricately linked to changes in both branching morphogenesis and the commitment of nephron progenitor cells. The commitment of progenitor cells is pivotal for nephron formation and is significantly influenced by nutritional factors, with a low protein diet driving alterations in this program which directly results in a reduced nephron endowment.

Significance Statement

While a mother's diet can negatively impact the number of nephrons in the kidneys of her offspring, the root cellular and molecular drivers of these deficits have not been rigorously explored. In this study we use advanced imaging and gene expression analysis in mouse models to define how a maternal low protein diet, analogous to that of impoverished communities, results in reduced nephron endowment. We find that low protein diet has pleiotropic effects on metabolism and the normal developmental programs of gene expression. These profoundly impact the process of branching morphogenesis necessary to establish niches for nephron generation and change cell behaviors which regulate how and when nephron progenitor cells commit to differentiation.

Associating transcription factors to single-cell trajectories with DREAMIT

Inferring gene regulatory networks from single-cell RNA-sequencing trajectories has been an active area of research yet methods are still needed to identify regulators governing cell transitions. We developed DREAMIT (Dynamic Regulation of Expression Across Modules in Inferred Trajectories) to annotate transcription-factor activity along single-cell trajectory branches, using ensembles of relations to target genes. Using a benchmark representing several different tissues, as well as external validation with ATAC-Seq and Perturb-Seq data on hematopoietic cells, the method was found to have higher tissue-specific sensitivity and specificity over competing approaches.

Multi-omic analysis of Huntington's disease reveals a compensatory astrocyte state

The mechanisms underlying the selective regional vulnerability to neurodegeneration in Huntington's disease (HD) have not been fully defined. To explore the role of astrocytes in this phenomenon, we used single-nucleus and bulk RNAseq, lipidomics, HTT gene CAG repeat-length measurements, and multiplexed immunofluorescence on HD and control post-mortem brains. We identified genes that correlated with CAG repeat length, which were enriched in astrocyte genes, and lipidomic signatures that implicated poly-unsaturated fatty acids in sensitizing neurons to cell death. Because astrocytes play essential roles in lipid metabolism, we explored the heterogeneity of astrocytic states in both protoplasmic and fibrous-like (CD44+) astrocytes. Significantly, one protoplasmic astrocyte state showed high levels of metallothioneins and was correlated with the selective vulnerability of distinct striatal neuronal populations. When modeled in vitro, this state improved the viability of HD-patient-derived spiny projection neurons. Our findings uncover key roles of astrocytic states in protecting against neurodegeneration in HD.

Fasting induces metabolic switches and spatial redistributions of lipid processing and neuronal interactions in tanycytes

The ependyma lining the third ventricle (3V) in the mediobasal hypothalamus plays a crucial role in energy balance and glucose homeostasis. It is characterized by a high functional heterogeneity and plasticity, but the underlying molecular mechanisms governing its features are not fully understood. Here, 5481 hypothalamic ependymocytes were cataloged using FACS-assisted scRNAseq from fed, 12h-fasted, and 24h-fasted adult male mice. With standard clustering analysis, typical ependymal cells and β2-tanycytes appear sharply defined, but other subpopulations, β1- and α-tanycytes, display fuzzy boundaries with few or no specific markers. Pseudospatial approaches, based on the 3V neuroanatomical distribution, enable the identification of specific versus shared tanycyte markers and subgroup-specific versus general tanycyte functions. We show that fasting dynamically shifts gene expression patterns along the 3V, leading to a spatial redistribution of cell type-specific responses. Altogether, we show that changes in energy status induce metabolic and functional switches in tanycyte subpopulations, providing insights into molecular and functional diversity and plasticity within the tanycyte population.

Recruited atypical Ly6G+ macrophages license alveolar regeneration after lung injury

The lung is constantly exposed to airborne pathogens and particles that can cause alveolar damage. Hence, appropriate repair responses are essential for gas exchange and life. Here, we deciphered the spatiotemporal trajectory and function of an atypical population of macrophages after lung injury. Post-influenza A virus (IAV) infection, short-lived monocyte-derived Ly6G-expressing macrophages (Ly6G+ Macs) were recruited to the alveoli of lung perilesional areas. Ly6G+ Macs engulfed immune cells, exhibited a high metabolic potential, and clustered with alveolar type 2 epithelial cells (AT2s) in zones of active epithelial regeneration. Ly6G+ Macs were partially dependent on granulocyte-macrophage colony-stimulating factor and interleukin-4 receptor signaling and were essential for AT2-dependent alveolar regeneration. Similar macrophages were recruited in other models of injury and in the airspaces of lungs from patients with suspected pneumonia. This study identifies perilesional alveolar Ly6G+ Macs as a spatially restricted, short-lived macrophage subset promoting epithelial regeneration postinjury, thus representing an attractive therapeutic target for treating lung damage.

Revealing the crucial roles of suppressive immune microenvironment in cardiac myxoma progression

Cardiac myxoma is a commonly encountered tumor within the heart that has the potential to be life-threatening. However, the cellular composition of this condition is still not well understood. To fill this gap, we analyzed 75,641 cells from cardiac myxoma tissues based on single-cell sequencing. We defined a population of myxoma cells, which exhibited a resemblance to fibroblasts, yet they were distinguished by an increased expression of phosphodiesterases and genes associated with cell proliferation, differentiation, and adhesion. The clinical relevance of the cell populations indicated a higher proportion of myxoma cells and M2-like macrophage infiltration, along with their enhanced spatial interaction, were found to significantly contribute to the occurrence of embolism. The immune cells surrounding the myxoma exhibit inhibitory characteristics, with impaired function of T cells characterized by the expression of GZMK and TOX, along with a substantial infiltration of tumor-promoting macrophages expressed growth factors such as PDGFC. Furthermore, in vitro co-culture experiments showed that macrophages promoted the growth of myxoma cells significantly. In summary, this study presents a comprehensive single-cell atlas of cardiac myxoma, highlighting the heterogeneity of myxoma cells and their collaborative impact on immune cells. These findings shed light on the complex pathobiology of cardiac myxoma and present potential targets for intervention.