- Open Access
Linc00659, a long noncoding RNA, acts as novel oncogene in regulating cancer cell growth in colorectal cancer
Molecular Cancer volume 17, Article number: 72 (2018)
Colorectal cancer (CRC) is one of the most common cancers and causes of cancer-related death worldwide. In patients with CRC, metastasis is a crucial problem that leads to treatment failure and is the primary cause of the lethality of colon cancer. Long noncoding RNAs (lncRNAs) have recently emerged as critical molecules in the development, cell growth, apoptosis, and metastasis of CRC.
We investigated the transcriptome profiles of human lncRNAs in the primary tumor tissues and in the corresponding normal mucosa of two patients with CRC by using a microarray approach. The expression levels of lncRNAs were verified in colon cancer by real-time PCR. Using bioinformatics approach to illustrate putative biological function of Linc00659 in colon cancer. The effects of Linc00659 on cell growth, proliferation, cell cycle and apoptosis were studies by in vitro assays.
Our data revealed that compared with adjacent normal tissues, 201 lncRNAs were deregulated (fold change ≥ 4 or ≤ 0.25) in CRC tissues. Among them, the expression levels of Linc00659 were significantly increased in colon cancer, and high expression levels were correlated with poor survival in patients with CRC. Bioinformatics analysis results indicated that Linc00659 was significantly coexpressed with cycle-related genes in CRC. Linc00659 expression knockdown could significantly suppress colon cancer cell growth by impairing cell cycle progression. In addition, our results showed that Linc00659 expression knockdown could accelerate cell apoptosis in colon cancer cells treated with chemotherapy drugs. Meanwhile, our results also demonstrated that silencing of Linc00659 expression leads to cell growth inhibition and induced apoptosis, possibly by suppressing PI3K-AKT signaling in colon cancer.
Linc00659 is a novel oncogenic lncRNA involved in colon cancer cell growth by modulating the cell cycle. Our findings give an insight into lncRNA regulation and provide an application for colon cancer therapy.
Colorectal cancer (CRC) is one of the most common cancers and causes of cancer-related death worldwide [1, 2]. The major causes of the lethality of colon cancer are distant metastasis and drug resistance, which usually lead to a poor survival rate . In patients diagnosed as having advanced-stage CRC, distant metastasis is frequently observed, leading to poor therapeutic effects of chemotherapy, radiotherapy, and surgery . Therefore, developing a diagnostic biomarker for early detection of CRC or as a therapeutic target will improve the survival rate of patients with CRC.
High-throughput genome-scale approaches have revealed that a large amount of transcripts are produced from more than 80% of human genomic DNA, and only 1–2% of the human genome is composed of protein-coding genes . Therefore, the human genome comprehensively harbors numerous transcripts with no protein-coding function. Recent studies have shown that numerous noncoding RNAs (ncRNAs) play a critical role in regulating various cellular processes . Long noncoding RNAs (lncRNAs) are a group of long RNA transcripts that are more than 200 nucleotides in length; however, they lack protein translation ability. An increasing body of evidence indicates that lncRNAs play a critical role in regulating various cellular processes, including cell development and growth, the cell cycle, and cancer metastasis, by (1) regulating protein-coding gene expression, (2) modifying epigenetic regulation, (3) modulating the alternative splicing process, and (4) acting as decoys that titrate microRNAs [7, 8]. To date, some studies have reported that numerous lncRNAs, including HOTAIR, PCAT-1, PVT-1, H19, BANCR, 91H, CCAT1, and MALAT1, were dysregulated in colon cancer [9,10,11,12,13,14,15,16,17,18,19,20,21]. These dysfunctional lncRNAs have been reported to have tumor-suppressive or oncogenic roles in regulating cell growth, the cell cycle, and cell migration in CRC. However, the function of most lncRNAs remains largely unknown in CRC. In this study, our objective was to identify dysregulated lncRNAs in CRC. For this purpose, we compared the transcriptome profiles of primary tumor with adjacent normal tissue from two patients by using a microarray approach. We further evaluated the expression levels and biological role of Linc00659 in CRC through bioinformatics and an experimental approach. Our finding revealed that Linc00659 is a novel oncogenic lncRNA that regulates CRC cell growth by impairing cell cycle progression, and thus may be a potential therapeutic target for gene therapy.
Eight human CRC cell lines, namely DLD-1, colo205, colo320DM, LS174T, HCT116, SW620, LOVO, and HT29, were obtained from the American Type Culture Collection. Cells were cultured in high-glucose Dulbecco’s modified Eagle’s medium (Invitrogen, Grand Island, NY, USA) with 10% fetal bovine serum (Hyclone Laboratories, Inc., South Logan, UT, USA) and 1% penicillin/streptomycin (Gibco, Thermo Fisher Scientific Inc., Waltham, MA, USA) in plastic tissue culture plates in a humidified atmosphere containing 5% CO2 at 37 °C.
Colon cancer tissue and the corresponding adjacent normal mucosa samples were obtained from 91 patients with CRC who underwent surgery at the Department of Surgery, Taipei Veterans General Hospital, Taiwan. Informed consent was obtained from all patients.
Extraction of RNA
Total RNA from tissue was extracted using TRIzol reagent (Invitrogen, Grand Island, NY, USA), according to instructions in the user’s manual. In brief, tissue samples were homogenized in 1 mL TRIzol reagent and mixed with 0.2 mL chloroform to extract protein, before RNA was precipitated using 0.5 mL isopropanol. The concentration, purity, and amount of total RNA were determined using a Nanodrop 1000 spectrophotometer (Nanodrop Technologies Inc., Wilmington, DE, USA).
LncRNA profiling and pathway enrichment analysis
Total RNAs of four clinical samples were extracted, involving two primary CRC tissues and two corresponding adjacent normal tissues obtained from two patients with CRC, respectively. The microarray experiments and data analysis were performed by Welgene Biotech (Taipei, Taiwan) using the Agilent Oligo Chip (Agilent SurePrint G3 Human V2 GE; including protein-coding genes and lncRNAs). Differential expression protein-coding genes (tumor versus normal ≥ 2 and ≤ 0.5) were selected from microarray data and candidate genes were mapped onto the Kyoto Encyclopedia of Genes and Genomes pathways based on the Enzyme Commission numbers by using the R package SubpathwayMiner v.3.1. Subsequently, a hypergeometric test was performed to identify significantly enriched pathways and calculate the false-positive discovery rate in the false discovery rate-corrected q value.
Expression data from the cancer genome atlas
The transcriptome expression profiles of colon cancer were downloaded from The Cancer Genome Atlas (TCGA) data portal (https://cancergenome.nih.gov). The expression profiles of 616 colon cancer tissues and 51 adjacent normal tissues were obtained from TCGA data portal. In this study, the transcriptome profiles of 29 N-T pairs were used for coexpression analysis and 616 cases were included in the survival analysis.
Reverse transcription and real-time polymerase chain reaction
In this reaction, 2 μg of total RNA was reverse transcribed with random primers (Invitrogen; Thermo Fisher Scientific Inc., Waltham, MA, USA) and SuperScript IV Reverse Transcriptase according to the user’s manual (Invitrogen; Thermo Fisher Scientific Inc., Waltham, MA, USA). The reaction was performed with incubation at 42 °C for 1 h, and the enzyme was subsequently inactivated by incubation at 85 °C for 5 min. cDNA was used for real-time polymerase chain reaction (PCR) analysis with gene-specific primers, and gene expression was detected using a Fast SYBR Green Master Mix (Applied Biosystems; Thermo Fisher Scientific Inc., Waltham, MA, USA). The expression of lncRNA was normalized to that of glyceraldehyde 3-phosphate dehydrogenase (GAPDH; △Ct = target lncRNA Ct – GAPDH Ct). The individual primers used in this study were as follows:
Small interfering RNA transfection
Small interfering RNA (siRNA) oligonucleotides targeting Linc00659 (si-Linc00659, sense: 3′-TTGGGAGGAACACGAAGUCUU-5′ and antisense: 5′-UUCUGAAGCACAAGGAGGGTT-3′) and a scrambled oligo as a negative control were designed and synthesized by GenDiscovery Biotechnology (Taipei, Taiwan). Cells were transfected with a final concentration (10 mM) of individual siRNA or control using Lipofectamine RNAiMAX (Invitrogen; Thermo Fisher Scientific Inc., Waltham, MA, USA). After transfection for 24 h, the RNA was extracted and knockdown efficiency was evaluated by real-time PCR.
Stable Linc00659 knockdown with short hairpin RNA
Stable HCT116 cells with Linc00659 knockdown were generated by infecting colon cancer HCT116 cells with lentiviruses expressing sh-Linc00659 in the presence of 8 μg/mL polybrene for 24 h followed by puromycin (4 μg/mL) selection for 3–5 days. The shLuc vector targeting the luciferase gene provided puromycin resistance and was used as the control. Linc00659 expression was confirmed by a real-time PCR assay. In this study, we designed two shRNA sequence targeted Linc00659, and individual sequence of shRNA used for constructs in this study were as follows:
For cell proliferation analysis, 2500 living cells were transfected with either siLinc00659 or a scrambled control and plated onto 96-well plates. Cell growth was determined at 0, 1, 2, 3, and 4 days. Cell viability was determined using a CellTiter-Glo One Solution cell proliferation assay (Promega Corporation, Madison, WI, USA), which was performed according to the manufacturer’s instructions.
Soft agar assay
In brief, this assay was performed as follows: The base agar layer was prepared with 1% agar (Laboratorios CONDA, Torrejón de Ardoz, Madrid, Spain) dissolved in 1.5 mL phosphate-buffered saline on 6-well plates; subsequently, 1.5 mL of 0.7% agarose (Invitrogen, Grand Island, NY, USA) solution containing the cells (10,000 cells per well) was inoculated on top of the base agar layer. After allowing the solution to harden, 0.5 mL of fresh medium was added to the top of the hard agar layer. After 3–4 weeks, agar plates were stained with iodonitrotetrazolium chloride (Santa Cruz Biotechnology, Santa Cruz, CA, USA) and colonies were counted.
Colony formation assay
A total of 4000 cells were seeded on each well of a 6-well plate, and the cells were transfected with individual siRNA oligonucleotides by using Lipofectamine RNAiMAX (Invitrogen; Thermo Fisher Scientific Inc., Waltham, MA, USA) in the 6-well plate. After incubation at 37 °C for 3 days, the cultured medium was replaced with a new medium. The cells were allowed to incubate at 37 °C for 10 days. Cell culture plates containing colonies were fixed with 4% formaldehyde for 2 min and colonies were stained with crystal violet solution (containing 0.5% crystal violet, 5% formaldehyde, 50% ethanol, and 0.85% sodium chloride) for 2 h. Wells were rinsed with H2O after air drying. The crystal violet staining of cells from each well was solubilized using 1 mL per well of 10% acetic acid, and the absorbance (optical density) of the solution was measured on a spectrophotometer at a wavelength of 595 nm.
Cell cycle analysis
A total of 1 × 106 cells were collected and mixed with 70% ethanol in a fixative at − 20 °C overnight. Cells were then stained with 4′,6-diamidino-2-phenylindole (ChemoMetec, Gydevang, Lillerød, Denmark) and detected using NucleoView NC-3000 software (ChemoMetec, Gydevang, Lillerød, Denmark).
Apoptotic cells stained with Annexin V, propidium iodide, and Hoechst 33,342
For the image flow cytometry assay, CF488A Annexin V and PI Apoptosis kit (Biotium, Fremont, CA, USA) and Hoechst 33,342 (ChemoMetec, Gydevang, Lillerød, Denmark) were used to detect cells in early apoptotic and late apoptotic stages. After staining, the cells were analyzed in NucleoCounter NC-3000 (ChemoMetec, Gydevang, Lillerød, Denmark) and NucleoView NC-3000 for automated image flow cytometry analysis.
The cells were harvested 24 h after transient transfection, washed with phosphate-buffered saline, and subjected to lysis in radioimmunoprecipitation assay buffer (50 mM Tris–HCl pH 8.0, 150 mM NaCl, 1% NP40, 0.5% deoxycholate, 0.1% sodium dodecyl sulfate) at 4 °C for 30 min. Lysed cells were collected and centrifuged to remove cell debris. Protein assays were performed with the Bio-Rad Protein Assay kit based on the Bradford dye-binding procedure (Bio-Rad). Protein samples (40 μg) were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis in 10% or 12% resolving gel using a Mighty Small II Deluxe Mini Vertical Electrophoresis Unit (Hoefer, Inc., Holliston, MA, USA). Proteins were then electrotransferred to polyvinylidene difluoride membranes (PerkinElmer, Inc., Waltham, MA, USA) by Mighty Small Transfer Tank (Hoefer, Inc., Holliston, MA, USA). Subsequently, the membranes were blocked with a blocking buffer (50 mM Tris–HCl pH 7.6, 150 mM NaCl, 0.1% Tween 20, 5% nonfat dry milk, 0.05% sodium azide) for 1 h at room temperature and incubated with the following primary antibodies overnight at 4 °C: CCNA2 (1:1000; 18,202–1-AP, Proteintech Group, Inc., Rosemont, IL, USA), CCNB1(1:1000; 55,004–1-AP, Proteintech Group, Inc., Rosemont, IL, USA), CCND1 (1:1000; RM-9104-S, Thermo Fisher Scientific Inc., Waltham, MA, USA), CDK4 (1:1000; MS-299-P, Thermo Fisher Scientific Inc., Waltham, MA, USA), CDKN1B (1:1000; #3686, Cell Signaling Technology, Inc., Beverly, MA, USA), CDKN1A (1:1000; #2947, Cell Signaling Technology, Inc., Beverly, MA, USA), PARP (1:1000; sc-8007, Santa Cruz Biotechnology, Inc. Santa Cruz, CA, USA), CASP3 (1:1000; #9662, Cell Signaling Technology, Inc., Beverly, MA, USA), E2F1 (1:200, sc-251, Santa Cruz Biotechnology, Inc., Santa Cruz, CA, USA), PI3K p85 (1:1000; #4292, Cell Signaling Technology, Inc., MA, USA), AKT (1:1000; #4691, Cell Signaling Technology, Inc., MA, USA), Phospho-AKT (Thr308) (1:1000; #4056, Cell Signaling Technology, Inc., MA, USA), Phospho-AKT (Ser473) (1:1000; #4060, Cell Signaling Technology, Inc., MA, USA), GSK3β (1:1000; #9315, Cell Signaling Technology, Inc., MA, USA), Bcl-2 (1:1000; #4223, Cell Signaling Technology, Inc., MA, USA), BAX (1:1000; #5023, Cell Signaling Technology, Inc., MA, USA), Bad (1:1000; #9239, Cell Signaling Technology, Inc., MA, USA), Phospho-Bad(Ser136) (1:1000; #4366, Cell Signaling Technology, Inc., MA, USA), β-catenin (1:4000; C2206, Sigma-Aldrich, Inc., Beverly, MO, USA) and ACTB (1:2000, MAB1501, EMD Millipore, Billerica, MA, USA). The membranes were then incubated with anti-rabbit (sc-2004) or anti-mouse (sc-2005) immunoglobulin G horseradish peroxidase-conjugated secondary antibodies (1:10000, Santa Cruz Biotechnology, Inc., Santa Cruz, CA, USA) for 1 h at room temperature. After three washes with Tris-buffered saline with Tween-20 buffer (50 mM Tris–HCl pH 7.6, 150 mM NaCl, 0.1% Tween-20), immunoreactive bands were detected using a WesternBright ECL (Advansta, Menlo Park, CA, USA).
The expression levels of Linc00659 in colon cancer from TCGA database were analyzed using Student t tests. The correlation of Linc00659 with the protein-coding genes in colon cancer was determined through Pearson coefficient analysis, with r and p values as indicated. Cumulative survival curves were estimated using the Kaplan–Meier method, and comparison between survival curves was performed using the log-rank test. The difference was considered significant when p < 0.05. The data on the expression levels of lncRNAs in paired colon tissues were analyzed using a paired t test. Cell proliferation, colony formation, and soft agar assay experiments were performed in triplicate. Histograms present the mean values, and the error bars indicate the standard deviation. These data were analyzed using Student t tests.
Generation of lncRNA expression profiles in CRC
We determined the gene expression profiles of the primary tumor and the corresponding adjacent normal tissues in two patients with CRC by using a microarray approach, in which protein-coding genes and lncRNAs were included (Agilent SurePrint G3 Human V2 GE; Fig. 1a). Most genes were consistently expressed between the CRC and the corresponding adjacent normal tissues in the two patients (R2 > 0.86; Fig. 1b and c). In total, 5790 protein-coding genes were detected as differentially expressed (fold change ≤ 2 or ≥ 0.5) between CRC and the corresponding adjacent normal mucosa, comprising 4866 upregulated genes and 924 downregulated genes, respectively (Fig. 1d, upper panels). We further examined the putative biological function of these differentially expressed genes by using pathway enrichment analysis. Our data revealed that these differentially expressed genes were significantly enriched in cell cycle, DNA replication, RNA transport, and pyrimidine metabolism signaling pathways (Fig. 1d, bottom panel). Abnormal lncRNAs in CRC were also identified by analyzing the same microarray data. Consistent with protein data, most lncRNAs were consistently expressed between the CRC and the corresponding adjacent normal tissues in the two patients with CRC (R2 > 0.83; Fig. 2a and b). Only a small subset of lncRNAs was differentially expressed between CRC and the corresponding adjacent normal tissues. A total of 201 lncRNAs were either upregulated (n = 165; fold change ≥ 4) or downregulated (n = 36; fold change ≤ 0.25) in the two CRC cases.
LncRNA candidates were confirmed in colon cancer
The top 10 of these abnormal lncRNAs upregulated and downregulated are presented in Table 1. To validate the microarray analysis results, 4 of the top 10 differentially expressed lncRNAs were selected for validation in 89 matched CRC samples by using reverse transcription-quantitative PCR. Our results indicated that the expression levels of Linc00659 and MNX1-AS1 were significantly increased, whereas those of Loc339524 and Linc00657 were significantly decreased in colon cancer (Fig. 2d–i), which were consistent with our microarray results.
Linc00659 coexpressed with cycle-related genes
The biological and clinical effects of Linc00659 expression were further examined by analyzing TCGA database. A total of 667 transcriptome data items were downloaded from TCGA database, comprising 51 normal and 616 CRC tissues. As shown in Fig. 3a, the expression levels of Linc00659 significantly increased in colon cancer compared with normal mucosa. Furthermore, high expression levels of Linc00659 were significantly correlated with a short survival curve of patients with CRC (Fig. 3b), implying that Linc00659 might play an oncogenic role in colon cancer progression. We further explored the biological role of Linc00659 in colon cancer tumorigenesis by identifying coexpression gene networks. A total of 1218 protein-coding genes were coexpressed with Linc00659 in colon cancer (r > 0.7 or <− 0.7; Fig. 4a). By determining these signaling pathways of the protein-coding genes coexpressed with Linc00659, we found that these highly interconnected genes were enriched in several cell proliferation-related signaling pathways (Fig. 4b). Notably, we also found that Linc00659 expression had a highly positive correlation with cell cycle-related gene expression in colon cancer, implying that Linc00659 might have an oncogenic role in colon cancer cell growth by accelerating cell cycle progression (Fig. 4c).
Linc00659 was involved in colon cancer growth
We examined the expression levels of Linc00659 in eight colon cancer cell lines, namely colo320dm, colo205, DLD1, Lovo, SW620, HCT116, LS174T, and HT29, by using real-time PCR. Compared with adjacent normal mucosa, the expression levels of Linc00659 were clearly overexpressed in eight colon cancer cell lines by more than 20-fold (Fig. 5a). Based on this result, we selected a loss-of-function approach to study the biological function of Linc00659 in colon cancer cells. As demonstrated in Fig. 5b, Linc00659 expression significantly decreased by 70% in HCT116 cells after siRNA transfection.. Knockdown of Linc00659 expression significantly suppressed cancer cell proliferation and colony formation in HCT116 (Fig. 5c–e). Using a soft agar assay, we indicated that Linc00659 knockdown could significantly suppress anchorage-independent growth in HCT116 cells (Fig. 5f and g). We also observed that the cell proliferation of Lovo, LS174T, DLD-1 and SW620 could be suppressed after Linc00659 knockdown for 4 days (Additional file 1: Figure S1). We further established stable knockdown of Linc00659 expression through sh-Linc00659 transfection and puromycin selection for 2 weeks. As illustrated in Fig. 6a and Additional file 2: Figure S2A, the expression levels of Linc00659 were significantly decreased by 60% in sh-Linc00659 compared with the control. The stable knockdown of Linc00659 expression exhibited similar effects on cell growth in accordance with si-Linc00659 treatment. The cell proliferation, colony formation, and anchorage-independent growth abilities were significantly decreased in HCT116 cells with stable Linc00659 knockdown (Fig. 6b–f and Additional file 2: Figure S2B-F).
Knockdown of Linc00659 suppressed cell growth by impairing cell cycle progression and inducing cell apoptosis
Our results revealed that Linc00659 could regulate colon cancer growth. Coexpression data indicated that Linc00659 had a high correlation with cell cycle-related genes; therefore, we explored the mechanism underlying the regulation of cell growth by examining the cell cycle distribution and apoptosis after Linc00659 knockdown. Image flow assay results indicated that the cell cycle distribution was disrupted after Linc00659 knockdown. Transient knockdown of Linc00659 cells also indicated that the G1 phase decreased in HCT116 cells (Fig. 7a and b). Furthermore, an examination of the cell cycle-related protein expression levels revealed that the expression levels of cyclin A2, cyclin B1, cyclin D1, CDK4, and E2F1 decreased in HCT116 cells with Linc00659 knockdown (Fig. 7c). Notably, the sub-G1 population cells were obviously increased in Linc00659 cells subjected to transient knockdown, including HCT116 (Fig. 7a and b). Apoptosis assay results indicated that the apoptotic cells significantly increased and that PARP and caspase-3 cleavage was induced in HCT116 cells with si-linc00659 transfection (Fig. 7d–f). These data indicate that Linc00659 knockdown could impair cell cycle progression and induce cell apoptosis.
Linc00659 knockdown confers drug sensitivity to colon cancer cells
To explore whether Linc00659 expression is involved in the regulation of drug responsiveness, we assessed the effects of Linc00659 knockdown on cell apoptosis in HCT116 with drug treatment, including oxaliplatin, 5-fluorouracil (5-FU), and irinotecan treatment. Linc00659 was stably knocked down and the respective control cells were exposed to various concentration of oxaliplatin, 5-FU, or irinotecan for 48 h, and cell viability was examined. As presented in Fig. 8a, Linc00659 knockdown could obviously accelerate the reduction of cell viability in HCT116 following oxaliplatin or 5-FU, and slightly reduced cell viability with irinotecan treatment. An examination of the cell cycle progression revealed that Linc00659 knockdown resulted in a significant increase of the sub-G1 phase of the cell cycle distribution after oxaliplatin and 5-FU treatment (Fig. 8b and c). In addition, Linc00659 knockdown increased apoptotic cells after oxaliplatin treatment compared to control group (Fig. 8d and e). Overall, the expression levels of Linc00659 conferred drug sensitivity in colon cancer by inducing apoptosis.
Effects of Linc00659 knockdown on PI3K-AKT-GSK3β signaling
Studies have indicated that AKT activation plays a crucial role in the cell-cycle progression and drug resistance of colon cancer [22,23,24]. Therefore, we sought to determine whether the P13K-AKT signaling pathway is involved in suppressing growth in HCT116 with Linc00659 knockdown. As shown in Fig. 9a, the protein levels of PI3K and p-AKT were decreased in HCT116 after Linc00659 knockdown. In addition, studies have indicated that AKT signaling acts has a crucial role in colon cancer cell survival by modulating BCl2 and Bad [25, 26]. Therefore, we examined these apoptosis-related proteins, revealing that the levels of anti-apoptosis proteins BCl2 and p-Bad both decreased after Linc00659 knockdown (Fig. 9b); we also observed that GSK3β expression was increased in Linc00659 knockdown compared with the control group. We further examined the expression levels of B-catenin, revealing that total protein and nuclear-B-catenin were decreased in HCT116 with Linc00659 knockdown (Fig. 9c). Taken together, our results demonstrate a putative mechanism of Linc00659 knockdown that impairs cell cycles and induces apoptosis; this may be a result of suppressing the PI3K-AKT-GSK3B axis in colon cancer.
Increasing evidence reveals that numerous lncRNAs play a critical role in the pathogenesis and progression of CRC. In the present study, transcriptome profiles were analyzed to identify differentially expressed lncRNA in colon cancer by using a microarray approach, which revealed that the expression levels of 201 lncRNAs were significantly different between colon cancer and the corresponding adjacent normal tissues. As expected, our microarray profiles identified several lncRNAs, which were aberrantly expressed in human cancer, such as LOC100190940, LOC100506178, H19, PVT1, VPS9D1-AS1, Linc00520, Linc00675, GNG12-AS1 CDKN2B-AS1, and CCAT1 [19, 27,28,29,30,31,32,33,34,35,36].
Chen et al. systematically examined global lncRNA dysregulation by analyzing RNA-seq data of colon cancer . They identified that LOC100190940 was upregulated in primary tumor tissues compared with adjacent normal tissues, and that high LOC100506178 expression levels were correlated with poor survival in patients with colon cancer. In addition, they reported that Linc00659 was highly expressed in primary tumor tissues with a higher H3K4me3 signal in its promoter region, compared with normal tissues. This result implies that the expression levels of Linc00659 might be regulated by the epigenetic mechanism in CRC . A previous study indicated that c-Myc could promote CCAT1 transcription by activating its promoter activity, and that ectopic CCAT1 expression could accelerate growth and invasion in colon cancer . Furthermore, high CCAT1 expression levels were significantly associated with poor prognosis and survival in patients with colon cancer [19, 29]. In this study, we identified hundreds of novel abnormal lncRNA genes to provide valuable information for lncRNA research in colon cancer.
The mechanisms through which lncRNAs participate in carcinogenesis and CRC progression are various and complex. To date, the mechanism underlying lncRNA function remains largely unknown, and the molecules with which they interact or how their expression influences signaling is unclear. Although studies have attempted to predict the putative mechanism of dysfunctional lncRNAs by using the bioinformatics approach, the lack of conserved sequences and the absence of functional motifs have made it difficult to determine the function of lncRNAs by analyzing their sequence or secondary structures. Using transcriptome data to determine which types of protein-coding genes are coexpressed with lncRNAs may facilitate exploring lncRNA function . In this study, we used TCGA database to determine coexpressed protein-coding genes and demonstrate the biological function of Linc00659 involved in modulating cancer cell growth. Using an experimental approach, we further confirmed that Linc00659 knockdown could significantly suppress colon cancer growth by impairing cell cycle progression and inducing cell apoptosis. Moreover, our data reveal that Linc00659 conferred drug sensitivity to colon cancer cells. Studies have increasingly indicated that some lncRNAs were associated with drug resistance in human cancer [38,39,40]. Lee et al. identified differentially expressed long intergenic noncoding RNAs (lincRNAs) in 5-FU-resistant colon cancer cells . A total of 42 lincRNA candidates were identified with differential expression in 5-FU-resistant cells. Among them, colon cancer cells with snaR loss could decrease in vitro sensitivity to 5-FU . In addition, 46 lncRNAs were identified with differential expression in vincristine-resistant colon cancer cells, compared with control cells, using next-generation sequencing profiling . The expression levels of Linc00152 conferred colon cancer cell resistance to oxaliplatin by modulating the miR-193a-3p/ERBB4/AKT signaling axis . In our study, we demonstrated that Linc00659 knockdown could accelerate cell apoptosis in colon cancer cells following oxaliplatin treatment. An analysis of the Gene Expression Omnibus database showed that the expression levels of Linc00659 were increased in drug-resistant breast and colon cancer cells [42, 43]. According to our present study, we determined a putative mechanism by which Linc00659 may play a novel oncogenic role in cell growth and drug resistance by modulating PI3K-AKT signaling.
Unlike protein-coding genes, whose transcripts are mainly exported to the cytoplasm for translation, lncRNAs are distributed in either the cytoplasm or the nucleus, depending on their function. Therefore, further examination of subcellular localization of lncRNA would facilitate determining the putative mechanism. We analyzed the subcellular localization of Linc00659, revealing that Linc00659 expression occurred both in the nucleus and cytoplasm in most colon cancer cells (Additional file 3: Figure S3). According to its location, Linc00659 might regulate colon cancer cell growth by binding with microRNA sponge, regulate transcription factors, or modulate epigenetic modification to prevent tumor suppressor gene expression. Notably, we observed that Linc00659 expression predominantly enriches in the nucleus compared with the cytoplasm in HT29 cells (Additional file 3: Figure S3). Studies have reported that the siRNA approach was less efficient for knockdown nuclear-lncRNA expression . In the present study, we were unable to explore the role of Linc00659 in HT29 cells using the same siRNA oligo. According to the varying localization of Linc00659, we suggest that Linc00659 might play an oncogenic role through the difference mechanism in HT29 cells. However, the detailed mechanism of Linc00659 remains unclear.
In this study, we identified lncRNAs with differential expression in colon cancer. We first reported that the expression levels of Linc00659 were significantly increased in colon cancer and that Linc00659 loss suppressed colon cancer cell growth by impairing cell cycle progression. Our finding also establishes Linc00659 as a prognostic biomarker for indicating CRC survival and to be used as adjuvant therapy for chemotherapy.
Shike M, Winawer SJ, Greenwald PH, Bloch A, Hill MJ, Swaroop SV. Primary prevention of colorectal cancer. The WHO collaborating Centre for the Prevention of colorectal cancer. Bull World Health Organ. 1990;68:377–85.
Eheman C, Henley SJ, Ballard-Barbash R, Jacobs EJ, Schymura MJ, Noone AM, Pan L, Anderson RN, Fulton JE, Kohler BA, et al. Annual report to the nation on the status of cancer, 1975-2008, featuring cancers associated with excess weight and lack of sufficient physical activity. Cancer. 2012;118:2338–66.
Sasaki H, Miura K, Horii A, Kaneko N, Fujibuchi W, Kiseleva L, Gu Z, Murata Y, Karasawa H, Mizoi T, et al. Orthotopic implantation mouse model and cDNA microarray analysis indicates several genes potentially involved in lymph node metastasis of colorectal cancer. Cancer Sci. 2008;99:711–9.
Cunningham D, Atkin W, Lenz HJ, Lynch HT, Minsky B, Nordlinger B, Starling N. Colorectal cancer. Lancet. 2010;375:1030–47.
Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, Margulies EH, Weng Z, Snyder M, Dermitzakis ET, Thurman RE, et al. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007;447:799–816.
Wu HH, Lin WC, Tsai KW. Advances in molecular biomarkers for gastric cancer: miRNAs as emerging novel cancer markers. Expert Rev Mol Med. 2014;16:e1.
Spizzo R, Almeida MI, Colombatti A, Calin GA. Long non-coding RNAs and cancer: a new frontier of translational research? Oncogene. 2012;31:4577–87.
Ma H, Hao Y, Dong X, Gong Q, Chen J, Zhang J, Tian W. Molecular mechanisms and function prediction of long noncoding RNA. ScientificWorldJournal. 2012;2012:541786.
Kogo R, Shimamura T, Mimori K, Kawahara K, Imoto S, Sudo T, Tanaka F, Shibata K, Suzuki A, Komune S, et al. Long noncoding RNA HOTAIR regulates polycomb-dependent chromatin modification and is associated with poor prognosis in colorectal cancers. Cancer Res. 2011;71:6320–6.
Wu ZH, Wang XL, Tang HM, Jiang T, Chen J, Lu S, Qiu GQ, Peng ZH, Yan DW. Long non-coding RNA HOTAIR is a powerful predictor of metastasis and poor prognosis and is associated with epithelial-mesenchymal transition in colon cancer. Oncol Rep. 2014;32:395–402.
Ge X, Chen Y, Liao X, Liu D, Li F, Ruan H, Jia W. Overexpression of long noncoding RNA PCAT-1 is a novel biomarker of poor prognosis in patients with colorectal cancer. Med Oncol. 2013;30:588.
Qi P, Xu MD, Ni SJ, Huang D, Wei P, Tan C, Zhou XY, Du X. Low expression of LOC285194 is associated with poor prognosis in colorectal cancer. J Transl Med. 2013;11:122.
Shi D, Zheng H, Zhuo C, Peng J, Li D, Xu Y, Li X, Cai G, Cai S. Low expression of novel lncRNA RP11-462C24.1 suggests a biomarker of poor prognosis in colorectal cancer. Med Oncol. 2014;31:31.
Zhai H, Fesler A, Schee K, Fodstad O, Flatmark K, Ju J. Clinical significance of long intergenic noncoding RNA-p21 in colorectal cancer. Clin Colorectal Cancer. 2013;12:261–6.
Yan B, Gu W, Yang Z, Gu Z, Yue X, Gu Q, Liu L. Downregulation of a long noncoding RNA-ncRuPAR contributes to tumor inhibition in colorectal cancer. Tumour Biol. 2014;35:11329–35.
Takahashi Y, Sawada G, Kurashige J, Uchi R, Matsumura T, Ueo H, Takano Y, Eguchi H, Sudo T, Sugimachi K, et al. Amplification of PVT-1 is involved in poor prognosis via apoptosis inhibition in colorectal cancers. Br J Cancer. 2013;110:164–71.
Guo Q, Zhao Y, Chen J, Hu J, Wang S, Zhang D, Sun Y. BRAF-activated long non-coding RNA contributes to colorectal cancer migration by inducing epithelial-mesenchymal transition. Oncol Lett. 2014;8:869–75.
Deng Q, He B, Gao T, Pan Y, Sun H, Xu Y, Li R, Ying H, Wang F, Liu X, et al. Up-regulation of 91H promotes tumor metastasis and predicts poor prognosis for patients with colorectal cancer. PLoS One. 2014;9:e103022.
He X, Tan X, Wang X, Jin H, Liu L, Ma L, Yu H, Fan Z. C-Myc-activated long noncoding RNA CCAT1 promotes colon cancer cell proliferation and invasion. Tumour Biol. 2014; 35:12181–8.
Ji P, Diederichs S, Wang W, Boing S, Metzger R, Schneider PM, Tidow N, Brandt B, Buerger H, Bulk E, et al. MALAT-1, a novel noncoding RNA, and thymosin beta4 predict metastasis and survival in early-stage non-small cell lung cancer. Oncogene. 2003;22:8031–41.
Yin D, He X, Zhang E, Kong R, De W, Zhang Z. Long noncoding RNA GAS5 affects cell proliferation and predicts a poor prognosis in patients with colorectal cancer. Med Oncol. 2014;31:253.
Yang F, Gao JY, Chen H, Du ZH, Zhang XQ, Gao W. Targeted inhibition of the phosphoinositide 3-kinase impairs cell proliferation, survival, and invasion in colon cancer. Onco Targets Ther. 2017;10:4413–22.
Yu Y, Yu X, Ma J, Tong Y, Yao J. Effects of NVP-BEZ235 on the proliferation, migration, apoptosis and autophagy in HT-29 human colorectal adenocarcinoma cells. Int J Oncol. 2016;49:285–93.
Wang H, Zhang L, Yang X, Jin Y, Pei S, Zhang D, Zhang H, Zhou B, Zhang Y, Lin D. PUMA mediates the combinational therapy of 5-FU and NVP-BEZ235 in colon cancer. Oncotarget. 2015;6:14385–98.
Virdee K, Parone PA, Tolkovsky AM. Phosphorylation of the pro-apoptotic protein BAD on serine 155, a novel site, contributes to cell survival. Curr Biol. 2000;10:1151–4.
Agarwal E, Brattain MG, Chowdhury S. Cell survival and metastasis regulation by Akt signaling in colorectal cancer. Cell Signal. 2013;25:1711–9.
Yang Q, Wang X, Tang C, Chen X, He J. H19 promotes the migration and invasion of colon cancer by sponging miR-138 to upregulate the expression of HMGA1. Int J Oncol. 2017;50:1801–9.
Lv Y, Li H, Li F, Liu P, Zhao X. Long noncoding RNA MNX1-AS1 knockdown inhibits cell proliferation and migration in ovarian cancer. Cancer Biother Radiopharm. 2017;32:91–9.
McCleland ML, Mesh K, Lorenzana E, Chopra VS, Segal E, Watanabe C, Haley B, Mayba O, Yaylaoglu M, Gnad F, Firestein R. CCAT1 is an enhancer-templated RNA that predicts BET sensitivity in colorectal cancer. J Clin Invest. 2016;126:639–52.
Tseng YY, Moriarity BS, Gong W, Akiyama R, Tiwari A, Kawakami H, Ronning P, Reuland B, Guenther K, Beadnell TC, et al. PVT1 dependence in cancer with MYC copy-number increase. Nature. 2014;512:82–6.
Yang L, Xu L, Wang Q, Wang M, An G. Dysregulation of long non-coding RNA profiles in human colorectal cancer and its association with overall survival. Oncol Lett. 2016;12:4068–74.
Henry WS, Hendrickson DG, Beca F, Glass B, Lindahl-Allen M, He L, Ji Z, Struhl K, Beck AH, Rinn JL, Toker A. LINC00520 is induced by Src, STAT3, and PI3K and plays a functional role in breast cancer. Oncotarget. 2016;7:81981–94.
Li DD, Fu ZQ, Lin Q, Zhou Y, Zhou QB, Li ZH, Tan LP, Chen RF, Liu YM. Linc00675 is a novel marker of short survival and recurrence in patients with pancreatic ductal adenocarcinoma. World J Gastroenterol. 2015;21:9348–57.
Stojic L, Niemczyk M, Orjalo A, Ito Y, Ruijter AE, Uribe-Lewis S, Joseph N, Weston S, Menon S, Odom DT, et al. Transcriptional silencing of long noncoding RNA GNG12-AS1 uncouples its transcriptional and product-related functions. Nat Commun. 2016;7:10406.
Tian Z, Wen S, Zhang Y, Shi X, Zhu Y, Xu Y, Lv H, Wang G. Identification of dysregulated long non-coding RNAs/microRNAs/mRNAs in TNM I stage lung adenocarcinoma. Oncotarget. Oncotarget. 2017;8:51703–18.
Chen X, Liu B, Yang R, Guo Y, Li F, Wang L, Hu H. Integrated analysis of long non-coding RNAs in human colorectal cancer. Oncotarget. 2016;7:23897–908.
Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, Huarte M, Zuk O, Carey BW, Cassady JP, et al. Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009;458:223–7.
Tong YS, Zhou XL, Wang XW, Wu QQ, Yang TX, Lv J, Yang JS, Zhu B, Cao XF. Association of decreased expression of long non-coding RNA LOC285194 with chemoradiotherapy resistance and poor prognosis in esophageal squamous cell carcinoma. J Transl Med. 2014;12:233.
Lee H, Kim C, Ku JL, Kim W, Yoon SK, Kuh HJ, Lee JH, Nam SW, Lee EK. A long non-coding RNA snaR contributes to 5-fluorouracil resistance in human colon cancer cells. Mol Cells. 2014;37:540–6.
Sun QL, Zhao CP, Wang TY, Hao XB, Wang XY, Zhang X, Li YC. Expression profile analysis of long non-coding RNA associated with vincristine resistance in colon cancer cells by next-generation sequencing. Gene. 2015;572:79–86.
Yue B, Cai D, Liu C, Fang C, Yan D. Linc00152 functions as a competing endogenous RNA to confer Oxaliplatin resistance and holds prognostic values in colon cancer. Mol Ther. 2016;24:2064–77.
Selga E, Morales C, Noe V, Peinado MA, Ciudad CJ. Role of caveolin 1, E-cadherin, enolase 2 and PKCalpha on resistance to methotrexate in human HT29 colon cancer cells. BMC Med Genet. 2008;1:35.
Liu L, Greger J, Shi H, Liu Y, Greshock J, Annan R, Halsey W, Sathe GM, Martin AM, Gilmer TM. Novel mechanism of lapatinib resistance in HER2-positive breast tumor cells: activation of AXL. Cancer Res. 2009;69:6871–8.
Lennox KA, Behlke MA. Cellular localization of long non-coding RNAs affects silencing by RNAi more than by antisense oligonucleotides. Nucleic Acids Res. 2016;44:863–77.
We acknowledge Wallace Academic Editing for editing this manuscript.
This work was supported by funding from the Ministry of Science and Technology (MOST105–2314-B-010-051), Kaohsiung Veterans General Hospital (VGHKS 105–135 and VGHKS105–092) and Zuoying Branch of Kaohsiung Armed Forces General Hospital (grant no. ZBH-105-11).
Ethics approval and consent to participate
This study was approved by the ethics committee of the Taipei Veterans General Hospital and written informed consents were obtained from all patients (IRB number: 2015–06-012CC).
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Knockdown of Linc00659 expression with siRNA suppressed growth of colon cancer cells (A)-(D) The expression levels of Linc00659 were examined after transfection of Lovo, LS174T, DLD-1 and SW620 with siRNA. (E)-(H) After transient transfection of Lovo, LS174T, DLD-1 and SW620 cells with si-Linc00659, cell proliferation was examined; the graph shows quantified values. (DOCX 124 kb)
Figure S2. Stable knockdown of Linc00659 expression with shLinc00659#2 suppressed HCT116 cell growth. (A) The expression levels of Linc00659 were examined in HCT116 cells with stable shLinc00659#2 knockdown. (B) and (C) The cell proliferation and colony formation assays were performed: (D) graph illustrating quantified values. (E) The effects of Linc00659 knockdown on anchorage-independent growth of HCT116 were examined by a soft agar formation assay: (F) graph illustrating quantified values. Data are reported as the number of colonies relative to the control (means ± standard deviation). (DOCX 205 kb)
Figure S3. Examination of the subcellular fractionation localization of Linc00659 in CRC cell lines. After nuclear and cytosolic separation, total RNA from Lovo, HCT116, HT29, SW620 and DLD-1 cells underwent RT and real-time PCR. GAPDH was used as a cytosol marker (A) and U6 was used as a nucleus marker (B). (C) RNA expression levels of Linc00659 candidates in the nucleus and cytoplasm were measured by real-time PCR, respectively. CRC, colorectal cancer; GAPDH, glyceraldehyde 3-phosphate dehydrogenase; PCR, polymerase chain reaction. (DOCX 109 kb)
About this article
Cite this article
Tsai, KW., Lo, YH., Liu, H. et al. Linc00659, a long noncoding RNA, acts as novel oncogene in regulating cancer cell growth in colorectal cancer. Mol Cancer 17, 72 (2018). https://doi.org/10.1186/s12943-018-0821-1
- Colorectal carcinoma
- Noncoding RNA