Genetic mapping of cell type specificity for complex traits

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.

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.

BackgroundThe rapid accumulation of single-cell RNA sequencing (scRNA-seq) data presents unique opportunities to decode the genetically mediated cell-type specificity in complex diseases. Here, we develop a new method, scGWAS, which effectively leverages scRNA-seq data to achieve two goals: (1) to infer the cell types in which the disease-associated genes manifest and (2) to construct cellular modules which imply disease-specific activation of different processes.ResultsscGWAS only utilizes the average gene expression for each cell type followed by virtual search processes to construct the null distributions of module scores, making it scalable to large scRNA-seq datasets. We demonstrated scGWAS in 40 genome-wide association studies (GWAS) datasets (average sample size N ≈ 154,000) using 18 scRNA-seq datasets from nine major human/mouse tissues (totaling 1.08 million cells) and identified 2533 trait and cell-type associations, each with significant modules for further investigation. The module genes were validated using disease or clinically annotated references from ClinVar, OMIM, and pLI variants.ConclusionsWe showed that the trait-cell type associations identified by scGWAS, while generally constrained to trait-tissue associations, could recapitulate many well-studied relationships and also reveal novel relationships, providing insights into the unsolved trait-tissue associations. Moreover, in each specific cell type, the associations with different traits were often mediated by different sets of risk genes, implying disease-specific activation of driving processes. In summary, scGWAS is a powerful tool for exploring the genetic basis of complex diseases at the cell type level using single-cell expression 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.

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.

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

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.

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.

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.

Based on benchmarking various methods and analyzing multiple real scRNA-seq datasets, a computational platform/workflow and a set of tips for best practices are developed for analyzing causal interactions or relationships in single cells.

Hypertrophic cardiomyopathy (HCM) is a common inherited heart disease characterized by left ventricular hypertrophy and excessive myocardial contraction. Endothelial cells (ECs) play a key role in the pathogenesis of HCM, and their secreted growth factors regulate the growth and function of cardiomyocytes. However, the specific molecular mechanisms and biomarkers associated with endothelial dysfunction in HCM remain unclear. Single-cell RNA sequencing (scRNA-seq) data were downloaded from the Gene Expression Omnibus (GEO) database (GSE255296), and bulk RNA sequencing data were obtained from GSE249925 and GSE180313. The scRNA-seq data were processed using the Seurat package for cell type annotation and batch effect removal. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed, and pathway activity was evaluated using the AUCell algorithm. SCENIC was applied to construct transcription factor-specific regulatory networks. Differential gene expression analysis was conducted using the limma package, and key genes were determined by intersecting differentially expressed genes (DEGs) with EC-specific highly expressed genes. The diagnostic efficacy of these key genes was assessed via receiver operating characteristic (ROC) curve analysis. The mRNA expression of genes in HUVECs treated with or without angiotensin II (Ang II) was assessed using qPCR. Cell viability was determined using the CCK-8 assay, and cell invasive ability was evaluated through Transwell assays to investigate the regulatory effect of TMTC1 overexpression on Ang II-induced injury in HUVECs. Seven major cell types were identified, with cardiomyocytes and ECs being the main cell populations in HCM. KEGG enrichment analysis revealed significant activation of pathways such as the transforming growth factor-beta (TGF-β) signaling pathway in ECs. High-activity transcription factors NFIB, MECOM, and FLI1 were identified in ECs, whose target genes were enriched in processes related to cell morphogenesis. Three key genes-TMTC1, RPS6KA2, and F8-were selected based on differential expression and EC specificity. These genes were lowly expressed in HCM and exhibited high diagnostic efficacy, with area under the curve (AUC) values > 0.7. Ultimately, the CCK-8 assay results indicated thatAng II inhibited the viability of HUVECs, and this inhibitory effect could be reversed by the overexpression of TMTC1 (oe-TMTC1). Furthermore, Transwell assays indicated that Ang II impaired cellular migration, and this impairment could be restored by oe-TMTC1. The identified key genes and pathways provide new targets for therapeutic intervention and contribute to a deeper understanding of the pathogenesis of HCM.

Chronic pain is a complex clinical problem comprising multiple conditions that may share a common genetic profile. Genome-wide association studies (GWAS) have identified many risk loci whose cell-type context remains unclear. Here, we integrated GWAS data on chronic pain (N = 1,235,695) with single-cell RNA sequencing (scRNA-seq) data from human brain and dorsal root ganglia (hDRG), and single-cell chromatin accessibility data from human brain and mouse dorsal horn. Pain-associated variants were enriched in glutamatergic neurons; mainly in prefrontal cortex, hippocampal CA1–3, and amygdala. In hDRG, the hPEP.TRPV1/A1.2 neuronal subtype showed robust enrichment. Chromatin accessibility analyses revealed variant enrichment in excitatory and inhibitory neocortical neurons in brain and in midventral neurons and oligodendrocyte precursor cells in the mouse dorsal horn. Gene-level heritability in the brain highlighted roles for kinase activity, GABAergic synapses, axon guidance, and neuron projection development. In hDRG, implicated genes related to glutamatergic signaling and neuronal projection. In cervical DRG of patients with acute or chronic pain (N = 12), scRNA-seq data from neuronal or non-neuronal cells were enriched for chronic pain-associated genes (e.g., EFNB2, GABBR1, NCAM1, SCN11A). This cell-type-specific genetic architecture of chronic pain across central and peripheral nervous system circuits provides a foundation for targeted translational research.

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.

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.

Childhood cancer survivors (CCSs) face lifelong side effects related to their treatment with chemotherapy. Anthracycline agents, such as doxorubicin (DOX), are important in the treatment of childhood cancers but are associated with cardiotoxicity. Cardiac toxicities represent a significant source of chronic disability that cancer survivors face; despite this, the chronic cardiotoxicity phenotype and how it relates to acute toxicity remains poorly defined. To address this critical knowledge gap, we studied the acute effect of DOX on murine cardiac nonmyocytes in vivo. Determination of the acute cellular effects of DOX on nonmyocytes, a cell pool with finite replicative capacity, provides a basis for understanding the pathogenesis of the chronic heart disease that CCSs face. To investigate the acute cellular effects of DOX, we present single-cell RNA sequencing (scRNAseq) data from homeostatic cardiac nonmyocytes and compare it with preexisting datasets, as well as a novel CyTOF datasets. SCANPY, a python-based single-cell analysis, was used to assess the heterogeneity of cells detected in scRNAseq and CyTOF. To further assist in CyTOF data annotation, joint analyses of scRNAseq and CyTOF data using an artificial neural network known as sparse autoencoder for clustering, imputation, and embedding (SAUCIE) are performed. Lastly, the panel is tested on a mouse model of acute DOX exposure at two time points (24 and 72 h) after the last dose of doxorubicin and examined with joint clustering. In sum, we report the first ever CyTOF study of cardiac nonmyocytes and characterize the effect of acute DOX exposure with scRNAseq and CyTOF.NEW & NOTEWORTHY We describe the first mass cytometry studies of murine cardiac nonmyocytes. The mass cytometry panel is compared with single-cell RNA sequencing data. Homeostatic cardiac nonmyocytes are characterized by mass cytometry to identify and quantify four major cell populations: endothelial cells, fibroblasts, leukocytes, and pericytes. The single-cell acute nonmyocyte response to doxorubicin is studied at 24 and 72 h after doxorubicin exposure given daily for 5 days at a dose of 4 mg/kg/day.