Establishment and Comprehensive Analysis of Underlying microRNA-mRNA Interactive Networks in Ovarian Cancer

Ye, Gengchen; Feng, Shuyue; Yang, Yufei; Cao, Zhengzhi; Zhang, Beilei; Wang, Fu

doi:https://doi.org/10.1155/2022/5120342

Journal of Oncology

On this page

Abstract Introduction Methods Results Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Bioinformatics in Cancer and Immune Microenvironment

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 5120342 | https://doi.org/10.1155/2022/5120342

Establishment and Comprehensive Analysis of Underlying microRNA-mRNA Interactive Networks in Ovarian Cancer

Gengchen Ye,¹Shuyue Feng,¹Yufei Yang,¹Zhengzhi Cao,¹Beilei Zhang,²and Fu Wang^1,3

Academic Editor: Feng Jiang

Received07 Jan 2022

Accepted15 Feb 2022

Published10 Mar 2022

Abstract

Background. The rate of ovarian cancer (OC) is one of the highest in women's reproductive systems. An improperly expressed microRNA (miRNA) has been discovered to have a vital role in the pathophysiology of OC. However, more research into OC's miRNA-message RNA (mRNA) gene interaction network is required. Methods. Firstly, the microarray data sets GSE25405 and GSE119055 from the GEO (Gene Expression Omnibus) database were downloaded and then analyzed with the GEO2R tool aiming at identifying DEMs (differential expressed miRNAs) between ovarian malignant tissue and ovarian normal tissue. The whole consistently changed miRNAs were then screened out to be candidate DEMs. For estimating underlying upstream transcription factors, FunRich was employed. miRNet was utilized to determine putative DEMs' downstream target genes. The R program was then used to do the GO annotation as well as the analysis of KEGG pathway enrichment for target genes. The PPI (protein-protein interaction), as well as the DEM-hub gene networks, were created by the Cytoscape software and STRING database. Finally, we chose the GSE74448 dataset to test the precision of hub gene expressions. Results. We have screened out six (five upregulated and one downregulated) DEMs. The majority of upregulated and downregulated DEMs are likely regulated by SP1 (specificity protein 1). SP4 (s protein 4), POU2F1 (POU class 2 homeobox 1), MEF2A (myocyte-specific enhancer factor 2A), ARID3A (AT-rich interaction domain 3A), and EGR1 (early growth response 1) can regulate upregulated and downregulated DEMs. We have found 807 target genes (656 upregulated and 151 downregulated DEM), being generally enriched in focal adhesion and proteoglycans in cancer, gastric cancer, hepatocellular carcinoma, as well as breast cancer. The majority of hub genes are projected to be controlled by hsa-miR-429, hsa-miR-140-5p, hsa-miR-199a-5p, and hsa-miR-199a-3p after the DEM-hub gene network was built. VEGFA (vascular endothelial growth factor A), EZH2 (enhancer of zeste 2 polycomb repressive complex 2 subunit), and HIF1A (hypoxia inducible factor 1 subunit alpha) expressions are consistent with the GSE74448 dataset in the first 18 hub genes. Conclusion. We have built an underlying miRNA-mRNA interacting network in OC, giving us unparalleled insight into the disease's diagnosis and treatment.

1. Introduction

OC is a malignancy in the female reproductive system, which is prevalent in patients over the age of 50. OC incidence has grown dramatically in Eastern Europe/Southern Europe and Asia in recent years [1]. Patients often show nonspecific pelvic or abdominal symptoms. Nowadays, carbohydrate antigen 125 (CA125) is the most widely utilized and concerning biomarker in OC. Nevertheless, CA125 is less reliable in the diagnosis and not capable of screening early ovarian cancer. In 2009, a randomized clinical trial in China showed that the routine surveillance of OC using the CA125 assay will increase the usage of chemotherapy, reduce the living quality of patients, and can’t improve the survival rate of patients [2]. The diagnosis of OC is so difficult as symptoms, medical history, and tests are usually not specific for OC. Because of its asymptomatic development, the disease is often diagnosed as advanced with poor treatment effects, and most patients will develop relapse, acquired drug resistance, and chemotherapy resistance. New early diagnosis and treatment methods are in urgent need.

miRNA, whose size is approximately 22 nucleotides, belongs to the group of small ncRNAs (noncoding RNAs). It regulates one-third of the genomes and downregulates the expression level of its target gene at the posttranscriptional stage as an antisense RNA. miRNA has been linked to a variety of human disorders, and it is thought to be the underlying target for cancer and other diseases detection and treatment [3]. Serum miRNA profiles are a potential and highly reliable biomarker for early detection and diagnosis of OC which has been proven in recent research [4]. miRNA might also operate as a wide regulator of gene expression, influencing particular target mRNA while also modulating the expression of several other genes. In tumor tissues, miRNA can be upregulated or downregulated, while a large proportion of miRNA seems to be upregulated in cancer [5]; for instance, researchers have demonstrated that members of the miR-200 family play an important role in the stress response of the human body and ovarian cancer formation [6]. Exosomal has-miR-21 inhibits ovarian cancer cell apoptosis in ovarian cancer cells while also conferring chemoresistance to cells by binding to the novel target APAF1 (apoptotic peptidase activating factor 1) [7]. Thus, miRNA may be useful as a noninvasive biomarker to improve the accuracy of ovarian cancer screening.

Although the miRNA-mRNA network has been explored in chronic obstructive pulmonary disease (COPD), there has not been an extensively elucidated miRNA-mRNA gene regulatory network in OC[8]. In our study, we analyzed 47 OC samples and 10 normal samples from two datasets (GSE25405 and GSE119055) from GEO in order to find DEMs in OC tissue with normal tissue, and 6 invariably changing miRNAs were filtered out. miRNet has successfully predicted the downstream target genes. The DEM target genes, DEM transcription factors, the DEM-hub gene network, as well as target gene function analysis have all been identified. Furthermore, in the GSE74448 dataset, the tissue level of the hub genes was evaluated and verified. In summary, we set out to discover particular targets for the diagnosis and treatment of OC by establishing an underlying miRNA-mRNA regulation network in the disease.

2. Methods

2.1. miRNA Microarray

To screen out appropriate miRNA datasets from the GEO (https://www.ncbi.nlm.nih.gov/geo/), our team has set four keywords: “Home sapiens” (organism), “noncoding RNA profiles by array” (study type), “OC” (study keywords), and “tissue” (attribute name). Eventually, two miRNA datasets (GSE25405 and GSE119055) were chosen for subsequent analysis, respectively, on the strength of the GPL7731 platform and the GPL21572 platform. The study design is depicted in Figure 1, and additional details on all the datasets are listed in Table 1.

2.2. The Screening of DEMs

GEO2R (http://www.ncbi.nlm.nih.gov/geo/geo2r/), an interactive web program for comparing at least two sets of samples in a GEO series in order to find DEMs in the same experimental conditions, was our choice. To detect DEMs in our investigation, we used log fold change (FC) <−2 or log fold change (FC) > 1.5 as well as a value <0.05 as criteria. The overlap of DEMs in the two datasets (GSE25405 and GSE119055) was analyzed using a Venn diagram, and these invariably changing miRNAs were identified as DEM candidates. The heatmaps of DEM candidates, as well as DEMs, are all created by the R program [9–18].

2.3. The Prediction of Transcription Factors Upstream of DEMs

A stand-alone software application FunRich (http://www.funrich.org/) was employed for forecasting underlying upstream transcription factors for potential DEMs. We chose it to study functional enrichment and protein-gene interaction networks. Statistical significance is determined by a value below 0.05.

2.4. The Prediction of Downstream Target Genes for DEMs

miRNet (https://www.mirnet.ca/) is a sophisticated toolkit that offers details on miRNA-target interactions as well as the ability to display the interactions in graphical networks. It is frequently used to anticipate downstream target genes of potential DEMs.

2.5. The GO and KEGG Analysis of Target Genes

For determining the biological function of DEM target genes, we employed R software to perform GO and KEGG pathway analysis. Statistical significance is determined by an adjusted value below 0.05.

2.6. The Construction of PPI Network and Screening of Hub Genes

We create a PPI (protein-protein interaction) network of DEM target genes employing STRING (http://stringdb.org/), then visualize it by using the Cytoscape software (version 3.9.0) and Cytoscape plug-in components, and sort nodes per the characteristics of the preloaded PPI network utilizing various topology algorithms. The hub gene (the top 10 nodes in the PPI network) is then screened out by using MCC (maximal clique centrality), which is substantially better than other methods for estimating the precision of essential proteins. Where S (v) represents the set of the biggest group including and the product of all positive integers is smaller, the proportion of MCC on node is defined as.

2.7. Expression Analysis of Hub Gene by GSE74448 Data Set

GSE74448, based on the platform of GPL18348 from the GEO database, was chosen to perform the analysis of the expressions of hub genes. In this dataset, we have picked up 29 OC and 11 normal tissues. Then, we performed the Student's t-test to screen out DEGs between OC and normal tissue. We have set to criteria, which significant hub genes must meet: one is that the upregulated DEM target genes are downregulated, or the downregulated DEM target genes are upregulated; the other is that the value must be below 0.05.

3. Results

3.1. Identification of DEMs in OC Tissue

In the GSE25405 dataset, 52 DEMs, including 37 upregulated and 15 downregulated, were identified after screening out under the standard of a value 0.05 and log fold change (FC) > 1.5 or log fold change (FC) -2. In the GSE119055 dataset, 140 DEMs, including 106 upregulated and 34 downregulated, were identified after screening out under the standard of a value 0.05 and log fold change (FC) > 1.5 or log fold change (FC) −2. Figures 2(a) and 2(b) show the heatmaps of DEMs, and Supplementary Tables S1 and S2 show the specific information about DEMs. Among these DEMs, five upregulated (has-miR-199a-3p, has-miR-140-3p, has-miR-199a-5p, has-miR-140-5p, and has-miR-455-5p) and one downregulated (has-miR-429) were invariably changed in both two datasets (Figures 3(a)–3(c)).

(a)

(b)

(a)

(b)

(c)

3.2. Prediction of Upstream Transcription Factors of DEMs

We employed a software named FunRich to forecast upstream transcription factors among candidate DEMs. Figures 4(a) and 4(b) are transcription factors that upregulate and downregulate DEM. The transcription factors were SP1, RORA (RAR related orphan receptor A), SP4, GATA1 (GATA binding protein 1), GFI1 (growth factor independent 1 transcriptional repressor), NKX2-1(NK2 homeobox 1), MEF2A, POU2F1, ARID3A, and EGR1 for downregulated DEMs. Concerning upregulated DEMs, transcription factors were NKX6-1 (NK6 homeobox 1), SP1 (specificity protein 1), SP4, EGR1, SOX1 (SRY-Box transcription factor 1), ARID3A, FOXA1 (forkhead box A1), POU2F1, FOXD1 (forkhead box D1), and MEF2A.

(a)

(b)

3.3. Prediction of Downstream Target Genes of DEMs

In all 865 candidate DEM target genes were found by the miRNet database (704 upregulated DEM target genes as well as 151 downregulated DEM target genes). Figure 5(a) and 5(b), respectively, show the upregulated, as well as downregulated DEMs, and the target genes of them, are represented in the DEM target gene network. Moreover, Table 2 shows the amount of target genes per DEM. Supplementary Table S3 shows the predicted target genes.

(a)

(b)

3.4. GO (Gene Ontology) and KEGG Analysis of Target Genes

On 865 DEM target genes, we have carried out GO and KEGG analysis. DEM target genes are abundant in epithelial cell proliferation, gland development, and positive regulation of kinase activity, according to biological process analysis (Figures 6(a)–6(f)), particularly; DEM target genes are mostly found at the cell-substrate junction and focal adhesion, according to cellular component analysis, and transcription regulator complex; DEM target genes are considerably enriched in protein serine/threonine kinase activity, protein serine kinase activity, and ubiquitin-like protein ligase binding, according to molecular function study.

(a)

(b)

(c)

(d)

(e)

(f)

KEGG pathway analysis shows that DEM target genes are significantly enriched in focal adhesion, proteoglycans in cancer, gastric cancer, hepatocellular carcinoma, and breast cancer (Figures 7(s) and 7(b)).

(a)

(b)

3.5. The Construction of PPI and DEM-Hub Gene Network

We used the STRING database to create a PPI (protein-protein interaction) network comprising DEM-targeted genes, which are upregulated and downregulated. Then, to screen out Hub genes, we used the Cytoscape and cytoHubba plugins. The first 10 hub genes, which are upregulated and downregulated in DEMs, are all shown in Figures 8(a) and 8(b). Concerning upregulated DEMs, the top 10 hub genes are MET (MET proto-oncogene and receptor tyrosine kinase), VEGFA, HIF1A, MTOR (mammalian target of rapamycin), GSK3B (glycogen synthase kinase 3 beta), SOX2 (SRY-box transcription factor 2), CDH2 (cadherin 2), IGF1R (insulin-like growth factor 1 receptor), FGF2 (fibroblast growth factor 2) and KDR (kinase insert domain receptor); and the top 10 hub genes for downregulated DEMs were MYC (MYC proto-oncogene and BHLH transcription factor), HIF1A, PTEN (MMAC1, mutated in multiple advanced cancers 1), VEGFA, KRAS (V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog), EZH2, JUN (Jun proto-oncogene and AP-1 transcription factor subunit), SOX2, EP300 (E1A binding protein P300) and ZEB1 (zinc finger E-box binding homeobox 1) (Table 3).

(a)

(b)

We used Cytoscape software to create a DEM-hub gene network to comprehend the molecular mechanism of the DEMs in OC tissue. The following are the relations between upregulated DEMs and hub genes: VEGFA, GSK3B, HIF1A, and CDH2 were associated with hsa-miR-199a-5p; hsa-miR-140-5p interacted with two hub genes, SOX2 and IGF1R; and hsa-miR-199a-3p interacted with MET, FGF2, KDR, and MTOR. Downregulated DEMs made interactions with 10 hub genes, including JUN, MYC, HIF1A, PTEN, SOX2, EZH2, KRAS, VEGFA, ZEB1, and EP300 (Figure 9).

3.6. Hub Genes Expression Validation

The GSE74448 dataset was employed for verifying tissue levels of the top four hub genes (EZH2, VEGFA, GSK3B, and HIF1A) based on the DEM-hub gene network. For downregulated DEMs, only the HIF1A expression surged in OC tissue (Figure 10(d)). Compared with normal tissue, the expression of upregulated DEMs (VEGFA and EZH2) was increased in OC tissue (Figures 10(a) and 10(b)). Therefore, hsa-miR-429-VEGFA, hsa-miR-429-EZH2, and hsa-miR-199a-5p-HIF1A were verified as three underlying regulatory pathways in OC tissues.

(a)

(b)

(c)

(d)

4. Discussion

In recent years, many studies on the diagnosis and treatment of OC have been carried out. However, the understanding of the underlying molecular biology pathogenesis of OC is still unclear, so the treatment effect of ovarian cancer patients is still poor. Through microarray technology, thousands of gene changes in the process of occurrence and development of various diseases can be revealed more easily, so as to construct an effective miRNA-mRNA gene regulatory network. In our study, we extracted data from two data sets, namely, GSE25405 and GSE119055, to distinguish DEMs between OC tissue and normal tissue.

Five upregulated DEMs (hsa-miR-455-5p, hsa-miR-199a-5p, hsa-miR-140-5p, hsa-miR-199a-3p, and hsa-miR-140-3p) and one downregulated DEM (hsa-miR-429) are constantly changing in the two data sets and have been identified as the candidate DEMs for further analysis. Both strands (-5p and -3p) of has-miR-455 target the same genes to regulate articular cartilage homeostasis, with underlying treatment value for osteoarthritis (OA), are upregulated by Sox9 (SRY-box transcription factor 9), a crucial transcription factor of cartilage differentiation and function [19]. has-miR-199a-5p is the key regulator of the abnormal α1-antitrypsin (AAT)-deficient monocyte unfolded protein response, which can be regulated by epigenetic silencing in chronic obstructive pulmonary disease [20]. has-miR-140-5p is an underlying prognostic factor for gastric cancer patients, being capable of mediating YES1 (YES proto-oncogene 1 and Src family tyrosine kinase) suppression [21]. It has appeared that the expression level of has-miR-199a-3p is concerned with acute heart failure (AHF) [22]. has-miR-140-3p downregulates the oncogene BRD9 (bromodomain containing 9) expression in squamous cell lung cancer (SqCLC) and inhibits SqCLC tumorigenesis, thus may become a new inhibitor of SqCLC [23]. As a member of the has-miR-200 family (composed of 5 miRNAs: has-miR-141, has-miR-200a, has-miR-200b, has-miR-200c, and has-miR-429), the expression of has-mir-429 is in connection with the survival rate of liver cancer [24]. It is worth noting that many of these miRNAs have been reported as underlying treatment targets for many human diseases, and their mechanism of action in OC requires further research to clarify more clearly.

In recent years, more and more research has revealed the fact that transcription factors are capable of regulating the expression of miRNA [25]. Therefore, we forecast that this process could potentially take place. SP1 (specific protein 1), a promoter-specific binding factor, has been proven to be overexpressed in a variety of cancers, and in connection with poor prognosis, downregulation of SP1 regulatory genes has also been shown to be drug-dependent. It is expected to account for the highest proportion of upregulated as well as downregulated DEM and has been widely studied in OC. Studies have shown that the new signal axis of has-miR-141/KLF12 (Kruppel-like factor 12)/SP1/survivin can enhance anoikis drug resistance, which is promising as a feasible therapeutic target for transmutative OC [26]. EGR1 is expected to regulate the upregulation and downregulation of DEM expression. Studies related to clinical cases have indicated that the expression levels of EGR1, as well as has-miR-152 in OC tissues, are significantly reduced, while EGR1-has-miR-152 downregulation of ATG14 (autophagy related 14) can make OC cells obtain cisplatin-induced apoptosis sensitivity by restraining cytoprotective autophagy [27].

The conclusion that DEMs are substantially enriched in focal adhesion, proteoglycans in cancer, and endocrine resistance is mainly indicated by the structure analyzed by KEGG. During cell migration, the dynamic regulation of focal adhesion formation and disassembly plays a crucial part in the metastasis of vast cancers including OC. SOCE (store-operated Ca entry) inhibitors significantly weaken the assembly and disassembly of focal adhesions. Therefore, SOCE may be a feasible medicinal target for OC [28]. The syndecans are a protein family formed by the joint combination of four single transmembrane domain proteins. Members can interact with plenty of ligands, including fibroblast growth factors, TGF-β (transforming growth factor-β), VEGF (vascular endothelial growth factor), and other molecules, and usually carry 3 to 5 heparan sulfate and chondroitin sulfate chains. The CD44 (CD44 molecule (Indian blood group)) receptor, a characteristic biomarker of ovarian cancer stem cells, is a group of intact membrane proteins with high sugar content, wide distribution, and a molecular weight of (85–160) × 10kD [29]. Endocrine therapy has been recognized as an effective way to treat chemically resistant ovarian cancer. Studies have shown that targeting the autocrine IL-6 (interleukin-6)-SPINK1 (serine peptidase inhibitor Kazal type) signal axis can inhibit the metastasis and spread of ovarian clear cell carcinoma, which provides a possibility while replicating the therapy. Therefore, inhibition of SOCE or syndecan-1, promoting CD44 receptor expression, and induction of the targeting autocrine IL-6-SPINK1 signal axis may represent a new treatment strategy for OC.

GO analysis results also have important reference values. In terms of molecular function (MF), ubiquitin-like protein ligase binding kinase regulatory activity, protein C-terminal binding, phosphatase binding, and growth factor binding have been shown to exhibit a significant correlation with OC. BRCA1 (breast cancer gene 1) is a nuclear phospholipid, which inhibits breast cancer and ovarian cancer tumors. It is usually combined with BARD1 (BRCA1 associated RING domain 1) to form an annular isoprenoid, acting as a ubiquitin (Ub) ligase. Its activity is greatly increased, and it contains a special domain that is interrupted by mutations that are susceptible to cancer [30]. CTBP2 (C-terminal binding protein 2) was found to be increased in epithelial ovarian carcinoma (EOC). Many research studies have demonstrated that CTBP2 and CTBP1 (C-terminal binding protein 1) have an extremely important role in suppressing the expression and activity of DR4/5 (death-receptors 4/5), and may maintain the possibility of CTBP as a high-grade serous ovarian cancer (HGSOC) treatment strategy [31]. In the pathological angiogenesis of female ovarian cancer, the expression level of vascular endothelial cadherin (VEC) and the expression of claudin-5 and vascular endothelial-protein tyrosine phosphatase are lowered in parallel [32]. TGF-β signaling is a considerable component in gene regulation in cancer. The TGF-β signal tends to increase in the late tumor malignant stage. Studies have shown that the expression of nc886 (a 101-nucleotide (nt)-long ncRNA) is induced by TGF-β, then nc886 is combined with dicer (an RNase III endonuclease) to inhibit miRNA maturation [33].

In terms of cellular components (CC), transcriptional regulatory complexes have been elucidated to play an essential part in the incidence and progression of OC. The overexpression of the oncogene MYBL2 (also known as B-Myb, MYB roto-oncogene like 2) has been confirmed to be associated with an increased rate of cell proliferation and can also become a valuable biomarker for poor cancer prognosis. However, the interaction of B-Myb with the MuvB (multivulval B, also known as LINC (LIN complex)) core complex in proliferating cells culminates in the formation of the product MMB (Myb-MuvB) complex, a process that is also critical for the transcription of genes required for mitosis promotion. Decreased or even loss of DREAM (drosophila, RB, E2F, and Myb) target gene repression in breast and ovarian cancer has been likewise demonstrated to be associated with overexpression of B-Myb. Based on this mechanism, targeting B-Myb to restore the body's control of the cell cycle may serve as a more effective and novel approach to treating not only ovarian cancer but other cancers as well [34].

In the biological process (BP), the positive regulation of kinase activity and hypoxia in the incidence and progression of OC has also been clarified to play a very essential role. To begin with, GSK-3 (Glycogen synthase kinase 3) can participate in the proliferation of female ovarian cancer cells by accelerating the expression level of cyclin D1, which is expected to be an explorable and therapeutic direction for OC [35]. The oxygen dependence of HIF1A protein works in the degradation and regulation of the activity of the transcription factor HIF (hypoxia-inducible factor), which is an important regulatory process in the hypoxia response. Past studies have demonstrated that the gene expression characteristics of hypoxia response are firm predictors of the clinical outcome of ovarian cancer and breast cancer and have an extremely important relationship with the significantly poorer prognosis of ovarian cancer and breast cancer [36].

By constructing the DEM-hub gene network, most of the hub genes may be targeted by hsa-miR-429, hsa-miR-199a-3p, hsa-miR-199a-5p, and hsa-miR-140-5p. In the first 18 central genes, the expression of only three genes (EZH2, VEGFA, and HIF1A) is consistent with the expression in the GSE74448 dataset, which may be due to the fact that they come from different sample sources. One of the epigenetic modifiers, EZH2 has been studied to promote ovarian cancer chemoresistance and recurrence [37]. VEGFA is a widely recognized treatment target of OC [38]. What we can know from the research is that has-miR-6086 can directly or indirectly downregulate the OC2 (onecut2)/VEGFA/EGFL6 (EGF-like domain multiple 6) axes to inhibit the angiogenesis network in OC, and may become a new research hotspot as an underlying target for tumor therapy [39]. HIF1A is a transcription factor that responds to hypoxia. By studying the transcription of inducible HIF1A-AS2 (hypoxia inducible factor 1 subunit alpha antisense RNA 2), there is hope for a new method for ovarian cancer treatment that combines effective anticancer therapy with OC standard-of-care chemotherapy [40].

Although we have uncovered an underlying miRNA-mRNA regulatory network in OC tissue, there were still some limitations in the current study. First, we concentrated on miRNAs and mRNAs only between OC and normal tissues, whereas a few of them may be distinct in different phases of OC. Second, through the prediction of the public database, we found miRNA-mRNA interactions.Further in vivo experiments are required to verify the analysis.

5. Conclusion

In conclusion, we have developed three underlying miRNA-mRNA pathways (hsa-miR-429-VEGFA, hsa-miR-429-EZH2, and hsa-miR-199a-5p-HIF1A) in OC tissues as underlying biomarkers for diagnosis and therapy of OC patients based on bioinformatics analysis and the GEO database. Our team truly hopes that all of the discoveries would be useful in further research and in improving the prognosis of patients with OC.

Data Availability

The supplementary information files provide the data used to support the findings of this investigation.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This paper was promoted by Natural Science Basic Research Program of Shaanxi (program no. 2022JC-56), the National Natural Science Foundation of China (no. 81772010), the “Young Talent Support Plan” of Xi’an Jiaotong University, and the Social Talent Grant of Tangdu Hospital of Air Force Medical University (grant no. 2021SHRC026).

Supplementary Materials

Table S1. DEMs between OC and normal tissue from the GSE25405 dataset. Table S2. DEMs between OC and normal tissue from the GSE119055 dataset. Table S3. Target genes of DEMs predicted by miRNet. (Supplementary Materials)

References

S. B. Coburn, F. Bray, M. E. Sherman, and B. Trabert, “International patterns and trends in ovarian cancer incidence, overall and by histologic subtype,” International Journal of Cancer, vol. 140, no. 11, 2017.
View at: Google Scholar
K. M. Esselen, A. M. Cronin, and K. Bixel, “Use of CA-125 tests and computed tomographic scans for surveillance in ovarian cancer,” JAMA Oncology, vol. 2, no. 11, 2016.
View at: Publisher Site | Google Scholar
S. M. Hammond, “An overview of microRNAs,” Advanced Drug Delivery Reviews, vol. 87, 2015.
View at: Publisher Site | Google Scholar
A. Yokoi, J. Matsuzaki, and Y. Yamamoto, “Integrated extracellular microRNA profiling for ovarian cancer screening,” Nature Communications, vol. 9, no. 1, 2018.
View at: Publisher Site | Google Scholar
J. Gerhard, H.-I. Eva, and M. Leticia, “Epigenetics of colorectal cancer: biomarker and therapeutic potential,” Nature Reviews Gastroenterology & Hepatology, vol. 17, no. 2, 2020.
View at: Google Scholar
L. Batista, B. Bourachot, B. Mateescu, F. Reyal, and F. Mechta-Grigoriou, “Regulation of miR-200c/141 expression by intergenic DNA-looping and transcriptional read-through,” Nature Communications, vol. 7, no. 1, p. 8959, 2016.
View at: Publisher Site | Google Scholar
L. Chi, N. N. Co, and T. Tsuruga, “Exosomal transfer of stroma-derived miR21 confers paclitaxel resistance in ovarian cancer cells through targeting APAF1,” Nature Communications, vol. 7, p. 11150, 2016.
View at: Google Scholar
M. Zhu, M. Ye, J. Wang, L. Ye, and M. Jin, “Construction of potential miRNA-mRNA regulatory network in COPD plasma by bioinformatics analysis,” International Journal of Chronic Obstructive Pulmonary Disease, vol. 15, pp. 2135–2145, 2020.
View at: Publisher Site | Google Scholar
R. Becker, The New S Language, CRC Press, Boca Raton, FL, USA, 1988.
L. Sabine and L. Morven, Cluster Analysis, John Wiley & Sons, NJ, USA, 2011.
J. A. Hartigan, Clustering Algorithms, John Wiley & Sons, , NJ, USA, 1918.
R. E. Anderberg, Cluster Analysis for Applications.
A. D. Gordon, Classification, Chapman & Hall, London, UK, 2nd edition, 1999.
L. L. Mcquitty, “Similarity analysis by reciprocal pairs for discrete and continuous data,” Educational and Psychological Measurement, vol. 26, no. 4, 1966.
View at: Publisher Site | Google Scholar
F. Murtagh and P. Legendre, “Ward's hierarchical agglomerative clustering method: which algorithms implement ward's criterion?” Journal of Classification, vol. 31, no. 3, 2014.
View at: Publisher Site | Google Scholar
R. R. Sokal, “Numerical taxonomy,” Scientific American, vol. 215, no. 6, pp. 106–116, 1966.
View at: Publisher Site | Google Scholar
F. Murtagh, Multidimensional Clustering Algorithms, Compstat Lectures Vienna Physika Verlag, Heidelberg, Germany, 1985.
C. F. Dormann, Basic and Applied Ecology, 2013.
Y. Ito, T. Matsuzaki, F. Ayabe et al., “Both microRNA-455-5p and -3p repress hypoxia-inducible factor-2α expression and coordinately regulate cartilage homeostasis,” Nature Communications, vol. 12, no. 1, p. 4148, 2021.
View at: Publisher Site | Google Scholar
T. Hassan, T. P. Carroll, P. G. Buckley et al., “miR-199a-5p silencing regulates the unfolded protein response in chronic obstructive pulmonary disease and α1-antitrypsin deficiency,” American Journal of Respiratory and Critical Care Medicine, vol. 189, no. 3, pp. 263–273, 2014.
View at: Publisher Site | Google Scholar
Z. Fang, S. Yin, R. Sun et al., “miR-140-5p suppresses the proliferation, migration and invasion of gastric cancer by regulating YES1,” Molecular Cancer, vol. 16, no. 1, p. 139, 2017.
View at: Publisher Site | Google Scholar
S. O. Ekaterina, S. Daniela, and L. V. Eline, “Signature of circulating microRNAs in patients with acute heart failure,” European Journal of Heart Failure, vol. 18, no. 4, 2016.
View at: Google Scholar
H. Huang, Y. Wang, and Q. Li, “miR-140-3p functions as a tumor suppressor in squamous cell lung cancer by regulating BRD9,” Cancer Letters, vol. 446, 2019.
View at: Google Scholar
G.-L. Huang, J. Sun, Y. Lu et al., “MiR-200 family and cancer: from a meta-analysis view,” Molecular Aspects of Medicine, vol. 70, pp. 57–71, 2019.
View at: Publisher Site | Google Scholar
L. Mingyang, Z. Yuqing, and Y. Jingxuan, “Zinc dependent regulation of ZEB1 and YAP1 Co-activation promotes EMT plasticity and metastasis in pancreatic cancer,” Gastroenterology, vol. 160, no. 5, 2021.
View at: Google Scholar
M. C. S. Mak, Y. M. M. Yung, and H. L. M. David, “MicroRNA-141 enhances anoikis resistance in metastatic progression of ovarian cancer through targeting KLF12/Sp1/survivin axis,” Molecular Cancer, vol. 16, no. 1, 2017.
View at: Google Scholar
H. Jun, Y. Jing-Jie, and X. Qing, “Downregulation of ATG14 by EGR1-MIR152 sensitizes ovarian cancer cells to cisplatin-induced apoptosis by inhibiting cyto-protective autophagy,” Autophagy, vol. 11, no. 2, 2015.
View at: Google Scholar
H. Hokai, L. Yihsin, and C. Hengai, “Chemoresistant ovarian cancer enhances its migration abilities by increasing store-operated Ca2+ entry-mediated turnover of focal adhesions,” Journal of Biomedical Science, vol. 27, no. 1, 2020.
View at: Google Scholar
C. D. Erik, W. M. Stephen, and E. S. Kevin, “Bimodal tumor-targeting from microenvironment responsive hyaluronan layer-by-layer (LbL) nanoparticles,” ACS Nano, vol. 8, no. 8, 2014.
View at: Google Scholar
S. B. Peter, R. K. Jennifer, and N. Hiroyuki, “Binding and recognition in the assembly of an active BRCA1/BARD1 ubiquitin-ligase complex,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 10, 2003.
View at: Google Scholar
D. Boxiao and Y. Fang, “CtBP determines ovarian cancer cell fate through repression of death receptors,” Cell Death & Disease, vol. 11, no. 4, 2020.
View at: Google Scholar
M. F. Morini, C. Giampietro, M. Corada et al., “VE-Cadherin-Mediated epigenetic regulation of endothelial gene expression,” Circulation Research, vol. 122, no. 2, pp. 231–245, 2018.
View at: Publisher Site | Google Scholar
A. Ji-Hye, L. Hyun-Sung, and L. Ju-Seog, “nc886 is induced by TGF-β and suppresses the microRNA pathway in ovarian cancer,” Nature Communications, vol. 9, no. 1, 2018.
View at: Google Scholar
N. I. Audra, F. Jessica, and A. Varsha, “The cell cycle regulatory DREAM complex is disrupted by high expression of oncogenic B-Myb,” Oncogene, vol. 38, no. 7, 2019.
View at: Google Scholar
Q. Cao, X. Lu, and Y.-J. Feng, “Glycogen synthase kinase-3β positively regulates the proliferation of human ovarian cancer cells,” Cell Research, vol. 16, no. 7, pp. 671–677, 2006.
View at: Publisher Site | Google Scholar
J.-T. Chi, Z. Wang, and D. S. A. Nuyten, “Gene expression programs in response to hypoxia: cell type specificity and prognostic significance in human cancers,” PLoS Medicine, vol. 3, no. 3, 2006.
View at: Publisher Site | Google Scholar
W. Yiping, H. Yaya, and Y. Xiaoqing, “EZH2 activates CHK1 signaling to promote ovarian cancer chemoresistance by maintaining the properties of cancer stem cells,” Theranostics, vol. 11, no. 4, 2021.
View at: Google Scholar
J. Kibeom, K. Minsoon, and A. C. Gilbert, “VEGFA activates an epigenetic pathway upregulating ovarian cancer-initiating cells,” EMBO Molecular Medicine, vol. 9, no. 3, 2017.
View at: Google Scholar
B. Wu, L. Zhang, Y. Yu et al., “miR-6086 inhibits ovarian cancer angiogenesis by downregulating the OC2/VEGFA/EGFL6 axis,” Cell Death & Disease, vol. 11, no. 5, p. 345, 2020.
View at: Publisher Site | Google Scholar
L. Tiangong, T. Jianming, and S. Binita, “Up-regulation of hypoxia-inducible factor antisense as a novel approach to treat ovarian cancer,” Theranostics, vol. 10, no. 15, 2020.
View at: Google Scholar

Copyright

Copyright © 2022 Gengchen Ye et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

461

Downloads

753

Citations