- Open Access
TEAD4 modulated LncRNA MNX1-AS1 contributes to gastric cancer progression partly through suppressing BTG2 and activating BCL2
Molecular Cancer volume 19, Article number: 6 (2020)
Gastric cancer (GC) is the third leading cause of cancer-related mortality globally. Long noncoding RNAs (lncRNAs) are dysregulated in obvious malignancies including GC and exploring the regulatory mechanisms underlying their expression is an attractive research area. However, these molecular mechanisms require further clarification, especially upstream mechanisms.
LncRNA MNX1-AS1 expression in GC tissue samples was investigated via microarray analysis and further determined in a cohort of GC tissues via quantitative reverse transcription polymerase chain reaction (qRT-PCR) assays. Cell proliferation and flow cytometry assays were performed to confirm the roles of MNX1-AS1 in GC proliferation, cell cycle regulation, and apoptosis. The influence of MNX1-AS1 on GC cell migration and invasion was explored with Transwell assays. A xenograft tumour model was established to verify the effects of MNX1-AS1 on in vivo tumourigenesis. The TEAD4-involved upstream regulatory mechanism of MNX1-AS1 was explored through ChIP and luciferase reporter assays. The mechanistic model of MNX1-AS1 in regulating gene expression was further detected by subcellular fractionation, FISH, RIP, ChIP and luciferase reporter assays.
It was found that MNX1-AS1 displayed obvious upregulation in GC tissue samples and cell lines, and ectopic expression of MNX1-AS1 predicted poor clinical outcomes for patients with GC. Overexpressed MNX1-AS1 expression promoted proliferation, migration and invasion of GC cells markedly, whereas decreased MNX1-AS1 expression elicited the opposite effects. Consistent with the in vitro results, MNX1-AS1 depletion effectively inhibited the growth of xenograft tumour in vivo. Mechanistically, TEAD4 directly bound the promoter region of MNX1-AS1 and stimulated the transcription of MNX1-AS1. Furthermore, MNX1-AS1 can sponge miR-6785-5p to upregulate the expression of BCL2 in GC cells. Meanwhile, MNX1-AS1 suppressed the transcription of BTG2 by recruiting polycomb repressive complex 2 to BTG2 promoter regions.
Our findings demonstrate that MNX1-AS1 may be able to serve as a prognostic indicator in GC patients and that TEAD4-activatd MNX1-AS1 can promote GC progression through EZH2/BTG2 and miR-6785-5p/BCL2 axes, implicating it as a novel and potent target for the treatment of GC.
Gastric cancer (GC) is the third leading cause of cancer-related death and the fifth most frequently diagnosed cancer worldwide, with 1,000,000 new gastric cancer cases and 783,000 death predicted in 2018 [1, 2]. Notably, the incidence rates are dramatically increased in Eastern Asia, which has the highest incidence level worldwide in both sexes [2, 3]. The majority of newly diagnosed GC patients have reached an advanced stage due to absence of sensitive and specific biomarkers [4, 5]. Consequently, there is an urgent need to further GC pathogenesis and identify useful markers implicated in GC.
Long non-coding RNAs (lncRNAs) are newly-discovered transcripts larger than 200 nucleotides [6,7,8]. Recent research has reported that lncRNAs participate in numerous biological activities, such as apoptosis, cell differentiation and metastasis, thus contributing to cancer progression [9,10,11]. Also, the interplay between lncRNAs and cellular molecules, such as protein, RNA and DNA, can lead to ectopic expression of critical genes [12,13,14,15]. Therefore, studies of potential molecules associated with gastric tumourigenesis represent an attempt to identify the possible application of lncRNAs as effective markers in GC.
The novel lncRNA MNX1-AS1 originally identified as a highly-expressed gene in colorectal cancer (CRC) is also known as CCAT5 [16, 17]. Wang et al. demonstrated that CCAT5 can facilitate CRC progression by interacting with STAT3 . Another research by Ye et al. reported that overexpression of MNX1-AS1 is activated by E2F1 in CRC cells and that MNX1-AS1 can increase SEC61A1 expression by sponging miR-218-5p to facilitate CRC development . In addition, CRC patients in high-MNX1-AS1 expression group exhibited remarkably worse prognostic outcome than those in low-MNX1-AS1 expression group, displaying the robust predictive value of MNX1-AS1 in CRC . Mounting evidence has also shown the presence of MNX1-AS1 overexpression in multiple cancers, including glioblastoma, cervical cancer, GC, ovarian cancer and breast cancer [20,21,22,23]. It is reported that MNX1-AS1 exerts promoting effects on the proliferative capacities of cancer cells in both glioblastoma and ovarian cancer [20, 21] In breast cancer, MNX1-AS1 has been characterized as an oncogene to induce epithelial-mesenchymal transition (EMT) and activate the AKT/mTOR signal pathways . Moreover, MNX1-AS1 facilitated cervical cancer progression through activating mitogen-activated protein kinase (MAPK) pathway . An earlier research reported by Zhang and his colleagues revealed that increased MNX1-AS1 level predicts poor prognosis in patients with GC . Recently, Ma et al. demonstrated that MNX1-AS1 can exert promotive effects on GC migration and invasion by repressing CDKN1A . However, diverse mechanisms behind MNX1-AS1 dysregulation remain largely unknown in GC. The present study is designed to comprehensively investigate MNX1-AS1-mediated mechanistic models in GC progression.
For the present research, RNA sequencing data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets were analysed to find deregulated expression of lncRNAs, which were correlated with GC progression. Among the significantly overexpressed lncRNAs, MNX1-AS1 was selected for further investigation, and the results demonstrated that MNX1-AS1 was significantly amplified in GC and its upregulation frequently predicted poor clinical outcomes. Also, TEA domain family member 4 (TEAD4) activated MNX1-AS1 can mediate the expression level of BTG anti-proliferation factor 2 (BTG2) and BCL2 apoptosis regulator (BCL2) in GC cells. Our findings may illuminate the effectors involved in MNX1-AS1 overexpression and further clarify the underlying mechanism by which MNX1-AS1 mediates the malignant phenotype of GC cells.
LncRNA expression profile analysis
Gene expression data of GC were obtained from TCGA and GEO (GSE62254) datasets. The BAM files and normalized probe-level intensity files were downloaded from TCGA and GEO databases. The probe sequences were downloaded from GEO or microarray manufacturers.
Collection of GC tissues
We collected tissue specimens from 174 GC patients at the People’s Hospital of Jiangsu Province (Nanjing, Jiangsu, China). All the patients signed an informed consents and received no chemotherapy before the operation. 2010 tumour–node–metastasis (TNM) staging recommended by American Joint Committee on Cancer system (AJCC 7th edition) was used to classify tumour stages. This study was approved by the Human Ethics Committee of Nanjing Medical University (Nanjing, Jiangsu, China) and implemented in accordance with the Helsinki Declaration of Principles.
Cell lines culture
We purchased the MGC803, SGC7901, BGC823, HEK-293 T and GES-1 cell lines from Shanghai Cell Bank Library of the Chinese Academy of Sciences (Shanghai, China). SGC7901, SGC803, GES-1 and HEK-293 T cells were cultured in RPMI DMEM (Invitrogen, Shanghai, China) with 10% foetal bovine serum (FBS) at 37 °C in 5% CO2. BGC823 and GES-1 cells were cultured in RPMI1640 medium supplemented with 10% FBS at 37 °C in 5% CO2.
RNA extraction and qRT-PCR assays
Total RNA was isolated with TRIZOL reagent (Invitrogen, Grand Island, NY, USA) according to the product descritption. One microliter total of RNA was reverse-transcribed into cDNA using Primescript RT reagent kit (Takara, Dalian, China). cDNA was quantified by RT-PCR and the data were acquired with SYBR Green (Takara, Dalian, China) using Applied Biosystems 7500 instrument. GAPDH was used as an internal control. The primers are listed in Additional file 4: Table S3.
The plasmid vectors were purified with a DNA Midiprep kit (Qiagen, Hilden, Germany). Three siRNAs and a scrambled RNA construct were designed to prevent the off-target effects and ensure the efficiency of interference. MiRNA mimics and their corresponding miRNA inhibitors were purchased from Ribobio (Guangzhou, China). Cells were transfected using Lipofectamine 3000 (Invitrogen). The sequences of siRNA and shRNA are presented in Additional file 4: Table S3.
Cell migration and invasion assays
The details of cell migration and invasion assays are described in Additional file 5.
Cell proliferation assay
Cell proliferation ability was examined using a Proliferation Reagent Kit I (MTT) (Roche, Basel, Switzerland). Details are available in Supplementary Materials and Methods.
Flow cytometry assay
We analysed cell cycle and apoptosis using flow cytometry assays based on the manufacturer’s protocol. Details of Flow cytometry assay are supplied in Additional file 5.
Tumour formation assay
The programme was approved by the Animal Experimental Ethics Committee of Nanjing Medical University. All experimental procedures involving animals were performed in strict accordance with the NIH Guide for the Care and Use of Laboratory Animals. Details of Tumour formation assay are described in Additional file 5.
Immunohistochemical (IHC) analysis
For hematoxylin-eosin stain (HE) staining, tumour tissues from mice were embedded and sectioned, and then stained with haematoxylin and eosin. For immunohistochemical studies, we incubated the samples with antibodies against Ki67. After washing the sample with PBS, the second antibody was added. The signal was amplified and visualized with 3′-diaminobenzidine chromogen (DAB), and then counterstained with haematoxylin.
An EZMagna RNA Immunoprecipitation (RIP) Kit (Millipore) was used according to the manufacturer’s protocol. Details of the RIP experiment are obtained in Additional file 5. The detailed information of antibodies are listed in Additional file 4: Table S3.
Chromatin immunoprecipitation assays(ChIP)
The TEAD4 binding sequence on the promoter region of MNX1-AS1 was determined by JASPAR (http://jaspar.genereg.net/). The details of luciferase assays are described in Additional file 5. Each experiment was repeated three times.
Fluorescence in situ hybridization (FISH) and subcellular separation
The FISH and subcellular separation assay of SGC-7901 and MGC-803 cells were conducted according to the method described in Additional file 5. The probe sequences are listed in the Additional file 4: Table S3.
The significance of the differences between the groups was evaluated using Student’s t test, Wilcoxon test or χ2 test. The data were analysed with SPSS 17.0 software (IBM), and the double-sided p value was calculated, p values less than 0.05 were considered to indicate a significant difference.
LncRNA MNX1-AS1 is significantly overexpressed and functions AS an unfavourable prognostic factor in patients with GC
To detect lncRNAs that potentially drove gastric tumourigenesis, we analysed the RNA-sequencing data from TCGA and GEO datasets (Fig. 1a, b). Then, lncRNA MNX1-AS1 was chosen for subsequent research. To determine MNX1-AS1 expression, we performed qRT-PCR assays in 174 GC tissues and matched non-tumour tissues (Fig. 1c). Compared with non-cancerous GC tissues, MNX1-AS1 displayed prominent increase in GC tissue samples. Moreover, a significantly higher MNX1-AS1 expression can be seen in GC cells than in the GES-1 cell line (Additional file 1: Figure S1A). Taken together, these results demonstrate a novel dysregulated lncRNA MNX1-AS1 in GC.
Further, we analysed the relationships between MNX1-AS1 level and clinical factors of GC patients. The high-MNX1-AS1 group (n = 87 > median) showed higher tumour stages than the low MNX1-AS1 group (n = 87 < median) (Fig. 1c, d). Additionally, overexpressed MNX1-AS1 expression was obviously associated with tumour size, depth of invasion, histologic grade, TNM stage, lymph node metastasis and distant metastasis in GC patients (Table 1).
As shown in Kaplan-Meier survival curve, GC patients in the high-MNX1-AS1 group had markedly shorter overall survival (OS) and disease-free survival (DFS) rates than those in the low-MNX1-AS1 group (Fig. 1e, f). Besides, factors associated with OS and DFS were evaluated using the univariate and multivariate cox regression models. It was found that tumour size, depth of tumour, lymphatic metastasis, TNM stage and MNX1-AS1 expression appeared to correlate with survival period of GC patients (Additional file 2: Table S1, Additional file 3: Table S2). Importantly, multivariate analysis showed that MNX1-AS1 is an independent prognostic factor for worse OS and DFS among GC patients (Additional file 2: Table S1, Additional file 3: Table S2).
Transcription factor TEAD4 activates lncRNA MNX1-AS1 transcription
Recent research has demonstrated that transcription factors (TFs) can be involved in activating the transcription of some lncRNAs [25,26,27]. To find the transcription factors closely associated with lncRNA MNX1-AS1 overexpression, we further explored the expression data from the GEO dataset (GSE62254) and conducted a correlation analysis between lncRNA MNX1-AS1 and transcription factors. As shown in Fig. 2b and c, TEAD4 was significantly increased in GC tissues and exhibited a positive correlation with lncRNA MNX1-AS1 in GSE62254 (Fig. 2b, c). Thus, we hypothesized that TEAD4 is involved in the transcription that leads to lncRNA MNX1-AS1 overexpression. To confirm this hypothesis, we made analysis of lncRNA MNX1-AS1 promoter using the JASPAR algorithm and found TEAD4-binding site regions (Fig. 2a).
To clarify the upstream regulatory mechanism, we conducted chromatin immunoprecipitation (ChIP) assays and found that TEAD4 directly interacted with the TEAD4 binding sites within the MNX1-AS1 promoter in SGC7901 and MGC803 cells (Fig. 2g). Next, we performed loss- and gain-of-function assays, and detected that the knockdown or overexpression (Fig. 2d, e), decreased or increased TEAD4 enrichment on the MNX1-AS1 promoter respectively (Fig. 2i, h), contributing to MNX1-AS1 downregulation or upregulation in GC cells, respectively (Fig. 2f).
As shown in Fig. 2a, three promoter regions, designated as Luc-MNX1-AS1-F (− 1300 to − 2000), Luc-MNX1-AS1-P1 (− 1500 to − 2000), and Luc-MNX1-AS1-P2 (− 1700 to − 2000), were cloned into a luciferase reporter plasmid to ensure the TEAD4-binding sites (Fig. 2a). Dual-luciferase reporter assays revealed that TEAD4 binds to the E2 (− 1500 to − 1700) region (Fig. 2j). Then, we cloned the E2 promoter region of MNX1-AS1 into a luciferase reporter plasmid and respectively made mutations at 5 putative binding sites to further determine the exact binding sites. HEK293T cells were co-transfected with the pcDNA-TEAD4/Empty-vector and Luc-F, Luc-M1, Luc-M2, Luc-M3, Luc-M4, or Luc-M5. Both Luc-M2 and Luc-M3 were found to markedly reduce promoter activity compared with Luc-F, Luc-M1, Luc-M4, and Luc-M5 (Fig. 2k). These results demonstrated that the TEAD4-binding site cluster is responsible for lncRNA MNX1-AS1 transcription.
lncRNA MNX1-AS1 enhances GC cell proliferation, migration, and invasion and inhibits apoptosis
To clarify the effects of lncRNA MNX1-AS1 dysregulation on GC cell phenotype, we conducted loss-of function and gain-of function assays in GC cells (Fig. 3a, b). MTT and colony formation assays indicated that knockdown of lncRNA MNX1-AS1 impaired GC cell proliferation, whereas overexpression of lncRNA MNX1-AS1 elicited opposite effects (Fig. 3c-g).
To determine the effects of lncRNA MNX1-AS1 on cell cycle regulation, flow cytometry assays were conducted in GC cells. An obvious increase in the G1/G0 phase population of MNX1-AS1-depleted GC cells was observed compared with the respective control groups (Fig. 4a, b). Moreover, SGC7901 and MGC803 cells with MNX1-AS1 knockdown exhibited higher levels of apoptosis than control group cells (Fig. 4c, d). Additionally, Transwell assays indicated that the suppression of MNX1-AS1 attenuated the invasive and migratory capacities of GC cells (Fig. 4e). Conversely, overexpression of MNX1-AS1 exhibited promotion effects on the invasive and migratory capacities of GC cells (Additional file 1: Figure S1B). These findings revealed that lncRNA MNX1-AS1 enhances GC cell proliferation, migration and invasion, induces G1/G0 phase arrest, and activates apoptosis, thus contributing to GC progression.
Inhibition of lncRNA MNX1-AS1 suppressed GC tumourigenesis in vivo
To explore whether lncRNA MNX1-AS1 affects tumour growth in vivo, we injected lncRNA MNX1-AS1 stable knockdown MGC803 cells, or control cells into male nude mice and constructed a xenograft tumour model. As shown in Fig. 5a, tumours formed in the sh-MNX1-AS1 group were dramatically smaller than those formed in the control group (Fig. 5a). In addition, the efficient knockdown of MNX1-AS1 significantly decreased tumour volume and weight compared with control group (Fig. 5b, c, d). Moreover, immunohistochemistry (IHC) staining assays revealed that lncRNA MNX1-AS1 knockdown group had lower ki-67 protein levels (Fig. 5 e). These findings indicated that lncRNA MNX1-AS1 may inhibit tumour growth in vivo.
lncRNA MNX1-AS1 mediated GC progression by repressing the expression of BTG2
To further explore the molecular mechanism about how lncRNA MNX1-AS1 contribute to the progression phenotype of GC cells, we performed RNA-sequencing assays from control or siRNAs against MNX1-AS1 in SGC7901 cells (Fig. 6a, b). GO analysis demonstrated that the potential targets of MNX1-AS1 were involved in various biological processes, such as signal transduction, cell growth and death (Fig. 6c). Then, qRT-PCR assays were used to confirm the changes in the level of possible targets mediated by MNX1-AS1 (Fig. 6d). We chose BTG2 for subsequent research, knowing that BTG2 exhibited the highest expression level in MNX1-AS1-knockdown cells (Fig. 6d). As shown in Fig. 6e, BTG2 was obviously increased in MNX1-AS1 depleted GC cells at mRNA and protein levels (Fig. 6e).
To further clarify the regulatory mode of MNX1-AS1 in GC, we performed FISH and subcellular fractionation assays to study MNX1-AS1 distribution in GC cells. It was shown that MNX1-AS1 mainly localized in the nucleus, implying the role of MNX1-AS1 in transcriptional processing (Fig. 7a). Current evidence has revealed that lncRNAs contribute to cancer phenotypes by interacting with RNA-binding proteins (RBPs) [27,28,29,30]. Interestingly, the prediction results by bioinformatics indicate possible binding between MNX1-AS1 and RBPs, such as Enhancer of zeste 2 (EZH2), SUZ12, EED, AGO2, LSD1 and SIRT (http://pridb.gdcb.iastate.edu//RPISeq/) (Fig. 7b). Subsequently, RNA immunoprecipitation (RIP) experiments demonstrated that MNX1-AS1 bound with polycomb repressive complex 2 (PRC2) and argonaute 2 (Ago2) in both SGC7901 and MGC803 cells (Fig. 7c).
EZH2, a critical member of PRC2, catalyses the demethylation and trimethylation of H3 lysine 27 trimethylation (H3K27me3), thus suppressing expression levels at the transcriptional levels [31, 32]. As shown in Fig. 7d and e, depletion of EZH2 markedly inhibited GC cell viability (Fig. 7d, e). Importantly, BTG2 was remarkably upregulated in EZH2-depleted SGC7901 and MGC803 cells, suggesting that BTG2 may become a prospective downstream target induced by MNX1-AS1 in GC cells (Fig. 7f). To further study whether lncRNA MNX1-AS1 modulates transcription repression of BTG2 through EZH2-mediated H3K27me3, we performed ChIP assays and found that EZH2 bound the promoter region of BTG2 and induced H3K27me3 modification in GC cells. In addition, the binding activity and H3K27me3 levels were obviously decreased in MNX1-AS1-depleted GC cells than those in control cells (Fig. 7g). These findings revealed that lncRNA MNX1-AS1 can suppress the expression level of BTG2 through EZH2-mediated H3K27me3 in GC.
Moreover, rescue experiments were used to study BTG2 involvement in MNX1-AS1-induced contributions to GC cell proliferation, migration and invasion. Therefore, we performed cotransfection in MGC803 cells to reach the purpose of upregulating both MNX1-AS1 and BTG2. MTT assays showed that the cotransfection could partly reverse pcDNA-MNX1-AS1-induced growth (Fig. 7h). The results of transwell experiments revealed that cotransfection with pcDNA-MNX1-AS1 and pcDNA-BTG2 could significantly rescued pcDNA-MNX1-AS1-mediated promoting effects on GC cell invasion and migration (Fig. 7i).
LncRNA MNX1-AS1 functions AS a ceRNA and sponges miR-6785-5p in GC cells
Previous studies have indicated that lncRNAs exert critical roles in regulating gene expression by acting as ceRNAs for miRNAs [33, 34]. As shown in Fig. 7a, approximately 35% of lncRNA MNX1-AS1 was distributed in the cytoplasm in SGC7901 cells and about 21% of MNX1-AS1 was distributed in the cytoplasm in MGC803 cells (Fig. 7a). These findings suggest that lncRNA MNX1-AS1 might also have a post-transcriptional regulation function that contributes to GC progression.
To confirm our hypothesis, we performed RIP assays in GC cells and found that MNX1-AS1 could directly interact with Ago2 (Fig. 7c), which has been identified as a critical member of RNA-induced silencing complex (RISC) and closely correlates with mRNA repression mediated by miRNA. Then, we predicted potential miRNA targeting sites on lncRNA MNX1-AS1 using miRanda, pita, and RNAhybrid (Fig. 8a). Among the miRNAs, we mainly focus on miR-125a-5p [35, 36], miR-760 [37, 38], miR-4728-5p [39, 40] and miR-6785-5p , knowing that they are closely correlated with other malignancies.
To determine the interaction between lncRNA MNX1-AS1 and the listed miRNAs, we performed dual-luciferase reporter assays in HEK-293 T cells transfected with a luciferase plasmid harbouring the sequence of MNX1-AS1 together with plasmids encoding the miRNAs or a control sequence. The results revealed that miR-760, miR-4728-5p and miR-6785-5p significantly inhibited luciferase activity driven by lncRNA MNX1-AS1, and among the listed miRNAs, the suppression ability of miR-6785-5p was the strongest (Fig. 7b). Thus, luciferase reporter constructs were generated to detect binding between MNX1-AS1 and miR-6785-5p (Fig. 8c). It can be observed that the luciferase activities were attenuated in HEK293T cells cotransfected with miR-6785-5p and MNX1-AS1-wile-type (MNX1-AS1-WT), but not the MNX1-AS1-mutation (MNX1-AS1-Mut) (Fig. 8d).
Moreover, qRT-PCR assays showed that miR-6785-5p expression was obviously increased in MNX1-AS1-depleted GC cells (Fig. 8e). Meanwhile, overexpression of MNX1-AS1 significantly decreased miR-6785-5p expression in GC cells (Fig. 8f), whereas miR-6785-5p upregulation failed to cause alteration of MNX1-AS1 expression (Fig. 8g). To further investigate the potential direct binding between lncRNA MNX1-AS1 and miR-6785-5p, RIP assays were performed in GC cells. It was notable that MNX1-AS1 and miR-6785-5p were enriched in immunoprecipitates of Ago2 than those of IgG in both SGC7901 and MGC803 cells (Fig. 8h).
LncRNA MNX1-AS1 modulates the expression levels of BCL2 through post-transcriptional regulation of miR-6785-5p
To further illuminate the biological function of miR-6785-5p in GC, miR-6785-5p mimic or miR-6785-5p inhibitor was used to upregulate or downregulate miR-6785-5p expression in SGC7901 and MGC803 cells (Fig. 9a). MTT assays showed that repression of miR-6785-5p obviously enhanced GC cell viability, while the overexpression of miR-6785-5p led to inverse effects (Fig. 9b). Furthermore, flow cytometric experiments showed that increased miR-6785-5p expression obviously activated apoptosis in GC cells (Fig. 9c).
Interestingly, through RNA-sequencing data analysis, BCL2 was found to be one of the most commonly altered genes (Fig. 6a, b). Consistent with the RNA-sequencing results, the expression level of BCL2 was significantly decreased in GC cells with MNX1-AS1 knockdown at both mRNA and protein levels (Fig. 9d, g). In addition, overexpression or downregulation of miR-6785-5p caused a significant decrease or increase in BCL2 expression at both mRNA and protein levels (Fig. 9e, g). According to these findings, we hypothesized that lncRNA MNX1-AS1, miR-6785-5p and BCL2 might form a common ceRNA network in GC.
To test our hypothesis, we predicted the potential target genes of miR-6785-5p using miRanda and found that BCL2 was likely to bind with miR-6785-5p (Fig. 9f). Indeed, luciferase reporter assays showed that WT-BCL2-driven luciferase expression was significantly inhibited by cotransfection with the miR-6785-5p mimic compared with the control. However, this inhibitory effect was abolished by mutation of the putative miR-6785-5p-binding site in the BCL2 3′ UTR. Taken together, these results indicate that miR-6785-5p can bind to the predicted site in the 3’UTR of BCL2 mRNA and mediated BCL2 expression levels in GC (Fig. 9f). Additionally, we performed rescue assays to evaluate whether lncRNA MNX1-ASl regulates BCL2 by competing for miR-6785-5p. As shown in Fig. 9g, the results showed that overexpression of lncRNA MNX1-AS1 increased BCL2 levels and ectopic expression of miR-6785-5p repressed this increase, whereas suppression of lncRNA MNX1-AS1 decreased BCL2 levels and effects of miR-6785-5p inhibitor impaired this downregulation (Fig. 9g). These findings illuminated the existence of a lncRNA MNX1-AS1-miR-6785-5p-BCL2 regulatory axis. Collectively, BCL2 is a target gene of miR-6785-5p and is indirectly mediated by lncRNA MNX1-AS1 (Fig. 10).
A growing body of research has illuminated that dysregulated lncRNAs may mediate cancer progression [41,42,43]. Recent evidence of the roles of lncRNAs in cancer pathogenesis has added to our knowledge of cancer biology. Although dysregulation of certain lncRNAs in gastric tumourigenesis is a recognized phenomenon, the functional mechanisms of most lncRNAs still remain undetermined in human GC.
Currently, the expression pattern and functional role of lncRNA MNX1-AS1 in diverse tumours have been gradually identified [22, 24, 44, 45] Yang et al. found that overexpression of MNX1-AS1 exerted tumour-promoting roles in lung adenocarcinoma . A current study by Ji et al. showed that enrichment of MNX1-AS1 is strongly linked to TNM stage in hepatocellular carcinoma (HCC) patients . More interestingly, MNX1-AS1 can sponge miR-218-5p to activate the expression of COMMD8, thus facilitating malignant properties of HCC . Previous data provided by Ma et al. showed that impaired MNX1-AS1 expression can regulate CDKN1A expression to suppress invasive capacities of GC cells . To the best of our knowledge, few studies have reported the use of in vitro (let alone in vivo) assays to explore the regulatory correlation between lncRNA MNX1-AS1 and GC, or tried to clarify the upstream regulation of lncRNAs. More studies in epigenetics may shed new light on cancer pathogenesis.
In this study, MNX1-AS1 expression was firstly detected via microarray analysis, and further determined in 174 paired GC tissue samples. Overexpressed MNX1-AS1 promoted GC cell growth and invasion, suggesting its critical value in GC patients. Current evidence has shown that TFs contribute to lncRNA dysregulation in various cancers [25, 47]. TEAD4, a critical member of the TEA domain (TEAD) family, has been identified as an oncogenic regulator in many cancers, including GC [48, 49]. Our data showed that MNX1-AS1 overexpression could be activated by TEAD4.
To clarify the regulatory mechanism of lncRNA MNX1-AS1 in GC, we used RNA-sequencing experiments to investigate potential target genes. It was reported in previous studies that epigenetic alterations such as aberrant histone modifications can participate in cancer development and progression [50,51,52]. Moreover, lncRNAs have been shown to cooperate with chromatin-modifying enzymes, thus activating epigenetic activation or gene silencing . EZH2, an important member of PRC2, could mediate H3K27me3 and exert oncogenic functions through repressing tumour suppressors [53, 54]. For this study, we found that MNX1-AS1 mainly localized in the nucleus, indicating that MNX1-AS1 may exert its oncogenic effects at transcriptional level. Mechanistic investigations elucidated that lncRNA MNX1-AS1 can recruit EZH2 and H3K27me3 to the promoter of BTG2 in the nucleus, thus partially silencing BTG2 expression and mediating oncogenic properties in GC.
Additionally, we also attach great importance to lncRNA-mediated post-transcriptional processing. For instance, our previous research has reported that linc00346 can function as a ceRNA to sponge miR-34a-5p in GC cells . In the present research, both bioinformatic analyses and luciferase reporter assays demonstrated that lncRNA MNX1-AS1 functions as a ceRNA for miR-6785-5p, thereby contributing to the depression of its endogenous target gene. In addition, miR-6785-5p upregulation impaired GC cell proliferation and induce cell apoptosis, while miR-6785-5p knockdown elicited the opposite effects. Interestingly, both RNA-seq data and subsequent verification assays revealed that BCL2 was obviously downregulated in MNX1-AS1-depleted-GC cells. Prediction results from miRanda also revealed that BCL2 was likely to bind with miR-6785-5p. Thus, we proposed the hypothesis that BCL2 may be a candidate target of miR-6785-5p, To confirm this hypothesis, we used luciferase reporter assays and showed that miR-6785-5p targeted BCL2 mRNA at its 3’UTR.
Collectively, we have uncovered that MNX1-AS1 is a novel lncRNA correlated with GC tumourigenesis and progression. TEAD4, a critical oncogenic transcription factor, may be involved in promoting the transcription of lncRNA MNX1-AS1 via binding to the promoter region of MNX1-AS1. In addition, we identified a molecular mechanism underlying GC development that involves lncRNA MNX1-AS1/EZH2/BTG2 and MNX1-AS1/miR-6785-5p/BCL2 axes (Fig. 10).
Our findings may facilitate better understanding of GC pathogenesis and shed new light on lncRNA-based diagnosis and treatment of GC and other human malignancies.
Availability of data and materials
The microarray data used in the study repsitory from TCGA and GEO (GSE62254) datasets. The RNAsequencing raw data and normalized results were submitted to GEO database.
B-cell lymphoma 2
BTG anti proliferation factor 2
Chromatin immunoprecipitation assays
Enhancer of zeste 2
Fetal bovine serum
Fluorescence in situ hybridization
Gene Expression Omnibus
Trimethylation of histone H3 lysine 27
Long non-coding RNAs
Mitogen-activated protein kinase
Mammalian target of rapamycin
Polycomb Repressive Complex
Quantitative reverse transcription polymerase chain reaction
RNA-induced silencing complex
The Cancer Genome Atlas
TEA domain family member4
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68:394–424.
Torre LA, Bray F, Siegel RL, Ferlay J, Lortet-Tieulent J, Jemal A. Global cancer statistics, 2012. CA Cancer J Clin. 2015;65:87–108.
Sugano K. Screening of gastric cancer in Asia. Best Pract Res Clin Gastroenterol. 2015;29:895–905.
Bernards N, Creemers GJ, Nieuwenhuijzen GA, Bosscha K, Pruijt JF, Lemmens VE. No improvement in median survival for patients with metastatic gastric cancer despite increased use of chemotherapy. Ann Oncol. 2013;24:3056–60.
Jiang Y, Xie J, Han Z, Liu W, Xi S, Huang L, Huang W, Lin T, Zhao L, Hu Y, et al. Immunomarker support vector machine classifier for prediction of gastric Cancer survival and adjuvant chemotherapeutic benefit. Clin Cancer Res. 2018;24:5574–84.
Tsai MC, Spitale RC, Chang HY. Long intergenic noncoding RNAs: new links in cancer progression. Cancer Res. 2011;71:3–7.
Schmidt K, Carroll JS, Yee E, Thomas DD, Wert-Lamas L, Neier SC, Sheynkman G, Ritz J, Novina CD. The lncRNA SLNCR recruits the androgen receptor to EGR1-bound genes in melanoma and inhibits expression of tumor suppressor p21. Cell Rep. 2019;27:2493–507.
Lan T, Yuan K, Yan X, Xu L, Liao H, Hao X, Wang J, Liu H, Xie K, Li J, et al. LncRNA SNHG10 facilitates hepatocarcinogenesis and metastasis by modulating its homolog SCARNA13 via a positive feedback loop. Cancer Res. 2019;79:3220-34.
Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154:26–46.
Geisler S, Coller J. RNA in unexpected places: long non-coding RNA functions in diverse cellular contexts. Nat Rev Mol Cell Biol. 2013;14:699–712.
Batista PJ, Chang HY. Long noncoding RNAs: cellular address codes in development and disease. Cell. 2013;152:1298–307.
Ma M, Xu H, Liu G, Wu J, Li C, Wang X, Zhang S, Xu H, Ju S, Cheng W, et al. Metabolism-induced tumor activator 1 (MITA1), an energy stress-inducible long noncoding RNA, promotes hepatocellular carcinoma metastasis. Hepatology. 2019;70:215-30.
Shahabi S, Kumaran V, Castillo J, Cong Z, Nandagopal G, Mullen DJ, Alvarado A, Correa MR, Saizan A, Goel R, et al. LINC00261 is an epigenetically regulated tumor suppressor essential for activation of the DNA damage response. Cancer Res. 2019;79:3050–62.
Qu D, Sun WW, Li L, Ma L, Sun L, Jin X, Li T, Hou W, Wang JH. Long noncoding RNA MALAT1 releases epigenetic silencing of HIV-1 replication by displacing the polycomb repressive complex 2 from binding to the LTR promoter. Nucleic Acids Res. 2019;47:3013–27.
Wang M, Mao C, Ouyang L, Liu Y, Lai W, Liu N, Shi Y, Chen L, Xiao D, Yu F, et al. Long noncoding RNA LINC00336 inhibits ferroptosis in lung cancer by functioning as a competing endogenous RNA. Cell Death Differ. 2019;26:2329-43.
Kim T, Jeon YJ, Cui R, Lee JH, Peng Y, Kim SH, Tili E, Alder H, Croce CM. Role of MYC-regulated long noncoding RNAs in cell cycle regulation and tumorigenesis. J Natl Cancer Inst. 2015;107:505-15.
Kim T, Croce CM. Long noncoding RNAs: Undeciphered cellular codes encrypting keys of colorectal cancer pathogenesis. Cancer Lett. 2018;417:89–95.
Wang Y, Yan XL, Tian SK. Downregulating long non-coding RNA CCAT5 inhibits tumor growth, invasion and metastasis in colorectal cancer through suppressing STAT3. Eur Rev Med Pharmacol Sci. 2019;23:7899–904.
Ye Y, Gu B, Wang Y, Shen S, Huang W. E2F1-mediated MNX1-AS1-miR-218-5p-SEC61A1 feedback loop contributes to the progression of colon adenocarcinoma. J Cell Biochem. 2019;120:6145–53.
Li AH, Zhang HH. Overexpression of lncRNA MNX1-AS1 is associated with poor clinical outcome in epithelial ovarian cancer. Eur Rev Med Pharmacol Sci. 2017;21:5618–23.
Gao Y, Xu Y, Wang J, Yang X, Wen L. Feng J: lncRNA MNX1-AS1 promotes Glioblastoma progression through inhibition of miR-4443. Oncol Res. 2019;27:341–7.
Cheng Y, Pan Y, Pan Y, Wang O. MNX1-AS1 is a functional oncogene that induces EMT and activates the AKT/mTOR pathway and MNX1 in breast cancer. Cancer Manag Res. 2019;11:803–12.
Liu X, Yang Q, Yan J, Zhang X, Zheng M. LncRNA MNX1-AS1 promotes the progression of cervical cancer through activating MAPK pathway. J Cell Biochem. 2019;120:4268–77.
Ma JX, Yang YL, He XY, Pan XM, Wang Z, Qian YW. Long noncoding RNA MNX1-AS1 overexpression promotes the invasion and metastasis of gastric cancer through repressing CDKN1A. Eur Rev Med Pharmacol Sci. 2019;23:4756–62.
Liu HT, Liu S, Liu L, Ma RR, Gao P. EGR1-mediated transcription of lncRNA-HNF1A-AS1 promotes cell-cycle progression in gastric Cancer. Cancer Res. 2018;78:5877–90.
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:136.
Xu TP, Ma P, Wang WY, Shuai Y, Wang YF, Yu T, Xia R, Shu YQ. KLF5 and MYC modulated LINC00346 contributes to gastric cancer progression through acting as a competing endogeous RNA and indicates poor outcome. Cell Death Differ. 2019;26:2179-93.
He F, Song Z, Chen H, Chen Z, Yang P, Li W, Yang Z, Zhang T, Wang F, Wei J, et al. Long noncoding RNA PVT1-214 promotes proliferation and invasion of colorectal cancer by stabilizing Lin28 and interacting with miR-128. Oncogene. 2019;38:164–79.
Sun M, Nie F, Wang Y, Zhang Z, Hou J, He D, Xie M, Xu L, De W, Wang Z, Wang J. LncRNA HOXA11-AS promotes proliferation and invasion of gastric Cancer by scaffolding the chromatin modification factors PRC2, LSD1, and DNMT1. Cancer Res. 2016;76:6299–310.
Klingenberg M, Gross M, Goyal A, Polycarpou-Schwarz M, Miersch T, Ernst AS, Leupold J, Patil N, Warnken U, Allgayer H, et al. The long noncoding RNA Cancer susceptibility 9 and RNA binding protein heterogeneous nuclear Ribonucleoprotein L form a complex and Coregulate genes linked to AKT signaling. Hepatology. 2018;68:1817–32.
Dardenne E, Beltran H, Benelli M, Gayvert K, Berger A, Puca L, Cyrta J, Sboner A, Noorzad Z, MacDonald T, et al. N-Myc induces an EZH2-mediated transcriptional program driving neuroendocrine prostate Cancer. Cancer Cell. 2016;30:563–77.
Xu K, Wu ZJ, Groner AC, He HH, Cai C, Lis RT, Wu X, Stack EC, Loda M, Liu T, et al. EZH2 oncogenic activity in castration-resistant prostate cancer cells is Polycomb-independent. Science. 2012;338:1465–9.
Huang Y, Xu Y, Lu Y, Zhu S, Guo Y, Sun C, Xu L, Chen X, Zhao Y. Yu B, et al: lncRNA Gm10451 regulates PTIP to facilitate iPSCs-derived beta-like cell differentiation by targeting miR-338-3p as a ceRNA. Biomaterials. 2019;216:119266.
Huang X, Xie X, Liu P, Yang L, Chen B, Song C, Tang H, Xie X. Adam12 and lnc015192 act as ceRNAs in breast cancer by regulating miR-34a. Oncogene. 2018;37:6316–26.
Nishida N, Mimori K, Fabbri M, Yokobori T, Sudo T, Tanaka F, Shibata K, Ishii H, Doki Y, Mori M. MicroRNA-125a-5p is an independent prognostic factor in gastric cancer and inhibits the proliferation of human gastric cancer cells in combination with trastuzumab. Clin Cancer Res. 2011;17:2725–33.
Yan L, Yu MC, Gao GL, Liang HW, Zhou XY, Zhu ZT, Zhang CY, Wang YB, Chen X. MiR-125a-5p functions as a tumour suppressor in breast cancer by downregulating BAP1. J Cell Biochem. 2018;119:8773–83.
Cao L, Liu Y, Wang D, Huang L, Li F, Liu J, Zhang C, Shen Z, Gao Q, Yuan W, Zhang Y. MiR-760 suppresses human colorectal cancer growth by targeting BATF3/AP-1/cyclinD1 signaling. J Exp Clin Cancer Res. 2018;37:83.
Zhu L, Xue F, Cui Y, Liu S, Li G, Li J, Guan B, Zeng H, Bian W, Yang C. Zhao C: miR-155-5p and miR-760 mediate radiation therapy suppressed malignancy of non-small cell lung cancer cells. Biofactors. 2019;45:393–400.
Pekow J, Hutchison AL, Meckel K, Harrington K, Deng Z, Talasila N, Rubin DT, Hanauer SB, Hurst R. Umanskiy K, et al: miR-4728-3p functions as a tumor suppressor in ulcerative colitis-associated colorectal Neoplasia through regulation of focal adhesion signaling. Inflamm Bowel Dis. 2017;23:1328–37.
Lulla AR, Slifker MJ, Zhou Y, Lev A, Einarson MB, Dicker DT. El-Deiry WS: miR-6883 family miRNAs target CDK4/6 to induce G1 phase cell-cycle arrest in Colon Cancer cells. Cancer Res. 2017;77:6902–13.
Hu Q, Ye Y, Chan LC, Li Y, Liang K, Lin A, Egranov SD, Zhang Y, Xia W, Gong J, et al. Oncogenic lncRNA downregulates cancer cell antigen presentation and intrinsic tumor suppression. Nat Immunol. 2019;20:835–51.
Zhuo W, Liu Y, Li S, Guo D, Sun Q, Jin J, Rao X, Li M, Sun M, Jiang M, et al. Long noncoding RNA GMAN, up-regulated in gastric Cancer tissues, is associated with metastasis in patients and promotes translation of Ephrin A1 by competitively binding GMAN-AS. Gastroenterology. 2019;156:676–91.
Zhang Y, Pitchiaya S, Cieslik M, Niknafs YS, Tien JC, Hosono Y, et al. Analysis of the androgen receptor-regulated lncRNA landscape identifies a role for ARLNC1 in prostate cancer progression. Nat Genet. 2018;50:814-24.
Zhang W, Huang L, Lu X, Wang K, Ning X, Liu Z. Upregulated expression of MNX1-AS1 long noncoding RNA predicts poor prognosis in gastric cancer. Bosn J Basic Med Sci. 2019;19:164–71.
Yang R, Wang L, Han M. MNX1-AS1 is a novel biomarker for predicting clinical progression and poor prognosis in lung adenocarcinoma. J Cell Biochem. 2018:120:7222-28.
Ji D, Wang Y, Sun B, Yang J, Luo X. Long non-coding RNA MNX1-AS1 promotes hepatocellular carcinoma proliferation and invasion through targeting miR-218-5p/COMMD8 axis. Biochem Biophys Res Commun. 2019;513:669–74.
Melo CA, Leveille N, Rooijers K, Wijchers PJ, Geeven G, Tal A, Melo SA, de Laat W, Agami R. A p53-bound enhancer region controls a long intergenic noncoding RNA required for p53 stress response. Oncogene. 2016;35:4399–406.
Liu Y, Wang G, Yang Y, Mei Z, Liang Z, Cui A, Wu T, Liu CY, Cui L. Increased TEAD4 expression and nuclear localization in colorectal cancer promote epithelial-mesenchymal transition and metastasis in a YAP-independent manner. Oncogene. 2016;35:2789–800.
Zhou Y, Huang T, Zhang J, Wong CC, Zhang B, Dong Y, Wu F, Tong J, Wu W, Cheng A, et al. TEAD1/4 exerts oncogenic role and is negatively regulated by miR-4269 in gastric tumorigenesis. Oncogene. 2017;36:6518–30.
Okugawa Y, Grady WM, Goel A. Epigenetic alterations in colorectal Cancer: emerging biomarkers. Gastroenterology. 2015;149:1204–25.
Herman JG, Baylin SB. Gene silencing in cancer in association with promoter hypermethylation. N Engl J Med. 2003;349:2042–54.
Weber M, Hellmann I, Stadler MB, Ramos L, Paabo S, Rebhan M, Schubeler D. Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007;39:457–66.
Frapporti A, Miro PC, Arnaiz O, Holoch D, Kawaguchi T, Humbert A, Eleftheriou E, Lombard B, Loew D, Sperling L, et al. The Polycomb protein Ezl1 mediates H3K9 and H3K27 methylation to repress transposable elements in paramecium. Nat Commun. 2019;10:2710.
Gentile C, Berlivet S, Mayran A, Paquette D, Guerard-Millet F, Bajon E, Dostie J, Kmita M. PRC2-associated chromatin contacts in the developing limb reveal a possible mechanism for the atypical role of PRC2 in HoxA gene expression. Dev Cell. 2019;50:184-96.
The work was supported by the National Natural Science Foundation of China (Grant number: 81672896/81602071) and Innovation Fund for Outstanding Doctoral Candidates of Peking University Health Science Center (Grant number: 71013Y2029).
Ethics approval and consent to participate
This study was approved by the Ethics Committee of Nanjing Medical University (Nanjing, Jiangsu, China) and implemented in accordance with the Helsinki Declaration of Principles.
Consent for publication
All the authors agree to publish this paper.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. LncRNA MNX1-AS1 was significantly upregulated in GC cells, and its ectopic expression promoted GC cell migration and invasion. A.The expression level of MNX1-AS1 in GC cells and GES-1 cells. B. The effects of MNX1-AS1 overexpression on GC cell migration and invasion.
. Univariate and multivariate analysis of clinicopathological factors for overall-survival in gastric cancer patients (n = 174).
. Univariate and multivariate analysis of clinicopathological factors for disease-free survival in gastric cancer patients (n = 174).
. Primers for qRT-PCR, siRNAs oligonucleotides and the company for antibody.
Materials and Methods.
About this article
Cite this article
Shuai, Y., Ma, Z., Liu, W. et al. TEAD4 modulated LncRNA MNX1-AS1 contributes to gastric cancer progression partly through suppressing BTG2 and activating BCL2. Mol Cancer 19, 6 (2020). https://doi.org/10.1186/s12943-019-1104-1
- Gastric cancer