Identifying gene expression programs of cell-type identity and cellular activity with single-cell RNA-Seq

Spatial oncology: Translating contextual biology to the clinic

Microscopic examination of cells in their tissue context has been the driving force behind diagnostic histopathology over the past two centuries. Recently, the rise of advanced molecular biomarkers identified through single cell profiling has increased our understanding of cellular heterogeneity in cancer but have yet to significantly impact clinical care. Spatial technologies integrating molecular profiling with microenvironmental features are poised to bridge this translational gap by providing critical in situ context for understanding cellular interactions and organization. Here, we review how spatial tools have been used to study tumor ecosystems and their clinical applications. We detail findings in cell-cell interactions, microenvironment composition, and tissue remodeling for immune evasion and therapeutic resistance. Additionally, we highlight the emerging role of multi-omic spatial profiling for characterizing clinically relevant features including perineural invasion, tertiary lymphoid structures, and the tumor-stroma interface. Finally, we explore strategies for clinical integration and their augmentation of therapeutic and diagnostic approaches.

The conserved genetic program of male germ cells uncovers ancient regulators of human spermatogenesis

Male germ cells share a common origin across animal species, therefore they likely retain a conserved genetic program that defines their cellular identity. However, the unique evolutionary dynamics of male germ cells coupled with their widespread leaky transcription pose significant obstacles to the identification of the core spermatogenic program. Through network analysis of the spermatocyte transcriptome of vertebrate and invertebrate species, we describe the conserved evolutionary origin of metazoan male germ cells at the molecular level. We estimate the average functional requirement of a metazoan male germ cell to correspond to the expression of approximately 10,000 protein-coding genes, a third of which defines a genetic scaffold of deeply conserved genes that has been retained throughout evolution. Such scaffold contains a set of 79 functional associations between 104 gene expression regulators that represent a core component of the conserved genetic program of metazoan spermatogenesis. By genetically interfering with the acquisition and maintenance of male germ cell identity, we uncover 161 previously unknown spermatogenesis genes and three new potential genetic causes of human infertility. These findings emphasize the importance of evolutionary history on human reproductive disease and establish a cross-species analytical pipeline that can be repurposed to other cell types and pathologies.

SIMVI reveals intrinsic and spatial-induced states in spatial omics data

Spatial omics technologies enable the analysis of gene expression and interaction dynamics in relation to tissue structure and function. However, existing computational methods may not properly distinguish cellular intrinsic variability and intercellular interactions, and may thus fail to capture spatial regulations for further biological discoveries. Here, we present Spatial Interaction Modeling using Variational Inference (SIMVI), an annotation-free framework that disentangles cell intrinsic and spatial-induced latent variables for modeling gene expression in spatial omics data. We derive theoretical support for SIMVI in disentangling intrinsic and spatial-induced variations. By this disentanglement, SIMVI enables estimation of spatial effects (SE) at a single-cell resolution, and opens up various opportunities for novel downstream analyses. To demonstrate the potential of SIMVI, we applied SIMVI to spatial omics data from diverse platforms and tissues (MERFISH human cortex, Slide-seqv2 mouse hippocampus, Slide-tags human tonsil, spatial multiome human melanoma, cohort-level CosMx melanoma). In all tested datasets, SIMVI effectively disentangles variations and infers accurate spatial effects compared with alternative methods. Moreover, on these datasets, SIMVI uniquely uncovers complex spatial regulations and dynamics of biological significance. In the human tonsil data, SIMVI illuminates the cyclical spatial dynamics of germinal center B cells during maturation. Applying SIMVI to both RNA and ATAC modalities of the multiome melanoma data reveals potential tumor epigenetic reprogramming states. Application of SIMVI on our newly-collected cohort-level CosMx melanoma dataset uncovers space-and-outcome-dependent macrophage states and the underlying cellular communication machinery in the tumor microenvironments.

Spatiotemporal modeling reveals high-resolution invasion states in glioblastoma

BACKGROUND

Diffuse invasion of glioblastoma cells through normal brain tissue is a key contributor to tumor aggressiveness, resistance to conventional therapies, and dismal prognosis in patients. A deeper understanding of how components of the tumor microenvironment (TME) contribute to overall tumor organization and to programs of invasion may reveal opportunities for improved therapeutic strategies.

RESULTS

Towards this goal, we apply a novel computational workflow to a spatiotemporally profiled GBM xenograft cohort, leveraging the ability to distinguish human tumor from mouse TME to overcome previous limitations in the analysis of diffuse invasion. Our analytic approach, based on unsupervised deconvolution, performs reference-free discovery of cell types and cell activities within the complete GBM ecosystem. We present a comprehensive catalogue of 15 tumor cell programs set within the spatiotemporal context of 90 mouse brain and TME cell types, cell activities, and anatomic structures. Distinct tumor programs related to invasion align with routes of perivascular, white matter, and parenchymal invasion. Furthermore, sub-modules of genes serving as program network hubs are highly prognostic in GBM patients.

CONCLUSION

The compendium of programs presented here provides a basis for rational targeting of tumor and/or TME components. We anticipate that our approach will facilitate an ecosystem-level understanding of the immediate and long-term consequences of such perturbations, including the identification of compensatory programs that will inform improved combinatorial therapies.

Scalable, compressed phenotypic screening using pooled perturbations

High-throughput phenotypic screens using biochemical perturbations and high-content readouts are constrained by limitations of scale. To address this, we establish a method of pooling exogenous perturbations followed by computational deconvolution to reduce required sample size, labor and cost. We demonstrate the increased efficiency of compressed experimental designs compared to conventional approaches through benchmarking with a bioactive small-molecule library and a high-content imaging readout. We then apply compressed screening in two biological discovery campaigns. In the first, we use early-passage pancreatic cancer organoids to map transcriptional responses to a library of recombinant tumor microenvironment protein ligands, uncovering reproducible phenotypic shifts induced by specific ligands distinct from canonical reference signatures and correlated with clinical outcome. In the second, we identify the pleotropic modulatory effects of a chemical compound library with known mechanisms of action on primary human peripheral blood mononuclear cell immune responses. In sum, our approach empowers phenotypic screens with information-rich readouts to advance drug discovery efforts and basic biological inquiry.

C-ziptf: stable tensor factorization for zero-inflated multi-dimensional genomics data

In the past two decades, genomics has advanced significantly, with single-cell RNA-sequencing (scRNA-seq) marking a pivotal milestone. ScRNA-seq provides unparalleled insights into cellular diversity and has spurred diverse studies across multiple conditions and samples, resulting in an influx of complex multidimensional genomics data. This highlights the need for robust methodologies capable of handling the complexity and multidimensionality of such genomics data. Furthermore, single-cell data grapples with sparsity due to issues like low capture efficiency and dropout effects. Tensor factorizations (TF) have emerged as powerful tools to unravel the complex patterns from multi-dimensional genomics data. Classic TF methods, based on maximum likelihood estimation, struggle with zero-inflated count data, while the inherent stochasticity in TFs further complicates result interpretation and reproducibility. Our paper introduces Zero Inflated Poisson Tensor Factorization (ZIPTF), a novel method for high-dimensional zero-inflated count data factorization. We also present Consensus-ZIPTF (C-ZIPTF), merging ZIPTF with a consensus-based approach to address stochasticity. We evaluate our proposed methods on synthetic zero-inflated count data, simulated scRNA-seq data, and real multi-sample multi-condition scRNA-seq datasets. ZIPTF consistently outperforms baseline matrix and tensor factorization methods, displaying enhanced reconstruction accuracy for zero-inflated data. When dealing with high probabilities of excess zeros, ZIPTF achieves up to 2.4 × better accuracy. Moreover, C-ZIPTF notably enhances the factorization's consistency. When tested on synthetic and real scRNA-seq data, ZIPTF and C-ZIPTF consistently uncover known and biologically meaningful gene expression programs. Access our data and code at: https://github.com/klarman-cell-observatory/scBTF and https://github.com/klarman-cell-observatory/scbtf_experiments .

Network modeling links kidney developmental programs and the cancer type-specificity of VHL mutations

Elucidating the molecular dependencies behind the cancer-type specificity of driver mutations may reveal new therapeutic opportunities. We hypothesized that developmental programs would impact the transduction of oncogenic signaling activated by a driver mutation and shape its cancer-type specificity. Therefore, we designed a computational analysis framework by combining single-cell gene expression profiles during fetal organ development, latent factor discovery, and information theory-based differential network analysis to systematically identify transcription factors that selectively respond to driver mutations under the influence of organ-specific developmental programs. After applying this approach to VHL mutations, which are highly specific to clear cell renal cell carcinoma (ccRCC), we revealed important regulators downstream of VHL mutations in ccRCC and used their activities to cluster patients with ccRCC into three subtypes. This classification revealed a more significant difference in prognosis than the previous mRNA profile-based method and was validated in an independent cohort. Moreover, we found that EP300, a key epigenetic factor maintaining the regulatory network of the subtype with the worst prognosis, can be targeted by a small inhibitor, suggesting a potential treatment option for a subset of patients with ccRCC. This work demonstrated an intimate relationship between organ development and oncogenesis from the perspective of systems biology, and the method can be generalized to study the influence of other biological processes on cancer driver mutations.

Organ Evolution: Emergence of Multicellular Function

Instances of multicellularity across the tree of life have fostered the evolution of complex organs composed of distinct cell types that cooperate, producing emergent biological functions. How organs originate is a fundamental evolutionary problem that has eluded deep mechanistic and conceptual understanding. Here I propose a cell- to organ-level transitions framework, whereby cooperative division of labor originates and becomes entrenched between cell types through a process of functional niche creation, cell-type subfunctionalization, and irreversible ratcheting of cell interdependencies. Comprehending this transition hinges on explaining how these processes unfold molecularly in evolving populations. Recent single-cell transcriptomic studies and analyses of terminal fate specification indicate that cellular functions are conferred by modular gene expression programs. These discrete components of functional variation may be deployed or combined within cells to introduce new properties into multicellular niches, or partitioned across cells to establish division of labor. Tracing gene expression program evolution at the level of single cells in populations may reveal transitions toward organ complexity.

Inferring pattern-driving intercellular flows from single-cell and spatial transcriptomics

From single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), one can extract high-dimensional gene expression patterns that can be described by intercellular communication networks or decoupled gene modules. These two descriptions of information flow are often assumed to occur independently. However, intercellular communication drives directed flows of information that are mediated by intracellular gene modules, in turn triggering outflows of other signals. Methodologies to describe such intercellular flows are lacking. We present FlowSig, a method that infers communication-driven intercellular flows from scRNA-seq or ST data using graphical causal modeling and conditional independence. We benchmark FlowSig using newly generated experimental cortical organoid data and synthetic data generated from mathematical modeling. We demonstrate FlowSig's utility by applying it to various studies, showing that FlowSig can capture stimulation-induced changes to paracrine signaling in pancreatic islets, demonstrate shifts in intercellular flows due to increasing COVID-19 severity and reconstruct morphogen-driven activator-inhibitor patterns in mouse embryogenesis.

INSPIRE: interpretable, flexible and spatially-aware integration of multiple spatial transcriptomics datasets from diverse sources

Recent advances in spatial transcriptomics technologies have led to a growing number of diverse datasets, offering unprecedented opportunities to explore tissue organizations and functions within spatial contexts. However, it remains a significant challenge to effectively integrate and interpret these data, often originating from different samples, technologies, and developmental stages. In this paper, we present INSPIRE, a deep learning method for integrative analyses of multiple spatial transcriptomics datasets to address this challenge. With designs of graph neural networks and an adversarial learning mechanism, INSPIRE enables spatially informed and adaptable integration of data from varying sources. By incorporating non-negative matrix factorization, INSPIRE uncovers interpretable spatial factors with corresponding gene programs, revealing tissue architectures, cell type distributions and biological processes. We demonstrate the capabilities of INSPIRE by applying it to human cortex slices from different samples, mouse brain slices with complementary views, mouse hippocampus and embryo slices generated through different technologies, and spatiotemporal organogenesis atlases containing half a million spatial spots. INSPIRE shows superior performance in identifying detailed biological signals, effectively borrowing information across distinct profiling technologies, and elucidating dynamical changes during embryonic development. Furthermore, we utilize INSPIRE to build 3D models of tissues and whole organisms from multiple slices, demonstrating its power and versatility.

Intercellular extrachromosomal DNA copy-number heterogeneity drives neuroblastoma cell state diversity

Neuroblastoma exhibits significant inter- and intra-tumor genetic heterogeneity and varying clinical outcomes. Extrachromosomal DNAs (ecDNAs) may drive this heterogeneity by independently segregating during cell division, leading to rapid oncogene amplification. While ecDNA-mediated oncogene amplification is linked to poor prognosis in various cancers, the effects of ecDNA copy-number heterogeneity on intermediate phenotypes are poorly understood. Here, we leverage DNA and RNA sequencing from the same single cells in cell lines and neuroblastoma patients to investigate these effects. By analyzing ecDNA amplicon structures, we reveal extensive intercellular ecDNA copy-number heterogeneity. We also provide direct evidence of how this heterogeneity influences the expression of cargo genes, including MYCN and its downstream targets, and the overall transcriptional state of neuroblastoma cells. Our findings highlight the role of ecDNA copy number in promoting rapid adaptability of cellular states within tumors, underscoring the need for ecDNA-specific treatment strategies to address tumor formation and adaptation.

Unraveling the phenotypic states of human innate-like T cells: Comparative insights with conventional T cells and mouse models

The "innate-like" T cell compartment, known as Tinn, represents a diverse group of T cells that straddle the boundary between innate and adaptive immunity. We explore the transcriptional landscape of Tinn compared to conventional T cells (Tconv) in the human thymus and blood using single-cell RNA sequencing (scRNA-seq) and flow cytometry. In human blood, the majority of Tinn cells share an effector program driven by specific transcription factors, distinct from those governing Tconv cells. Conversely, only a fraction of thymic Tinn cells displays an effector phenotype, while others share transcriptional features with developing Tconv cells, indicating potential divergent developmental pathways. Unlike the mouse, human Tinn cells do not differentiate into multiple effector subsets but develop a mixed type 1/type 17 effector potential. Cross-species analysis uncovers species-specific distinctions, including the absence of type 2 Tinn cells in humans, which implies distinct immune regulatory mechanisms across species.

Cell type and cell signaling innovations underlying mammalian pregnancy

How fetal and maternal cell types have co-evolved to enable mammalian placentation poses a unique evolutionary puzzle. Here, we present a multi-species atlas integrating single-cell transcriptomes from six species bracketing therian mammal diversity. We find that invasive trophoblasts share a gene-expression signature across eutherians, and evidence that endocrine decidual cells evolved stepwise from an immunomodulatory cell type retained in Tenrec with affinity to human decidua of menstruation. We recover evolutionary patterns in ligand-receptor signaling: fetal and maternal cells show a pronounced tendency towards disambiguation, but a predicted arms race dynamic between them is limited. We reconstruct cell communication networks of extinct mammalian ancestors, finding strong integration of fetal trophoblast into maternal networks. Together, our results reveal a dynamic history of cell type and signaling evolution.

Synopsis

The fetal-maternal interface is one of the most intense loci of cell-cell signaling in the human body. Invasion of cells from the fetal placenta into the uterus, and the corresponding transformation of maternal tissues called decidualization, first evolved in the stem lineage of eutherian mammals( 1 , 2 ). Single-cell studies of the human fetal-maternal interface have provided new insight into the cell type diversity and cell-cell interactions governing this chimeric organ( 3-5 ). However, the fetal-maternal interface is also one of the most rapidly evolving, and hence most diverse, characters among mammals( 6 ), and an evolutionary analysis is missing. Here, we present and compare single-cell data from the fetal-maternal interface of species bracketing key events in mammal phylogeny: a marsupial (opossum, Monodelphis domestica ), the afrotherian Tenrec ecaudatus, and four Euarchontoglires - guinea pig and mouse (Rodentia) together with recent macaque and human data (primates) ( 4 , 5 , 7 ). We infer cell type homologies, identify a gene-expression signature of eutherian invasive trophoblast conserved over 99 million years, and discover a predecidual cell in the tenrec which suggests stepwise evolution of the decidual stromal cell. We reconstruct ancestral cell signaling networks, revealing the integration of fetal cell types into the interface. Finally, we test two long-standing theoretical predictions, the disambiguation hypothesis( 8 ) and escalation hypothesis( 9 ), at transcriptome-wide scale, finding divergence between fetal and maternal signaling repertoires but arms race dynamics restricted to a small subset of ligand-receptor pairs. In so doing, we trace the co-evolutionary history of cell types and their signaling across mammalian viviparity.

Unsupervised pattern identification in spatial gene expression atlas reveals mouse brain regions beyond established ontology

The rapid growth of large-scale spatial gene expression data demands efficient and reliable computational tools to extract major trends of gene expression in their native spatial context. Here, we used stability-driven unsupervised learning (i.e., staNMF) to identify principal patterns (PPs) of 3D gene expression profiles and understand spatial gene distribution and anatomical localization at the whole mouse brain level. Our subsequent spatial correlation analysis systematically compared the PPs to known anatomical regions and ontology from the Allen Mouse Brain Atlas using spatial neighborhoods. We demonstrate that our stable and spatially coherent PPs, whose linear combinations accurately approximate the spatial gene data, are highly correlated with combinations of expert-annotated brain regions. These PPs yield a brain ontology based purely on spatial gene expression. Our PP identification approach outperforms principal component analysis and typical clustering algorithms on the same task. Moreover, we show that the stable PPs reveal marked regional imbalance of brainwide genetic architecture, leading to region-specific marker genes and gene coexpression networks. Our findings highlight the advantages of stability-driven machine learning for plausible biological discovery from dense spatial gene expression data, streamlining tasks that are infeasible by conventional manual approaches.

Neoadjuvant PARPi or chemotherapy in ovarian cancer informs targeting effector Treg cells for homologous-recombination-deficient tumors

Homologous recombination deficiency (HRD) is prevalent in cancer, sensitizing tumor cells to poly (ADP-ribose) polymerase (PARP) inhibition. However, the impact of HRD and related therapies on the tumor microenvironment (TME) remains elusive. Our study generates single-cell gene expression and T cell receptor profiles, along with validatory multimodal datasets from >100 high-grade serous ovarian cancer (HGSOC) samples, primarily from a phase II clinical trial (NCT04507841). Neoadjuvant monotherapy with the PARP inhibitor (PARPi) niraparib achieves impressive 62.5% and 73.6% response rates per RECIST v.1.1 and GCIG CA125, respectively. We identify effector regulatory T cells (eTregs) as key responders to HRD and neoadjuvant therapies, co-occurring with other tumor-reactive T cells, particularly terminally exhausted CD8+ T cells (Tex). TME-wide interferon signaling correlates with cancer cells upregulating MHC class II and co-inhibitory ligands, potentially driving Treg and Tex fates. Depleting eTregs in HRD mouse models, with or without PARP inhibition, significantly suppresses tumor growth without observable toxicities, underscoring the potential of eTreg-focused therapeutics for HGSOC and other HRD-related tumors.

SPP1+ TAM Regulates the Metastatic Colonization of CXCR4+ Metastasis-Associated Tumor Cells by Remodeling the Lymph Node Microenvironment

Lymph node metastasis, the initial step in distant metastasis, represents a primary contributor to mortality in patients diagnosed with oral squamous cell carcinoma (OSCC). However, the underlying mechanisms of lymph node metastasis in OSCC remain incompletely understood. Here, the transcriptomes of 56 383 single cells derived from paired tissues of six OSCC patients are analyzed. This study founds that CXCR4+ epithelial cells, identified as highly malignant disseminated tumor cells (DTCs), exhibited a propensity for lymph node metastasis. Importantly, a distinct subset of tumor-associated macrophages (TAMs) characterized by exclusive expression of phosphoprotein 1 (SPP1) is discovered. These TAMs may remodel the metastatic lymph node microenvironment by potentially activating fibroblasts and promoting T cell exhaustion through SPP1-CD44 and CD155-CD226 ligand-receptor interactions, thereby facilitating colonization and proliferation of disseminated tumor cells. The research advanced the mechanistic understanding of metastatic tumor microenvironment (TME) and provided a foundation for the development of personalized treatments for OSCC patients with metastasis.

Deciphering the impact of genomic variation on function

Our genomes influence nearly every aspect of human biology-from molecular and cellular functions to phenotypes in health and disease. Studying the differences in DNA sequence between individuals (genomic variation) could reveal previously unknown mechanisms of human biology, uncover the basis of genetic predispositions to diseases, and guide the development of new diagnostic tools and therapeutic agents. Yet, understanding how genomic variation alters genome function to influence phenotype has proved challenging. To unlock these insights, we need a systematic and comprehensive catalogue of genome function and the molecular and cellular effects of genomic variants. Towards this goal, the Impact of Genomic Variation on Function (IGVF) Consortium will combine approaches in single-cell mapping, genomic perturbations and predictive modelling to investigate the relationships among genomic variation, genome function and phenotypes. IGVF will create maps across hundreds of cell types and states describing how coding variants alter protein activity, how noncoding variants change the regulation of gene expression, and how such effects connect through gene-regulatory and protein-interaction networks. These experimental data, computational predictions and accompanying standards and pipelines will be integrated into an open resource that will catalyse community efforts to explore how our genomes influence biology and disease across populations.

Single-cell transcriptomic analysis of immune cell dynamics in the healthy human endometrium

The microenvironment of the endometrial immune system is crucial to the success of placental implantation and healthy pregnancy. However, the functionalities of immune cells across various stages of the reproductive cycle have yet to be fully comprehended. To address this, we conducted advanced bioinformatic analysis on 230,049 high-quality single-cell transcriptomes from healthy endometrial samples obtained during the proliferative, secretory, early pregnancy, and late pregnancy stages. Our investigation has unveiled that proliferative natural killer (NK) cells, a potential source of endometrial NK cells, exhibit the most robust proliferative and differentiation potential during non-pregnant stages. We have also identified similar differentiation trajectories of NK cells originating from proliferative NK cells across four stages. Notably, during early pregnancy, NK cells demonstrate the highest oxidative phosphorylation metabolism activity, and, in conjunction with macrophages and T cells, exhibit the strongest type II interferon response. With spatial transcriptome data, we have discerned that the most robust immune-non-immune interactions are associated with the promotion and inhibition of cell proliferation, differentiation and migration during four stages. Furthermore, we have compiled lists of stage-specific risk genes implicated in reproductive diseases, which hold promise as potential disease biomarkers. Our study provides insights into the dynamics of the endometrial immune microenvironment during different reproductive cycle stages, thus serving as a reference for detecting pathological changes during pregnancy.

Dissecting the regulatory logic of specification and differentiation during vertebrate embryogenesis

The interplay between transcription factors and chromatin accessibility regulates cell type diversification during vertebrate embryogenesis. To systematically decipher the gene regulatory logic guiding this process, we generated a single-cell multi-omics atlas of RNA expression and chromatin accessibility during early zebrafish embryogenesis. We developed a deep learning model to predict chromatin accessibility based on DNA sequence and found that a small number of transcription factors underlie cell-type-specific chromatin landscapes. While Nanog is well-established in promoting pluripotency, we discovered a new function in priming the enhancer accessibility of mesendodermal genes. In addition to the classical stepwise mode of differentiation, we describe instant differentiation, where pluripotent cells skip intermediate fate transitions and terminally differentiate. Reconstruction of gene regulatory interactions reveals that this process is driven by a shallow network in which maternally deposited regulators activate a small set of transcription factors that co-regulate hundreds of differentiation genes. Notably, misexpression of these transcription factors in pluripotent cells is sufficient to ectopically activate their targets. This study provides a rich resource for analyzing embryonic gene regulation and reveals the regulatory logic of instant differentiation.

Cell states and neighborhoods in distinct clinical stages of primary and metastatic esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is a highly lethal cancer of the upper gastrointestinal tract with rising incidence in western populations. To decipher EAC disease progression and therapeutic response, we performed multiomic analyses of a cohort of primary and metastatic EAC tumors, incorporating single-nuclei transcriptomic and chromatin accessibility sequencing, along with spatial profiling. We identified tumor microenvironmental features previously described to associate with therapy response. We identified five malignant cell programs, including undifferentiated, intermediate, differentiated, epithelial-to-mesenchymal transition, and cycling programs, which were associated with differential epigenetic plasticity and clinical outcomes, and for which we inferred candidate transcription factor regulons. Furthermore, we revealed diverse spatial localizations of malignant cells expressing their associated transcriptional programs and predicted their significant interactions with microenvironmental cell types. We validated our findings in three external single-cell RNA-seq and three bulk RNA-seq studies. Altogether, our findings advance the understanding of EAC heterogeneity, disease progression, and therapeutic response.

Applying 12 machine learning algorithms and Non-negative Matrix Factorization for robust prediction of lupus nephritis

Lupus nephritis (LN) is a challenging condition with limited diagnostic and treatment options. In this study, we applied 12 distinct machine learning algorithms along with Non-negative Matrix Factorization (NMF) to analyze single-cell datasets from kidney biopsies, aiming to provide a comprehensive profile of LN. Through this analysis, we identified various immune cell populations and their roles in LN progression and constructed 102 machine learning-based immune-related gene (IRG) predictive models. The most effective models demonstrated high predictive accuracy, evidenced by Area Under the Curve (AUC) values, and were further validated in external cohorts. These models highlight six hub IRGs (CD14, CYBB, IFNGR1, IL1B, MSR1, and PLAUR) as key diagnostic markers for LN, showing remarkable diagnostic performance in both renal and peripheral blood cohorts, thus offering a novel approach for noninvasive LN diagnosis. Further clinical correlation analysis revealed that expressions of IFNGR1, PLAUR, and CYBB were negatively correlated with the glomerular filtration rate (GFR), while CYBB also positively correlated with proteinuria and serum creatinine levels, highlighting their roles in LN pathophysiology. Additionally, protein-protein interaction (PPI) analysis revealed significant networks involving hub IRGs, emphasizing the importance of the interleukin family and chemokines in LN pathogenesis. This study highlights the potential of integrating advanced genomic tools and machine learning algorithms to improve diagnosis and personalize management of complex autoimmune diseases like LN.

A multi-layered integrative analysis reveals a cholesterol metabolic program in outer radial glia with implications for human brain evolution

The definition of molecular and cellular mechanisms contributing to brain ontogenetic trajectories is essential to investigate the evolution of our species. Yet their functional dissection at an appropriate level of granularity remains challenging. Capitalizing on recent efforts that have extensively profiled neural stem cells from the developing human cortex, we develop an integrative computational framework to perform trajectory inference and gene regulatory network reconstruction, (pseudo)time-informed non-negative matrix factorization for learning the dynamics of gene expression programs, and paleogenomic analysis for a higher-resolution mapping of derived regulatory variants in our species in comparison with our closest relatives. We provide evidence for cell type-specific regulation of gene expression programs during indirect neurogenesis. In particular, our analysis uncovers a key role for a cholesterol program in outer radial glia, regulated by zinc-finger transcription factor KLF6. A cartography of the regulatory landscape impacted by Homo sapiens-derived variants reveals signals of selection clustering around regulatory regions associated with GLI3, a well-known regulator of radial glial cell cycle, and impacting KLF6 regulation. Our study contributes to the evidence of significant changes in metabolic pathways in recent human brain evolution.

Comprehensive single-cell analysis deciphered microenvironmental dynamics and immune regulator olfactomedin 4 in pathogenesis of gallbladder cancer

OBJECTIVE

Elucidating complex ecosystems and molecular features of gallbladder cancer (GBC) and benign gallbladder diseases is pivotal to proactive cancer prevention and optimal therapeutic intervention.

DESIGN

We performed single-cell transcriptome analysis on 230 737 cells from 15 GBCs, 4 cholecystitis samples, 3 gallbladder polyps, 5 gallbladder adenomas and 16 adjacent normal tissues. Findings were validated through large-scale histological assays, digital spatial profiler multiplexed immunofluorescence (GeoMx), etc. Further molecular mechanism was demonstrated with in vitro and in vivo studies.

RESULTS

The cell atlas unveiled an altered immune landscape across different pathological states of gallbladder diseases. GBC featured a more suppressive immune microenvironment with distinct T-cell proliferation patterns and macrophage attributions in different GBC subtypes. Notably, mutual exclusivity between stromal and immune cells was identified and remarkable stromal ecosystem (SC) heterogeneity during GBC progression was unveiled. Specifically, SC1 demonstrated active interaction between Fibro-iCAF and Endo-Tip cells, correlating with poor prognosis. Moreover, epithelium genetic variations within adenocarcinoma (AC) indicated an evolutionary similarity between adenoma and AC. Importantly, our study identified elevated olfactomedin 4 (OLFM4) in epithelial cells as a central player in GBC progression. OLFM4 was related to T-cell malfunction and tumour-associated macrophage infiltration, leading to a worse prognosis in GBC. Further investigations revealed that OLFM4 upregulated programmed death-ligand 1 (PD-L1) expression through the MAPK-AP1 axis, facilitating tumour cell immune evasion.

CONCLUSION

These findings offer a valuable resource for understanding the pathogenesis of gallbladder diseases and indicate OLFM4 as a potential biomarker and therapeutic target for GBC.

Interpretable single-cell factor decomposition using sciRED

Single-cell RNA sequencing (scRNA-seq) maps gene expression heterogeneity within a tissue. However, identifying biological signals in this data is challenging due to confounding technical factors, sparsity, and high dimensionality. Data factorization methods address this by separating and identifying signals in the data, such as gene expression programs, but the resulting factors must be manually interpreted. We developed Single-Cell Interpretable Residual Decomposition (sciRED) to improve the interpretation of scRNA-seq factor analysis. sciRED removes known confounding effects, uses rotations to improve factor interpretability, maps factors to known covariates, identifies unexplained factors that may capture hidden biological phenomena and determines the genes and biological processes represented by the resulting factors. We apply sciRED to multiple scRNA-seq datasets and identify sex-specific variation in a kidney map, discern strong and weak immune stimulation signals in a PBMC dataset, reduce ambient RNA contamination in a rat liver atlas to help identify strain variation, and reveal rare cell type signatures and anatomical zonation gene programs in a healthy human liver map. These demonstrate that sciRED is useful in characterizing diverse biological signals within scRNA-seq data.

Interpretable single-cell factor decomposition using sciRED

Single-cell RNA sequencing (scRNA-seq) maps gene expression heterogeneity within a tissue. However, identifying biological signals in this data is challenging due to confounding technical factors, sparsity, and high dimensionality. Data factorization methods address this by separating and identifying signals in the data, such as gene expression programs, but the resulting factors must be manually interpreted. We developed Single-Cell Interpretable Residual Decomposition (sciRED) to improve the interpretation of scRNA-seq factor analysis. sciRED removes known confounding effects, uses rotations to improve factor interpretability, maps factors to known covariates, identifies unexplained factors that may capture hidden biological phenomena and determines the genes and biological processes represented by the resulting factors. We apply sciRED to multiple scRNA-seq datasets and identify sex-specific variation in a kidney map, discern strong and weak immune stimulation signals in a PBMC dataset, reduce ambient RNA contamination in a rat liver atlas to help identify strain variation, and reveal rare cell type signatures and anatomical zonation gene programs in a healthy human liver map. These demonstrate that sciRED is useful in characterizing diverse biological signals within scRNA-seq data.

FASN-mediated fatty acid biosynthesis remodels immune environment in Clonorchis sinensis infection-related intrahepatic cholangiocarcinoma

BACKGROUND & AIMS

Intrahepatic cholangiocarcinoma (iCCA) is the second most common primary liver cancer and is highly lethal. Clonorchis sinensis (C. sinensis) infection is an important risk factor for iCCA. Here we investigated the clinical impact and underlying molecular characteristics of C. sinensis infection-related iCCA.

METHODS

We performed single-cell RNA sequencing, whole-exome sequencing, RNA sequencing, metabolomics and spatial transcriptomics in 251 patients with iCCA from three medical centers. Alterations in metabolism and the immune microenvironment of C. sinensis-related iCCAs were validated through an in vitro co-culture system and in a mouse model of iCCA.

RESULTS

We revealed that C. sinensis infection was significantly associated with iCCA patients' overall survival and response to immunotherapy. Fatty acid biosynthesis and the expression of fatty acid synthase (FASN), a key enzyme catalyzing long-chain fatty acid synthesis, were significantly enriched in C. sinensis-related iCCAs. iCCA cell lines treated with excretory/secretory products of C. sinensis displayed elevated FASN and free fatty acids. The metabolic alteration of tumor cells was closely correlated with the enrichment of tumor-associated macrophage (TAM)-like macrophages and the impaired function of T cells, which led to formation of an immunosuppressive microenvironment and tumor progression. Spatial transcriptomics analysis revealed that malignant cells were in closer juxtaposition with TAM-like macrophages in C. sinensis-related iCCAs than non-C. sinensis-related iCCAs. Importantly, treatment with a FASN inhibitor significantly reversed the immunosuppressive microenvironment and enhanced anti-PD-1 efficacy in iCCA mouse models treated with excretory/secretory products from C. sinensis.

CONCLUSIONS

We provide novel insights into metabolic alterations and the immune microenvironment in C. sinensis infection-related iCCAs. We also demonstrate that the combination of a FASN inhibitor with immunotherapy could be a promising strategy for the treatment of C. sinensis-related iCCAs.

IMPACT AND IMPLICATIONS

Clonorchis sinensis (C. sinensis)-infected patients with intrahepatic cholangiocarcinoma (iCCA) have a worse prognosis and response to immunotherapy than non-C. sinensis-infected patients with iCCA. The underlying molecular characteristics of C. sinensis infection-related iCCAs remain unclear. Herein, we demonstrate that upregulation of FASN (fatty acid synthase) and free fatty acids in C. sinensis-related iCCAs leads to formation of an immunosuppressive microenvironment and tumor progression. Thus, administration of FASN inhibitors could significantly reverse the immunosuppressive microenvironment and further enhance the efficacy of anti-PD-1 against C. sinensis-related iCCAs.

Single-cell multiregion dissection of Alzheimer's disease

Alzheimer's disease is the leading cause of dementia worldwide, but the cellular pathways that underlie its pathological progression across brain regions remain poorly understood1-3. Here we report a single-cell transcriptomic atlas of six different brain regions in the aged human brain, covering 1.3 million cells from 283 post-mortem human brain samples across 48 individuals with and without Alzheimer's disease. We identify 76 cell types, including region-specific subtypes of astrocytes and excitatory neurons and an inhibitory interneuron population unique to the thalamus and distinct from canonical inhibitory subclasses. We identify vulnerable populations of excitatory and inhibitory neurons that are depleted in specific brain regions in Alzheimer's disease, and provide evidence that the Reelin signalling pathway is involved in modulating the vulnerability of these neurons. We develop a scalable method for discovering gene modules, which we use to identify cell-type-specific and region-specific modules that are altered in Alzheimer's disease and to annotate transcriptomic differences associated with diverse pathological variables. We identify an astrocyte program that is associated with cognitive resilience to Alzheimer's disease pathology, tying choline metabolism and polyamine biosynthesis in astrocytes to preserved cognitive function late in life. Together, our study develops a regional atlas of the ageing human brain and provides insights into cellular vulnerability, response and resilience to Alzheimer's disease pathology.

Interface-guided phenotyping of coding variants in the transcription factor RUNX1

Single-gene missense mutations remain challenging to interpret. Here, we deploy scalable functional screening by sequencing (SEUSS), a Perturb-seq method, to generate mutations at protein interfaces of RUNX1 and quantify their effect on activities of downstream cellular programs. We evaluate single-cell RNA profiles of 115 mutations in myelogenous leukemia cells and categorize them into three functionally distinct groups, wild-type (WT)-like, loss-of-function (LoF)-like, and hypomorphic, that we validate in orthogonal assays. LoF-like variants dominate the DNA-binding site and are recurrent in cancer; however, recurrence alone does not predict functional impact. Hypomorphic variants share characteristics with LoF-like but favor protein interactions, promoting gene expression indicative of nerve growth factor (NGF) response and cytokine recruitment of neutrophils. Accessible DNA near differentially expressed genes frequently contains RUNX1-binding motifs. Finally, we reclassify 16 variants of uncertain significance and train a classifier to predict 103 more. Our work demonstrates the potential of targeting protein interactions to better define the landscape of phenotypes reachable by missense mutations.

The neurotransmitter calcitonin gene-related peptide shapes an immunosuppressive microenvironment in medullary thyroid cancer

Neurotransmitters are key modulators in neuro-immune circuits and have been linked to tumor progression. Medullary thyroid cancer (MTC), an aggressive neuroendocrine tumor, expresses neurotransmitter calcitonin gene-related peptide (CGRP), is insensitive to chemo- and radiotherapies, and the effectiveness of immunotherapies remains unknown. Thus, a comprehensive analysis of the tumor microenvironment would facilitate effective therapies and provide evidence on CGRP's function outside the nervous system. Here, we compare the single-cell landscape of MTC and papillary thyroid cancer (PTC) and find that expression of CGRP in MTC is associated with dendritic cell (DC) abnormal development characterized by activation of cAMP related pathways and high levels of Kruppel Like Factor 2 (KLF2), correlated with an impaired activity of tumor infiltrating T cells. A CGRP receptor antagonist could offset CGRP detrimental impact on DC development in vitro. Our study provides insights of the MTC immunosuppressive microenvironment, and proposes CGRP receptor as a potential therapeutic target.

SingleCellGGM enables gene expression program identification from single-cell transcriptomes and facilitates universal cell label transfer

Gene co-expression analysis of single-cell transcriptomes, aiming to define functional relationships between genes, is challenging due to excessive dropout values. Here, we developed a single-cell graphical Gaussian model (SingleCellGGM) algorithm to conduct single-cell gene co-expression network analysis. When applied to mouse single-cell datasets, SingleCellGGM constructed networks from which gene co-expression modules with highly significant functional enrichment were identified. We considered the modules as gene expression programs (GEPs). These GEPs enable direct cell-type annotation of individual cells without cell clustering, and they are enriched with genes required for the functions of the corresponding cells, sometimes at levels greater than 10-fold. The GEPs are conserved across datasets and enable universal cell-type label transfer across different studies. We also proposed a dimension-reduction method through averaging by GEPs for single-cell analysis, enhancing the interpretability of results. Thus, SingleCellGGM offers a unique GEP-based perspective to analyze single-cell transcriptomes and reveals biological insights shared by different single-cell datasets.

mosaicMPI: a framework for modular data integration across cohorts and -omics modalities

Advances in molecular profiling have facilitated generation of large multi-modal datasets that can potentially reveal critical axes of biological variation underlying complex diseases. Distilling biological meaning, however, requires computational strategies that can perform mosaic integration across diverse cohorts and datatypes. Here, we present mosaicMPI, a framework for discovery of low to high-resolution molecular programs representing both cell types and states, and integration within and across datasets into a network representing biological themes. Using existing datasets in glioblastoma, we demonstrate that this approach robustly integrates single cell and bulk programs across multiple platforms. Clinical and molecular annotations from cohorts are statistically propagated onto this network of programs, yielding a richly characterized landscape of biological themes. This enables deep understanding of individual tumor samples, systematic exploration of relationships between modalities, and generation of a reference map onto which new datasets can rapidly be mapped. mosaicMPI is available at https://github.com/MorrissyLab/mosaicMPI.

Gene regulatory network structure informs the distribution of perturbation effects

Gene regulatory networks (GRNs) govern many core developmental and biological processes underlying human complex traits. Even with broad-scale efforts to characterize the effects of molecular perturbations and interpret gene coexpression, it remains challenging to infer the architecture of gene regulation in a precise and efficient manner. Key properties of GRNs, like hierarchical structure, modular organization, and sparsity, provide both challenges and opportunities for this objective. Here, we seek to better understand properties of GRNs using a new approach to simulate their structure and model their function. We produce realistic network structures with a novel generating algorithm based on insights from small-world network theory, and we model gene expression regulation using stochastic differential equations formulated to accommodate modeling molecular perturbations. With these tools, we systematically describe the effects of gene knockouts within and across GRNs, finding a subset of networks that recapitulate features of a recent genome-scale perturbation study. With deeper analysis of these exemplar networks, we consider future avenues to map the architecture of gene expression regulation using data from cells in perturbed and unperturbed states, finding that while perturbation data are critical to discover specific regulatory interactions, data from unperturbed cells may be sufficient to reveal regulatory programs.

Spatial transcriptomics reveals influence of microenvironment on intrinsic fates in melanoma therapy resistance

Resistance to cancer therapy is driven by both cell-intrinsic and microenvironmental factors. Previous work has revealed that multiple resistant cell fates emerge in melanoma following treatment with targeted therapy and that, in vitro, these resistant fates are determined by the transcriptional state of individual cells prior to exposure to treatment. What remains unclear is whether these resistant fates are shared across different genetic backgrounds and how, if at all, these resistant fates interact with the tumor microenvironment. Through spatial transcriptomics and single-cell RNA sequencing, we uncovered distinct resistance programs in melanoma cells shaped by both intrinsic cellular states and the tumor microenvironment. Consensus non-negative matrix factorization revealed shared intrinsic resistance programs across different cell lines, highlighting the presence of universal and unique resistance pathways. In patient samples, we demonstrated that these resistance programs coexist within individual tumors and associate with diverse immune signatures, suggesting that the tumor microenvironment and distribution of resistant fates are closely connected. Single-cell resolution spatial transcriptomics in xenograft models revealed both intrinsically determined and extrinsically influenced resistant fates. Overall, this work demonstrates that each therapy resistant fate coexists with a distinct immune microenvironment in tumors and that, in vivo, tissue features, such as regions of necrosis, can influence which resistant fate is adopted.

An oncolytic virus delivering tumor-irrelevant bystander T cell epitopes induces anti-tumor immunity and potentiates cancer immunotherapy

Tumor-specific T cells are crucial in anti-tumor immunity and act as targets for cancer immunotherapies. However, these cells are numerically scarce and functionally exhausted in the tumor microenvironment (TME), leading to inefficacious immunotherapies in most patients with cancer. By contrast, emerging evidence suggested that tumor-irrelevant bystander T (TBYS) cells are abundant and preserve functional memory properties in the TME. To leverage TBYS cells in the TME to eliminate tumor cells, we engineered oncolytic virus (OV) encoding TBYS epitopes (OV-BYTE) to redirect the antigen specificity of tumor cells to pre-existing TBYS cells, leading to effective tumor inhibition in multiple preclinical models. Mechanistically, OV-BYTE induced epitope spreading of tumor antigens to elicit more diverse tumor-specific T cell responses. Remarkably, the OV-BYTE strategy targeting human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific T cell memory efficiently inhibited tumor progression in a human tumor cell-derived xenograft model, providing important insights into the improvement of cancer immunotherapies in a large population with a history of SARS-CoV-2 infection or coronavirus disease 2019 (COVID-19) vaccination.

Unfolding the mysteries of heterogeneity from a high-resolution perspective: integration analysis of single-cell multi-omics and spatial omics revealed functionally heterogeneous cancer cells in ccRCC

The genomic landscape of clear cell renal cell carcinoma (ccRCC) has a considerable intra-tumor heterogeneity, which is a significant obstacle in the field of precision oncology and plays a pivotal role in metastasis, recurrence, and therapeutic resistance of cancer. The mechanisms of intra-tumor heterogeneity in ccRCC have yet to be fully established. We integrated single-cell RNA sequencing (scRNA-seq) and transposase-accessible chromatin sequencing (scATAC-seq) data from a single-cell multi-omics perspective. Based on consensus non-negative matrix factorization (cNMF) algorithm, functionally heterogeneous cancer cells were classified into metabolism, inflammatory, and EMT meta programs, with spatial transcriptomics sequencing (stRNA-seq) providing spatial information of such disparate meta programs of cancer cells. The bulk RNA sequencing (RNA-seq) data revealed high clinical prognostic values of functionally heterogeneous cancer cells of three meta programs, with transcription factor regulatory network and motif activities revealing the key transcription factors that regulate functionally heterogeneous ccRCC cells. The interactions between varying meta programs and other cell subpopulations in the microenvironment were investigated. Finally, we assessed the sensitivity of cancer cells of disparate meta programs to different anti-cancer agents. Our findings inform on the intra-tumor heterogeneity of ccRCC and its regulatory networks and offers new perspectives to facilitate the designs of rational therapeutic strategies.

Characterization of transcriptional heterogeneity and novel therapeutic targets using single cell RNA-sequencing of primary and circulating Ewing sarcoma cells

Ewing sarcoma is the second most common bone cancer in children, accounting for 2% of pediatric cancer diagnoses. Patients who present with metastatic disease at the time of diagnosis have a dismal prognosis, compared to the >70% 5-year survival of those with localized disease. Here, we utilized single cell RNA-sequencing to characterize the transcriptional landscape of primary Ewing sarcoma tumors and surrounding tumor microenvironment (TME). Copy-number analysis identified subclonal evolution within patients prior to treatment. Primary tumor samples demonstrate a heterogenous transcriptional landscape with several conserved gene expression programs, including those composed of genes related to proliferation and EWS targets. Single cell RNA-sequencing and immunofluorescence of circulating tumor cells at the time of diagnosis identified TSPAN8 as a novel therapeutic target.

An integrated single-cell RNA-seq map of human neuroblastoma tumors and preclinical models uncovers divergent mesenchymal-like gene expression programs

BACKGROUND

Neuroblastoma is a common pediatric cancer, where preclinical studies suggest that a mesenchymal-like gene expression program contributes to chemotherapy resistance. However, clinical outcomes remain poor, implying we need a better understanding of the relationship between patient tumor heterogeneity and preclinical models.

RESULTS

Here, we generate single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We develop an unsupervised machine learning approach ("automatic consensus nonnegative matrix factorization" (acNMF)) to compare the gene expression programs found in preclinical models to a large cohort of patient tumors. We confirm a weakly expressed, mesenchymal-like program in otherwise adrenergic cancer cells in some pre-treated high-risk patient tumors, but this appears distinct from the presumptive drug-resistance mesenchymal programs evident in cell lines. Surprisingly, however, this weak-mesenchymal-like program is maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 h, suggesting an uncharacterized therapy-escape mechanism.

CONCLUSIONS

Collectively, our findings improve the understanding of how neuroblastoma patient tumor heterogeneity is reflected in preclinical models, provides a comprehensive integrated resource, and a generalizable set of computational methodologies for the joint analysis of clinical and pre-clinical single-cell RNA-seq datasets.

Computational analysis of super-resolved in situ sequencing data reveals genes modified by immune-tumor contact events

Cancer cells can manipulate immune cells and escape from the immune system response. Quantifying the molecular changes that occur when an immune cell touches a tumor cell can increase our understanding of the underlying mechanisms. Recently, it became possible to perform such measurements in situ-for example, using expansion sequencing, which enabled in situ sequencing of genes with super-resolution. We systematically examined whether individual immune cells from specific cell types express genes differently when in physical proximity to individual tumor cells. First, we demonstrated that a dense mapping of genes in situ can be used for the segmentation of cell bodies in 3D, thus improving our ability to detect likely touching cells. Next, we used three different computational approaches to detect the molecular changes that are triggered by proximity: differential expression analysis, tree-based machine learning classifiers, and matrix factorization analysis. This systematic analysis revealed tens of genes, in specific cell types, whose expression separates immune cells that are proximal to tumor cells from those that are not proximal, with a significant overlap between the different detection methods. Remarkably, an order of magnitude more genes are triggered by proximity to tumor cells in CD8 T cells compared to CD4 T cells, in line with the ability of CD8 T cells to directly bind major histocompatibility complex (MHC) class I on tumor cells. Thus, in situ sequencing of an individual biopsy can be used to detect genes likely involved in immune-tumor cell-cell interactions. The data used in this manuscript and the code of the InSituSeg, machine learning, cNMF, and Moran's I methods are publicly available at doi:10.5281/zenodo.7497981.

Decoding temporal heterogeneity in NSCLC through machine learning and prognostic model construction

BACKGROUND

Non-small cell lung cancer (NSCLC) is a prevalent and heterogeneous disease with significant genomic variations between the early and advanced stages. The identification of key genes and pathways driving NSCLC tumor progression is critical for improving the diagnosis and treatment outcomes of this disease.

METHODS

In this study, we conducted single-cell transcriptome analysis on 93,406 cells from 22 NSCLC patients to characterize malignant NSCLC cancer cells. Utilizing cNMF, we classified these cells into distinct modules, thus identifying the diverse molecular profiles within NSCLC. Through pseudotime analysis, we delineated temporal gene expression changes during NSCLC evolution, thus demonstrating genes associated with disease progression. Using the XGBoost model, we assessed the significance of these genes in the pseudotime trajectory. Our findings were validated by using transcriptome sequencing data from The Cancer Genome Atlas (TCGA), supplemented via LASSO regression to refine the selection of characteristic genes. Subsequently, we established a risk score model based on these genes, thus providing a potential tool for cancer risk assessment and personalized treatment strategies.

RESULTS

We used cNMF to classify malignant NSCLC cells into three functional modules, including the metabolic reprogramming module, cell cycle module, and cell stemness module, which can be used for the functional classification of malignant tumor cells in NSCLC. These findings also indicate that metabolism, the cell cycle, and tumor stemness play important driving roles in the malignant evolution of NSCLC. We integrated cNMF and XGBoost to select marker genes that are indicative of both early and advanced NSCLC stages. The expression of genes such as CHCHD2, GAPDH, and CD24 was strongly correlated with the malignant evolution of NSCLC at the single-cell data level. These genes have been validated via histological data. The risk score model that we established (represented by eight genes) was ultimately validated with GEO data.

CONCLUSION

In summary, our study contributes to the identification of temporal heterogeneous biomarkers in NSCLC, thus offering insights into disease progression mechanisms and potential therapeutic targets. The developed workflow demonstrates promise for future applications in clinical practice.

Mutation and cell state compatibility is required and targetable in Ph+ acute lymphoblastic leukemia minimal residual disease

Efforts to cure BCR::ABL1 B cell acute lymphoblastic leukemia (Ph+ ALL) solely through inhibition of ABL1 kinase activity have thus far been insufficient despite the availability of tyrosine kinase inhibitors (TKIs) with broad activity against resistance mutants. The mechanisms that drive persistence within minimal residual disease (MRD) remain poorly understood and therefore untargeted. Utilizing 13 patient-derived xenograft (PDX) models and clinical trial specimens of Ph+ ALL, we examined how genetic and transcriptional features co-evolve to drive progression during prolonged TKI response. Our work reveals a landscape of cooperative mutational and transcriptional escape mechanisms that differ from those causing resistance to first generation TKIs. By analyzing MRD during remission, we show that the same resistance mutation can either increase or decrease cellular fitness depending on transcriptional state. We further demonstrate that directly targeting transcriptional state-associated vulnerabilities at MRD can overcome BCR::ABL1 independence, suggesting a new paradigm for rationally eradicating MRD prior to relapse. Finally, we illustrate how cell mass measurements of leukemia cells can be used to rapidly monitor dominant transcriptional features of Ph+ ALL to help rationally guide therapeutic selection from low-input samples.

CD63 + tumor-associated macrophages drive the progression of hepatocellular carcinoma through the induction of epithelial-mesenchymal transition and lipid reprogramming

BACKGROUND

Tumor-associated macrophages (TAMs) constitute a substantial part of human hepatocellular carcinoma (HCC). The present study was devised to explore TAM diversity and their roles in HCC progression.

METHODS

Through the integration of multiple 10 × single-cell transcriptomic data derived from HCC samples and the use of consensus nonnegative matrix factorization (an unsupervised clustering algorithm), TAM molecular subtypes and expression programs were evaluated in detail. The roles played by these TAM subtypes in HCC were further probed through pseudotime, enrichment, and intercellular communication analyses. Lastly, vitro experiments were performed to validate the relationship between CD63, which is an inflammatory TAM expression program marker, and tumor cell lines.

RESULTS

We found that the inflammatory expression program in TAMs had a more obvious interaction with HCC cells, and CD63, as a marker gene of the inflammatory expression program, was associated with poor prognosis of HCC patients. Both bulk RNA-seq and vitro experiments confirmed that higher TAM CD63 expression was associated with the growth of HCC cells as well as their epithelial-mesenchymal transition, metastasis, invasion, and the reprogramming of lipid metabolism.

CONCLUSIONS

These analyses revealed that the TAM inflammatory expression program in HCC is closely associated with malignant tumor cells, with the hub gene CD63 thus representing an ideal target for therapeutic intervention in this cancer type.

Patterning and folding of intestinal villi by active mesenchymal dewetting

Tissue folds are structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, finger-like protrusions that enable nutrient absorption. However, the molecular and mechanical processes driving villus morphogenesis remain unclear. Here, we identify an active mechanical mechanism that simultaneously patterns and folds the intestinal epithelium to initiate villus formation. At the cellular level, we find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. This symmetry-breaking process requires altered cell and extracellular matrix interactions that are enabled by matrix metalloproteinase-mediated tissue fluidization. Computational models, together with in vitro and in vivo experiments, revealed that these cellular features manifest at the tissue level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active dewetting of a thin liquid film.

The Exonuclease TREX1 Constitutes an Innate Immune Checkpoint Limiting cGAS/STING-Mediated Antitumor Immunity

The DNA exonuclease three-prime repair exonuclease 1 (TREX1) is critical for preventing autoimmunity in mice and humans by degrading endogenous cytosolic DNA, which otherwise triggers activation of the innate cGAS/STING pathway leading to the production of type I IFNs. As tumor cells are prone to aberrant cytosolic DNA accumulation, we hypothesized that they are critically dependent on TREX1 activity to limit their immunogenicity. Here, we show that in tumor cells, TREX1 restricts spontaneous activation of the cGAS/STING pathway, and the subsequent induction of a type I IFN response. As a result, TREX1 deficiency compromised in vivo tumor growth in mice. This delay in tumor growth depended on a functional immune system, systemic type I IFN signaling, and tumor-intrinsic cGAS expression. Mechanistically, we show that tumor TREX1 loss drove activation of CD8+ T cells and NK cells, prevented CD8+ T-cell exhaustion, and remodeled an immunosuppressive myeloid compartment. Consequently, TREX1 deficiency combined with T-cell-directed immune checkpoint blockade. Collectively, we conclude that TREX1 is essential to limit tumor immunogenicity, and that targeting this innate immune checkpoint remodels the tumor microenvironment and enhances antitumor immunity by itself and in combination with T-cell-targeted therapies. See related article by Toufektchan et al., p. 673.

Identifying gene expression programs in single-cell RNA-seq data using linear correlation explanation

OBJECTIVE

Gene expression analysis through single-cell RNA sequencing (scRNA-seq) has revolutionized our understanding of gene regulation in diverse cell types, tissues, and organisms. While existing methods primarily focus on identifying cell type-specific gene expression programs (GEPs), the characterization of GEPs associated with biological processes and stimuli responses remains limited. In this study, we aim to infer biologically meaningful GEPs that are associated with both cellular phenotypes and activity programs directly from scRNA-seq data.

METHODS

We applied linear CorEx, a machine-learning-based approach, to infer GEPs by grouping genes based on total correlation optimization function in simulated and real-world scRNA-seq datasets. Additionally, we utilized a transfer learning approach to project CorEx-inferred GEPs to other scRNA-seq datasets.

RESULTS

By leveraging total correlation optimization, linear CorEx groups genes and demonstrates superior performance in identifying cell types and activity programs compared to similar methods using simulated data. Furthermore, we apply this same approach to real-world scRNA-seq data from the mouse dentate gyrus and embryonic colon development, uncovering biologically relevant GEPs related to cell types, developmental ages, and cell cycle programs. We also demonstrate the potential for transfer learning by evaluating similar datasets, showcasing the cross-species sensitivity of linear CorEx.

CONCLUSION

Our findings validate linear CorEx as a valuable tool for comprehensively analyzing complex signals in scRNA-seq data, leading to deeper insights into gene expression dynamics, cellular heterogeneity, and regulatory mechanisms.

The molecular basis of scale development highlighted by a single-cell atlas of Bicyclus anynana butterfly pupal forewings

Butterfly wings display a diversity of cell types, including large polyploid scale cells, yet the molecular basis of such diversity is poorly understood. To explore scale cell diversity at a transcriptomic level, we employ single-cell RNA sequencing of ∼5,200 large cells (>6 μm) from 22.5- to 25-h male pupal forewings of the butterfly Bicyclus anynana. Using unsupervised clustering, followed by in situ hybridization, immunofluorescence, and CRISPR-Cas9 editing of candidate genes, we annotate various cell types on the wing. We identify genes marking non-innervated scale cells, pheromone-producing glandular cells, and innervated sensory cell types. We show that senseless, a zinc-finger transcription factor, and HR38, a hormone receptor, determine the identity, size, and color of different scale cell types and are important regulators of scale cell differentiation. This dataset and the identification of various wing cell-type markers provide a foundation to compare and explore scale cell-type diversification across arthropod species.

Dissecting tumor transcriptional heterogeneity from single-cell RNA-seq data by generalized binary covariance decomposition

Profiling tumors with single-cell RNA sequencing (scRNA-seq) has the potential to identify recurrent patterns of transcription variation related to cancer progression, and produce new therapeutically relevant insights. However, the presence of strong inter-tumor heterogeneity often obscures more subtle patterns that are shared across tumors, some of which may characterize clinically relevant disease subtypes. Here we introduce a new statistical method, generalized binary covariance decomposition (GBCD), to address this problem. We show that GBCD can help decompose transcriptional heterogeneity into interpretable components - including patient-specific, dataset-specific and shared components relevant to disease subtypes - and that, in the presence of strong inter-tumor heterogeneity, it can produce more interpretable results than existing methods. Applied to data from three studies on pancreatic cancer adenocarcinoma (PDAC), GBCD produces a refined characterization of existing tumor subtypes (e.g., classical vs. basal), and identifies a new gene expression program (GEP) that is prognostic of poor survival independent of established prognostic factors such as tumor stage and subtype. The new GEP is enriched for genes involved in a variety of stress responses, and suggests a potentially important role for the integrated stress response in PDAC development and prognosis.

Coordinated immune dysregulation in juvenile dermatomyositis revealed by single-cell genomics

Juvenile dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I IFN response and autoantibodies. Treatment options are limited due to an incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of patients with JDM at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment toward an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I IFN response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK cells, CD8+ T cells, and γδ T cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell Spatio- and Neighborhood-informed Annotation and Patterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

Unveiling the molecular complexity of proliferative diabetic retinopathy through scRNA-seq, AlphaFold 2, and machine learning

Background

Proliferative diabetic retinopathy (PDR), a major cause of blindness, is characterized by complex pathogenesis. This study integrates single-cell RNA sequencing (scRNA-seq), Non-negative Matrix Factorization (NMF), machine learning, and AlphaFold 2 methods to explore the molecular level of PDR.

Methods

We analyzed scRNA-seq data from PDR patients and healthy controls to identify distinct cellular subtypes and gene expression patterns. NMF was used to define specific transcriptional programs in PDR. The oxidative stress-related genes (ORGs) identified within Meta-Program 1 were utilized to construct a predictive model using twelve machine learning algorithms. Furthermore, we employed AlphaFold 2 for the prediction of protein structures, complementing this with molecular docking to validate the structural foundation of potential therapeutic targets. We also analyzed protein-protein interaction (PPI) networks and the interplay among key ORGs.

Results

Our scRNA-seq analysis revealed five major cell types and 14 subcell types in PDR patients, with significant differences in gene expression compared to those in controls. We identified three key meta-programs underscoring the role of microglia in the pathogenesis of PDR. Three critical ORGs (ALKBH1, PSIP1, and ATP13A2) were identified, with the best-performing predictive model demonstrating high accuracy (AUC of 0.989 in the training cohort and 0.833 in the validation cohort). Moreover, AlphaFold 2 predictions combined with molecular docking revealed that resveratrol has a strong affinity for ALKBH1, indicating its potential as a targeted therapeutic agent. PPI network analysis, revealed a complex network of interactions among the hub ORGs and other genes, suggesting a collective role in PDR pathogenesis.

Conclusion

This study provides insights into the cellular and molecular aspects of PDR, identifying potential biomarkers and therapeutic targets using advanced technological approaches.

Durable lymph-node expansion is associated with the efficacy of therapeutic vaccination

Following immunization, lymph nodes dynamically expand and contract. The mechanical and cellular changes enabling the early-stage expansion of lymph nodes have been characterized, yet the durability of such responses and their implications for adaptive immunity and vaccine efficacy are unknown. Here, by leveraging high-frequency ultrasound imaging of the lymph nodes of mice, we report more potent and persistent lymph-node expansion for animals immunized with a mesoporous silica vaccine incorporating a model antigen than for animals given bolus immunization or standard vaccine formulations such as alum, and that durable and robust lymph-node expansion was associated with vaccine efficacy and adaptive immunity for 100 days post-vaccination in a mouse model of melanoma. Immunization altered the mechanical and extracellular-matrix properties of the lymph nodes, drove antigen-dependent proliferation of immune and stromal cells, and altered the transcriptional features of dendritic cells and inflammatory monocytes. Strategies that robustly maintain lymph-node expansion may result in enhanced vaccination outcomes.