Ovarian cancer-specific dysregulated genes with prognostic significance: scRNA-Seq with bulk RNA-Seq data and experimental validation

Identification of a G-protein coupled receptor-related gene signature through bioinformatics analysis to construct a risk model for ovarian cancer prognosis

BACKGROUND

Ovarian cancer (OV) is a common malignant tumor of the female reproductive system with a 5-year survival rate of ∼30 %. Inefficient early diagnosis and prognosis leads to poor survival in most patients. G protein-coupled receptors (GPCRs, the largest family of human cell surface receptors) are associated with OV. We aimed to identify GPCR-related gene (GPCRRG) signatures and develop a novel model to predict OV prognosis.

METHOD

We downloaded data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Prognostic GPCRRGs were screened using least absolute shrinkage and selection operator (LASSO) Cox regression analysis, and a prognostic model was constructed. The predictive ability of the model was evaluated by Kaplan-Meier (K-M) survival analysis. The levels of GPCRRGs were examined in normal and OV cell lines using quantitative reverse-Etranscription polymerase chain reaction. The immunological characteristics of the high- and low-risk groups were analyzed using single-sample gene set enrichment analysis (ssGSEA) and CIBERSORT.

RESULTS

Based on the risks scores, 17 GPCRRGs were associated with OV prognosis. CXCR4, GPR34, LGR6, LPAR3, and RGS2 were significantly expressed in three OV datasets and enabled accurate OV diagnosis. K-M analysis of the prognostic model showed that it could differentiate high- and low-risk patients, which correspond to poorer and better prognoses, respectively. GPCRRG expression was correlated with immune infiltration rates.

CONCLUSIONS

Our prognostic model elaborates on the roles of GPCRRGs in OV and provides a new tool for prognosis and immune response prediction in patients with OV.

The expression of immune related genes and potential regulatory mechanisms in schizophrenia

OBJECTIVE

This study aimed to investigate the role of immune dysfunction in the pathogenesis of schizophrenia through single-cell transcriptome and bulk RNA data analyses.

METHODS

The single-cell RNA sequencing (scRNA-seq) was selected to assess the cellular composition and gene expression profiles of the brain tissue. Further, bulk RNA sequencing data was utilized to corroborate findings from the single-cell analyses and provide additional insights into the molecular changes associated with the disease. Gene-drug interaction data was also included to identify potential therapeutic drugs targeting these dysregulated immune-related genes in schizophrenia.

RESULTS

We discovered significant differences in cellular composition within schizophrenia tissue, including increased infiltration of fibroblasts, horizontal basal cells, monocytes, mesenchymal cells, and smooth muscle cells. The investigation of immune-related genes revealed significantly different expression of genes such as S100A2, CCL14, IGHA1, BPIFA1, GDF15, IL32, BPIFB2, HLA-DRA, S100A8, PTX3, TPM2, TNFRSF12A, GREM1 and others. These genes possibly contribute to the progression of schizophrenia through various pathways such as humoral immune response, IL-17 signaling pathway, adaptive immune response, antigen processing and presentation, and gut IgA production. Our findings also suggest possible transcriptional regulation in schizophrenia's immune dysfunction by transcription factors in monocytes, neutrophils, endothelial cells, and epithelial cells. Lastly, potential therapeutic drugs were identified through gene-drug interaction data, such as those targeting HLA-A and HLAB.

CONCLUSION

The cellular heterogeneity and immune-related gene dysregulation play important roles in schizophrenia, which provides a foundation for understanding the pathogenesis and developing new treatment methods.

A novel defined programmed cell death related gene signature for predicting the prognosis of serous ovarian cancer

PURPOSE

This study aims to explore the contribution of differentially expressed programmed cell death genes (DEPCDGs) to the heterogeneity of serous ovarian cancer (SOC) through single-cell RNA sequencing (scRNA-seq) and assess their potential as predictors for clinical prognosis.

METHODS

SOC scRNA-seq data were extracted from the Gene Expression Omnibus database, and the principal component analysis was used for cell clustering. Bulk RNA-seq data were employed to analyze SOC-associated immune cell subsets key genes. CIBERSORT and single-sample gene set enrichment analysis (ssGSEA) were utilized to calculate immune cell scores. Prognostic models and nomograms were developed through univariate and multivariate Cox analyses.

RESULTS

Our analysis revealed that 48 DEPCDGs are significantly correlated with apoptotic signaling and oxidative stress pathways and identified seven key DEPCDGs (CASP3, GADD45B, GNA15, GZMB, IL1B, ISG20, and RHOB) through survival analysis. Furthermore, eight distinct cell subtypes were characterized using scRNA-seq. It was found that G protein subunit alpha 15 (GNA15) exhibited low expression across these subtypes and a strong association with immune cells. Based on the DEGs identified by the GNA15 high- and low-expression groups, a prognostic model comprising eight genes with significant prognostic value was constructed, effectively predicting patient overall survival. Additionally, a nomogram incorporating the RS signature, age, grade, and stage was developed and validated using two large SOC datasets.

CONCLUSION

GNA15 emerged as an independent and excellent prognostic marker for SOC patients. This study provides valuable insights into the prognostic potential of DEPCDGs in SOC, presenting new avenues for personalized treatment strategies.

SLAMF7 predicts prognosis and correlates with immune infiltration in serous ovarian carcinoma

OBJECTIVE

Signaling lymphocytic activation molecule family members (SLAMFs) play a critical role in immune regulation of malignancies. This study aims to investigate the prognostic value and function of SLAMFs in ovarian cancer (OC).

METHODS

The expression analysis of SLAMFs was conducted based on The Cancer Genome Atlas Ovarian Cancer Collection (TCGA-OV) and Gene Expression Omnibus (GEO) databases. Immunohistochemistry (IHC) was further performed on tissue arrays (n=98) to determine the expression of SLAMF7. Kaplan-Meier plotter and multivariate Cox regression model were used to evaluate the correlation of SLAMF7 expression with survival outcomes of patients. The molecular function of SLAMF7 in OC was further investigated using Gene Set Enrichment Analysis (GSEA).

RESULTS

SLAMF7 mRNA expression were significantly upregulated in OC tumor tissue compared to normal tissue. IHC revealed that SLAMF7 expression was located in the interstitial parts of tumor tissue, and higher SLAMF7 expression was associated with favorable survival outcomes. GSEA demonstrated that SLAMF7 is involved immune-related pathways. Further analysis showed that SLAMF7 had a strong correlation with the T cell-specific biomarker (CD3) but not with the B cell (CD19, CD22, and CD23) and natural killer cell-specific biomarkers (CD85C, CD336, and CD337). Furthermore, IHC analysis confirmed that SLAMF7 was expressed in tumor-infiltrating T cells, and the IHC score of SLAMF7 was positively correlated with CD3 (r=0.85, p<0.001).

CONCLUSION

SLAMF7 is expressed in the interstitial components of clinical OC tissue, and higher SLAMF7 expression indicated a favorable prognosis for patients with OC. Additionally, SLAMF7 is involved in T-cell immune infiltration in OC.

Identification of lncRNAs Deregulated in Epithelial Ovarian Cancer Based on a Gene Expression Profiling Meta-Analysis

Epithelial ovarian cancer (EOC) is one of the deadliest gynecological cancers worldwide, mainly because of its initially asymptomatic nature and consequently late diagnosis. Long non-coding RNAs (lncRNA) are non-coding transcripts of more than 200 nucleotides, whose deregulation is involved in pathologies such as EOC, and are therefore envisaged as future biomarkers. We present a meta-analysis of available gene expression profiling (microarray and RNA sequencing) studies from EOC patients to identify lncRNA genes with diagnostic and prognostic value. In this meta-analysis, we include 46 independent cohorts, along with available expression profiling data from EOC cell lines. Differential expression analyses were conducted to identify those lncRNAs that are deregulated in (i) EOC versus healthy ovary tissue, (ii) unfavorable versus more favorable prognosis, (iii) metastatic versus primary tumors, (iv) chemoresistant versus chemosensitive EOC, and (v) correlation to specific histological subtypes of EOC. From the results of this meta-analysis, we established a panel of lncRNAs that are highly correlated with EOC. The panel includes several lncRNAs that are already known and even functionally characterized in EOC, but also lncRNAs that have not been previously correlated with this cancer, and which are discussed in relation to their putative role in EOC and their potential use as clinically relevant tools.

Molecular cluster mining of high-grade serous ovarian cancer via multi-omics data analysis aids precise medicine

PURPOSE

HGSOC is a kind of gynecological cancer with high mortality and strong heterogeneity. The study used multi-omics and multiple algorithms to identify novel molecular subtypes, which can help patients obtain more personalized treatments.

METHODS

Firstly, the consensus clustering result was obtained using a consensus ensemble of ten classical clustering algorithms, based on mRNA, lncRNA, DNA methylation, and mutation data. The difference in signaling pathways was evaluated using the single-sample gene set enrichment analysis (ssGSEA). Meanwhile, the relationship between genetic alteration, response to immunotherapy, drug sensitivity, prognosis, and subtypes was further analyzed. Finally, the reliability of the new subtype was verified in three external datasets.

RESULTS

Three molecular subtypes were identified. Immune desert subtype (CS1) had little enrichment in the immune microenvironment and metabolic pathways. Immune/non-stromal subtype (CS2) was enriched in the immune microenvironment and metabolism of polyamines. Immune/stromal subtype (CS3) not only enriched anti-tumor immune microenvironment characteristics but also enriched pro-tumor stroma characteristics, glycosaminoglycan metabolism, and sphingolipid metabolism. The CS2 had the best overall survival and the highest response rate to immunotherapy. The CS3 had the worst prognosis and the lowest response rate to immunotherapy but was more sensitive to PARP and VEGFR molecular-targeted therapy. The similar differences among three subtypes were successfully validated in three external cohorts.

CONCLUSION

We used ten clustering algorithms to comprehensively analyze four types of omics data, identified three biologically significant subtypes of HGSOC patients, and provided personalized treatment recommendations for each subtype. Our findings provided novel views into the HGSOC subtypes and could provide potential clinical treatment strategies.

Development and Validation of a Prognostic Risk Model Based on Nature Killer Cells for Serous Ovarian Cancer

Nature killer (NK) cells are increasingly considered important in tumor microenvironment, but their role in predicting the prognosis of ovarian cancer has not been revealed. This study aimed to develop a prognostic risk model for ovarian cancer based on NK cells. Firstly, differentially expressed genes (DEGs) of NK cells were found by single-cell RNA-sequencing dataset analysis. Based on six NK-cell DEGs identified by univariable, Lasso and multivariable Cox regression analyses, a prognostic risk model for serous ovarian cancer was developed in the TCGA cohort. This model was then validated in three external cohorts, and evaluated as an independent prognostic factor by multivariable Cox regression analysis together with clinical characteristics. With the investigation of the underlying mechanism, a relation between a higher risk score of this model and more immune activities in tumor microenvironment was revealed. Furthermore, a detailed inspection of infiltrated immunocytes indicated that not only quantity, but also the functional state of these immunocytes might affect prognostic risk. Additionally, the potential of this model to predict immunotherapeutic response was exhibited by evaluating the functional state of cytotoxic T lymphocytes. To conclude, this study introduced a novel prognostic risk model based on NK-cell DEGs, which might provide assistance for the personalized management of serous ovarian cancer patients.

Targeting tumor-associated macrophages for successful immunotherapy of ovarian carcinoma

Epithelial ovarian cancer (EOC) is among the top five causes of cancer-related death in women, largely reflecting early, prediagnosis dissemination of malignant cells to the peritoneum. Despite improvements in medical therapies, particularly with the implementation of novel drugs targeting homologous recombination deficiency, the survival rates of patients with EOC remain low. Unlike other neoplasms, EOC remains relatively insensitive to immune checkpoint inhibitors, which is correlated with a tumor microenvironment (TME) characterized by poor infiltration by immune cells and active immunosuppression dominated by immune components with tumor-promoting properties, especially tumor-associated macrophages (TAMs). In recent years, TAMs have attracted interest as potential therapeutic targets by seeking to reverse the immunosuppression in the TME and enhance the clinical efficacy of immunotherapy. Here, we review the key biological features of TAMs that affect tumor progression and their relevance as potential targets for treating EOC. We especially focus on the therapies that might modulate the recruitment, polarization, survival, and functional properties of TAMs in the TME of EOC that can be harnessed to develop superior combinatorial regimens with immunotherapy for the clinical care of patients with EOC.

Focus on the molecular mechanisms of cisplatin resistance based on multi-omics approaches

Cisplatin is commonly used in combination with other cytotoxic agents as a standard treatment regimen for a variety of solid tumors, such as lung, ovarian, testicular, and head and neck cancers. However, the effectiveness of cisplatin is accompanied by toxic side effects, for instance, nephrotoxicity and neurotoxicity. The response of tumors to cisplatin treatment involves multiple physiological processes, and the efficacy of chemotherapy is limited by the intrinsic and acquired resistance of tumor cells. Although enormous efforts have been made toward molecular mechanisms of cisplatin resistance, the development of omics provides new insights into the understanding of cisplatin resistance at genome, transcriptome, proteome, metabolome and epigenome levels. Mechanism studies using different omics approaches revealed the necessity of multi-omics applications, which provide information at different cellular function levels and expand our recognition of the peculiar genetic and phenotypic heterogeneity of cancer. The present work systematically describes the underlying mechanisms of cisplatin resistance in different tumor types using multi-omics approaches. In addition to the classical mechanisms such as enhanced drug efflux, increased DNA damage repair and changes in the cell cycle and apoptotic pathways, other changes like increased protein damage clearance, increased protein glycosylation, enhanced glycolytic process, dysregulation of the oxidative phosphorylation pathway, ferroptosis suppression and mRNA m6A methylation modification can also induce cisplatin resistance. Therefore, utilizing the integrated omics to identify key signaling pathways, target genes and biomarkers that regulate chemoresistance are essential for the development of new drugs or strategies to restore tumor sensitivity to cisplatin.

Integration of bulk RNA sequencing and single-cell analysis reveals a global landscape of DNA damage response in the immune environment of Alzheimer's disease

Background

We developed a novel system for quantifying DNA damage response (DDR) to help diagnose and predict the risk of Alzheimer's disease (AD).

Methods

We thoroughly estimated the DDR patterns in AD patients Using 179 DDR regulators. Single-cell techniques were conducted to validate the DDR levels and intercellular communications in cognitively impaired patients. The consensus clustering algorithm was utilized to group 167 AD patients into diverse subgroups after a WGCNA approach was employed to discover DDR-related lncRNAs. The distinctions between the categories in terms of clinical characteristics, DDR levels, biological behaviors, and immunological characteristics were evaluated. For the purpose of choosing distinctive lncRNAs associated with DDR, four machine learning algorithms, including LASSO, SVM-RFE, RF, and XGBoost, were utilized. A risk model was established based on the characteristic lncRNAs.

Results

The progression of AD was highly correlated with DDR levels. Single-cell studies confirmed that DDR activity was lower in cognitively impaired patients and was mainly enriched in T cells and B cells. DDR-related lncRNAs were discovered based on gene expression, and two different heterogeneous subtypes (C1 and C2) were identified. DDR C1 belonged to the non-immune phenotype, while DDR C2 was regarded as the immune phenotype. Based on various machine learning techniques, four distinctive lncRNAs associated with DDR, including FBXO30-DT, TBX2-AS1, ADAMTS9-AS2, and MEG3 were discovered. The 4-lncRNA based riskScore demonstrated acceptable efficacy in the diagnosis of AD and offered significant clinical advantages to AD patients. The riskScore ultimately divided AD patients into low- and high-risk categories. In comparison to the low-risk group, high-risk patients showed lower DDR activity, accompanied by higher levels of immune infiltration and immunological score. The prospective medications for the treatment of AD patients with low and high risk also included arachidonyltrifluoromethane and TTNPB, respectively.

Conclusions

In conclusion, immunological microenvironment and disease progression in AD patients were significantly predicted by DDR-associated genes and lncRNAs. A theoretical underpinning for the individualized treatment of AD patients was provided by the suggested genetic subtypes and risk model based on DDR.

Establishment of an ovarian cancer omentum metastasis-related prognostic model by integrated analysis of scRNA-seq and bulk RNA-seq

OBJECTIVE

Ovarian cancer has the highest mortality rate among gynecological malignant tumors, and it preferentially metastasizes to omental tissue, leading to intestinal obstruction and death. scRNA-seq is a powerful technique to reveal tumor heterogeneity. Analyzing omentum metastasis of ovarian cancer at the single-cell level may be more conducive to exploring and understanding omentum metastasis and prognosis of ovarian cancer at the cellular function and genetic levels.

METHODS

The omentum metastasis site scRNA-seq data of GSE147082 were acquired from the GEO (Gene Expression Omnibus) database, and single cells were clustered by the Seruat package and annotated by the SingleR package. Cell differentiation trajectories were reconstructed through the monocle package. The ovarian cancer microarray data of GSE132342 were downloaded from GEO and were clustered by using the ConsensusClusterPlus package into omentum metastasis-associated clusters according to the marker genes gained from single-cell differentiation trajectory analysis. The tumor microenvironment (TME) and immune infiltration differences between clusters were analyzed by the estimate and CIBERSORT packages. The expression matrix of genes used to cluster GSE132342 patients was extracted from bulk RNA-seq data of TCGA-OV (The Cancer Genome Atlas ovarian cancer), and least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression were performed to establish an omentum metastasis-associated gene (OMAG) signature. The signature was then tested by GSE132342 data. Finally, the clinicopathological characteristics of TCGA-OV were screened by univariate and multivariate Cox regression analysis to draw the nomogram.

RESULTS

A total of 9885 cells from 6 patients were clustered into 18 cell clusters and annotated into 14 cell types. Reconstruction of differentiation trajectories divided the cells into 5 branches, and a total of 781 cell trajectory-related characteristic genes were obtained. A total of 3769 patients in GSE132342 were subtyped into 3 clusters by 74 cell trajectory-related characteristic genes. Kaplan-Meier (K-M) survival analysis showed that the prognosis of cluster 2 was the worst, P < 0.001. The TME analysis showed that the ESTIMATE score and stromal score in cluster 2 were significantly higher than those in the other two clusters, P < 0.001. The immune infiltration analysis showed differences in the fraction of 8 immune cells among the 3 clusters, P < 0.05. The expression data of 74 genes used for GEO clustering were extracted from 379 patients in TCGA-OV, and combined with survival information, 10 candidates for OMAGs were filtered by LASSO. By using multivariate Cox regression, the 6-OMAGs signature was established as RiskScore = 0.307*TIMP3 + 3.516*FBN1-0.109*IGKC + 0.209*RPL21 + 0.870*UCHL1 + 0.365*RARRES1. Taking TCGA-OV as the training set and GSE132342 as the test set, receiver operating characteristic (ROC) curves were drawn to verify the prognostic value of 6-OMAGs. Screened by univariate and multivariate Cox regression analysis, 3 (age, cancer status, primary therapy outcome) of 5 clinicopathological characteristics were used to construct the nomogram combined with risk score.

CONCLUSION

We constructed an ovarian cancer prognostic model related to omentum metastasis composed of 6-OMAGs and 3 clinicopathological features and analyzed the potential mechanism of these 6-OMAGs in ovarian cancer omental metastasis.

Identification of a Recurrence Gene Signature for Ovarian Cancer Prognosis by Integrating Single-Cell RNA Sequencing and Bulk Expression Datasets

Ovarian cancer is one of the most common gynecological malignancies in women, with a poor prognosis and high mortality. With the expansion of single-cell RNA sequencing technologies, the inner biological mechanism involved in tumor recurrence should be explored at the single-cell level, and novel prognostic signatures derived from recurrence events were urgently identified. In this study, we identified recurrence-related genes for ovarian cancer by integrating two Gene Expression Omnibus datasets, including an ovarian cancer single-cell RNA sequencing dataset (GSE146026) and a bulk expression dataset (GSE44104). Based on these recurrence genes, we further utilized the merged expression dataset containing a total of 524 ovarian cancer samples to identify prognostic signatures and constructed a 13-gene risk model, named RMGS (recurrence marker gene signature). Based on the RMGS score, the samples were stratified into high-risk and low-risk groups, and these two groups displayed significant survival difference in two independent validation cohorts including The Cancer Genome Atlas (TCGA). Also, the RMGS score remained significantly independent in multivariate analysis after adjusting for clinical factors, including the tumor grade and stage. Furthermore, there existed close associations between the RMGS score and immune characterizations, including checkpoint inhibition, EMT signature, and T-cell infiltration. Finally, the associations between RMGS scores and molecular subtypes revealed that samples with mesenchymal subtypes displayed higher RMGS scores. In the meanwhile, the genomics characterization from these two risk groups was also identified. In conclusion, the recurrence-related RMGS model we identified could provide a new understanding of ovarian cancer prognosis at the single-cell level and offer a reference for therapy decisions for patient treatment.