Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics

Selective utilization of glucose metabolism guides mammalian gastrulation

The prevailing dogma for morphological patterning in developing organisms argues that the combined inputs of transcription factor networks and signalling morphogens alone generate spatially and temporally distinct expression patterns. However, metabolism has also emerged as a critical developmental regulator1-10, independent of its functions in energy production and growth. The mechanistic role of nutrient utilization in instructing cellular programmes to shape the in vivo developing mammalian embryo remains unknown. Here we reveal two spatially resolved, cell-type- and stage-specific waves of glucose metabolism during mammalian gastrulation by using single-cell-resolution quantitative imaging of developing mouse embryos, stem cell models and embryo-derived tissue explants. We identify that the first spatiotemporal wave of glucose metabolism occurs through the hexosamine biosynthetic pathway to drive fate acquisition in the epiblast, and the second wave uses glycolysis to guide mesoderm migration and lateral expansion. Furthermore, we demonstrate that glucose exerts its influence on these developmental processes through cellular signalling pathways, with distinct mechanisms connecting glucose with the ERK activity in each wave. Our findings underscore that-in synergy with genetic mechanisms and morphogenic gradients-compartmentalized cellular metabolism is integral in guiding cell fate and specialized functions during development. This study challenges the view of the generic and housekeeping nature of cellular metabolism, offering valuable insights into its roles in various developmental contexts.

Age-associated clonal B cells drive B cell lymphoma in mice

Although cancer is an age-related disease, how the processes of aging contribute to cancer progression is not well understood. In this study, we uncovered how mouse B cell lymphoma develops as a consequence of a naturally aged system. We show here that this malignancy is associated with an age-associated clonal B cell (ACBC) population that likely originates from age-associated B cells. Driven by c-Myc activation, promoter hypermethylation and somatic mutations, IgM+ ACBCs clonally expand independently of germinal centers and show increased biological age. ACBCs become self-sufficient and support malignancy when transferred into young recipients. Inhibition of mTOR or c-Myc in old mice attenuates pre-malignant changes in B cells during aging. Although the etiology of mouse and human B cell lymphomas is considered distinct, epigenetic changes in transformed mouse B cells are enriched for changes observed in human B cell lymphomas. Together, our findings characterize the spontaneous progression of cancer during aging through both cell-intrinsic and microenvironmental changes and suggest interventions for its prevention.

Combined expansion and STED microscopy reveals altered fingerprints of postsynaptic nanostructure across brain regions in ASD-related SHANK3-deficiency

Synaptic dysfunction is a key feature of SHANK-associated disorders such as autism spectrum disorder, schizophrenia, and Phelan-McDermid syndrome. Since detailed knowledge of their effect on synaptic nanostructure remains limited, we aimed to investigate such alterations in ex11|SH3 SHANK3-KO mice combining expansion and STED microscopy. This enabled high-resolution imaging of mosaic-like arrangements formed by synaptic proteins in both human and murine brain tissue. We found distinct shape-profiles as fingerprints of the murine postsynaptic scaffold across brain regions and genotypes, as well as alterations in the spatial and molecular organization of subsynaptic domains under SHANK3-deficient conditions. These results provide insights into synaptic nanostructure in situ and advance our understanding of molecular mechanisms underlying synaptic dysfunction in neuropsychiatric disorders.

Septo-dentate gyrus cholinergic circuits modulate function and morphogenesis of adult neural stem cells through granule cell intermediaries

Cholinergic neurons in the basal forebrain play a crucial role in regulating adult hippocampal neurogenesis (AHN). However, the circuit and molecular mechanisms underlying cholinergic modulation of AHN, especially the initial stages of this process related to the generation of newborn progeny from quiescent radial neural stem cells (rNSCs), remain unclear. Here, we report that stimulation of the cholinergic circuits projected from the diagonal band of Broca (DB) to the dentate gyrus (DG) neurogenic niche promotes proliferation and morphological development of rNSCs, resulting in increased neural stem/progenitor pool and rNSCs with longer radial processes and larger busy heads. Interestingly, DG granule cells (GCs) are required for DB-DG cholinergic circuit-dependent modulation of proliferation and morphogenesis of rNSCs. Furthermore, single-nucleus RNA sequencing of DG reveals cell type-specific transcriptional changes in response to cholinergic circuit stimulation, with GCs (among all the DG niche cells) exhibiting the most extensive transcriptional changes. Our findings shed light on how the DB-DG cholinergic circuits orchestrate the key niche components to support neurogenic function and morphogenesis of rNSCs at the circuit and molecular levels.

Next-generation spatial transcriptomics: unleashing the power to gear up translational oncology

The growing advances in spatial transcriptomics (ST) stand as the new frontier bringing unprecedented influences in the realm of translational oncology. This has triggered systemic experimental design, analytical scope, and depth alongside with thorough bioinformatics approaches being constantly developed in the last few years. However, harnessing the power of spatial biology and streamlining an array of ST tools to achieve designated research goals are fundamental and require real-world experiences. We present a systemic review by updating the technical scope of ST across different principal basis in a timeline manner hinting on the generally adopted ST techniques used within the community. We also review the current progress of bioinformatic tools and propose in a pipelined workflow with a toolbox available for ST data exploration. With particular interests in tumor microenvironment where ST is being broadly utilized, we summarize the up-to-date progress made via ST-based technologies by narrating studies categorized into either mechanistic elucidation or biomarker profiling (translational oncology) across multiple cancer types and their ways of deploying the research through ST. This updated review offers as a guidance with forward-looking viewpoints endorsed by many high-resolution ST tools being utilized to disentangle biological questions that may lead to clinical significance in the future.

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.

DNA methylation controls stemness of astrocytes in health and ischaemia

Astrocytes are the most abundant cell type in the mammalian brain and provide structural and metabolic support to neurons, regulate synapses and become reactive after injury and disease. However, a small subset of astrocytes settles in specialized areas of the adult brain where these astrocytes instead actively generate differentiated neuronal and glial progeny and are therefore referred to as neural stem cells1-3. Common parenchymal astrocytes and quiescent neural stem cells share similar transcriptomes despite their very distinct functions4-6. Thus, how stem cell activity is molecularly encoded remains unknown. Here we examine the transcriptome, chromatin accessibility and methylome of neural stem cells and their progeny, and of astrocytes from the striatum and cortex in the healthy and ischaemic adult mouse brain. We identify distinct methylation profiles associated with either astrocyte or stem cell function. Stem cell function is mediated by methylation of astrocyte genes and demethylation of stem cell genes that are expressed later. Ischaemic injury to the brain induces gain of stemness in striatal astrocytes7. We show that this response involves reprogramming the astrocyte methylome to a stem cell methylome and is absent if the de novo methyltransferase DNMT3A is missing. Overall, we unveil DNA methylation as a promising target for regenerative medicine.

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.

Alternative Balance between Transcriptional and Epigenetic Regulation during Developmental Proliferation of Human Cranial Neural Crest Cells

Cranial neural crest cells are implicated in multiple transcriptional events at the different stages of differentiation during development. The alteration of some transcription factors expressed during neural crest development, like PAX7, could be implicated in the etiology of face malformation in murine models. Epigenetic regulation has been shown to be an important mechanistic actor in the control of timing and the level of gene expression at different stages of neural crest development. During this work, we investigated the interconnection between epigenetics and transcription factors across a diversity of human development cranial neural crest cells. Across a diversity of neural cells from human developing cranial tissues, in accordance with their proliferation stage, an alternative balance of regulation between transcription factors and epigenetic factors was identified.

The pericardium forms as a distinct structure during heart formation

The heart integrates diverse cell lineages into a functional unit, including the pericardium, a mesothelial sac that supports heart movement, homeostasis, and immune responses. However, despite its critical roles, the developmental origins of the pericardium remain uncertain due to disparate models. Here, using live imaging, lineage tracking, and single-cell transcriptomics in zebrafish, we find the pericardium forms within the lateral plate mesoderm from dedicated anterior mesothelial progenitors and distinct from the classic heart field. Imaging of transgenic reporters in zebrafish documents lateral plate mesoderm cells that emerge lateral of the classic heart field and among a continuous mesothelial progenitor field. Single-cell transcriptomics and trajectories of hand2-expressing lateral plate mesoderm reveal distinct populations of mesothelial and cardiac precursors, including pericardial precursors that are distinct from the cardiomyocyte lineage. The mesothelial gene expression signature is conserved in mammals and carries over to post-natal development. Light sheet-based live-imaging and machine learning-supported cell tracking documents that during heart tube formation, pericardial precursors that reside at the anterior edge of the heart field migrate anteriorly and medially before fusing, enclosing the embryonic heart to form a single pericardial cavity. Pericardium formation proceeds even upon genetic disruption of heart tube formation, uncoupling the two structures. Canonical Wnt/β-catenin signaling modulates pericardial cell number, resulting in a stretched pericardial epithelium with reduced cell number upon canonical Wnt inhibition. We connect the pathological expression of secreted Wnt antagonists of the SFRP family found in pediatric dilated cardiomyopathy to increased pericardial stiffness: sFRP1 in the presence of increased catecholamines causes cardiomyocyte stiffness in neonatal rats as measured by atomic force microscopy. Altogether, our data integrate pericardium formation as an independent process into heart morphogenesis and connect disrupted pericardial tissue properties such as pericardial stiffness to pediatric cardiomyopathies.

Lipid Peroxidation and Type I Interferon Coupling Fuels Pathogenic Macrophage Activation Causing Tuberculosis Susceptibility

A quarter of human population is infected with Mycobacterium tuberculosis, but less than 10% of those infected develop pulmonary TB. We developed a genetically defined sst1-susceptible mouse model that uniquely reproduces a defining feature of human TB: the development of necrotic lung granulomas and determined that the sst1-susceptible phenotype was driven by the aberrant macrophage activation. This study demonstrates that the aberrant response of the sst1-susceptible macrophages to prolonged stimulation with TNF is primarily driven by conflicting Myc and antioxidant response pathways leading to a coordinated failure 1) to properly sequester intracellular iron and 2) to activate ferroptosis inhibitor enzymes. Consequently, iron-mediated lipid peroxidation fueled IFNβ superinduction and sustained the Type I Interferon (IFN-I) pathway hyperactivity that locked the sst1-susceptible macrophages in a state of unresolving stress and compromised their resistance to Mtb. The accumulation of the aberrantly activated, stressed, macrophages within granuloma microenvironment led to the local failure of anti-tuberculosis immunity and tissue necrosis. The upregulation of Myc pathway in peripheral blood cells of human TB patients was significantly associated with poor outcomes of TB treatment. Thus, Myc dysregulation in activated macrophages results in an aberrant macrophage activation and represents a novel target for host-directed TB therapies.

Olfactory ensheathing cells from adult female rats are hybrid glia that promote neural repair

Olfactory ensheathing cells (OECs) are unique glial cells found in both central and peripheral nervous systems where they support continuous axonal outgrowth of olfactory sensory neurons to their targets. Previously we reported that following severe spinal cord injury, OECs transplanted near the injury site modify the inhibitory glial scar and facilitate axon regeneration past the scar border and into the lesion. To better understand the mechanisms underlying the reparative properties of OECs, we used single-cell RNA-sequencing of OECs from adult rats to study their gene expression programs. Our analyses revealed five diverse OEC subtypes, each expressing novel marker genes and pathways indicative of progenitor, axonal regeneration, secreted molecules, or microglia-like functions. We found substantial overlap of OEC genes with those of Schwann cells, but also with microglia, astrocytes, and oligodendrocytes. We confirmed established markers on cultured OECs, and localized select top genes of OEC subtypes in olfactory bulb tissue. We also show that OECs secrete Reelin and Connective tissue growth factor, extracellular matrix molecules which are important for neural repair and axonal outgrowth. Our results support that OECs are a unique hybrid glia, some with progenitor characteristics, and that their gene expression patterns indicate functions related to wound healing, injury repair and axonal regeneration.

Deconvolving organogenesis in space and time via spatial transcriptomics in thick tissues

Organ development is guided by a space-time landscape that constraints cell behavior. This landscape is challenging to characterize for the hair follicle - the most abundant mini organ - due to its complex microscopic structure and asynchronous development. We developed 3DEEP, a tissue clearing and spatial transcriptomic strategy for characterizing tissue blocks up to 400 µm in thickness. We captured 371 hair follicles at different stages of organogenesis in 1 mm3 of skin of a 12-hour-old mouse with 6 million transcripts from 81 genes. From this single time point, we deconvoluted follicles by age based on whole-organ molecular pseudotimes to animate a stop-motion 3D atlas of follicle development along its trajectory. We defined molecular stages for hair follicle organogenesis and characterized the order of emergence for its structures, differential signaling dynamics at its top and bottom, morphogen shifts preceding and accompanying structural changes, and series of structural changes leading to the formation of its canal and opening. We further found that hair follicle stem cells and their niche are established and stratified early in organogenesis, before the formation of the hair bulb. Overall, this work demonstrates the power of increased depth of spatial transcriptomics to provide a four-dimensional analysis of organogenesis.

Cell cycle expression heterogeneity predicts degree of differentiation

Methods that predict fate potential or degree of differentiation from transcriptomic data have identified rare progenitor populations and uncovered developmental regulatory mechanisms. However, some state-of-the-art methods are too computationally burdensome for emerging large-scale data and all methods make inaccurate predictions in certain biological systems. We developed a method in R (stemFinder) that predicts single cell differentiation time based on heterogeneity in cell cycle gene expression. Our method is computationally tractable and is as good as or superior to competitors. As part of our benchmarking, we implemented four different performance metrics to assist potential users in selecting the tool that is most apt for their application. Finally, we explore the relationship between differentiation time and cell fate potential by analyzing a lineage tracing dataset with clonally labelled hematopoietic cells, revealing that metrics of differentiation time are correlated with the number of downstream lineages.

COFFEE: consensus single cell-type specific inference for gene regulatory networks

The inference of gene regulatory networks (GRNs) is crucial to understanding the regulatory mechanisms that govern biological processes. GRNs may be represented as edges in a graph, and hence, it have been inferred computationally for scRNA-seq data. A wisdom of crowds approach to integrate edges from several GRNs to create one composite GRN has demonstrated improved performance when compared with individual algorithm implementations on bulk RNA-seq and microarray data. In an effort to extend this approach to scRNA-seq data, we present COFFEE (COnsensus single cell-type speciFic inFerence for gEnE regulatory networks), a Borda voting-based consensus algorithm that integrates information from 10 established GRN inference methods. We conclude that COFFEE has improved performance across synthetic, curated, and experimental datasets when compared with baseline methods. Additionally, we show that a modified version of COFFEE can be leveraged to improve performance on newer cell-type specific GRN inference methods. Overall, our results demonstrate that consensus-based methods with pertinent modifications continue to be valuable for GRN inference at the single cell level. While COFFEE is benchmarked on 10 algorithms, it is a flexible strategy that can incorporate any set of GRN inference algorithms according to user preference. A Python implementation of COFFEE may be found on GitHub: https://github.com/lodimk2/coffee.

Dynamic regulation of tissue fluidity controls skin repair during wound healing

During wound healing, different pools of stem cells (SCs) contribute to skin repair. However, how SCs become activated and drive the tissue remodeling essential for skin repair is still poorly understood. Here, by developing a mouse model allowing lineage tracing and basal cell lineage ablation, we monitor SC fate and tissue dynamics during regeneration using confocal and intravital imaging. Analysis of basal cell rearrangements shows dynamic transitions from a solid-like homeostatic state to a fluid-like state allowing tissue remodeling during repair, as predicted by a minimal mathematical modeling of the spatiotemporal dynamics and fate behavior of basal cells. The basal cell layer progressively returns to a solid-like state with re-epithelialization. Bulk, single-cell RNA, and epigenetic profiling of SCs, together with functional experiments, uncover a common regenerative state regulated by the EGFR/AP1 axis activated during tissue fluidization that is essential for skin SC activation and tissue repair.

Tick extracellular vesicles undermine epidermal wound healing during hematophagy

Wound healing has been extensively studied through the lens of inflammatory disorders and cancer, but limited attention has been given to hematophagy and arthropod-borne diseases. Hematophagous ectoparasites, including ticks, subvert the wound healing response to maintain prolonged attachment and facilitate blood-feeding. Here, we unveil a strategy by which extracellular vesicles (EVs) ensure blood-feeding and arthropod survival in three medically relevant tick species. We demonstrate through single cell RNA sequencing and murine genetics that wildtype animals infested with EV-deficient Ixodes scapularis display a unique population of keratinocytes with an overrepresentation of pathways connected to wound healing. Tick feeding affected keratinocyte proliferation in a density-dependent manner, which relied on EVs and dendritic epidermal T cells (DETCs). This occurrence was linked to phosphoinositide 3-kinase activity, keratinocyte growth factor (KGF) and transforming growth factor β (TGF-β) levels. Collectively, we uncovered a strategy employed by a blood-feeding arthropod that impairs the integrity of the epithelial barrier, contributing to ectoparasite fitness.

From single-cell to spatial transcriptomics: decoding the glioma stem cell niche and its clinical implications

Background

Gliomas are aggressive brain tumors associated with a poor prognosis. Cancer stem cells (CSCs) play a significant role in tumor recurrence and resistance to therapy. This study aimed to identify and characterize glioma stem cells (GSCs), analyze their interactions with various cell types, and develop a prognostic signature.

Methods

Single-cell RNA sequencing data from 44 primary glioma samples were analyzed to identify GSC populations. Spatial transcriptomics and gene regulatory network analyses were performed to investigate GSC localization and transcription factor activity. CellChat analysis was conducted to infer cell-cell communication patterns. A GSC signature (GSCS) was developed using machine learning algorithms applied to bulk RNA sequencing data from multiple cohorts. In vitro and in vivo experiments were conducted to validate the role of TUBA1C, a key gene within the signature.

Results

A distinct GSC population was identified, characterized by high proliferative potential and an enrichment of E2F1, E2F2, E2F7, and BRCA1 regulons. GSCs exhibited spatial proximity to myeloid-derived suppressor cells (MDSCs). CellChat analysis revealed an active MIF signaling pathway between GSCs and MDSCs. A 26-gene GSCS demonstrated superior performance compared to existing prognostic models. Knockdown of TUBA1C significantly inhibited glioma cell migration, and invasion in vitro, and reduced tumor growth in vivo.

Conclusion

This study offers a comprehensive characterization of GSCs and their interactions with MDSCs, while presenting a robust GSCS. The findings offer new insights into glioma biology and identify potential therapeutic targets, particularly TUBA1C, aimed at improving patient outcomes.

High ratio of resident to exhausted CD4 + T cells predicts favorable prognosis and potentially better immunotherapeutic efficacy in hepatocellular carcinoma

BACKGROUND

Tumor-infiltrating lymphocytes (TILs) are significantly implicated in regulating the tumor immune microenvironment (TIME) and immunotherapeutic response. However, little is known about the impact of the resident and exhausted status of TILs in hepatocellular carcinoma (HCC).

METHODS

Single-cell RNA sequencing data was applied to discover resident and exhausted signatures of TILs. Survival outcomes, biological function, immune infiltration, genomic variation, immunotherapeutic efficacy, and sorafenib response were further explored the clinical significance and molecular association of TILs in HCC. Moreover, a candidate gene with predictive capability for the dismal subtype was identified through univariate Cox regression analysis, survival analysis, and the BEST website.

RESULTS

Single-cell analysis revealed that CD8 + T, CD4 + T, and NK cells were strongly associated with resident and exhausted patterns. Specific resident and exhausted signatures for each subpopulation were extracted in HCC. Further multivariate Cox analysis revealed that the ratio of resident to exhausted CD4 + T cells in TIME was an independent prognostic factor. After incorporating tumor purity with the ratio of resident to exhausted CD4 + T cells, we stratified HCC patients into three subtypes and found that (i) CD4 residencyhighexhaustionlow subtype was endowed with favorable prognosis, immune activation, and sensitivity to immunotherapy; (ii) CD4 exhaustionhighresidencylow subtype was characterized by genome instability and sensitivity to sorafenib; (iii) Immune-desert subtype was associated with malignant-related pathways and poor prognosis. Furthermore, spindle assembly abnormal protein 6 homolog (SASS6) was identified as a key gene, which accurately predicted the immune-desert subtype. Prognostic analysis as well as in vitro and in vivo experiments further demonstrated that SASS6 was closely associated with tumor prognosis, proliferation, and migration.

CONCLUSIONS

The ratio of resident to exhausted CD4 + T cells shows promise as a potential biomarker for HCC prognosis and immunotherapy response and SASS6 may serve as a biomarker and therapeutic target for prognostic assessment of HCC.

Optimal transport reveals dynamic gene regulatory networks via gene velocity estimation

Inferring gene regulatory networks from gene expression data is an important and challenging problem in the biology community. We propose OTVelo, a methodology that takes time-stamped single-cell gene expression data as input and predicts gene regulation across two time points. It is known that the rate of change of gene expression, which we will refer to as gene velocity, provides crucial information that enhances such inference; however, this information is not always available due to the limitations in sequencing depth. Our algorithm overcomes this limitation by estimating gene velocities using optimal transport. We then infer gene regulation using time-lagged correlation and Granger causality via regularized linear regression. Instead of providing an aggregated network across all time points, our method uncovers the underlying dynamical mechanism across time points. We validate our algorithm on 13 simulated datasets with both synthetic and curated networks and demonstrate its efficacy on 4 experimental data sets.

CD4+ Effective Memory T Cell Markers GBP2 and LAG3 Are Risk Factors for PTB and COVID-19 Infection: A Study Integrating Single-Cell Expression Quantitative Trait Locus and Mendelian Randomization Analyses

Observational studies indicate that variations in peripheral blood mononuclear cell (PBMC) subsets are associated with an increased risk of pulmonary tuberculosis (PTB) and coronavirus disease 2019 (COVID-19), but causal validation is lacking. Here, we combined single-cell expression quantitative trait locus (sc-eQTL) and two-sample mendelian randomization (MR) analyses to elucidate the causal relationship between PBMC subsets and the occurrence of PTB and COVID-19 and verified by RT-qPCR. We observed an increase in the CD4+ Effective Memory T Cell (CD4+ TEM) cluster in both PTB and COVID-19 patients according to the single-cell transcriptional landscape of PBMC. Through MR analysis using an inverse variance weighted (IVW) method, we found strong evidence of positive correlations between CD4+ TEM cell markers (GBP2, TRAV1-2, and ODF2L) and PTB, and between markers (LAG3 and SLFN5) and COVID-19, especially highlighted by lead eQTL-SNPs of GBP2 (rs2256752, p = 4.76321 × 10-15) and LAG3 (rs67706382, p = 6.16× 10-16). Similar results were observed in validation sets, and no pleiotropy was detected in sensitivity analyses including weighted median (WM), MR-Egger, MR-pleiotropy residual sum and outlier, and leave-one-out analyses (all p > 0.05). We visualized the colocalization of marker-eQTLs and markers of PTB and COVID-19 genome-wide association study (GWAS) associations. Based on CellChat analyses, monocytes communicated predominantly with CD4+ TEM cells positively expressing PTB markers (GBP2, TRAV1-2, and ODF2L) and COVID-19 markers (LAG3 and SLFN5) in both PTB and COVID-19. Our data suggest a causal effect between two key CD4+ TEM cell markers (GBP2 and LAG3) and the risk for PTB and COVID-19 infection. Our findings provide novel insights into the biological mechanism for PTB and COVID-19 infection, but future single-cell studies are necessary to further enhance understanding of this find.

Single-Cell RNA-Sequencing Identifies Bone Marrow-Derived Progenitor Cells as a Main Source of Extracellular Matrix-Producing Cells Across Multiple Organ-Based Fibrotic Diseases

Fibrosis is characterized by the aberrant deposition of extracellular matrix (ECM) due to dysregulated tissue repair responses, imposing a significant global burden on fibrosis-related diseases. Although alpha-smooth muscle actin (α-SMA/ACTA2)-expressing myofibroblasts are considered as key player in fibrogenesis, the origin of ECM-producing cells remains controversial. To address this issue, we integrated and analyzed large-scale single-cell transcriptomic datasets from patients with distinct fibrotic diseases involving the heart, lung, liver, or kidney. Unexpectedly, not all ACTA2-expressing cells were ECM-producing cells identified by expressing collagen genes; instead, the majority of ECM-producing cells were myofibroblasts and fibroblasts derived from circulating bone marrow precursor, and to a lesser extent from local pericytes and vascular smooth cells in all fibrotic diseases. This was confirmed in sex-mismatched kidney transplants by the discovery that ECM-producing cells originated from recipient, not donor, bone marrow-derived progenitor cells (BMPCs). Moreover, these BMPCs-derived ECM-producing cells exhibited a proinflammatory phenotype. Thus, bone marrow-derived proinflammatory and profibrotic fibroblasts/myofibroblasts with stem cell properties serve as a major source of ECM-producing cells and may play a driving role in tissue fibrosis across a wide range of human fibrotic diseases. Targeting these ECM-producing cells may provide a novel therapy for diseases with fibrosis.

CD8+ T-cell differentiation and dysfunction inform treatment response in acute myeloid leukemia

ABSTRACT

The interplay between T-cell states of differentiation, dysfunction, and treatment response in acute myeloid leukemia (AML) remains unclear. Here, we leveraged a multimodal approach encompassing high-dimensional flow cytometry and single-cell transcriptomics and found that early memory CD8+ T cells are associated with therapy response and exhibit a bifurcation into 2 distinct terminal end states. One state is enriched for markers of activation, whereas the other expresses natural killer (NK)-like and senescence markers. The skewed clonal differentiation trajectory toward CD8+ senescence was also a hallmark indicative of therapy resistance. We validated these findings by generating an AML CD8+ single-cell atlas integrating our data and other independent data sets. Finally, our analysis revealed that an imbalance between CD8+ early memory and senescent-like cells is linked to AML treatment refractoriness and poor survival. Our study provides crucial insights into the dynamics of CD8+ T-cell differentiation and advances our understanding of CD8+ T-cell dysfunction in AML.

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

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

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.

Setdb1 protects genome integrity in murine muscle stem cells to allow for regenerative myogenesis and inflammation

The histone H3 lysine 9 methyltransferase SETDB1 controls transcriptional repression to direct stem cell fate. Here, we show that Setdb1 expression by adult muscle stem cells (MuSCs) is required for skeletal muscle regeneration. We find that SETDB1 represses the expression of endogenous retroviruses (ERVs) in MuSCs. ERV de-repression in Setdb1-null MuSCs prevents their amplification following exit from quiescence and promotes cell death. Multi-omics profiling shows that chromatin decompaction at ERV loci activates the DNA-sensing cGAS-STING pathway, entailing cytokine expression by Setdb1-null MuSCs. This is followed by aberrant infiltration of inflammatory cells, including pathological macrophages. The ensuing histiocytosis is accompanied by myofiber necrosis, which, in addition to progressive MuSCs depletion, completely abolishes tissue repair. In contrast, loss of Setdb1 in fibro-adipogenic progenitors (FAPs) does not impact immune cells. In conclusion, genome maintenance by SETDB1 in an adult somatic stem cell is necessary for both its regenerative potential and adequate reparative inflammation.

Cell identity revealed by precise cell cycle state mapping links data modalities

Several methods for cell cycle inference from sequencing data exist and are widely adopted. In contrast, methods for classification of cell cycle state from imaging data are scarce. We have for the first time integrated sequencing and imaging derived cell cycle pseudo-times for assigning 449 imaged cells to 693 sequenced cells at an average resolution of 3.4 and 2.4 cells for sequencing and imaging data respectively. Data integration revealed thousands of pathways and organelle features that are correlated with each other, including several previously known interactions and novel associations. The ability to assign the transcriptome state of a profiled cell to its closest living relative, which is still actively growing and expanding opens the door for genotype-phenotype mapping at single cell resolution forward in time.

A primer on single-cell RNA-seq analysis using dendritic cells as a case study

Recent advances in single-cell (sc) transcriptomics have revolutionized our understanding of dendritic cells (DCs), pivotal players of the immune system. ScRNA-sequencing (scRNA-seq) has unraveled a previously unrecognized complexity and heterogeneity of DC subsets, shedding light on their ontogeny and specialized roles. However, navigating the rapid technological progress and computational methods can be daunting for researchers unfamiliar with the field. This review aims to provide immunologists with a comprehensive introduction to sc transcriptomic analysis, offering insights into recent developments in DC biology. Addressing common analytical queries, we guide readers through popular tools and methodologies, supplemented with references to benchmarks and tutorials for in-depth understanding. By examining findings from pioneering studies, we illustrate how computational techniques have expanded our knowledge of DC biology. Through this synthesis, we aim to equip researchers with the necessary tools and knowledge to navigate and leverage scRNA-seq for unraveling the intricacies of DC biology and advancing immunological research.

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.

Convergent evolution of monocyte differentiation in adult skin instructs Langerhans cell identity

Langerhans cells (LCs) are distinct among phagocytes, functioning both as embryo-derived, tissue-resident macrophages in skin innervation and repair and as migrating professional antigen-presenting cells, a function classically assigned to dendritic cells (DCs). Here, we demonstrate that both intrinsic and extrinsic factors imprint this dual identity. Using ablation of embryo-derived LCs in the murine adult skin and tracking differentiation of incoming monocyte-derived replacements, we found intrinsic intraepidermal heterogeneity. We observed that ontogenically distinct monocytes give rise to LCs. Within the epidermis, Jagged-dependent activation of Notch signaling, likely within the hair follicle niche, provided an initial site of LC commitment before metabolic adaptation and survival of monocyte-derived LCs. In the human skin, embryo-derived LCs in newborns retained transcriptional evidence of their macrophage origin, but this was superseded by DC-like immune modules after postnatal expansion. Thus, adaptation to adult skin niches replicates conditioning of LC at birth, permitting repair of the embryo-derived LC network.

Scoping Review: Methods and Applications of Spatial Transcriptomics in Tumor Research

Spatial transcriptomics (ST) examines gene expression within its spatial context on tissue, linking morphology and function. Advances in ST resolution and throughput have led to an increase in scientific interest, notably in cancer research. This scoping study reviews the challenges and practical applications of ST, summarizing current methods, trends, and data analysis techniques for ST in neoplasm research. We analyzed 41 articles published by the end of 2023 alongside public data repositories. The findings indicate cancer biology is an important focus of ST research, with a rising number of studies each year. Visium (10x Genomics, Pleasanton, CA, USA) is the leading ST platform, and SCTransform from Seurat R library is the preferred method for data normalization and integration. Many studies incorporate additional data types like single-cell sequencing and immunohistochemistry. Common ST applications include discovering the composition and function of tumor tissues in the context of their heterogeneity, characterizing the tumor microenvironment, or identifying interactions between cells, including spatial patterns of expression and co-occurrence. However, nearly half of the studies lacked comprehensive data processing protocols, hindering their reproducibility. By recommending greater transparency in sharing analysis methods and adapting single-cell analysis techniques with caution, this review aims to improve the reproducibility and reliability of future studies in cancer research.

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.

Integrated multi-omic high-throughput strategies across-species identified potential key diagnostic, prognostic, and therapeutic targets for atherosclerosis under high glucose conditions

Diabetes is a well-known risk factor for atherosclerosis (AS), but the underlying molecular mechanism remains unknown. The dysregulated immune response is an important reason. High glucose is proven to induce foam cell formation under lipidemia situations in clinical patients. Exploring the potential regulatory programs of accelerated foam cell formation stimulated by high glucose is meaningful. Macrophage-derived foam cells were induced in vitro, and high-throughput sequencing was performed. Coexpression gene modules were constructed using weighted gene co-expression network analysis (WGCNA). Highly related modules were identified. Hub genes were identified by multiple integrative strategies. The potential roles of selected genes were further validated in bulk-RNA and scRNA datasets of human plaques. By transfection of the siRNA, the role of the screened gene during foam cell formation was further explored. Two modules were found to be both positively related to high glucose and ox-LDL. Further enrichment analyses confirmed the association between the brown module and AS. The high correlation between the brown module and macrophages was identified and 4 hub genes (Aldoa, Creg1, Lgmn, and Pkm) were screened. Further validation in external bulk-RNA and scRNA revealed the potential diagnostic and therapeutic value of selected genes. In addition, the survival analysis confirmed the prognostic value of Aldoa while knocking down Aldoa expression alleviated the foam cell formation in vitro. We systematically investigated the synergetic effects of high glucose and ox-LDL during macrophage-derived foam cell formation and identified that ALDOA might be an important diagnostic, prognostic, and therapeutic target in these patients.

Decoding changes in tumor-infiltrating leukocytes through dynamic experimental models and single-cell technologies

The ability to characterize immune cells and explore the molecular interactions that govern their functions has never been greater, fueled in recent years by the revolutionary advance of single-cell analysis platforms. However, precisely how immune cells respond to different stimuli and where differentiation processes and effector functions operate remain incompletely understood. Inferring cellular fate within single-cell transcriptomic analyses is now omnipresent, despite the assumptions typically required in such analyses. Recently developed experimental models support dynamic analyses of the immune response, providing insights into the temporal changes that occur within cells and the tissues in which such transitions occur. Here we will review these approaches and discuss how these can be combined with single-cell technologies to develop a deeper understanding of the immune responses that should support the development of better therapeutic options for patients.

BootCellNet, a resampling-based procedure, promotes unsupervised identification of cell populations via robust inference of gene regulatory networks

Recent advances in measurement technologies, particularly single-cell RNA sequencing (scRNA-seq), have revolutionized our ability to acquire large amounts of omics-level data on cellular states. As measurement techniques evolve, there has been an increasing need for data analysis methodologies, especially those focused on cell-type identification and inference of gene regulatory networks (GRNs). We have developed a new method named BootCellNet, which employs smoothing and resampling to infer GRNs. Using the inferred GRNs, BootCellNet further infers the minimum dominating set (MDS), a set of genes that determines the dynamics of the entire network. We have demonstrated that BootCellNet robustly infers GRNs and their MDSs from scRNA-seq data and facilitates unsupervised identification of cell clusters using scRNA-seq datasets of peripheral blood mononuclear cells and hematopoiesis. It has also identified COVID-19 patient-specific cells and their potential regulatory transcription factors. BootCellNet not only identifies cell types in an unsupervised and explainable way but also provides insights into the characteristics of identified cell types through the inference of GRNs and MDS.

Cellular landscape of adrenocortical carcinoma at single-nuclei resolution

Adrenocortical carcinoma (ACC) is a rare yet devastating tumour of the adrenal gland with a molecular pathology that remains incompletely understood. To gain novel insights into the cellular landscape of ACC, we generated single-nuclei RNA sequencing (snRNA-seq) data sets from twelve ACC tumour samples and analysed these alongside snRNA-seq data sets from normal adrenal glands (NAGs). We find the ACC tumour microenvironment to be relatively devoid of immune cells compared to NAG tissues, consistent with known high tumour purity values for ACC as an immunologically "cold" tumour. Our analysis identifies three separate groups of ACC samples that are characterised by different relative compositions of adrenocortical cell types. These include cell populations that are specifically enriched in the most clinically aggressive and hormonally active tumours, displaying hallmarks of reorganised cell mechanobiology and dysregulated steroidogenesis, respectively. We also identified and validated a population of mitotically active adrenocortical cells that strongly overexpress genes POLQ, DIAPH3 and EZH2 to support tumour expansion alongside an LGR4+ progenitor-like or cell-of-origin candidate for adrenocortical carcinogenesis. Trajectory inference suggests the fate adopted by malignant adrenocortical cells upon differentiation is associated with the copy number or allelic balance state of the imprinted DLK1/MEG3 genomic locus, which we verified by assessing bulk tumour DNA methylation status. In conclusion, our results therefore provide new insights into the clinical and cellular heterogeneity of ACC, revealing how genetic perturbations to healthy adrenocortical renewal and zonation provide a molecular basis for disease pathogenesis.

Heterogeneity effects of bisphenol A and its substitute, fluorene-9-bisphenol, on intestinal homeostasis

Bisphenol A (BPA) and its substitute fluorene-9-bisphenol (BHPF) are used in consumer products; however, their toxic effects on intestinal epithelium remain largely unknown. In this study, we combined intestinal organoids and single-cell RNA sequencing to investigate the impact of BPA and BHPF exposure on intestinal cell composition, differentiation, and function. Both compounds inhibited the growth of small intestinal organoids, with BHPF exhibiting a more potent inhibitory effect. BPA and BHPF did not significantly alter the overall cell type composition; however, they led to different alterations in cell-cell communications. Gene Ontology enrichment analysis showed that BPA and BHPF exposures affected various biological processes, such as glutathione transferase activity, antioxidant activity, and lipid metabolism, in cell type-specific and compound-dependent manners. Trajectory analysis demonstrated that BPA and BHPF altered the differentiation trajectory of the intestinal cells. To further connect the cellular mechanism to the phenotypic impact in vivo, we constructed a mouse model exposed to BPA or BHPF and observed significant alterations in intestinal morphology, including reduced crypt depth and villus length and impaired stem cell proliferation and self-renewal. These results provide novel insights into the cell type-specific effects of BPA and BHPF on the intestinal epithelium and highlight the potential risks of exposure to these compounds. Our findings underscore the importance of evaluating the safety of BPA substitutes and contribute to a better understanding of the effects of environmental chemicals on gut health.

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.

Tools to study neural and glioma stem cell quiescence

Quiescence is a prolonged but reversible state of cell-cycle arrest that is an adaptive feature of most adult stem cell populations. In the brain, quiescence helps to protect adult neural stem cells from stress and supports lifelong neurogenesis. Unfortunately however, entry into a quiescent or a slow-cycling state is also a malignant feature of brain cancer stem cells. In glioblastoma, where the process has been best characterised, quiescent glioma stem cells preferentially survive chemoradiation, and after therapy, reactivate to regrow the tumour and drive recurrence. In this Review, we discuss the in vitro and in vivo models that have been developed for studying neural stem cell quiescence and how these tools may be used to deepen biological understanding and to develop novel therapies targeting quiescent glioma stem cells.

Cancer marker TNFRSF1A: From single‑cell heterogeneity of renal cell carcinoma to functional validation

During the progression of renal cell carcinoma (RCC), tumor growth, metastasis and treatment response heterogeneity are regulated by both the tumor itself and the tumor microenvironment (TME). The aim of the present study was to investigate the role of the TME in RCC and construct a crosstalk network for clear cell RCC (ccRCC). An additional aim was to evaluate whether TNF receptor superfamily member 1A (TNFRSF1A) is a potential therapeutic target for ccRCC. Single-cell data analysis of RCC was performed using the GSE152938 dataset, focusing on key cellular components and their involvement in the ccRCC TME. Additionally, cell-cell communication was analyzed to elucidate the complex network of the ccRCC microenvironment. Analyses of data from The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium databases were performed to further mine the key TNF receptor genes, with a particular focus on the prediction and assessment of the cancer-associated features of TNFRSF1A. In addition, following the silencing of TNFRSF1A using small interfering RNA in the 786-O ccRCC cell line, a number of in vitro experiments were conducted to further investigate the cancer-promoting characteristics of TNFRSF1A. These included 5-ethynyl-2'-deoxyuridine incorporation, Cell Counting Kit-8, colony formation, Transwell, cell cycle and apoptosis assays. The TNF signaling pathway was found to have a critical role in the development of ccRCC. Based on the specific crosstalk identified between TNF and TNFRSF1A, the communication of this signaling pathway within the TME was elucidated. The results of the cellular phenotype experiments indicated that TNFRSF1A promotes the proliferation, migration and invasion of ccRCC cells. Consequently, it is proposed that targeting TNFRSF1A may disrupt tumor progression and serve as a therapeutic strategy. In conclusion, by understanding the TME and identifying significant crosstalk within the TNF signaling pathway, the potential of TNFRSF1A as a therapeutic target is highlighted. This may facilitate an advance in precision medicine and improve the prognosis for patients with RCC.

An intestinal TH17 cell-derived subset can initiate cancer

Approximately 25% of cancers are preceded by chronic inflammation that occurs at the site of tumor development. However, whether this multifactorial oncogenic process, which commonly occurs in the intestines, can be initiated by a specific immune cell population is unclear. Here, we show that an intestinal T cell subset, derived from interleukin-17 (IL-17)-producing helper T (TH17) cells, induces the spontaneous transformation of the intestinal epithelium. This subset produces inflammatory cytokines, and its tumorigenic potential is not dependent on IL-17 production but on the transcription factors KLF6 and T-BET and interferon-γ. The development of this cell type is inhibited by transforming growth factor-β1 (TGFβ1) produced by intestinal epithelial cells. TGFβ signaling acts on the pretumorigenic TH17 cell subset, preventing its progression to the tumorigenic stage by inhibiting KLF6-dependent T-BET expression. This study therefore identifies an intestinal T cell subset initiating cancer.

γδ T-cells in human malignancies: insights from single-cell studies and analytical considerations

γδ T-cells are a rare population of T-cells with both adaptive and innate-like properties. Despite their low prevalence, they have been found to be implicated various human diseases. γδ T-cell infiltration has been associated with improved clinical outcomes in solid cancers, prompting renewed interest in understanding their biology. To date, their biology remains elusive due to their low prevalence. The introduction of high-resolution single-cell sequencing has allowed various groups to characterize key effector subsets in various contexts, as well as begin to elucidate key regulatory mechanisms directing the differentiation and activity of these cells. In this review, we will review some of insights obtained from single-cell studies of γδ T-cells across various malignancies and highlight some important questions that remain unaddressed.

Monocyte Single-Cell Multimodal Profiling in Cardiovascular Disease Risk States

BACKGROUND

Monocytes are a critical innate immune system cell type that serves homeostatic and immunoregulatory functions. They have been identified historically by the cell surface expression of CD14 and CD16. However, recent single-cell studies have revealed that they are much more heterogeneous than previously realized.

METHODS

We utilized cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) and single-cell RNA sequencing to describe the comprehensive transcriptional and phenotypic landscape of 437 126 monocytes.

RESULTS

This high-dimensional multimodal approach identified vast phenotypic diversity and functionally distinct subsets, including IFN-responsive, MHCIIhi (major histocompatibility complex class II), monocyte-platelet aggregates, as well as nonclassical, and several subpopulations of classical monocytes. Using flow cytometry, we validated the existence of MHCII+CD275+ MHCIIhi, CD42b+ monocyte-platelet aggregates, CD16+CD99- nonclassical monocytes, and CD99+ classical monocytes. Each subpopulation exhibited unique characteristics, developmental trajectories, transcriptional regulation, and tissue distribution. In addition, alterations associated with cardiovascular disease risk factors, including race, smoking, and hyperlipidemia were identified. Moreover, the effect of hyperlipidemia was recapitulated in mouse models of elevated cholesterol.

CONCLUSIONS

This integrative and cross-species comparative analysis provides a new perspective on the comparison of alterations in monocytes in pathological conditions and offers insights into monocyte-driven mechanisms in cardiovascular disease and the potential for monocyte subpopulation targeted therapies.

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.

GABAergic neuronal lineage development determines clinically actionable targets in diffuse hemispheric glioma, H3G34-mutant

Diffuse hemispheric gliomas, H3G34R/V-mutant (DHG-H3G34), are lethal brain tumors lacking targeted therapies. They originate from interneuronal precursors; however, leveraging this origin for therapeutic insights remains unexplored. Here, we delineate a cellular hierarchy along the interneuron lineage development continuum, revealing that DHG-H3G34 mirror spatial patterns of progenitor streams surrounding interneuron nests, as seen during human brain development. Integrating these findings with genome-wide CRISPR-Cas9 screens identifies genes upregulated in interneuron lineage progenitors as major dependencies. Among these, CDK6 emerges as a targetable vulnerability: DHG-H3G34 tumor cells show enhanced sensitivity to CDK4/6 inhibitors and a CDK6-specific degrader, promoting a shift toward more mature interneuron-like states, reducing tumor growth, and prolonging xenograft survival. Notably, a patient with progressive DHG-H3G34 treated with a CDK4/6 inhibitor achieved 17 months of stable disease. This study underscores interneuronal progenitor-like states, organized in characteristic niches, as a distinct vulnerability in DHG-H3G34, highlighting CDK6 as a promising clinically actionable target.

Spatiotemporal transcriptomic mapping of regenerative inflammation in skeletal muscle reveals a dynamic multilayered tissue architecture

Tissue regeneration is orchestrated by macrophages that clear damaged cells and promote regenerative inflammation. How macrophages spatially adapt and diversify their functions to support the architectural requirements of actively regenerating tissue remains unknown. In this study, we reconstructed the dynamic trajectories of myeloid cells isolated from acutely injured and early stage dystrophic muscles. We identified divergent subsets of monocytes/macrophages and DCs and validated markers (e.g., glycoprotein NMB [GPNMB]) and transcriptional regulators associated with defined functional states. In dystrophic muscle, specialized repair-associated subsets exhibited distinct macrophage diversity and reduced DC heterogeneity. Integrating spatial transcriptomics analyses with immunofluorescence uncovered the ordered distribution of subpopulations and multilayered regenerative inflammation zones (RIZs) where distinct macrophage subsets are organized in functional zones around damaged myofibers supporting all phases of regeneration. Importantly, intermittent glucocorticoid treatment disrupted the RIZs. Our findings suggest that macrophage subtypes mediated the development of the highly ordered architecture of regenerative tissues, unveiling the principles of the structured yet dynamic nature of regenerative inflammation supporting effective tissue repair.

IRF5 mediates adaptive immunity via altered glutamine metabolism, mTORC1 signaling and post-transcriptional regulation following T cell receptor activation

Although dynamic alterations in transcriptional, translational, and metabolic programs have been described in T cells, the factors and pathways guiding these molecular shifts are poorly understood, with recent studies revealing a disassociation between transcriptional responses and protein expression following T cell receptor (TCR) stimulation. Previous studies identified interferon regulatory factor 5 (IRF5) in the transcriptional regulation of cytokines, chemotactic molecules and T effector transcription factors following TCR signaling. In this study, we identified T cell intrinsic IRF5 regulation of mTORC1 activity as a key modulator of CD40L protein expression. We further demonstrated a global shift in T cell metabolism, with alterations in glutamine metabolism accompanied by shifts in T cell populations at the single cell level due to loss of Irf5. T cell conditional Irf5 knockout mice in a murine model of experimental autoimmune encephalomyelitis (EAE) demonstrated protection from clinical disease with conserved defects in mTORC1 activity and glutamine regulation. Together, these findings expand our mechanistic understanding of IRF5 as an intrinsic regulator of T effector function(s) and support the therapeutic targeting of IRF5 in multiple sclerosis.

Baf155 controls hematopoietic differentiation and regeneration through chromatin priming

Chromatin priming promotes cell-type-specific gene expression, lineage differentiation, and development. The mechanism of chromatin priming has not been fully understood. Here, we report that mouse hematopoietic stem and progenitor cells (HSPCs) lacking the Baf155 subunit of the BAF (BRG1/BRM-associated factor) chromatin remodeling complex produce a significantly reduced number of mature blood cells, leading to a failure of hematopoietic regeneration upon transplantation and 5-fluorouracil (5-FU) injury. Baf155-deficient HSPCs generate particularly fewer neutrophils, B cells, and CD8+ T cells at homeostasis, supporting a more immune-suppressive tumor microenvironment and enhanced tumor growth. Single-nucleus multiomics analysis reveals that Baf155-deficient HSPCs fail to establish accessible chromatin in selected regions that are enriched for putative enhancers and binding motifs of hematopoietic lineage transcription factors. Our study provides a fundamental mechanistic understanding of the role of Baf155 in hematopoietic lineage chromatin priming and the functional consequences of Baf155 deficiency in regeneration and tumor immunity.

Accurate long-read transcript discovery and quantification at single-cell, pseudo-bulk and bulk resolution with Isosceles

Accurate detection and quantification of mRNA isoforms from nanopore long-read sequencing remains challenged by technical noise, particularly in single cells. To address this, we introduce Isosceles, a computational toolkit that outperforms other methods in isoform detection sensitivity and quantification accuracy across single-cell, pseudo-bulk and bulk resolution levels, as demonstrated using synthetic and biologically-derived datasets. Here we show Isosceles improves the fidelity of single-cell transcriptome quantification at the isoform-level, and enables flexible downstream analysis. As a case study, we apply Isosceles, uncovering coordinated splicing within and between neuronal differentiation lineages. Isosceles is suitable to be applied in diverse biological systems, facilitating studies of cellular heterogeneity across biomedical research applications.