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

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.

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.

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.

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.

Reproducible single cell annotation of programs underlying T-cell subsets, activation states, and functions

T-cells recognize antigens and induce specialized gene expression programs (GEPs) enabling functions including proliferation, cytotoxicity, and cytokine production. Traditionally, different classes of helper T-cells express mutually exclusive responses - for example, Th1, Th2, and Th17 programs. However, new single-cell RNA sequencing (scRNA-Seq) experiments have revealed a continuum of T-cell states without discrete clusters corresponding to these subsets, implying the need for new analytical frameworks. Here, we advance the characterization of T-cells with T-CellAnnoTator (TCAT), a pipeline that simultaneously quantifies pre-defined GEPs capturing activation states and cellular subsets. From 1,700,000 T-cells from 700 individuals across 38 tissues and five diverse disease contexts, we discover 46 reproducible GEPs reflecting the known core functions of T-cells including proliferation, cytotoxicity, exhaustion, and T helper effector states. We experimentally characterize several novel activation programs and apply TCAT to describe T-cell activation and exhaustion in Covid-19 and cancer, providing insight into T-cell function in these diseases.

Comparative single cell analysis reveals complex patterns of cell type and cell signaling innovations at the fetal-maternal interface

The decidual-placental interface is one of the most diverse and rapidly evolving tissues in mammals. Its origin as a chimeric fetal-maternal tissue poses a unique evolutionary puzzle. We present single-cell RNA sequencing atlases from the fetal-maternal interfaces of the opossum, a marsupial, the Malagasy common tenrec, an afrotherian with primitive reproductive features, and mouse, guinea pig, and human. Invasive trophoblast shares a common transcriptomic signature across eutherians, which we argue represents a cell type family that radiated following the evolution of hemochorial placentation. We find evidence that the eutherian decidual stromal cell evolved stepwise from a predecidual state retained in Tenrec , followed by a second decidual cell type originating in Boreoeutheria with endocrine characteristics. We reconstruct ligand-receptor signaling to test evolutionary hypotheses at scale. Novel trophoblast and decidual cell types display strong integration into signaling networks compared to other cells. Additionally, we find consistent disambiguation between fetal and maternal signaling. Using phylogenetic analysis, we infer the cell-cell signaling network of the Placental common ancestor, and identify increased rates of signaling evolution in Euarchontoglires. Together, our findings reveal novel cell type identities and cell signaling dynamics at the mammalian fetal-maternal interface.

Gene networks and systems biology in Alzheimer's disease: Insights from multi-omics approaches

Despite numerous studies in the field of dementia and Alzheimer's disease (AD), a comprehensive understanding of this devastating disease remains elusive. Bulk transcriptomics have provided insights into the underlying genetic factors at a high level. Subsequent technological advancements have focused on single-cell omics, encompassing techniques such as single-cell RNA sequencing and epigenomics, enabling the capture of RNA transcripts and chromatin states at a single cell or nucleus resolution. Furthermore, the emergence of spatial omics has allowed the study of gene responses in the vicinity of amyloid beta plaques or across various brain regions. With the vast amount of data generated, utilizing gene regulatory networks to comprehensively study this disease has become essential. This review delves into some techniques employed in the field of AD, explores the discoveries made using these techniques, and provides insights into the future of the field.

Early Immune Remodeling Steers Clinical Response to First-Line Chemoimmunotherapy in Advanced Gastric Cancer

Adding anti-programmed cell death protein 1 (anti-PD-1) to 5-fluorouracil (5-FU)/platinum improves survival in some advanced gastroesophageal adenocarcinomas (GEA). To understand the effects of chemotherapy and immunotherapy, we conducted a phase II first-line trial (n = 47) sequentially adding pembrolizumab to 5-FU/platinum in advanced GEA. Using serial biopsy of the primary tumor at baseline, after one cycle of 5-FU/platinum, and after the addition of pembrolizumab, we transcriptionally profiled 358,067 single cells to identify evolving multicellular tumor microenvironment (TME) networks. Chemotherapy induced early on-treatment multicellular hubs with tumor-reactive T-cell and M1-like macrophage interactions in slow progressors. Faster progression featured increased MUC5A and MSLN containing treatment resistance programs in tumor cells and M2-like macrophages with immunosuppressive stromal interactions. After pembrolizumab, we observed increased CD8 T-cell infiltration and development of an immunity hub involving tumor-reactive CXCL13 T-cell program and epithelial interferon-stimulated gene programs. Strategies to drive increases in antitumor immune hub formation could expand the portion of patients benefiting from anti-PD-1 approaches.

SIGNIFICANCE

The benefit of 5-FU/platinum with anti-PD-1 in first-line advanced gastric cancer is limited to patient subgroups. Using a trial with sequential anti-PD-1, we show coordinated induction of multicellular TME hubs informs the ability of anti-PD-1 to potentiate T cell-driven responses. Differential TME hub development highlights features that underlie clinical outcomes. This article is featured in Selected Articles from This Issue, p. 695.

Characterizing cell type specific transcriptional differences between the living and postmortem human brain

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

The role of quiescent thymic progenitors in TAL/LMO2-induced T-ALL chemotolerance

Relapse in T-cell acute lymphoblastic leukemia (T-ALL) may signify the persistence of leukemia-initiating cells (L-ICs). Ectopic TAL1/LMO expression defines the largest subset of T-ALL, but its role in leukemic transformation and its impact on relapse-driving L-ICs remain poorly understood. In TAL1/LMO mouse models, double negative-3 (DN3; CD4-CD8-CD25+CD44-) thymic progenitors harbored L-ICs. However, only a subset of DN3 leukemic cells exhibited L-IC activity, and studies linking L-ICs and chemotolerance are needed. To investigate L-IC heterogeneity, we used mouse models and applied single-cell RNA-sequencing and nucleosome labeling techniques in vivo. We identified a DN3 subpopulation with a cell cycle-restricted profile and heightened TAL1/LMO2 activity, that expressed genes associated with stemness and quiescence. This dormant DN3 subset progressively expanded throughout leukemogenesis, displaying intrinsic chemotolerance and enrichment in genes linked to minimal residual disease. Examination of TAL/LMO patient samples revealed a similar pattern in CD7+CD1a- thymic progenitors, previously recognized for their L-IC activity, demonstrating cell cycle restriction and chemotolerance. Our findings substantiate the emergence of dormant, chemotolerant L-ICs during leukemogenesis, and demonstrate that Tal1 and Lmo2 cooperate to promote DN3 quiescence during the transformation process. This study provides a deeper understanding of TAL1/LMO-induced T-ALL and its clinical implications in therapy failure.

Identifying somatic fingerprints of cancers defined by germline and environmental risk factors

Numerous studies over the past generation have identified germline variants that increase specific cancer risks. Simultaneously, a revolution in sequencing technology has permitted high-throughput annotations of somatic genomes characterizing individual tumors. However, examining the relationship between germline variants and somatic alteration patterns is hugely challenged by the large numbers of variants in a typical tumor, the rarity of most individual variants, and the heterogeneity of tumor somatic fingerprints. In this article, we propose statistical methodology that frames the investigation of germline-somatic relationships in an interpretable manner. The method uses meta-features embodying biological contexts of individual somatic alterations to implicitly group rare mutations. Our team has used this technique previously through a multilevel regression model to diagnose with high accuracy tumor site of origin. Herein, we further leverage topic models from computational linguistics to achieve interpretable lower-dimensional embeddings of the meta-features. We demonstrate how the method can identify distinctive somatic profiles linked to specific germline variants or environmental risk factors. We illustrate the method using The Cancer Genome Atlas whole-exome sequencing data to characterize somatic tumor fingerprints in breast cancer patients with germline BRCA1/2 mutations and in head and neck cancer patients exposed to human papillomavirus.

An integrated single-nucleus and spatial transcriptomics atlas reveals the molecular landscape of the human hippocampus

The hippocampus contains many unique cell types, which serve the structure's specialized functions, including learning, memory and cognition. These cells have distinct spatial topography, morphology, physiology, and connectivity, highlighting the need for transcriptome-wide profiling strategies that retain cytoarchitectural organization. Here, we generated spatially-resolved transcriptomics (SRT) and single-nucleus RNA-sequencing (snRNA-seq) data from adjacent tissue sections of the anterior human hippocampus across ten adult neurotypical donors. We defined molecular profiles for hippocampal cell types and spatial domains. Using non-negative matrix factorization and transfer learning, we integrated these data to define gene expression patterns within the snRNA-seq data and infer the expression of these patterns in the SRT data. With this approach, we leveraged existing rodent datasets that feature information on circuit connectivity and neural activity induction to make predictions about axonal projection targets and likelihood of ensemble recruitment in spatially-defined cellular populations of the human hippocampus. Finally, we integrated genome-wide association studies with transcriptomic data to identify enrichment of genetic components for neurodevelopmental, neuropsychiatric, and neurodegenerative disorders across cell types, spatial domains, and gene expression patterns of the human hippocampus. To make this comprehensive molecular atlas accessible to the scientific community, both raw and processed data are freely available, including through interactive web applications.

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

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. Here, we generated single-cell RNA-seq maps of neuroblastoma cell lines, patient-derived xenograft models (PDX), and a genetically engineered mouse model (GEMM). We developed 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 confirmed 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 was maintained in PDX and could be chemotherapy-induced in our GEMM after only 24 hours, suggesting an uncharacterized therapy-escape mechanism. 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.

A map of signaling responses in the human airway epithelium

Receptor-mediated signaling plays a central role in tissue regeneration, and it is dysregulated in disease. Here, we build a signaling-response map for a model regenerative human tissue: the airway epithelium. We analyzed the effect of 17 receptor-mediated signaling pathways on organotypic cultures to determine changes in abundance and phenotype of epithelial cell types. This map recapitulates the gamut of known airway epithelial signaling responses to these pathways. It defines convergent states induced by multiple ligands and diverse, ligand-specific responses in basal cell and secretory cell metaplasia. We show that loss of canonical differentiation induced by multiple pathways is associated with cell-cycle arrest, but that arrest is not sufficient to block differentiation. Using the signaling-response map, we show that a TGFB1-mediated response underlies specific aberrant cells found in multiple lung diseases and identify interferon responses in COVID-19 patient samples. Thus, we offer a framework enabling systematic evaluation of tissue signaling responses. A record of this paper's transparent peer review process is included in the supplemental information.

q-Diffusion leverages the full dimensionality of gene coexpression in single-cell transcriptomics

Unlocking the full dimensionality of single-cell RNA sequencing data (scRNAseq) is the next frontier to a richer, fuller understanding of cell biology. We introduce q-diffusion, a framework for capturing the coexpression structure of an entire library of genes, improving on state-of-the-art analysis tools. The method is demonstrated via three case studies. In the first, q-diffusion helps gain statistical significance for differential effects on patient outcomes when analyzing the CALGB/SWOG 80405 randomized phase III clinical trial, suggesting precision guidance for the treatment of metastatic colorectal cancer. Secondly, q-diffusion is benchmarked against existing scRNAseq classification methods using an in vitro PBMC dataset, in which the proposed method discriminates IFN-γ stimulation more accurately. The same case study demonstrates improvements in unsupervised cell clustering with the recent Tabula Sapiens human atlas. Finally, a local distributional segmentation approach for spatial scRNAseq, driven by q-diffusion, yields interpretable structures of human cortical tissue.

Topical hidden genome: discovering latent cancer mutational topics using a Bayesian multilevel context-learning approach

Inferring the cancer-type specificities of ultra-rare, genome-wide somatic mutations is an open problem. Traditional statistical methods cannot handle such data due to their ultra-high dimensionality and extreme data sparsity. To harness information in rare mutations, we have recently proposed a formal multilevel multilogistic "hidden genome" model. Through its hierarchical layers, the model condenses information in ultra-rare mutations through meta-features embodying mutation contexts to characterize cancer types. Consistent, scalable point estimation of the model can incorporate 10s of millions of variants across thousands of tumors and permit impressive prediction and attribution. However, principled statistical inference is infeasible due to the volume, correlation, and noninterpretability of mutation contexts. In this paper, we propose a novel framework that leverages topic models from computational linguistics to effectuate dimension reduction of mutation contexts producing interpretable, decorrelated meta-feature topics. We propose an efficient MCMC algorithm for implementation that permits rigorous full Bayesian inference at a scale that is orders of magnitude beyond the capability of existing out-of-the-box inferential high-dimensional multi-class regression methods and software. Applying our model to the Pan Cancer Analysis of Whole Genomes dataset reveals interesting biological insights including somatic mutational topics associated with UV exposure in skin cancer, aging in colorectal cancer, and strong influence of epigenome organization in liver cancer. Under cross-validation, our model demonstrates highly competitive predictive performance against blackbox methods of random forest and deep learning.

Interpretable Spatial Gradient Analysis for Spatial Transcriptomics Data

Cellular anatomy and signaling vary across niches, which can induce gradated gene expressions in subpopulations of cells. Such spatial transcriptomic gradient (STG) makes a significant source of intratumor heterogeneity and can influence tumor invasion, progression, and response to treatment. Here we report Local Spatial Gradient Inference (LSGI), a computational framework that systematically identifies spatial locations with prominent, interpretable STGs from spatial transcriptomic (ST) data. To achieve so, LSGI scrutinizes each sliding window employing non-negative matrix factorization (NMF) combined with linear regression. With LSGI, we demonstrated the identification of spatially proximal yet opposite directed pathway gradients in a glioblastoma dataset. We further applied LSGI to 87 tumor ST datasets reported from nine published studies and identified both pan-cancer and tumor-type specific pathways with gradated expression patterns, such as epithelial mesenchymal transition, MHC complex, and hypoxia. The local gradients were further categorized according to their association to tumor-TME (tumor microenvironment) interface, highlighting the pathways related to spatial transcriptional intratumoral heterogeneity. We conclude that LSGI enables highly interpretable STG analysis which can reveal novel insights in tumor biology from the increasingly reported tumor ST datasets.

Highly sensitive spatial transcriptomics using FISHnCHIPs of multiple co-expressed genes

High-dimensional, spatially resolved analysis of intact tissue samples promises to transform biomedical research and diagnostics, but existing spatial omics technologies are costly and labor-intensive. We present Fluorescence In Situ Hybridization of Cellular HeterogeneIty and gene expression Programs (FISHnCHIPs) for highly sensitive in situ profiling of cell types and gene expression programs. FISHnCHIPs achieves this by simultaneously imaging ~2-35 co-expressed genes (clustered into modules) that are spatially co-localized in tissues, resulting in similar spatial information as single-gene Fluorescence In Situ Hybridization (FISH), but with ~2-20-fold higher sensitivity. Using FISHnCHIPs, we image up to 53 modules from the mouse kidney and mouse brain, and demonstrate high-speed, large field-of-view profiling of a whole tissue section. FISHnCHIPs also reveals spatially restricted localizations of cancer-associated fibroblasts in a human colorectal cancer biopsy. Overall, FISHnCHIPs enables fast, robust, and scalable cell typing of tissues with normal physiology or undergoing pathogenesis.

A Spatial Multi-Modal Dissection of Host-Microbiome Interactions within the Colitis Tissue Microenvironment

The intricate and dynamic interactions between the host immune system and its microbiome constituents undergo dynamic shifts in response to perturbations to the intestinal tissue environment. Our ability to study these events on the systems level is significantly limited by in situ approaches capable of generating simultaneous insights from both host and microbial communities. Here, we introduce Microbiome Cartography (MicroCart), a framework for simultaneous in situ probing of host features and its microbiome across multiple spatial modalities. We demonstrate MicroCart by comprehensively investigating the alterations in both gut host and microbiome components in a murine model of colitis by coupling MicroCart with spatial proteomics, transcriptomics, and glycomics platforms. Our findings reveal a global but systematic transformation in tissue immune responses, encompassing tissue-level remodeling in response to host immune and epithelial cell state perturbations, and bacterial population shifts, localized inflammatory responses, and metabolic process alterations during colitis. MicroCart enables a deep investigation of the intricate interplay between the host tissue and its microbiome with spatial multiomics.

Spatial transcriptomics data and analytical methods: An updated perspective

Spatial transcriptomics (ST) is a newly emerging field that integrates high-resolution imaging and transcriptomic data to enable the high-throughput analysis of the spatial localization of transcripts in diverse biological systems. The rapid progress in this field necessitates the development of innovative computational methods to effectively tackle the distinct challenges posed by the analysis of ST data. These platforms, integrating AI techniques, offer a promising avenue for understanding disease mechanisms and expediting drug discovery. Despite significant advances in the development of ST data analysis techniques, there is an ongoing need to enhance these models for increased biological relevance. In this review, we briefly discuss the ST-related databases and current deep-learning-based models for spatial transcriptome data analyses and highlight their roles and future perspectives in biomedical applications.

Single cell RNA-seq: a novel tool to unravel virus-host interplay

Over the last decade, single cell RNA sequencing (scRNA-seq) technology has caught the momentum of being a vital revolutionary tool to unfold cellular heterogeneity by high resolution assessment. It evades the inadequacies of conventional sequencing technology which was able to detect only average expression level among cell populations. In the era of twenty-first century, several epidemic and pandemic viruses have emerged. Being an intracellular entity, viruses totally rely on host. Complex virus-host dynamics result when the virus tend to obtain factors from host cell required for its replication and establishment of infection. As a prevailing tool, scRNA-seq is able to understand virus-host interplay by comprehensive transcriptome profiling. Because of technological and methodological advancement, this technology is capable to recognize viral genome and host cell response heterogeneity. Further development in analytical methods with multiomics approach and increased availability of accessible scRNA-seq datasets will improve the understanding of viral pathogenesis that can be helpful for development of novel antiviral therapeutic strategies.

A concerted neuron-astrocyte program declines in ageing and schizophrenia

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Single-cell transcriptome landscape of circulating CD4+ T cell populations in autoimmune diseases

CD4+ T cells are key mediators of various autoimmune diseases; however, their role in disease progression remains unclear due to cellular heterogeneity. Here, we evaluated CD4+ T cell subpopulations using decomposition-based transcriptome characterization and canonical clustering strategies. This approach identified 12 independent gene programs governing whole CD4+ T cell heterogeneity, which can explain the ambiguity of canonical clustering. In addition, we performed a meta-analysis using public single-cell datasets of over 1.8 million peripheral CD4+ T cells from 953 individuals by projecting cells onto the reference and cataloging cell frequency and qualitative alterations of the populations in 20 diseases. The analyses revealed that the 12 transcriptional programs were useful in characterizing each autoimmune disease and predicting its clinical status. Moreover, genetic variants associated with autoimmune diseases showed disease-specific enrichment within the 12 gene programs. The results collectively provide a landscape of single-cell transcriptomes of CD4+ T cell subpopulations involved in autoimmune disease.

Residual ANTXR1+ myofibroblasts after chemotherapy inhibit anti-tumor immunity via YAP1 signaling pathway

Although cancer-associated fibroblast (CAF) heterogeneity is well-established, the impact of chemotherapy on CAF populations remains poorly understood. Here we address this question in high-grade serous ovarian cancer (HGSOC), in which we previously identified 4 CAF populations. While the global content in stroma increases in HGSOC after chemotherapy, the proportion of FAP+ CAF (also called CAF-S1) decreases. Still, maintenance of high residual CAF-S1 content after chemotherapy is associated with reduced CD8+ T lymphocyte density and poor patient prognosis, emphasizing the importance of CAF-S1 reduction upon treatment. Single cell analysis, spatial transcriptomics and immunohistochemistry reveal that the content in the ECM-producing ANTXR1+ CAF-S1 cluster (ECM-myCAF) is the most affected by chemotherapy. Moreover, functional assays demonstrate that ECM-myCAF isolated from HGSOC reduce CD8+ T-cell cytotoxicity through a Yes Associated Protein 1 (YAP1)-dependent mechanism. Thus, efficient inhibition after treatment of YAP1-signaling pathway in the ECM-myCAF cluster could enhance CD8+ T-cell cytotoxicity. Altogether, these data pave the way for therapy targeting YAP1 in ECM-myCAF in HGSOC.

Single-cell RNA sequencing of cervical exfoliated cells reveals potential biomarkers and cellular pathogenesis in cervical carcinogenesis

Cervical cancer (CC) is a common gynecological malignancy. Despite the current screening methods have been proved effectively and significantly decreased CC morbidity and mortality, deficiencies still exist. Single-cell RNA sequencing (scRNA-seq) approach can identify the complex and rare cell populations at single-cell resolution. By scRNA-seq, the heterogeneity of tumor microenvironment across cervical carcinogenesis has been mapped and described. Whether these alterations could be detected and applied to CC screening is unclear. Herein, we performed scRNA-seq of 56,173 cervical exfoliated cells from 15 samples, including normal cervix, low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), and malignancy. The present study delineated the alteration of immune and epithelial cells derived during the cervical lesion progression. A subset of lipid-associated macrophage was identified as a tumor-promoting element and could serve as a biomarker for predicting the progression of LSIL into HSIL, which was then verified by immunofluorescence. Furthermore, cell-cell communication analysis indicated the SPP1-CD44 axis might exhibit a protumor interaction between epithelial cell and macrophage. In this study, we investigated the cervical multicellular ecosystem in cervical carcinogenesis and identified potential biomarkers for early detection.

Convergence of coronary artery disease genes onto endothelial cell programs

Linking variants from genome-wide association studies (GWAS) to underlying mechanisms of disease remains a challenge1-3. For some diseases, a successful strategy has been to look for cases in which multiple GWAS loci contain genes that act in the same biological pathway1-6. However, our knowledge of which genes act in which pathways is incomplete, particularly for cell-type-specific pathways or understudied genes. Here we introduce a method to connect GWAS variants to functions. This method links variants to genes using epigenomics data, links genes to pathways de novo using Perturb-seq and integrates these data to identify convergence of GWAS loci onto pathways. We apply this approach to study the role of endothelial cells in genetic risk for coronary artery disease (CAD), and discover 43 CAD GWAS signals that converge on the cerebral cavernous malformation (CCM) signalling pathway. Two regulators of this pathway, CCM2 and TLNRD1, are each linked to a CAD risk variant, regulate other CAD risk genes and affect atheroprotective processes in endothelial cells. These results suggest a model whereby CAD risk is driven in part by the convergence of causal genes onto a particular transcriptional pathway in endothelial cells. They highlight shared genes between common and rare vascular diseases (CAD and CCM), and identify TLNRD1 as a new, previously uncharacterized member of the CCM signalling pathway. This approach will be widely useful for linking variants to functions for other common polygenic diseases.

Rat microbial biogeography and age-dependent lactic acid bacteria in healthy lungs

The laboratory rat emerges as a useful tool for studying the interaction between the host and its microbiome. To advance principles relevant to the human microbiome, we systematically investigated and defined the multitissue microbial biogeography of healthy Fischer 344 rats across their lifespan. Microbial community profiling data were extracted and integrated with host transcriptomic data from the Sequencing Quality Control consortium. Unsupervised machine learning, correlation, taxonomic diversity and abundance analyses were performed to determine and characterize the rat microbial biogeography and identify four intertissue microbial heterogeneity patterns (P1-P4). We found that the 11 body habitats harbored a greater diversity of microbes than previously suspected. Lactic acid bacteria (LAB) abundance progressively declined in lungs from breastfed newborn to adolescence/adult, and was below detectable levels in elderly rats. Bioinformatics analyses indicate that the abundance of LAB may be modulated by the lung-immune axis. The presence and levels of LAB in lungs were further evaluated by PCR in two validation datasets. The lung, testes, thymus, kidney, adrenal and muscle niches were found to have age-dependent alterations in microbial abundance. The 357 microbial signatures were positively correlated with host genes in cell proliferation (P1), DNA damage repair (P2) and DNA transcription (P3). Our study established a link between the metabolic properties of LAB with lung microbiota maturation and development. Breastfeeding and environmental exposure influence microbiome composition and host health and longevity. The inferred rat microbial biogeography and pattern-specific microbial signatures could be useful for microbiome therapeutic approaches to human health and life quality enhancement.