- 1Department of Gastrointestinal Surgery, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 2Laboratory of General Surgery, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 3Digestive Disease Center, The Seventh Affiliated Hospital, Sun Yat-sen University, Shenzhen, China
Long noncoding RNAs (lncRNAs) have been reported to regulate diverse tumorigenic processes. However, little is known about long intergenic non-protein coding RNA 00893 (LINC00893) and its role in gastric cancer (GC). Herein we investigated its biological functions and molecular mechanism in GC. LINC00893 was decreased in GC tissues but significantly elevated in AGS cells after treatment with Nutlin-3. In GC patients, it was found that low expression of LINC00893 was correlated with tumor growth, metastasis and poor survival. Functionally, overexpression of LINC00893 suppressed the proliferation, migration and invasion of GC cells. Mechanistically, LINC00893 regulated the expression of epithelial-mesenchymal transition (EMT)-related proteins by binding to RNA binding fox-1 homolog 2 (RBFOX2) and promoting its ubiquitin-mediated degradation, thus suppressing the EMT and related functions of GC. In addition, the transcription factor p53 can regulate the expression of LINC00893 in an indirect way. Taken together, these results suggested that LINC00893 regulated by p53 repressed GC proliferation, migration and invasion by functioning as a binding site for RBFOX2 to regulate its stability and the expression of EMT-related proteins. LINC00893 acts as a tumor-inhibiting lncRNA that is induced by p53 in GC and regulates EMT by binding to RBFOX2, thus providing a novel experimental basis for the clinical treatment of GC.
Introduction
Gastric cancer (GC) is one of the leading causes of cancer-related human mortalities, and remains the fourth leading cause of cancer-related deaths in 2020 (Sung et al., 2021). Strikingly, GC cases in China account for about 42% of all cases worldwide, mainly owing to the high prevalence of Helicobacter pylori infection (Cai et al., 2021). Moreover, although massive efforts have been made by researchers to develop noninvasive biomarkers and therapeutic targets to reduce GC mortality, little success has been achieved in increasing the overall survival rate. Nonresectable or metastatic GC is always associated with poor prognosis, and patients can get minimal benefit from systemic chemotherapeutic approaches (Mehta et al., 2020). Therefore, novel biomarkers for improving the early diagnosis, tumor grading and prognostic evaluation of GC are urgently needed.
Owing to the rapid development in high-throughput sequencing, it has been determined that nonprotein coding RNA accounts for the vast majority of RNA transcribed from human genome, which indicates that a relatively small number of effectors is dedicated to be regulated by a large group of RNA regulators (Katayama et al., 2005; Djebali et al., 2012). Among these newly discovered nonprotein RNA elements, long noncoding RNAs (lncRNAs) have been widely reported to function as key regulators in multiple cellular processes, such as cell differentiation, development, and disease pathogenesis (Kung et al., 2013). To date, the expression profile of lncRNAs in GC have been screened by different groups and several disordered lncRNAs related to GC carcinogenesis have been investigated, such as MALAT1, HOTAIR, and UCA1(Liu et al., 2014; Li et al., 2017; Wang et al., 2019). LncRNAs are critical regulators of GC initiation, progression, dissemination and immune evasion (Dragomir et al., 2020). It has also been widely proved that there is a close connection between lncRNAs and gastric cancer metastasis. Some of them are reported to promote the metastasis of gastric cancer, such as lncRNA UBE2CP3, HOTAIR and GMAN (Wei et al., 2020; Xu et al., 2020; Li et al., 2021). Meanwhile, lncRNAs such as lnc-CTSLP4 and lnc-LEMGC can inhibit the metastasis of gastric cancer through different ways (Pan et al., 2021; Zhang et al., 2022).
The gene TP53, encoding p53 protein, is the most frequently mutated tumor suppressor gene in diverse human cancers and, therefore, one of the most broadly investigated regulatory networks (Vousden and Prives, 2009; Allen et al., 2014). For decades, when studying the regulatory mechanisms and downstream targets which execute the biological functions of p53, researchers mainly focused on the products of protein-coding genes (Donehower et al., 2019). Recently, it was discovered that lncRNAs were key regulatory players which could shape p53 activity and biological outcomes. The p53 regulatory networks are now reported to be intervened by multitudinous p53-regulated lncRNAs either directly or indirectly in cancers (Jain, 2020). LincRNA-p21, GUARDIN, NEAT1 and TP53TG1 are transcriptionally regulated by p53 (Huarte et al., 2010; Diaz-Lagares et al., 2016; Mello et al., 2017; Hu et al., 2018), while DDSR1 and LINP1 are induced by p53 in indirect ways (Sharma et al., 2015; Zhang et al., 2016).
In this study, with bioinformatic and statistical analysis of high-throughput sequencing data, we found that the expression of LINC00893 was significantly reduced in GC tissues but increased in AGS cells after increasing p53, so we further investigated LINC00893 in GC. We discovered that LINC00893 was regulated by p53 in an indirect way, while low expression of LINC00893 was positively associated with GC growth, metastasis and indicated a poor survival of GC patients. We further revealed that LINC00893 inhibited proliferation, metastasis and EMT processes of GC cells through binding to the protein RBFOX2 and promoting its degradation. Our new findings on the p53/LINC00893/RBFOX2 regulatory axis illustrate its potential as a target for GC diagnostic and therapeutic development.
Materials and Methods
Cell Culture and Cell Transfection
The GC cell lines AGS, MKN28 and MKN45 were obtained from the Type Culture Collection of the Chinese Academy of Sciences (China). The normal gastric epithelial cell line GES-1 was purchased from the Lab Animal Center of the Fourth Military Medical University (China). All cell lines were routinely tested and had negative results for mycoplasma. All cells except AGS were maintained in DMEM (Invitrogen) and AGS in F12 (Invitrogen) supplemented with 10% fetal bovine serum (GIBCO) in an incubator with 5% CO2 and 95% room air at 37°C. The LINC00893 gene was cloned into the pEZ-Lv201 lentivirus vector (GeneCopoeia), RBFOX2 gene was cloned into the pcDNA3.1 vector (GeneCopoeia), and wild-type TP53 were cloned into the pEnter lentivirus vector (Vigene Biosciences). All the small interfering RNAs (siRNA) and oligonucleotides (ASO) were purchased from RiboBio (China). All plasmids, siRNA and ASO were transfected with Lipofectamine 3000 reagent (Invitrogen) following the instructions. The sequences of siRNAs and ASO used are listed in Supplementary Table S5.
RNA Sequencing and Bioinformatics Analysis
RNA sequencing was performed using an Illumina HiSeqTM 2500. Sample preparation was performed based on the manufacturer’s standard protocols. The cDNA/DNA/small RNA libraries were sequenced on the Illumina sequencing platform by Genedenovo Biotechnology Co., Ltd. (Guangzhou, China). The public datasets analyzed in the current study are available in TCGA (https://cancergenome.nih.gov/), GEPIA (http://gepia.cancer-pku.cn/), and KM Plotter (http://kmplot.com/analysis/). GSEA 3.0 (http://www.gsea-msigdb.org/gsea/index.jsp) was used to perform Gene set enrichment analysis (GSEA).
LncRNA Coding Capacity Prediction
Strand-specific genomic coordinates for all exons of human LINC00893, PINCR, NEAT1, GAPDH, SDHA and UBC genes were downloaded from the UCSC Genome Browser (GRCh38/hg38) in BED format. PhyloCSF was applied to generate FASTA alignments to assess the coding potential (the codon substitution frequencies score - CSF) of individual exons and of mature transcripts of the above genes. The Coding Potential Assessment Tool (CAPT, http://lilab. research.bcm.edu/cpat/) constructs a logistic regression model built with four sequence features: open reading frame coverage, open reading frame size, hexamer usage bias and the Fickett TESTCODE statistic. The cutoff value for human coding probability (CP) was set at 0.364. CP < 0.364 was considered a noncoding transcript, whereas CP ≥ 0.364 was defined as a coding transcript.
Real-Time Quantitative Reverse Transcription PCR (RT-qPCR)
RNA isolation plus (TaKaRa, Japan) was used to extract total RNA from tumor tissues or cells. Single-stranded cDNA was generated from 1 μg of total RNA in a 20 μL reaction volume system using oligodT primers according to the protocol supplied with Primer ScriptTM RT Reagent (TaKaRa, Japan). The relative expression levels were measured by qPCR using an ABI 7900HT instrument (Applied Biosystems, United States) in a total volume of 10 µL with the SYBR green detection system (Takara, Japan), and GAPDH was used as an endogenous control. The primers used for the qPCR analysis are listed in Supplementary Table S5.
Proliferation Assay
Cell viability was examined by Cell Counting Kit-8 (CCK-8; Dojindo Laboratories, Japan). Briefly, after treatments, CCK-8 solution (10 μL) was added to each well, and the plates were incubated for an additional 2 h. The absorbance of each well was measured at 450 nm using a microplate reader (MultiskanEX, Lab systems, Finland). In the colony formation assay, treated cells (1×103 cells/well) were seeded into 6-well plates and cultured in DMEM/F12. Then the medium was replaced every 3 days. After 14 days, colonies were fixed with paraformaldehyde and stained with 0.1% crystal. After that, colonies were counted.
Transwell Assays
Firstly, GC cells were transfected as described above for 48 h. Next, cells were harvested and counted. Then a certain number of transfected cells (5×104 cells/well) were cultured with serum-free medium in the upper chamber of Transwell plates (Corning, United States). Medium containing fifteen percent FBS was added into the lower chamber of the Transwell plates. For the invasion experiments, the upper chamber was covered with DMEM (F12) and Matrigel (BD Biosciences, CA) mixture in advance. Finally, cells on the top of the chamber were wiped with cotton swabs, while cells that went through the membrane were stained with 0.1% crystal violet, observed and counted under microscope at ×100 magnification.
In vivo Assays
Four-week-old female athymic BALB/c nude mice were fostered under specific pathogen-free conditions and treated according to protocols. Lv201-LINC00893 or control vector were stably transfected into MKN45 cells. To assess the cell proliferation viability in different groups in vivo, the stable overexpressing cells were harvested, washed with cold PBS and resuspended to 5 × 107 cells/ml. Then, 100 μL cell suspension was subcutaneously injected into the armpits of each mouse. Tumor growth was examined and recorded every week, and tumor volumes were calculated using equation V = 0.5 × D × d2 (V, volume; D, longitudinal diameter; d, latitudinal diameter). Four weeks after injection, mice were sacrificed, and two groups of subcutaneous tumors were obtained and imaged separately. In the cell metastasis assay in vivo, the transfected cells were resuspended at a concentration of 1 × 107 cells/ml, and then 100 μL of resuspended cells was injected into the tail veins of the two groups of mice. Eight weeks after injection, the mice were euthanized. Lungs of the mice in two groups were removed and photographed. We then performed hematoxylin and eosin (HE) staining after tissue fixation and sectioning to evaluate the metastatic nodes, besides, the number of metastatic nodes in each lung was counted. The animal experiments were conducted in accordance with the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. Our whole protocol was approved by the Committee on the Ethics of Animal Experiments of the First Affiliated Hospital of Sun Yat-sen University.
Nuclear and Cytoplasmic Extraction
Cytoplasmic and nuclear fractions of RNA were isolated in accordance with the instructions of PARIS™ Kit (AM1556, Thermo Fisher Scientific, United States). Briefly, AGS and MKN28 cells were washed with cold PBS and lysed in Cell Fraction Buffer on ice for 10 min. Then gathering the fractions into new tubes and centrifugating at 500 × g for 3 min at 4°C. After that, the supernatant was collected as the cytoplasmic fraction. Finally, washing the rest pellet with Cell Fraction Buffer and collecting the nuclei.
Chromatin Isolation by RNA Purification (CHIRP)
RNA was extracted from AGS cells or transcribed in vitro and purified. Purified RNA was then incubated with specific prewashed beads at 37°C for 30 min. After centrifugation, the supernatant was hybridized with denatured biotin-labeled RNA probes at 55°C for 2 h. Next, we added prewashed beads, incubated them at 37°C for 1 h, centrifuged them, and washed the beads. The pulldown RNA was extracted and analyzed by quantitative reverse transcription PCR (RT-qPCR). LINC00893 RNAs were purified, and equal amounts of transcripts were used in the binding assays. The binding components were pulled down by biotin-labeled RNA probes and analyzed by western blotting.
RNA Immunoprecipitation (RIP)
RNA immunoprecipitation (RIP) assays were performed using a Magna RIP RNA-Binding Protein Immunoprecipitation Kit (Millipore, United States) according to its instructions. Anti-RBFOX2 (Bethyl, United States) and normal Rabbit IgG (CST, United States) were used to immunoprecipitated with target RNA or as a negative control. Finally, RT-qPCR analysis was performed to detect the enrichment of the coprecipitated RNA binding to RBFOX2.
Fluorescence in situ Hybridization (FISH)
Cy3 labeled LINC00893 probe was designed and synthesized by GenePharma (Shanghai, China). For FISH assay, AGS or MKN28 cells (1 × 105) were fixed in 4% formaldehyde and permeabilized with 0.3% Triton X-100 for 15 min, washed with PBS three times and once in 2× SSC buffer. Hybridization was carried out using DNA probe sets at 37 °C for 16 h. Images were obtained with confocal laser scanning microscope (Zeiss, Germany) and processed using the ZEN imaging software.
Western Blotting Analysis
Total proteins were extracted using RIPA supplemented with protease and phosphatase inhibitor reagents (Thermo-Fisher Scientific, United States). Equal amounts of protein samples were separated by SDS/PAGE and then electrophoretically transferred to PVDF membranes (Millipore, United States). Next, the membranes were probed with specific primary antibodies against p53 (dilution, 1:1,000; Abcam, United States), p21 (dilution, 1:1,000; CST, United States), RBFOX2 (dilution, 1:1,000; Bethyl, United States), E-cadherin (dilution, 1:1,000; CST, United States), N-cadherin (dilution, 1:1,000; CST, United States), vimentin (dilution, 1:1,000; CST, United States) and GAPDH (dilution, 1:5,000; Proteintech, China) at 4°C overnight, followed by incubation with HRP-conjugated secondary antibodies (dilution, 1:5,000; CST, United States) at room temperature for 1 h. Immunoreactive bands were visualized using enhanced chemiluminescence reagents (Bio-Rad, United States), and GAPDH was considered an inner loading control.
Statistical Analysis
All data are presented as the mean ± standard deviation, appropriate statistical methods including Student’s t-test, Wilcoxon signed-rank test, Mann-Whitney test, Pearson chi-square test were used to calculate differences between groups, and p < 0.05 indicates statistical significance. Statistical analyses were performed using SPSS Statistics software (version 18.0) or GraphPad Prism 7 software.
Results
Decreased LINC00893 Predicts Poorer Survival in GC Patients
First, to screen out the potential tumor suppressor lncRNAs in GC, we performed RNA sequencing of three pairs of human GC tissues and matched nontumor tissues. Differentially expressed lncRNA transcripts were identified. A total of 409 lncRNA transcripts were found to have significant decreases in GC tissues vs nontumor tissues (Fold change<0.5, p < 0.05), while 385 were found to have significant increases. (Fold change>2, p < 0.05) (Figure 1A; Supplementary Table S1). Nutlin-3 is an inhibitor of human homolog of murine double minute 2 (MDM2), which is a negative regulator of p53 (Endo et al., 2011). AGS is a TP53-wild type GC cell line, in which nutlin-3 has been proved to significantly elevate the expression of p53 (Wei et al., 2010; Vaseva et al., 2011). To determine the optimal concentration and time of nutlin-3 stimulation, we treated AGS cells with gradient concentrations and time periods of nutlin-3 and then detected the expression of p53 and its downstream effector p21 (Supplementary Figure S1A–D). Afterwards, we treated AGS cells with nutlin-3 at the most appropriate concentration (10 µM) and time (12 h) and then conducted RNA sequencing to identify the differentially expressed lncRNA transcripts (Figures 1B,C). As the cluster map shows, 152 transcripts were downregulated when AGS was stimulated with nutlin-3 (Fold change<0.5, p < 0.05), while 146 transcripts were upregulated, which were considered potential targets of p53 (Fold change>2, p < 0.05) (Figure 1D; Supplementary Table S2). Next, we cross-analyzed the lncRNAs differentially expressed in gastric tissues and AGS cells elevated after the accumulation of p53, and we found there were 40 lncRNAs differentially expressed in both sequencing data (Figure 1E). Intriguingly, there were only four lncRNAs decreased in GC tissues (Fold change<0.5, p < 0.05) and, at the same time, upregulated in AGS cells treated with nutlin-3 (Fold change>2, p < 0.05) (Figure 1F). This might be due to the highly heterogeneous nature of GC tissues and huge differences of the lncRNA expression profiles between GC tissues and AGS cells. We used GEPIA (Tang et al., 2017) to analyze the expression of these four candidate lncRNAs (Figure 2A, Supplementary Figure S1E). In addition, RT-qPCR confirmed that the level of LINC00893 was significantly decreased in fresh gastric tumor tissues compared with nontumor tissues, while there were no differences in the other three lncRNAs (Figure 2B, Supplementary Figure S1F). We also found there were no significant differences between the expression of these three lncRNAs and the survival of GC patients (Supplementary Figure S1G). Therefore, we conducted further research on LINC00893. The noncoding capability of LINC00893 was confirmed by analysis with PhyloCSF (Lin et al., 2011) and the Coding Potential Assessing Tool (CPAT) (Duan et al., 2021) (Figure 2C, Supplementary Figure S1H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g001.jpg)
FIGURE 1. Identification of p53-regulated tumor suppressor lncRNAs in gastric cancer (A) Heat map is shown for the differentially expressed lncRNAs identified by RNA sequencing between three gastric cancer tissues and paired normal tissues (B,C) Isogenic TP53-WT AGS cells were untreated or treated with Nutlin-3 for 12 h, RT-qPCR and immunoblotting was performed for p53 and the loading control GAPDH (D) Heat map is shown for the differentially expressed lncRNAs identified by RNA sequencing in duplicate from AGS cells untreated or treated with nutlin-3 for 12 h. Upregulated genes are shown in red and downregulated genes in green (E,F) Venn diagram showing the overlap between the transcriptomes downregulated in gastric cancer tissues and the transcriptomes upregulated after nutlin-3 treatment of AGS cells.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g002.jpg)
FIGURE 2. LINC00893 is significantly low-expressed in gastric cancer and regulated by p53 (A) Analysis of LINC00893 expression in unpaired GC (T = 408) and normal tissues (N = 211) in the GEPIA (p < 0.05) (B) Relative LINC00893 expression in 50 paired fresh frozen normal gastric tissues and gastric cancer tissues (C) Maximum CSF scores of LINC00893 as well as other coding and noncoding RNAs determined by analysis with PhyloCSF (D) Kaplan-Meier survival curve indicated the lower expression of LINC00893 associated with poorer survival rates (E,F) p53 and p21 expression in AGS treated with gradient nutlin-3 or 5-FU (G,H) p53 and p21 expression in AGS after transfection of TP53 siRNA or overexpression plasmid (I,J) LINC00893 expression in AGS treated with gradient nutlin-3 or 5-FU (K,L) LINC00893 expression in AGS after transfection of TP53 siRNA or overexpression plasmid. *p < 0.05, ****p < 0.0001, all RT-qPCR results were normalized to GAPDH.
In addition, the correlation of the expression of LINC00893 and the clinicopathological characteristics of patients was analyzed. As shown in Table 1, LINC00893 expression was negatively associated with tumor size, primary tumor (T) stage, lymph node metastasis, distant metastasis and tumor-node-metastasis (TNM) stage. Further analysis indicated that low expression of LINC00893 was negatively correlated with worse survival of GC patients (Figure 2D), suggesting that LINC00893 may inhibit GC progression.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-t001.jpg)
TABLE 1. Correlation between LINC00893 expression and clinicopathological characteristics in GC patients.
LINC00893 Is Regulated by p53 Indirectly in GC Cells
To delineate the possible link to p53, we looked for a correlation between LINC00893 and p53. First, we found that nutlin-3 treatment significantly increased the expression of LINC00893 in AGS cells, which express wild-type (WT) p53 (Figures 2E,I). Next, we treated AGS cells with the chemotherapy drug 5-fluorouracil (5-FU), which could cause DNA damage and increase the level of p53. We found that LINC00893 increased with higher concentration gradients, consistent with p53 (Figures 2F,J). To rule out other effects of nutlin-3 and 5-FU, we transfected TP53 plasmids or siRNAs and verified the expression of p53 (Figures 2G,H). With the overexpression of p53, LINC00893 was significantly elevated, while it was reduced after the knockdown of p53 (Figures 2K,L). Despite its p53-dependent induction, we were unable to identify an obvious p53 binding site in the vicinity of the LINC00893 promoter, nor was one identified in a recent p53 ChIP-seq study performed in TP53 WT cells (Sanchez et al., 2014), thus suggesting that p53 may regulate LINC00893 expression via an indirect pathway, possibly through other p53 targets associated with transcription or p53 binding to enhancer regions. To identify promoter elements involved in p53-mediated regulation of LINC00893, we then analyzed the promoter of LINC00893 and predicted the potential transcription factors (TFs). Among the predicted TFs, Sp1 transcription factor (SP1) was predicted both in GeneCards and PROMO (Farre et al., 2003). Intriguingly, SP1 has been proved to be regulated by p53 and act as an intermediary TF between p53 and its target gene (Innocente and Lee, 2005; Li et al., 2014a). Therefore, we speculated that p53 regulates LINC00893 expression through p53/SP1/LINC00893 axis.
LINC00893 Inhibits GC Proliferation, Migration and Invasion in vitro
According to the correlation between LINC00893 expression and GC, we focused on the biological functions of LINC00893 in GC proliferation and metastasis. AGS is an epithelial cell line derived from gastric adenocarcinoma and MKN28 is a highly metastatic GC cell line, which are widely used as models to detect the roles of oncogene or tumor suppressor gene in GC in vitro. The expression vector Lv201-LINC00893 or ASO-LINC00893 was transfected into GC cells, and the efficiency of LINC00893 overexpression or knockdown was confirmed by RT-qPCR (Figure 3A). CCK8 and colony formation assays were conducted, and they showed strikingly inhibited proliferation ability when LINC00893 was overexpressed or a promoting proliferation ability when LINC00893 was knocked down (Figures 3B,C). Transwell assays showed that LINC00893 overexpression significantly decreased the migration and invasion potential of GC cells, while LINC00893 knockdown showed the opposite effects (Figures 3D,E). To explore the downstream signaling pathways involved in LINC00893, we performed RNA-seq and bioinformatic analysis in AGS cells treated with empty vector or LINC00893 overexpression vector. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment and Gene Set Enrichment Analysis (GSEA) results indicated that LINC00893 was correlated with pathways directly involved in the regulation of cell adherens junctions (Supplementary Table S3; Figure 3F; Supplementary Figure S1I). Therefore, we detected the key representative proteins of EMT. We found that overexpression of LINC00893 increased the level of E-cadherin and decreased the levels of N-cadherin and vimentin. The opposite effect was confirmed after knockdown of LINC00893 (Figure 3G). This suggests that LINC00893 can inhibit the process of EMT in GC cells.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g003.jpg)
FIGURE 3. LINC00893 inhibits GC cell proliferation, migration, invasion and EMT in vitro (A) LINC00893 was overexpressed or knocked down in AGS and MKN28 cells (B,C) The effects of LINC00893 overexpression or knockdown on proliferation and plate colony-forming ability were measured in GC cells (D,E) The effects of LINC00893 overexpression or knockdown on migration and invasion were detected in GC cells (F) RNA-seq and Gene Set Enrichment Analysis (GSEA) of AGS treated with vector or LINC00893 overexpressed plasmid (G) Relative expression of E-cadherin, N-cadherin and Vimentin in GC cells after LINC00893 were either overexpressed or knocked down. *p < 0.05, **p < 0.01, ***p < 0.001.
LINC00893 Suppresses GC Progression in vivo
To determine the roles of LINC00893 in GC tumorigenesis in vivo, MKN45, a GC cell line derived from gastric lymph nodes in a female patient with signet ring cell carcinoma, was used to construct model cells stably overexpressing LINC00893 or the control vector. Then, cells were subcutaneously injected into nude mice (Figure 4A). We found mice in the LINC00893-overexpressing group generated smaller tumors than the control group (Figure 4B), and HE staining was performed to identify the tumors (Figure 4C left panel). Tumor volume and weight also indicated that the overexpression of LINC00893 suppressed tumor growth in vivo (Figures 4D,E). Besides, a significant difference in Ki-67 staining was found between the LINC00893 and control group (Figure 4C right panel). Furthermore, by injecting the stably transfected MKN45 cells into tail veins of nude mice, it was found that LINC00893 dramatically inhibited the development of pulmonary metastasis (Figures 4F–H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g004.jpg)
FIGURE 4. LINC00893 inhibits the growth and metastasis of gastric cancer cells in vivo (A) LINC00893 stable-overexpressed MKN45 cells was constructed (B) Representative images of subcutaneous xenograft tumors from the two group (C) Representative images of H&E and IHC with anti-Ki67 of the subcutaneous xenograft tumors (D) Quantification of tumor volumes in NC or LINC00893 overexpressed MKN45 cells in xenograft mouse models (E) Tumor weights were analyzed (F) Representative images of lung metastasis of LINC00893 overexpression groups and control group (G) Representative hematoxylin and eosin (H&E) staining results of corresponding lung metastatic nodules (H) Statistical analysis of numbers of metastatic nodules in the lung. ***p < 0.001.
LINC00893 Directly Binds to RBFOX2 and Decreases Its Expression
Several recent studies have revealed the closely related regulatory network between lncRNAs and proteins. GO (Gene Ontology) enrichment indicated that LINC00893 might play a role though binding to proteins (Supplementary Figure S2B). First, to identify the subcellular localization of LINC00893, we detected LINC00893 expression in cytoplasmic and nuclear fractions by RT-qPCR analysis (Figure 5A). The results showed that LINC00893 was localized predominantly in the nucleus, with some localization in the cytoplasm, which was further verified by RNA fluorescence in situ hybridization assays (Figure 5B). The secondary structure of LINC00893 was also predicted (Figure 5C). By using the StarBase database (Li et al., 2014), a large number of proteins were predicted to directly bind to LINC00893 (Supplementary Table S4). The splicing factor RNA Binding Protein Fox-1, Homolog 2 (RBFOX2) was included, which is a classic nuclear regulator of Polycomb Repressive Complex 2 (PRC2) recruitment to active genes, particularly those with bivalent histone modifications, to mediate dynamic transcriptional control and regulate autoregulation in RNA-binding protein networks (Venables et al., 2013; Braeutigam et al., 2014; Jangi et al., 2014; Wei et al., 2016). We found that nucleotides 753-848 of the LINC00893 sequence and amino acids 126-177 of the RBFOX2 sequence were the most likely binding areas (Figure 5D). To verify the direct interaction between LINC00893 and RBFOX2, western blotting was carried out using antibodies specific for RBFOX2 to test the protein samples precipitated by the LINC00893 CHIRP assay. This result indicated that RBFOX2 was pulled down by LINC00893 (Figures 5E,F). To further confirm this interaction, we carried out an RIP assay using a RBFOX2-specific antibody, followed by RT-qPCR using primers specific for LINC00893. As expected, LINC00893 was enriched in the anti-RBFOX2 group compared to the control IgG group (Figure 5G). These results indicated that RBFOX2 directly interacts with LINC00893.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g005.jpg)
FIGURE 5. Identification of RBFOX2 as a binding partner for LINC00893 (A) RT-qPCR analysis from nuclear and cytoplasmic fractions of GC cells; the cytoplasmic GAPDH mRNA and the nuclear lncRNA MALAT1 were used as controls (B) Subcellular localization of LINC00893 detected by FISH. Scale bar: 20 μm (C) The secondary structure of LINC00893 was predicted (http://rna.tbi.univie.ac.at). The red color indicates strong confidence for the prediction of each base (D) The potential binding areas between LINC00893 and RBFOX2 were predicted using the catRAPID database (E) SDS-PAGE of proteins purified from CHIRP assay using biotinylated LINC00893 or antisense RNA (F) Western blotting analyses following CHIRP assays in AGS cells confirmed the interaction between LINC00893 and RBFOX2 (G) RIP assay followed by RT-qPCR suggested LINC00893 binds to RBFOX2 (H) Representative images of immunofluorescence staining for RBFOX2 expression in MKN28 cells after transfection with LINC00893 overexpression or empty vectors. Scale bar: 20 μm (I) Western blotting to analyze total level of RBFOX2 after LINC00893 were overexpressed or knocked down in GC cells (J) A CHX treatment was administered to assess RBFOX2 degradation in AGS and MKN28 cells (K) MG-132 abolished the downregulation of RBFOX2 protein expression induced by LINC00893 overexpression in AGS and MKN28 cells (L) Ubiquitination assays revealed that LINC00893 overexpression increased the level of the ubiquitinated RBFOX2 protein in AGS and MKN28 cells.
Then, we examined RBFOX2 expression when LINC00893 levels were altered to determine whether LINC00893 regulated RBFOX2 expression. Interestingly, LINC00893 overexpression reduced the levels of the RBFOX2 protein (Figures 5H,I). According to our findings above, we inferred that LINC00893 regulated RBFOX2 protein stability. To elucidate the underlying mechanism by which LINC00893 regulated the RBFOX2 protein, we performed cycloheximide chase assays and observed an obvious shorter half-life of RBFOX2 in LINC00893-overexpressing cells than in empty control cells (Figure 5J). The ubiquitin-proteasome degradation pathway has been reported to be a primary mechanism regulating protein levels in cells. The proteasome inhibitor MG132 antagonized the decrease in RBFOX2 protein levels caused by LINC00893, indicating that LINC00893 downregulated the expression of RBFOX2 by promoting its proteasomal degradation (Figure 5K). Subsequent ubiquitination assays revealed that LINC00893 overexpression increased the level of the ubiquitinated RBFOX2 protein in AGS and MKN28 cells (Figure 5L). Taken together, LINC00893 directly interacts with RBFOX2 and promotes its ubiquitin-mediated degradation in GC cells.
RBFOX2 Is Overexpressed in GC Cells and Promotes Its Progression
RBFOX2 promotes tumor progression by acting as a splicer of EMT-related genes in several types of cancers (Shapiro et al., 2011; Braeutigam et al., 2014); however, the status and role of RBFOX2 in GC is unclear. Therefore, we focused on how RBFOX2 affects EMT-related genes and its role in GC cells. First, we knocked down RBFOX2 with siRNA or increased it by transfecting an overexpression plasmid. We found that the proliferation ability of GC cells was greatly decreased when RBFOX2 was downregulated but increased when RBFOX2 was upregulated (Figures 6A,B). The migration and invasion abilities of GC cells showed the same tendency (Figures 6C,D). Additionally, we found that when RBFOX2 was knocked down, the expression of E-cadherin increased, while N-cadherin and vimentin decreased, in accordance with the results when RBFOX2 was overexpressed (Figure 6E). In addition, analysis from KM Plotter (http://kmplot.com/analysis/) showed that high expression of RBFOX2 was correlated with a worse survival of GC patients, which is the opposite tendency with LINC00893 in terms of the overall survival rate (Figure 6F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g006.jpg)
FIGURE 6. RBFOX2 promotes malignant behavior and EMT in GC cells in vitro (A,B) The effects of RBFOX2 knockdown or overexpression on proliferation and plate colony-forming ability were measured in GC cells (C,D) The effects of RBFOX2 knockdown or overexpression on migration and invasion were detected (E) Relative expression of E-cadherin, N-cadherin and Vimentin was detected by western blotting in gastric cancer cells after RBFOX2 were either knocked down or overexpressed (F) Survival rates analysis of between RBFOX2 and GC were analyzed by KM Plotter database (https://kmplot.com). *p < 0.05, **p < 0.01, ***p < 0.001.
The Effects of LINC00893 on GC Cells Are Mediated by RBFOX2
Subsequently, we determined whether LINC00893 exerts its function though RBFOX2. We found that the LINC00893-mediated suppression of proliferation, migration and invasion was remarkably abolished when we overexpressed LINC00893 and RBFOX2 in GC cells at the same time. Consistent with these results, transfecting siRNA of RBFOX2 and ASO of LINC00893 simultaneously strikingly eliminated the promotion of GC cell proliferation, migration and invasion induced by LINC00893 knockdown, suggesting that LINC00893-mediated RBFOX2 downregulation contributes to the suppression of LINC00893 on GC cell proliferation, migration and invasion (Figures 7A–D). In addition, the changes of the EMT-related proteins N-cadherin, vimentin and E-cadherin caused by the upregulation or downregulation of LINC00893 was rescued by the overexpression or knockdown of RBFOX2, respectively (Figure7E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/796451/fcell-09-796451-HTML/image_m/fcell-09-796451-g007.jpg)
FIGURE 7. LINC00893 suppresses the proliferation and metastasis of GC cells through RBFOX2 (A,B) The effects of LINC00893 overexpression or knockdown on proliferation and plate colony-forming ability were rescued by RBFOX2 in GC cells (C,D) The effects of LINC00893 overexpression or knockdown on migration and invasion ability were rescued by RBFOX2 in GC cells (E,F) The effects of LINC00893 overexpression or knockdown on E-cadherin, N-cadherin and Vimentin were rescued by RBFOX2 in GC cells.
Discussion
LncRNAs are classified as a large sort of transcripts with unknown biological roles and functions. Emerging researches have implied that dysregulation of certain lncRNAs could play a vital role in the initiation, tumorigenesis, and progression of GC (Tan et al., 2020; Xie et al., 2020; Yuan et al., 2020). LncRNA can crosstalk with chromatin, DNA, RNA, and proteins, regulating GC cell progression via transcriptional and posttranscriptional mechanisms (Kong et al., 2018; Zhuo et al., 2019). The cellular distribution of lncRNAs is informative for lncRNA functions (Shin et al., 2020). A large number of lncRNAs function by binding to a protein partner (Castello et al., 2016; Chen, 2016). Cytoplasmic lncRNAs mainly play roles in modulating mRNA stability or translation and thus influencing cellular signaling cascades, while nuclear lncRNAs are enriched for functionality involving chromatin interactions, transcriptional regulation, and RNA processing (Huang et al., 2018). A growing body of evidence has shown that lncRNAs can interact with proteins to modulate protein function, regulate protein-protein interactions, or directly localize within cellular compartments (Batista and Chang, 2013).
In this study, we analyzed tumor suppressor lncRNAs induced by p53 through RNA sequencing of GC tissues and p53-elevated AGS cells, and 40 lncRNAs were differentially expressed in both RNA sequence data. Some of these lncRNAs have been reported to play vital roles in GC and other cancers. LINC01140 promotes the progression and tumor immune escape in lung cancer (Xia et al., 2021). LncRNA GATA6-AS1 Inhibits LNM and EMT via FZD4 through the Wnt/β-Catenin Signaling Pathway in GC (Li et al., 2020). We screened out and explored an interesting lncRNA, LINC00893, which was downregulated in GC tissues but increased after the elevation of p53 in AGS cells. Moreover, LINC00893 regulated the proliferation, migration and invasion of GC cells in vitro and in vivo. Mechanistically, our results indicated that LINC00893 binds to RBFOX2 and regulates the subsequent ubiquitin-mediated degradation. Finally, we proved that LINC00893 suppresses the progression of GC cells through RBFOX2.
RBFOX2 is a regulator of alternative splicing involved in the EMT process and contributes to the increased invasiveness of cancer cells that have undergone EMT (Yeo et al., 2009). RBFOX2 has been linked to a proepithelial to EMT phenotype in several cancer types (Murphy et al., 2018; Choi et al., 2019). In the current study, we revealed the expression status and oncogenic role of RBFOX2 in GC cells in vitro. We also found that its expression can be regulated by LINC00893 through direct binding. KEGG and GO analyses showed that LINC00893 was related to ubiquitin-mediated proteolysis and ubiquitin-protein transferase activity. We verified that the RNA-protein binding of LINC00893 and RBFOX2 can influence the ubiquitination of RBFOX2 to reduce its expression, which is a mechanism often reported in the interactions of lncRNAs and proteins.
TP53 is the most important tumor suppressor and the most frequently somatically mutated gene in human cancer; therefore, it is essential to figure out the functions of lncRNA components in the p53 network. Despite decades of studies on its significance in tumor suppression, the downstream effectors of the p53 response remain unclear. In recent years, increasing evidence has revealed that both proteins and noncoding genes including lncRNAs are involved in p53 regulatory networks. LncRNAs regulated by p53 in indirect ways play critical roles in genomic stability, DNA damage repair, cell cycle arrest, and apoptosis, which contributes to the tumor-suppressor functions of p53 (Chaudhary and Lal, 2017). Here, we discovered a lncRNA, LINC00893, regulated by p53 indirectly in GC cells, which enriched the p53 regulatory network, yet the specific mechanism by which p53 regulates LINC00893 remains to be verified. As we mentioned above, some other lncRNAs are also regulated by p53 in different ways. We infer that there may be connections between LINC00893 and them in the p53 regulatory network. It is possible that there is a synergic effect between LINC00893 and those regulated by p53 positively, such as lincRNA-p21, GUARDIN and TP53TG1 (Bao et al., 2015; Diaz-Lagares et al., 2016; Hu et al., 2018). On the other side, an antagonic effect may be possible between LINC00893 and those negatively regulated by the WT p53 pathway, such as LINP1(Zhang et al., 2016). However, further investigations are needed to verify these hypotheses.
In conclusion, we demonstrate that LINC00893 promotes the degradation of RBFOX2 by direct binding, thereby suppressing GC progression by decreasing the EMT and functional ability of GC cells. LINC00893 could serve as a prognostic indicator and potential therapeutic target for GC patients.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Sequence Read Archive (SRA) of RNA-seq data, https://www.ncbi.nlm.nih.gov/sra. Accession numbers: PRJNA792194 and PRJNA792206.
Ethics Statement
The studies involving human participants were reviewed and approved by The Ethical Committee Review Board of the First Affiliated Hospital of Sun Yat-sen University. The patients/participants provided their written informed consent to participate in this study.The animal study was reviewed and approved by The Ethical Committee Review Board of the First Affiliated Hospital of Sun Yat-sen University (No. 2020035).
Author Contributions
JX designed the research; XO, XZ, JL, HL, DF and QC performed the research; XO and XZ analyzed the data; XO wrote the manuscript; JX, JY, SC and YH reviewed the manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (81871915, 81672343 and 81472260), Natural Science Foundation of Guangdong, China (2017A030313570 and 2015A030313053). Science and Technology Program of Guangzhou, China (201607010050). The Fundamental Research Funds for the Central Universities (15ykpy17).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.796451/full#supplementary-material
References
Allen, M. A., Andrysik, Z., Dengler, V. L., Mellert, H. S., Guarnieri, A., Freeman, J. A., et al. (2014). Global Analysis of P53-Regulated Transcription Identifies its Direct Targets and Unexpected Regulatory Mechanisms. Elife 3, e02200. doi:10.7554/eLife.02200
Bao, X., Wu, H., Zhu, X., Guo, X., Hutchins, A. P., Luo, Z., et al. (2015). The P53-Induced lincRNA-P21 Derails Somatic Cell Reprogramming by Sustaining H3K9me3 and CpG Methylation at Pluripotency Gene Promoters. Cell Res 25 (1), 80–92. doi:10.1038/cr.2014.165
Batista, P. J., and Chang, H. Y. (2013). Long Noncoding RNAs: Cellular Address Codes in Development and Disease. Cell 152 (6), 1298–1307. doi:10.1016/j.cell.2013.02.012
Braeutigam, C., Rago, L., Rolke, A., Waldmeier, L., Christofori, G., and Winter, J. (2014). The RNA-Binding Protein Rbfox2: an Essential Regulator of EMT-Driven Alternative Splicing and a Mediator of Cellular Invasion. Oncogene 33 (9), 1082–1092. doi:10.1038/onc.2013.50
Cai, Q., Shi, P., Yuan, Y., Peng, J., Ou, X., Zhou, W., et al. (2021). Inflammation-Associated Senescence Promotes Helicobacter Pylori-Induced Atrophic Gastritis. Cell Mol. Gastroenterol. Hepatol. 11 (3), 857–880. doi:10.1016/j.jcmgh.2020.10.015
Castello, A., Fischer, B., Frese, C. K., Horos, R., Alleaume, A.-M., Foehr, S., et al. (2016). Comprehensive Identification of RNA-Binding Domains in Human Cells. Mol. Cel 63 (4), 696–710. doi:10.1016/j.molcel.2016.06.029
Chaudhary, R., and Lal, A. (2017). Long Noncoding RNAs in the P53 Network. WIREs RNA 8 (3), e1410. doi:10.1002/wrna.1410
Chen, L.-L. (2016). Linking Long Noncoding RNA Localization and Function. Trends Biochem. Sci. 41 (9), 761–772. doi:10.1016/j.tibs.2016.07.003
Choi, S., Sa, M., Cho, N., Kim, K. K., and Park, S.-H. (2019). Rbfox2 Dissociation from Stress Granules Suppresses Cancer Progression. Exp. Mol. Med. 51 (4), 1–12. doi:10.1038/s12276-019-0246-y
Diaz-Lagares, A., Crujeiras, A. B., Lopez-Serra, P., Soler, M., Setien, F., Goyal, A., et al. (2016). Epigenetic Inactivation of the P53-Induced Long Noncoding RNA TP53 Target 1 in Human Cancer. Proc. Natl. Acad. Sci. USA 113 (47), E7535–E7544. doi:10.1073/pnas.1608585113
Djebali, S., Davis, C. A., Merkel, A., Dobin, A., Lassmann, T., Mortazavi, A., et al. (2012). Landscape of Transcription in Human Cells. Nature 489 (7414), 101–108. doi:10.1038/nature11233
Donehower, L. A., Soussi, T., Korkut, A., Liu, Y., Schultz, A., Cardenas, M., et al. (2019). Integrated Analysis of TP53 Gene and Pathway Alterations in the Cancer Genome Atlas. Cel Rep. 28 (5), 1370–1384. doi:10.1016/j.celrep.2019.07.001
Dragomir, M. P., Kopetz, S., Ajani, J. A., and Calin, G. A. (2020). Non-coding RNAs in GI Cancers: from Cancer Hallmarks to Clinical Utility. Gut 69, 748–763. doi:10.1136/gutjnl-2019-318279
Duan, Y., Zhang, W., Cheng, Y., Shi, M., and Xia, X.-Q. (2021). A Systematic Evaluation of Bioinformatics Tools for Identification of Long Noncoding RNAs. RNA 27 (1), 80–98. doi:10.1261/rna.074724.120
Endo, S., Yamato, K., Hirai, S., Moriwaki, T., Fukuda, K., Suzuki, H., et al. (2011). Potent In Vitro and In Vivo Antitumor Effects of MDM2 Inhibitor Nutlin-3 in Gastric Cancer Cells. Cancer Sci. 102 (3), 605–613. doi:10.1111/j.1349-7006.2010.01821.x
Farre, D., Roset, R., Huerta, M., Adsuara, J. E., Rosello, L., Alba, M. M., et al. (2003). Identification of Patterns in Biological Sequences at the ALGGEN Server: PROMO and MALGEN. Nucleic Acids Res. 31 (13), 3651–3653. doi:10.1093/nar/gkg605
Hu, W. L., Jin, L., Xu, A., Wang, Y. F., Thorne, R. F., Zhang, X. D., et al. (2018). GUARDIN Is a P53-Responsive Long Non-coding RNA that Is Essential for Genomic Stability. Nat. Cel Biol 20 (4), 492–502. doi:10.1038/s41556-018-0066-7
Huang, R., Han, M., Meng, L., and Chen, X. (2018). Transcriptome-wide Discovery of Coding and Noncoding RNA-Binding Proteins. Proc. Natl. Acad. Sci. USA 115 (17), E3879–e3887. doi:10.1073/pnas.1718406115
Huarte, M., Guttman, M., Feldser, D., Garber, M., Koziol, M. J., Kenzelmann-Broz, D., et al. (2010). A Large Intergenic Noncoding RNA Induced by P53 Mediates Global Gene Repression in the P53 Response. Cell 142 (3), 409–419. doi:10.1016/j.cell.2010.06.040
Innocente, S. A., and Lee, J. M. (2005). p53 Is a NF-Y- and P21-independent, Sp1-dependent Repressor of Cyclin B1 Transcription. FEBS Lett. 579 (5), 1001–1007. doi:10.1016/j.febslet.2004.12.073
Jain, A. K. (2020). Emerging Roles of Long Non-coding RNAs in the P53 Network. RNA Biol. 17 (11), 1648–1656. doi:10.1080/15476286.2020.1770981
Jangi, M., Boutz, P. L., Paul, P., and Sharp, P. A. (2014). Rbfox2 Controls Autoregulation in RNA-Binding Protein Networks. Genes Dev. 28 (6), 637–651. doi:10.1101/gad.235770.113
Katayama, S., Tomaru, Y., Kasukawa, T., Waki, K., Nakanishi, M., Nakamura, M., et al. (2005). Antisense Transcription in the Mammalian Transcriptome. Science 309 (5740), 1564–1566. doi:10.1126/science.1112009
Kong, F., Deng, X., Kong, X., Du, Y., Li, L., Zhu, H., et al. (2018). ZFPM2-AS1, a Novel lncRNA, Attenuates the P53 Pathway and Promotes Gastric Carcinogenesis by Stabilizing MIF. Oncogene 37 (45), 5982–5996. doi:10.1038/s41388-018-0387-9
Kung, J. T. Y., Colognori, D., and Lee, J. T. (2013). Long Noncoding RNAs: Past, Present, and Future. Genetics 193 (3), 651–669. doi:10.1534/genetics.112.146704
Li, D., She, J., Hu, X., Zhang, M., Sun, R., and Qin, S. (2021). The ELF3-Regulated lncRNA UBE2CP3 Is Over-stabilized by RNA-RNA Interactions and Drives Gastric Cancer Metastasis via miR-138-5p/ITGA2 axis. Oncogene 40 (35), 5403–5415. doi:10.1038/s41388-021-01948-6
Li, J.-H., Liu, S., Zhou, H., Qu, L.-H., and Yang, J.-H. (2014). starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and Protein-RNA Interaction Networks from Large-Scale CLIP-Seq Data. Nucl. Acids Res. 42, D92–D97. (Database issue). doi:10.1093/nar/gkt1248
Li, H., Zhang, Y., Strose, A., Tedesco, D., Gurova, K., and Selivanova, G. (2014a). Integrated High-Throughput Analysis Identifies Sp1 as a Crucial Determinant of p53-Mediated Apoptosis. Cell Death Differ. 21 (9), 1493–1502. doi:10.1038/cdd.2014.69
Li, Y., Wu, Z., Yuan, J., Sun, L., Lin, L., Huang, N., et al. (2017). Long Non-coding RNA MALAT1 Promotes Gastric Cancer Tumorigenicity and Metastasis by Regulating Vasculogenic Mimicry and Angiogenesis. Cancer Lett. 395, 31–44. doi:10.1016/j.canlet.2017.02.035
Li, Z.-T., Zhang, X., Wang, D.-W., Xu, J., Kou, K.-J., Wang, Z.-W., et al. (2020). Overexpressed lncRNA GATA6-AS1 Inhibits LNM and EMT via FZD4 through the Wnt/β-Catenin Signaling Pathway in GC. Mol. Ther. - Nucleic Acids 19, 827–840. doi:10.1016/j.omtn.2019.09.034
Lin, M. F., Jungreis, I., and Kellis, M. (2011). PhyloCSF: a Comparative Genomics Method to Distinguish Protein Coding and Non-coding Regions. Bioinformatics 27 (13), i275–i282. doi:10.1093/bioinformatics/btr209
Liu, X.-h., Sun, M., Nie, F.-q., Ge, Y.-b., Zhang, E.-b., Yin, D.-d., et al. (2014). Lnc RNA HOTAIR Functions as a Competing Endogenous RNA to Regulate HER2 Expression by Sponging miR-331-3p in Gastric Cancer. Mol. Cancer 13, 92. doi:10.1186/1476-4598-13-92
Mehta, R., Kommalapati, A., and Kim, R. D. (2020). The Impact of Ramucirumab Treatment on Survival and Quality of Life in Patients with Gastric Cancer. Cancer Manag. Res. 12, 51–57. doi:10.2147/cmar.S199827
Mello, S. S., Sinow, C., Raj, N., Mazur, P. K., Bieging-Rolett, K., Broz, D. K., et al. (2017). Neat1 Is a P53-Inducible lincRNA Essential for Transformation Suppression. Genes Dev. 31 (11), 1095–1108. doi:10.1101/gad.284661.116
Murphy, P. A., Butty, V. L., Boutz, P. L., Begum, S., Kimble, A. L., Sharp, P. A., et al. (2018). Alternative RNA Splicing in the Endothelium Mediated in Part by Rbfox2 Regulates the Arterial Response to Low Flow. Elife 7, e29494. doi:10.7554/eLife.29494
Pan, T., Yu, Z., Jin, Z., Wu, X., Wu, A., Hou, J., et al. (2021). Tumor Suppressor Lnc-CTSLP4 Inhibits EMT and Metastasis of Gastric Cancer by Attenuating HNRNPAB-dependent Snail Transcription. Mol. Ther. - Nucleic Acids 23, 1288–1303. doi:10.1016/j.omtn.2021.02.003
Sánchez, Y., Segura, V., Marín-Béjar, O., Athie, A., Marchese, F. P., González, J., et al. (2014). Genome-wide Analysis of the Human P53 Transcriptional Network Unveils a lncRNA Tumour Suppressor Signature. Nat. Commun. 5, 5812. doi:10.1038/ncomms6812
Shapiro, I. M., Cheng, A. W., Flytzanis, N. C., Balsamo, M., Condeelis, J. S., Oktay, M. H., et al. (2011). An EMT-Driven Alternative Splicing Program Occurs in Human Breast Cancer and Modulates Cellular Phenotype. Plos Genet. 7 (8), e1002218. doi:10.1371/journal.pgen.1002218
Sharma, V., Khurana, S., Kubben, N., Abdelmohsen, K., Oberdoerffer, P., Gorospe, M., et al. (2015). A BRCA 1-interacting Lnc RNA Regulates Homologous Recombination. EMBO Rep. 16 (11), 1520–1534. doi:10.15252/embr.201540437
Shin, T.-J., Lee, K.-H., and Cho, J.-Y. (2020). Epigenetic Mechanisms of LncRNAs Binding to Protein in Carcinogenesis. Cancers 12 (10), 2925. doi:10.3390/cancers12102925
Sung, H., Ferlay, J., Siegel, R. L., Laversanne, M., Soerjomataram, I., Jemal, A., et al. (2021). Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A. Cancer J. Clin. 71, 209–249. doi:10.3322/caac.21660
Tan, H., Zhang, S., Zhang, J., Zhu, L., Chen, Y., Yang, H., et al. (2020). Long Non-coding RNAs in Gastric Cancer: New Emerging Biological Functions and Therapeutic Implications. Theranostics 10 (19), 8880–8902. doi:10.7150/thno.47548
Tang, Z., Li, C., Kang, B., Gao, G., Li, C., and Zhang, Z. (2017). GEPIA: a Web Server for Cancer and normal Gene Expression Profiling and Interactive Analyses. Nucleic Acids Res. 45 (W1), W98–w102. doi:10.1093/nar/gkx247
Vaseva, A. V., Yallowitz, A. R., Marchenko, N. D., Xu, S., and Moll, U. M. (2011). Blockade of Hsp90 by 17AAG Antagonizes MDMX and Synergizes with Nutlin to Induce P53-Mediated Apoptosis in Solid Tumors. Cell Death Dis 2, e156. doi:10.1038/cddis.2011.39
Venables, J. P., Brosseau, J.-P., Gadea, G., Klinck, R., Prinos, P., Beaulieu, J.-F., et al. (2013). RBFOX2 Is an Important Regulator of Mesenchymal Tissue-specific Splicing in Both normal and Cancer Tissues. Mol. Cell Biol. 33 (2), 396–405. doi:10.1128/MCB.01174-12
Vousden, K. H., and Prives, C. (2009). Blinded by the Light: The Growing Complexity of P53. Cell 137 (3), 413–431. doi:10.1016/j.cell.2009.04.037
Wang, C.-J., Zhu, C.-C., Xu, J., Wang, M., Zhao, W.-Y., Liu, Q., et al. (2019). The lncRNA UCA1 Promotes Proliferation, Migration, Immune Escape and Inhibits Apoptosis in Gastric Cancer by Sponging Anti-tumor miRNAs. Mol. Cancer 18 (1), 115. doi:10.1186/s12943-019-1032-0
Wei, C., Xiao, R., Chen, L., Cui, H., Zhou, Y., Xue, Y., et al. (2016). RBFox2 Binds Nascent RNA to Globally Regulate Polycomb Complex 2 Targeting in Mammalian Genomes. Mol. Cel 62 (6), 982. doi:10.1016/j.molcel.2016.06.003
Wei, J., Nagy, T. A., Vilgelm, A., Zaika, E., Ogden, S. R., Romero–Gallo, J., et al. (2010). Regulation of P53 Tumor Suppressor by Helicobacter pylori in Gastric Epithelial Cells. Gastroenterology 139 (4), 1333–1343. doi:10.1053/j.gastro.2010.06.018
Xia, R., Geng, G., Yu, X., Xu, Z., Guo, J., Liu, H., et al. (2021). LINC01140 Promotes the Progression and Tumor Immune Escape in Lung Cancer by Sponging Multiple microRNAs. J. Immunother. Cancer 9 (8), e002746. doi:10.1136/jitc-2021-002746
Xie, S., Chang, Y., Jin, H., Yang, F., Xu, Y., Yan, X., et al. (2020). Non-coding RNAs in Gastric Cancer. Cancer Lett. 493, 55–70. doi:10.1016/j.canlet.2020.06.022
Xu, J., Lu, Y., Liu, Q., Xia, A., Zhao, J., Xu, X., et al. (2020). Long Noncoding RNA GMAN Promotes Hepatocellular Carcinoma Progression by Interacting with eIF4B. Cancer Lett. 473, 1–12. doi:10.1016/j.canlet.2019.12.032
Yeo, G. W., Coufal, N. G., Liang, T. Y., Peng, G. E., Fu, X.-D., and Gage, F. H. (2009). An RNA Code for the FOX2 Splicing Regulator Revealed by Mapping RNA-Protein Interactions in Stem Cells. Nat. Struct. Mol. Biol. 16 (2), 130–137. doi:10.1038/nsmb.1545
Yuan, L., Xu, Z.-Y., Ruan, S.-M., Mo, S., Qin, J.-J., and Cheng, X.-D. (2020). Long Non-coding RNAs towards Precision Medicine in Gastric Cancer: Early Diagnosis, Treatment, and Drug Resistance. Mol. Cancer 19 (1), 96. doi:10.1186/s12943-020-01219-0
Zhang, H., Ma, R.-R., Zhang, G., Dong, Y., Duan, M., Sun, Y., et al. (2022). Long Noncoding RNA Lnc-LEMGC Combines with DNA-PKcs to Suppress Gastric Cancer Metastasis. Cancer Lett. 524, 82–90. doi:10.1016/j.canlet.2021.09.042
Zhang, Y., He, Q., Hu, Z., Feng, Y., Fan, L., Tang, Z., et al. (2016). Long Noncoding RNA LINP1 Regulates Repair of DNA Double-Strand Breaks in Triple-Negative Breast Cancer. Nat. Struct. Mol. Biol. 23 (6), 522–530. doi:10.1038/nsmb.3211
Zhuo, W., Liu, Y., Li, S., Guo, D., Sun, Q., Jin, J., et al. (2019). 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 156 (3), 676–691. e611. doi:10.1053/j.gastro.2018.10.054
Keywords: p53, LINC00893, gastric cancer, RBFOX2, epithelial-mesenchymal transition
Citation: Ou X, Zhou X, Li J, Ye J, Liu H, Fang D, Cai Q, Cai S, He Y and Xu J (2022) p53-Induced LINC00893 Regulates RBFOX2 Stability to Suppress Gastric Cancer Progression. Front. Cell Dev. Biol. 9:796451. doi: 10.3389/fcell.2021.796451
Received: 16 October 2021; Accepted: 29 December 2021;
Published: 19 January 2022.
Edited by:
Wen G. Jiang, Cardiff University, United KingdomReviewed by:
Jing Zhang, Shanghai Jiao Tong University, ChinaJose Roberto Sotelo-Silveira, Instituto de Investigaciones Biológicas Clemente Estable, Uruguay
Copyright © 2022 Ou, Zhou, Li, Ye, Liu, Fang, Cai, Cai, He and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yulong He, heyulong@mail.sysu.edu.cn; Jianbo Xu, xjianb@mail.sysu.edu.cn
†These authors have contributed equally to this work and share first authorship