- Open Access
HOXD-AS1 promotes the epithelial to mesenchymal transition of ovarian cancer cells by regulating miR-186-5p and PIK3R3
Journal of Experimental & Clinical Cancer Researchvolume 38, Article number: 110 (2019)
Epithelial ovarian cancer (EOC) is one of the most malignant gynecological tumors worldwide. Deregulation of long non-coding RNAs (lncRNAs) has been implicated in various oncogenic processes in multiple cancers. In this study, we aim to identify and characterize clinically relevant lncRNA deregulation in EOC.
LncRNAs, mRNAs and miRNAs were profiled using expression microarrays and validated using reverse transcription quantitative PCR in EOC cells and tissues. siRNAs targeting either HOXD-AS1 or PIK3R3 together with miR-186-5p inhibitors were used to modulate endogenous target expression in EOC cell lines in vitro. In vitro wound healing assay, trans-well assay, Western-blot assay,and Dual-luciferase reporter assay were used to explore the biological roles and molecular function underlying HOXD-AS1 in the EOC cells. Progression-free survival (PFS) and overall survival (OS) were statistically analyzed by Kaplan-Meier method test.
HOXD-AS1 was found to be significantly over-expressed in EOC tumors. High HOXD-AS1 expression significantly correlated with poorer PFS and OS of EOC patients. Multivariate Cox proportional hazards modeling indicated that HOXD-AS1 was an independent risk predictor of EOC patients (HR = 1.92, p = 0.004). SiRNA inhibition of HOXD-AS1 reduced cell migration, invasion, and epithelial-mesenchymal transition (EMT) in EOC cells in vitro by preventing HOXD-AS1 directly binding to miR-186-5p, and resulting in down-regulating of PIK3R3. The novel HOXD-AS1/miR-186-5p/PIK3R3 pathway was clinically relevant as we observed a significantly inverse correlation between HOXD-AS1/miR-186-5p and between miR-186-5p/PIK3R3 in an independent cohort of 200 EOC tissues.
HOXD-AS1/miR-186-5p/PIK3R3 is a novel pathway to promote cell migration, invasion, and EMT in EOC.
EOC is the third most common and second most lethal gynecological malignancy worldwide . The standard treatment is maximal cytoreductive surgery followed by platinum-based chemotherapy regimen . However, the occurrence of metastasis and acquired drug resistance often lead to poor prognosis in patients with advanced ovarian cancer [3, 4]. Therefore, there is an urgent clinical need to understand the molecular mechanisms leading to the occurrence of tumor metastasis and drug resistance in EOC. In addition, novel biomarkers which can accurately predict tumor metastasis and drug resistance are needed to provide valuable insight into the clinical management of EOC.
EMT is a biological process in which polarized epithelial cells adopt mesenchymal cell phenotype which includes enhanced migratory capacity, invasiveness, elevated resistance to apoptosis, and greatly increased production of extracellular matrix components [5, 6]. It has been well demonstrated that the activation of EMT program endows cancer cells with the multiple features that are required to execute most steps of the invasion-metastasis cascade [7, 8]. In addition, activation of EMT pathway induces chemoresistance in multiple cancers such as pancreatic cancer , lung cancer [10, 11], gastric cancer , liver cancer  and EOC .
LncRNA and miRNAs both belong to non-protein encoding RNAs which are implicated to play important roles in a variety of cellular processes such as cell proliferation, apoptosis, and migration [15, 16]. The aberrant expression of lncRNAs in tumors can induce EMT and increase the invasion and migration of tumor cells. For example, XIST, which is a lncRNA over-expressed in colorectal cancer tissues and cells, can promote invasion, proliferation, and EMT progress through XIST/miR-200b-3p/ZEB1 pathway . Another example is PVT1, a highly expressed lncRNA in prostate cancer cells that can co-regulate Twist with miR-186 and promote tumor invasion, migration, and the progress of EMT . This lncRNA/miRNA/mRNA pathway has been proposed as a new model by which lncRNAs regulate target gene expression through miRNA response elements . Salmena et al. have termed this phenomenon as the competing endogenous RNA (ceRNA) theory .
In this study, we identified a HOXD-AS1/miR-186-5p/PIK3R3 pathway in EOC as another example exhibiting the ceRNA property. The HOXD-AS1/miR-186-5p/PIK3R3 pathway plays an important role in invasion, migration and EMT of EOC.
Patients and EOC tissue samples
We collected EOC (n = 36) and normal ovarian (n = 14) tissue samples from patients who underwent surgery at the Hunan Cancer Hospital (Affiliated Cancer Hospital of Xiangya School of Medicine, Central South University, Changsha, Hunan, China) from 2012 to 2014. The resected tissues were separately identified and confirmed by two pathologists and kept in liquid nitrogen immediately for subsequent use. In addition we also collected formalin-fixed paraffin-embedded (FFPE) specimens from an independent cohort of 200 patients to analyze the relationship of PFS and OS with the expression levels of HOXD-AS1. The clinicopathologic characteristics of patients are shown in Table 1. The ethics committee of Hunan Cancer Hospital approved the protocol, and all patients signed informed consents before sample collection.
Differentially expressed gene analysis
Differentially expressed genes and lncRNAs were identified using the threshold of fold change ≥2 and Student’s t-test p-value < 0.01 between EOC tumor versus normal tissues. The R package, Significance analysis of microarray (SAM) was used to identify differentially expressed miRNA at a threshold of fold change≥1.5, false discovery rate (FDR) q-value< 0.05. Hierarchical clustering was performed on the pheatmap package of R language and heatmap was constructed. Volcano maps were drawn using the ggplot2 packages.
Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analysis
DAVID (The Database for Annotation, Visualization and Integrated Discovery, https://david.ncifcrf.gov/) was used to annotate the function of differentially expressed mRNA and enrichment of the signaling pathway. GO included cell component, biological process, and molecular function. Aberrant mRNAs were enriched into different pathways according to the KEGG database. TOP 10 terms were presented. The recommended p-value cutoff of p < 0.05 was used.
Target gene prediction
The potential miRNAs that may bind to HOXD-AS1 were predicted by DIANA (http://carolina.imis.athena-innovation.gr/diana_tools/web/index.php?r=lncbasev2%2Findex-predicted). The candidate mRNA that may bind to miR-186-5p were predicted by starBase v2.0 (http://starbase.sysu.edu.cn/) database, which integrates five programs: targetScan, picTar, RNA22, PITA, and miRanda/mirSVR. Candidate mRNA are defined as the genes which are predicted to be target mRNA of miR-186-5p by at least three programs.
EOC cell line A2780 was cultured in Dulbecco’s modified Eagle’s medium (TransGen, Beijing, China), and EOC cell line SKOV3 and normal ovarian cell line IOSE80 were cultured in RPMI-1640 (TransGen, Beijing, China) with 10% fetal bovine serum (Gibco) and 1% penicillin/streptomycin (TransGen, Beijing, China). All of the cells were cultured in an incubator at 37 °C with 5% CO2.
RNA extraction and RT-qPCR
Total RNAs from tissues and cells were extracted using Trizol reagent (Invitrogen). Total RNAs from FFPE were extracted using RNeasy®FFPE Kit (Qiagen, Germany). cDNA was synthesized with the TransScript® II One-Step gDNA Removal and cDNA Synthesis SuperMix Kit (TransGen, Beijing, China) according to the manufacturer’s instructions. RT-qPCR was performed on the LightCycler®96 system (Roche) using a standard protocol from the TransStart® Tip Green qPCRSuperMix Kit (TransGen, Beijing, China). The results were normalized to GAPDH. The sequences of primers used for RT-qPCR assays are listed in Additional file 1.
For the detection of miRNA expression, reverse transcription and qPCR were performed using stem-loop primers and Bulge-Loop™ miRNA RT-qPCR Starter Kit from Ribobio (Guangzhou, China). U6 snoRNA was used as the endogenous control. All experiments were performed in triplicate according to the manufacturer’s manual. Expression fold changes were calculated using 2-∆∆Ct methods.
In situ hybridization (ISH)
In situ hybridization was performed to detect the expression of HOXD-AS1, miR-186-5p, and PIK3R3. The probes of HOXD-AS1, miR-186-5p, and PIK3R3 were designed by BOSTER (Wuhan, China). RNA ISH was performed using the in situ hybridization detection kit (BOSTER) and DAB chromogenic kit (ZSGB-BIO) on 3-μm FFPE tissue sections according to the manufacturer’s instructions. The staining scores were determined based on both the intensity and proportion of target gene-positive cells under 20X and 40X objective.
Small interfering RNA (siRNA) and miRNA inhibitor transfections
LncRNA HOXD-AS1, siRNA (si-HOXD-AS1), PIK3R3 siRNA (si-PIK3R3), miR-186-5p inhibitor, and negative control (si-NC, miR-NC) were all obtained from Ribobio (Guangzhou, China). RiboFECT™ CP transfection Kit (Guangzhou, China) was used for cell transfection based on manufacturer’s instructions. SKOV3 and A2780 were seeded in a six-well plate and transfection would be performed when cell’s confluence reached 50%. Cells were transfected with a final concentration of 100 nM of individual siRNA and 200 nM of miR-186-5p inhibitor. All of the steps were carried out according to the manufacturer’s specification. Experiments were performed three times.
Wound healing assay
The cell monolayer was scraped in straight lines using a 10ul pipette tip and washed with PBS two times when transfected cells were up to 80% confluence in a six-well culture dish. The cells were cultured in culture medium with 3% FBS and 1% penicillin/streptomycin (TransGen, Beijing, China). Images were captured at 0, 12, 24, 36, and 48 h following the initial scratch from five different areas for each wound. We used Image J (National Institutes of Health, Bethesda, MD, USA) software to calculate cell wound healing rate as: (the original wound areas - the actual wound areas at different times)/ (the original wound areas).
The invasion ability was measured using transwell chambers coated with 1:30 diluted Matrigel Basement Membrane Matrix (Corning Life Sciences, NY, USA). Migration assays were performed without matrigel. 1 × 105 cells were seeded into upper chambers (8 μM pore size, Corning Life Sciences, NY, USA) with 200ul serum-free media. A total of 600 ul of 20% FBS media was added into the lower chambers. After incubation for 48 h at 37 °C, the cells in the upper chamber that did not invade through the pores were wiped out with cotton wool. The cells in the lower chamber were fixed with methanol and stained with 0.5% crystal violet. Cells were counted and imaged at × 200 magnification. The experiments were repeated three separate times.
After transfection for 48 h, cells were lysed using RIPA buffer (Thermo, USA) with proteinase inhibitor cocktail (Roche), then incubated at 4 °C for 15 min. The lysate was centrifuged and protein concentration was measured by BCA Protein Assay Kit (Thermo, USA). The protein was separated on 10% sodium dodecyl sulfate- polyacrylamide gel electrophoresis (SDS-PAGE) and then transferred to 0.22 μm PVDF membrane (Thermo, USA). After blocking with 5% fat-free milk for 1 h, the membranes were incubated with primary antibodies at 4 °C overnight with gentle shaking. Next, the membranes were incubated with secondary antibody at room temperature for 2 h. Finally, the immunofluorescence of protein bands was visualized using Pierce™ ECL Western Blotting Substrate (32,109, Thermo Scientific, USA) by Image J 2 software (Madison, WI, USA). Primary antibodies are as following: E-cadherin (1: 500, ab15148, Abcam company), Vimentin (1:500, #7431, Cell Signaling Technology (CST)), GAPDH (1:5000, 60,004-l-lg, Proteintech company), Rac1 (1:500, ab33186, Abcam), Twist (1500, GTX60776, Gene Tex).
Dual-luciferase reporter assay
The dual luciferase reporter plasmids (LncRNA-HOXD-AS1-WT, LncRNA-HOXD-AS1-Mut, PIK3R3-WT, PIK3R3-Mut, miR-186-5p, and negative control) were synthesized by GenePharma (Shanghai, China). The SKOV3 cells were seeded into 96-well plates and cultured to 80% confluence. The plasmid vectors or empty plasmid were transiently transfected into SKOV3 cells using lipofectamine 3000 (Invitrogen). After 48 h incubation, the luciferase activity was examined by Dual-Glo® Luciferase Assay System (Promega, USA) according to the manufacturer’s protocol.
All experiments were performed three separate times. Statistical analysis was performed with SPSS 13.0 and graphs were constructed with GraphPad Prism 5 (La Jolla, CA, USA). The data were presented as the mean ± Standard Deviation (SD). The significance of differences between groups was assessed by Student’s t-test and Chi-square test. The Kaplan-Meier method test was used for PFS analysis and OS analysis. Pearson correlation analysis was used to determine the relationship between two genes. We used the univariate and multivariate Cox proportional hazards modeling to determine the effects of variables on survival. A p-value of < 0.05 was considered statistically significant.
Transcriptome landscape of EOC
We performed transcriptome profiling on 6 human epithelial ovarian cancer tissues and three matched normal ovary samples using the Agilent Human lncRNArray + mRNA Array V4.0. at a threshold of absolute fold change of greater than 2 and p-value of less than 0.01. We observed a total of 2552 mRNAs in EOC to be differentially expressed compared to normal ovary samples, of which 1496 were up-regulated, and 1056 were down-regulated (Fig. 1a and b). These differentially expressed mRNAs were enriched in a number of KEGG pathways such as Straphylococcus aureus infection and cell adhesion molecules which are relevant pathways linking EOC for either increased risk or aggressiveness of the cancer (Fig. 1c). In the 1922 well-annotated lncRNAs on the array, we observed 288 lncRNAs to be differentially expressed (Fig. 1d). Among which, 51 were significantly up-regulated and 237 were significantly down-regulated in EOC tumors versus the normal ovary tissues. These 288 lncRNAs can stratify the samples into distinct clusters based on expression values (Fig. 1e). We validated the robustness of the array data by randomly choosing ten lncRNAs from these 288 differentially expressed lncRNAs in an independent cohort of 50 samples comprising of 36 EOC tissues and 14 normal ovary tissues. All ten lncRNAs showed differential expression between EOC and normal ovary tissues which was in concordance with what we observed in the microarrays. In conclusion, the transcriptome profile constructed from high quality data using the Agilent Human lncRNArray+mRNA Array were reproducible in an independent cohort. Table 1 summarized the top ten most up-regulated and down-regulated lncRNAs in EOC tumors and we selected HOXD-AS1 for further characterization.
HOXD-AS1 is highly expressed in EOC tissues which negatively impacts patients survival
We performed RT-qPCR to validate HOXD-AS1 up-regulation in an independent cohort comprising of 36 EOC tissues and 14 normal ovarian tissues and found HOXD-AS1 to be consistently and significantly up-regulated in EOC tissues (p < 0.05, Fig. 2a). We further examined HOXD-AS1 expression in another independent cohort of 200 EOC FFPE samples and correlated HOXD-AS1 expression with the clinicopathological features of EOC patients. As shown in Table 2, high HOXD-AS1 expression (above the median value = 0.000565) was positively associated more aggressive EOC tumors characterized by higher tumor grade (Grade G2/3), advanced stage (FIGO stage III/IV), and the presence of lymph node metastasis (p < 0.05). We observed that high HOXD-AS1 level was significantly associated with poorer PFS and OS of EOC patients through Kaplan-Meier analysis (p < 0.05, Fig. 2b and c). Univariate Cox regression analysis showed that the poor prognosis was observed in EOC patients with high-grade tumor, FIGO stage III / IV, lymph node metastasis, and high HOXD-AS1 expression. Multivariate Cox regression analysis suggested that tumor grade, FIGO stage, lymph node metastasis status, and HOXD-AS1 expression are independent risk predictors of prognosis for EOC patients (Table 3). Our data suggested that HOXD-AS1 may play an important functional role in EOC tumor.
HOXD-AS1 promotes EOC cell migration, invasion, and EMT transition
To understand the molecular function of HOXD-AS1 on EOC cells, we designed HOXD-AS1-specific siRNAs to inhibit the endogenous HOXD-AS1 expression in SKOV3 and A2780 EOC cell lines, both of which showed higher HOXD-AS1 expression than normal ovarian cell line IOSE80 (Fig. 3a). SKOV3 and A2780 cells showed significantly reduced ability to migrate in wound healing assays (Fig. 3c and d) when HOXD-AS1 was knocked-down (Fig. 3b). EOC cells in trans-well assays transfected with siRNA against HOXD-AS1 consistently demonstrated significantly reduced capacity to migrate or invade when compared to the mock-transfected control cells (Fig. 3e-h). EMT process is closely related to the invasion and migration of cancer cells therefore we examined the EMT markers between the HOXD-AS1 knocked-down cells and control cells. As shown in Fig. 3i, we observed a significant increase in the expression of epithelial marker E-cadherin and a corresponding decrease in mesenchymal marker vimentin in the HOXD-AS1 knock-down cells as compared to controlls. These results suggested that the down-regulation of HOXD-AS1 inhibited the invasion, migration, and EMT process of EOC cells in vitro.
HOXD-AS1 binds miR-186-5p to promote cell migration, invasion, and EMT
We performed Affymetrix miRNA 4.0 chips to evaluate if HOXD-AS1 can increases EOC cell migration, invasion, and EMT through its interaction with miRNAs in a ceRNA manner. We identified a total of 164 miRNAs to be significantly down-regulated in the same nine samples profiled for miRNAs of which 73 were found to be consistently reported by an independent GEO dataset (GSE53829). We further shortlisted a total of 6 miRNAs by overlapping these 73 down-regulated miRNAs in EOC tissues with 293 miRNAs predicted to interact with HOXD-AS1 by DIANA algorithm (Fig. 4a and Table 4). The increased miR-186-5p expression amongst these 6 miRNAs were confirmed upon HOXD-AS1 inhibition (Fig. 4b). We selected miR-186-5p for further research because miR-186-5p consistently showed an inverse correlation with HOXD-AS1 in both patient samples and cell lines. HOXD-AS1 with 6 putative miR-186-5p binding sites further enforced our interest. Inhibition of miR-186-5p promote cell migration, invasion and EMT in EOC cells (Additional file 2: Figure S1).
To demonstrate that HOXD-AS1 interacts with miR-186-5p through its putative miR-186-5p binding sites, we cloned the wildtype and a mutant HOXD-AS1 in which all six putative miR-186-5p binding sites were mutated and inserted downstream of a firefly luciferase gene (Fig. 4c). As shown in Fig. 4d, we observed significantly reduced reporter activity in the wildtype HOXD-AS1 construct when the cells were co-transfected with miR-186-5p compared to wildtype HOXD-AS1 construct co-transfected with the miRNA controls. However, such difference was abrogated when the putative miR-186-5p binding sites were mutated, indicating that HOXD-AS1 physically interacts with miR-186-5p at its putative binding sites to regulate reporter gene activity. It is further evidenced by concurrent increase in miR-186-5p expression when EOC cells were transfected with siRNAs targeting HOXD-AS1. The HOXD-AS1 knocked-down cells exhibited more epithelial and less mesenchymal phenotype (Fig. 4e, f, g and h, middle panel) which lead to reduced ability to migrate or invade. We observed a corresponding reversal of the above phenotype in cell migration, invasion, and EMT (Fig. 4e and f right panel) when miR-186-5p inhibitors were co-transfected with si-HOXD-AS1 to partly negate the increase in miR-186-5p expression. Therefore, our data demonstrated HOXD-AS1 promotes cell migration, invasion, and EMT through inhibiting miR-186-5p.
miR-186-5p targets PIK3R3 to negatively regulate cell migration, invasion, and EMT
In order to investigate how miRNAs regulate cellular functions through its target genes we queried starBase v2.0 to identify a total of 284 predicted targets of miR-186-5p, among which 33 were significantly up-regulated in EOC tissues with low miR-186-5p expression (Fig. 5a). KEGG pathway enrichment analysis identified four pathways such as focal adhesion all are important for cell migration and invasion. PIK3R3 was involved in all four pathways which suggestes PIK3R3 might be a direct miR-186-5p target. To test this hypothesis, we cloned the wildtype 3′ untranslated region (3’UTR) of PIK3R3 and inserted into the downstream region of the luciferase reporter gene. We mutated the two putative miR-186-5p binding sites along the PIK3R3 3’UTR to construct a mutant clone (Fig. 5b). Specific reduction in luciferase activity was only observed in EOC cells co-transfected with miR-186-5p and wildtype PIK3R3 3’UTR but not the mutant PIK3R3 3’UTR, confirming our hypothesis that miR-186-5p interacts with putative binding sites along PIK3R3 3’UTR to down-regulated luciferase reporter gene expression (Fig. 5c). Furthermore, EOC cells transfected with siRNA against PIK3R3 exhibited greater epithelial phenotypes and reduced ability to migrate and invade. This effect could be rescued by the co-transfection of miR-186-5p inhibitors as our data suggested that miR-186-5p binds to PIK3R3 to negatively regulate cell migration, invasion, and EMT (Fig. 5d, e, f & g). As shown in Fig. 5h, in EOC cells transfected with siRNA against HOXD-AS1, we observed increased miR-186-5p and decreased PIK3R3 along with Twist and Rac1 which contributed to reduced cell migration, invasion and EMT.
HOXD-AS1/miR-186-5p/PIK3R3 is a clinically relevant pathway regulating EMT in EOC tissues
We proceeded to evaluate the clinical relevance of the in vitro HOXD-AS1/miR-186-5p/PIK3R3 pathway in 200 FFPE EOC samples. We observed significant inverse correlations between HOXD-AS1/miR-186-5-p (Fig. 6a) and between miR-186-5p/PIK3R3 (Fig. 6b) and significant positive correlation between HOXD-AS1 and PIK3R3 (Fig. 6c) using RT-qPCR quantitation. Additional in situ hybridization (ISH) data were shown to be consistent with RT-qPCR data (Fig. 6d). Furthermore, patients with high HOXD-AS1 (above the median value = 0.000565), low miR-186-5p (below the median value = 0.00000413), and high PIK3R3 (above the median value = 0.001908) expression showed significantly poorer PFS and OS compared to patients with low HOXD-AS1, high miR-186-5p, and low PIK3R3 (Fig. 6e and f). As summarized in Fig. 6g, in normal ovary cells, banlanced HOXD-AS1and miR-186-5p expression keeps EMT in check by regulating PIK3R3. In EOC cells, significantly elevated HOXD-AS1 level suppresses cellular miR-186-5p leading to increased PIK3R3 level which in turn leads to EMT.
Tumor metastasis is a major factor contributing to the poor prognosis of EOC patients [21, 22]. EMT is a crucial cellular process associated with metastasis in tumor cells [23, 24]. Our group has previously demonstrated that platinum-based drugs could induce the activation of Rac1 and initiate EMT . Understanding the crucial players in the EMT process in EOC will provide insights into the development of new therapeutic strategies to treat the disease.
LncRNAs are an class of non-coding RNAs which plays an important regulatory role in the biological processes of tumorigenesis, tumor cell proliferation, invasion, and migration . In EOC, aberrant expression of lncRNA regulates target genes through different mechanisms affecting biological behaviour of tumor cells. For example, lncRNA ABHD11-AS1 through regulating RhoC can promote A2780 and OVCAR3 proliferation and inhibit apoptosis ; lncRNA PTAR through miR-101 -3p/ZEB1 pathway behave in ceRNA model promotes invasion, migration, and EMT in EOC ; lncRNA PANDAR through SFRS2/p53 phosphorylation involved in the acquired Cisplatin resistance . lncRNA aberrant expression in EOC can affect tumor proliferation, migration, and drug sensitivity. It also plays an important regulatory role in cancer initiation and progression.
lncRNAs have been divided into 5 categories based on the location of the lncRNA gene in the genome: (1) sense lncRNA, (2) antisense lncRNA, (3) bidirectional lncRNA, (4) intronic lncRNA, and (5) intergenic lncRNA . In this study, we systematically investigated lncRNAs deregulation in EOC tissues versus normal ovary tissues and identified HOXD-AS1 to be significantly up-regulated in EOC tumors. HOXD-AS1 is an antisense lncRNA located on chromosome 2q31.2. Lu et al. have previously shown that HOXD-AS1 is over-expressed in metastatic hepatocellular carcinoma and promotes hepatocellular carcinoma metastasis through miR19a/ARHGAP11A signaling pathway . Furthermore HOXD-AS1 has been demonstrated to cooperate with miR-130a to regulate SOX4 and E2F8, and promote metastasis in hepatocellular carcinoma and glioma respectively [31, 32]. In this study, we demonstrated that HOXD-AS1 acting as a cellular sponge in a ceRNA mechanism inhibits miR-186-5p and resulted in the corresponding up-regulation of PIK3R3 in EOC cells. This novel HOXD-AS1/miR-186-5p/PIK3R3 pathway plays an important role in EMT in EOC cells. In addition, high HOXD-AS1 level is significantly associated with poorer PFS and OS of EOC patients (Figs. 2 & 6), and hence demonstrating the clinical relevance of this novel pathway.
MiRNA-186-5p can bind to HOXD-AS1 and affect the expression of downstream target genes in the form of ceRNA. Abnormal expression of miRNA-186-5p in various tumors plays a regulatory role in the process of EMT. For example, high expression of miR-186-5p in colorectal cancer cells can inhibit the proliferation, metastasis, and EMT process of tumor cells; Expression changes of miR-186-5p are closely related to the proliferation and invasion ability of non-small cell lung cancer . Therefore, the study further focused on whether HOXD-AS1 can regulate the expression of miR-186-5p through a ceRNA sponge model. Data of the present study suggests that HOXD-AS1 can negatively regulate miR-186-5p.
PIK3R3 is a regulatory subunit of PI3K, whose differential expression can promote the proliferation, invasion, and migration process of NSCLC cells, cervical carcinoma cells, and pancreatic cancer cells. Our studies have suggested that PIK3R3 is the target gene regulated by HOXD-AS1/miRNA-186-5p in the mode of ceRNA. PIK3R3 was down-regulated in SKOV3 cells transfected with si-PIK3R3 accompanied with the down-regulation of Twist and Rac1, which suggests Twist and Rac1 may mediate changes in PIK3R3 expression and affecting EMT in tumor cells [34,35,36].
In this study, we used lncRNA, microRNA, and mRNA biochips as carriers to discover a number of differentially expressed genes in ovarian cancer tissues. We utilized GEO, KEGG, DIANA and other common tumor databases for bioinformatics analysis. The results suggest that HOXD-AS1, as a ceRNA, upregulates PIK3R3 expression by sponging miR-186-5p. The higher expression of HOXD-AS1 is associated with shorter PFS and OS of EOC patients. Collectively, lncRNA-HOXD-AS1 may function as a promising biomarker and therapeutic target for human EOC treatment.
While this paper is in the process of preparation, Wang et al. reported HOXD-AS1 could regulate cell migration and invasion through miR-608 and FZD4 in ovarian cancer . Therefore, these data suggest that HOXD-AS1 may function as a master regulator of metastasis through its interaction with a number of miRNAs and their respective target genes. HOXD-AS1 functions as a cellular sponge to modulate the availability of specific miRNAs at the 3’UTR of the miRNA target genes in ceRNA theory .
In summary, we profiled the transcriptome landscape to report significant differentially expressed mRNA, lncRNA and miRNAs in EOC tissues compared with normal ovary tissues. We delineated a novel HOXD-AS1/miR-186-5p/PIK3R3 molecular pathway that regulates cell migration, invasion, and EMT process in EOC. Suppressing HOXD-AS1 can be a possible strategy to improve the clinical outcome of EOC patients.
Rho GTPase activating protein 11A
Carbohydrate antigen 125
Competing endogenous RNA
The Database for Annotation, Visualization and Integrated Discovery
E2F transcription factor 8
Epithelial ovarian cancer
Fetal bovine serum
False discovery rate
International Federation of Gynecology and Obstetrics
Gene Expression Omnibus
HOXD antisense growth-associated long non-coding RNA
Kyoto Encyclopedia of Genes and Genomes
Long non-coding RNAs
The promoter of CDKN1A antisense DNA damage activated RNA
Phosphate buffer saline
A regulatory subunit of phosphatidylinositol 3-kinase (PI3K)
Pro-transition associated RNA
Plasmacytoma variant translocation 1
Rac family small GTPase 1
Ras homolog gene family member C
Reverse transcription quantitative PCR
Significance analysis of microarray
Sodium dodecyl sulfate- polyacrylamide gel electrophoresis
A member of the SOX (SRY-related HMG-box) gene family
Twist family bHLH transcription factor 1
X-inactive specific transcript
Zinc finger E-box binding homeobox 1
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65(2):87–108.
Fujiwara K, McAlpine JN, Lheureux S, Matsumura N, Oza AM. Paradigm shift in the management strategy for epithelial ovarian Cancer. Am Soc Clin Oncol Educ Book. 2016;35:e247–57.
Nasioudis D, Sisti G, Kanninen TT, Holcomb K, Di Tommaso M, Fambrini M, Witkin SS. Epidemiology and outcomes of squamous ovarian carcinoma; a population-based study. Gynecol Oncol. 2016;141(1):128–33.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66(1):7–30.
Kalluri R, Weinberg RA. The basics of epithelial-mesenchymal transition. J Clin Invest. 2009;119(6):1420–8.
Nieto MA, Huang RY, Jackson RA, Thiery JP. Emt: 2016. Cell. 2016;166(1):21–45.
Polyak K, Weinberg RA. Transitions between epithelial and mesenchymal states: acquisition of malignant and stem cell traits. Nat Rev Cancer. 2009;9(4):265–73.
Acloque H, Adams MS, Fishwick K, Bronner-Fraser M, Nieto MA. Epithelial-mesenchymal transitions: the importance of changing cell state in development and disease. J Clin Invest. 2009;119(6):1438–49.
Zheng X, Carstens JL, Kim J, Scheible M, Kaye J, Sugimoto H, Wu CC, LeBleu VS, Kalluri R. Epithelial-to-mesenchymal transition is dispensable for metastasis but induces chemoresistance in pancreatic cancer. Nature. 2015;527(7579):525–30.
Fischer KR, Durrans A, Lee S, Sheng J, Li F, Wong ST, Choi H, El Rayes T, Ryu S, Troeger J, et al. Epithelial-to-mesenchymal transition is not required for lung metastasis but contributes to chemoresistance. Nature. 2015;527(7579):472–6.
Liang SQ, Marti TM, Dorn P, Froment L, Hall SR, Berezowska S, Kocher G, Schmid RA, Peng RW. Blocking the epithelial-to-mesenchymal transition pathway abrogates resistance to anti-folate chemotherapy in lung cancer. Cell Death Dis. 2015;6:e1824.
Huang L, Wu RL, Xu AM. Epithelial-mesenchymal transition in gastric cancer. Am J Transl Res. 2015;7(11):2141–58.
Zhang PF, Li KS, Shen YH, Gao PT, Dong ZR, Cai JB, Zhang C, Huang XY, Tian MX, Hu ZQ, et al. Galectin-1 induces hepatocellular carcinoma EMT and sorafenib resistance by activating FAK/PI3K/AKT signaling. Cell Death Dis. 2016;7:e2201.
Liu S, Sun J, Cai B, Xi X, Yang L, Zhang Z, Feng Y, Sun Y. NANOG regulates epithelial-mesenchymal transition and chemoresistance through activation of the STAT3 pathway in epithelial ovarian cancer. Tumour Biol. 2016;37(7):9671–80.
Huarte M. The emerging role of lncRNAs in cancer. Nat Med. 2015;21(11):1253–61.
Wang Y, Lee CG. MicroRNA and cancer--focus on apoptosis. J Cell Mol Med. 2009;13(1):12–23.
Chen DL, Chen LZ, Lu YX, Zhang DS, Zeng ZL, Pan ZZ, Huang P, Wang FH, Li YH, Ju HQ, et al. Long noncoding RNA XIST expedites metastasis and modulates epithelial-mesenchymal transition in colorectal cancer. Cell Death Dis. 2017;8(8):e3011.
Chang Z, Cui J, Song Y. Long noncoding RNA PVT1 promotes EMT via mediating microRNA-186 targeting of Twist1 in prostate cancer. Gene. 2018;654:36–42.
Barabasi AL, Oltvai ZN. Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004;5(2):101–13.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Liu J, Matulonis UA. New strategies in ovarian cancer: translating the molecular complexity of ovarian cancer into treatment advances. Clin Cancer Res. 2014;20(20):5150–6.
Davidson B, Trope CG. Ovarian cancer: diagnostic, biological and prognostic aspects. Womens Health (Lond). 2014;10(5):519–33.
Ombrato L, Malanchi I. The EMT universe: space between cancer cell dissemination and metastasis initiation. Crit Rev Oncog. 2014;19(5):349–61.
Singh A, Settleman J. EMT, cancer stem cells and drug resistance: an emerging axis of evil in the war on cancer. Oncogene. 2010;29(34):4741–51.
Kun ZYX. Chemoresistance to CDDP induces EMT through activation of RAC1 in SKOV3 human ovarian Cancer cells. Prog Biochem Biophys. 2015;42(2):195–6.
Wu DD, Chen X, Sun KX, Wang LL, Chen S, Zhao Y. Role of the lncRNA ABHD11-AS1 in the tumorigenesis and progression of epithelial ovarian cancer through targeted regulation of RhoC. Mol Cancer. 2017;16(1):138.
Liang H, Yu T, Han Y, Jiang H, Wang C, You T, Zhao X, Shan H, Yang R, Yang L, et al. LncRNA PTAR promotes EMT and invasion-metastasis in serous ovarian cancer by competitively binding miR-101-3p to regulate ZEB1 expression. Mol Cancer. 2018;17(1):119.
Wang H, Fang L, Jiang J, Kuang Y, Wang B, Shang X, Han P, Li Y, Liu M, Zhang Z, et al. The cisplatin-induced lncRNA PANDAR dictates the chemoresistance of ovarian cancer via regulating SFRS2-mediated p53 phosphorylation. Cell Death Dis. 2018;9(11):1103.
Zhong Y, Gao D, He S, Shuai C, Peng S. Dysregulated expression of long noncoding RNAs in ovarian Cancer. Int J Gynecol Cancer. 2016;26(9):1564–70.
Lu S, Zhou J, Sun Y, Li N, Miao M, Jiao B, Chen H. The noncoding RNA HOXD-AS1 is a critical regulator of the metastasis and apoptosis phenotype in human hepatocellular carcinoma. Mol Cancer. 2017;16(1):125.
Wang H, Huo X, Yang XR, He J, Cheng L, Wang N, Deng X, Jin H, Wang N, Wang C, et al. STAT3-mediated upregulation of lncRNA HOXD-AS1 as a ceRNA facilitates liver cancer metastasis by regulating SOX4. Mol Cancer. 2017;16(1):136.
Chen Y, Zhao F, Cui D, Jiang R, Chen J, Huang Q, Shi J. HOXD-AS1/miR-130a sponge regulates glioma development by targeting E2F8. Int J Cancer. 2018;142(11):2313–22.
Wang H, Shen Q, Zhang X, Yang C, Cui S, Sun Y, Wang L, Fan X, Xu S. The long non-coding RNA XIST controls non-small cell lung Cancer proliferation and invasion by modulating miR-186-5p. Cell Physiol Biochem. 2017;41(6):2221–9.
Mei Q, Li X, Zhang K, Wu Z, Li X, Meng Y, Guo M, Luo G, Fu X, Han W. Genetic and methylation-induced loss of miR-181a2/181b2 within chr9q33.3 facilitates tumor growth of cervical Cancer through the PIK3R3/Akt/FoxO signaling pathway. Clin Cancer Res. 2017;23(2):575–86.
Yu T, Li J, Yan M, Liu L, Lin H, Zhao F, Sun L, Zhang Y, Cui Y, Zhang F, et al. MicroRNA-193a-3p and -5p suppress the metastasis of human non-small-cell lung cancer by downregulating the ERBB4/PIK3R3/mTOR/S6K2 signaling pathway. Oncogene. 2015;34(4):413–23.
Peng YP, Zhu Y, Yin LD, Wei JS, Liu XC, Zhu XL, Miao Y. PIK3R3 promotes metastasis of pancreatic Cancer via ZEB1 induced epithelial-mesenchymal transition. Cell Physiol Biochem. 2018;46(5):1930–8.
Wang Y, Zhang W, Wang Y, Wang S. HOXD-AS1 promotes cell proliferation, migration and invasion through miR-608/FZD4 axis in ovarian cancer. Am J Cancer Res. 2018;8(1):170–82.
We would like to thank our colleagues especially Dr. Junjun Li from the department of pathology for their assistance in this project; and thank the staff of Hunan Cancer Hospital / The Affiliated Cancer Hospital of Xiangya School of Medicine; Central South University for their assistance to the study.
This study was supported by the National Natural Science Foundation (No. 81201730, No. 81703005 and No.81602480), Key Research and Development Project of Hunan Science and Technology Department (No. 2018SK2126), Key Project of Changsha Science and Technology Bureau (No. kq1706046), Research Project of National Cancer Center Cancer (No.NCC2017A21), Research Project of Health and Family Planning Commission of Hunan Province (No. B2017098), the Natural Science Foundation of Hunan Province (No. 12JJ5073, No. 2017JJ3195, and No. 2018JJ3311), CBRG-NIG grant (NMRC/BNIG/2037/2015) and OF-YIRG (NMRC/OFYIRG/0076/2018) from National Medical Research Council of Singapore and SingHealth ONCO-ACP grant to YW (NCCRF-OACPCCS-YR2015-Aug-1).
Availability of data and materials
The lncRNA+mRNA arrays and miRNA arrays data have been deposited into GEO database under the accession number GSE119056. The reviewer’s link is provided below. (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE119056).
An additional set of miRNA expression microarray data, GSE53829, was retrieved from GEO database.
Ethics approval and consent to participate
The research was approved by the Ethics Committee of Hunan Cancer Hospital / The Affiliated Cancer Hospital of Xiangya School of Medicine; Central South University.
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.