scRNA-seq Datasets Research Articles

Tumor necrosis factor receptor 2 promotes endothelial cell-mediated suppression of CD8+ T cells through tuning glycolysis in chemoresistance of breast cancer

BackgroundT cells play a pivotal role in chemotherapy-triggered anti-tumor effects. Emerging evidence underscores the link between impaired anti-tumor immune responses and resistance to paclitaxel therapy in triple-negative breast cancer (TNBC). Tumor-related endothelial cells (ECs) have potential immunoregulatory activity. However, how ECs regulate T cell activity during TNBC chemotherapy remains poorly understood.MethodsSingle-cell analysis of ECs in patients with TNBC receiving paclitaxel therapy was performed using an accessible single-cell RNA sequencing (scRNA-seq) dataset to identify key EC subtypes and their immune characteristics. An integrated analysis of a tumor-bearing mouse model, immunofluorescence, and a spatial transcriptome dataset revealed the spatial relationship between ECs, especially Tumor necrosis factor receptor (TNFR) 2+ ECs, and CD8+ T cells. RNA sequencing, CD8+ T cell proliferation assays, flow cytometry, and bioinformatic analyses were performed to explore the immunosuppressive function of TNFR2 in ECs. The downstream metabolic mechanism of TNFR2 was further investigated using RNA sequencing, cellular glycolysis assays, and western blotting.ResultsIn this study, we identified an immunoregulatory EC subtype, characterized by enhanced TNFR2 expression in non-responders. By a mouse model of TNBC, we revealed a dynamic reduction in the proportion of the CD8+ T cell-contacting tumor vessels that could co-localize spatially with CD8+ T cells during chemotherapy and an increased expression of TNFR2 by ECs. TNFR2 suppresses glycolytic activity in ECs by activating NF-κB signaling in vitro. Tuning endothelial glycolysis enhances programmed death-ligand (PD-L) 1-dependent inhibitory capacity, thereby inducing CD8+ T cell suppression. In addition, TNFR2+ ECs showed a greater spatial affinity for exhausted CD8+ T cells than for non-exhausted CD8+ T cells. TNFR2 blockade restores impaired anti-tumor immunity in vivo, leading to the loss of PD-L1 expression by ECs and enhancement of CD8+ T cell infiltration into the tumors.ConclusionsThese findings reveal the suppression of CD8+ T cells by ECs in chemoresistance and indicate the critical role of TNFR2 in driving the immunosuppressive capacity of ECs via tuning glycolysis. Targeting endothelial TNFR2 may serve as a potent strategy for treating TNBC with paclitaxel.

Computing cell state discriminates the aberrant hematopoiesis and activated microenvironment in Myelodysplastic syndrome (MDS) through a single cell genomic study

BackgroundMyelodysplastic syndrome (MDS) is a complicated hematopoietic malignancy characterized by bone marrow (BM) dysplasia with symptoms like anemia, neutropenia, or thrombocytopenia. MDS exhibits considerable heterogeneity in prognosis, with approximately 30% of patients progressing to acute myeloid leukemia (AML). Single cell RNA-sequencing (scRNA-seq) is a new and powerful technique to profile disease landscapes. However, the current available scRNA-seq datasets for MDS are only focused on CD34+ hematopoietic progenitor cells. We argue that using entire BM cell for MDS studies probably will be more informative for understanding the pathophysiology of MDS.MethodsFive MDS patients and four healthy donors were enrolled in the study. Unsorted cells from BM aspiration were collected for scRNA-seq analysis to profile overall alteration in hematopoiesis.ResultsStandard scRNA-seq analysis of unsorted BM cells successfully profiles deficient hematopoiesis in all five MDS patients, with three classified as high-risk and two as low-risk. While no significant increase in mutation burden was observed, high-risk MDS patients exhibited T-cell activation and abnormal myelogenesis at the stages between hematopoietic stem and progenitor cells (HSPC) and granulocyte–macrophage progenitors (GMP). Transcriptional factor analysis on the aberrant myelogenesis suggests that the epigenetic regulator chromatin structural protein-encoding gene HMGA1 is highly activated in the high-risk MDS group and moderately activated in the low-risk MDS group. Perturbation of HMGA1 by CellOracle simulated deficient hematopoiesis in mouse Lineage-negative (Lin-) BM cells. Projecting MDS and AML cells on a BM cell reference by our newly developed MarcoPolo pipeline intuitively visualizes a connection for myeloid leukemia development and abnormalities of hematopoietic hierarchy, indicating that it is technically feasible to integrate all diseased bone marrow cells on a common reference map even when the size of the cohort reaches to 1,000 patients or more.ConclusionThrough scRNA-seq analysis on unsorted cells from BM aspiration samples of MDS patients, this study systematically profiled the development abnormalities in hematopoiesis, heterogeneity of risk, and T-cell microenvironment at the single cell level.

Systematic Analysis of Human Colorectal Cancer scRNA-seq Revealed Limited Pro-tumoral IL-17 Production Potential in Gamma Delta T Cells.

Gamma delta (γδ) T cells play a crucial role in anti-tumor immunity due to their cytotoxic properties. However, the role and extent of γδ T cells in production of pro-tumorigenic interleukin- 17 (IL-17) within the tumor microenvironment (TME) of colorectal cancer (CRC) remains controversial. In this study, we re-analyzed nine published human CRC whole-tissue single-cell RNA sequencing (scRNA-seq) datasets, identifying 18,483 γδ T cells out of 951,785 total cells, in the neoplastic or adjacent normal tissue of 165 human CRC patients. Our results confirm that tumor-infiltrating γδ T cells exhibit high cytotoxicity-related transcription in both tumor and adjacent normal tissues, but critically, none of the γδ T cell clusters showed IL-17 production potential. We also identified various γδ T cell subsets, including Teff, TRM, Tpex, and Tex, and noted an increased expression of cytotoxic molecules in tumor-infiltrating γδ T cells compared to their normal area counterparts. Our work demonstrates that γδ T cells in CRC primarily function as cytotoxic effector cells rather than IL-17 producers, mitigating the concerns about their potential pro-tumorigenic roles in CRC, highlighting the importance of accurately characterizing these cells for cancer immunotherapy research and the unneglectable cross-species discrepancy between the mouse and human immune system in the study of cancer immunology.

Accurate identification of single-cell types via correntropy-based Sparse PCA combining hypergraph and fusion similarity

ABSTRACT The advent of single-cell RNA sequencing (scRNA-seq) technology enables researchers to gain deep insights into cellular heterogeneity. However, the high dimensionality and noise of scRNA-seq data pose significant challenges to clustering. Therefore, we propose a new single-cell type identification method, called CHLSPCA, to address these challenges. In this model, we innovatively combine correntropy with PCA to address the noise and outliers inherent in scRNA-seq data. Meanwhile, we integrate the hypergraph into the model to extract more valuable information from the local structure of the original data. Subsequently, to capture crucial similarity information not considered by the PCA model, we employ the Gaussian kernel function and the Euclidean metric to mine the similarity information between cells, and incorporate this information into the model as the similarity constraint. Furthermore, the principal components (PCs) of PCA are very dense. A new sparse constraint is introduced into the model to gain sparse PCs. Finally, based on the principal direction matrix learned from CHLSPCA, we conduct extensive downstream analyses on real scRNA-seq datasets. The experimental results show that CHLSPCA performs better than many popular clustering methods and is expected to promote the understanding of cellular heterogeneity in scRNA-seq data analysis and support biomedical research.

Expression patterns of keratin family members during tooth development and the role of keratin 17 in cytodifferentiation of stratum intermedium and stellate reticulum.

Keratins are typical intermediate filament proteins of the epithelium that exhibit highly specific expression patterns related to the epithelial type and stage of cellular differentiation. They are important for cytoplasmic stability and epithelial integrity and are involved in various intracellular signaling pathways. Several keratins are associated with enamel formation. However, information on their expression patterns during tooth development remains lacking. In this study, we analyzed the spatiotemporal expression of keratin family members during tooth development using single-cell RNA-sequencing (scRNA-seq) and microarray analysis. scRNA-seq datasets from postnatal Day 1 mouse molars revealed that several keratins are highly expressed in the dental epithelium, indicating the involvement of keratin family members in cellular functions. Among various keratins, keratin 5 (Krt5), keratin 14 (Krt14), and keratin 17 (Krt17) are highly expressed in the tooth germ; KRT17 is specifically expressed in the stratum intermedium (SI) and stellate reticulum (SR). Depletion of Krt17 did not affect cell proliferation in the dental epithelial cell line SF2 but suppressed their differentiation ability. These results suggest that Krt17 is essential for SI cell differentiation. Furthermore, scRNA-seq results indicated that Krt5, Krt14, and Krt17 exhibited distinct expression patterns in ameloblast, SI, and SR cells. Our findings contribute to the elucidation of novel mechanisms underlying tooth development.

Graph attention network with convolutional layer for predicting gene regulations from single-cell ribonucleic acid sequence data

Reconstructing gene regulatory networks (GRNs) is an important task to reveal the regulatory relationship between genes and understand the mechanism of intracellular gene expression regulation. With the development of single-cell ribonucleic acid sequencing (scRNA-seq) technology, researchers begin to attempt to infer GRN within cells. In this paper, we propose graph attention network with convolutional layer (GATCL) to infer the latent interactions between transcription factors (TFs) and target genes in GRN. Firstly, GATCL uses graph attention network (GAT), which can effectively extract information about genes and TFs. Secondly, we combine multi-head attention layer with one-head attention layer and propose a new method of using convolution instead of weight matrix. Thirdly, we use the exponential linear unit (ELU) activation function to replace the leaky rectified linear unit (LReLU) commonly used in GAT, which further improves the accuracy of GATCL. The AUROC of our method on seven scRNA-seq datasets with four types of ground-truth networks reached an average of 0.827, which is higher than other state-of-the-art models. GATCL applies a supervised deep learning algorithm to solve the problems existing in the inference of GRN from scRNA-seq in the engineering field, and the effectiveness of this method is verified by a large number of experiments.

Diffusion on PCA-UMAP Manifold: The Impact of Data Structure Preservation to Denoise High-Dimensional Single-Cell RNA Sequencing Data.

Single-cell transcriptomics (scRNA-seq) is revolutionizing biological research, yet it faces challenges such as inefficient transcript capture and noise. To address these challenges, methods like neighbor averaging or graph diffusion are used. These methods often rely on k-nearest neighbor graphs from low-dimensional manifolds. However, scRNA-seq data suffer from the 'curse of dimensionality', leading to the over-smoothing of data when using imputation methods. To overcome this, sc-PHENIX employs a PCA-UMAP diffusion method, which enhances the preservation of data structures and allows for a refined use of PCA dimensions and diffusion parameters (e.g., k-nearest neighbors, exponentiation of the Markov matrix) to minimize noise introduction. This approach enables a more accurate construction of the exponentiated Markov matrix (cell neighborhood graph), surpassing methods like MAGIC. sc-PHENIX significantly mitigates over-smoothing, as validated through various scRNA-seq datasets, demonstrating improved cell phenotype representation. Applied to a multicellular tumor spheroid dataset, sc-PHENIX identified known extreme phenotype states, showcasing its effectiveness. sc-PHENIX is open-source and available for use and modification.

T-distributed Stochastic Neighbor Network for unsupervised representation learning

Unsupervised representation learning (URL) is still lack of a reasonable operator (e.g. convolution kernel) for exploring meaningful structural information from generic data including vector, image and tabular data. In this paper, we propose a simple end-to-end T-distributed Stochastic Neighbor Network (TsNet) for URL with clustering downstream task. Concretely, our TsNet model has three major components: (1) an adaptive connectivity distribution learning module is presented to construct a pairwise graph for preserving the local structure of generic data; (2) a T-distributed stochastic neighbor embedding based loss function is designed to learn a transformation between embeddings and original data, which improves the discrimination of representations; (3) a nonlinear parametric mapping is learned via our TsNet on an unsupervised generalized manner, which can address the “out-of-sample” issue. By combining these components, our method is able to considerably outperform previous related unsupervised learning approaches on visualization and clustering of generic data. A simple deep neural network equipped on our model respectively achieves 74.90%, 76.56% ACC and NMI, which is 8% relative improvement over previous state-of-the-art on real single-cell RNA-sequencing (scRNA-seq) datasets clustering.

A novel CAF-cancer cell crosstalk-related gene prognostic index based on machine learning: prognostic significance and prediction of therapeutic response in head and neck squamous cell carcinoma

BackgroundCancer-associated fibroblast (CAF)-cancer cell crosstalk (CCCT) plays an important role in tumor microenvironment shaping and immunotherapy response. Current prognostic indexes are insufficient to accurately assess immunotherapy response in patients with head and neck squamous cell carcinoma (HNSCC). This study aimed to develop a CCCT-related gene prognostic index (CCRGPI) for assessing the prognosis and response to immune checkpoint inhibitor (ICI) therapy of HNSCC patients.MethodsTwo cellular models, the fibroblast-cancer cell indirect coculture (FCICC) model, and the fibroblast-cancer cell organoid (FC-organoid) model, were constructed to visualize the crosstalk between fibroblasts and cancer cells. Based on a HNSCC scRNA-seq dataset, the R package CellChat was used to perform cell communication analysis to identify gene pairs involved in CCCT. Least absolute shrinkage and selection operator (LASSO) regression was then applied to further refine the selection of these gene pairs. The selected gene pairs were subsequently subjected to stepwise regression to develop CCRGPI. We further performed a comprehensive analysis to determine the molecular and immune characteristics, and prognosis associated with ICI therapy in different CCRGPI subgroups. Finally, the connectivity map (CMap) analysis and molecular docking were used to screen potential therapeutic drugs.ResultsFCICC and FC-organoid models showed that cancer cells promoted the activation of fibroblasts into CAFs, that CAFs enhanced the invasion of cancer cells, and that CCCT was somewhat heterogeneous. The CCRGPI was developed based on 4 gene pairs: IGF1-IGF1R, LGALS9-CD44, SEMA5A-PLXNA1, and TNXB-SDC1. Furthermore, a high CCRGPI score was identified as an adverse prognostic factor for overall survival (OS). Additionally, a high CCRGPI was positively correlated with the activation of the P53 pathway, a high TP53 mutation rate, and decreased benefit from ICI therapy but was inversely associated with the abundance of various immune cells, such as CD4+ T cells, CD8+ T cells, and B cells. Moreover, Ganetespib was identified as a potential drug for HNSCC combination therapy.ConclusionsThe CCRGPI is reliable for predicting the prognosis and immunotherapy response of HSNCC patients and may be useful for guiding the individualized treatment of HNSCC patients.

ScBoolSeq: Linking scRNA-seq statistics and Boolean dynamics.

Boolean networks are largely employed to model the qualitative dynamics of cell fate processes by describing the change of binary activation states of genes and transcription factors with time. Being able to bridge such qualitative states with quantitative measurements of gene expression in cells, as scRNA-seq, is a cornerstone for data-driven model construction and validation. On one hand, scRNA-seq binarisation is a key step for inferring and validating Boolean models. On the other hand, the generation of synthetic scRNA-seq data from baseline Boolean models provides an important asset to benchmark inference methods. However, linking characteristics of scRNA-seq datasets, including dropout events, with Boolean states is a challenging task. We present scBoolSeq, a method for the bidirectional linking of scRNA-seq data and Boolean activation state of genes. Given a reference scRNA-seq dataset, scBoolSeq computes statistical criteria to classify the empirical gene pseudocount distributions as either unimodal, bimodal, or zero-inflated, and fit a probabilistic model of dropouts, with gene-dependent parameters. From these learnt distributions, scBoolSeq can perform both binarisation of scRNA-seq datasets, and generate synthetic scRNA-seq datasets from Boolean traces, as issued from Boolean networks, using biased sampling and dropout simulation. We present a case study demonstrating the application of scBoolSeq's binarisation scheme in data-driven model inference. Furthermore, we compare synthetic scRNA-seq data generated by scBoolSeq with BoolODE's, data for the same Boolean Network model. The comparison shows that our method better reproduces the statistics of real scRNA-seq datasets, such as the mean-variance and mean-dropout relationships while exhibiting clearly defined trajectories in two-dimensional projections of the data.

CrossMP: Enabling Cross-Modality Translation between Single-Cell RNA-Seq and Single-Cell ATAC-Seq through Web-Based Portal.

In recent years, there has been a growing interest in profiling multiomic modalities within individual cells simultaneously. One such example is integrating combined single-cell RNA sequencing (scRNA-seq) data and single-cell transposase-accessible chromatin sequencing (scATAC-seq) data. Integrated analysis of diverse modalities has helped researchers make more accurate predictions and gain a more comprehensive understanding than with single-modality analysis. However, generating such multimodal data is technically challenging and expensive, leading to limited availability of single-cell co-assay data. Here, we propose a model for cross-modal prediction between the transcriptome and chromatin profiles in single cells. Our model is based on a deep neural network architecture that learns the latent representations from the source modality and then predicts the target modality. It demonstrates reliable performance in accurately translating between these modalities across multiple paired human scATAC-seq and scRNA-seq datasets. Additionally, we developed CrossMP, a web-based portal allowing researchers to upload their single-cell modality data through an interactive web interface and predict the other type of modality data, using high-performance computing resources plugged at the backend.

Establishment of an ovarian cancer exhausted CD8+T cells-related genes model by integrated analysis of scRNA-seq and bulk RNA-seq

Ovarian cancer (OC) was the fifth leading cause of cancer death and the deadliest gynecological cancer in women. This was largely attributed to its late diagnosis, high therapeutic resistance, and a dearth of effective treatments. Clinical and preclinical studies have revealed that tumor-infiltrating CD8+T cells often lost their effector function, the dysfunctional state of CD8+T cells was known as exhaustion. Our objective was to identify genes associated with exhausted CD8+T cells (CD8TEXGs) and their prognostic significance in OC. We downloaded the RNA-seq and clinical data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CD8TEXGs were initially identified from single-cell RNA-seq (scRNA-seq) datasets, then univariate Cox regression, the least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression were utilized to calculate risk score and to develop the CD8TEXGs risk signature. Kaplan–Meier analysis, univariate Cox regression, multivariate Cox regression, time-dependent receiver operating characteristics (ROC), nomogram, and calibration were conducted to verify and evaluate the risk signature. Gene set enrichment analyses (GSEA) in the risk groups were used to figure out the closely correlated pathways with the risk group. The role of risk score has been further explored in the homologous recombination repair deficiency (HRD), BRAC1/2 gene mutations and tumor mutation burden (TMB). A risk signature with 4 CD8TEXGs in OC was finally built in the TCGA database and further validated in large GEO cohorts. The signature also demonstrated broad applicability across various types of cancer in the pan-cancer analysis. The high-risk score was significantly associated with a worse prognosis and the risk score was proven to be an independent prognostic biomarker. The 1-, 3-, and 5-years ROC values, nomogram, calibration, and comparison with the previously published models confirmed the excellent prediction power of this model. The low-risk group patients tended to exhibit a higher HRD score, BRCA1/2 gene mutation ratio and TMB. The low-risk group patients were more sensitive to Poly-ADP-ribose polymerase inhibitors (PARPi). Our findings of the prognostic value of CD8TEXGs in prognosis and drug response provided valuable insights into the molecular mechanisms and clinical management of OC.

Disease-specific T cell receptors maintain pathogenic T helper cell responses in postinfectious Lyme arthritis.

Antibiotic-Refractory Lyme Arthritis (ARLA) involves a complex interplay of T cell responses targeting Borrelia burgdorferi antigens succeeding towards autoantigens by epitope spreading. However, the precise molecular mechanisms driving the pathogenic T cell response in ARLA remain unclear. Our aim was to elucidate the molecular program of disease-specific Th cells. Using flow cytometry, high-throughput T cell receptor (TCR) sequencing and scRNA-seq of CD4+ Th cells isolated from the joints of European ARLA patients, we aimed at inferring antigen specificity through unbiased analysis of TCR repertoire patterns, identifying surrogate markers for disease-specific TCRs and connecting TCR specificity to transcriptional patterns. PD-1hiHLA-DR+CD4+ effector T cells were clonally expanded within the inflamed joints and persisted throughout disease course. Among these cells, we identified a distinct TCRβ motif restricted to HLA-DRB1*11 or *13 alleles. These alleles, being underrepresented in North American ARLA patients, were unexpectedly prevalent in our European cohort. The identified TCRβ motif served as surrogate marker for a convergent TCR response specific to ARLA, distinguishing it from other rheumatic diseases. In the scRNA-seq dataset, the TCRβ motif particularly mapped to peripheral T helper (TPH) cells displaying signs of sustained proliferation, continuous TCR signaling, and expressing CXCL13 and IFN-γ. By inferring disease-specific TCRs from synovial T cells we identified a convergent TCR response in the joints of ARLA patients that continuously fueled the expansion of TPH cells expressing a pathogenic cytokine effector program. The identified TCRs will aid in uncovering the major antigen targets of the maladaptive immune response. Supported by the German Research Foundation (DFG) MO 2160/4-1; the Federal Ministry of Education and Research (BMBF; Advanced Clinician Scientist-Program INTERACT; 01EO2108) embedded in the Interdisciplinary Center for Clinical Research (IZKF) of the University Hospital Würzburg; the German Center for Infection Research (DZIF; Clinical Leave Program; TI07.001_007) and the Interdisciplinary Center for Clinical Research (IZKF) Würzburg (Clinician Scientist Program, Z-2/CSP-30).

EzSingleCell: an integrated one-stop single-cell and spatial omics analysis platform for bench scientists

ezSingleCell is an interactive and easy-to-use application for analysing various single-cell and spatial omics data types without requiring prior programing knowledge. It combines the best-performing publicly available methods for in-depth data analysis, integration, and interactive data visualization. ezSingleCell consists of five modules, each designed to be a comprehensive workflow for one data type or task. In addition, ezSingleCell allows crosstalk between different modules within a unified interface. Acceptable input data can be in a variety of formats while the output consists of publication ready figures and tables. In-depth manuals and video tutorials are available to guide users on the analysis workflows and parameter adjustments to suit their study aims. ezSingleCell’s streamlined interface can analyse a standard scRNA-seq dataset of 3000 cells in less than five minutes. ezSingleCell is available in two forms: an installation-free web application (https://immunesinglecell.org/ezsc/) or a software package with a shinyApp interface (https://github.com/JinmiaoChenLab/ezSingleCell2) for offline analysis.

S CIRCLE-An interactive visual exploration tool for single cell RNA-Seq data.

sCIRCLE (single-Cell Interactive Real-time Computer visualization for Low-dimensional Exploration) is a tool for exploratory analysis of single cell RNA-seq (scRNA-seq) data sets, with a focus on bacterial scRNA-seq. The software takes an information design perspective to re-envision visually and interactively exploring low dimensional representations of scRNA-Seq data. Users can project cells in various 3D and 2D spaces and interactively query and paint cells using rich metadata sets reporting on cell cluster, gene function, and gene expression. As a standalone application it contains, among other features, options for dimensionality reduction, navigation and interaction with data in 3d and 2d space, gene filtering, fold change and metacell computation as well as various capabilities for visualization, data import and export.

Constructing an immune-related prognostic signature for predicting prognosis and immune response in hepatocellular carcinoma

BackgroundCurrently, there are few studies on immune-related prognostic analysis of hepatocellular carcinoma (HCC). Our aim was to establish an immune-correlated prognostic model for HCC. MethodsImmune-associated cells were obtained from the scRNA-seq dataset (GSE149614) of HCC. Differentially expressed genes between normal and tumor cells from immune-associated cells and the immune-related genes from the ImmPort database were used to identify immune-related differentially expressed genes (IRDEGs). Subsequently, the risk model was established in the TCGA-LIHC cohort (n = 438) from the Cancer Genome Atlas (TCGA) database by using Kaplan-Meier (K-M) survival curve, univariate/multivariate Cox regression analysis. Subsequently, we further analyzed tumor immune microenvironment characteristics, somatic mutation, immune checkpoint and its ligand expression levels between high- and low-risk groups, as well as drug sensitivity prediction. ICGC cohort was set as the validation cohort. TCGA-LIHC cohort and three independent the Gene Expression Omnibus (GEO) datasets (GSE54236, GSE14520, and GSE64041) was used to verify IRDEGs expression, as well as PCR assays using clinical samples. ResultsThe IRDEGs was composed of 4 genes, namely B2M, SPP1, PPIA, and HRG. The 438 HCC patients were divided into high- and low-risk group. The high-risk group was associated with poor prognosis, including higher T stage, advanced pathological stages, less immune cell infiltration, higher TP53 mutation rate, the high expression of CTLA4 and HAVCR2. Besides, high-risk populations benefit from most chemotherapy drugs. Similarly, the performance of the risk model was validated in the ICGC. All four datasets (TCGA-LIHC cohort, GSE54236, GSE14520, and GSE64041) and clinical q-PCR results demonstrated that, compared with normal samples, the expressions of B2M and HRG were lower in tumor samples, and the expression of SPP1 was higher. ConclusionIn summary, the immune-related prognostic signature had a good predictive performance on prognosis and immunotherapy for HCC patients.

Single cell transcriptomic analysis reveals tumor immune infiltration by NK cells gene signature in lung adenocarcinoma

BackgroundNatural Killer (NK) cells are vital components of the innate immune system, crucial for combating infections and tumor growth, making them pivotal in cancer prognosis and immunotherapy. We sought to understand the diverse characteristics of NK cells within lung adenocarcinoma (LUAD) by conducting single-cell RNA sequencing analyses. MethodsUsing the scRNA-seq dataset for multiple primary lung cancers (MPLCs), we examined two major NK cell groups, NK1 and NK2, comparing the expression profiles of 422 differentially expressed NK signature genes. We identified eight genes (SPON2, PLEKHG3, CAMK2N1, RAB27B, CTBP2, EFHD2, GOLM1, and PLOD1) that distinguish NK1 from NK2 cells. A prognostic signature, the NK gene signature (NKGS) score, was established through LASSO Cox regression. High NKGS scores were linked to poorer overall survival in TCGA-LUAD patients and consistently validated in other datasets (GSE31210 and GSE14814). ResultsFunctional analysis revealed an enrichment of genes related to the TGF-β signaling pathway in the high NKGS score group. Moreover, a high NKGS score correlated with an immunosuppressive tumor microenvironment (TME) driven by immune evasion mechanisms. We also observed reduced T-cell receptor (TCR) repertoire diversity in the high-risk NKGS group, indicating a negative association between inflammation and risk score. ConclusionThis study introduced the innovative NKGS score, differentiating NK1 from NK2 cells. High NKGS scores were associated with the TGF-β pathway and provided insights into LUAD prognosis and immune activities.

DeepGSEA: explainable deep gene set enrichment analysis for single-cell transcriptomic data.

Gene set enrichment (GSE) analysis allows for an interpretation of gene expression through pre-defined gene set databases and is a critical step in understanding different phenotypes. With the rapid development of single-cell RNA sequencing (scRNA-seq) technology, GSE analysis can be performed on fine-grained gene expression data to gain a nuanced understanding of phenotypes of interest. However, with the cellular heterogeneity in single-cell gene profiles, current statistical GSE analysis methods sometimes fail to identify enriched gene sets. Meanwhile, deep learning has gained traction in applications like clustering and trajectory inference in single-cell studies due to its prowess in capturing complex data patterns. However, its use in GSE analysis remains limited, due to interpretability challenges. In this paper, we present DeepGSEA, an explainable deep gene set enrichment analysis approach which leverages the expressiveness of interpretable, prototype-based neural networks to provide an in-depth analysis of GSE. DeepGSEA learns the ability to capture GSE information through our designed classification tasks, and significance tests can be performed on each gene set, enabling the identification of enriched sets. The underlying distribution of a gene set learned by DeepGSEA can be explicitly visualized using the encoded cell and cellular prototype embeddings. We demonstrate the performance of DeepGSEA over commonly used GSE analysis methods by examining their sensitivity and specificity with four simulation studies. In addition, we test our model on three real scRNA-seq datasets and illustrate the interpretability of DeepGSEA by showing how its results can be explained. https://github.com/Teddy-XiongGZ/DeepGSEA.

Identifying cell states in single-cell RNA-seq data at statistically maximal resolution.

Single-cell RNA sequencing (scRNA-seq) has become a popular experimental method to study variation of gene expression within a population of cells. However, obtaining an accurate picture of the diversity of distinct gene expression states that are present in a given dataset is highly challenging because of the sparsity of the scRNA-seq data and its inhomogeneous measurement noise properties. Although a vast number of different methods is applied in the literature for clustering cells into subsets with 'similar' expression profiles, these methods generally lack rigorously specified objectives, involve multiple complex layers of normalization, filtering, feature selection, dimensionality-reduction, employ ad hoc measures of distance or similarity between cells, often ignore the known measurement noise properties of scRNA-seq measurements, and include a large number of tunable parameters. Consequently, it is virtually impossible to assign concrete biophysical meaning to the clusterings that result from these methods. Here we address the following problem: Given raw unique molecule identifier (UMI) counts of an scRNA-seq dataset, partition the cells into subsets such that the gene expression states of the cells in each subset are statistically indistinguishable, and each subset corresponds to a distinct gene expression state. That is, we aim to partition cells so as to maximally reduce the complexity of the dataset without removing any of its meaningful structure. We show that, given the known measurement noise structure of scRNA-seq data, this problem is mathematically well-defined and derive its unique solution from first principles. We have implemented this solution in a tool called Cellstates which operates directly on the raw data and automatically determines the optimal partition and cluster number, with zero tunable parameters. We show that, on synthetic datasets, Cellstates almost perfectly recovers optimal partitions. On real data, Cellstates robustly identifies subtle substructure within groups of cells that are traditionally annotated as a common cell type. Moreover, we show that the diversity of gene expression states that Cellstates identifies systematically depends on the tissue of origin and not on technical features of the experiments such as the total number of cells and total UMI count per cell. In addition to the Cellstates tool we also provide a small toolbox of software to place the identified cellstates into a hierarchical tree of higher-order clusters, to identify the most important differentially expressed genes at each branch of this hierarchy, and to visualize these results.

Role of TAP1 in the identification of immune-hot tumor microenvironment and its prognostic significance for immunotherapeutic efficacy in gastric carcinoma.

Gastric cancer (GC), a multifaceted gastrointestinal malignancy, is the fourth most prevalent contributor to cancer-related fatalities globally. As a member of the ATP-binding cassette (ABC) family, transporter associated with antigen processing 1 (TAP1) is crucial for conveying antigen peptides from the cytoplasm to the lumen of the endoplasmic reticulum and subsequently loading them onto the major histocompatibility complex (MHC) class I molecules. Recent studies have established the biological significance of TAP1 in upholding tumor survival and facilitating immune evasion by remodeling the tumor microenvironment (TME) and orchestrating immune infiltration. The study was conducted to elucidate the association of TAP1 expression with immunological characteristics, and sought to exploit the value of TAP1 as a biomarker reflecting the inflamed TME and immunotherapeutic response. RNA-sequencing profiles and clinical annotations were obtained from The Cancer Genome Atlas-stomach adenocarcinoma (TCGA-STAD) cohort and Gene Expression Omnibus (GEO) portal. Preprocessing was conducting using the limma package. Weighted gene co-expression network analysis (WGCNA) was used to identify gene modules and TAP1 co-expressed genes (CEGs) based on correlation patterns. Consensus clustering and silhouette analysis determined the optimal number of TAP1-related groups. Gene expression profiles were integrated and classified using the pamr package. The Estimation of STromal and Immune cells in MAlignant Tumors using Expression data (ESTIMATE) algorithm and single-sample gene set enrichment analysis (ssGSEA) were used to evaluate immunological characteristics. Differential expression analysis was conducted using the limma package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Single-cell RNA sequencing (scRNA-seq) datasets were analyzed using the Seurat toolkit to characterize cell types. Within this investigation, no significant differences in TAP1 expression were observed among patients exhibiting various clinicopathological features, indicating that TAP1 expression was not specific to molecular subtypes. Subsequent analysis revealed a positive correlation between TAP1 and diverse immunological traits, encompassing immunomodulators, tumor-infiltrating immune cells, as well as immune checkpoints across multiple datasets. Besides, within a GC immunotherapy cohort, individuals displaying high TAP1 expression demonstrated an increased likelihood of achieving complete remission (CR) post-treatment, suggesting heightened sensitivity to immunotherapy. In the clinical cohort, TAP1 overexpression in GC patients was positively correlated with CD8. TAP1 appears linked to an inflamed TME and serves as a prospective biomarker for discerning immunological attributes and gauging immunotherapeutic responses in GC, particularly in identifying immune-reactive tumors.