- Open Access
A risk model of gene signatures for predicting platinum response and survival in ovarian cancer
Journal of Ovarian Research volume 15, Article number: 39 (2022)
Ovarian cancer (OC) is the deadliest tumor in the female reproductive tract. And increased resistance to platinum-based chemotherapy represents the major obstacle in the treatment of OC currently. Robust and accurate gene expression models are crucial tools in distinguishing platinum therapy response and evaluating the prognosis of OC patients.
In this study, 230 samples from The Cancer Genome Atlas (TCGA) OV dataset were subjected to mRNA expression profiling, single nucleotide polymorphism (SNP), and copy number variation (CNV) analysis comprehensively to screen out the differentially expressed genes (DEGs). An SVM classifier and a prognostic model were constructed using the Random Forest algorithm and LASSO Cox regression model respectively via R. The Gene Expression Omnibus (GEO) database was applied as the validation set.
Forty-eight differentially expressed genes (DEGs) were figured out through integrated analysis of gene expression, single nucleotide polymorphism (SNP), and copy number variation (CNV) data. A 10-gene classifier was constructed which could discriminate platinum-sensitive samples precisely with an AUC of 0.971 in the training set and of 0.926 in the GEO dataset (GSE638855). In addition, 8 optimal genes were further selected to construct the prognostic risk model whose predictions were consistent with the actual survival outcomes in the training cohort (p = 9.613e-05) and validated in GSE638855 (p = 0.04862). PNLDC1, SLC5A1, and SYNM were then identified as hub genes that were associated with both platinum response status and prognosis, which was further validated by the Fudan University Shanghai cancer center (FUSCC) cohort.
These findings reveal a specific risk model that could serve as effective biomarkers to identify patients’ platinum response status and predict survival outcomes for OC patients. PNLDC1, SLC5A1, and SYNM are the hub genes that may serve as potential biomarkers in OC treatment.
Ovarian cancer (OC), the most lethal gynecological cancer, is one of the main causes of cancer-related death among females worldwide . The five-year overall survival rate of epithelial OC patients ranges from 20% at stage IV to 89% at stage I, however, 80% of OC cases can not be diagnosed timely until the tumor has progressed to advanced stages with severe clinical outcomes due to its insidious onset without specific clinical manifestations and the lack of mature early diagnosis methods . Cytoreductive surgery followed by chemotherapy based on platinum or combined with taxanes is the standard treatment for OC . Although most patients with OC show initially highly response to platinum therapy, tumors demonstrate increasing resistance during treatment inevitably. Reportedly, about 70% of patients suffer from tumor relapse a few months after treatment and develop resistance to therapy eventually, no matter primary or secondary resistance, representing the major challenge in OC treatment [4, 5]. Identification of nonresponders is an important step toward greater life expectancy for OC patients . Meanwhile, the specific biomarkers predicting platinum therapy response remain obscure, therefore, it is of vital importance to figure out the potential indicators, which could aid clinical decisions and improve prognosis.
Nowadays, the rapid development of next-generation sequencing (NGS) has revolutionized and renewed how we comprehend cancer treatment and promoted the progress of precision medicine . Increasing evidence has authenticated that molecular biomarkers contribute to the prognosis evaluation and prediction of tumors . Besides, researchers found that rather than conventional single-gene biomarkers, gene signatures containing several genes can provide stronger evidence to prognosis and survival . For example, based on the public database, a six-gene model (TGFBI, SFRP1, COL16A1, THY1, PPIB, BGN) was built and serves as an independent prognostic biomarker for overall survival . Bi et al. identified eight glycolysis-related prognostic genes that effectively predicted survival in ovarian cancer ; A tumor mutation burden (TMB) associated immune risk score signature was built by Cui et al. for TMB and prognosis evaluation ; Salinas et al. applied SNP data from TCGA to find SNPs associated with chemo-response in ovarian cancer ; And another study developed and validated an immune-related gene signature that was significantly associated with survival . Despite encouraging developments, no biomarkers for the prediction of response to therapy and prognosis are applied into clinical practice yet.
In this study, through bioinformatics data analysis, we integrated the gene expression profiles of the transcriptome, single nucleotide polymorphism (SNP), and copy number variation (CNV) to identify differentially expressed genes (DEGs) firstly. Then a support vector machine classifier was constructed to distinguish patients’ responses to platinum therapy. Next, combined with L1-LASSO and Cox-Proportional Hazards regression, we constructed a prognostic risk model based on 8 optimum genes to predict prognosis which could mirror the prognosis related to platinum response status as well. Finally, after the intersection of the classifier and prognostic model, PNLDC1, SLC5A1, and SYNM were identified as the hub genes related to both platinum response and prognosis, which was further verified by IHC analysis. The flowchart of this study was displayed in Fig. 1. In summary, the comprehensive analysis of gene expression level, SNP and CNV in our study could provide more accurate and robust molecular markers for diagnosis, prediction and bring new insights into clinical treatment strategies for OC patients.
Materials and methods
OC datasets extraction
The OC datasets used in this study were derived from TCGA and GEO databases. mRNA-seq expression profile data (platform: Illumina HiSeq 2000 RNA Sequencing) and SNP, CNV information (platform: Affymetrix Genome-Wide Human SNP Array 6.0) were downloaded from TCGA. Collectively, 419 OC tumor tissue samples with expression profiles and 481 samples with SNP information were included. After barcode matching, 230 OC samples possessing clinical platinum response status were obtained, comprised of 69 resistant and 161 sensitive samples respectively, which serve as our training dataset. As for validation dataset, the gene expression profiles of GSE63885(n = 101, platform: GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array) from the GEO database (https://www.ncbi.nlm.nih.gov/geo/) was acquired, containing a total of 75 samples having platinum response status, consisting of 34 resistant and 41sensitive samples, respectively. Both the expression profiles and clinical characteristics can be obtained publicly, so there was no need to acquire ethics committee approval. The abovementioned data were displayed in Table 1.
For the original gene expression profile FPKM data downloaded from TCGA, preprocessCore Version 1.40.0  (http://bioconductor.org/packages/release/bioc/html/preprocessCore.html) in R was used to perform standardization based on the quantiles algorithm. For the SNP 6.0 chip data, PICNIC software  (ftp://ftp.sanger.ac.uk/pub/cancer) was applied to convert and process the data of CEL format to obtain CN segment data. (The segment data indicated the copy value in the detection region. Usually, the segment value of the diploid was 0, implying there was no copy number variation, and the other non-zero signal indicated the region was missing or amplificated). The human gene annotation file (Release 27 (GRCh38.p10)) from the GENCODE database (http://www.gencodegenes.org/releases/current.html) was extracted and gene annotation was employed. And at the same time, the oligo package Version 1.42.0  (http://www.bioconductor.org/packages/release/bioc/html) was used to convert to the original data of the GSE63885 data set, fill the missing data (median method), conduct background correction (MAS method) and data standardization (quantiles).
After normalization, in the light of platinum response status, we divide the patients into two groups: chemotherapy-sensitive (n = 161) and chemotherapy-resistant (n = 69). Differently expressed genes and genes with different CN signals were selected based on the ‘limma’ package . Significant DEGs were defined as those with adjusted P < 0 .05 and |log FC| ≥1. The overlapping genes, that is, genes with significant differences in both expression level and CN signal value between the resistant and sensitive groups were further screened out. Next, only genes containing CN variant sites were retained (variant types include SNP, INS, DEL, etc.). We used DAVID version6.8  (https://david.ncifcrf.gov/) to analyze the DEGs for molecular function (MF), cellular component (CC), biology procedure (BP) enrichment by studying the Gene Ontology (GO) terms. And the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of DEGs was carried out as well. A p-value < 0.05 was set as the threshold for significant enrichment. And these genes with significantly different CNV expressing were used for further analysis.
Construction and validation of the classifier
Selection of best representative gene features using Random Forest
We identified the best combination of representative genes using a Random Forest machine learning method. The Random Forest method is an ensemble algorithm comprised of a series of decision trees. Each tree randomly selects several features in the sample zone to make a prediction. These predictions will be aggregated, and the final prediction will be decided using a voting method which refers to the category having the highest votes.
We implemented the Random Forest model using the randomForest R package Version 4.6–12  (https://cran.r-project.org/web/packages/randomForest/index.html). The model was built on expression levels of genes identified in TCGA samples. Details of the algorithm were as follows:
We randomly sampled k samples from TCGA samples with replacement using a bootstrap method to construct k regression trees for classification. The unselected samples in each round constructed k out-of-bag (OOB) sets (k was iteratively set from 1 to the total number of samples N).
For a total number of n features, we randomly selected m features at each splitting node of each tree and calculated the predicting power of each feature. We then exploited the most powerful feature to assign samples at that node. N was set from 1 to the total number of variables and m was set to the secondary square root of the total number of variables.
We let each tree grow to the maximum without any pruning.
We aggregated all decision trees to construct a Random Forest (RF) model. The RF model adopted a voting method that defined the category with the highest votes as the final classification.
We evaluated the RF model using the OOB error rate and features of the model with the lowest OOB error rate were selected as the optimal combination.
Development and validation of the SVM classifier
We used the e1071 R package (https://cran.r-project.org/web/packages/e1071) to develop an SVM model  (Support Vector Machine). SVM is a supervised machine learning method for classification. Using representative features of each sample, the model predicts the possibility belonging to a certain category to implement classification. We developed the SVM model on the TCGA training set, using “Sigmoid” as the Kernel and the optimal signature genes as features. The parameters were selected by 100-fold cross-validation. To test the performance of our model, we measured five indicators on a validation set GSE63885, including sensitivity, specificity, positive prediction value (PPV), negative prediction value (NPV) and area under curves (AUC) of the receiver operating characteristic curve (ROC) . Calculation methods of these indicators are listed below:
Construction of the prognostic risk model
Selection of prognostic genes and clinical factors
We performed univariate Cox-PH (proportional hazards) regression model to select prognostic genes and clinical factors. Based on the expression levels of genes in 2.3 of TCGA samples, the model was built by the survival R package (Version 2.41–1) . Prognostic genes and clinical factors were identified using P < 0.05 as the threshold (log-rank test).
Selection of the optimal genes
Based on expression levels of the prognostic genes identified prior, we implemented L1-Regularized Cox-PH regression analysis on TCGA samples to select the optimal combination of these genes . The model was developed using the penalized R package (Version0.9–50, http://bioconductor.org/packages/penalized/) . The “lambda” parameter was identified by the 1000 cross-validation likelihood (cvl) method.
Development of a Cox-PH model using optimal genes
We constructed a prognosis score (PS) using coefficients of optimal genes from the Cox-PH regression model. Using the median of PS as the threshold, we separate training samples into high and low-risk groups. The prognostic value of PS was then evaluated by the Kaplan-Meier survival curve  using log-rank test in the training set and then validated in the GSE63885 dataset.
Hub genes screening and validation
Selection of the hub genes related to platinum therapy and prognosis
Further comparison between feature genes comprised in SVM classifier and included in the prognostic risk prediction model, crucial genes were screened out.
Validation via immunohistochemistry
A total of 80 OC patients from FUSCC who received platinum-based chemotherapy after surgery were selected randomly, and tissue microarrays composed the representative cores from each specimen. Immunohistochemistry proceeded as described before . In brief, specimens were incubated first with an anti- PNLDC1 antibody (1:2000, Proteintech, China), anti-SLC5A1 antibody (1:1000, Abcam, UK) or anti-SYNM antibody (1:200, Proteintech, China) overnight at 4 °C and then with a biotinylated secondary antibody (1:100, goat anti-rabbit IgG) for 30 min at 37 °C. A well-established IRS was then used to calculate the protein expression level of these three hub genes . Firstly, the staining intensity (SI) was scored using a 4-point scale from 0 to 3, with 0 if there was no staining. For weak, moderate, and strong staining, the scores were 1, 2 and 3, respectively. Secondly, the percentage of positive cells was scored into five categories: no staining, 1–10, 11–50, 51–80, 81–100 percentage positive cells. And the scores were 0, 1, 2, 3 and 4, respectively. An IRS was calculated by multiplying the percentage of hub genes by the SI score, resulting in a scale from 0 to 12. The IRS was divided into four groups: 0 (IRS 0–1), 1 (IRS 2–3), 2 (IRS 4–8) and 3 (IRS 8–12). Then, 0 and 1 were stratified into low expression group and 2 and 3 into high expression group and performed survival analyses. The expression of each hub gene was quantified by using an Image-Pro Plus Image Analysis Software and the IOD (Integral optic density) was measured as reported previously . To find the optimal cut-off points, the X-tile program was used .
All statistical analyses in this study were performed by R (version 3.4.1). The statistical significance threshold was set at 0.05 if not explicitly mentioned. In our study, progression-free survival (PFS) is defined as the time from operation to relapse or progression, whichever occurred first. And patients were divided into platinum-sensitive and platinum-resistant subgroups according to the platinum-free interval (PFI). PFI is defined as the time from the end of the first chemotherapy course to disease recurrence, and PFI > 6 months was regarded as platinum-sensitive, whereas PFI < 6 months was the platinum-resistant group.
Data pre-processing and DEGs screening
The gene expression profiles obtained from TCGA and GEO datasets were firstly normalized and the box diagram before and after standardization was shown in Supplementary Figure 1. Concerning the data of CN signal, gene annotation was performed, followed by depicting the distribution of chromosomes. And we found that, in different samples, but the same sites, the CN signals distributed similarly (Fig. 2A). Via Limma package in R, 1144 DEGs differently expressed (524 downregulated and 620 upregulated) and 1864 DEGs with diversified CN signals (727 downregulated and 1137 upregulated) between platinum-sensitive and the resistant group were obtained from the TCGA database (Fig. 2B and C, Supplementary Table 1).
Further analysis uncovered 108 genes as overlapping genes in both expression and CN signal levels (Fig. 2D, Supplementary Table 2). Integrated with SNP information, we found 48 genes had variant sites, including 94 SNP,1 INS, and 1 DEL (Supplementary Table 3), indicating they were differentially expressed and had diversified CNV signals between platinum-resistant and sensitive groups simultaneously. To reveal the biological functions of the 48 DEGs, the GO and KEGG enrichment analyses conducted by DAVID were employed (Fig. 2E). Regarding biological process (BP), the GO analysis results showed that the intersecting DEGs were mainly enriched in terms related to cell adhesion (Supplementary Table 4). As for KEGG pathway analysis, the DEGs were enriched in Cell adhesion molecules (Table 2).
Construction and validation of the classifier
To extract the most representative and feature genes in these 48 genes, the RandomForest algorithm was performed. And 10 genes were figured out optimally when the OOB error is minimum (Fig. 3A). And we found that the CNV of the 10 genes were all SNPs, comprising one known SNP site in CD209 and 9 unreported SNP sites (Table 3).
Based on expression profiles of the above 10 feature genes in the TCGA dataset, we constructed an SVM classifier to determine the platinum-sensitive and resistant samples. It was able to accurately distinguish 212 out of 230 samples (161 sensitivity vs. 51 resistance), with a precision rate of 92.17% and an average AUC of 0.971 (Fig. 3B, solid line). The sensitivity and specificity were 1 and 0.839, respectively, and the PPV and NPV are 0.899 and 1 as well. To further verify and evaluate the predictive effects of this model, GSE63885 was used as an independent external validation dataset. The result of the validation cohort showed that 71 samples (41sensitivity vs. 30 resistance) out of 75 samples could be discriminated precisely, and the accuracy rate was 94.76% with an AUC of 0.926 (Fig. 3B, dotted line). The sensitivity, specificity, PPV, and NPV were 1, 0.882, 0.899, and 1, respectively. To sum up, this model could accurately classify the drug-sensitive samples and moderately found the drug-resistant patients, indicating these 10 genes had strong correlations with drug sensitivity. The gene profiles of 10 feature genes in the TCGA and GSE638855 datasets were displayed in Supplementary Table 5.
Construction and validation of the prognostic risk model
Univariate cox regression analysis
Combined with the clinical information, the overlapping 48 DEGs were filtered via univariate cox regression analysis in the TCGA training cohort to acquire genes significantly related to prognosis. Consequently, 34 genes were differentially expressed and 29 genes with different CN signals were obtained separately (Supplementary Table 6). After the intersection, 20 crossed genes were left (Fig. 4A). Meanwhile, the clinicopathological factors related to prognosis identified by univariate analysis were merely platinum response status (Table 4). And a conspicuous OS difference was noted between sensitive and resistant Kaplan-Meier curve (HR = 0.22, p = 5.55e-16), reflecting a better survival in the platinum-sensitive group (Fig. 4B). It also verified that the genes we screened were indeed related to the platinum response status to some extent.
Selection of the optimal genes
According to the intersectional 20 genes related to prognosis, the Cox-PH model based on the L1-penalized regularized regression algorithm was exploited to further select the optimal genes. The maximum value of cvl − 771.2244 was obtained when the lambda value is 20.88803 after 1000 cycles of cvl algorithm calculations (Fig. 5A). Under this Circumstance, 8 optimum genes were received (Table 5) and the gene prognostic coefficients are shown in Fig. 5B.
Construction and validation of the prognostic risk model
Based on the Cox-PH prognostic coefficients of the 8 optimized genes, a risk model was constructed by the following formula: Prognosis score (PS) = (− 0.42542) × ExpGJA8+ (0.430375) × ExpPNLDC1 + (− 0.20707) × Exp SLC5A1 + (1.169891) × ExpVSTM2L +(1.195075) × Exp CACNA1C +(− 1.64918) × ExpSEZ6L+(0.442726) × ExpGDF3 + (− 1.78725) × Exp SYNM.
To validate the survival-predicting performance of the model, the prognostic score (PS) of each sample was calculated and the median PS was applied as the threshold to subdivide the training cohort into a high-risk group (HRG) and a low-risk group (LRG). First, in the training set, the correlation between the model’s predictions and the actual prognosis was evaluated through the Kaplan-Meier curve. We discovered that LRG had a longer median OS time than HRG. In detail, the median OS of HRG (115 samples) was 38.56 ± 21.31 months, while the average OS of the LRG (115 samples) was 47.30 ± 26.11 months (Fig. 5C). And the correlation between the groupings predicted by the model and the actual survival outcome was significant and consistent (p = 9.613e-05). Concurrently, the results of the validation set GSE63885 saw identical results, showing that the average survival time of the HRG (35 samples) was 33.78 ± 18.33 months, whereas the LRG (35 samples) had a longer median OS of 46.11 ± 30.49 months (Fig. 5D). Model predictions and actual results had a significant correlation (p = 0.04862). The detailed survival information of the TCGA and GSE63885 datasets and the PS information of the samples were shown in Supplementary Table 7.
Additionally, to further clarify the correlation between the prediction of the prognostic model and the platinum response status, subgroup survival analysis dividing into the sensitive and resistant groups was carried out (Fig. 5E). On the whole, the prognosis predicted by the prognostic model is significantly correlated with the prognosis that depends on the platinum response state (p = 1.694e-08). Specifically, in the sensitive subgroup, the prognosis prediction of the model was remarkably correlated with the actual survival. It is supposed that after being determined by the SVM classifier, patients who identified as the sensitive group could accept platinum-based chemotherapy continually and the prognostic model could foresee the prognosis accurately. Conversely, patients who were considered to be resistant could alter and optimize therapeutic strategies as early as possible. In a nutshell, our model had an important role in different ways for the two groups which helped with the clinical decision to some degree.
Identification of hub genes related to platinum response and prognosis
To find hub genes, genes used to establish the SVM classifier and the prognostic model were consolidated, and 3 intersecting genes named PNLDC1, SLC5A1, and SYNM were obtained (Fig. 6A). And we found that, both in training and validation sets, PNLDC1 was down-regulated in the sensitive group while SLC5A1 and SYNM were up-regulated (Fig. 6B and C; Supplementary Table 8), hence we assumed that these 3 genes reflect both platinum response status and prognosis, and their CNV information is shown in Table 6.
To confirm their relationships with platinum response status and prognosis, IHC of 80 patients from the FUSCC cohort was applied (Fig. 6D). And in the FUSCC cohort, 20 patients were resistant to platinum-based therapy and the other 60 were in the sensitive group (Supplementary Table 9). In line with our former findings, the intensity and quantity of PNLDC1’s expression were remarkably higher in the resistant group (p = 0.0096), while SLC5A1 (p = 0.0058) and SYNM (p = 0.0022) were significantly amplified in the sensitive patients (Fig. 6E). Since the deficient number of events regarding the OS analysis, PFS was analyzed instead. X-tile was exploited to find the best cut-off points with minimum p values. The K-M curve uncovered that those patients with high expression of PNLDC1 had shorter PFS (Log-rank p = 0.009, Fig. 7A). Conversely, survival analysis demonstrated that the SLC5A1 and SYNM high group showed adverse outcomes (Log-rank p < 0.001, Fig. 7B; Log-rank p = 0.0015, Fig. 7C). These results further supported that PNLDC1, SLC5A1 and SYNM were the hub genes associated with the platinum response status and prognosis indicating they could be used as potential biomarkers in clinical practice.
Nowadays, the common treatment regimen for OC consists of tumor debulking, followed by administration of platinum-based chemotherapy , however, resistance to platinum therapy limits therapeutic options, and makes platinum-resistant patients the most challenging to treat. Apart from PFI, regarding as a predictive factor to subsequent platinum therapy , biomarkers reflecting platinum response status are urgently needed. Therefore, to begin with, the 230 samples in the TCGA database were stratified into the sensitive and resistant groups, and through comprehensive analysis from 3 aspects, including gene expression level, CN signal, and SNP data, we obtained 48 overlapping DEGs. Considering a burgeoning number of researches concerning gene signature of tumors show up, which was supported by the development of RNA-sequencing and microarray, as well as available public databases [14, 32]. We constructed a classifier via the RandomForest algorithm aiming at dividing patients into sensitivity and resistance groups. The classifier comprised of 10 genes (CD209, CD274, HIST1H3I, HIST1H4L, NLGN1, NTRK3, PNLDC1, SLC22A3, SLC5A1 and SYNM), and its validation in the GEO dataset showed a satisfying consistency, which indicating the classifier could differentiate the sensitive group accurately and aid the resistant patients to receive other effective therapy as soon as possible.
Liu et al. used the TCGA dataset to validate a seven genes-based model which can predict the survival of FIGO stage IIIc serous ovarian carcinoma (HG3cSOC) and served as a valuable marker for the response to platinum-based chemotherapy , however, they chose HG3cSOC to analysis and confined to gene expression analysis only; Zhao et al. identified AGGF1 and MFAP4 as potential predictors of primary platinum-based chemoresistance , nonetheless, they just focused on the gene expression level and lack of experimental-level validation, such as immunohistochemistry; Dugo et al. mainly focused on HGSOC patients who received complete cytoreduction (R0) and analyzed focal copy number alterations ; A qualitative transcriptional signature for predicting recurrence risk for high-grade serous ovarian cancer patients was constructed by Liu et al. ; Salinas et al. found 19 SNPs were associated with chemo-response . Although these studies have similarities to ours, we conducted a comprehensive analysis from diversified levels to figure out the most robust markers indicating platinum treatment response and prognosis.
Detailedly, the current treatment of OC could prolong the interval between recurrences but does not benefit overall survival . In the study, via univariate Cox regression analyses, genes and clinicopathologic parameters related to prognosis were obtained. And we found that sensitivity to platinum-based therapy was the only clinical factor that contribute to better survival. To make survival predictions, 8 optimum genes were applied to establish the prognostic model, and we found that upregulation of PNLDC1, VSTM2L, CACNA1C, and GDF3 were related to worse clinical outcomes, however, high expressions of GJA8, SEZ6L, SLC5A1, and SYNM were associated with better prognoses. Specifically, PNLDC1, a PARN-like 3′-to-5′ exonuclease located at the membrane of the mitochondria in a mouse, is critical to the processing of piRNAs and PNLDC1 disrupted in mice would lead to azoospermia and male infertility ultimately [38, 39]. Former studies pointed that PNLDC1 was related to survival in CRC patients , and PNLDC1 expressed higher in normal colorectal tissues than in cancer tissues . VSTM2L could induce adverse survival outcomes in rectal cancer . CACNA1C, as a voltage-gated calcium channel, is up-regulated in brain tumors, leukemia, breast cancer, and other tumors  and plays as an oncogene in OC tumors . GDF3 is widely accepted as a pluripotency marker and expressed in several cancer types such as breast carcinoma  melanoma . GJA8 is amplified in Wilms tumors . SEZ6L plays a role in signal transduction and protein-protein interaction and increases in lung cancer . SLC5A1, a member of the GLUT family, encodes SGLT1 facilitating glucose transport at the basolateral membrane of the cells [49, 50]. Aberrant expression of SLC5A1 in different types of human cancers is observed, including ovarian cancer , cervical cancer , colorectal cancer [53, 54], hepatocellular carcinoma ,prostate cancer . SYNM is a type IV intermediate filament , which is supposed to modulate biological processes such as cell adhesion, cell motility. Reportedly, the silencing of synemin results in the suppression of tumor proliferation , whereas its hypomethylation is associated with aggressiveness in breast cancer . Given all these data, the contradictory function of a single gene may ascribe to the “dual roles” of genes in different cancers, even in diversified stages of one cancer type.
Apart from predicting survival, the potential of prognostic classifiers lies in the ability to recognize patients that are more likely to respond to particular therapies . Interestingly, in our study, the predictions of the prognostic risk model were significantly associated with the platinum response status, especially in the sensitive population, indicating the two models we built had a prominent correlation and may serve the clinical practice.
Finally, we screened out 3 hub genes, namely PNLDC1, SLC5A1 and SYNM, and explored their relationship with platinum therapy response and prognosis. The expression level of these genes detected by IHC in the FUSCC cohort indicates that PNLDC1 expresses higher in the resistant group, whereas overexpression of the SLC5A1 and SYNM were detected in the sensitive patients. Despite the insufficient number of events in OS analysis, all these 3 genes show significant correlations with PFS. Longer and more detailed follow-up data is required in further study. To date, the role of PNLDC1 and SYNM in OC has not been reported yet, and seldom do SLC5A1 as well, indicating they might be treated as new potential biomarkers in platinum-based chemotherapy and prognosis in OC, even as the therapeutic targets to some extent. For example, the development of individual SGLT1 inhibitors which target SLC5A1 is on the way .
Although our study brings new insights into OC treatment and survival, there are limitations to this study. Firstly, despite we included two completely independent datasets as training and validation cohorts and included an IHC cohort to validate our findings, this is a retrospective study. Next, the sample sizes in the training set, validation set, and IHC cohort were relatively small. Larger cohorts are required to prove our findings. Lastly, although bioinformatic analysis is a powerful tool for exploring potential biomarkers in tumors, further corroboration is needed.
Taken together, based on the public databases, considering multiple aspects, we selected the feature DEGs to construct the SVM classifier to determine the patients’ responses to platinum-based chemotherapy. Meanwhile, a prognostic risk model was established to help predict patients’ prognoses. And 3 hub genes, PNLDC1, SLC5A1 and SYNM, related to the platinum therapy response and prognosis were screened out, which could be used as new biomarkers in OC treatment.
Availability of data and materials
The datasets used and/or analyzed in this study are available from the corresponding author on reasonable request.
The Cancer Genome Atlas
Gene Expression Omnibus
Single nucleotide polymorphism
Copy number variation
Differentially expressed genes
Negative prediction value
Area under curves
Receiver operating characteristic curve
Integral optic density
Kyoto Encyclopedia of Genes and Genomes
Tumor Mutation Burden
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Torre LA, Trabert B, DeSantis CE, Miller KD, Samimi G, Runowicz CD, et al. Ovarian cancer statistics, 2018. CA Cancer J Clin. 2018;68:284–96.
Pignata S, Cannella L, Leopardo D, Pisano C, Bruni GS, Facchini G. Chemotherapy in epithelial ovarian cancer. Cancer Lett. 2011;303:73–83.
Davis A, Tinker AV, Friedlander M. "platinum resistant" ovarian cancer: what is it, who to treat and how to measure benefit? Gynecol Oncol. 2014;133:624–31.
Matsuo K, Lin YG, Roman LD, Sood AK. Overcoming platinum resistance in ovarian carcinoma. Expert Opin Investig Drugs. 2010;19:1339–54.
Kelland L. The resurgence of platinum-based cancer chemotherapy. Nat Rev Cancer. 2007;7:573–84.
Nakagawa H, Fujita M. Whole genome sequencing analysis for cancer genomics and precision medicine. Cancer Sci. 2018;109:513–22.
Bonome T, Levine DA, Shih J, Randonovich M, Pise-Masison CA, Bogomolniy F, et al. A gene signature predicting for survival in suboptimally debulked patients with ovarian cancer. Cancer Res. 2008;68:5478–86.
Raman P, Maddipati R, Lim KH, Tozeren A. Pancreatic cancer survival analysis defines a signature that predicts outcome. PLoS One. 2018;13:e0201751.
Pan X, Ma X. A novel six-gene signature for prognosis prediction in ovarian Cancer. Front Genet. 2020;11:1006.
Bi J, Bi F, Pan X, Yang Q. Establishment of a novel glycolysis-related prognostic gene signature for ovarian cancer and its relationships with immune infiltration of the tumor microenvironment. J Transl Med. 2021;19:382.
Cui M, Xia Q, Zhang X, Yan W, Meng D, Xie S, et al. Development and validation of a tumor mutation burden-related immune prognostic signature for ovarian cancers. Front Genet. 2021;12:688207.
Salinas EA, Newtson AM, Leslie KK, Gonzalez-Bosquet J. Association analysis of a chemo-response signature identified within the Cancer genome atlas aimed at predicting genetic risk for chemo-response in ovarian cancer. Int J Mol Epidemiol Genet. 2016;7:41–4.
Zhang B, Nie X, Miao X, Wang S, Li J, Wang S. Development and verification of an immune-related gene pairs prognostic signature in ovarian cancer. J Cell Mol Med. 2021;25:2918–30.
Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–93.
Greenman CD, Bignell G, Butler A, Edkins S, Hinton J, Beare D, et al. PICNIC: an algorithm to predict absolute allelic copy number variation with microarray cancer data. Biostatistics. 2010;11:164–75.
Parrish RS, Spencer HJ 3rd. Effect of normalization on significance testing for oligonucleotide microarrays. J Biopharm Stat. 2004;14:575–89.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43:e47.
Huang da W, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 2009;37:1–13.
Gardner IA, Greiner M. Receiver-operating characteristic curves and likelihood ratios: improvements over traditional methods for the evaluation and application of veterinary clinical pathology tests. Vet Clin Pathol. 2006;35:8–17.
Wang Q, Liu X. Screening of feature genes in distinguishing different types of breast cancer using support vector machine. Onco Targets Ther. 2015;8:2311–7.
Wang P, Wang Y, Hang B, Zou X, Mao JH. A novel gene expression-based prognostic scoring system to predict survival in gastric cancer. Oncotarget. 2016;7:55343–51.
Tibshirani R. The lasso method for variable selection in the cox model. Stat Med. 1997;16:385–95.
Goeman JJ. L1 penalized estimation in the cox proportional hazards model. Biom J. 2010;52:70–84.
Goel MK, Khanna P, Kishore J. Understanding survival analysis: Kaplan-Meier estimate. Int J Ayurveda Res. 2010;1:274–8.
Camp RL, Chung GG, Rimm DL. Automated subcellular localization and quantification of protein expression in tissue microarrays. Nat Med. 2002;8:1323–7.
Specht E, Kaemmerer D, Sanger J, Wirtz RM, Schulz S, Lupp A. Comparison of immunoreactive score, HER2/neu score and H score for the immunohistochemical evaluation of somatostatin receptors in bronchopulmonary neuroendocrine neoplasms. Histopathology. 2015;67:368–77.
Jin GZ, Dong H, Yu WL, Li Y, Lu XY, Yu H, et al. A novel panel of biomarkers in distinction of small well-differentiated HCC from dysplastic nodules and outcome values. BMC Cancer. 2013;13:161.
Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10:7252–9.
Markman M, Liu PY, Wilczynski S, Monk B, Copeland LJ, Alvarez RD, et al. Phase III randomized trial of 12 versus 3 months of maintenance paclitaxel in patients with advanced ovarian cancer after complete response to platinum and paclitaxel-based chemotherapy: a southwest oncology group and gynecologic oncology group trial. J Clin Oncol. 2003;21:2460–5.
Bolis G, Scarfone G, Luchini L, Ferraris C, Zanaboni F, Presti M, et al. Response to second-line weekly cisplatin chemotherapy in ovarian cancer previously treated with a cisplatin- or carboplatin-based regimen. Eur J Cancer. 1994;30A:1764–8.
D'Alessandris N, Travaglino A, Santoro A, Arciuolo D, Scaglione G, Raffone A, et al. TCGA molecular subgroups of endometrial carcinoma in ovarian endometrioid carcinoma: a quantitative systematic review. Gynecol Oncol. 2021.
Liu G, Chen L, Ren H, Liu F, Dong C, Wu A, et al. Seven genes based novel signature predicts clinical outcome and platinum sensitivity of high grade IIIc serous ovarian carcinoma. Int J Biol Sci. 2018;14:2012–22.
Zhao H, Sun Q, Li L, Zhou J, Zhang C, Hu T, et al. High expression levels of AGGF1 and MFAP4 predict primary platinum-based Chemoresistance and are associated with adverse prognosis in patients with serous ovarian Cancer. J Cancer. 2019;10:397–407.
Dugo M, Devecchi A, De Cecco L, Cecchin E, Mezzanzanica D, Sensi M, et al. Focal recurrent copy number alterations characterize disease relapse in high grade serous ovarian Cancer patients with good clinical prognosis: a pilot study. Genes. 2019;10.
Liu Y, Zhang Z, Li T, Li X, Zhang S, Li Y, et al. A qualitative transcriptional signature for predicting recurrence risk for high-grade serous ovarian Cancer patients treated with platinum-Taxane adjuvant chemotherapy. Front Oncol. 2019;9:1094.
Eisenhauer EA. Real-world evidence in the treatment of ovarian cancer. Ann Oncol. 2017;28: viii61-viii5.
Izumi N, Shoji K, Suzuki Y, Katsuma S, Tomari Y. Zucchini consensus motifs determine the mechanism of pre-piRNA production. Nature. 2020;578:311–6.
Nagirnaja L, Morup N, Nielsen JE, Stakaitis R, Golubickaite I, Oud MS, et al. Variant PNLDC1, defective piRNA processing, and Azoospermia. N Engl J Med. 2021;385:707–19.
Li T, Hui W, Halike H, Gao F. RNA binding protein-based model for prognostic prediction of colorectal Cancer. Technol Cancer Res Treat. 2021;20:15330338211019504.
Miao Y, Zhang H, Su B, Wang J, Quan W, Li Q, et al. Construction and validation of an RNA-binding protein-associated prognostic model for colorectal cancer. PeerJ. 2021;9:e11219.
Liu H, Zhang Z, Zhen P, Zhou M. High expression of VSTM2L induced resistance to Chemoradiotherapy in rectal Cancer through downstream IL-4 signaling. J Immunol Res. 2021;2021:6657012.
Wang CY, Lai MD, Phan NN, Sun Z, Lin YC. Meta-analysis of public microarray datasets reveals voltage-gated calcium gene signatures in clinical Cancer patients. PLoS One. 2015;10:e0125766.
Chang X, Dong Y. CACNA1C is a prognostic predictor for patients with ovarian cancer. J Ovarian Res. 2021;14:88.
Li Q, Ling Y, Yu L. GDF3 inhibits the growth of breast cancer cells and promotes the apoptosis induced by Taxol. J Cancer Res Clin Oncol. 2012;138:1073–9.
Ehira N, Oshiumi H, Matsumoto M, Kondo T, Asaka M, Seya T. An embryo-specific expressing TGF-beta family protein, growth-differentiation factor 3 (GDF3), augments progression of B16 melanoma. J Exp Clin Cancer Res. 2010;29:135.
Krepischi ACV, Maschietto M, Ferreira EN, Silva AG, Costa SS, da Cunha IW, et al. Genomic imbalances pinpoint potential oncogenes and tumor suppressors in Wilms tumors. Mol Cytogenet. 2016;9:20.
Gorlov IP, Meyer P, Liloglou T, Myles J, Boettger MB, Cassidy A, et al. Seizure 6-like (SEZ6L) gene and risk for lung cancer. Cancer Res. 2007;67:8406–11.
Mueckler M, Thorens B. The SLC2 (GLUT) family of membrane transporters. Mol Asp Med. 2013;34:121–38.
Roder PV, Geillinger KE, Zietek TS, Thorens B, Koepsell H, Daniel H. The role of SGLT1 and GLUT2 in intestinal glucose transport and sensing. PLoS One. 2014;9:e89977.
Lai B, Xiao Y, Pu H, Cao Q, Jing H, Liu X. Overexpression of SGLT1 is correlated with tumor development and poor prognosis of ovarian carcinoma. Arch Gynecol Obstet. 2012;285:1455–61.
Perez M, Praena-Fernandez JM, Felipe-Abrio B, Lopez-Garcia MA, Lucena-Cacace A, Garcia A, et al. MAP 17 and SGLT1 protein expression levels as prognostic markers for cervical tumor patient survival. PLoS One. 2013;8:e56169.
Mojica L, Luna-Vital DA. Gonzalez de Mejia E. black bean peptides inhibit glucose uptake in Caco-2 adenocarcinoma cells by blocking the expression and translocation pathway of glucose transporters. Toxicol Rep. 2018;5:552–60.
Guo GF, Cai YC, Zhang B, Xu RH, Qiu HJ, Xia LP, et al. Overexpression of SGLT1 and EGFR in colorectal cancer showing a correlation with the prognosis. Med Oncol (Northwood, London, England). 2011;28(Suppl 1):S197–203.
Lei S, Yang J, Chen C, Sun J, Yang L, Tang H, et al. FLIP(L) is critical for aerobic glycolysis in hepatocellular carcinoma. J Exp Clin Cancer Res. 2016;35:79.
Ren J, Bollu LR, Su F, Gao G, Xu L, Huang WC, et al. EGFR-SGLT1 interaction does not respond to EGFR modulators, but inhibition of SGLT1 sensitizes prostate cancer cells to EGFR tyrosine kinase inhibitors. Prostate. 2013;73:1453–61.
Granger BL, Lazarides E. Synemin: a new high molecular weight protein associated with desmin and vimentin filaments in muscle. Cell. 1980;22:727–38.
Kapoor S. Synemin: an evolving role in tumor growth and progression. J Cachexia Sarcopenia Muscle. 2014;5:347–8.
Noetzel E, Rose M, Sevinc E, Hilgers RD, Hartmann A, Naami A, et al. Intermediate filament dynamics and breast cancer: aberrant promoter methylation of the Synemin gene is associated with early tumor relapse. Oncogene. 2010;29:4814–25.
Verhaak RG, Tamayo P, Yang JY, Hubbard D, Zhang H, Creighton CJ, et al. Prognostically relevant gene signatures of high-grade serous ovarian carcinoma. J Clin Invest. 2013;123:517–25.
Rieg T, Vallon V. Development of SGLT1 and SGLT2 inhibitors. Diabetologia. 2018;61:2079–86.
This work was supported by the National Natural Science Foundation of China (grant number: 82002747); the Chinese Society of Clinical Oncology (CSCO)-Rhoche oncologic research foundation (grant numbers: Y-2019Roche-077); Chinese Society of Clinical Oncology (CSCO)-BMS oncologic research foundation (grant numbers: Y-BMS2019–018); Pandeng Foundation of China National Cancer Center (grant numbers: NCC201809B032); Shanghai Sailing Program (grant number: 20YF1408000); and Shanghai Anticancer Association EYAS PROJECT (grant number: SACA-CY19A07).
Ethics approval and consent to participate
The research program was approved by the Ethics Committee of Fudan University Shanghai Cancer Center and with consent from all patients. All procedures were performed in accordance with the Declaration of Helsinki and relevant policies in China.
Consent for publication
The authors declare that there are no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, S., Wu, Y., Wang, S. et al. A risk model of gene signatures for predicting platinum response and survival in ovarian cancer. J Ovarian Res 15, 39 (2022). https://doi.org/10.1186/s13048-022-00969-3
- Ovarian cancer
- Platinum response
- Prognostic model