High-throughput sequencing identification of differentially expressed microRNAs in metastatic ovarian cancer with experimental validations
Cancer Cell International volume 20, Article number: 517 (2020)
Ovarian cancer (OC) is a common gynecological cancer and characterized by high metastatic potential. MicroRNAs (miRNAs, miRs) have the promise to be harnessed as prognostic and therapeutic biomarkers for OC. Herein, we sought to identify differentially expressed miRNAs and mRNAs in metastatic OC, and to validate them with functional experiments.
Differentially expressed miRNAs and mRNAs were screened from six pairs of primary OC tissues and metastatic tissues using a miRStar™ Human Cancer Focus miRNA and Target mRNA PCR Array. Then, gene expression profiling results were verified by reverse transcription quantitative polymerase chain reaction (RT-qPCR) and western blot assays. The binding affinity between miR-7-5p and TGFβ2 was validated by dual-luciferase reporter assay. Expression of miR-7-5p and TGFβ2 was manipulated to assess their roles in malignant phenotypes of highly metastatic HO-8910PM cells.
MiRNA profiling and sequencing identified 12 miRNAs and 10 mRNAs that were differentially expressed in metastatic tissues. Gene ontology and Pathway analyses determined that 3 differentially expressed mRNAs (ITGB3, TGFβ2 and TNC) were related to OC metastasis. The results of RT-qPCR confirmed that the decrease of miR-7-5p was most significant in OC metastasis, while TGFβ2 was up-regulated in OC metastasis. Moreover, miR-7-5p targeted and negatively regulated TGFβ2. MiR-7-5p overexpression accelerated HO-8910PM cell viability and invasion, and TGFβ2 overexpression reversed the results. Meanwhile, simultaneous miR-7-5p and TGFβ2 overexpression rescued the cell activities.
This study characterizes differentially expressed miRNAs and mRNAs in metastatic OC, where miR-7-5p and its downstream target were most closely associated with metastatic OC. Overexpression of miR-7-5p targets and inhibits TGFβ2 expression, thereby inhibiting the growth and metastasis of OC.
Ovarian cancer (OC) is the deadliest gynecological cancer, due to the absence of symptoms at the early stage [1, 2]. The non-specific symptoms of OC contribute to a delayed diagnosis at advanced stages of cancer, which causes poor 5-year overall survival rate that is about 30–40% . Moreover, the widespread and distant metastases are observed in 59% OC patients, which are clearly correlated with unfavorable prognosis of OC patients . Therefore, it is pivotal to unveil invasive and metastatic mechanisms of OC. Of note, accurate prognostic strategies with the assistance of biomarker targets would aid in the diagnosis and prognosis of OC patients . However, the understanding of mechanisms behind biomarkers and signaling pathways controlling invasive and metastatic potential of OC should been thoroughly probed .
MicroRNAs (miRNAs), small RNA molecules [21–23 nucleotides (nt)], function as potent modulators of gene expression through mRNA translation blockade or RNA interference . MiRNAs have become a hotspot of tumor researches over the recent years. Interestingly, the differential expression of miRNAs has been reported in OC. For example, in a previous study conducted by Chong et al.  differential miRNA expression profiles showed that miR-551b, miR-19b, miR-196b and miR-3198 were upregulated and miR-8084, miR-3201, miR-3613 and miR-7515 were downregulated in OC. Besides, another study revealed that there were 17 differentially expressed miRNAs in metastatic OC, among which miR-21, miR-150, and miR-146a were upregulated in metastatic OC tumors .
Moreover, the dysregulated miRNAs may act as novel oncogenes or anti-oncogenes in OC by regulating target genes . For instance, miR-7 was found to block the EGFR/ERK pathway to decrease OC cell invasion . Stated thus, a hypothesis can be made that miRNAs may participate in OC metastasis via modulating their target genes. Notably, emerging evidence has demonstrated the importance of the miRNA-mRNA interplay in the pathogenesis and metastasis of OC [12, 13]. Therefore, differentially expressed miRNAs and mRNAs were screened from primary and metastatic tissues of OC using high-throughput sequencing technique. We found 12 miRNAs and 10 mRNAs that were differentially expressed in and related to metastatic potential of OC. The expression profiling results were validated by functional experimental validations.
All 31 pairs of metastatic tissues (lymph node metastasis tissue) and primary tissues of OC were harvested from patients receiving surgery at Shengjing Hospital of China Medical University from 2014 to 2016. Of these, 6 pairs were used for expression profiling experiments, and 25 pairs were employed for RT-qPCR validation. Tissue samples resected during surgery were immediately placed into 1.5-mL centrifuge tubes with RNase-removal high-pressure treatment, and rapidly placed in liquid nitrogen for quick freezing, and then transferred to a − 80 ˚C cryogenic freezer for preservation. All patients were diagnosed with ovarian serous cystadenocarcinoma by postoperative paraffin pathology, and their clinical staging was made according to the International Federation of Gynecology and Obstetrics (FIGO) staging system for OC (2009). None of the patients had undergone preoperative radiotherapy, chemotherapy or other special treatments. Patients who provided tissue specimens for expression profiling experiments were aged 48–64 years with the mean age of 57.5 years, and those providing samples for RT-qPCR validation were aged 38–80 years with the mean age of 57 years. Data and tissues were harvested upon receiving the informed consent of patients and approval by the ethics committee of Shengjing Hospital of China Medical University (Approval number: 2016PS040K).
RNA extraction from tissues
A Trizol Kit (Invitrogen Inc., Carlsbad, CA, USA) was applied for isolation of total RNA from metastatic and primary tissues of OC. RNA purity was estimated by adopting a NanoDrop® ND-1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA), and assessed by the ratio of optical density (OD) 260 nm/280 nm. A ratio between 1.8 and 2.0 signified the sample to be highly pure.
MiRNA profiling and sequencing
Six pairs of metastatic and primary tissues of OC were obtained from Shengjing Hospital of China Medical University. Total RNA was extracted from tissues using a Trizol Kit. cDNA was synthesized, purified, and hybridized based on the instructions of an Arraystar miRNA First-Strand cDNA Synthesis Kit (Cat# AS-MR-004), which was fully compatible with an miRStar™ Human Cancer Focus miRNA and Target mRNA PCR Array (Arraystar, Rockville, MD, USA). After standardization of the original data from fluorescence quantitative PCR, fold changes in miRNAs and mRNAs were calculated based on a comparison between metastatic and primary tissues of OC. A fold change > 1.5 and p < 0.05 were considered as a noteworthy up-regulation, while fold change < − 1.5 and p < 0.05 were considered as a striking down-regulation.
Gene ontology (GO) functional enrichment and pathway analyses
GO analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were employed to analyze differentially expressed mRNAs in OC based on a David database (https://david.ncifcrf.gov/). GO analysis provided annotations and function classifications for differentially expressed mRNAs by analyzing the classifications of Biological Processes, Cellular Components and Molecular Function. Pathway analysis was based on KEGG results combined with Fisher's exact test and a t test. A p < 0.05 was regarded statistically distinct.
The miRNA-mRNA network analysis
Using a miRStar™ Human Cancer Focus miRNA and Target mRNA PCR Array expression profiling, we carried out the screening of 6 pairs of metastatic and primary tissues of OC in 184 cancer-related miRNAs and 178 target mRNAs. The miRNA-mRNA networks with a potential targeting relationship were screened out from miRNA microarray analysis and then plotted with the assistance of online tools targetscan (https://www.targetscan.org/vert_71/) and miranda (https://18.104.22.168/microrna/home.do).
The metastatic and primary tissues of OC frozen at − 80 °C were collected. TRIzol reagent (Invitrogen) was applied for total RNA isolation from tissues and cells. A mixture (14.5 µL) of 2.0 µg RNA, 1 µL Oligo, 1.6 µL deoxy-ribonucleoside triphosphate Mix and H2O was incubated at 65 °C in water for 5 min and then put on ice for 2 min. After centrifugation, the mixture was mixed with reverse transcription reaction liquid, and placed at 37 °C for 1 min. The mRNA was reversely transcribed into cDNA, and the samples were then temporarily placed on ice. All cDNA template samples were prepared in a PCR amplification system. The fluorescence value of SYBR Green was detected. PCR primer sequences (Table 1) were designed and synthesized by Sangon Biotech (Shanghai, China). The relative expression of mRNAs or miRNAs, normalized to β-actin or U6, was calculated using 2−ΔΔCt method.
Tissues were immersed in radioimmunoprecipitation assay lysis buffer (P0013B; Beyotime, Shanghai, China) containing phenylmethyl sulfonylfluoride to extract total protein. The samples were allowed to stand on ice for 30 min, followed by 20-min centrifugation at 12,000×g. The supernatant was collected and transferred to an EP tube to determine the total protein concentration according to the instructions of a bicinchoninic acid kit (KC-430, Shanghai Kangchen Biological Engineering Co., Ltd., Shanghai, China). The samples were separated at 100 V for 1.5 h, and electroblotted onto a polyvinylidene fluoride membrane using a sandwich method at 80 V for 1.5 h. After 1-h blocking in 5% bovine serum albumin-Tris-buffered saline (TBST), the membrane was incubated with diluted (1:1000) rabbit anti-human TGFβ2 primary antibody overnight at 4 °C. Afterwards, the membrane was then re-probed with goat anti-rabbit secondary antibody (1:5000; Shanghai Kangchen Biological Engineering Co., Ltd., Shanghai, China) complexed to horseradish peroxidase at RT for 1 h. After rinsing with TBST, the membrane was soaked in electrogenerated chemiluminescence (P0018, Beyotime) for 1 min, developed, photographed and recorded. β-actin (1:10,000; Cell Signaling Technology, Danvers, MA, USA) served as an endogenous control for TGFβ2. A Gel-Pro Analyzer 4.0 program (Media Cybernetics, Silver Spring, MD, USA) was adopted for gray value analysis. The gray value ratio of each target protein band to the endogenous control band was calculated as the relative expression of each protein. Experiments were independently conducted three times.
Dual-luciferase reporter assay
Binding sites for miR-7-5p and TGFβ2-3′untranslated region (UTR) were predicted using Genebank and Targetscan combined with expression profiling analysis results. TGFβ2-3′-UTR wild type (WT), TGFβ2-3′-UTR mutant type (MUT) luciferase plasmid and miR-7-5p mimic/negative control (NC) were constructed by Shanghai GeneChem Co., Ltd (Shanghai, China). After plasmid transfection, commonly used tool HEK293T cells in the logarithmic growth phase were prepared into a suspension that was then cultured in a 24-well plate so as to achieve 85% confluence at 24 h. Four groups were used in the experiment (NC mimic + TGFβ2-WT, miR-7-5p mimic + TGFβ2-WT, NC mimic + TGFβ2-MUT, and miR-7-5p mimic + TGFβ2-MUT). After 48-h transfection, the cells were lysed to obtain the supernatant that was centrifuged for 3–5 min. The luciferase activity of TGFβ2 3′UTR in response to miR-7-5p mimic was detected by utilizing a Dual-Luciferase® Reporter Assay System (E2910; Promega, Fitchburg, WI, USA). The luminescence intensity was measured by GloMax®20/20 luminometer (Promega, Madison, WI, USA).
Cell treatment and grouping
The cells used for the cell phenotype detection were the human highly metastatic OC cell line HO-8910PM, which has also been previously reported to establish a highly metastatic animal model [14, 15]. Besides, the human normal ovarian cell line HUM-CELL-0088 was selected as a control. Both cell lines were purchased from Shanghai Institute of Biochemistry and Cell Biology (Shanghai, China). HEK293T cells (used for dual-luciferase reporter assay) were cultured in Dulbecco minimum essential medium (DMEM; 12800017; Gibco, Carlsbad, CA, USA) containing 10% fetal bovine serum (FBS; FBS500-S; Ausgenex, Brisbane, Australia) in an incubator (BB15; Thermo Fisher Scientific, Waltham, MA, USA) with saturated humidity and 5% CO2 at 37 °C. HO-8910-PM cells and HUM-CELL-0088 cells were cultured in RPMI-1640 medium (Gibco-BRL, Life Technologies, Gaithersburg, MD, USA) containing 10% FBS in an incubator with saturated humidity and 5% CO2 at 37 °C. The medium was renewed every 24 h, and subculture was performed every 72 h. After removal of medium, cells were treated with 0.25% trypsin for 3 min, and then the reaction was ended by adding RPMI-1640 medium containing 10% FBS. Afterwards, cells were pipetted into single cell suspension. Then HO-8910-PM cells were transduced with lentiviral vectors containing miR-7-5p mimic, NC mimic, miR-7-5p inhibitor, NC inhibitor, TGFβ2 overexpression plasmid (oe-TGFβ2) or oe-NC. Lentiviruses were purchased from Shanghai Genechem Co., Ltd. (Shanghai, China). The negative control and miR-7-5p mimic and inhibitor were obtained from GenePharma (Shanghai, China). The sequence of the negative control, miR-7-5p mimic and inhibitor were provided in Additional file 1: Table S1. Cells in logarithmic growth phase were treated with trypsin to prepare into 3 × 104 cells/mL cell suspension. Cell suspension (2 mL) was seeded into 6-well plates for further culture. The infection medium (Eni.S + polybrene) was renewed, and the cell density was 20%/well at the time of infection. The optimal amount (infection MOI = 10) was added for infection. Twelve hours after infection, the conventional medium was renewed to continue culture cells.
Cell counting kit (CCK-8) assay
The 80% confluent cells were dispersed into single cell suspension with the reaction of 0.25% trypsin. After counting, the cells were separately seeded in a 96-well plate with 100 μL (2 × 104 cells/mL) per well, and set with 6 duplicated wells. Cells were cultured in the incubator, and 10 μL CCK-8 solutions (Sigma-Aldrich, St. Louis, MO, USA) were supplemented to each well after culture for 24 h, 48 h, 72 h, and 96 h. After further 2-h culture, OD value at 450 nm was measured using an enzyme-linked immunometric meter. Cell viability curve was drawn with time as X-axis and OD value as Y-axis.
After transfection, cells were suspended in 0.25% trypsin (C0205, Beyotime), followed by centrifugation at 400 rpm for 5 min with the supernatant discarded. Subsequently, cells were resuspended in RPMI-1640 medium. Cell density was adjusted to 105 cells/mL. Before the experiment, Transwell chambers (3413, Millipore, Massachusetts, USA) were coated with Matrigel. The 100 μL cell suspension was supplemented into the 24-well apical chamber (105 cells/well). RPMI-1640 medium (600 μL) containing 30% FBS was added into the basolateral chamber. Each group was set 3 duplicated wells. After 48-h culture, cells not invading through Matrigel were discarded using cotton swabs, and cells invading through Matrigel were stained with 0.1% Giemsa (32884, Sigma-Aldrich) for 10 min. The number of cells invading through Matrigel was counted in five random fields under the microscope. The number of cells invading through Matrigel was considered as an indicator of the invasion ability.
Even horizontal linear scratches were made in the back of 6-well plate by a marker pen with 3 scratches in each well. Six pairs of cells at logarithmic growth phase were seeded in the culture plate (5 × 105 cells/well). Following 24-h culture, when cells settled at the bottom of the well, 200 μL pipette tips was adopted to make even vertical linear scratches in the back of the plate with 3 scratches in each well. Afterwards, cell incubation was conducted in RPMI-1640 medium with 5% CO2 at 37 °C. The migration distance was observed and photographed under inverted microscope at 0 h and 24 h after scratch. Three replicates were made in each group. The number of cells passing through scratches was considered as an indicator of the migration ability.
Overnight cell fixative was conducted in 70% chilled ethanol at 4 °C. After centrifugation at 800×g and 4 °C with the removal of supernatant, cells were washed with phosphate buffer saline (PBS) containing 1% FBS. Cells were resuspended with 400 μL binding buffer, followed by 30-min incubation with 50 μL RNAase at 37 °C. After that, cells were added with 50 μL propidium iodide (PI) (50 mg/L, Sigma-Aldrich) for 30-min incubation in dark. Flow cytometry was used to detect cell cycle.
All data were summarized as mean ± standard deviation and analyzed using SPSS 21.0 statistical software (IBM Corp., Armonk, New York, USA). Data between two groups were compared using unpaired t-test. Statistical analysis in relation to time-based measurements within each group was realized using repeated measures analysis of variance (ANOVA), followed by a Bonferroni’s post-hoc test. Pearson correlation coefficient was adopted for correlation analyses. Significance was defined as p < 0.05.
Differentially expressed miRNAs and mRNAs with metastatic OC
Differentially expressed miRNAs and mRNAs in metastatic OC were screened using miRStar™ Human Cancer Focus miRNA and Target mRNA PCR Array analysis. As shown in Table 2, 12 differentially expressed miRNAs were screened out from metastatic and primary tissues of OC (Fig. 1a): (in sequence of absolute fold change) miR-200a-3p, miR-141-3p, miR-200b-3p, miR-10a-5p, miR-15a-5p, miR-187-5p, miR-16-5p, miR-9-3p, miR-195-5p, miR-7-5p, miR-584-5p and miR-27a-3p (p < 0.05). The expression of miR-10a-5p, miR-15a-5p, miR-16-5p, miR-195-5p and miR-27a-3p in metastatic tissues from OC was increased, as compared with primary OC tissues (p < 0.05, fold change > 1.5), while that of miR-141-3p, miR-7-5p, miR-187-5p, miR-200a-3p, miR-200b-3p, miR-584-5p and miR-9-3p was lower than that in primary OC tissues (p < 0.05, fold change < − 1.5). As shown in Table 3, a total of 10 differentially expressed target mRNAs were identified. Relative to primary OC tissues, the gene expression of BCL2L1, TGFβ2, TNC, HOXB5, HOXB7, HOXB8, ITGB3 and MTSS1 in metastatic tissues from OC was substantially up-regulated (p < 0.05, fold change > 2; Fig. 1b), while that of PAK1 and FOXA2 was down-regulated (p < 0.05, fold change < − 2).
GO functional enrichment and Pathway analyses for the identified differentially expressed mRNAs
To determine which mRNAs were involved in the regulatory effect of miRNAs in metastatic OC, GO analysis (Fig. 2) was performed to investigate the function of 10 differentially expressed mRNAs. For Biological Processes, 10 differentially expressed mRNAs mostly participated in anatomical structure morphogenesis, embryo development, cell differentiation and cellular developmental processes (p < 0.001, Enrichment scores > 8; Table 4). For cellular components, the differentially expressed mRNAs mostly participated in ruffle, adherens junction, anchoring junction and cell junction (p < 0.001, Enrichment scores > 3.3; Table 5). For molecular functions, the differentially expressed mRNAs mostly participated in RNA polymerase II distal enhancer sequence-specific DNA binding transcription factor activity involved in positive regulation of transcription (p < 0.001, Enrichment scores > 3; Table 6). Furthermore, the Pathway analysis based on KEGG method (Fig. 3) indicated that, among the 10 differentially expressed mRNAs, ITGB3, TGFβ2 and TNC were linked to the pathway of MicroRNAs in Cancer.
The down-regulation of miR-7-5p is the most significant in metastatic tissue of OC
Next, we moved to assess expression of the screened differentially expressed miRNAs as mentioned above, so as to determine the research focus for this present study. As illustrated in Fig. 4, when compared with primary OC tissues, there were 5 miRNAs miR-7-5p, miR-141-3p, miR-187-5p, miR-200a-3p and miR-200b-3p were appreciably decreased in metastatic tissues, where miR-7-5p expression presented with the most significant decrease. Thus, miR-7-5p was selected as a research focus for following experiments.
The up-regulation of TGFβ2 is the most significant in metastatic tissue of OC
Further, we proceeded to probe the downstream target genes of aforementioned miR-7-5p. A differentially expressed miRNA-mRNA network (Fig. 5a) was plotted, where we found that the target genes of miR-7-5p were TGFβ2 and HOXB5. Therefore, mRNA expression of TGFβ2 and HOXB5 was measured by RT-qPCR (Fig. 5b, c). Relative to the primary OC tissues, expression of TGFβ2 mRNA was enhanced in metastatic tissues (p < 0.05), while no significant difference was witnessed in the slightly up-regulated expression of HOXB5 mRNA (p > 0.05). This may be implicated in factors, such as sample size and ethnic differences. In addition, western blot assay (Fig. 5d) further confirmed that the protein level of TGFβ2 was elevated in metastatic tissues versus that in primary OC tissues (p = 0.0014). Accordingly, we selected miR-7-5p and its downstream target TGFβ2 for further investigations.
miR-7-5p targets and inhibits TGFβ2
We used the commonly used tool cell HEK293T cells to conduct a dual luciferase reporter gene assay to verify the targeting relationship between miR-7-5p and TGFβ2. It was found that the luciferase activity of TGFβ2-WT was reduced in response to miR-7-5p mimic transfection, while there was no significant difference in luciferase activity of TGFβ2-MUT (Fig. 6a, b), suggesting that miR-7-5p could target TGFβ2. Furthermore, we transfected miR-7-5p mimic or miR-7-5p inhibitor into highly metastatic OC cell line HO-8910PM to validate its regulatory effect on TGFβ2. As shown in Fig. 6c, we have also tested the transfection efficiency of miR-7-5p inhibitor and miR-7-5p mimic. The RT-qPCR detection presented that miR-7-5p expression was upregulated in the presence of miR-7-5p mimic and downregulated in the presence of miR-7-5p inhibitor. In addition, Western blot analysis manifested that protein levels of TGFβ2 was diminished in response to miR-7-5p mimic and elevated in response to miR-7-5p inhibitor in highly metastatic HO-8910PM cells (Fig. 6c, d). As shown in Fig. 6c, we have also tested the transfection efficiency of miR-7-5p inhibitor and miR-7-5p mimic. Besides, miR-7-5p expression was inversely correlated with TGFβ2 expression in metastatic tissues from OC (Fig. 6e).
Additionally, we also examined the expression of miR-7-5p and TGFβ2 in normal ovarian cell line HUM-CELL-0088 and highly metastatic OC cell line HO-8910PM. The RT-qPCR results (Fig. 6f) showed that, compared with normal ovarian cells, miR-7-5p was poorly expressed in highly metastatic OC cells, while TGFβ2 was highly expressed in highly metastatic OC cells. In conclusion, miR-7-5p directly bound to TGFβ2 3′UTR and negatively regulated its expression.
Overexpression of miR-7-5p inhibits OC cell viability, migration and invasion by targeting TGFβ2
After determining the negative correlation between miR-7-5p and TGFβ2, we shifted to validate their roles functionally in malignant phenotypes of highly metastatic OC cell line HO-8910PM by using CCK-8 assay (Figs. 7a, 8b), Scratch test (Fig. 7b), Transwell invasion assay (Figs. 7c, 8c), flow cytometric detection (Fig. 7d), Annexin V-APC staining (Fig. 7e) and western blots (Fig. 7f). In the presence of miR-7-5p mimic, HO-8910PM cells presented with attenuated viability (Fig. 7a), migration (Fig. 7b) after 24 h (no significant difference observed at 8 h), and invasion (Fig. 7c). These results indicated that up-regulated miR-7-5p appreciably restricted the proliferative, migratory and invasive capacities of OC cells in vitro.
Moreover, the quantitative flow cytometric analysis (Fig. 7d) showed that the percentages of cells in G1 phase of the NC mimic group and the miR-7-5p mimic group were (63.62 ± 0.85)% and (68.32 ± 0.62)%, respectively; the percentages of cells in S phase of the NC mimic group and the miR-7-5p mimic group were (15.82 ± 0.79)% and (14.12 ± 0.52)%, respectively; and the percentages of cells in G2/M phase of the NC mimic group and the miR-7-5p mimic group were (20.56 ± 0.11)% and (17.57 ± 0.11)%, respectively. Accordingly, the miR-7-5p mimic treatment prolonged G1 phase (increased cell percentage) and shortened S and G2/M phases (reduced cell percentage). These data suggested that up-regulated miR-7-5p blocked OC cell cycle progression.
In terms of cell apoptosis detected by Annexin V-APC staining method (Fig. 7e), the apoptosis rate of HO-8910PM cells in the NC mimic group and miR-7-5p mimic group was lower than 5%, and there was no significant difference in the apoptosis rate between the two groups. Western blot assay further verified the effect of miR-7-5p mimic on cyclin D1 protein expression and caspase-3 activity in HO-8910PM cells. The results showed that overexpression of miR-7-5p up-regulated cyclin D1 protein expression, but had no effect on caspase-3 (Fig. 7f). Moreover, as shown in Fig. 7g, HO-8910PM cells transfected with miR-7-5p mimic presented with elevated expression of E-cadherin and lowered expression of Vimentin. This result indicated that overexpression of miR-7-5p suppressed epithelial–mesenchymal transition (EMT) of highly metastatic OC cell line HO-8910PM.
Additionally, we conducted rescue experiments (Fig. 8) on HO-8910PM cells overexpressing miR-7-5p and TGFβ2 simultaneously. Firstly, western blot analysis (Fig. 8a) suggested that expression of TGFβ2 protein was elevated in response to oe-TGFβ2 alone and reduced to miR-7-5p mimic alone, while the combined treatment of oe-TGFβ2 and miR-7-5p mimic rescued its expression. Furthermore, miR-7-5p mimic alone diminished viability and invasion of HO-8910PM cells, while TGFβ2 overexpression alone facilitated viability and invasion of HO-8910PM cells. However, simultaneous miR-7-5p and TGFβ2 overexpression rescued cell viability and invasion in response to miR-7-5p mimic or TGFβ2 overexpression alone (Fig. 8b, c).
MiRNAs confer pivotal roles in multiple cellular functions, such as proliferation, apoptosis and differentiation. However, aberrant expression of miRNAs appears in carcinogenesis and metastasis during cancer progression . Mounting evidence highlights a crucial role of miRNAs in OC progression. For instance, alterations in miR-101, miR-206, miR-200a and miR-203 are involved with OC cell proliferation and invasion [17,18,19]. Therefore, this current study presented differentially expressed miRNAs and mRNAs related to metastatic OC and revealed a tumor-inhibiting role of miR-7-5p in the restriction of OC metastasis, which was associated with the regulation of TGFβ2.
Different types of cancer samples show varying gene expression patterns, and even the same miRNA often regulate various target genes in different types of cancers. In studies on OC, the application of a miRNA expression profiling technique is helpful in identifying differentially expressed miRNAs affecting malignant phenotypes and metastasis in OC in a more rapid and comprehensive manner. In the present study, we analyzed the miRNA expression profiles of primary OC tissues and metastatic tissues using a miRStar™ Human Cancer Focus miRNA and Target mRNA PCR Array expression profiling. With the application of a miRCURY LNA™ microRNA Array expression profiling technique, Liu et al.  identified 31 differentially expressed miRNAs from another two pairs of cisplatin-resistant OC cell lines, of which the expression of 21 miRNAs was up-regulated and that of 10 miRNAs was down-regulated. Moreover, Cheng et al.  found by using an Affymetrix miRNA 3.0 Array expression profiling that 37 miRNAs were differentially expressed in adult and juvenile granulosa cell tumors, and that these could be used as markers for the diagnosis and recurrence of ovarian granulosa cell tumors.
The high-throughput sequencing screening of this study found 12 differentially expressed miRNAs (miR-195-5p, miR-27a-3p, miR-16-5p, miR-15a-5p, miR-10a-5p, miR-584-5p, miR-7-5p, miR-187-5p, miR-141-3p, miR-200b-3p, miR-200a-3p and miR-9-3p) and 10 differentially expressed mRNAs (BCL2L1, TGFβ2, TNC, HOXB5, HOXB7, HOXB8, ITGB3, MTSS1, PAK1 and FOXA2) in OC. GO analysis and Pathway analysis finally determined that 3 differentially expressed mRNAs (ITGB3, TGFβ2 and TNC) were related to OC metastasis. The results of RT-qPCR confirmed that the decrease of miR-7-5p was most significant in OC metastases, while TGFβ2 was significantly up-regulated in OC metastases. Therefore, miR-7-5p and its downstream target gene TGFβ2 were selected as research focuses. A wide array of miRNAs have been reported to present differential expression in OC . For instance, miR-187 has been found to exert a dual role in OC by regulating the disabled homolog-2 gene . Gao et al.  also found that miR-141 acted as a potential diagnostic and prognostic biomarker for OC, which was consistent with our in silico prediction results. Of interest, underexpressed miR-7-5p has been documented in metastatic breast cancer  and invasive pancreatic cancer . Another study also claimed that miR-7-5p deficiency was associated with recurrence in glioblastoma patients, and that its overexpression decreased glioblastoma cell stemness . Furthermore, it was also elucidated that ectopic expression of miR-7 functioned as an anti-oncogene in OC by repressing cell invasion and proliferation .
Recently, a growing number of studies elucidated that miRNAs are implicated in the invasive and metastatic potential of tumors. One miRNA can regulate nearly 100 different mRNAs, and a single mRNA can bind to several different miRNAs, thus forming a miRNA-target gene network regulatory relationship [28, 29]. Of note, miRNAs can post-transcriptionally mediate a number of genes through the binding of specific sequences in target mRNA molecules . TGFβ2 as a putative target of miR-7-5p was verified using a luciferase activity assay of the present study. Moreover, there was a negative correlation between miR-7-5p expression and TGFβ2 expression in OC tissues, and miR-7-5p targeted and negatively regulated TGFβ2. MiR-7-5p has been found to target and regulate target genes as SATB1 and PARP1 in some malignancies [31, 32]. For instance, Zhu et al.  demonstrated that miR-7-5p suppresses cell migratory and invasive potential by targeting SOX18 in pancreatic ductal adenocarcinoma.
In addition, through gain- and loss-of-function experiments in highly metastatic OC cells, our results revealed that miR-7-5p overexpression reduced OC cell viability and invasion, and blocked cell cycle entry by targeting TGFβ2. In addition, we found that OC metastatic tissues presented higher levels of TGFβ2 relative to primary OC tissues. TGFβ2 is an isoform of TGFβ in mammals  and is abnormally expressed in various cancers, such as human melanoma and hepatocellular carcinoma [35, 36]. Furthermore, TGFβ2 and the TGFβ pathway could interact with miRNAs, such as miR-187 and miR-141 identified in this study, to affect growth of tumors. Ectopic miR-187 expression could negate the contribution of TGFβ to aggressiveness in metastatic colorectal cancer . Besides, the interplay of miR-141 and oncogenic TGFβ2 was found to orchestrate malignancy in gastric cancer . Additionally, many researchers have also explored the role of TGFβ in OC. Cao et al.  reported that TGFβ-induced transglutaminase accelerated EMT and a cancer stem cell phenotype that consequently enhanced ovarian tumor metastasis. TGFβ induced EMT and a more invasive phenotype in epithelial OC cells in collaboration with the EGF pathway, indicating TGFβ may be a promising target candidate for the treatment of metastatic OC in future .
In conclusion, our study suggested several differentially expressed miRNAs and mRNAs for metastatic OC, and uncovered the tumor-suppressive function of miR-7-5p in highly metastatic OC cells. miR-7-5p directly bound to TGFβ2 3′UTR to inhibit its expression, thus restricting the invasive and metastatic potential of OC. miR-7-5p may be harnessed as a potential target for preventive and therapeutic strategies against metastatic OC. Although bioinformatic analyses and in vitro findings shed new light on the mechanistic understanding of tumor-suppressive miR-7-5p and its target TGFβ2 in OC metastasis, the conclusions of this study need to be further confirmed in future studies based on in vivo experimental validations. Moreover, it is necessary to conduct further cohort studies with larger sample sizes to validate our results.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
- miR or miRNA:
Reverse transcription quantitative polymerase chain reaction
Kyoto Encyclopedia of Genes and Genomes
Fetal bovine serum
Cell counting kit 8
Phosphate buffer saline
Moufarrij S, Dandapani M, Arthofer E, Gomez S, Srivastava A, Lopez-Acevedo M, Villagra A, Chiappinelli KB. Epigenetic therapy for ovarian cancer: promise and progress. Clin Epigenet. 2019;11(1):7.
Franier B, Thompson M. Early stage detection and screening of ovarian cancer: a research opportunity and significant challenge for biosensor technology. Biosens Bioelectron. 2019;135:71–81.
Emmings E, Mullany S, Chang Z, Landen CN Jr, Linder S, Bazzaro M. Targeting mitochondria for treatment of chemoresistant ovarian cancer. Int J Mol Sci. 2019;20(1):229.
Hoter A, Naim HY. Heat shock proteins and ovarian cancer: important roles and therapeutic opportunities. Cancers. 2019;11(9):1389.
Riester M, Wei W, Waldron L, Culhane AC, Trippa L, Oliva E, et al. Risk prediction for late-stage ovarian cancer by meta-analysis of 1525 patient samples. J Natl Cancer Inst. 2014;106(5):dju048.
Zhang H, Lu B. microRNAs as biomarkers of ovarian cancer. Expert Rev Anticancer Ther. 2020;20(5):373–85.
Gebert LFR, MacRae IJ. Regulation of microRNA function in animals. Nat Rev Mol Cell Biol. 2019;20(1):21–37.
Chong GO, Jeon HS, Han HS, Son JW, Lee YH, Hong DG, Lee YS, Cho YL. Differential MicroRNA expression profiles in primary and recurrent epithelial ovarian cancer. Anticancer Res. 2015;35(5):2611–7.
Vang S, Wu HT, Fischer A, Miller DH, MacLaughlan S, Douglass E, Comisar L, Steinhoff M, Collins C, Smith PJ, et al. Identification of ovarian cancer metastatic miRNAs. PLoS ONE. 2013;8(3):e58226.
Chen SN, Chang R, Lin LT, Chern CU, Tsai HW, Wen ZH, Li YH, Li CJ, Tsui KH. MicroRNA in ovarian cancer: biology, pathogenesis, and therapeutic opportunities. Int J Environ Res Public Health. 2019;16(9):1510.
Cui X, Sun Y, Shen M, Song K, Yin X, Di W, Duan Y. Enhanced chemotherapeutic efficacy of paclitaxel nanoparticles co-delivered with microRNA-7 by inhibiting paclitaxel-induced EGFR/ERK pathway activation for ovarian cancer therapy. ACS Appl Mater Interfaces. 2018;10(9):7821–31.
Liu M, Zhang X, Hu CF, Xu Q, Zhu HX, Xu NZ. MicroRNA-mRNA functional pairs for cisplatin resistance in ovarian cancer cells. Chin J Cancer. 2014;33(6):285–94.
Dou YD, Huang T, Wang Q, Shu X, Zhao SG, Li L, et al. Integrated microRNA and mRNA signatures in peripheral blood lymphocytes of familial epithelial ovarian cancer. Biochem Biophys Res Commun. 2018;496(1):191–8.
Ruibin J, Guoping C, Zhiguo Z, Maowei N, Danying W, Jianguo F, et al. Establishment and characterization of a highly metastatic ovarian cancer cell line. Biomed Res Int. 2018;2018:3972534.
Dong PX, Jia N, Xu ZJ, Liu YT, Li DJ, Feng YJ. Silencing of IQGAP1 by shRNA inhibits the invasion of ovarian carcinoma HO-8910PM cells in vitro. J Exp Clin Cancer Res. 2008;27:77.
Lin CW, Chang YL, Chang YC, Lin JC, Chen CC, Pan SH, Wu CT, Chen HY, Yang SC, Hong TM, et al. MicroRNA-135b promotes lung cancer metastasis by regulating multiple targets in the Hippo pathway and LZTS1. Nat Commun. 2013;4:1877.
Semaan A, Qazi AM, Seward S, Chamala S, Bryant CS, Kumar S, Morris R, Steffes CP, Bouwman DL, Munkarah AR, et al. MicroRNA-101 inhibits growth of epithelial ovarian cancer by relieving chromatin-mediated transcriptional repression of p21(waf(1)/cip(1)). Pharm Res. 2011;28(12):3079–90.
Li S, Li Y, Wen Z, Kong F, Guan X, Liu W. microRNA-206 overexpression inhibits cellular proliferation and invasion of estrogen receptor alpha-positive ovarian cancer cells. Mol Med Rep. 2014;9(5):1703–8.
Wang S, Zhao X, Wang J, Wen Y, Zhang L, Wang D, Chen H, Chen Q, Xiang W. Upregulation of microRNA-203 is associated with advanced tumor progression and poor prognosis in epithelial ovarian cancer. Med Oncol. 2013;30(3):681.
Liu MX, Siu MK, Liu SS, Yam JW, Ngan HY, Chan DW. Epigenetic silencing of microRNA-199b-5p is associated with acquired chemoresistance via activation of JAG1-Notch1 signaling in ovarian cancer. Oncotarget. 2014;5(4):944–58.
Cheng WT, Rosario R, Muthukaruppan A, Wilson MK, Payne K, Fong PC, Shelling AN, Blenkiron C. MicroRNA profiling of ovarian granulosa cell tumours reveals novel diagnostic and prognostic markers. Clin Epigenetics. 2017;9:72.
Taylor DD, Gercel-Taylor C. MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol. 2008;110(1):13–21.
Chao A, Lin CY, Lee YS, Tsai CL, Wei PC, Hsueh S, Wu TI, Tsai CN, Wang CJ, Chao AS, et al. Regulation of ovarian cancer progression by microRNA-187 through targeting disabled homolog-2. Oncogene. 2012;31(6):764–75.
Gao YC, Wu J. MicroRNA-200c and microRNA-141 as potential diagnostic and prognostic biomarkers for ovarian cancer. Tumour Biol. 2015;36(6):4843–50.
Okuda H, Xing F, Pandey PR, Sharma S, Watabe M, Pai SK, et al. miR-7 suppresses brain metastasis of breast cancer stem-like cells by modulating KLF4. Cancer Res. 2013;73(4):1434–44.
Ma J, Fang B, Zeng F, Pang H, Zhang J, Shi Y, et al. Curcumin inhibits cell growth and invasion through up-regulation of miR-7 in pancreatic cancer cells. Toxicol Lett. 2014;231(1):82–91.
Jia B, Liu W, Gu J, Wang J, Lv W, Zhang W, Hao Q, Pang Z, Mu N, Zhang W, et al. MiR-7-5p suppresses stemness and enhances temozolomide sensitivity of drug-resistant glioblastoma cells by targeting Yin Yang 1. Exp Cell Res. 2019;375(1):73–81.
Shin C, Nam JW, Farh KK, Chiang HR, Shkumatava A, Bartel DP. Expanding the microRNA targeting code: functional sites with centered pairing. Mol Cell. 2010;38(6):789–802.
de Giorgio A, Krell J, Harding V, Stebbing J, Castellano L. Emerging roles of competing endogenous RNAs in cancer: insights from the regulation of PTEN. Mol Cell Biol. 2013;33(20):3976–82.
Elton TS, Yalowich JC. Experimental procedures to identify and validate specific mRNA targets of miRNAs. EXCLI J. 2015;14:758–90.
Yin CY, Kong W, Jiang J, Xu H, Zhao W. miR-7-5p inhibits cell migration and invasion in glioblastoma through targeting SATB1. Oncol Lett. 2019;17(2):1819–25.
Lai J, Yang H, Zhu Y, Ruan M, Huang Y, Zhang Q. MiR-7-5p-mediated downregulation of PARP1 impacts DNA homologous recombination repair and resistance to doxorubicin in small cell lung cancer. BMC Cancer. 2019;19(1):602.
Zhu W, Wang Y, Zhang D, Yu X, Leng X. MiR-7-5p functions as a tumor suppressor by targeting SOX18 in pancreatic ductal adenocarcinoma. Biochem Biophys Res Commun. 2018;497(4):963–70.
Xie B, Zhang C, Kang K, Jiang S. miR-599 inhibits vascular smooth muscle cells proliferation and migration by targeting TGFB2. PLoS ONE. 2015;10(11):e0141512.
Li JR, Wang JQ, Gong Q, Fang RH, Guo YL. MicroRNA-328 inhibits proliferation of human melanoma cells by targeting TGFbeta2. Asian Pac J Cancer Prev. 2015;16(4):1575–9.
Dropmann A, Dediulia T, Breitkopf-Heinlein K, Korhonen H, Janicot M, Weber SN, Thomas M, Piiper A, Bertran E, Fabregat I, et al. TGF-beta1 and TGF-beta2 abundance in liver diseases of mice and men. Oncotarget. 2016;7(15):19499–518.
Zhang F, Luo Y, Shao Z, Xu L, Liu X, Niu Y, et al. MicroRNA-187, a downstream effector of TGFbeta pathway, suppresses Smad-mediated epithelial–mesenchymal transition in colorectal cancer. Cancer Lett. 2016;373(2):203–13.
Lei K, Liang X, Gao Y, Xu B, Xu Y, Li Y, et al. Lnc-ATB contributes to gastric cancer growth through a MiR-141-3p/TGFbeta2 feedback loop. Biochem Biophys Res Commun. 2017;484(3):514–21.
Cao L, Shao M, Schilder J, Guise T, Mohammad KS, Matei D. Tissue transglutaminase links TGF-beta, epithelial to mesenchymal transition and a stem cell phenotype in ovarian cancer. Oncogene. 2012;31(20):2521–34.
Xu Z, Jiang Y, Steed H, Davidge S, Fu Y. TGFbeta and EGF synergistically induce a more invasive phenotype of epithelial ovarian cancer cells. Biochem Biophys Res Commun. 2010;401(3):376–81.
We acknowledge and appreciate our colleagues for their valuable efforts and comments on this paper.
This study was supported by the Free Researcher Project of Shengjing Hospital (Grant No. 201302).
Ethics approval and consent to participate
Data and tissues were harvested upon receiving the informed consent of patients and approval by the ethics committee of Shengjing Hospital of China Medical University (Approval number: 2016PS040K).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gu, Y., Zhang, S. High-throughput sequencing identification of differentially expressed microRNAs in metastatic ovarian cancer with experimental validations. Cancer Cell Int 20, 517 (2020). https://doi.org/10.1186/s12935-020-01601-4