- Open Access
Long noncoding RNA NORAD, a novel competing endogenous RNA, enhances the hypoxia-induced epithelial-mesenchymal transition to promote metastasis in pancreatic cancer
Molecular Cancervolume 16, Article number: 169 (2017)
Pancreatic cancer, one of the top two most fatal cancers, is characterized by a desmoplastic reaction that creates a dense microenvironment, promoting hypoxia and inducing the epithelial-to-mesenchymal transition (EMT) to facilitate invasion and metastasis. Recent evidence indicates that the long noncoding RNA NORAD may be a potential oncogenic gene and that this lncRNA is significantly upregulated during hypoxia. However, the overall biological role and clinical significance of NORAD remains largely unknown.
NORAD expression was measured in 33 paired cancerous and noncancerous tissue samples by real-time PCR. The effects of NORAD on pancreatic cancer cells were studied by overexpression and knockdown in vitro. Insights into the mechanism of competitive endogenous RNAs (ceRNAs) were gained from bioinformatics analyses and luciferase assays. In vivo, metastatic potential was identified using an orthotopic model of PDAC and quantified using bioluminescent signals. Alterations in RhoA expression and EMT levels were identified and verified by immunohistochemistry and Western blotting.
NORAD is highly expressed in pancreatic cancer tissues and upregulated in hypoxic conditions. NORAD upregulation is correlated with shorter overall survival in pancreatic cancer patients. Furthermore, NORAD overexpression promoted the migration and invasion of pancreatic carcinoma cells, while NORAD depletion inhibited EMT and metastasis in vitro and in vivo. In particular, NORAD may function as a ceRNA to regulate the expression of the small GTP binding protein RhoA through competition for hsa-miR-125a-3p, thereby promoting EMT.
Elevated expression of NORAD in pancreatic cancer tissues is linked to poor prognosis and may confer a malignant phenotype upon tumor cells. NORAD may function as a ceRNA to regulate the expression of the small GTP binding protein RhoA through competition for hsa-miR-125a-3p. This finding may contribute to a better understanding of the role played by lncRNAs in hypoxia-induced EMT and provide a potential novel diagnostic and therapeutic target for pancreatic cancer.
In 2012, an estimated 338 000 new cases of pancreatic cancer were diagnosed worldwide, while the disease was responsible for 331 000 deaths . In contrast to the steady increase in survival rates for most cancers, advances have been slow for pancreatic cancers, for which the 5-year relative survival remained at approximately 8% for 2005–2011, only 4–5% higher than 1975–1989 . Both morbidity and mortality continue to increase. Pancreatic cancer is projected to be one of the top two cancers in terms of fatalities in the next ten years [2, 3]. At present, surgical resection offers the only chance for a cure. Unfortunately, 80–85% of patients present with advanced metastasis and unresectable tumors when diagnosed . Metastasis is the leading cause of death among pancreatic cancer patients . However, the underlying molecular mechanisms that trigger metastasis remain unclear. It is believed that the hypoxic microenvironment and hypoxia-induced epithelial-to-mesenchymal transition (EMT) are critical drivers of pancreatic cancer metastasis and progression [6,7,8]. These processes occur via activation of various signaling pathways such as Hedgehog signaling [7, 9], PI3K/Akt signaling  and Notch signaling .
Recent studies have revealed that the aberrant expression of noncoding RNAs plays an important role in gene control and may play a role in cancer biology as oncogenes or tumor suppressor genes . Existing evidence shows that several noncoding RNAs, such as miR-1236 , miR-30c  and miRNA-34a , also participate in regulating hypoxia-induced EMT. However, it remains to be shown whether long noncoding RNAs (lncRNAs) participate in modulating hypoxia-induced EMT.
It has been reported that lncRNAs can act as competing endogenous RNAs (ceRNAs) or “RNA sponges”, interacting with microRNAs in a manner that can sequester these molecules and reduce their regulatory effect on target mRNAs [16, 17]. In other words, lncRNAs, mRNA transcripts and false gene transcripts can affect each other by competitively combining with a miRNA response element (MRE) to influence post-transcriptional regulation .
For instance, the long non-coding RNA Unigene56159 promotes metastasis by acting as a ceRNA for miR-140-5p in hepatocellular carcinoma cells . The lncRNA HOTAIR may act as a ceRNA, effectively becoming a sink for miR-331-3p, thereby modulating the derepression of HER2 and promoting the proliferation, migration and invasion of gastric carcinoma cells .
NORAD (annotated as LINC00657 in RefSeq) is a highly abundant, conserved mammalian noncoding RNA that can directly regulate both ploidy and chromosomal stability by sequestering PUMILIO proteins . Furthermore, NORAD is reported to be a potential oncogenic gene and is associated with overall survival in breast cancer . NORAD is more highly expressed in the cytoplasm and is significantly upregulated by hypoxia in endothelial cells . However, the mechanisms behind these events have not been revealed.
In the present study, we demonstrated that the lncRNA NORAD is upregulated during hypoxia and promotes EMT in tumor cells. NORAD facilitates invasion and metastasis by increasing the expression of the small GTP binding protein RhoA in vitro and in vivo. Mechanistically, NORAD may function as a ceRNA to regulate the expression of RhoA through competition for miR-125a-3p, thus playing an oncogenic role in the pathogenesis of pancreatic cancer. We also found that NORAD is upregulated in pancreatic ductal adenocarcinomas (PDACs). Elevated NORAD expression in tumor tissues is linked to poor prognosis and recurrence, demonstrating the clinical significance of this molecule during cancer progression.
Microarray analysis and bioinformatics methods
Expression profiles were obtained from the GEO database (www.ncbi.nim.nih.gov/geo/) as mentioned previously . Expression data from 55 pancreatic tumor tissues and their matched non-tumor tissues from two human microarray datasets GSE15471  and GSE16515  (39 pairs from GSE15471 and 16 pairs from GSE16515) were obtained and compared using the robust multiarray average (RMA) algorithm. Next, we used Venn analysis to analyze aberrant lncRNAs. Afterward, GSE32688  was analyzed for the expression of miR-125a-3p and RhoA. In addition, box plots (Additional file 1: Figure S1) were made to graphically view the value distributions among datasets. By analyzing the boxplots, these databases were confirmed to be median-centered and cross-comparable.
The potential microRNA binding sites within NORAD were predicted by computer-aided algorithms obtained from Starbase v2.0 (http://starbase.sysu.edu.cn) [28, 29]. MicroRNA-mRNA binding sites were obtained from miRTarBase (http://mirtarbase.mbc.nctu.edu.tw) .
PDAC specimens and the corresponding adjacent non-tumor tissues were obtained from patients who had undergone curative surgery at Shanghai Jiao Tong University School of Medicine Affiliated Ruijin Hospital between 2013 and 2015. None of them had received radiotherapy or chemotherapy before their surgery. Staging had been performed clinically, radiographically, and pathologically. The fresh primary PDAC tissues and adjacent non-tumor pancreas tissue samples were collected and snap-frozen in liquid nitrogen immediately after resection, then stored at −80 °C for further use.
Cell culture and hypoxia model
HEK-293Ts and the pancreatic cancer cell lines SW1990, Capan-1, PANC-1, AsPC-1, CFPAC-1, MIAPaCa-2, and BxPC-3 were purchased from ATCC (www.atcc.org). Premalignant HPDE cells (HPV16-immortalized normal pancreatic ductal epithelium) were generously provided by Shanghai Cell Bank of the Chinese Academy of Sciences. All cell lines were tested and authenticated by DNA typing at the Shanghai Jiao Tong University Analysis Core. All cells were cultured in DMEM, IMDM or RPMI 1640 (Gibco) with 10% Fetal Bovine Serum (FBS, Gibco) at 37 °C under 5% CO2 in a humidified chamber.
During hypoxic treatment conditions , cells were cultured in a modular incubator chamber (Billups-Rothenberger) that was flushed with 1% O2, 5% CO2, and 94% N2. Cells were harvested after 48 h. During normoxic control treatments, cells were cultured in an incubator with a humidified atmosphere containing 5% CO2 and 20% O2.
Quantitative real-time PCR
Total RNA was isolated from cells with TRIzol reagent (Invitrogen, Carlsbad, CA, USA). Next, cDNA was synthesized with a reverse transcription kit (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. qRT-PCR was performed on the 7500 Fast Real-time PCR system (Applied Biosystems) using a standard protocol from the SYBR Green PCR Kit (Toyobo). The results were normalized to the expression of GAPDH. Primers used for qRT-PCR assays are listed in Additional file 2: Table S1. For the detection of miRNA expression, reverse transcription was performed and microRNAs were detected with stem-loop primers purchased from RiboBio (Guangzhou, China). U6 snoRNA was used as the endogenous control. Relative fold changes were calculated using the 2-△△Ct method. All PCR assays were repeated three times.
NORAD cDNA was cloned into the mammalian expression vector pcDNA3.1 (Invitrogen). NORAD cDNA fragments containing either the predicted potential microRNA binding sites (wildtype, wt) or scrambled microRNA binding site sequences (mutation, mut) were amplified by PCR. Moreover, plasmids containing the firefly luciferase reporter were constructed with RhoA-3′-UTR-wildtype ( RhoA-3′-UTR-wt-luc) and RhoA-3′-UTR-mutation (mutated miR-125a-3p binding region; RhoA-3′-UTR-mut-luc). Mimics and inhibitors of hsa-miR-125a-3p and rno-miR-239b-5p were purchased from RiboBio (Guangzhou, China). rno-miR-239b-5p was used as a negative control.
Virus packaging and infection
shRNA sequences targeting NORAD were designed, inserted in lentiviral plasmids and then transfected into HEK-293Ts with virus packaging plasmids to produce the lentivirus. The viral supernatants were collected at 48 and 72 h after transfection and added into PANC-1 and SW1990 cells to construct the stably transfected cell lines. Puromycin was added to the media 48 h after infection and maintained for at least 1 week to select stably transfected cell lines (PANC-1/SW1990-sh-NORAD or PANC-1/SW1990-sh-NC).
Wound healing assay
PDAC cells (1 × 106 cells/well) were treated with different reagents, seeded in six-well plates and cultured until they reached confluence. Wounds were made in the cell monolayer by making a scratch with a 20 μL pipette tip. Plates were washed once with fresh medium after 24 h in culture to remove non-adherent cells. Following this wash, plates were photographed.
Transwell migration and invasion assay
Approximately 1 × 105 PDAC cells treated with different reagents were suspended in serum-free medium and seeded in either Matrigel-coated chambers with 8-μm pores (BD Biosciences, Franklin Lakes, NY, USA) for migration assays or in chambers (Corning Costar, NY, USA) that were not coated with Matrigel for invasion assays. For both assays, medium containing 10% FBS was added to the lower chamber as a chemoattractant. After 12 h, the cells that had invaded through the membrane to the lower surface were fixed with methanol, stained with 0.5% crystal violet solution and counted in 5 random fields of view (200×). RhoA pathway specific inhibitor CCG-1423  was purchased from selleck and used as recommended.
In vivo orthotopic model
Athymic BALB/C mice (male, 4–6 weeks old) were purchased from the Chinese Academy of Sciences (Shanghai, China) and were maintained in a specific pathogen-free facility. Orthotopic implantation of pancreatic cancer cell lines was performed as described previously . For experimental metastasis assays, the SW1990-sh-NORAD and SW1990-sh-NC cell lines were labeled with firefly luciferase. Each mouse was then injected with 100 μL of cell suspension (1 × 106 cells) into the pancreas. Animal health conditions and metastatic progression were monitored. Metastasis was quantified using a noninvasive bioluminescence In Vivo Imaging System (IVIS, Xenogen) as described previously . All mice were sacrificed and both the primary tumor and metastatic foci were observed by general observation, H & E staining and immunohistochemical staining.
Dual-luciferase reporter assay
Recombinant plasmids or an empty plasmid encoding the firefly luciferase reporter were transiently transfected into SW1990-sh-NORAD and SW1990-sh-NC cells using Lipofectamine 2000 (Invitrogen, USA) The renilla luciferase reporter pRL-CMV plasmid (Promega) was also transfected into these cells as an internal control. Each group above was transfected with miR-125a-3p mimic or a negative control. Luciferase activities in the indicated cells were measured using the Dual-Luciferase Reporter Assay System (Promega, USA) 48 h after transfection according to the manufacturer’s instructions.
Western blot assay
Cells were lysed with RIPA cell lysis buffer in the presence of a protease inhibitor cocktail (Sigma, USA) in order to harvest total cell protein. Standard Western blotting was carried out as described previously . In brief, the BCA Protein Assay Kit (Pierce, USA) was used to quantify the protein concentration of each cell lysate. Equal amounts of protein samples were loaded onto a 10% SDS-PAGE gel and then transferred onto PVDF membranes. After blocking with skim milk, the membranes were incubated with primary antibodies diluted in TBST buffer overnight at 4 °C and then with the HRP-conjugated secondary antibody for 2–3 h at room temperature. Finally, the images of protein bands were captured by a Tanon detection system using ECL reagent (Thermo).
Data were presented as the mean ± one standard deviation (SD) from at least three separate experiments. The Student t test, X2 test, log-rank test, nonparametric Mann–Whitney U test, or Fisher exact test was used for comparisons between groups. The Kaplan–Meier method was used to estimate overall survival (os) and recurrence rate. Differences were deemed statistically significant at P < 0.05. The results are shown as the mean ± SEM.
This study was approved by the Clinical Research Ethics Committee of Shanghai Jiao-Tong University. All patients were fully informed of the experimental procedures before tissue acquisition. All animal experimental procedures were performed in compliance with the institutional ethical requirements and were approved by the Shanghai Jiao-Tong University School of Medicine Committee for the Use and Care of Animals.
LncRNA expression profiles in PDAC
To identify eligible lncRNAs in PDAC, we analyzed microarray expression profiles. Two human PDAC microarray datasets (GSE15471 and GSE16515) were selected and analyzed for consistently aberrant lncRNAs between PDAC and normal pancreas tissues using GEO2R. Fold changes in robust multiarray average (RMA) signal intensities (which represent expression differences between tumor and normal tissue) were calculated and cutoff values were set at |RMA fold change| > 1.5 and a p value < 0.05. Subsequently, we identified 272 consistent non-coding RNAs that were significantly upregulated or downregulated by Venn analysis (Fig. 1a).
NORAD was notably increased during hypoxia
As is already known, intratumoral hypoxia plays a crucial role in pancreatic cancer metastasis. After the above analyses, aberrant lncRNAs (Fig. 1b) with potential roles in the hypoxia pathway were chosen for further study.
To evaluate hypoxia-induced alterations in lncRNA expression, the pancreatic cancer cell line SW1990 was cultured under hypoxic or normoxic conditions for 48 h. We found that NORAD was significantly increased after hypoxic stimulation (Fig. 1c). Additionally, we found that NEAT1, PVT1, NBR2, HOTAIR, MALAT1 and NORAD increased during hypoxia, while MEG3 decreased during hypoxia.
NORAD promotes tumor metastasis in vitro and in vivo
As NORAD was upregulated under hypoxia-stimulating conditions, we hypothesized that superfluous NORAD may participate in the development of hypoxia-induced malignant phenotypes.
To investigate further, we selected SW1990 and PANC-1, two pancreatic cell lines with highest NORAD expression among 8 cell lines (shown in Additional file 3: Figure S2), and knocked down NORAD using shRNAs and then assessed tumor cell migration and invasion abilities using wound-healing assays as well as transwell migration and invasion assays. As shown in Fig. 2a-d, knockdown of NORAD in PANC-1 or SW1990 cells markedly suppressed cell migration and invasion compared with their respective control cells.
To further quantify metastatic and dissemination potential in vivo, we established an orthotopic metastatic model, as mentioned above. SW1990-sh-NORAD and SW1990-sh-NC cells were labeled with firefly luciferase and injected into the pancreas of nude mice. After 3 weeks, the scope of metastasis was quantified by non-invasively detecting the bioluminescent signal on a bioluminescence In Vivo Imaging System (Fig. 2e-f). Signals from the group injected with SW1990-sh-NC cells were significantly stronger than those from the group injected with SW1990-sh-NORAD cells (p = 0.043). Moreover, more distant micrometastatic foci were observed in the negative control group than the sh-NORAD group by H & E staining. It is very obvious that metastasis was repressed when knocking down NORAD (Fig. 2g).
NORAD inhibits miR-125a-3p via direct binding
Accumulating evidence has suggested that lncRNAs might function as ceRNAs by binding to miRNAs and functionally liberating other RNA transcripts . To examine whether NORAD operates by a similar mechanism in PDAC cells, we used Starbase v2.0 to predict the potential miRNA binding sites in NORAD and found a sequence complementary to miRNA-125a-3p (Fig. 3a). Furthermore, we found that the expression of NORAD was negatively associated with the expression of miR-125a-3p in 33 PDAC tissues (Fig. 3b). We then tested the effects of shRNA-mediated downregulation of NORAD on the expression of miR-125a-3p. Notably, miR-125a-3p was upregulated in the sh-NORAD-treated SW1990 and PANC-1 cells compared with negative control-treated cells (Fig. 3c). Moreover, we selected BxPC-3 and Canpan-1, two pancreatic cell lines with lowest NORAD expression among 8 cell lines (shown in Additional file 3: Figue S2), and overexpress NORAD with plasmids containing wildtype NORAD (NORAD-wt) or NORAD with mutation in its MRE (NORAD-mut) respectively. After we overexpressed NORAD in BxPC-3 and Canpan-1, we found a significant decrease in miR-125a-3p. This phenomenon disappeared when the MRE of NORAD was mutated (Fig. 3d), suggesting that the lncRNA NORAD can repress the expression of miR-125a-3p via direct binding at the MRE.
NORAD increases the expression of RhoA acting as a ceRNA of miR-125a-3p
Using miRTarBase, we found that RhoA was a potential target of miR-125a-3p. Additionally, we found that RhoA expression is negatively correlated with miR-125a-3p in GSE32688  (data provided in Additional file 4: Table S2) and that miR-125a-3p can complement the seed region of RhoA. Interestingly, the MRE of NORAD was found to be similar to the seed region of RhoA and the expression levels of NORAD and RhoA were found to be positively correlated (Fig. 4a-c). RhoA is downregulated after NORAD knockdown in vitro and vivo (Fig. 4d). And we also found that RhoA is increased by NORAD-wt, while NORAD-mut failed to increase RhoA (Fig. 4e). This suggests that NORAD may enhance RhoA by repressing miR-125a-3p. To confirm the binding in predicted sites of RhoA-3′-UTR, we mutated binding sites in RhoA and then performed reporter assays. We found that miR-125a-3p could inhibit the activity of RhoA-3′-UTR-wt, but had no influence on the activity of RhoA-3′-UTR-mut (Fig. 4f). In addition, we found that downregulation of NORAD could significantly inhibit the activity of RhoA-3′-UTR. This effect could be partially retrieved using a miR-125a-3p inhibitor (rather than miR-239b-5p), which suggests that NORAD competes with the RhoA-3′-UTR for miR-125a-3p (Fig. 4g).
NORAD promotes EMT via regulating RhoA expression
As we know, hypoxic stimuli can trigger the expression of RhoA and its downstream effector ROCK1, leading to cell metastasis, and RhoA plays a crucial role in mediating the epithelial-mesenchymal transition (EMT) . To further investigate the molecular mechanisms underlying these events, we explored whether NORAD regulates EMT in cancer cells. As shown in Fig. 5a-b, knockdown of NORAD significantly decreased the expression levels of the mesenchymal cell markers N-cadherin, Vimentin, and ZEB1 but increased the expression levels of the epithelial cell marker E-cadherin. And overexpression of NORAD resulted in RhoA increasing and E-cadherin decreasing. This suggests that the more NORAD expressed the more RhoA protein will be produced, and then triggering more aggressive EMT and cell migration and invasion. Furthermore, as Fig. 5b-d shown, when we treated cell lines with RhoA pathway specific inhibitor CCG-1423 , RhoA downstream protein ROCK1 decreased and aggressive EMT and invasive behaviors caused by NORAD overexpression were partly redressed. In another words, blocking RhoA activity with the RhoA inhibitor can impede the flow of EMT induced by NORAD, indicating that NORAD promotes EMT via regulating RhoA.
Excess NORAD is associated with poor prognosis in PDAC
The expression of NORAD was significantly higher in 75 pancreatic cancer samples than in 55 normal pancreas tissues (GSE15471 and GSE16515, p < 0.001, Fig. 6a). We then assayed NORAD expression in a panel of 33 PDAC with complete clinicopathological information (Additional file 5: Table S3), along with 33 paired normal pancreas tissues. NORAD was significantly upregulated in PDAC specimens compared with normal pancreas tissues (Fig. 6b, p<0.01). The expression of NORAD was also upregulated in most tumor tissues compared with the non-tumor tissues from the same donor (Fig. 6c), thus linking the upregulation of NORAD to pathological changes in pancreatic tissue.
Next, we evaluated the clinical significance of NORAD in PDAC. Thirty-three PDAC patients from our hospital were divided into two groups (high: n = 16, low: n = 17) according to the level of NORAD expression.
The Kaplan–Meier estimator was used to estimate recurrence-free survival (RFS) and overall survival (OS) rates and p values were calculated using the log-rank test. As shown in Fig. 6d-e, PDAC patients with higher NORAD expression have shorter overall survival (p = 0.0335) and recurrence-free survival (p = 0.0419) rates. The median overall and recurrence-free survival rates were sharply decreased (RFS: 9.35 versus 19.8 months; OS: 14.35 versus 22.8 months) in the high-NORAD group, which indicates that NORAD upregulation is associated with poor prognosis.
Pancreatic cancer is a devastating disease in which the majority of patients succumb within one year of diagnosis. The disease is characterized by a desmoplastic reaction that creates a dense microenvironment, promoting hypoxia and limiting cancer drug delivery due to decreased blood perfusion . Stromal fibrosis and tissue necrosis are major histological sequelae of hypoxia . Standard chemotherapeutic agents fail because they are unsuccessful at targeting cells within the hypoxic microenvironment. Out of all therapeutic obstacles, hypoxia might be the most important to target for elimination .
Responses to hypoxia are thought to operate via hypoxia-inducible factors (HIFs). Both hypoxia and overexpression of hypoxia-inducible factor (HIF)-1α and/or HIF-2α have been shown to promote EMT and metastatic phenotypes . Mechanisms possibly connecting hypoxia with EMT include intracellular reactive oxygen species (ROS)-dependent HIF accumulation, Snail translocation and HIF1α-dependent accumulation of the transcription factors Snail, Twist, ZEB1 and ZEB2, key regulators of EMT . Numerous studies have suggested that EMT both contributes to the early-stage dissemination of cancer cells in PDAC [5, 41] and induces chemoresistance, which may be more important to the unsatisfactory prognosis in pancreatic cancer .
Although various potential candidates emerge for disrupting EMT or treating advanced-stage pancreatic cancer [43,44,45], this cancer remains resistant to conventional therapies and is susceptible to early distant metastasis . We believe that the mechanism of hypoxia-induced EMT is much more complicated than our current understanding and that the prognosis of pancreatic cancer will not be improved until we increase our understanding.
Recent evidence indicates that lncRNAs play essential roles in numerous biological processes, such as epigenetic regulation , cell cycle control , RNA decay and transcription [49, 50]. MiRNAs can function as master gene regulators, impacting a variety of cellular pathways involved in normal cellular functions as well as cancer development and progression . LncRNAs that act as competing endogenous RNAs (ceRNAs) or natural microRNA sponges are important post-transcriptional regulators of gene expression that communicate with and co-regulate each other by competing for binding to shared microRNAs. Understanding this novel RNA crosstalk will lead to significant insights into gene regulatory networks and human development and disease .
NORAD is highly conserved and abundant, with expression levels of approximately 500–1000 copies per cell . Notably, NORAD loss-of-function results in chromosomal instability and triggers dramatic aneuploidy in previously karyotypically stable cell lines. NORAD can directly regulate both ploidy and chromosomal stability by sequestering PUMILIO proteins, which repress key mitotic, DNA repair, and DNA replication factors . NORAD is a cytoplasmic multivalent PUMILIO binding platform that acts as a negative regulator of PUMILIO activity . In breast cancer, NORAD knockout significantly suppressed tumor cell growth and proliferation, suggesting that it plays an oncogenic role . However, this study was unable to shed light on its mechanism.
Previous studies have revealed that miR-125a-3p plays a tumor suppressive role in many cancers, including non-small-cell lung cancer , gastric cancer , hepatocellular carcinoma  and myeloma . Moreover, miR-125a-3p is widely involved in cell proliferation , apoptosis , differentiation, as well as migration and invasion. However, little is known about its function in PDAC or its link to NORAD. Only one study reports that assessing the levels of miR-125a-3p in combination with 7 other miRNAs can be clinically valuable for identifying patients with pancreatic and biliary-tract cancers who could benefit from surgical intervention .
To the best of our knowledge, we are the first to demonstrate that NORAD acts as a ceRNA of miR-125a-3p and that it is involved in regulating pancreatic cancer metastasis and invasion. In addition, miR-125a-3p has also been reported to participate in several other pathways such as the p53 pathway , the p38/ MAPK pathway  and the PI3K/AKT/mTOR pathway , all of which may be associated with the malignant phenotype mediated by NORAD/miR-125a-3p. In the present study, we focused on RhoA, which is not only an important downstream factor in hypoxia but also a crucial regulator of EMT. According to bioinformatics analyses, RhoA and NORAD share the same response elements for miR-125a-3p. We further confirmed that the expression level of miR-125a-3p can be downregulated by NORAD and that miR-125a-3p can bind to NORAD directly at its MRE. We found that NORAD-wt increased RhoA, NORAD knockdown downregulated RhoA, and NORAD-mut failed to increase RhoA levels. Knock-down of NORAD could partially rescue the hypoxia-induced suppression of RhoA and EMT. Thus, NORAD may act as a ceRNA to regulate the expression of RhoA by sequestering its inhibitor miR-125a-3p, thus bridging the response to hypoxia and EMT (Fig. 7).
In conclusion, we determined that NORAD is an independent prognostic factor for PDAC patients and that it promotes the metastasis and invasion of PDAC both in vitro and in vivo. We also shed light on the role of NORAD in hypoxia-induced EMT. We found that hypoxia can increase NORAD expression and stimulate EMT in pancreatic cancer cells and that NORAD promotes EMT and metastasis by regulating RhoA in a miR-125a-3p-dependent manner. Knockdown of NORAD can alleviate malignant phenomena caused by hypoxia. Thus, we believe that NORAD may be a potential novel therapeutic target for the treatment of PDAC.
Neoptolemos JP, Palmer DH, Ghaneh P, Psarelli EE, Valle JW, Halloran CM, Faluyi O, O'Reilly DA, Cunningham D, Wadsley J, et al. Comparison of adjuvant gemcitabine and capecitabine with gemcitabine monotherapy in patients with resected pancreatic cancer (ESPAC-4): a multicentre, open-label, randomised, phase 3 trial. Lancet. London: England; 2017.
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2016. CA Cancer J Clin. 2016;66(1):7–30.
Rahib L, Smith BD, Aizenberg R, Rosenzweig AB, Fleshman JM, Matrisian LM. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States. Cancer Res. 2014;74(11):2913–21.
Vincent A, Herman J, Schulick R, Hruban RH, Goggins M. Pancreatic cancer. Lancet (London, England). 2011;378(9791):607–20.
Rhim AD, Mirek ET, Aiello NM, Maitra A, Bailey JM, McAllister F, Reichert M, Beatty GL, Rustgi AK, Vonderheide RH, et al. EMT and dissemination precede pancreatic tumor formation. Cell. 2012;148(1–2):349–61.
Lei J, Ma J, Ma Q, Li X, Liu H, Xu Q, Duan W, Sun Q, Xu J, Wu Z, et al. Hedgehog signaling regulates hypoxia induced epithelial to mesenchymal transition and invasion in pancreatic cancer cells via a ligand-independent manner. Mol Cancer. 2013;12:66.
Spivak-Kroizman TR, Hostetter G, Posner R, Aziz M, Hu C, Demeure MJ, Von Hoff D, Hingorani SR, Palculict TB, Izzo J, et al. Hypoxia triggers hedgehog-mediated tumor-stromal interactions in pancreatic cancer. Cancer Res. 2013;73(11):3235–47.
Liu L, Salnikov AV, Bauer N, Aleksandrowicz E, Labsch S, Nwaeburu C, Mattern J, Gladkich J, Schemmer P, Werner J, et al. Triptolide reverses hypoxia-induced epithelial-mesenchymal transition and stem-like features in pancreatic cancer by NF-kappaB downregulation. Int J Cancer. 2014;134(10):2489–503.
Cao L, Xiao X, Lei J, Duan W, Ma Q, Li W. Curcumin inhibits hypoxia-induced epithelialmesenchymal transition in pancreatic cancer cells via suppression of the hedgehog signaling pathway. Oncol Rep. 2016;35(6):3728–34.
Du R, Xia L, Ning X, Liu L, Sun W, Huang C, Wang H, Sun S. Hypoxia-induced Bmi1 promotes renal tubular epithelial cell-mesenchymal transition and renal fibrosis via PI3K/Akt signal. Mol Biol Cell. 2014;25(17):2650–9.
Sahlgren C, Gustafsson MV, Jin S, Poellinger L, Lendahl U. Notch signaling mediates hypoxia-induced tumor cell migration and invasion. Proc Natl Acad Sci U S A. 2008;105(17):6392–7.
Schmitt AM, Chang HY. Long noncoding RNAs in cancer pathways. Cancer Cell. 2016;29(4):452–63.
Chen SY, Teng SC, Cheng TH, KJ W. miR-1236 regulates hypoxia-induced epithelial-mesenchymal transition and cell migration/invasion through repressing SENP1 and HDAC3. Cancer Lett. 2016;378(1):59–67.
Huang J, Yao X, Zhang J, Dong B, Chen Q, Xue W, Liu D, Huang Y. Hypoxia-induced downregulation of miR-30c promotes epithelial-mesenchymal transition in human renal cell carcinoma. Cancer Sci. 2013;104(12):1609–17.
Li H, Rokavec M, Jiang L, Horst D, Hermeking H. Antagonistic effects of p53 and HIF1A on microRNA-34a regulation of PPP1R11 and STAT3 and hypoxia-induced epithelial to Mesenchymal transition in colorectal cancer cells. Gastroenterology. 2017;153(2):505–20.
Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, Tramontano A, Bozzoni I. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011;147(2):358–69.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52.
MH L, Tang B, Zeng S, CJ H, Xie R, YY W, Wang SM, He FT, Yang SM. Long noncoding RNA BC032469, a novel competing endogenous RNA, upregulates hTERT expression by sponging miR-1207-5p and promotes proliferation in gastric cancer. Oncogene. 2016;35(27):3524–34.
Lv J, Fan HX, Zhao XP, Lv P, Fan JY, Zhang Y, Liu M, Tang H. Long non-coding RNA Unigene56159 promotes epithelial-mesenchymal transition by acting as a ceRNA of miR-140-5p in hepatocellular carcinoma cells. Cancer Lett. 2016;382(2):166–75.
Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, Kong R, Xia R, KH L, Li JH, et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 2014;13:92.
Lee S, Kopp F, Chang TC, Sataluri A, Chen B, Sivakumar S, Yu H, Xie Y, Mendell JT. Noncoding RNA NORAD regulates genomic stability by sequestering PUMILIO proteins. Cell. 2016;164(1–2):69–80.
Liu H, Li J, Koirala P, Ding X, Chen B, Wang Y, Wang Z, Wang C, Zhang X, Mo YY. Long non-coding RNAs as prognostic markers in human breast cancer. Oncotarget. 2016;7(15):20584–96.
Michalik KM, You X, Manavski Y, Doddaballapur A, Zornig M, Braun T, John D, Ponomareva Y, Chen W, Uchida S, et al. Long noncoding RNA MALAT1 regulates endothelial cell function and vessel growth. Circ Res. 2014;114(9):1389–97.
Li H, Wang X, Fang Y, Huo Z, Lu X, Zhan X, Deng X, Peng C, Shen B. Integrated expression profiles analysis reveals novel predictive biomarker in pancreatic ductal adenocarcinoma. Oncotarget. 2017;
Badea L, Herlea V, Dima SO, Dumitrascu T, Popescu I. Combined gene expression analysis of whole-tissue and microdissected pancreatic ductal adenocarcinoma identifies genes specifically overexpressed in tumor epithelia. Hepato-Gastroenterology. 2008;55(88):2016–27.
Pei H, Li L, Fridley BL, Jenkins GD, Kalari KR, Lingle W, Petersen G, Lou Z, Wang L. FKBP51 affects cancer cell response to chemotherapy by negatively regulating Akt. Cancer Cell. 2009;16(3):259–66.
Donahue TR, Tran LM, Hill R, Li Y, Kovochich A, Calvopina JH, Patel SG, Wu N, Hindoyan A, Farrell JJ, et al. Integrative survival-based molecular profiling of human pancreatic cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2012;18(5):1352–63.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.
Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH. StarBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2011;39(Database issue):D202–9.
Chou CH, Chang NW, Shrestha S, Hsu SD, Lin YL, Lee WH, Yang CD, Hong HC, Wei TY, Tu SJ, et al. miRTarBase 2016: updates to the experimentally validated miRNA-target interactions database. Nucleic Acids Res. 2016;44(D1):D239–47.
Aguilera KY, Rivera LB, Hur H, Carbon JG, Toombs JE, Goldstein CD, Dellinger MT, Castrillon DH, Brekken RA. Collagen signaling enhances tumor progression after anti-VEGF therapy in a murine model of pancreatic ductal adenocarcinoma. Cancer Res. 2014;74(4):1032–44.
Evelyn CR, Wade SM, Wang Q, Wu M, Iniguez-Lluhi JA, Merajver SD, Neubig RR. CCG-1423: a small-molecule inhibitor of RhoA transcriptional signaling. Mol Cancer Ther. 2007;6(8):2249–60.
Li X, Wang Z, Ma Q, Xu Q, Liu H, Duan W, Lei J, Ma J, Wang X, Lv S, et al. Sonic hedgehog paracrine signaling activates stromal cells to promote perineural invasion in pancreatic cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2014;20(16):4326–38.
Li T, Xie J, Shen C, Cheng D, Shi Y, Wu Z, Deng X, Chen H, Shen B, Peng C, et al. Amplification of long noncoding RNA ZFAS1 promotes metastasis in hepatocellular carcinoma. Cancer Res. 2015;75(15):3181–91.
Zhu Y, Shi M, Chen H, Gu J, Zhang J, Shen B, Deng X, Xie J, Zhan X, Peng C. NPM1 activates metabolic changes by inhibiting FBP1 while promoting the tumorigenicity of pancreatic cancer cells. Oncotarget. 2015;6(25):21443–51.
Gulhati P, Bowen KA, Liu J, Stevens PD, Rychahou PG, Chen M, Lee EY, Weiss HL, O'Connor KL, Gao T, et al. mTORC1 and mTORC2 regulate EMT, motility, and metastasis of colorectal cancer via RhoA and Rac1 signaling pathways. Cancer Res. 2011;71(9):3246–56.
Kong B, Cheng T, Wu W, Regel I, Raulefs S, Friess H, Erkan M, Esposito I, Kleeff J, Michalski CW. Hypoxia-induced endoplasmic reticulum stress characterizes a necrotic phenotype of pancreatic cancer. Oncotarget. 2015;6(31):32154–60.
Feig C, Gopinathan A, Neesse A, Chan DS, Cook N, Tuveson DA. The pancreas cancer microenvironment. Clinical cancer research : an official journal of the American Association for Cancer Research. 2012;18(16):4266–76.
Lu X, Kang Y. Hypoxia and hypoxia-inducible factors: master regulators of metastasis. Clinical cancer research : an official journal of the American Association for Cancer Research. 2010;16(24):5928–35.
Xi Y, Wei Y, Sennino B, Ulsamer A, Kwan I, Brumwell AN, Tan K, Aghi MK, McDonald DM, Jablons DM, et al. Identification of pY654-beta-catenin as a critical co-factor in hypoxia-inducible factor-1alpha signaling and tumor responses to hypoxia. Oncogene. 2013;32(42):5048–57.
Hotz B, Arndt M, Dullat S, Bhargava S, Buhr HJ, Hotz HG. Epithelial to mesenchymal transition: expression of the regulators snail, slug, and twist in pancreatic cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2007;13(16):4769–76.
Zheng X, Carstens JL, Kim J, Scheible M, Kaye J, Sugimoto H, CC W, LeBleu VS, Kalluri R. Epithelial-to-mesenchymal transition is dispensable for metastasis but induces chemoresistance in pancreatic cancer. Nature. 2015;527(7579):525–30.
Rumman M, Jung KH, Fang Z, Yan HH, Son MK, Kim SJ, Kim J, Park JH, Lim JH, Hong S, et al. HS-173, a novel PI3K inhibitor suppresses EMT and metastasis in pancreatic cancer. Oncotarget. 2016;7(47):78029–47.
Kashiwaya K, Hosokawa M, Eguchi H, Ohigashi H, Ishikawa O, Shinomura Y, Nakamura Y, Nakagawa H. Identification of C2orf18, termed ANT2BP (ANT2-binding protein), as one of the key molecules involved in pancreatic carcinogenesis. Cancer Sci. 2009;100(3):457–64.
Han Y, GQ R, Mou X, Wang HJ, Ma Y, He XL, Yan Z, Huang D. AUTS2 is a potential therapeutic target for pancreatic cancer patients with liver metastases. Med Hypotheses. 2015;85(2):203–6.
Zhan HX, JW X, Wu D, Zhang TP, SY H. Pancreatic cancer stem cells: new insight into a stubborn disease. Cancer Lett. 2015;357(2):429–37.
Mercer TR, Mattick JS. Structure and function of long noncoding RNAs in epigenetic regulation. Nat Struct Mol Biol. 2013;20(3):300–7.
Cheng Y, Jutooru I, Chadalapaka G, Corton JC, Safe S. The long non-coding RNA HOTTIP enhances pancreatic cancer cell proliferation, survival and migration. Oncotarget. 2015;6(13):10840–52.
Hirata H, Hinoda Y, Shahryari V, Deng G, Nakajima K, Tabatabai ZL, Ishii N, Dahiya R. Long noncoding RNA MALAT1 promotes aggressive renal cell carcinoma through Ezh2 and interacts with miR-205. Cancer Res. 2015;75(7):1322–31.
Zheng S, Chen H, Wang Y, Gao W, Fu Z, Zhou Q, Jiang Y, Lin Q, Tan L, Ye H, et al. Long non-coding RNA LOC389641 promotes progression of pancreatic ductal adenocarcinoma and increases cell invasion by regulating E-cadherin in a TNFRSF10A-related manner. Cancer Lett. 2016;371(2):354–65.
Chen L, Kang C. miRNA interventions serve as 'magic bullets' in the reversal of glioblastoma hallmarks. Oncotarget. 2015;6(36):38628–42.
Tichon A, Gil N, Lubelsky Y, Havkin Solomon T, Lemze D, Itzkovitz S, Stern-Ginossar N, Ulitsky I. A conserved abundant cytoplasmic long noncoding RNA modulates repression by Pumilio proteins in human. Cell. 2016;7:12209.
Zhang H, Zhu X, Li N, Li D, Sha Z, Zheng X, Wang H. miR-125a-3p targets MTA1 to suppress NSCLC cell proliferation, migration, and invasion. Acta Biochim Biophys Sin. 2015;47(7):496–503.
Hashiguchi Y, Nishida N, Mimori K, Sudo T, Tanaka F, Shibata K, Ishii H, Mochizuki H, Hase K, Doki Y, et al. Down-regulation of miR-125a-3p in human gastric cancer and its clinicopathological significance. Int J Oncol. 2012;40(5):1477–82.
Tang H, Li RP, Liang P, Zhou YL, Wang GW. miR-125a inhibits the migration and invasion of liver cancer cells via suppression of the PI3K/AKT/mTOR signaling pathway. Oncol Lett. 2015;10(2):681–6.
Alzrigat M, Parraga AA, Agarwal P, Zureigat H, Osterborg A, Nahi H, Ma A, Jin J, Nilsson K, Oberg F, et al. EZH2 inhibition in multiple myeloma downregulates myeloma associated oncogenes and upregulates microRNAs with potential tumor suppressor functions. Oncotarget. 2017;8(6):10213–24.
Ninio-Many L, Grossman H, Levi M, Zilber S, Tsarfaty I, Shomron N, Tuvar A, Chuderland D, Stemmer SM, Ben-Aharon I, et al. MicroRNA miR-125a-3p modulates molecular pathway of motility and migration in prostate cancer cells. Oncoscience. 2014;1(4):250–61.
Kojima M, Sudo H, Kawauchi J, Takizawa S, Kondou S, Nobumasa H, Ochiai A. MicroRNA markers for the diagnosis of pancreatic and biliary-tract cancers. PLoS One. 2015;10(2):e0118220.
Jiang L, Chang J, Zhang Q, Sun L, Qiu X. MicroRNA hsa-miR-125a-3p activates p53 and induces apoptosis in lung cancer cells. Cancer Investig. 2013;31(8):538–44.
Dong Y, Li P, Ni Y, Zhao J, Liu Z. Decreased microRNA-125a-3p contributes to upregulation of p38 MAPK in rat trigeminal ganglions with orofacial inflammatory pain. PLoS One. 2014;9(11):e111594.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Box plots of GSE15471 and GSE16515. (TIFF 154 kb)
Primers used for qRT-PCR assays. (XLSX 9 kb)
NORAD expression among 8 cell lines. (TIFF 728 kb)
The expression of RhoA and miR-125a-3p in GSE32688. (XLSX 11 kb)
Clinicopathological information of 33 PDAC. (XLSX 9 kb)