Interpretable modeling of time-resolved single-cell gene-protein expression with CrossmodalNet.

Background: Cancer stem cells (CSCs) play a key role in tumor initiation, progression, and resistance to conventional therapies. These cells possess the ability to self-renew and differentiate, contributing to tumor heterogeneity and recurrence. Preclinical data indicate that targeting cancer stem cell-specific pathways could lead to more effective treatments and prevent relapse, thereby improving outcomes. However, there is limited clinical data to support this, in part due to the challenges of measuring CSCs in the clinical setting. Single-cell RNA sequencing (scRNAseq) data in combination with stemness gene expression signatures were investigated as a novel approach to detect and assess changes in CSC numbers in response to treatment. Methods: Single-cell RNAseq datasets from two clinical trials were interrogated – OPPORTUNE, a window-of-opportunity trial evaluating anastrozole vs anastrozole plus pictilisib (PI3K inhibitor) in 75 patients with ER+ HER2- early breast cancer, and FELINE, a neoadjuvant trial which compared letrozole plus placebo with letrozole plus ribociclib (CDK4/6 inhibitor) in 120 patients with ER+ HER2- early breast cancer. The primary endpoint for OPPORTUNE was the inhibition of tumor cell proliferation as measured by Ki67. The primary endpoint for FELINE was the rate of preoperative endocrine prognostic index (PEPI) score 0 after neoadjuvant endocrine therapy. Thirty-four patients from the FELINE study had available scRNAseq data for analysis; this number was 62 patients for the OPPORTUNE study. Stemness gene expression signatures with published evidence of selectivity for breast CSCs were identified from the literature and their utility in detecting CSCs compared in the datasets. Subsequently, changes in CSC fraction with treatment was assessed. Results: Eight different stemness gene signatures were identified from the published literature. Among these, four (Kim_myc, benporath_es2, Bhattacharya_hESC, and Shats_consensus) exhibited selective expression in a minority of tumor cells, with the Shats consensus signature demonstrating the highest specificity to tumor cells over non-malignant cells. Utilizing a Gaussian mixture model, we estimated that approximately 4.5% of cells within a tumor are likely stem cells. In the OPPORTUNE study, higher stemness scores were observed in luminal B compared to luminal A tumors, but there was no association with PIK3CA mutation status. Changes to Ki67 were inversely associated with the CSC fraction – tumors with higher stemness scores were less likely to achieve a complete cell cycle arrest versus endocrine-sensitive tumors. Similarly, in the FELINE trial, non-responders in the letrozole arm showed a trend towards elevated stemness scores. Conclusion: The use of stemness gene expression signatures in scRNAseq data is a feasible method to assess changes in putative CSC fraction with treatment in breast cancer. Increased stemness scores were associated with more aggressive subtypes and resistance to treatment. Citation Format: Peter Hall,Alejandro Chibly, Peter Schmid, Sarah E Pinder, Arnie Purushotham, Alastair Thompson, Steven Gendreau. Evaluation of breast cancer stem cell gene expression signatures in single-cell RNA sequencing (scRNAseq) data from the OPPORTUNE and FELINE trials, and the association with treatment resistance [abstract]. In: Proceedings of the San Antonio Breast Cancer Symposium 2024; 2024 Dec 10-13; San Antonio, TX. Philadelphia (PA): AACR; Clin Cancer Res 2025;31(12 Suppl):Abstract nr PS13-01.

Single-cell RNA sequencing (scRNA-seq) data are widely used to study cancer cell states and their heterogeneity. However, the tumour microenvironment is usually a mixture of healthy and cancerous cells and it can be difficult to fully separate these two populations based on transcriptomics alone. If available, somatic single-nucleotide variants (SNVs) observed in the scRNA-seq data could be used to identify the cancer population and match that information with the single cells' expression profile. However, calling somatic SNVs in scRNA-seq data is a challenging task, as most variants seen in the short-read data are not somatic, but can instead be germline variants, RNA edits or transcription, sequencing, or processing errors. In addition, only variants present in actively transcribed regions for each individual cell will be seen in the data. To address these challenges, we develop CCLONE (Cancer Cell Labelling On Noisy Expression), an interpretable tool adapted to handle the uncertainty and sparsity of SNVs called from scRNA-seq data. CCLONE jointly identifies cancer clonal populations, and their associated variants. We apply CCLONE on two acute myeloid leukaemia datasets and one lung adenocarcinoma dataset and show that CCLONE captures both genetic clones and somatic events for multiple patients. These results show how CCLONE can be used to gather insight into the course of the disease and the origin of cancer cells in scRNA-seq data. Source code is available at github.com/HaghverdiLab/CCLONE.

Cell type prediction is one of the most challenging goals in single-cell RNA sequencing (scRNA-seq) data. Existing methods use unsupervised learning to identify signature genes in each cluster, followed by a literature survey to look up those genes for assigning cell types. However, finding potential marker genes in each cluster is cumbersome, which impedes the systematic analysis of single-cell RNA sequencing data. To address this challenge, we proposed a framework based on regularized multi-task learning (RMTL) that enables us to simultaneously learn the subpopulation associated with a particular cell type. Learning the structure of subpopulations is treated as a separate task in the multi-task learner. Regularization is used to modulate the multi-task model (e.g., W1, W2, … Wt) jointly, according to the specific prior. For validating our model, we trained it with reference data constructed from a single-cell RNA sequencing experiment and applied it to a query dataset. We also predicted completely independent data (the query dataset) from the reference data which are used for training. We have checked the efficacy of the proposed method by comparing it with other state-of-the-art techniques well known for cell type detection. Results revealed that the proposed method performed accurately in detecting the cell type in scRNA-seq data and thus can be utilized as a useful tool in the scRNA-seq pipeline.

Anchored-fusion enables targeted fusion search in bulk and single-cell RNA sequencing data

Transcriptomic analysis across species is increasingly used to reveal conserved gene regulations which implicate crucial regulators. Cross-species analysis of single-cell RNA sequencing (scRNA-seq) data provides new opportunities to identify the cellular and molecular conservations, especially for cell types and cell type-specific gene regulations. However, few methods have been developed to analyze cross-species scRNA-seq data to uncover both molecular and cellular conservations. Here, we built a tool called CACIMAR, which can perform cross-species analysis of cell identities, markers, regulations, and interactions using scRNA-seq profiles. Based on the weighted sum models of the conserved features, we developed different conservation scores to measure the conservation of cell types, regulatory networks, and intercellular interactions. Using publicly available scRNA-seq data on retinal regeneration in mice, zebrafish, and chick, we demonstrated four main functions of CACIMAR. First, CACIMAR allows to identify conserved cell types even in evolutionarily distant species. Second, the tool facilitates the identification of evolutionarily conserved or species-specific marker genes. Third, CACIMAR enables the identification of conserved intracellular regulations, including cell type-specific regulatory subnetworks and regulators. Lastly, CACIMAR provides a unique feature for identifying conserved intercellular interactions. Overall, CACIMAR facilitates the identification of evolutionarily conserved cell types, marker genes, intracellular regulations, and intercellular interactions, providing insights into the cellular and molecular mechanisms of species evolution.

Osteoporosis is a common disease, especially among the elderly. This study aimed to comprehensively examine the roles of immune microenvironment in osteoporosis pathogenesis. Expression profiles of GSE35959, GSE7158, and GSE13850 datasets were used to analyze differential expression and identify hub genes related to immune features. Based on the single-cell RNA sequencing (scRNA-seq) data of an osteoporosis patient, different cell types were classified and the relation between immune environment and osteoporosis was explored. Twelve hub genes significantly associated with immune features were selected and 11 subgroups were defined using scRNA-seq data. The expression of two hub genes (CDKN1A and TEFM) was greatly altered during the transformation from mesenchymal stem cells (MSCs) to osteoblasts. Chemokines and chemokine receptors were differentially enriched in different cell types. CXCL12 was high-expressed in MSCs. This study emphasized that immune microenvironment played a critical role in the pathogenesis of osteoporosis. Chemokines and chemokine receptors can modify cell development and affect the interactions among different cell types, leading to unbalanced bone remodeling.

Analysis of single-cell RNA sequencing (scRNA-seq) data of immune cells from the tumor microenvironment (TME) may identify tumor progression biomarkers. This study was designed to investigate the prognostic value of differentially expressed genes (DEGs) in intrahepatic cholangiocarcinoma (ICC) using scRNA-seq. We downloaded the scRNA-seq data of 33,991 cell samples, including 17,090 ICC cell samples and 16,901 ICC adjacent tissue cell samples regarded as normal cells. scRNA-seq data were processed and classified into 20 clusters. The immune cell clusters were extracted and processed again in the same way, and each type of immune cells was divided into several subclusters. In total, 337 marker genes of macrophages and 427 marker genes of B cells were identified by comparing ICC subclusters with normal subclusters. Finally, 659 DEGs were obtained by merging B cell and macrophage marker genes. ICC sample clinical information and gene expression data were downloaded. A nine-prognosis-related-gene (PRG) signature was established by analyzing the correlation between DEGs and overall survival in ICC. The robustness and validity of the signature were verified. Functional enrichment analysis revealed that the nine PRGs were mainly involved in tumor immune mechanisms. In conclusion, we established a PRG signature based on scRNA-seq data from immune cells of patients with ICC. This PRG signature not only reflects the TME immune status but also provides new biomarkers for ICC prognosis.

Most sequencing-based spatial transcriptomics (ST) technologies do not achieve single-cell resolution where each captured location (spot) may contain a mixture of cells from heterogeneous cell types, and several cell-type decomposition methods have been proposed to estimate cell type proportions of each spot by integrating with single-cell RNA sequencing (scRNA-seq) data. However, these existing methods did not fully consider the effect of distribution difference between scRNA-seq and ST data for decomposition, leading to biased cell-type-specific genes derived from scRNA-seq for ST data. To address this issue, we develop an instance-based transfer learning framework to adjust scRNA-seq data by ST data to correctly match cell-type-specific gene expression. We evaluate the effect of raw and adjusted scRNA-seq data on cell-type decomposition by eight leading decomposition methods using both simulated and real datasets. Experimental results show that data adjustment can effectively reduce distribution difference and improve decomposition, thus enabling for a more precise depiction on spatial organization of cell types. We highlight the importance of data adjustment in integrative analysis of scRNA-seq with ST data and provide guidance for improved cell-type decomposition.

The incidence and mortality rates of lung cancer are exceptionally high. Many patients are diagnosed with early stage lung cancer but experience rapid recurrence post-surgery. Many research studies have shown that the unfavorable prognosis of patients may be associated with micro-metastasis in the lymph nodes. Our research aimed to develop a nomogram to predict the prognosis of lung adenocarcinoma (LUAD). Single-cell RNA sequencing (scRNA-seq) data were analyzed to identify 11 cell clusters. Patterns of incoming and outgoing signals were identified across the entire cell population. A weighted gene co-expression network analysis (WGCNA) was conducted to uncover critical genes in LUAD. The intersecting marker genes were used to construct the prognostic model. scRNA-seq data were analyzed to identify 19 cell clusters. We identified 3,464 marker genes from the scRNA-seq dataset, 1,994 differentially expressed genes from the bulk RNA sequencing (RNA-seq) dataset, and 1,863 genes associated with a key module identified by the WGCNA. After performing the intersection, univariate Cox, and least absolute shrinkage and selection operator analyses, a prognostic model was established based on the expression levels of 13 signature genes. Subsequent functional experiments confirmed the role of selected regulated genes. Through the integration of scRNA-seq data and bulk RNA-seq data, we developed an innovative model to predict the prognosis of patients. The risk score was found to be a significant independent predictor and clinical-pathological features of LUAD.

Joint learning of multiple gene networks from single-cell gene expression data

Unsupervised clustering is a fundamental step of single-cell RNA-sequencing (scRNA-seq) data analysis. This issue has inspired several clustering methods to classify cells in scRNA-seq data. However, accurate prediction of the cell clusters remains a substantial challenge. In this study, we propose a new algorithm for scRNA-seq data clustering based on Sparse Optimization and low-rank matrix factorization (scSO). We applied our scSO algorithm to analyze multiple benchmark datasets and showed that the cluster number predicted by scSO was close to the number of reference cell types and that most cells were correctly classified. Our scSO algorithm is available at https://github.com/QuKunLab/scSO. Overall, this study demonstrates a potent cell clustering approach that can help researchers distinguish cell types in single- scRNA-seq data.

In this paper, we introduce Gene Knockout Inference (GenKI), a virtual knockout (KO) tool for gene function prediction using single-cell RNA sequencing (scRNA-seq) data in the absence of KO samples when only wild-type (WT) samples are available. Without using any information from real KO samples, GenKI is designed to capture shifting patterns in gene regulation caused by the KO perturbation in an unsupervised manner and provide a robust and scalable framework for gene function studies. To achieve this goal, GenKI adapts a variational graph autoencoder (VGAE) model to learn latent representations of genes and interactions between genes from the input WT scRNA-seq data and a derived single-cell gene regulatory network (scGRN). The virtual KO data is then generated by computationally removing all edges of the KO gene—the gene to be knocked out for functional study—from the scGRN. The differences between WT and virtual KO data are discerned by using their corresponding latent parameters derived from the trained VGAE model. Our simulations show that GenKI accurately approximates the perturbation profiles upon gene KO and outperforms the state-of-the-art under a series of evaluation conditions. Using publicly available scRNA-seq data sets, we demonstrate that GenKI recapitulates discoveries of real-animal KO experiments and accurately predicts cell type-specific functions of KO genes. Thus, GenKI provides an in-silico alternative to KO experiments that may partially replace the need for genetically modified animals or other genetically perturbed systems.

Modeling cellular dynamics from single-cell RNA sequencing (scRNA-seq) data is critical for understanding cell development and underlying gene regulatory relationships. Many current methods rely on single-cell velocity to obtain pseudotime, which can lead to inconsistencies between pseudotime and velocity. It is challenging to simultaneously infer cell pseudotime and gene interaction networks, especially in multi-branch differentiation scenarios. We present single-cell Piecewise Network (scPN), a novel high-dimensional dynamical modeling approach that iteratively extracts temporal patterns and inter-gene relationships from scRNA-seq data. To tackle multi-branch differentiation challenges, scPN models gene regulatory dynamics using piecewise gene–gene interaction networks, offering an interpretable framework for deciphering complex gene regulation patterns over time. Results on synthetic data and multiple scRNA-seq datasets demonstrate the superior performance of scPN in reconstructing cellular dynamics and identifying key transcription factors involved in development compared to existing methods. To the best of our knowledge, scPN is the first attempt at modeling that can recover pseudotime, velocity fields, and gene interactions all at once on multi-branch datasets.

Tumor-associated macrophages (TAMs) are crucial in hepatocellular carcinoma (HCC) development and invasion. This study explores monocyte/ macrophage-associated gene expression profiles in HCC, constructs a prognostic model based on these genes, and examines its relationship with drug resistance and immune therapy responses. Single-cell RNA sequencing(scRNA-seq) data from 10 HCC tissue biopsy samples, totaling 24,597 cells, were obtained from the GEO database to identify monocyte/macrophage-associated genes. A prognostic model was constructed and validated using external datasets and Western blot. Relationships between the model, clinical correlates, drug sensitivity, and immune therapy responses were investigated. From scRNA-seq data, 2,799 monocyte/macrophage marker genes were identified. Using the TCGA dataset, a prognostic model based on the single-gene UQCRH was constructed, stratifying patients into high-risk and low-risk groups based on overall survival rates. High-risk group patients showed reduced survival rates and higher UQCRH expression in tumor tissues. Western blot analysis further confirmed the elevated expression of UQCRH in HCC cell lines. Spatial transcriptomics analysis revealed that high UQCRH expression co-localized with malignant cells in the tumor tissue. Drug sensitivity analysis revealed that the high-risk group had lower sensitivity to sorafenib and axitinib. Immune therapy response analysis indicated poorer outcomes in the high-risk group, with more pronounced APC inhibition and a weaker IFN-II response. Clinical indicator analysis showed a positive correlation between high UQCRH expression and tumor invasion. Enrichment analysis of UQCRH and associated molecules indicated involvement in oxidative phosphorylation and mitochondrial electron transport. This study introduces a prognostic model for HCC patients based on monocyte/macrophage marker genes. The single-gene model predicts HCC patient survival and treatment outcomes, identifying high-risk individuals with varying drug sensitivities and immune suppression states.

Ischemic stroke (IS) is a leading cause of adult disability worldwide. The inflammatory processes involved are complex, making it challenging to fully understand the pathological mechanisms of IS. Phagocytosis plays an important role in eliminating neurotoxic or damaged neurons resulting from inflammatory responses. This study employed bioinformatics methods to analyze single-cell RNA sequencing (scRNA-seq) data to investigate the cell types and molecular biological processes involved in IS. scRNA-seq data for IS were obtained from the Gene Expression Omnibus (GEO). Following sample screening and reprocessing, 5,582 single cells were identified from healthy controls and patients with IS. Uniform manifold approximation and projection (UMAP) was utilized to further explore the cellular composition in IS. Functional enrichment analysis of differentially expressed genes was conducted to identify transcriptional regulators, whereas cell developmental trajectories were predicted to uncover potential cell fate decisions. iTALK was employed to identify potential ligand-receptor axes within the cell-type immune microenvironment of IS. Based on scRNA-seq data analysis, we identified four cell types and their associated subclusters, along with genes exhibiting significant differential expression within these subclusters. Phagocytosis was significantly enriched in cell types linked to IS, while the differentiation trajectories of subpopulations in IS was different. Additionally, multiple receptor-ligand axes were identified, indicating diverse interactions within the immune microenvironment of IS. This study demonstrated that phagocytosis in IS cell types critically influences disease progression. It also predicted the trajectories of infarct cells. These findings provide valuable insights into the molecular and cellular mechanisms underlying IS and highlight potential pathways for therapeutic intervention.