101
Views
14
CrossRef citations to date
0
Altmetric
Original Research

Identification of key candidate genes and small molecule drugs in cervical cancer by bioinformatics strategy

, , , , , & show all
Pages 3533-3549 | Published online: 13 Sep 2018

Abstract

Purpose

Cervical cancer (CC) is one of the most common malignant tumors among women. The present study aimed at integrating two expression profile datasets to identify critical genes and potential drugs in CC.

Materials and methods

Expression profiles, GSE7803 and GSE9750, were integrated using bioinformatics methods, including differentially expressed genes analysis, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and protein–protein interaction (PPI) network construction. Subsequently, survival analysis was performed among the key genes using Gene Expression Profiling Interactive Analysis websites. Connectivity Map (CMap) was used to query potential drugs for CC.

Results

A total of 145 upregulated genes and 135 downregulated genes in CC were identified. The functional changes of these differentially expressed genes related to CC were mainly associated with cell cycle, DNA replication, p53 signaling pathway, and oocyte meiosis. A PPI network was identified by STRING with 220 nodes and 2,111 edges. Thirteen key genes were identified as the intersecting genes of the enrichment pathways and the top 20 nodes in PPI network. Survival analysis revealed that high mRNA expression of MCM2, PCNA, and RFC4 was significantly associated with longer overall survival, and the survival was significantly better in the low-expression RRM2 group. Moreover, CMap predicted nine small molecules as possible adjuvant drugs to treat CC.

Conclusion

Our study found key dysregulated genes involved in CC and potential drugs to combat it, which might provide insights into CC pathogenesis and might shed light on potential CC treatments.

Introduction

Cervical cancer (CC) is the second most common malignant tumor among women, responsible for ~527,600 new cases and >265,700 deaths annually.Citation1 Despite advances in screening detection and new treatment strategies, CC is one of the leading causes of cancer death among females in many developing countries.Citation2,Citation3 Although most patients can be cured if diagnosed at an early stage, poor prognosis is observed with secondary metastatic cancer and tumor relapse.

Although human papillomavirus (HPV) is a prerequisite for CC, only a small number of women infected by this virus develop cancer. Thus, other risk factors should be considered as cofactors contributing to the progression of CC.Citation4 Dysregulated genes play important roles in CC development.Citation5 Several studies have used gene expression profiling to identify key genes between CC samples and normal cervix.Citation6Citation9 Hundreds of differentially expressed genes (DEGs) were detected. However, DEGs reported in different studies vary enormously with only some of them consistently detected. Therefore, the discovery of novel effective therapeutic targets against CC is urgently required.

A number of chemotherapeutic agents have shown activity against CC, including cisplatin,Citation10 bevacizumab,Citation11 carboplatin,Citation12 paclitaxel,Citation13 ifosfamide,Citation14 and topotecan.Citation15 Various combinations of these agents are recommended as therapies.Citation16 A recent systematic literature review found that carboplatin–paclitaxel is equally effective and less toxic than cisplatin–paclitaxel as the first-line therapy for metastatic CC.Citation17 However, patients overall survival (OS) times remains short, indicating an urgent need to discover some molecular drugs that are more efficient and selective. Based on bioinformatics approaches, several studies found small molecules as potential anticancer agents.Citation18Citation20

In this study, we selected the following microarray datasets GSE7803 and GSE9750 from the Gene Expression Omnibus (GEO) database to identify DEGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the identified DEGs was investigated. A protein–protein interaction (PPI) network was constructed to elucidate the significant relationships among DEGs and to identify key genes. Furthermore, the Kaplan–Meier estimator was used on the Gene Expression Profiling Interactive Analysis (GEPIA) website. Candidate small molecules were identified for their potential use in the treatment of CC.

Materials and methods

Data collection

Two CC microarray datasets were downloaded from the GEO website (http://www.ncbi.nlm.nih.gov/geo/). GSE7803 microarray data contained 21 CC tissues and 10 normal cervical epithelia tissues.Citation6 GSE9750 included 33 tumors samples and 24 healthy cervical samples.Citation7 Both the profile datasets were based on the Affymetrix GPL96 platform (Affymetrix Human Genome U133A Array). Because Connectivity Map (CMap) strictly required all probesets obtained from the Affymetrix Human Genome U133A Array,Citation21 we predicted the drugs for the DEGs measured only in this platform with high accuracy. GSE63514 data included 28 cancer cases and 24 normal casesCitation8 and were chosen to validate RRM2 mRNA expression in our analysis.

Data preprocessing and DEGs screening

The raw data were standardized and transformed into expression values using the affy package of Bioconductor (http://www.bioconductor.org/).Citation22 DEGs between cancer and normal samples were selected by significance analysis using the empirical Bayes methods within limma package.Citation23 False discovery rate (FDR) <0.05 and |log2 (fold change)| >1 were set as the cutoff criteria for the identification of DEGs. Common dysregulated probesets between GSE7803 and GSE9750 were selected for subsequent analyzes.

KEGG pathway analysis

Pathway enrichment analysis was performed using the clusterProfiler package and a pathway with an adjusted P-value <0.05 was considered significantly enriched.Citation24 DEGs that we identified could be involved in multiple pathways, Thus, some overlap was observed among the pathways. We identified the significant pathways that shared the same DEGs and used Cytoscape (version 3.5.1) to construct graphical representations of the interactive relationships among the pathways.Citation25

PPI network construction and analysis

The PPI pairs of the screened DEGs were analyzed using the online database STRING version 10.5 (https://string-db.org/).Citation26 The pairs with combined scores >0.4 were used for the PPI network construction, then the Cytoscape software was used to construct the network and analyze the interaction relationship of the candidate DEGs encoding proteins in CC.

Validation of key genes

Key genes were identified as the intersecting genes of the enrichment pathways and top 20 nodes in PPI network. To confirm the reliability of these genes from our detection, we analyzed their prognostic and expression in CC using GEPIA.Citation27 GEPIA is an interactive web application for gene expression analysis based on 9,736 tumors and 8,587 normal samples from the Cancer Genome Atlas (TCGA) and the Genotype-Tissue Expression databases.Citation28,Citation29 We evaluated the expression of key genes in CC tissues and normal tissues. Then the survival curve and boxplot were performed to visualize the relationships.

Identification of candidate small molecules

The CC gene signature was used to query CMap to find potential drugs for use in patients.Citation21 CMap is an in silico method to predict potential drugs that could possibly reverse, or induce, the biological state encoded in particular gene expression signatures. The common differently expressed probesets in GSE7803 and GSE9750 between CC samples and healthy controls were divided into upregulated and downregulated groups. Then, these probesets were used to query the CMap database. Finally, the enrichment score representing similarity was calculated, ranging from −1 to 1. A positive connectivity score indicates that a drug is able to induce the input signature in human cell lines. Conversely, a negative connectivity score indicates that a drug is able to reverse the input signature. Negative connectivity scores were investigated, which indicate potential therapeutic value. After rank ordering all instances, the connectivity score of various instances were filter by the number of instances (N>10) and P-value (<0.05).

Results

DEGs identification

The two mRNA expression profiles, including 54 patients with CC and 34 healthy individuals, were included in our study. Using a FDR <0.05 and |logFC| >1 as cutoff criteria, we extracted 443 and 848 differentially expressed probesets from the expression profile datasets GSE7803 and GSE9750, respectively. In GSE7803, 212 unregulated probes and 231 downregulated probes were identified. A total of 376 unregulated probes and 472 downregulated probes were identified in GSE9750. After being overlapped, the common 149 upregulated and 146 downregulated probesets corresponding to 145 upregulated and 135 downregulated genes were identified from the two profile datasets ().

CC significant pathways evaluation

A total of 16 pathways with adjusted P-value <0.05 were found enriched including 10 upregulated and 6 downregulated pathways (). The most significant upregulated pathway was cell cycle; the other significant pathways included DNA replication, oocyte meiosis, p53 signaling pathway, microRNAs in cancer, and cellular senescence. The downregulated pathways included arachidonic acid metabolism, serotonergic synapse, gap junction, estrogen signaling pathway, signaling pathways regulating pluripotency of stem cells, and proteoglycans in cancer (). In order to consider the potentially biological complexities in which a gene may belong to multiple pathways and provide information of numeric changes, we constructed pathway–gene networks to extract the complex association (). Cell cycle pathway contained the most significant genes in the network.

Table 1 Pathway enrichment analysis of DEGs function in CC

Figure 1 Significantly enriched pathway terms associated to DEGs in CC.

Notes: (A) KEGG pathways in CC DEGs enrichment analysis. (B) Upregulated pathway–gene network including 35 upregulated genes and 10 pathways. (C) Downregulated pathway–gene network including 26 downregulated genes and 6 pathways.

Abbreviations: CC, cervical cancer; DEGs, differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes; HTLV-I, human T-lymphotropic virus type 1.

Figure 1 Significantly enriched pathway terms associated to DEGs in CC.Notes: (A) KEGG pathways in CC DEGs enrichment analysis. (B) Upregulated pathway–gene network including 35 upregulated genes and 10 pathways. (C) Downregulated pathway–gene network including 26 downregulated genes and 6 pathways.Abbreviations: CC, cervical cancer; DEGs, differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes; HTLV-I, human T-lymphotropic virus type 1.

PPI network construction

STRING was used for mining proteins expressed by DEGs which can interact with others. At a combined score >0.4, a total of 222 DEGs (118 upregulated and 104 downregulated genes) among the 280 commonly altered DEGs were filtered into the DEGs PPI network, containing 222 nodes and 2,111 edges (). NetworkAnalyzer app in Cyto-scape was used to calculate the node degree.Citation25 The genes CDK1, PCNA, TOP2A, CCNB1, RFC4, MAD2L1, NDC80, CCNB2, AURKA, TYMS, MCM2, FEN1, RRM2, NCAPG, TTK, PRC1, MCM4, ZWINT, DTL, and MCM6 were the most significant 20 node degree genes and were selected as the hub nodes, since they might play important roles in CC progression ().

Figure 2 PPI network analysis.

Notes: (A) Using the STRING online database, a total of 222 DEGs (118 upregulated in red standing for upregulation and 104 downregulated genes in green standing for downregulation) were filtered into the DEGs PPI network. Bigger nodes represent genes with more links. (B) Degree of the top 20 nodes in the PPI network. All these nodes are upregulated genes.

Abbreviations: DEGS, differentially expressed genes; PPI, protein–protein interaction.

Figure 2 PPI network analysis.Notes: (A) Using the STRING online database, a total of 222 DEGs (118 upregulated in red standing for upregulation and 104 downregulated genes in green standing for downregulation) were filtered into the DEGs PPI network. Bigger nodes represent genes with more links. (B) Degree of the top 20 nodes in the PPI network. All these nodes are upregulated genes.Abbreviations: DEGS, differentially expressed genes; PPI, protein–protein interaction.

Key gene signatures identification in CC

Compared with KEGG enrichment genes, 13 of the top 20 nodes in the PPI network, including AURKA, CCNB1, CCNB2, CDK1, FEN1, MAD2L1, MCM2, MCM4, MCM6, PCNA, RFC4, RRM2, and TTK were found as key genes. Further survival analyses on these key genes were employed to evaluate their effects on CC patients’ survival using GEPIA. Expression levels of MCM2, PCNA, RFC4, and RRM2 were significantly related to the OS of patients with cervical squamous cancer (P<0.05). High expression of MCM2, PCNA, and RFC4 could result in a high survival rate, and increased RRM2 expression in CC was significantly associated with shorter patients’ survival (). The expression of these four genes was significantly higher in CC tissues compared to that of normal tissues (P<0.01; ). Together, the high level of these four genes might represent the important prognostic factor to predict the survival of CC. GSE63514 was used to validate RRM2 mRNA expression. The results showed that RRM2 expression was significantly higher in CC compared to that of normal tissues (P<0.01; ). The PPI network based on RRM2 found that PCNA and RFC4 have a close relationship with RRM2, and most of the proteins in the network were related to cell cycle ().

Figure 3 Survival curves and expression boxplots of key genes using GEPIA website.

Notes: (A–D) Expression level of MCM2, PCNA, RFC4, and RRM2 was significantly related to the overall survival of patients with cervical squamous cancer (P<0.05). (E–H) MCM2, PCNA, RFC4, and RRM2 were significantly upregulated in cervical squamous cancer compared with normal tissues (P<0.01).

Abbreviations: CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; GEPIA, Gene Expression Profiling Interactive Analysis; TPM, transcripts per million.

Figure 3 Survival curves and expression boxplots of key genes using GEPIA website.Notes: (A–D) Expression level of MCM2, PCNA, RFC4, and RRM2 was significantly related to the overall survival of patients with cervical squamous cancer (P<0.05). (E–H) MCM2, PCNA, RFC4, and RRM2 were significantly upregulated in cervical squamous cancer compared with normal tissues (P<0.01).Abbreviations: CESC, cervical squamous cell carcinoma and endocervical adenocarcinoma; GEPIA, Gene Expression Profiling Interactive Analysis; TPM, transcripts per million.

Figure 4 RRM2 validation using GSE63514 and PPI network.

Notes: (A) GSE63514 showed higher expression of RRM2 in CC tissues compared with normal cervical tissues (P<0.01). (B) RRM2 PPI network based on STRING.

Abbreviations: CC, cervical cancer; PPI, protein–protein interaction.

Figure 4 RRM2 validation using GSE63514 and PPI network.Notes: (A) GSE63514 showed higher expression of RRM2 in CC tissues compared with normal cervical tissues (P<0.01). (B) RRM2 PPI network based on STRING.Abbreviations: CC, cervical cancer; PPI, protein–protein interaction.

Related small molecule drugs screening

In order to screen out small molecule drugs, consistent differently expressed probesets between CC samples and healthy controls were analyzed with CMap. The related small molecules with highly significant correlations are listed in . Among these molecules, trichostatin A (TSA), tanespimycin, vorinostat, trifluoperazine, prochlorperazine, and thioridazine showed higher negative correlation and the potential to treat CC.

Table 2 Results of CMap analysis

Discussion

Driver genes play vital roles during stages of cancer progression. Although many studies on CC development are available, more efforts are needed to identify driver genes and candidate drugs that may shed light on CC treatments. This study integrated two gene profile datasets based on Affymetrix Human Genome U133A Array, utilized bioinformatics methods to analyze these datasets, and identified 280 commonly changed DEGs (145 upregulated and 135 downregulated). Pathway enrichment analysis indicated that cell cycle, DNA replication, oocyte meiosis, p53 signaling pathway, cellular senescence, and DNA repair-relevant biological pathways were overrepresented among the upregulated genes. The PPI network was constructed including 222 nodes/DEGs and 2,111 edges. Thirteen key genes were identified and chosen for survival analysis. MCM2, PCNA, RFC4, and RRM2 were clearly related to the prognosis of patients. In addition, small molecules that can provide new insights in CC therapeutic studies were identified.

Many researchers have found that four key genes were involved in cell cycle, participating in tumorigenesis and tumor proliferation. MCM2 has been studied in a wide range of human malignancies and is associated with tumor histopathological grade in several malignancies, including colon, oral cavity, ovarian, urothelial, and non-small cell lung carcinoma.Citation30Citation34 In cervical carcinoma and precancerous lesions, MCM2 is overexpressed and positively correlated with high risk types of HPV.Citation35 Amaro Filho et al also reported an increasing expression of MCM2 in invasive CC compared to control, but they suggested that MCM2 is not a good biomarker when comparing the different clinical stages of CC.Citation36 PCNA acts as a central coordinator of DNA transactions by providing a multivalent interaction surface for factors involved in DNA replication and cell cycle regulation. Owing to its function, PCNA has been widely used as a tumor marker for cancer cell progression and patient prognosis.Citation37Citation39 A recent systematic literature review found that the expression of PCNA is significantly associated with poor 5-year survival, International Federation of Gynecology and Obstetrics stage, or WHO grade, suggesting its use as a valuable prognostic and diagnostic biomarker in CC and gliomas.Citation40 RFC4 is involved in cancer. Knockdown of RFC4 in HepG2 cells induces apoptosis.Citation41 Similar results were discovered in breast carcinoma.Citation42 In colorectal cancer, overexpression of RFC4 is associated with tumor progression and poor survival outcome.Citation43 Additionally, with gene network reconstruction, RFC4 is regarded as one of the main drivers in cell cycle network in CC.Citation44 Together with our results, MCM2, PCNA, and RFC4 were significantly upregulated in CC compared with normal samples, and in CC patients, the survival rate was positively correlated with the high expression of these genes.

RRM2 is markedly upregulated in many patients’ cancer types and indeed acts as an oncogene.Citation45 RRM2 knockdown reduces cell proliferation and invasive ability in gastric cancer and pancreatic adenocarcinoma.Citation46,Citation47 Wang et al reported that RRM2 expression inhibition significantly increases apoptosis, promotes cell cycle arrest at the G1 phase, and inhibits tumor formation in CC nude mice transplant models.Citation48 Several studies showed that RRM2 is an independent prognostic factor and may predict poor survival in ovarian cancer, bladder cancer, breast cancer, and CC.Citation49Citation52 In this study, according to the PPI network, RRM2 closely interacts with PCNA and RFC4 involved in CC progression. Therefore, a further exploration of cell cycle and related genes was of enormous significance.

Consistent with our results, recent studies have also reported the identification of DEGs in CC. van Dam et al used three publicly available Affymetrix gene expression datasets (GSE5787, GSE7803, and GSE9750) and identified five cancer hallmarks enriched pathways in CC, showing that cell cycle deregulation is the major component of CC biology. They also identified seven probesets that were highly expressed in both CIN3 samples compared to normal samples and in cancer samples compared to CIN3 samples. From these probesets, six genes (AURKA, DTL, HMGB3, KIF2C, NEK2, and RFC4) were overexpressed in CC cell lines compared to cancer samples, suggesting their potential role as biomarkers in CC early diagnosis.Citation53 One of these genes, such as RFC4, was also identified in our study. Furthermore, our conclusion generated from both expression and survival analysis suggested that RFC4 might have a prognostic value. Another report from Li et al was based on TCGA data.Citation54 They found that MCM2, MCM4, MCM5, PCNA, and RNASEH2A participating in DNA replication pathway might be prognostic biomarkers in CC patients. MCM2 and PCNA were also found in our results.

Several small molecules with potential therapeutic efficacy against CC were identified. The most significant small molecules in our result have been reported to display anticancer activity. TSA, as a histone deacetylase (HDAC) inhibitor, shows a potential therapeutic effect in various types of cancer cells, when combined with radiotherapy or chemotherapy.Citation55,Citation56 In particular, TSA and its hydroxamate analogs can effectively and selectively induce tumor growth arrest at very low concentrations.Citation57 Additionally, TSA can inhibit HeLa cells growth via Bcl-2-mediated and caspase-dependent apoptosis.Citation58 Vorinostat is a hydroxamate-based pan-HDAC inhibitor also known as suberoylanilide hydroxamic acid used for the treatment of cutaneous T-cell lymphoma.Citation59 In HeLa cell, both mRNA and protein levels of HPV18 E6 and E7 were reduced after vorinostat treatment.Citation60 Furthermore, vorinostat promotes SiHa apoptosis through upregulation of p21 and Bax mRNA and protein, leading to cell cycle arrest in G0/G1 phase.Citation61 Thioridazine, a derivative of phenothiazine, displays anticancer abilities in a variety of cancer types and can reverse multidrug resistance.Citation62Citation64 Kang et al found that thioridazine can inhibit the PI3K/Akt/mTOR/p70S6K signaling pathway and exert cytotoxic effect on CC cells by inducing cell cycle arrest and apoptosis.Citation65 Thus, we might suppose that these identified drugs could play certain roles to combat CC.

Conclusion

Using bioinformatics analysis, 280 DEGs were identified, which were significantly enriched in several pathways, mainly associated with cell cycle, DNA replication, oocyte meiosis, p53 signaling pathway, and cellular senescence. We also identified key genes including MCM2, PCNA, RFC4, and RRM2 that might play important roles in CC and that might represent novel biomarkers in CC diagnosis, prognosis, and therapy. Additionally, a group of small molecules was identified that might be exploited as adjuvant drugs for improved therapeutics for CC. However, further investigations are required to validate the predicted drugs.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (grant no 81360336), the Joint Special Funds for the Department of Science and Technology of Yunnan Province – Kunming Medical University (grant no 2015FB017) and One Hundred Young and Middle - Aged Academic and Technical Backbone Project of Kunming Medical University (grant no. 60117190449).

Supplementary material

Table S1 Common dysregulated probes identified in GSE7803 and GSE9750

Disclosure

The authors report no conflicts of interest in this work.

References

  • TorreLABrayFSiegelRLFerlayJLortet-TieulentJJemalAGlobal cancer statistics, 2012CA Cancer J Clin20156528710825651787
  • ForouzanfarMHForemanKJDelossantosAMBreast and cervical cancer in 187 countries between 1980 and 2010: a systematic analysisLancet201137898011461148421924486
  • MaguireRKotronoulasGSimpsonMPatersonCA systematic review of the supportive care needs of women living with and beyond cervical cancerGynecol Oncol2015136347849025462200
  • zur HausenHHuman papillomaviruses in the pathogenesis of anogenital cancerVirology199118419131651607
  • OjesinaAILichtensteinLFreemanSSLandscape of genomic alterations in cervical carcinomasNature2014506748837137524390348
  • ZhaiYKuickRNanBGene expression analysis of prein-vasive and invasive cervical squamous cell carcinomas identifies HOXC10 as a key mediator of invasionCancer Res20076721101631017217974957
  • ScottoLNarayanGNandulaSVIdentification of copy number gain and overexpressed genes on chromosome arm 20q by an integrative genomic approach in cervical cancer: potential role in progressionGenes Chromosomes Cancer200847975576518506748
  • den BoonJAPyeonDWangSSMolecular transitions from papillomavirus infection to cervical precancer and cancer: role of stromal estrogen receptor signalingProc Natl Acad Sci U S A201511225E3255E326426056290
  • PappaKIPolyzosAJacob-HirschJProfiling of discrete gyne-cological cancers reveals novel transcriptional modules and common features shared by other cancer types and embryonic stem cellsPLoS One20151011e014222926559525
  • LeischingGRLoosBBothaMHEngelbrechtAMThe role of mTOR during cisplatin treatment in an in vitro and ex vivo model of cervical cancerToxicology2015335727826201060
  • TewariKSSillMWLongHJImproved survival with bevacizumab in advanced cervical cancerN Engl J Med2014370873474324552320
  • KatanyooKTangjitgamolSChongthanakornMTreatment outcomes of concurrent weekly carboplatin with radiation therapy in locally advanced cervical cancer patientsGynecol Oncol2011123357157621955483
  • WangXShenYZhaoYAdjuvant intensity-modulated radiotherapy (IMRT) with concurrent paclitaxel and cisplatin in cervical cancer patients with high risk factors: a phase II trialEur J Surg Oncol20154181082108826033265
  • DownsLSChuraJCArgentaPAIfosfamide, paclitaxel, and carboplatin, a novel triplet regimen for advanced, recurrent, or persistent carcinoma of the cervix: a phase II trialGynecol Oncol2011120226526921145100
  • MuderspachLIBlessingJALevenbackCMooreJLJrA phase II study of topotecan in patients with squamous cell carcinoma of the cervix: a gynecologic oncology group studyGynecol Oncol200181221321511354055
  • EskanderRNTewariKSChemotherapy in the treatment of metastatic, persistent, and recurrent cervical cancerCurr Opin Obstet Gynecol201426431432124979076
  • LorussoDPetrelliFCoinuARaspagliesiFBarniSA systematic review comparing cisplatin and carboplatin plus paclitaxel-based chemotherapy for recurrent or metastatic cervical cancerGynecol Oncol2014133111712324486604
  • YehCTWuATChangPMTrifluoperazine, an antipsychotic agent, inhibits cancer stem cell growth and overcomes drug resistance of lung cancerAm J Respir Crit Care Med2012186111180118823024022
  • ChenMHLinKJYangWLGene expression-based chemical genomics identifies heat-shock protein 90 inhibitors as potential therapeutic drugs in cholangiocarcinomaCancer2013119229330322810956
  • HassaneDCGuzmanMLCorbettCDiscovery of agents that eradicate leukemia stem cells using an in silico screen of public gene expression dataBlood2008111125654566218305216
  • LambJCrawfordEDPeckDThe Connectivity Map: using gene-expression signatures to connect small molecules, genes, and diseaseScience200631357951929193517008526
  • GautierLCopeLBolstadBMIrizarryRAaffy – analysis of Affyme-trix GeneChip data at the probe levelBioinformatics200420330731514960456
  • RitchieMEPhipsonBWuDlimma powers differential expression analyses for RNA-sequencing and microarray studiesNucleic Acids Res2015437e4725605792
  • YuGWangLGHanYHeQYclusterProfiler: an R package for comparing biological themes among gene clustersOMICS201216528428722455463
  • ShannonPMarkielAOzierOCytoscape: a software environment for integrated models of biomolecular interaction networksGenome Res200313112498250414597658
  • SzklarczykDMorrisJHCookHThe STRING database in 2017: quality-controlled protein–protein association networks, made broadly accessibleNucleic Acids Res201745D1D362D36827924014
  • TangZLiCKangBGaoGLiCZhangZGEPIA: a web server for cancer and normal gene expression profiling and interactive analysesNucleic Acids Res201745W1W98W10228407145
  • Cancer Genome Atlas Research NetworkWeinsteinJNCollissonEAMillsGBThe Cancer Genome Atlas Pan-Cancer analysis projectNat Genet201345101113112024071849
  • GTEx ConsortiumThe Genotype-Tissue Expression (GTEx) projectNat Genet201345658058523715323
  • WangYLiYZhangWYmRNA expression of minichromosome maintenance 2 in colonic adenoma and adenocarcinomaEur J Cancer Prev2009181404519077563
  • SzelachowskaJDziegielPJelen-KrzeszewskaJCorrelation of metallothionein expression with clinical progression of cancer in the oral cavityAnticancer Res200929258959519331208
  • GakiopoulouHKorkolopoulouPLevidouGMinichromosome maintenance proteins 2 and 5 in non-benign epithelial ovarian tumours: relationship with cell cycle regulators and prognostic implicationsBr J Cancer20079781124113417940502
  • BurgerMDenzingerSHartmannAWielandWFStoehrRObermannECMcm2 predicts recurrence hazard in stage Ta/T1 bladder cancer more accurately than CK20, Ki67 and histological gradeBr J Cancer200796111711171517505513
  • ZhangXTengYYangFMCM2 is a therapeutic target of lovastatin in human non-small cell lung carcinomasOncol Rep20153352599260525738322
  • ZhengJDiagnostic value of MCM2 immunocytochemical staining in cervical lesions and its relationship with HPV infectionInt J Clin Exp Pathol20158187588025755789
  • Amaro FilhoSMNuovoGJCunhaCBCorrelation of MCM2 detection with stage and virology of cervical cancerInt J Biol Markers2014294e363e37124706378
  • TachibanaKEGonzalezMAColemanNCell-cycle-dependent regulation of DNA replication and its relevance to cancer pathologyJ Pathol2005205212312915643673
  • ZhaoHLoYHMaLTargeting tyrosine phosphorylation of PCNA inhibits prostate cancer growthMol Cancer Ther2011101293621220489
  • WangLFChaiCYKuoWRTaiCFLeeKWHoKYCorrelation between proliferating cell nuclear antigen and p53 protein expression and 5-year survival rate in nasopharyngeal carcinomaAm J Otolaryngol200627210110516500472
  • LvQZhangJYiYProliferating cell nuclear antigen has an association with prognosis and risks factors of cancer patients: a systematic reviewMol Neurobiol20165396209621726558632
  • AraiMKondohNImazekiNThe knockdown of endogenous replication factor C4 decreases the growth and enhances the chemosensitivity of hepatocellular carcinoma cellsLiver Int2009291556218492021
  • SrihariSKalimuthoMLalSUnderstanding the functional impact of copy number alterations in breast cancer using a network modeling approachMol Biosyst201612396397226805938
  • XiangJFangLLuoYLevels of human replication factor C4, a clamp loader, correlate with tumor progression and predict the prognosis for colorectal cancerJ Transl Med20141232025407051
  • MineKLShulzhenkoNYambartsevAGene network reconstruction reveals cell cycle and antiviral genes as major drivers of cervical cancerNat Commun20134180623651994
  • XuXPageJLSurteesJABroad overexpression of ribonucleotide reductase genes in mice specifically induces lung neoplasmsCancer Res20086882652266018413732
  • KangWTongJHChanAWTargeting ribonucleotide reductase M2 subunit by small interfering RNA exerts anti-oncogenic effects in gastric adenocarcinomaOncol Rep20143162579258624756820
  • DuxburyMSWhangEERRM2 induces NF-kappaB-dependent MMP-9 activation and enhances cellular invasivenessBiochem Biophys Res Commun2007354119019617222798
  • WangNLiYZhouJDownregulation of ribonucleotide reductase subunits M2 induces apoptosis and G1 arrest of cervical cancer cellsOncol Lett20181533719372529556274
  • FerrandinaGMeyVNannizziSExpression of nucleoside transporters, deoxycitidine kinase, ribonucleotide reductase regulatory subunits, and gemcitabine catabolic enzymes in primary ovarian cancerCancer Chemother Pharmacol201065467968619639316
  • MorikawaTMaedaDKumeHHommaYFukayamaMRibonucleotide reductase M2 subunit is a novel diagnostic marker and a potential therapeutic target in bladder cancerHistopathology201057688589221166702
  • KretschmerCSterner-KockASiedentopfFSchoeneggWSchlagPMKemmnerWIdentification of early molecular markers for breast cancerMol Cancer20111011521314937
  • SuYFWuTFKoJLThe expression of ribonucleotide reductase M2 in the carcinogenesis of uterine cervix and its relationship with clinicopathological characteristics and prognosis of cancer patientsPLoS One201493e9164424637958
  • van DamPAvan DamPJRolfoCIn silico pathway analysis in cervical carcinoma reveals potential new targets for treatmentOncotarget2016732780279526701206
  • LiXTianRGaoHIdentification of significant gene signatures and prognostic biomarkers for patients with cervical cancer by integrated bioinformatic methodsTechnol Cancer Res Treat201817112
  • RanganathanPRangnekarVMExploiting the TSA connections to overcome apoptosis-resistanceCancer Biol Ther20054439139215908800
  • HajjiNWallenborgKVlachosPNymanUHermansonOJosephBCombinatorial action of the HDAC inhibitor trichostatin A and etoposide induces caspase-mediated AIF-dependent apoptotic cell death in non-small cell lung carcinoma cellsOncogene200827223134314418071312
  • VanhaeckeTPapeleuPElautGRogiersVTrichostatin A-like hydroxa-mate histone deacetylase inhibitors as therapeutic agents: toxicological point of viewCurr Med Chem200411121629164315180568
  • YouBRParkWHTrichostatin A induces apoptotic cell death of HeLa cells in a Bcl-2 and oxidative stress-dependent mannerInt J Oncol201342135936623165748
  • DuvicMVuJVorinostat: a new oral histone deacetylase inhibitor approved for cutaneous T-cell lymphomaExpert Opin Investig Drugs200716711111120
  • HeHLiuXWangDSAHA inhibits the transcription initiation of HPV18 E6/E7 genes in HeLa cervical cancer cellsGene201455329810425300249
  • XingJWangHXuSHanPXinMZhouJLSensitization of suberoylanilide hydroxamic acid (SAHA) on chemoradiation for human cervical cancer cells and its mechanismEur J Gynaecol Oncol201536211712226050346
  • Gil-AdIShtaifBLevkovitzYPhenothiazines induce apoptosis in a B16 mouse melanoma cell line and attenuate in vivo melanoma tumor growthOncol Rep200615110711216328041
  • LiJYaoQYXueJSDopamine D2 receptor antagonist sulpiride enhances dexamethasone responses in the treatment of drug-resistant and metastatic breast cancerActa Pharmacol Sin20173891282129628649130
  • SeoSUChoHKMinKJThioridazine enhances sensitivity to carboplatin in human head and neck cancer cells through downregulation of c-FLIP and Mcl-1 expressionCell Death Dis201782e259928182008
  • KangSDongSMKimBRThioridazine induces apoptosis by targeting the PI3K/Akt/mTOR pathway in cervical and endometrial cancer cellsApoptosis201217998999722460505