Gene expression analysis in ovarian cancer - faults and hints from DNA microarray study

Deciphering Folate Receptor alphaGene Expression and mRNA Signatures in Ovarian Cancer: Implications for Precision Therapies

Antibody-drug conjugates targeting folate receptor alpha (FRα) are a promising treatment for platinum-resistant ovarian cancer (OC) with high FRα expression. Challenges persist in accurately assessing FRα expression levels. Our study aimed to better elucidate FRα gene expression and identify mRNA signatures in OC. We pooled OC gene expression data from 16 public datasets, encompassing 1832 OC and 30 normal ovarian tissues. Additional data included DNA copy number and methylation data from TCGA and protein data from 363 cancer cell lines from the Broad Institute Cancer Cell Line Encyclopedia. FOLR1 mRNA expression was significantly correlated with protein expression in pan-cancer cell lines and ovarian cancer cell lines. FOLR1 expression was higher in OC samples than in normal ovarian tissues (OR = 3.88, p = 6.97 × 10-12). Patients with high FOLR1 expression were more likely to be diagnosed with serous histology, FIGO stage III-IV, and high-grade tumors; however, nearly similar percentages of patients with low FOLR1 expression were also diagnosed with these features. FOLR1 mRNA expression was not correlated with platinum sensitivity or complete surgery, nor with prognosis. However, we identified a 187-gene signature associated with high FOLR1 expression that was significantly associated with improved survival (HR = 0.71, p = 1.18 × 10-6), independently from clinicopathological features. We identified a gene expression signature correlated to high FRα expression and OC prognosis, which may be used to refine therapeutic strategies targeting FRα in OC. These findings warrant validation in larger cohorts.

Targeting the GPI transamidase subunit GPAA1 abrogates the CD24 immune checkpoint in ovarian cancer

CD24 is frequently overexpressed in ovarian cancer and promotes immune evasion by interacting with its receptor Siglec10, present on tumor-associated macrophages, providing a "don't eat me" signal that prevents targeting and phagocytosis by macrophages. Factors promoting CD24 expression could represent novel immunotherapeutic targets for ovarian cancer. Here, using a genome-wide CRISPR knockout screen, we identify GPAA1 (glycosylphosphatidylinositol anchor attachment 1), a factor that catalyzes the attachment of a glycosylphosphatidylinositol (GPI) lipid anchor to substrate proteins, as a positive regulator of CD24 cell surface expression. Genetic ablation of GPAA1 abolishes CD24 cell surface expression, enhances macrophage-mediated phagocytosis, and inhibits ovarian tumor growth in mice. GPAA1 shares structural similarities with aminopeptidases. Consequently, we show that bestatin, a clinically advanced aminopeptidase inhibitor, binds to GPAA1 and blocks GPI attachment, resulting in reduced CD24 cell surface expression, increased macrophage-mediated phagocytosis, and suppressed growth of ovarian tumors. Our study highlights the potential of targeting GPAA1 as an immunotherapeutic approach for CD24+ ovarian cancers.

Low FHL1 expression indicates a good prognosis and drug sensitivity in ovarian cancer

Chemotherapy resistance is the main reason for the poor prognosis of ovarian cancer (OC). FHL1 is an important tumour regulator, but its relationship with the prognosis, drug resistance, and tumour microenvironment of OC is unknown. Immunohistochemistry was used to determine FHL1 expression in OC. Kaplan‒Meier plotter was used for survival analysis. The value of gene expression in predicting drug resistance was estimated using the area under the curve (AUC). Bivariate correlation was used to determine the coexpression of two genes. Functional cluster and pathway enrichment were used to uncover hidden signalling pathways. The relationship between gene levels and the tumour microenvironment was visualised through the ggstatsplot and pheatmap packages. The mRNA and protein levels of FHL1 were downregulated in 426 and 100 OC tissues, respectively. Low FHL1 expression was correlated with good progression-free survival (PFS), postprogression survival, and overall survival (OS) in 1815 OC patients, and was further confirmed to be associated with good OS by immunohistochemistry in 152 OC tissues. Furthermore, FHL1 was downregulated in drug-sensitive tissues, while its high expression predicted drug resistance (AUC > 0.65). Mechanistically, FHL1 was coexpressed with FLNC, CAV1, PPP1R12B, and FLNA at the mRNA and protein levels in 558 and 174 OC tissues, respectively, and their expression was downregulated in OC. Additionally, very strong coexpression of FHL1 with the four genes was identified in at least 23 different tumours. Low expression of the four genes was associated with good PFS, and the combination of FHL1 with the four genes provided better prognostic power. Meanwhile, the expression of all five genes was strongly and positively associated with the abundance of macrophages. Low FHL1 expression acts as a favourable factor in OC, probably via positive coexpression with FLNC, CAV1, PPP1R12B, and FLNA.

Signatures of tumor-associated macrophages correlate with treatment response in ovarian cancer patients

Ovarian cancer (OC) ranks as the second leading cause of death among gynecological cancers. Numerous studies have indicated a correlation between the tumor microenvironment (TME) and the clinical response to treatment in OC patients. Tumor-associated macrophages (TAMs), a crucial component of the TME, exert influence on invasion, metastasis, and recurrence in OC patients. To delve deeper into the role of TAMs in OC, this study conducted an extensive analysis of single-cell data from OC patients. The aim is to develop a new risk score (RS) to characterize the response to treatment in OC patients to inform clinical treatment. We first identified TAM-associated genes (TAMGs) in OC patients and examined the protein and mRNA expression levels of TAMGs by Western blot and PCR experiments. Additionally, a scoring system for TAMGs was constructed, successfully categorizing patients into high and low RS subgroups. Remarkably, significant disparities were observed in immune cell infiltration and immunotherapy response between the high and low RS subgroups. The findings revealed that patients in the high RS group had a poorer prognosis but displayed greater sensitivity to immunotherapy. Another important finding was that patients in the high RS subgroup had a higher IC50 for chemotherapeutic agents. Furthermore, further experimental investigations led to the discovery that THEMIS2 could serve as a potential target in OC patients and is associated with EMT (epithelial-mesenchymal transition). Overall, the TAMGs-based scoring system holds promise for screening patients who would benefit from therapy and provides valuable information for the clinical treatment of OC.

Comprehensive analysis based on glycolytic and glutaminolytic pathways signature for predicting prognosis and immunotherapy in ovarian cancer

Background: Our study attempts to develop and identify an aerobic glycolysis and glutamine-related genes (AGGRGs) signature for estimating prognostic effectively of ovarian cancer (OV) patients. Materials & methods: OV related data were extracted from the multiple public databases, including TCGA-OV, GSE26193, GSE63885, and ICGC-OV. A consistent clustering approach was used to characterize the subtypes associated with AGGRGs. LASSO Cox regressions was utilized to construct the prognosis signatures of AGGRGs. In addition, GSE26193, GSE63885 and ICGC-OV served as independent external cohorts to assess the reliability of the model. In vitro and in vivo experiments were conducted to study the role of AAK1 in the malignant progression and glutamine metabolism of OV, and assessed its therapeutic potential for treating OV patients. Results: OV patients could be separated into four subtypes (quiescent, glycolysis, glutaminolytic, and mixed subtypes). The survival outcome of glutaminolytic subtype was notably worse than the glycolytic subtype. Besides, we identified eight AGGRGs (AAK1, GJB6, HMGN5, LPIN3, INTS6L, PPOX, SPAG4, and ZNF316) to establish a prognostic signature for OV patients. Comprehensive analysis revealed that the signature risk score served as an independent prognostic factor for OV. Additionally, high-risk OV patients were less sensitive to platinum and, conversely, were proved to be more responsive to immunotherapy than low-risk score. In cytological experiments, we found that AAK1 could promote cancer progression and glutamine metabolism via activating the Notch3 pathway in OV cells. Furthermore, knockdown of AAK1 significantly inhibited tumor growth and weight, decreased lung metastases, and ultimately extended the survival time of the nude mice. Conclusions: The prognostic signature of AGGRGs constructed could efficiently estimate the prognosis and immunotherapy effectiveness of OV patients. In addition, AAK1 may represent a promising therapeutic target for OV.

BDP1 as a biomarker in serous ovarian cancer

BACKGROUND

TFIIIB, an RNA polymerase III specific transcription factor has been found to be deregulated in human cancers with much of the research focused on the TBP, BRF1, and BRF2 subunits. To date, the TFIIIB specific subunit BDP1 has not been investigated in ovarian cancer but has previously been shown to be deregulated in neuroblastoma, breast cancer, and Non-Hodgkins lymphoma.

RESULTS

Using in silico analysis of clinically derived platforms, we report a decreased BDP1 expression as a result of deletion in serous ovarian cancer and a correlation with higher and advanced ovarian stages. Further analysis in the context of TP53 mutations, a major contributor to ovarian tumorigenesis, suggests that high BDP1 expression is unfavorable for overall survival and high BDP1 expression occurs in stages 2, 3 and 4 serous ovarian cancer. Additionally, high BDP1 expression is disadvantageous and unfavorable for progression-free survival. Lastly, BDP1 expression significantly decreased in patients treated with first-line chemotherapy, platin and taxane, at twelve-month relapse-free survival.

CONCLUSIONS

Taken together with a ROC analysis, the data suggest BDP1 could be of clinical relevance as a predictive biomarker in serous ovarian cancer. Lastly, this study further demonstrates that both the over- and under expression of BDP1 warrants further investigation and suggests BDP1 may exhibit dual function in the context of tumorigenesis.

CD4+ conventional T cells-related genes signature is a prognostic indicator for ovarian cancer

Introduction

It is believed that ovarian cancer (OC) is the most deadly form of gynecological cancer despite its infrequent occurrence, which makes it one of the most salient public health concerns. Clinical and preclinical studies have revealed that intratumoral CD4+ T cells possess cytotoxic capabilities and were capable of directly killing cancer cells. This study aimed to identify the CD4+ conventional T cells-related genes (CD4TGs) with respect to the prognosis in OC.

Methods

We obtained the transcriptome and clinical data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CD4TGs were first identified from single-cell datasets, then univariate Cox regression was used to screen prognosis-related genes, LASSO was conducted to remove genes with coefficient zero, and multivariate Cox regression was used to calculate riskscore and to construct the CD4TGs risk signature. Kaplan-Meier analysis, univariate Cox regression, multivariate Cox regression, time-dependent receiver operating characteristics (ROC), decision curve analysis (DCA), nomogram, and calibration were made to verify and evaluate the risk signature. Gene set enrichment analyses (GSEA) in risk groups were conducted to explore the tightly correlated pathways with the risk group. The role of riskscore has been further explored in the tumor microenvironment (TME), immunotherapy, and chemotherapy. A risk signature with 11 CD4TGs in OC was finally established in the TCGA database and furtherly validated in several GEO cohorts.

Results

High riskscore was significantly associated with a poorer prognosis and proven to be an independent prognostic biomarker by multivariate Cox regression. The 1-, 3-, and 5-year ROC values, DCA curve, nomogram, and calibration results confirmed the excellent prediction power of this model. Compared with the reported risk models, our model showed better performance. The patients were grouped into high-risk and low-risk subgroups according to the riskscore by the median value. The low-risk group patients tended to exhibit a higher immune infiltration, immune-related gene expression and were more sensitive to immunotherapy and chemotherapy.

Discussion

Collectively, our findings of the prognostic value of CD4TGs in prognosis and immune response, provided valuable insights into the molecular mechanisms and clinical management of OC.

Microfibril Associated Protein 5 (MFAP5) Is Related to Survival of Ovarian Cancer Patients but Not Useful as a Prognostic Biomarker

Ovarian cancer (OC) is usually diagnosed late due to its nonspecific symptoms and lack of reliable tools for early diagnostics and screening. OC studies concentrate on the search for new biomarkers and therapeutic targets. This study aimed to validate the MFAP5 gene, and its encoded protein, as a potential prognostic biomarker. In our previous study, we found that patients with high-grade serous OC who had higher MFAP5 mRNA levels had shorter survival, as compared with those with lower levels. Here, we used the Kaplan-Meier Plotter and CSIOVDB online tools to analyze possible associations of MFAP5 expression with survival and other clinico-pathological features. In these analyses, higher MFAP5 mRNA expression was observed in the more advanced FIGO stages and high-grade tumors, and was significantly associated with shorter overall and progression-free survival. Next, we analyzed the expression of the MFAP5 protein by immunohistochemistry (IHC) in 108 OC samples and tissue arrays. Stronger MFAP5 expression was associated with stronger desmoplastic reaction and serous vs. non-serous histology. We found no significant correlation between IHC results and survival, although there was a trend toward shorter survival in patients with the highest IHC scores. We searched for co-expressed genes/proteins using cBioPortal and analyzed potential MFAP5 interaction networks with the STRING tool. MFAP5 was shown to interact with many extracellular matrix proteins, and was connected to the Notch signaling pathway. Therefore, although not suitable as a prognostic biomarker for evaluation with a simple diagnostic tool like IHC, MFAP5 is worth further studies as a possible therapeutic target.

Identification and validation of an immune-related lncRNAs signature to predict the overall survival of ovarian cancer

Ovarian cancer (OC) is the most lethal gynecological cancer in women. Studies had reported that immune-related lncRNAs signatures were valuable in predicting the survival and prognosis of patients with various cancers. In our study, the prognostic value of immune-related lncRNAs was investigated in OC patients from TCGA-RNA-seq cohort (n=378) and HG-U133_Plus_2 cohort (n=590), respectively. Pearson correlation analysis was implemented to screen the immune-related lncRNA and then univariate Cox regression analysis was performed to explore their prognostic value in OC patients. Five prognostic immune-related lncRNAs were identified as prognostic lncRNAs. Besides, they were inputted into a LASSO Cox regression to establish and validate an immune-related lncRNA prognostic signature in TCGA-RNA-Seq cohort and HG-U133_Plus_2 cohort, respectively. Based on the best cut-off value of risk score, patients were divided into high- and low-risk groups. Survival analysis suggested that patients in the high-risk group had a worse overall survival (OS) than those in the low-risk group in both cohorts. The association between clinicopathological feathers and risk score was then evaluated by using stratification analysis. Moreover, we constructed a nomogram based on risk score, age and stage, which had a strong ability to forecast the OS of the OC patients. The influence of risk score on immune infiltration and immunotherapy response were assessed and the results suggested that patients with high-risk score might recruit multiple immune cells and stromal cells, leading to facilitating immune surveillance evasive. Ultimately, we demonstrated that the risk model was associated with chemotherapy response of multiple antitumor drugs, especially for paclitaxel, metformin and veliparib, which are commonly used in treating OC patients. In conclusion, we constructed a novel immune-related lncRNA signature, which had a potential prognostic value for OC patients and might facilitate personalized counselling for immunotherapy and chemotherapy.

A review of bioinformatics tools and web servers in different microarray platforms used in cancer research

Over the past decade, conventional lab work strategies have gradually shifted from being limited to a laboratory setting towards a bioinformatics era to help manage and process the vast amounts of data generated by omics technologies. The present work outlines the latest contributions of bioinformatics in analyzing microarray data and their application to cancer. We dissect different microarray platforms and their use in gene expression in cancer models. We highlight how computational advances empowered the microarray technology in gene expression analysis. The study on protein-protein interaction databases classified into primary, derived, meta-database, and prediction databases describes the strategies to curate and predict novel interaction networks in silico. In addition, we summarize the areas of bioinformatics where neural graph networks are currently being used, such as protein functions, protein interaction prediction, and in silico drug discovery and development. We also discuss the role of deep learning as a potential tool in the prognosis, diagnosis, and treatment of cancer. Integrating these resources efficiently, practically, and ethically is likely to be the most challenging task for the healthcare industry over the next decade; however, we believe that it is achievable in the long term.

Integrative genomic and transcriptomic analysis reveals immune subtypes and prognostic markers in ovarian clear cell carcinoma

BACKGROUND

We performed an integrative genomic and transcriptomic profiling to identify molecular subtypes and prognostic markers with special focus on immune-related pathways.

METHODS

Totally, 50 Chinese patients were subjected to targeted next-generation sequencing and transcriptomic sequencing.

RESULTS

Two distinct subgroups were identified as immune (22.0%) and non-immune (78.0%) based on the immune-pathway related hierarchical clustering. Surprisingly, patients with immune subtype had a significantly worse survival. The prognostic capacity was validated in external cohorts. The immune group had higher expression of genes involved in pro-inflammation and checkpoints. PD-1 signalling pathway was enriched in the immune subtype. Besides, the immune cluster presented enriched expression of genes involved in epithelial-mesenchymal transition, angiogenesis and PI3K-AKT-mTOR signalling, while the non-immune subtype had higher expression of metabolic pathways. The immune subtype had a higher mutation rate of PIK3CA though significance was not achieved. Lastly, we established a prognostic immune signature for overall survival. Interestingly, the immune signature could also be applied to renal clear cell carcinoma, but not to other histologic subtype of ovarian cancer.

CONCLUSIONS

An immune subtype of OCCC was identified with poor survival and enrichment of PD-1 and PI3K-AKT-mTOR signalling. We constructed and validated a robust prognostic immune signature of OCCC patients.

Glycosylation-Related Genes Predict the Prognosis and Immune Fraction of Ovarian Cancer Patients Based on Weighted Gene Coexpression Network Analysis (WGCNA) and Machine Learning

Background

Ovarian cancer (OC) is a malignancy exhibiting high mortality in female tumors. Glycosylation is a posttranslational modification of proteins but research has failed to demonstrate a systematic link between glycosylation-related signatures and tumor environment of OC.

Purpose

This study is aimed at developing a novel model with glycosylation-related messenger RNAs (GRmRNAs) to predict the prognosis and immune function in OC patients.

Methods

The transcriptional profiles and clinical phenotypes of OC patients were collected from the Gene Expression Omnibus and The Cancer Genome Atlas databases. A weighted gene coexpression network analysis and machine learning were performed to find the optimal survival-related GRmRNAs. Least absolute shrinkage and selection operator regression (LASSO) and Cox regression were carried out to calculate the coefficients of each GRmRNA and compute the risk score of each patient as well as develop a prognostic model. A nomogram model was constructed, and several algorithms were used to investigate the relationship between risk subtypes and immune-infiltrating levels.

Results

A total of four signatures (ALG8, DCTN4, DCTN6, and UBB) were determined to calculate the risk scores, classifying patients into the high-and low-risk groups. High-risk patients exhibited significantly poorer survival outcomes, and the established nomogram model had a promising prediction for OC patients' prognosis. Tumor purity and tumor mutation burden were negatively correlated with risk scores. In addition, risk scores held statistical associations with pathway signatures such as Wnt, Hippo, and reactive oxygen species, and nonsynonymous mutation counts.

Conclusion

The currently established risk scores based on GRmRNAs can accurately predict the prognosis, the immune microenvironment, and the immunotherapeutic efficacy of OC patients.

OSov: An Interactive Web Server to Evaluate Prognostic Biomarkers for Ovarian Cancer

Simple Summary

The OSov web server incorporates gene expression profiles with clinical risk factors to estimate the ovarian cancers patients’ survival, and provides a tool for multiple analysis, such as forest-plot, uni/multi-variate survival analysis, Kaplan-Meier plot and nomogram construction.

Abstract

Ovarian cancer is one of the most aggressive and highly lethal gynecological cancers. The purpose of our study is to build a free prognostic web server to help researchers discover potential prognostic biomarkers by integrating gene expression profiling data and clinical follow-up information of ovarian cancer. We construct a prognostic web server OSov (Online consensus Survival analysis for Ovarian cancer) based on RNA expression profiles. OSov is a user-friendly web server which could present a Kaplan–Meier plot, forest plot, nomogram and survival summary table of queried genes in each individual cohort to evaluate the prognostic potency of each queried gene. To assess the performance of OSov web server, 163 previously published prognostic biomarkers of ovarian cancer were tested and 72% of them had their prognostic values confirmed in OSov. It is a free and valuable prognostic web server to screen and assess survival-associated biomarkers for ovarian cancer.

Construction of a Macrophage Infiltration Regulatory Network and Related Prognostic Model of High-Grade Serous Ovarian Cancer

Background

High-grade serous ovarian cancer (HGSOC) carries the highest mortality in the gynecological cancers; however, therapeutic outcomes have not significantly improved in recent decades. Macrophages play an essential role in the occurrence and development of ovarian cancer, so the mechanisms of macrophage infiltration should be elucidated.

Method

We downloaded transcriptome data of ovarian cancers from the Gene Expression Omnibus and The Cancer Genome Atlas. After rigorous screening, 1566 HGSOC were used for data analysis. CIBERSORT was used to estimate the level of macrophage infiltration and WGCNA was used to identify macrophage-related modules. We constructed a macrophage-related prognostic model using machine learning LASSO algorithm and verified it using multiple HGSOC cohorts.

Results

In the GPL570-OV cohort, high infiltration level of M1 macrophages was associated with a good outcome, while high infiltration level of M2 macrophages was associated with poor outcomes. We used WGCNA to select genes correlated with macrophage infiltration. These genes were used to construct protein-protein interaction maps of macrophage infiltration. IFL44L, RSAD2, IFIT3, MX1, IFIH1, IFI44, and ISG15 were the hub genes in the network. We then constructed a macrophage-related prognostic model composed of CD38, ACE2, BATF2, HLA-DOB, and WARS. The model had the ability to predict the overall survival rate of HGSOC patients in GPL570-OV, GPL6480-OV, TCGA-OV, GSE50088, and GSE26712. In exploring the immune microenvironment, we found that CD4 memory T cells and activated mast cells showed that the degree of infiltration was higher in the high-risk group, while M1 macrophages were the opposite, and HLA molecules were overexpressed in the high-risk group.

Conclusion

We constructed a macrophage infiltration-related protein interaction network that provides a basis for studying macrophages in HGSOC. Our macrophage-related prognostic model is robust and widely applicable. It predicts overall survival in HGSOC patients and may improve HGSOC treatment.

Identification of a glycolysis-related gene signature for survival prediction of ovarian cancer patients

Background

Ovarian cancer (OV) is deemed the most lethal gynecological cancer in women. The aim of this study was to construct an effective gene prognostic model for predicting overall survival (OS) in patients with OV.

Methods

The expression profiles of glycolysis‐related genes (GRGs) and clinical data of patients with OV were extracted from The Cancer Genome Atlas (TCGA) database. Univariate, multivariate, and least absolute shrinkage and selection operator Cox regression analyses were conducted, and a prognostic signature based on GRGs was constructed. The predictive ability of the signature was analyzed using training and test sets.

Results

A gene risk signature based on nine GRGs (ISG20, CITED2, PYGB, IRS2, ANGPTL4, TGFBI, LHX9, PC, and DDIT4) was identified to predict the survival outcome of patients with OV. The signature showed a good prognostic ability for OV, particularly high‐grade OV, in the TCGA dataset, with areas under the curve (AUC) of 0.709 and 0.762 for 3‐ and 5‐year survival, respectively. Similar results were found in the test sets, and the AUCs of 3‐, 5‐year OS were 0.714 and 0.772 in the combined test set. And our signature was an independent prognostic factor. Moreover, a nomogram combining the prediction model and clinical factors was developed.

Conclusion

Our study established a nine‐GRG risk model and nomogram to better predict OS in patients with OV. The risk model represents a promising and independent prognostic predictor for patients with OV. Moreover, our study on GRGs could offer guidance for the elucidation of underlying mechanisms in future studies.

Landscape of Immune Microenvironment in Epithelial Ovarian Cancer and Establishing Risk Model by Machine Learning

Background

Epithelial ovarian cancer (EOC) is an extremely lethal gynecological malignancy and has the potential to benefit from the immune checkpoint blockade (ICB) therapy, whose efficacy highly depends on the complex tumor microenvironment (TME).

Method and Result

We comprehensively analyze the landscape of TME and its prognostic value through immune infiltration analysis, somatic mutation analysis, and survival analysis. The results showed that high infiltration of immune cells predicts favorable clinical outcomes in EOC. Then, the detailed TME landscape of the EOC had been investigated through “xCell” algorithm, Gene set variation analysis (GSVA), cytokines expression analysis, and correlation analysis. It is observed that EOC patients with high infiltrating immune cells have an antitumor phenotype and are highly correlated with immune checkpoints. We further found that dendritic cells (DCs) may play a dominant role in promoting the infiltration of immune cells into TME and forming an antitumor immune phenotype. Finally, we conducted machine-learning Lasso regression, support vector machines (SVMs), and random forest, identifying six DC-related prognostic genes (CXCL9, VSIG4, ALOX5AP, TGFBI, UBD, and CXCL11). And DC-related risk stratify model had been well established and validated.

Conclusion

High infiltration of immune cells predicted a better outcome and an antitumor phenotype in EOC, and the DCs might play a dominant role in the initiation of antitumor immune cells. The well-established risk model can be used for prognostic prediction in EOC.

Common miRNAs, candidate genes and their interaction network across four subtypes of epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is categorized into four major histological subtypes such as clear cell carcinoma (CCC), endometrioid carcinoma (EC), mucinous carcinoma (MC), and serous carcinoma (SC). Heterogeneity of the EOC leads to different clinical outcomes of the disease, although all the subtypes are originated from the same layer of tissue. Therefore, it is of interest to identify the common candidate genes, miRNA and their interaction network in four the subtypes of EOC. A comparative gene expression analysis identified 248 common differentially expressed genes (DEGs) in the four subtypes of EOC. Identified common DEGs were found to be enriched in cancer specific pathways. A protein-protein interaction (PPI) network of the common DEGs were constructed, and subsequent module and survival analyses identified seven key candidate genes (CCNB1, CENPM, CEP55, RACGAP1, TPX2, UBE2C, and ZWINT). We also documented 10 key candidate miRNAs (hsa-mir-16-5p, hsa-mir-23b-3p, hsa-mir-34a-5p, hsa-mir-103a-3p, hsa-mir-107, hsa-mir-124-3p, hsa-mir-129-2-3p, hsa-mir-147a, hsa-mir-205-5p, and hsa-mir-195-5p) linked to the candidate genes. These derived data find application in the understanding of EOC.

Differential expression analysis in ovarian cancer: A functional genomics and systems biology approach

BACKGROUND

Ovarian cancer is one of the rarest lethal oncologic diseases that have hardly any specific biomarkers. The availability of high-throughput genomic data and advancement in bioinformatics tools allow us to predict gene biomarkers and apply systems biology approaches to get better diagnosis, and prognosis of the disease with a tentative drug that may be repurposed.

OBJECTIVE

To perform genome-wide association studies using microarray gene expression of ovarian cancer and identify gene biomarkers, construction and analyze networks, perform survival analysis, and drug interaction studies for better diagnosis, prognosis, and treatment of ovarian cancer.

METHOD

The gene expression profiles of both healthy and serous ovarian cancer epithelial samples were considered. We applied a series of bioinformatics methods and tools, including fold-change statistics for differential expression analysis, DisGeNET and NCBI-Gene databases for gene-disease association mapping, DAVID 6.8 for GO enrichment analysis, GeneMANIA for network construction, Cytoscape 3.8 with its plugins for network visualization, analysis, and module detection, the UALCAN for patient survival analysis, and PubChem, DrugBank and DGIdb for gene-drug interaction.

RESULTS

We identified 8 seed genes that were subjected for drug-gene interaction studies. Because of over-expression in all the four stages of ovarian cancer, we discern that genes HMGA1 and PSAT1 are potential therapeutic biomarkers for its diagnosis at an early stage (stage I). Our analysis suggests that there are 11 drugs common in the seed genes. However, hypermethylated seed genes HMGA1 and PSAT1 showcased a good interaction affinity with drugs cisplatin, cyclosporin, bisphenol A, progesterone, and sunitinib, and are crucial in the proliferation of ovarian cancer.

CONCLUSION

Our study reveals that HMGA1 and PSAT1 can be deployed for initial screening of ovarian cancer and drugs cisplatin, bisphenol A, cyclosporin, progesterone, and sunitinib are effective in curbing the epigenetic alteration.

A risk prediction model of gene signatures in ovarian cancer through bagging of GA-XGBoost models

Introduction

Ovarian cancer (OC) is one of the most frequent gynecologic cancers among women, and high-accuracy risk prediction techniques are essential to effectively select the best intervention strategies and clinical management for OC patients at different risk levels. Current risk prediction models used in OC have low sensitivity, and few of them are able to identify OC patients at high risk of mortality, which would both optimize the treatment of high-risk patients and prevent unnecessary medical intervention in those at low risk.

Objectives

To this end, we have developed a bagging-based algorithm with GA-XGBoost models that predicts the risk of death from OC using gene expression profiles.

Methods

Four gene expression datasets from public sources were used as training (n = 1) or validation (n = 3) sets. The performance of our proposed algorithm was compared with fine-tuning and other existing methods. Moreover, the biological function of selected genetic features was further interpreted, and the response to a panel of approved drugs was predicted for different risk levels.

Results

The proposed algorithm showed good sensitivity (74-100%) in the validation sets, compared with two simple models whose sensitivity only reached 47% and 60%. The prognostic gene signature used in this study was highly connected to AKT, a key component of the PI3K/AKT/mTOR signaling pathway, which influences the tumorigenesis, proliferation, and progression of OC.

Conclusion

These findings demonstrated an improvement in the sensitivity of risk classification of OC patients with our risk prediction models compared with other methods. Ongoing effort is needed to validate the outcomes of this approach for precise clinical treatment.

Identification of the miRNA signature associated with survival in patients with ovarian cancer

Ovarian cancer is a major gynaecological malignant tumor associated with a high mortality rate. Identifying survival-related variants may improve treatment and survival in patients with ovarian cancer. In this work, we proposed a support vector regression (SVR)-based method called OV-SURV, which is incorporated with an inheritable bi-objective combinatorial genetic algorithm for feature selection to identify a miRNA signature associated with survival in patients with ovarian cancer. There were 209 patients with miRNA expression profiles and survival information of ovarian cancer retrieved from The Cancer Genome Atlas database. OV-SURV achieved a mean correlation coefficient of 0.77±0.01and a mean absolute error of 0.69±0.02 years using 10-fold cross-validation. Analysis of the top ranked miRNAs revealed that the miRNAs, hsa-let-7f, hsa-miR-1237, hsa-miR-98, hsa-miR-933, and hsa-miR-889, were significantly associated with the survival in patients with ovarian cancer. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that four of these miRNAs, hsa-miR-182, hsa-miR-34a, hsa-miR-342, and hsa-miR-1304, were highly enriched in fatty acid biosynthesis, and the five miRNAs, hsa-let-7f, hsa-miR-34a, hsa-miR-342, hsa-miR-1304, and hsa-miR-24, were highly enriched in fatty acid metabolism. The prediction model with the identified miRNA signature consisting of prognostic biomarkers can benefit therapeutic decision making of ovarian cancer.

Synthetic lethality-mediated precision oncology via the tumor transcriptome

Precision oncology has made significant advances, mainly by targeting actionable mutations in cancer driver genes. Aiming to expand treatment opportunities, recent studies have begun to explore the utility of tumor transcriptome to guide patient treatment. Here, we introduce SELECT (synthetic lethality and rescue-mediated precision oncology via the transcriptome), a precision oncology framework harnessing genetic interactions to predict patient response to cancer therapy from the tumor transcriptome. SELECT is tested on a broad collection of 35 published targeted and immunotherapy clinical trials from 10 different cancer types. It is predictive of patients' response in 80% of these clinical trials and in the recent multi-arm WINTHER trial. The predictive signatures and the code are made publicly available for academic use, laying a basis for future prospective clinical studies.

Development and Validation of an Immune-Related Gene-Pair Model of High-Grade Serous Ovarian Cancer After Platinum-Based Chemotherapy

Background

High-grade serous ovarian cancer (HGSOC) is a common cause of death from gynecological cancer, with an overall survival rate that has not significantly improved in decades. Reliable bio-markers are needed to identify high-risk HGSOC to assist in the selection and development of treatment options.

Method

The study included ten HGSOC cohorts, which were merged into four separate cohorts including a total of 1,526 samples. We used the relative expression of immune genes to construct the gene-pair matrix, and the least absolute shrinkage and selection operator regression was performed to build the prognosis model using the training set. The prognosis of the model was verified in the training set (363 cases) and three validation sets (of 251, 354, and 558 cases). Finally, the differences in immune cell infiltration and gene enrichment pathways between high and low score groups were identified.

Results

A prognosis model of HGSOC overall survival rate was constructed in the training set, and included data for 35 immune gene-related gene pairs and the regression coefficients. The risk stratification of HGSOC patients was successfully performed using the training set, with a p-value of Kaplan-Meier of < 0.001. A score from this model is an independent prognostic factor of HGSOC, and prognosis was evaluated in different clinical subgroups. This model was also successful for the other three validation sets, and the results of Kaplan-Meier analysis were statistically significant (p < 0.05). The model can also predict patient progression-free survival with HGSOC to reflect tumor growth status. There was a lower infiltration level of M1 macrophages in the high-risk group compared to that in the low-risk group (p < 0.001). Finally, the immune-related pathways were enriched in the low-risk group.

Conclusion

The prognostic model based on immune-related gene pairs developed is a potential prognostic marker for high-grade serous ovarian cancer treated with platinum. The model has robust prognostic ability and wide applicability. More prospective studies will be needed to assess the practical application of this model for precision therapy.

Suppression of ARID1A associated with decreased CD8 T cells improves cell survival of ovarian clear cell carcinoma

Objective

AT-rich interactive domain 1A (ARID1A) plays an important role as a tumor suppressor gene in ovarian clear cell carcinoma (OCCC), but the clinical application of ARID1A remains unclear. The aim of this study was to analyze clinicopathological parameters, molecular interactions and immune-infiltration in patients with low ARID1A expression and to provide candidate target drugs.

Methods

We investigated the clinicopathologic parameters, specific gene sets/genes, and immunological relevance according to ARID1A expression in 998 OCCC patients from 12 eligible studies (using meta-analyses); 30 OCCC patients from the Hanyang University Guri Hospital (HYGH) cohort; and 52 OCCC patients from gene set enrichment (GSE) 65986 (25 patients), 63885 (9 patients), and 54809 (6 patients and 12 healthy people) of the Gene Expression Omnibus (GEO). We analyzed network-based pathways based on gene set enrichment analysis (GSEA) and performed in vitro drug screening.

Results

Low ARID1A expression was associated with poor survival in OCCC from the meta-analysis, HYGH cohort and GEO data. In GSEA, low ARID1A expression was related to the tumor invasion process as well as a low immune-infiltration. In silico cytometry showed that CD8 T cells were decreased with low ARID1A expression. In pathway analysis, ARID1A was associated with angiogenic endothelial cell signaling. In vitro drug screening revealed that cabozantinib and bicalutamide effectively inhibited specific hub genes, such as vascular endothelial growth factor-A and androgen receptor, in OCCC cells with low ARID1A expression.

Conclusions

Therapeutic strategies making use of low ARID1A could contribute to better clinical management/research for patients with OCCC.

Development and Validation of the Gene Expression Predictor of High-grade Serous Ovarian Carcinoma Molecular SubTYPE (PrOTYPE)

PURPOSE

Gene expression-based molecular subtypes of high-grade serous tubo-ovarian cancer (HGSOC), demonstrated across multiple studies, may provide improved stratification for molecularly targeted trials. However, evaluation of clinical utility has been hindered by nonstandardized methods, which are not applicable in a clinical setting. We sought to generate a clinical grade minimal gene set assay for classification of individual tumor specimens into HGSOC subtypes and confirm previously reported subtype-associated features.

EXPERIMENTAL DESIGN

Adopting two independent approaches, we derived and internally validated algorithms for subtype prediction using published gene expression data from 1,650 tumors. We applied resulting models to NanoString data on 3,829 HGSOCs from the Ovarian Tumor Tissue Analysis consortium. We further developed, confirmed, and validated a reduced, minimal gene set predictor, with methods suitable for a single-patient setting.

RESULTS

Gene expression data were used to derive the predictor of high-grade serous ovarian carcinoma molecular subtype (PrOTYPE) assay. We established a de facto standard as a consensus of two parallel approaches. PrOTYPE subtypes are significantly associated with age, stage, residual disease, tumor-infiltrating lymphocytes, and outcome. The locked-down clinical grade PrOTYPE test includes a model with 55 genes that predicted gene expression subtype with >95% accuracy that was maintained in all analytic and biological validations.

CONCLUSIONS

We validated the PrOTYPE assay following the Institute of Medicine guidelines for the development of omics-based tests. This fully defined and locked-down clinical grade assay will enable trial design with molecular subtype stratification and allow for objective assessment of the predictive value of HGSOC molecular subtypes in precision medicine applications.See related commentary by McMullen et al., p. 5271.

Four genes relevant to pathological grade and prognosis in ovarian cancer

BACKGROUND: Ovarian cancer is the common tumor in female, the prognostic of which is influenced by a series of factors. In this study, 4 genes relevant to pathological grade in ovarian cancer were screened out by the construction of weighted gene co-expression network analysis. METHODS: GSE9891 with 298 ovarian cancer cases had been used to construct co-expression networks. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses was used to analyze the possible mechanism of genes involved in the malignant process of ovarian cancer. Hub genes were validated in other independent datasets, such as GSE63885, GSE26193 and GSE30161. Survival analysis based on the hub genes was performed by website of Kaplan Meier-plotter. RESULTS: The result based on weighted gene co-expression network analysis indicated that turquoise module has the highest association with pathological grade. Gene Ontology enrichment analysis revealed that the genes in turquoise module main enrichment in inflammatory response and immune response. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the genes in turquoise module main enrichment in cytokine-cytokine receptor interaction and chemokine signaling pathway. In turquoise module, a total of 4 hub genes (MS4A4A, CD163, CPR65, MS4A6A) were identified. Then, 4 hub genes were effectively verified in the test datasets (GSE63885, GSE26193 and GSE30161) and tissue samples from Shengjing Hospital of China Medical University. Survival analysis indicated that the 4 hub genes were associated with poor progression-free survival of ovarian cancer. CONCLUSIONS: In conclusion, 4 hub genes (MS4A4A, CD163, CPR65, MS4A6A) were verified associated with pathological grade of ovarian cancer. Moreover, MS4A4A, CD163, MS4A6A may serve as a surface marker for M2 macrophages. Targeting the 4 hub genes may can improve the prognosis of ovarian cancer.

Evaluation of the Role of ITGBL1 in Ovarian Cancer

In our previous microarray study we identified two subgroups of high-grade serous ovarian cancers with distinct gene expression and survival. Among differentially expressed genes was an Integrin beta-like 1 (ITGBL1), coding for a poorly characterized protein comprised of ten EGF-like repeats. Here, we have analyzed the influence of ITGBL1 on the phenotype of ovarian cancer (OC) cells. We analyzed expression of four putative ITGBL1 mRNA isoforms in five OC cell lines. OAW42 and SKOV3, having the lowest level of any ITGBL1 mRNA, were chosen to produce ITGBL1-overexpressing variants. In these cells, abundant ITGBL1 mRNA expression could be detected by RT-PCR. Immunodetection was successful only in the culture media, suggesting that ITGBL1 is efficiently secreted. We found that ITGBL1 overexpression affected cellular adhesion, migration and invasiveness, while it had no effect on proliferation rate and the cell cycle. ITGBL1-overexpressing cells were significantly more resistant to cisplatin and paclitaxel, major drugs used in OC treatment. Global gene expression analysis revealed that signaling pathways affected by ITGBL1 overexpression were mostly those related to extracellular matrix organization and function, integrin signaling, focal adhesion, cellular communication and motility; these results were consistent with the findings of our functional studies. Overall, our results indicate that higher expression of ITGBL1 in OC is associated with features that may worsen clinical course of the disease.

Differences in Gene Expression Profile of Primary Tumors in Metastatic and Non-Metastatic Papillary Thyroid Carcinoma-Do They Exist?

Molecular mechanisms of distant metastases (M1) in papillary thyroid cancer (PTC) are poorly understood. We attempted to analyze the gene expression profile in PTC primary tumors to seek the genes associated with M1 status and characterize their molecular function. One hundred and twenty-three patients, including 36 M1 cases, were subjected to transcriptome oligonucleotide microarray analyses: (set A-U133, set B-HG 1.0 ST) at transcript and gene group level (limma, gene set enrichment analysis (GSEA)). An additional independent set of 63 PTCs, including 9 M1 cases, was used to validate results by qPCR. The analysis on dataset A detected eleven transcripts showing significant differences in expression between metastatic and non-metastatic PTC. These genes were validated on microarray dataset B. The differential expression was positively confirmed for only two genes: IGFBP3, (most significant) and ECM1. However, when analyzed on an independent dataset by qPCR, the IGFBP3 gene showed no differences in expression. Gene group analysis showed differences mainly among immune-related transcripts, indicating the potential influence of tumor immune infiltration or signal within the primary tumor. The differences in gene expression profile between metastatic and non-metastatic PTC, if they exist, are subtle and potentially detectable only in large datasets.

IQGAP3 promotes cancer proliferation and metastasis in high-grade serous ovarian cancer

Ovarian cancer is a type of gynecological cancer with the highest mortality rate worldwide. Due to a lack of effective screening methods, most cases are diagnosed at later stages where the survival rates are poor. Thus, it is termed a ‘silent killer’ and is the most lethal of all the malignancies in women. IQ motif containing GTPase Activating Protein 3 (IQGAP3) is a member of the Rho family of GTPases, and plays a crucial role in the development and progression of several types of cancer. The aim of the present study was to investigate the oncogenic functions and mechanisms of IQGAP3 on the proliferation and metastasis of high-grade serous ovarian cancer (HGSOC). Therefore, the expression levels of IQGAP3 in HGSOC and normal tissue samples were compared, and IQGAP3 knockdown was performed to examine its functional role using various in vitro and in vivo experiments. It was demonstrated that the expression of IQGAP3 was upregulated in HGSOC tissues compared with the healthy tissues; this differential expression was also observed in the ovarian cancer cell lines. Functional experimental results suggested that IQGAP3 silencing significantly reduced proliferation, migration and invasion in ovarian cancer cell lines. Moreover, in vivo experimental findings validated the in vitro results, where the tumorigenic and metastatic capacities of IQGAP3-silenced cells were significantly lower in the nude mice compared with the mice implanted with the control cells. Furthermore, knockdown of IQGAP3 resulted in increased apoptosis, and the effects of IQGAP3 expression on various epithelial-mesenchymal transition markers were identified, suggesting a possible mechanism associated with the role of IQGAP3 in metastasis. The effect of IQGAP3 silencing on chemosensitivity towards olaparib was also assessed. Collectively, the present results indicated that IQGAP3 is a potential diagnostic and prognostic marker, and a putative therapeutic target of HGSOC.

Fibronectin and Periostin as Prognostic Markers in Ovarian Cancer

Previously, based on a DNA microarray experiment, we identified a 96-gene prognostic signature associated with the shorter survival of ovarian cancer patients. We hypothesized that some differentially expressed protein-coding genes from this signature could potentially serve as prognostic markers. The present study was aimed to validate two proteins, namely fibronectin (FN1) and periostin (POSTN), in the independent set of ovarian cancer samples. Both proteins are mainly known as extracellular matrix proteins with many important functions in physiology. However, there are also indications that they are implicated in cancer, including ovarian cancer. The expression of these proteins was immunohistochemically analyzed in 108 surgical samples of advanced ovarian cancer (majority: high-grade serous) and additionally on tissue arrays representing different stages of the progression of ovarian and fallopian tube epithelial tumors, from normal epithelia, through benign tumors, to adenocarcinomas of different stages. The correlation with clinical, pathological, and molecular features was evaluated. Kaplan-Meier survival analysis and Cox-proportional hazards models were used to estimate the correlation of the expression levels these proteins with survival. We observed that the higher expression of fibronectin in the tumor stroma was highly associated with shorter overall survival (OS) (Kaplan-Meier analysis, log-rank test p = 0.003). Periostin was also associated with shorter OS (p = 0.04). When we analyzed the combined score, calculated by adding together individual scores for stromal fibronectin and periostin expression, Cox regression demonstrated that this joint FN1&POSTN score was an independent prognostic factor for OS (HR = 2.16; 95% CI: 1.02-4.60; p = 0.044). The expression of fibronectin and periostin was also associated with the source of ovarian tumor sample: metastases showed higher expression of these proteins than primary tumor samples (χ² test, p = 0.024 and p = 0.032). Elevated expression of fibronectin and periostin was also more common in fallopian cancers than in ovarian cancers. Our results support some previous observations that fibronectin and periostin have a prognostic significance in ovarian cancer. In addition, we propose the joint FN1&POSTN score as an independent prognostic factor for OS. Based on our results, it may also be speculated that these proteins are related to tumor progression and/or may indicate fallopian-epithelial origin of the tumor.

Analysis of gene expression signatures identifies prognostic and functionally distinct ovarian clear cell carcinoma subtypes

Background

Ovarian clear cell carcinoma (OCCC) is a histological subtype of epithelial ovarian cancer (EOC) with distinct pathological, biological, and molecular features. OCCCs are more resistant to conventional treatment regimen of EOC and have the worst stage-adjusted prognosis amongst EOC subtypes. As the OCCC incidence rate in Asian populations has significantly increased in recent decades, it is critical to elucidate its molecular features that could lead to OCCC-tailored therapeutic strategies.

Methods

Gene expression profiles of 222 OCCC were analyzed by hierarchical clustering and statistical analyses.

Findings

We identified two OCCC gene expression subtypes: EpiCC—epithelial-like, which is associated with early-stage disease, with a relatively higher rate of gene mutations in the SWI/SNF complex; and MesCC—mesenchymal-like, associated with late-stage and higher enrichment of immune-related pathway activity. Genetic, copy number and transcriptomic analyses showed that both EpiCC and MesCC carried OCCC-associated aberrations. The EpiCC/MesCC classification was reproducible in validation cohorts and OCCC cell lines. MesCC tumors had a poorer progression-free survival (PFS) than EpiCC tumors (HR: 3·0, p = 0·0006). Functional assays in cell lines showed that the MesCC subtype was more proliferative and more anoikis-resistant than the EpiCC. By applying the EpiCC/MesCC classification to the TCGA renal clear cell carcinoma cohort, our results indicated interoperability of the subtyping scheme, and revealed preferential drug response of MesCC to bevacizumab.

Interpretation

The EpiCC/MesCC classification shows promise for prognostic and therapeutic stratification in OCCC patients and warrants further investigation in the context of OCCC gene expression subtype-tailored treatment strategies.

ATAD2 predicts poor outcomes in patients with ovarian cancer and is a marker of proliferation

The oncogene ATPase family AAA domain-containing protein 2 (ATAD2) has been demonstrated to promote malignancy in a number of different types of tumor; however, its expression and role in ovarian cancer (OC) remain unknown. In the present study, it was demonstrated that ATAD2 acts as both a marker and a driver of cell proliferation in OC. Immunohistochemistry (IHC) and bioinformatics analyses were used to evaluate ATAD2 expression in OC, and multi-omics integrated analyses were used to dissect which factor resulted in its upregulation. Multiplex IHC assay was used to reveal the specific expression of ATAD2 in proliferating OC cells. CRISPR-Cas9-mediated gene editing was performed to investigate the effect of ATAD2 deletion on OC proliferation. The results demonstrated that ATAD2 is elevated in primary OC tissues compared with the adjacent normal tissue and metastases from the stomach. Genetic copy number amplification is a primary cause resulting in upregulation of ATAD2, and this was most frequently observed in OC. High ATAD2 expression was associated with advanced progression and predicted an unfavorable prognosis. ATAD2 could be used to identify cases of OC with a high proliferation signature and could label proliferating cells in OC. CRISPR-Cas9-mediated ATAD2 deletion resulted in a significant decrease in both cell proliferation and colony formation ability. Mechanistically, ATAD2-knockdown resulted in deactivation of the mitogen-activated protein kinase (MAPK) pathways, particularly the JNK-MAPK pathway, resulting in suppression of proliferation. Collectively, the data from the present study demonstrated that the ATD2 gene was frequently amplified and protein expression levels were upregulated in OC. Therefore, ATAD2 may serve as an attractive diagnostic and prognostic OC marker, which may be used to identify patients with primary OC, whom are most likely to benefit from ATAD2 gene-targeted proliferation intervention therapies.

Integrated pan-cancer gene expression and drug sensitivity analysis reveals SLFN11 mRNA as a solid tumor biomarker predictive of sensitivity to DNA-damaging chemotherapy

Background

Precision oncology seeks to integrate multiple layers of data from a patient’s cancer to effectively tailor therapy. Conventional chemotherapies are sometimes effective but accompanied by adverse events, warranting the identification of a biomarker of chemosensitivity.

Objective

Identify an mRNA biomarker that predicts chemosensitivity across solid tumor subtypes.

Methods

We performed a pan-solid tumor analysis integrating gene expression and drug sensitivity profiles from 3 cancer cell line datasets to identify transcripts correlated with sensitivity to a panel of chemotherapeutics. We then tested the ability of an mRNA biomarker to predictive clinical outcomes in cohorts of patients with breast, lung, or ovarian cancer.

Results

Expression levels of several mRNA transcripts were significantly correlated with sensitivity or resistance chemotherapeutics in cancer cell line datasets. The only mRNA transcript significantly correlated with sensitization to multiple classes of DNA-damaging chemotherapeutics in all 3 cell line datasets was encoded by Schlafen Family Member 11 (SLFN11). Analyses of multiple breast, lung, and ovarian cancer patient cohorts treated with chemotherapy confirmed SLFN11 mRNA expression as a predictive biomarker of longer overall survival and improved tumor response.

Conclusions

Tumor SLFN11 mRNA expression is a biomarker of sensitivity to an array of DNA-damaging chemotherapeutics across solid tumor subtypes.

Increased expression of TET3 predicts unfavorable prognosis in patients with ovarian cancer-a bioinformatics integrative analysis

Ovarian carcinoma is a lethal gynecological malignancy. Women with ovarian cancer (OC) are highly recurrent and typically diagnosed at late stage. Ten-eleven translocation protein 3 (TET3) belongs to the family of ten-eleven translocations (TETs) which induce DNA demethylation and gene regulation in epigenetic level by converting 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC). Previous studies indicated that TET3 is overexpressed in ovarian cancer tissues. However, the clinic-pathological functions and prognostic values of TET3 remain unclear. Here we performed an integrative study to identify the role of TET3 by bioinformatics analysis. The TET3 expression in ovarian cancer was assessed with Oncomine database, and validated with TCGA and GTEx database. The correlation of TET3 gene alteration and clinic-pathological functions was addressed by integrative analysis of GEO datasets. Then we showed mainly TET3 gain and diploid but less deletion in ovarian cancer by copy number alteration (CNA) or mutation analysis with cBioPortal. Furthermore, by using Kaplan-Meier plotter (K-M plotter), we evaluated that high TET3 level was associated with poor survival in ovarian cancer patients, which was validated with analysis by PrognoScan database and gene differential analyses with TCGA and GTEx. This is the first study demonstrated that elevated expression of TET3 is associated with poor clinic-pathological functions, poor prognosis, wherein TET3, which presents epigenetic changes or methylation changes, might be served as a diagnostic marker or therapeutic target for ovarian cancer.

Novel Model for Comprehensive Assessment of Robust Prognostic Gene Signature in Ovarian Cancer Across Different Independent Datasets

Different analytical methods or models can often find completely different prognostic biomarkers for the same cancer. In the study of prognostic molecular biomarkers of ovarian cancer (OvCa), different studies have reported a variety of prognostic gene signatures. In the current study, based on geometric concepts, the linearity-clustering phase diagram with integrated P-value (LCP) method was used to comprehensively consider three indicators that are commonly employed to estimate the quality of a prognostic gene signature model. The three indicators, namely, concordance index, area under the curve, and level of the hazard ratio were determined via calculation of the prognostic index of various gene signatures from different datasets. As evaluation objects, we selected 13 gene signature models (Cox regression model) and 16 OvCa genomic datasets (including gene expression information and follow-up data) from published studies. The results of LCP showed that three models were universal and better than other models. In addition, combining the three models into one model showed the best performance in all datasets by LCP calculation. The combination gene signature model provides a more reliable model and could be validated in various datasets of OvCa. Thus, our method and findings can provide more accurate prognostic biomarkers and effective reference for the precise clinical treatment of OvCa.

The core genes involved in the promotion of depression in patients with ovarian cancer

The present study aimed to identify the core genes and pathways involved in depression in patients with ovarian cancer (OC) who suffer from high or low-grade depression. The dataset GSE9116 from Gene Expression Omnibus database was analyzed to identify differentially expressed genes (DEGs) in these patients. To elucidate how certain genes could promote depression in patients with OC, pathway crosstalk, protein-protein interaction (PPI) and comprehensive gene-pathway analyses were determined using WebGestalt, ToppGene and Search Tool for the Retrieval of Interacting Genes and gene ontology analysis. Key genes and pathways were extracted from the gene-pathway network, and gene expression and survival analysis were evaluated. A total of 93 DEGs were identified from GSE9116 dataset, including 84 upregulated genes and nine downregulated genes. The PPI, pathway crosstalk and comprehensive gene-pathway analyses highlighted C-C motif chemokine ligand 2 (CCL2), Fos proto-oncogene, AP-1 transcription factor subunit (FOS), serpin family E member 1 (SERPINE1) and serpin family G member 1 (SERPING1) as core genes involved in the promotion of depression in patients with OC. These core genes were involved in the following four pathways 'Ensemble of genes encoding ECM-associated proteins including ECM-affiliated proteins', 'ECM regulators and secreted factors', 'Ensemble of genes encoding extracellular matrix and extracellular matrix-associated proteins' and 'MAPK signaling pathway and IL-17 signaling pathway'. The results from gene expression and survival analysis demonstrated that these four key genes were upregulated in patients with OC and high-grade depression and could worsen patients' survival. These results suggested that CCL2, FOS, SERPINE1 and SERPING1 may serve a crucial role in the promotion of depression in patients with OC. This finding may provide novel markers for predicting and treating depression in patients with OC; however, the underlying mechanisms remain unknown and require further investigation.

The epigenetic factor BORIS (CTCFL) controls the androgen receptor regulatory network in ovarian cancer

The identification of prognostic biomarkers is a priority for patients suffering from high-grade serous ovarian cancer (SOC), which accounts for >70% of ovarian cancer (OC) deaths. Meanwhile, borderline ovarian cancer (BOC) is a low malignancy tumor and usually patients undergo surgery with low probabilities of recurrence. However, SOC remains the most lethal neoplasm due to the lack of biomarkers for early diagnosis and prognosis. In this regard, BORIS (CTCFL), a CTCF paralog, is a promising cancer biomarker that is overexpressed and controls transcription in several cancer types, mainly in OC. Studies suggest that BORIS has an important function in OC by altering gene expression, but the effect and extent to which BORIS influences transcription in OC from a genome-wide perspective is unclear. Here, we sought to identify BORIS target genes in an OC cell line (OVCAR3) with potential biomarker use in OC tumor samples. To achieve this, we performed in vitro knockout and knockdown experiments of BORIS in OVCAR3 cell line followed by expression microarrays and bioinformatics network enrichment analysis to identify relevant BORIS target genes. In addition, ex vivo expression data analysis of 373 ovarian cancer patients were evaluated to identify the expression patterns of BORIS target genes. In vitro, we uncovered 130 differentially expressed genes and obtained the BORIS-associated regulatory network, in which the androgen receptor (AR) acts as a major transcription factor. Also, FN1, FAM129A, and CD97 genes, which are related to chemoresistance and metastases in OC, were identified. In SOC patients, we observed that malignancy is associated with high levels of BORIS expression while BOC patients show lower levels. Our study suggests that BORIS acts as a main regulator, and has the potential to be used as a prognostic biomarker and to yield novel drug targets among the genes BORIS controls in SOC patients.

Genome-wide bioinformatics analysis reveals CTCFL is upregulated in high-grade epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is the most lethal gynecological malignancy that threatens the health of females. Previous studies have demonstrated that the survival outcomes of patients with different EOC grades varied. Therefore, the EOC grade is considered to serve as a distinctive prognostic factor. To date, the evaluation of ovarian cancer grade relies on pathological examination and a quantitative index for diagnosis is lacking. Furthermore, the dysregulation of genes has been demonstrated to exert pivotal functions in the carcinogenesis of EOCs. Therefore, the identification of effective biomarkers associated with EOC grade is of importance for the development of therapeutic regimens, and also contributes to the prediction of EOC prognosis. Microarrays have been increasingly applied for the identification of potential molecular biomarkers for numerous diseases including EOC. In the present study, four public microarray datasets (GSE26193, GSE63885, GSE30161 and GSE9891) were analyzed. A total of 6,103 upregulated probes corresponding to 5,766 genes, and 4,004 downregulated probes corresponding to 3,707 genes were identified in the GSE26193, GSE63885 and GSE30161 datasets. ALK and LTK ligand 2 was the most downregulated gene associated with the tumor grade, while CCCTC-binding factor like (CTCFL), EGF like domain multiple 6, radical S-adenosyl methionine domain containing 2 and SAM and HD domain containing deoxynucleoside triphosphate triphosphohydrolase 1 were the most upregulated genes associated with EOC grade. The GSE9891 dataset was added for further analysis. Only one probe (1552368_at) encoding for CTCFL was identified to be consistently upregulated in the four examined datasets. Immunohistochemical analysis was used to detect the expression of CTCFL between low- and high-grade EOC tissues and revealed that the EOC grade was closely associated with CTCFL level. This was corroborated via the reverse transcription-quantitative polymerase chain reaction. Taken together, the results of the present study suggested that CTCFL is upregulated in high-grade epithelial ovarian cancer.

Long non-coding RNA LINC01627 is a prognostic risk factor for epithelial ovarian cancer

Ovarian malignancies are commonly diagnosed cancers of the female reproductive system. Recent studies have revealed that long non-coding RNAs (lncRNAs) can regulate a variety of oncological processes. In the present study, ovarian cancer expression datasets were searched in the GEO database using the GPL570 platform. Differential lncRNA expression between normal ovarian tissues and ovarian tumors were analyzed using the R package 'limma', and patient prognosis was accessed using the package 'survival'. Four databases, GSE14001, GSE18520, GSE38666 and GSE40595, were used for the analysis. A total of 64 lncRNAs were highly expressed and 4 were downregulated within these four databases. Prognostic analysis of the 68 lncRNAs in the four databases was performed, and revealed that the expression of long intergenic non-protein coding RNA 1627 (LINC01627) was negatively associated with patient prognosis in GSE19829 and GSE62193; there was no association between LINC01627 expression and patient's prognosis in GSE18520 or GSE63885. To investigate the proposed association between LINC01627 and patient prognosis, meta-analysis revealed that the total hazard ratio was 1.38 and the 95% confidence interval was between 1.04 and 1.83. Subgroup analysis revealed that LINC01627 may predict patient prognosis in high-grade, advanced and serous epithelial ovarian cancer, which was a risk factor for prognosis. Further assessment was performed in clinical samples and ovarian cancer cells, where the knockdown of LINC01627 inhibited the proliferative and migratory capacities of HO8910 and HEY cells. Collectively, the present results suggested that lncRNA LINC01627 may serve an oncogenic role in the development of epithelial ovarian tumors.

Gene Expression Indicates Altered Immune Modulation and Signaling Pathway Activation in Ovarian Cancer Patients Resistant to Topotecan

Epithelial ovarian cancer (EOC) is one of the deadliest gynecological malignancies. Topotecan remains an essential tool in second-line therapy; even so, most patients develop resistance within a short period of time. We aimed to identify biomarkers of topotecan resistance by using gene expression signatures derived from patient specimens at surgery and available subsequent responses to therapy. Gene expression was collected for 1436 patients and 10,103 genes. Based on disease progression, patients were categorized as responders/nonresponders depending on their progression free survival (PFS) state at 9, 12, 15 and 18 months after surgery. For each gene, the median expression was compared between responders and nonresponders for two treatment regimens (chemotherapy including/excluding topotecan) with Mann-Whitney U test at each of the four different PFS cutoffs. Statistical significance was accepted in the case of p < 0.05 with a fold change (FC) ≥ 1.44. Four genes (EPB41L2, HLA-DQB1, LTF and SFRP1) were consistently overexpressed across multiple PFS cutoff times in initial tumor samples of patients with disease progression following topotecan treatment. A common theme linked to topotecan resistance was altered immune modulation. Genes associated with disease progression after systemic chemotherapy emphasize the role of the initial organization of the tumor microenvironment in therapy resistance. Our results uncover biomarkers with potential utility for patient stratification.

Six-mRNA risk score system and nomogram constructed for patients with ovarian cancer

Platinum is a commonly used drug for the treatment of ovarian cancer (OC). The aim of the current study was to design and construct a risk score system for predicting the prognosis of patients with OC receiving platinum chemotherapy. The mRNA sequencing data and copy number variation (CNV) information (training set) of patients with OC were downloaded from The Cancer Genome Atlas database. A validation set, GSE63885, was obtained from Gene Expression Omnibus database. The differentially expressed genes (DEGs) and CNV genes (DECNs) between platinum-resistant and platinum-sensitive groups were identified using the limma package. The intersection between DEGs and DECNs were selected. Cox regression analysis was used to identify the genes and clinical factors associated with prognosis. Risk score system assessment and nomogram analysis were performed using the survival and rms packages in R. Gene Set Enrichment Analysis was used to identify the enriched pathways in high and low risk score groups. From 1,144 DEGs and 1,864 DECNs, 48 genes that occurred in the two datasets were selected. A total of six independent prognostic genes (T-box transcription factor T, synemin, tektin 5, growth differentiation factor 3, solute carrier family 22 member 3 and calcium voltage-gated channel subunit α1 C) and platinum response status were revealed to be associated with prognosis. Based on the six independent prognostic genes, a risk score system was constructed and assessed. Nomogram analysis revealed that the patients with the sensitive status and low risk scores had an improved prognosis. Furthermore, the current study revealed that the 574 DEGs identified were involved in eight pathways, including chemokine signaling pathway, toll-like receptor signaling pathway, cytokine-cytokine receptor interaction, RIG I like receptor signaling pathway, natural killer cell mediated cytotoxicity, apoptosis, T cell receptor signaling pathway and Fc ε receptor 1 signaling pathway. The six-mRNA risk score system designed in the present study may be used as prognosis predictor in patients with OC, whereas the nomogram may be valuable for identifying patients with OC who may benefit from platinum chemotherapy.

Identification of pathological grade and prognosis-associated lncRNA for ovarian cancer

Ovarian carcinoma (OC) is one of the most common malignant tumors in female genitals. In recent years, the therapeutic effect of OC has been significantly improved through the application of effective chemotherapy regimen. However, the 5-year survival rate is also lower than 30% due to high rate of relapse. So, it is needed to screen reliable predictive and prognostic markers of OC. Ovarian cancer gene expression data and corresponding clinical data used were downloaded from Gene Expression Omnibus database. Weighted gene expression network analysis (WGCNA) and Cox proportional hazards regression (PHR) were used to screen Pathological Grade and Prognosis-associated long noncoding RNA (lncRNA). Kaplan-Meier analysis and receiver operating characteristic curves analysis were performed to evaluate the predictive ability of the selected lncRNA. Gene Ontology (GO) enrichment and Gene Set Enrichment Analysis (GSEA) enrichment analysis methods were used to explore the possible mechanisms of the selected lncRNA affecting the development of OC. Five reliably lncRNAs (LINC00664, LINC00667, LINC01139, LINC01419, and LOC286437) was identified through a series of bioinformatics methods. In testing cohorts, we found that the five lncRNAs in predicting the risk of OC recurrence is robustness, and multivariate Cox PHR analysis indicate that the five lncRNAs is an independent risk factor for OC recurrence. Moreover, GO and GSEA enrichment analysis showed that the five lncRNAs are involved in multiple ovarian cancer occurrence mechanism. In summary, all these findings indicated that the five lncRNAs can effectively predict the risk of recurrence of ovarian cancer.

High Expression Levels of AGGF1 and MFAP4 Predict Primary Platinum-Based Chemoresistance and are Associated with Adverse Prognosis in Patients with Serous Ovarian Cancer

Primary platinum-based chemoresistance occurs in approximately one-third of patients with serous ovarian cancer (SOC); however, traditional clinical indicators are poor predictors of chemoresistance. So we aimed to identify novel genes as predictors of primary platinum-based chemoresistance. Gene expression microarray analyses were performed to identify the genes related to primary platinum resistance in SOC on two discovery datasets (GSE51373, GSE63885) and one validation dataset (TCGA). Univariate and multivariate analyses with logistic regression were performed to evaluate the predictive values of the genes for platinum resistance. Machine learning algorithms (linear kernel support vector machine and artificial neural network) were applied to build prediction models. Univariate and multivariate analyses with Cox proportional hazards regression and log-rank tests were used to assess the effects of these gene signatures for platinum resistance on prognosis in two independent datasets (GSE9891, GSE32062). AGGF1 and MFAP4 were found highly expressed in patients with platinum-resistant SOC and independently predicted platinum resistance. Platinum resistance prediction models based on these targets had robust predictive power (highest AUC: 0.8056, 95% CI: 0.6338-0.9773; lowest AUC: 0.7245, 95% CI: 0.6052-0.8438). An AGGF1- and MFAP4-centered protein interaction network was built, and hypothetical regulatory pathways were identified. Enrichment analysis indicated that aberrations of extracellular matrix may play important roles in platinum resistance in SOC. High AGGF1 and MFAP4 expression levels were also related to shorter recurrence-free and overall survival in patients with SOC after adjustment for other clinical variables. Therefore, AGGF1 and MFAP4 are potential predictive biomarkers for response to platinum-based chemotherapy and survival outcomes in SOC.

WWOX, the FRA16D gene: A target of and a contributor to genomic instability

WWOX is one of the largest human genes spanning over 1.11 Mbp in length at chr16q23.1-q23.2 and containing FRA16D, the second most common chromosomal fragile site. FRA16D is a hot spot of genomic instability, prone to breakage and for causing germline and somatic copy number variations (CNVs). Consequentially WWOX is frequent target for deletions in cancer. Esophageal, stomach, colon, bladder, ovarian, and uterine cancers are those most commonly affected by WWOX deep focal deletions. WWOX deletions significantly correlate with various clinicopathological features in esophageal carcinoma. WWOX is also a common target for translocations in multiple myeloma. By mapping R-loop (RNA:DNA hybrid) forming sequences (RFLS) we observe this to be a consistent feature aligning with germline and somatic CNV break points at the edges and core of FRA16D spanning from introns 5 to 8 of WWOX. Germline CNV polymorphisms affecting WWOX are extremely common in humans across different ethnic groups. Importantly, structural variants datasets allowed us to identify a specific hot spot for germline duplications and deletions within intron 5 of WWOX coinciding with the 5' edge of the FRA16D core and various RFLS. Recently, multiple pathogenic CNVs spanning WWOX have been identified associated with neurological conditions such as autism spectrum disorder, infantile epileptic encephalopathies, and other developmental anomalies. Loss of WWOX function has recently been associated with DNA damage repair abnormalities, increased genomic instability, and resistance to chemoradiotherapy. The described observations place WWOX both as a target of and a contributor to genomic instability. Both of these aspects will be discussed in this review.

Prognostic Characteristics of MACC1 Expression in Epithelial Ovarian Cancer

Recent studies have shown that overexpression of metastasis-associated in colon cancer 1 (MACC1) is significantly associated with adverse prognoses of patients with different kinds of cancer. However, the exact survival effect of MACC1 on epithelial ovarian cancer (EOC) patients has not yet been established. Thus, the objective of this study was to explore the prognostic role of MACC1 mRNA in EOC by using Kaplan-Meier (KM) plotter and ONCOMINE database. Our results indicated that MACC1 mRNA high expression was significantly associated with unfavorable overall survival (hazard ratio (HR) = 1.51 (95% confidence interval (CI): 1.21 - 1.88), P = 0.00025) and progression-free survival (HR = 1.53 (95% CI: 1.24 - 1.89), P = 5.8e-05) in EOC patients. We also found that the expression of MACC1 mRNA in EOC was 2.5 times higher than that in normal surface ovarian epithelium, which was statistically significant (P = 2.86e-7). Our results suggest that MACC1 expression might be a biomarker for poor prognosis in individual EOC patients.

Novel tumor suppressor SPRYD4 inhibits tumor progression in hepatocellular carcinoma by inducing apoptotic cell death

BACKGROUND

Hepatocellular carcinoma (HCC) is one of the leading causes of cancer-associated deaths worldwide. Although recent studies have proposed different biomarkers for HCC progression and therapy resistance, a better understanding of the molecular mechanisms underlying HCC progression and recurrence, as well as the identification of molecular markers with a higher diagnostic accuracy, are necessary for the development of more effective clinical management strategies. Here, we aimed to identify novel players in HCC progression.

METHODS

SPRYD4 mRNA and protein expression analyses were carried out on a normal liver-derived cell line (HL-7702) and four HCC-derived cell lines (HepG2, SMMC7721, Huh-7, BEL-7402) using qRT-PCR and Western blotting, respectively. Cell proliferation Cell Counting Kit-8 (CCK-8) assays, protein expression analyses for apoptosis markers using Western blotting, and Caspase-Glo 3/7 apoptosis assays were carried out on the four HCC-derived cell lines. Expression comparison, functional annotation, gene set enrichment, correlation and survival analyses were carried out on patient data retrieved from the NCBI Gene module, the NCBI GEO database and the TCGA database.

RESULTS

Through a meta-analysis we found that the expression of SPRYD4 was downregulated in primary HCC tissues compared to non-tumor tissues. We also found that the expression of SPRYD4 was downregulated in HCC-derived cells compared to normal liver-derived cells. Subsequently, we found that the expression of SPRYD4 was inversely correlated with a gene signature associated with HCC cell proliferation. Exogenous SPRYD4 expression was found to inhibit HCC cell proliferation by inducing apoptotic cell death. We also found that SPRYD4 expression was associated with a good prognosis and that its expression became downregulated when HCCs progressed towards more aggressive stages and higher grades. Finally, we found that SPRYD4 expression may serve as a biomarker for a good overall and relapse-free survival in HCC patients.

CONCLUSIONS

Our data indicate that a decreased SPRYD4 expression may serve as an independent predictor for a poor prognosis in patients with HCC and that increased SPRYD4 expression may reduce HCC growth and progression through the induction of apoptotic cell death, thereby providing a potential therapeutic target.

A panel of 7 prognosis-related long non-coding RNAs to improve platinum-based chemoresistance prediction in ovarian cancer

In order to study the role of long non-coding RNAs (lncRNAs) in predicting platinum-based chemoresistance in patients with high-grade serous ovarian carcinoma (HGS-OvCa), a=7-lncRNA signature was developed by analyzing 561 microarrays and 136 specimens from RNA-sequencing (RNA-seq) obtained from online databases [odds ratio (OR), 2.859; P<0.0001]. The upregulated lncRNAs (RP11-126K1.6, ZBED3-AS1, RP11-439E19.10 and RP11‑348N5.7) and downregulated lncRNAs [RNF144A-AS1, growth arrest specific 5 (GAS5) and F11-AS1] exhibited high sensitivity and specificity in predicting chemoresistance in the Gene Expression Omnibus and the Cancer Genome Atlas (area under curve >0.8). The lncRNA signature was independent of clinical characteristics and 4 HGS-OvCa molecular subtypes. This signature was negatively associated with disease-free survival (n=47; log-rank, P<0.01). Furthermore, the expression of the 7 lncRNAs was consistent with microarray (GSE63885, GSE51373, GSE15372 and GSE9891) and RNA-seq data. In in vitro experiments, ZBED3-AS1, F11-AS1 and GAS5 were differentially expressed in cell lines that are known to be resistant and non-resistant to platinum-based drugs, which was consistent with the results in the present study. This lncRNA signature may be used as a prognostic marker for predicting resistance to platinum-based chemotherapeutics in HGS-OvCa. These findings may contribute to individualized therapies in patients with HGS-OvCa in the future.

Curated compendium of human transcriptional biomarker data

One important use of genome-wide transcriptional profiles is to identify relationships between transcription levels and patient outcomes. These translational insights can guide the development of biomarkers for clinical application. Data from thousands of translational-biomarker studies have been deposited in public repositories, enabling reuse. However, data-reuse efforts require considerable time and expertise because transcriptional data are generated using heterogeneous profiling technologies, preprocessed using diverse normalization procedures, and annotated in non-standard ways. To address this problem, we curated 45 publicly available, translational-biomarker datasets from a variety of human diseases. To increase the data's utility, we reprocessed the raw expression data using a uniform computational pipeline, addressed quality-control problems, mapped the clinical annotations to a controlled vocabulary, and prepared consistently structured, analysis-ready data files. These data, along with scripts we used to prepare the data, are available in a public repository. We believe these data will be particularly useful to researchers seeking to perform benchmarking studies—for example, to compare and optimize machine-learning algorithms' ability to predict biomedical outcomes.

Transcriptome-based signature predicts the effect of taxol in serous ovarian cancer

Taxol is a widely used chemotherapy drug used clinically for ovarian cancer, although the response to Taxol among individuals varies due to the heterogeneity among ovarian cancer patients. In this work, we analyzed differences in the prognostic effect of gene expression and Taxol usage in the Cancer Genome Atlas (TCGA) dataset and identified specific genes associated with the Taxol effect. Using the Cox regression model, a risk model (Taxol score) was developed to assess the outcome of ovarian cancer patients who underwent chemotherapy with Taxol. According to the results, survival was significantly associated with the Taxol score. Moreover, the patients in the high and low Taxol score group had different responses to Taxol. This result was further validated in another two independent datasets. The correlation between clinicopathological indicators was also analyzed, and we determined that the Taxol score is not associated with age, pathological stage, or Taxol treatment, while there was significant correlation with tumor size and grade. Gene Set Enrichment Analysis (GSEA) showed that various signaling pathways including ECM receptor, drug metabolism and ascorbate metabolism pathways were significantly enriched in the high Taxol score group. Collectively, these results indicate that the model is robust for predicting the effectiveness of Taxol by reflecting the various cell statuses of serous ovarian carcinoma.

Advances in ovarian cancer therapy

Epithelial ovarian cancer is typically diagnosed at an advanced stage. Current state-of-the-art surgery and chemotherapy result in the high incidence of complete remissions; however, the recurrence rate is also high. For most patients, the disease eventually becomes a continuum of symptom-free periods and recurrence episodes. Different targeted treatment approaches and biological drugs, currently under development, bring the promise of turning ovarian cancer into a manageable chronic disease. In this review, we discuss the current standard in the therapy for ovarian cancer, major recent studies on the new variants of conventional therapies, and new therapeutic approaches, recently approved and/or in clinical trials. The latter include anti-angiogenic therapies, polyADP-ribose polymerase (PARP) inhibitors, inhibitors of growth factor signaling, or folate receptor inhibitors, as well as several immunotherapeutic approaches. We also discuss cost-effectiveness of some novel therapies and the issue of better selection of patients for personalized treatment.