- Open Access
Bioinformatics analysis of mRNA and miRNA microarray to identify the key miRNA-mRNA pairs in cisplatin-resistant ovarian cancer
BMC Cancer volume 21, Article number: 452 (2021)
Ovarian cancer (OC) is a gynecological malignancy with the highest mortality rate. Cisplatin (DDP) based chemotherapy is a standard strategy for ovarian cancer. Despite good response rates for initial chemotherapy, almost 80% of the patients treated with DDP based chemotherapy will experience recurrence due to drug-resistant, which will ultimately result in fatality. The aim of the present study was to examine the pathogenesis and potential molecular markers of cisplatin-resistant OC by studying the differential expression of mRNAs and miRNAs between cisplatin resistant OC cell lines and normal cell lines.
Two mRNA datasets (GSE58470 and GSE45553) and two miRNA sequence datasets (GSE58469 and GSE148251) were downloaded from the Gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) were screened by the NetworkAnalyst. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to analyze the biological functions of DEGs. The protein-protein interaction network was constructed using STRING and Cytoscape software to identify the molecular mechanisms of key signaling pathways and cellular activities. FunRich and MiRNATip databases were used to identify the target genes of the DEMs.
A total of 380 DEGs, and 5 DEMs were identified. Protein–protein interaction (PPI) network of DEGs containing 379 nodes and 1049 edges was constructed, and 4 key modules and 24 hub genes related to cisplatin-resistant OC were screened. Two hundred ninety-nine target genes of the 5 DEMs were found out. Subsequently, one of these 299 target genes (UBB) belonging to the hub genes of GSE58470 and GSE45553 was identified by MCODE and CytoHubba,which was regulated by one miRNA (mir-454).
One miRNA–mRNA regulatory pairs (mir-454-UBB) was established. Taken together, our study provided evidence concerning the alteration genes involved in cisplatin-resistant OC, which will help to unravel the mechanisms underlying drug resistant.
Ovarian cancer (OC), originated in gynecological genital tract, is the ninth most common cancer in female with 313,959 newly diagnosed cases and 207,252 new death in 2020 worldwide . OC is more frequently diagnosed at an advanced stage because of the lack of efficient screening measures [2, 3]. And the advanced stages prognosis (FIGO stage III and IV) is extremely poor with 5-year survival rates of approximately 39 and 17%, respectively . The standard treatment of OC involves cytoreductive surgery and platinum-based chemotherapy. Despite the good response rates to initial surgery and chemotherapy, for over a decade, the median progression-free survival rate of patients remains low at about 18 months . Indeed, almost 80% of the patients experience recurrence as a drug resistant population in tumors that ultimately results in fatality .
Cisplatin is the first generation of platinum-based drugs that can directly interact with DNA of cancer cells to prevent DNA synthesis and RNA transcription . Numerous protein-coding genes were revealed that are connected with ovarian cancer cisplatin resistance. For example, PTGER3 overexpression confers to cisplatin resistance in OC through up-regulation of Ras-MAPK/Erk-ETS1-ELK1/CFTR1 axis . CpG island promoter hypermethylation of BRCA1, PTK6, PRKCE can enhance sensitivity to cisplatin in OC [9, 10]. Genes with coding ability only account for 2% of all genes in human genome, while the rest belong to noncoding RNAs including microRNA (miRNA), long noncoding RNA (lncRNA), and so on . MicroRNAs (miRNAs) are endogenous, highly conserved small RNAs, 20–24 nucleotides in length, and specifically bind to target mRNA to inhibit post-transcriptional gene expression. MiRNA were also reported to play an important role in cisplatin resistance. MiR-375 overexpression can enhance the cisplatin sensitivity in OC by targeting PAX2 . MiR-138-5p downregulation promotes overexpression of EZH2 and SIRT1 in OC cell, thereby regulating the cisplatin sensitivity .
In the present study, bioinformatics tools were used to analyze the cisplatin resistant OC expression profile chips in a public gene chip database for the purpose of identifying differentially expressed genes (DEGs), differentially expressed miRNAs (DEMs), and of constructing miRNA–mRNA regulatory networks involved in cisplatin-resistant OC.
Materials and methods
We screened the gene chip data by using the Gene expression omnibus (GEO, http://www.ncbi.nlm.nih.gov/geo) database . Cisplatin resistant ovarian cancer public gene expression data sets (GSE58470 and GSE45553) and miRNA expression data sets (GSE58469 and GSE148251) were selected and downloaded from GEO, with the keywords “ovarian cancer” and “cisplatin resistant” [15,16,17]. The GSE58470 and GSE45553 datasets, which are two mRNA datasets, contain 2 ovarian cancer cell lines (IGROV-1, OVCAR-8) and 14 samples, including 7 cisplatin resistant cell samples and 7 parental cell samples. The miRNA GSE58469 and GSE148251 datasets, which were analyzed using MicroRNA expression beadchips, contain 2 OC cell lines (IGROV-1, W1) and 12 samples, including 6 cisplatin resistant cell samples and 6 parental cell samples.
Differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs)
The raw data of GSE58470, GSE45553 and GSE58469, GSE148251 were effectively processed using the NetworkAnalyst, an online data-analytics platform, using correction, normalization and log2 conversion . The DEGs in cisplatin resistant OC cell lines compared with parental cell lines were determined using limma algorithm . DEGs were screened with a false discovery rate (FDR) corrected p-value< 0.05 and |log2 fold-change (FC)| > 1. The DEMs in cisplatin resistant OC cell lines compared with parental cell lines were also processed by the NetworkAnalyst. The FDR corrected p-value< 0.05 and |log2 fold-change (FC)| > 1 were used as the screening thresholds.
Functional enrichment analysis of DEGs
DAVID (https://david.ncifcrf.gov/), a widely used web-based genomic functional annotation tool, was used for data annotation analysis . In our study, we used the DAVID to perform Gene Ontology (GO) analysis including cellular component (CC), molecular function (MF), biological process (BP), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. A p-value that is smaller than 0.05 and enrichment score above 1.5 were considered as significant enrichment.
Construction of protein-protein interaction (PPI) networks and module research
The PPI network of the DEGs was constructed and visualized using the STRING (https://string-db.org) database to determine the molecular mechanisms of key signaling pathways and cellular activities in cisplatin resistant ovarian cancer . An interaction score > 0.4 was considered to identify the significant PPIs. We then used the Cytoscape software (version 3.8.2; www. cytoscape.org) to analyze the PPI network . Relationships among DEGs were analyzed by NetworkAnalyzer plug-in of Cytoscape software to characterize small-world network through calculating the network properties such as the clustering coefficient of the network, distribution of node degree and the shortest path . Molecular Complex Detection (MCODE) was used to identify key clusters of genes within PPI network using the cutoff criteria (MCODE score > 5) [24, 25] with the default parameters (degree cutoff = 5, node score cutoff = 0.2, K-core = 2, and Max depth = 100). Hub genes in the network were selected using CytoHubba through connection degree method (cutoff criteria: degree above 20). Finally, we summarized the overlapping genes between results of MCODE and CytoHubba to create a consensus of predictions to identify more accurate hub genes.
miRNA target prediction
MiRNATip and FunRich are two bioinformatics platform and analysis tools for predicting DEMs target genes and miRNA-gene pairs [26,27,28]. In the present study, the targets of the DEMs were predicted using MiRNATip and FunRich (version 3.1.3). The Venny 2.1 online Tool (http://bioinfogp.cnb. csic.es) was used to find overlapping genes between DEGs and predictive targeted genes of DEMs . The miRNA-gene regulatory network was depicted and visualized using Cytoscape software.
Construction of miRNA–mRNA regulatory pairs
MiRNA–mRNA regulatory pairs related to cisplatin-resistance in ovarian cancer were constructed according to the miRNA targets prediction results and DEGs to display the interaction among miRNA and mRNA. Furthermore, the K-M plotter database (https://kmplot.com/analysis/) and OncomiR database (http://www.oncomir.org/) were used to assess the potential prognostic significance of our selected mRNAs and miRNAs by performing survival analysis [30, 31].
Identification of DEGs, DEMs
The data was successfully normalized to ensure the accuracy. Dataset GSE58470 and GSE45553 were utilized to identify DEGs in IGROV-1 and OVCAR-8 cell lines. As presented in Fig. 1a, a total of 380 overlapping DEGs were identified between IGROV-1 and OVCAR-8 cell lines. Volcano plots were drawn to give a direct presentation of all the examined genes in the data sets. The red dots represent the significantly upregulated genes, while the blue dots represent the significantly downregulated genes. Volcano plots for gene expression in IGROV-1 and OVCAR-8 cell lines were presented in Fig. 1b and c. In miRNA datasets GES58469 and GSE148251,5 overlapping DEMs were detected between IGROV-1 and W1 cell lines which was shown in Fig. 1d. Volcano plots for miRNA expression in IGROV-1 and W1 cell lines are presented in Fig. 1e and f.
Functional enrichment analysis of DEGs
GO analysis including cellular components (CC), molecular function (MF), biological process (BF), and KEGG analysis was performed using DAVID database to understand the functions of DEGs. The GO functional enrichment analysis resulted in a total of 380 DEGs mapped to 265 GO terms. With the FDR corrected p-value< 0.05 and enrichment score > 1.5 as the cut-off value, 17 significant enriched functional clusters were screened (Fig. 2a). In total, 7 GO terms were significantly enriched in cellular components including ‘extracellular exosome’, ‘cytosol’, ‘cell surface’, ‘extracellular space’, ‘cell-cell adherens junction’, ‘focal adhesion’ and ‘perinuclear region of cytoplasm’. Enrichment of 4 GO terms, such as ‘heparin binding’, ‘GTPase activity’, ‘cadherin binding involved in cell-cell adhesion’ and ‘protein homodimerization activity’, belongs to molecular functions. A total of 6 biological processes were enriched, mainly involving ‘type I interferon signaling pathway’, ‘response to estradiol’, ‘response to virus’, ‘cell migration’, ‘positive regulation of apoptotic process’ and ‘negative regulation of apoptotic process’. A total of 380 DEGs were mapped into the KEGG database using DAVID, enrichment score > 1.5 and p-value< 0.05 were used as an enrichment screening standard. In total, 27 enriched functional clusters of the DEGs were obtained (Fig. 2b), such as ‘Pathways in cancer’ (20 genes), ‘PI3K-Akt signaling pathway’ (18 genes), ‘Proteoglycans in cancer’ (15 genes), ‘MAPK signaling pathway’ 13 genes), ‘Regulation of actin cytoskeleton’ (13 genes).
Establishment of PPI network and identification of hub gene
The PPI network of 380 DEGs was constructed and visualized using STRING database. The isolated nodes and partially loosely connected gene nodes were removed, and the remaining DEGs together constituted a complex multi-center interaction network map, which contained 379 nodes and 1049 edges. The average node degree was 5.54 and the average local clustering coefficient was 0.394. Among the 379 nodes, top 20 DEGs with the highest degree of nodes were screened based on the Cytoscape software analysis results (Fig. 3). The results of the top 10 DEGs were as follows: TNF, CXCL8, CD44, UBB, CDK1, ITGB3, RAC1, RELA, OASL, VCL. The key clusters of genes were obtained using MCODE, with 9 key modules and a false degree cutoff = 5. Four significant key modules including 40 key genes with the MCODE score > 5 were identified (Fig. 4). Subsequently, the CytoHubba was used to find the hub genes in the PPI network of the DEGs. In total, 25 hub genes were identified. At last, we summarized the overlapping genes between the MCODE and CytoHubba results (Table 1). Twenty-four hub genes belonging to the GSE58470 and GSE45553 were identified. Furthermore, the functional and pathway enrichment of these genes was also analyzed using DAVID online tools shown in Tables 2 and 3. The 24 significantly upregulated or downregulated hub genes were also presented in heatmaps (Fig. 5a,b).
miRNA-gene regulatory network
The FunRich (http://www.funrich.org) online tool searched and predicted 611 target genes of the 5 DEMs and the MiRNATip predicted 1534 target genes of 5 DEMs. And then, there were 299 overlapping genes between the predicted targeted genes of FunRich and MiRNATip. As presented in Fig. 6, the 299 overlapping genes were regulated by 3 DEMs (mir-146a, mir-708, mir-454), but 2 DEMs (mir-675, mir-1250) have no targeted genes. In addition, the target gene UBB among the 299 overlapping genes also belonged to the 24 hub genes of GSE58470 and GSE45553 and was regulated by mir-454 (Fig. 7a). The expression of the UBB was low and the expression of the mir-454 was high in cisplatin resistant cell lines (Fig. 7b,c). Furthermore, to validate the predicting results of our present study, we have performed the survival analysis by using K-M plotter database and OncomiR database. The results of Kaplan-Meier survival analysis illustrated that lower UBB expression predicted poor progression-free-survival and overall survival both in ovarian cancers (Fig. 8a,b) and cisplatin treated ovarian cancers (Fig. 8c,d), respectively. The Kaplan-Meier survival analysis of mir-454 revealed that high expression of the mir-454 predicted poor prognosis in ovarian cancer (Fig. 8e). On the basis of these results, the miRNA–mRNA regulatory pair was established, indicating the importance of this miRNA-mRNA pair in cisplatin-resistant ovarian cancer.
Ovarian cancer is a lethal malignant cancer with different histopathological and biological characteristics, and a < 40% overall remission rate at all stages . Over the past a few decades, systematic chemotherapy has brought plenty of benefits to patients with OC . Cisplatin-based chemotherapy is the most commonly used treatment regimen; however, acquisition of cisplatin resistance is a major clinical obstacle to treatment of OC . Thence, there is a compulsive need to identify more effective biomarkers to overcome cisplatin resistance and improve the prognosis of OC. Noncoding RNAs are reported to have vital roles in regulating human cell behaviors [35, 36]. MicroRNAs (miRNA) are a kind of noncoding RNAs, which are composed of 22–25 nucleotides and bind with the 3′-untranslated region of targeted mRNAs to regulate mRNA expression . miRNAs have been regarded as essential regulators of chemoresistance of many cancers to common chemotherapy agents [38, 39], including cisplatin-based chemotherapy [40, 41].
In our study, 380 DEGs were screened from the GSE58470 and GSE45553 datasets and processed using bioinformatics methods. The results of the KEGG and GO enrichment analysis of the DEGs revealed that the identified genes were enriched in various signaling pathways, such as ‘Pathway in cancer’, ‘Proteoglycans in cancer’, ‘RIG-I-like receptor signaling pathway’, ‘Hippo signaling pathway’. RIG-I-like receptors, a member of cytosolic pattern-recognition receptors (PRRs), can distinguish pathogen-associated molecular patterns and activate antiviral immune responses . The RIG-I-like receptor signaling pathway induce the expression of large numbers of IFN-β and IFN-α proteins . A recent research has illustrated that the downregulation of RIG-I expression was correlated with poorer prognosis and weakened response to IFNα therapy in hepatocellular carcinoma (HCC) patients . The KEGG analysis results of our study revealed that the aberrant activation of the RIG-I like receptor signaling pathway may be related with the cisplatin resistance in the OC patients, which needed to be further validated. Hippo signaling pathway is one of the eight major signaling pathways commonly altered in human cancers . Dysfunction of Hippo signaling pathway has been implicated in an increasing number of human cancers . For example, miR372 and miR373 mediated silencing of LATS2 expression, a Hippo pathway tumor suppressor, is related with testicular germ cell tumors . MiR-149-5p can aggravate chemoresistance in ovarian cancer cells by directly targets the core kinase components of the Hippo signaling pathway .
In the present study, 5 overlapping DEMs (mir-146a, mir-708, mir-454, mir-675, mir-1250) were identified from the GSE58469 and GSE148251 datasets. Subsequently, we predicted the target genes of these 5 DEMs by using the FunRich and MiRNATip analysis tools and got 299 different target genes regulated by 3 DEMs (mir-146a, mir-708, mir-454). Among these 299 target genes, UBB also belongs to the 24 hub genes of GSE58470 and GSE45553 which is regulated by mir-454. Stefanie et al. has demonstrated that mir-454 can regulate UBB expression in kidney tissue and HEK293 cell lines . But, the role of mir-454-UBB regulatory pair has not been revealed in cisplatin treated ovarian cancer before. So, we perform our bioinformatics analysis and the results of our present study indicate that the mir-454-UBB regulatory pair is associated with cisplatin resistance in OC. To validate the results of our study, the survival analysis was performed and illustrated that UBB and mir-454 expressions are associated with the prognosis of ovarian cancer. The K-M survival plots proven that low expression of UBB and high expression of mir-454 may not only predict poor prognosis of OC (Fig. 8a,b,e) but also poor prognosis of cisplatin treated OC (Fig. 8c,d), which indicated that mir-454 and UBB regulatory pair may be correlated with cisplatin resistance in OC.
MiR-454 has been reported to be implicated in the progression of many types of cancer and play important roles in chemotherapeutic drug resistance. Several studies show that miR-454 functions as an oncogene in colorectal cancer , hepatocellular carcinoma , non-small cell lung cancer  and induce the oxaliplatin resistance in gastric carcinoma cells by targeting CYLD . UBB also known as ubiquitin B is highly expressed in all eukaryotic cells and can mark some target proteins for ubiquitin-proteasome system degradation. Previous studies emphasize UBB and UBB dependent ubiquitin-proteasomal protein degradation are essential in histone deacetylase inhibitor-induced tumor selectivity. Tian et al. studied the role of UBB expression in cervical cancer and demonstrated that UBB can maintain cancer stem-like characteristics . The results of our current study showed that miR-454 -UBB regulatory pair was significant in cisplatin resistance OC cell lines. This study may advance the understanding of the mechanism of cisplatin resistance in ovarian cancer and suggest that miR-454 and UBB may be two novel biomarker and therapeutic targets for ovarian cancer patients.
There are some limitations in our study. First, the regulatory network we constructed using bioinformatic methods was not validated by experimental work. In our future work, we will validate these findings using experiment manners on cell lines and human tissues. Second, all the data analyzed in our study was retrieved from one online database, which might result in some biases. So, further studies consist of larger sample sizes are needed to validate our findings.
Our study analyzed gene and miRNA expression between cisplatin resistant OC cell lines and parental OC cell lines using mRNA data and non-coding RNA data from the GEO database, and identified aberrant expression of mRNAs and miRNAs in cisplatin resistant OC cell lines. Based on these analysis results, we found that mir-454-UBB regulatory pair is significant in cisplatin resistance ovarian cancer cell lines. These finds will no doubt help us to understand the mechanisms under the skin of ovarian cancer cisplatin-resistance.
Availability of data and materials
All data generated or analyzed during this study are included in Gene expression omnibus GEO) public database. The mRNA datasets (GSE58470, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE58470 and GSE45553, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE45553). The miRNA datasets (GSE58469, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE58469 and GSE148251, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE148251).
Gene expression omnibus
Differentially expressed genes
Differentially expressed miRNAs
Kyoto Encyclopedia of Genes and Genomes
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021. https://doi.org/10.3322/caac.21660.
Jayson GC, Kohn EC, Kitchener HC, Ledermann JA. Ovarian cancer. Lancet (London, England). 2014;384(9951):1376–88.
Vaughan S, Coward JI, Bast RC Jr, Berchuck A, Berek JS, Brenton JD, et al. Rethinking ovarian cancer: recommendations for improving outcomes. Nat Rev Cancer. 2011;11(10):719–25. https://doi.org/10.1038/nrc3144.
Goh J, Mohan GR, Ladwa R, Ananda S, Cohen PA, Baron-Hay S. Frontline treatment of epithelial ovarian cancer. Asia Pac J Clin Oncol. 2015;11(S6):1–16.
Binju M, Padilla MA, Singomat T, Kaur P, Suryo Rahmanto Y, Cohen PA, et al. Mechanisms underlying acquired platinum resistance in high grade serous ovarian cancer - a mini review. Biochim Biophys Acta, Gen Subj. 2019;1863(2):371–8. https://doi.org/10.1016/j.bbagen.2018.11.005.
Grisham RN, Hyman DM, Iyer G. Targeted therapies for treatment of recurrent ovarian cancer. Clin Adv Hematol Oncol. 2014;12(3):158–62.
Yamashita K, Tsunoda S, Gunji S, Murakami T, Suzuki T, Tabata Y, et al. Intraperitoneal chemotherapy for peritoneal metastases using sustained release formula of cisplatin-incorporated gelatin hydrogel granules. Surg Today. 2019;49(9):785–94. https://doi.org/10.1007/s00595-019-01792-y.
Rodriguez-Aguayo C, Bayraktar E, Ivan C, Aslan B, Mai J, He G, et al. PTGER3 induces ovary tumorigenesis and confers resistance to cisplatin therapy through up-regulation Ras-MAPK/Erk-ETS1-ELK1/CFTR1 axis. EBioMedicine. 2019;40:290–304. https://doi.org/10.1016/j.ebiom.2018.11.045.
Stefansson OA, Villanueva A, Vidal A, Martí L, Esteller M. BRCA1 epigenetic inactivation predicts sensitivity to platinum-based chemotherapy in breast and ovarian cancer. Epigenetics. 2012;7(11):1225–9. https://doi.org/10.4161/epi.22561.
Yu W, Jin C, Lou X, Han X, Li L, He Y, et al. Global analysis of DNA methylation by methyl-capture sequencing reveals epigenetic control of cisplatin resistance in ovarian cancer cell. PLoS One. 2011;6(12):e29450. https://doi.org/10.1371/journal.pone.0029450.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52. https://doi.org/10.1038/nature12986.
Yang S, Yang R, Lin R, Si L. MicroRNA-375 inhibits the growth, drug sensitivity and metastasis of human ovarian cancer cells by targeting PAX2. J BUON. 2019;24(6):2341–6.
Zhang Y, Ai H, Fan X, Chen S, Wang Y, Liu L. Knockdown of long non-coding RNA HOTAIR reverses cisplatin resistance of ovarian cancer cells through inhibiting miR-138-5p-regulated EZH2 and SIRT1. Biol Res. 2020;53(1):18. https://doi.org/10.1186/s40659-020-00286-3.
Clough E, Barrett T. The Gene Expression Omnibus Database. Methods Mol Biol (Clifton, NJ). 2016;1418:93–110.
Arrighetti N, Cossa G, De Cecco L, Stucchi S, Carenini N, Corna E, et al. PKC-alpha modulation by miR-483-3p in platinum-resistant ovarian carcinoma cells. Toxicol Appl Pharmacol. 2016;310:9–19. https://doi.org/10.1016/j.taap.2016.08.005.
Chowanadisai W, Messerli SM, Miller DH, Medina JE, Hamilton JW, Messerli MA, et al. Cisplatin resistant spheroids model clinically relevant survival mechanisms in ovarian tumors. PLoS One. 2016;11(3):e0151089. https://doi.org/10.1371/journal.pone.0151089.
Kazmierczak D, Jopek K, Sterzynska K, Ginter-Matuszewska B, Nowicki M, Rucinski M, et al. The Significance of MicroRNAs Expression in Regulation of Extracellular Matrix and Other Drug Resistant Genes in Drug Resistant Ovarian Cancer Cell Lines. Int J Mol Sci. 2020;21(7):2619.
Xia J, Gill EE, Hancock RE. NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data. Nat Protoc. 2015;10(6):823–44. https://doi.org/10.1038/nprot.2015.052.
Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019;47(W1):W234–w241. https://doi.org/10.1093/nar/gkz240.
Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, et al. DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35(Web Server issue):W169–75.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–d613. https://doi.org/10.1093/nar/gky1131.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Assenov Y, Ramírez F, Schelhorn SE, Lengauer T, Albrecht M. Computing topological parameters of biological networks. Bioinformatics. 2008;24(2):282–4. https://doi.org/10.1093/bioinformatics/btm554.
Zhao X, Tang DY, Zuo X, Zhang TD, Wang C. Identification of lncRNA-miRNA-mRNA regulatory network associated with epithelial ovarian cancer cisplatin-resistant. J Cell Physiol. 2019;234(11):19886–94. https://doi.org/10.1002/jcp.28587.
Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4(1):2. https://doi.org/10.1186/1471-2105-4-2.
Pathan M, Keerthikumar S, Chisanga D, Alessandro R, Ang CS, Askenase P, et al. A novel community driven software for functional enrichment analysis of extracellular vesicles data. J Extracellular Vesicles. 2017;6(1):1321455. https://doi.org/10.1080/20013078.2017.1321455.
Pathan M, Keerthikumar S, Ang CS, Gangoda L, Quek CY, Williamson NA, et al. FunRich: an open access standalone functional enrichment and interaction network analysis tool. Proteomics. 2015;15(15):2597–601. https://doi.org/10.1002/pmic.201400515.
Fiannaca A, Rosa M, Paglia L, Rizzo R, Urso A. MiRNATIP: a SOM-based miRNA-target interactions predictor. BMC Bioinformatics. 2016;17(Suppl 11):321. https://doi.org/10.1186/s12859-016-1171-x.
Mao Y, Xue P, Li L, Xu P, Cai Y, Chu X, et al. Bioinformatics analysis of mRNA and miRNA microarray to identify the key miRNA-gene pairs in small-cell lung cancer. Mol Med Rep. 2019;20(3):2199–208. https://doi.org/10.3892/mmr.2019.10441.
Nagy Á, Lánczky A, Menyhárt O, Győrffy B. Validation of miRNA prognostic power in hepatocellular carcinoma using expression data of independent datasets. Sci Rep. 2018;8(1):9227. https://doi.org/10.1038/s41598-018-27521-y.
Sarver AL, Sarver AE, Yuan C, Subramanian S. OMCD: OncomiR Cancer database. BMC Cancer. 2018;18(1):1223. https://doi.org/10.1186/s12885-018-5085-z.
Chen YN, Ren CC, Yang L, Nai MM, Xu YM, Zhang F, et al. MicroRNA let-7d-5p rescues ovarian cancer cell apoptosis and restores chemosensitivity by regulating the p53 signaling pathway via HMGA1. Int J Oncol. 2019;54(5):1771–84. https://doi.org/10.3892/ijo.2019.4731.
Lloyd KL, Cree IA, Savage RS. Prediction of resistance to chemotherapy in ovarian cancer: a systematic review. BMC Cancer. 2015;15(1):117. https://doi.org/10.1186/s12885-015-1101-8.
Xiao S, Li Y, Pan Q, Ye M, He S, Tian Q, et al. MiR-34c/SOX9 axis regulates the chemoresistance of ovarian cancer cell to cisplatin-based chemotherapy. J Cell Biochem. 2019;120(3):2940–53. https://doi.org/10.1002/jcb.26865.
Jonas S, Izaurralde E. Towards a molecular understanding of microRNA-mediated gene silencing. Nat Rev Genet. 2015;16(7):421–33. https://doi.org/10.1038/nrg3965.
Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17(1):47–62. https://doi.org/10.1038/nrg.2015.10.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97. https://doi.org/10.1016/S0092-8674(04)00045-5.
Bach DH, Hong JY, Park HJ, Lee SK. The role of exosomes and miRNAs in drug-resistance of cancer cells. Int J Cancer. 2017;141(2):220–30. https://doi.org/10.1002/ijc.30669.
MacDonagh L, Gray SG, Finn SP, Cuffe S, O'Byrne KJ, Barr MP. The emerging role of microRNAs in resistance to lung cancer treatments. Cancer Treat Rev. 2015;41(2):160–9. https://doi.org/10.1016/j.ctrv.2014.12.009.
Chen X, Lu P, Wu Y, Wang DD, Zhou S, Yang SJ, et al. MiRNAs-mediated cisplatin resistance in breast cancer. Tumour Biol. 2016;37(10):12905–13. https://doi.org/10.1007/s13277-016-5216-6.
Fadejeva I, Olschewski H, Hrzenjak A. MicroRNAs as regulators of cisplatin-resistance in non-small cell lung carcinomas. Oncotarget. 2017;8(70):115754–73. https://doi.org/10.18632/oncotarget.22975.
Tang XD, Zhang DD, Jia L, Ji W, Zhao YS. lncRNA AFAP1-AS1 promotes migration and invasion of non-small cell lung Cancer via up-regulating IRF7 and the RIG-I-like receptor signaling pathway. Cell Physiol Biochem. 2018;50(1):179–95. https://doi.org/10.1159/000493967.
Yoshino H, Iwabuchi M, Kazama Y, Furukawa M, Kashiwakura I. Effects of retinoic acid-inducible gene-I-like receptors activations and ionizing radiation cotreatment on cytotoxicity against human non-small cell lung cancer in vitro. Oncol Lett. 2018;15(4):4697–705. https://doi.org/10.3892/ol.2018.7867.
Hou J, Zhou Y, Zheng Y, Fan J, Zhou W, Ng IO, et al. Hepatic RIG-I predicts survival and interferon-α therapeutic response in hepatocellular carcinoma. Cancer Cell. 2014;25(1):49–63. https://doi.org/10.1016/j.ccr.2013.11.011.
Sanchez-Vega F, Mina M, Armenia J, Chatila WK, Luna A, La KC, et al. Oncogenic Signaling Pathways in The Cancer Genome Atlas. Cell. 2018;173(2):321–337.e310.
Zheng Y, Pan D. The hippo signaling pathway in development and disease. Dev Cell. 2019;50(3):264–82. https://doi.org/10.1016/j.devcel.2019.06.003.
Voorhoeve PM, le Sage C, Schrier M, Gillis AJ, Stoop H, Nagel R, et al. A genetic screen implicates miRNA-372 and miRNA-373 as oncogenes in testicular germ cell tumors. Cell. 2006;124(6):1169–81. https://doi.org/10.1016/j.cell.2006.02.037.
Xu M, Xiao J, Chen M, Yuan L, Li J, Shen H. Yao S: miR-149-5p promotes chemotherapeutic resistance in ovarian cancer via the inactivation of the hippo signaling pathway. Int J Oncol. 2018;52(3):815–27. https://doi.org/10.3892/ijo.2018.4252.
Grosswendt S, Filipchyk A, Manzano M, Klironomos F, Schilling M, Herzog M, et al. Unambiguous identification of miRNA:target site interactions by different types of ligation reactions. Mol Cell. 2014;54(6):1042–54. https://doi.org/10.1016/j.molcel.2014.03.049.
Liang HL, Hu AP, Li SL, Xie JP, Ma QZ, Liu JY. MiR-454 prompts cell proliferation of human colorectal cancer cells by repressing CYLD expression. Asian Pac J Cancer Prev. 2015;16(6):2397–402. https://doi.org/10.7314/APJCP.2015.16.6.2397.
Yu L, Gong X, Sun L, Yao H, Lu B. Zhu L: miR-454 functions as an oncogene by inhibiting CHD5 in hepatocellular carcinoma. Oncotarget. 2015;6(36):39225–34. https://doi.org/10.18632/oncotarget.4407.
Zhu DY, Li XN, Qi Y, Liu DL, Yang Y, Zhao J, et al. MiR-454 promotes the progression of human non-small cell lung cancer and directly targets PTEN. Biomed Pharmacother. 2016;81:79–85.
Huang C, Liu J, Pan X, Peng C, Xiong B, Feng M, et al. miR-454 promotes survival and induces oxaliplatin resistance in gastric carcinoma cells by targeting CYLD. Exp Ther Med. 2020;19(6):3604–10. https://doi.org/10.3892/etm.2020.8655.
Tian Y, Ding W, Wang Y, Ji T, Sun S, Mo Q, et al. Ubiquitin B in cervical cancer: critical for the maintenance of cancer stem-like cell characters. PLoS One. 2013;8(12):e84457. https://doi.org/10.1371/journal.pone.0084457.
We deeply appreciate Xianyu Jin and Lifeng Liu for kind suggestions and generous help in data extraction and analysis.
Ethics approval and consent to participate
Consent for publication
The authors declare that there is no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xue, B., Li, S., Jin, X. et al. Bioinformatics analysis of mRNA and miRNA microarray to identify the key miRNA-mRNA pairs in cisplatin-resistant ovarian cancer. BMC Cancer 21, 452 (2021). https://doi.org/10.1186/s12885-021-08166-z
- Ovarian cancer
- Bioinformatics analysis
- miRNA-mRNA pairs