A computationally constructed ceRNA interaction network based on a comparison of the SHEE and SHEEC cell lines
- Jiachun Sun†1,
- Junqiang Yan†2,
- Xiaozhi Yuan1,
- Ruina Yang1,
- Tanyou Dan1,
- Xinshuai Wang1,
- Guoqiang Kong1 and
- Shegan Gao1Email author
© The Author(s) 2016
Received: 17 August 2016
Accepted: 10 September 2016
Published: 26 September 2016
Long non-coding RNAs (lncRNAs) play critical and complicated roles in the regulation of various biological processes, including chromatin modification, transcription and post-transcriptional processing. Interestingly, some lncRNAs serve as miRNA “sponges” that inhibit interaction with miRNA targets in post-transcriptional regulation. We constructed a putative competing endogenous RNA (ceRNA) network by integrating lncRNA, miRNA and mRNA expression based on high-throughput RNA sequencing and microarray data to enable a comparison of the SHEE and SHEEC cell lines. Using Targetscan and miRanda bioinformatics algorithms and miRTarbase microRNA-target interactions database, we established that 51 miRNAs sharing 13,623 MREs with 2260 genes and 82 lncRNAs were involved in this ceRNA network. Through a biological function analysis, the ceRNA network appeared to be primarily involved in cell proliferation, apoptosis, the cell cycle, invasion and metastasis. Functional pathway analyses demonstrated that the ceRNA network potentially modulated multiple signaling pathways, such as the MAPK, Ras, HIF-1, Rap1, and PI3K/Akt signaling pathways. These results might provide new clues to better understand the regulation of the ceRNA network in cancer.
KeywordsSHEE SHEEC miRNA lncRNA ceRNA
Human esophageal cancer is the sixth leading cause of cancer-related death worldwide . The most frequent subtype, esophageal squamous cell carcinoma (ESCC), has an obvious geographic variation in its distribution, with eastern Asia being a high incidence area [1, 2]. Most ESCC patients are diagnosed at an advanced stage, past the best opportunity for surgical treatment, and their overall 5-year survival rate is less than 15 %. By contrast, patients diagnosed at an early stage can have 10-year overall post-operative survival rates up to 95 % .
The immortalized human esophageal epithelial cell line SHEE and the malignantly transformed esophageal carcinoma cell line SHEEC were established and cultured in consecutive passes by Zeng Yi of the National Institute for Viral Disease Control and Prevention and Shen Zhongying of the Medical College of Shantou University [4, 5]. From a biological perspective, SHEE is similar to a primary cell line and retains its proliferation ability and differentiation potential. It retains the phenotype of primary epithelial cells, such as growth in a monolayer and anchorage-dependent cell aggregation without colony formation in soft agar or tumor formation after transplantation. SHEEC was established through induction of the SHEE cell line using 12-otetradeanoy-lphorbol-13-acetate (TPA). These two cell lines are the most appropriate models for ESCC tumorigenesis.
To date, considerable research has indicated that differential expression of miRNAs has been observed in various types of human cancers, such as lung cancer, gastric cancer, breast cancer, colorectal cancer and ESCC [6–9]. The expression of miRNAs has been used as a useful biomarker for tumor prognosis and clinical treatment. MicroRNAs can also regulate gene expression at the epigenetic regulation level. They are involved in most aspects of cellular processes, including development, differentiation, growth and apoptosis [10, 11].
Unlike small ncRNAs, increasing evidence indicates that lncRNAs play critical and complicated roles in the regulation of various biological processes, including chromatin modification, transcription and post-transcriptional processing [12–15]. Interestingly, additional studies have revealed that some lncRNAs can serve as miRNA “sponges” that inhibit interaction with their miRNA targets in post-transcriptional regulation. This means they can be considered competing endogenous RNA (ceRNA). Through gene expression data analysis, Sumazin et al. found >7000 transcripts acting as ceRNAs in glioblastoma tissues . Cesana et al.  revealed the ceRNA function of linc-MD1, which liberates the differentiation factors MEF2C and MAML1 from repression by decoying miR-135 and -133 to control muscle differentiation. Another lncRNA, linc-RoR, acts as an endogenous sponge to mediate miR-145 regulation of self-renewal of human embryonic stem cells [18, 19]. Interestingly, CHRF, a lncRNA connected to cardiac hypertrophy-related factor, was found to act as an endogenous sponge of miR-489, directly binding to and downregulating miR-489 expression levels and, in turn, regulating its target Myd88 expression and hypertrophy .
The above studies identified a new regulatory mechanism in post-transcriptional regulation: ceRNA networks. The mechanism of the ceRNA network is that all types of RNA transcripts (lncRNA, pseudogenes and circular RNAs) could communicate with each other by competing for binding to shared miRNA-binding sites (MREs) .
In this study, we constructed a putative competing endogenous RNA (ceRNA) network by integrating lncRNA, miRNA and mRNA expression based on high-throughput RNA sequencing and microarray data to enable a comparison of the SHEE and SHEEC cell lines. Using Targetscan and miRanda bioinformatics algorithms and the miRTarbase microRNA–target interactions database, we established that 51 miRNAs sharing 13,623 MREs with 2260 genes and 82 lncRNAs were involved in this ceRNA network. Based on a biological function analysis, this ceRNA network may participate in the PI3K/Akt pathway and, consistent with previous reports, may play a modulating role in the regulation of stem-like cells in primary ESCC . These results might provide new clues to better understand the regulation of the ceRNA network in cancer.
Materials and methods
SHEE and SHEEC cells were cultured in Gibco MEM medium supplemented with 100 ml/l fetal bovine serum (containing 100 μg/ml penicillin and 100 μg/ml streptomycin) and incubated at 37 °C in a humidified atmosphere of 50 ml/l CO2. Cells were harvested after growth into a full monolayer and were kept at −70 °C until use.
Total RNA isolation
Total RNA from each sample was extracted using a TRK-1001 Total RNA Purification Kit (LC Sciences) according to the manufacturer’s protocol. Total RNA was quantified on a NanoDrop ND-2000 (Thermo Scientific), and the RNA integrity was assessed using an Agilent Bioanalyzer 2100 (Agilent Technologies).
μParaflo microRNA microarray assay
The microarray assay was performed by a service provider (LC Sciences). The assay started with a 2 to 5 μg total RNA sample that was 3′-extended with a poly(A) tail using poly(A) polymerase. An oligonucleotide tag was then ligated to the poly(A) tail for later fluorescent dye staining. Hybridization was performed overnight on a μParaflo microfluidic chip using a microcirculation pump (Atactic Technologies) [23, 24]. On the microfluidic chip, each detection probe consisted of a chemically modified nucleotide-coding segment complementary to the target microRNA (from miRBase 20.0, http://www.mirbase.org/) or other RNA (control or customer-defined sequences) and a spacer segment of polyethylene glycol to extend the coding segment away from the substrate.
The detection probes were made via in situ synthesis using PGR (photogenerated reagent) chemistry. The hybridization melting temperatures were balanced by chemical modifications of the detection probes. Hybridization used 100 l 6× SSPE buffer consisting of 0.90 M NaCl, 60 mM Na2HPO4, 6 mM EDTA (pH 6.8) and 25 % formamide at 34 °C. After RNA hybridization, tag-conjugating Cy3 dye was circulated through the microfluidic chip for dye staining. Fluorescence images were collected using a laser scanner (GenePix 4000B, Molecular Device) and digitized using Array-Pro image analysis software (Media Cybernetics).
Data were analyzed by first subtracting the background, and then normalizing the signals using a LOWESS filter (locally weighted regression) . The miRNA differential expression based on the normalized signal was analyzed via selective application of Student’s t test. The significance threshold was p < 0.05 and fold-change >2.
Construction of cDNA libraries and high-throughput sequencing
To construct the next-generation sequencing libraries, approximately 3 μg of total RNA was used to deplete ribosomal RNA according to the manufacturer’s protocol for the Human/Mouse/Rat Ribo-Zero rRNA Removal Kit (Epicentre/Illumina). Following purification, the poly(A)- or poly(A) + RNA fractions were fragmented into small pieces using divalent cations at elevated temperatures. The cleaved RNA fragments were reverse transcribed to construct a cDNA library using the dUTP method as described previously . The average insert size for the paired-end libraries was 300 bp (±50 bp). RNA libraries were then sequenced on the Illumina HiSeq 2500 platform using 125-bp paired-end reads at LC Biotech in Hangzhou, China.
The RNAseq data were aligned to hg19 using TopHat v2.0.9  with the default parameters. The mapped reads were assembled using Cufflinks v2.11 [28, 29]. All multiple assembled transcript files (GTF format) were then merged to produce a unique transcriptome set using the Cuffmerge utility provided in the Cufflinks package .
We filtered the assembled novel transcripts from the two pooled cell lines to obtain putative lncRNAs. First, identical and overlapping transcripts were merged to remove redundancy. Then, transcripts overlapping with known gene exons were removed. Only transcripts with a length >200 nt were retained. To identify and eliminate potential known lncRNA transcripts, we compared the merged transcriptome with lncRNA and protein-coding genes in the public authoritative database GENCODE .
To obtain a reliable dataset of putative lncRNAs, single exon models were filtered out. Next, we removed transcripts that were likely to be assembly artifacts or PCR run-on fragments according to their class code (annotated by Cuffmerge). Among the different classes defined by Cufflinks , only those annotated by “u”, “i”, “x”, “o” and “j” were retained. Extremely low gene expression is generally considered to be transcriptional noise . On average, 79 % of the initial reads with a quality score >30 could be aligned to the hg19 assembly of the human genome sequence. We used CPC (http://cpc.cbi.pku.edu.cn)  and the coding non-coding index (CNCI)  to assess the protein-coding potential of each novel transcript. Those putative transcripts with a CPC score < −1 and a CNCI score < −1 were retained as candidate lncRNAs for further analysis.
Transcript differential expression analysis
Expression levels of all of the transcripts, including putative lncRNAs and mRNAs, were quantified as fragments per kilobase of exon per million fragments mapped (FPKM) using the Cuffdiff program from the Cufflinks package. Differential gene expression was determined using Cuffdiff with a p-value of <0.05 and q-value of <0.05.
Construction of ceRNA network
The putative miRNA–lncRNA interactions were evaluated using the algorithms of Targetscan version 6.2 (http://www.targetscan.org/) and miRanda version 3.3a (http://www.microrna.org/microrna/home.do). The miRNA binding-site prediction in lncRNAs was based on their full-length sequence in consideration of their non-coding properties. High-confidence miRNA–lncRNA pairs had a Targetscan context+ score percentile >50 and miRanda max energy < −20. To reduce the false positives, at least two miRNA binding-sites were retained with each lncRNA. The mRNAs that were targeted by miRNAs with experimental support were from miRTarbase (http://mirtarbase.mbc.nctu.edu.tw/). The ceRNA relationships were integrated using an in-house Perl script. The information including all of the above interactions was imported into Cytoscape software version 3.3.0 (http://www.cytoscape.org) to construct a regulatory network.
Aberrantly expressed miRNAs
Aberrantly expressed lncRNAs and protein-coding genes
To systematically characterize the transcriptional changes between the SHEE and SHEEC cell lines, we performed high-throughput RNA sequencing according to previously published methods [36–38]. Briefly, we generated >17 million reads from strand-specific, paired-end, 125-bp RNAseq reads of total RNA depleted of rRNA from both cell populations. Using TopHat , an average of 95 % clean reads were mapped to the human GRCh38 genome (http://www.gencodegenes.org/releases/21.html). Using the lncRNA discovery pipeline described in the Methods section above, we detected 11,557 unique protein-coding genes and 15,932 unique lncRNAs that met the criteria of FPKM > 1 and length >200 nucleotides.
The expression levels of all transcripts were then quantified as fragments per kilobase of exon per million fragments mapped (FPKM) using the Cuffdiff program from the Cufflinks package . To examine the aberrantly expressed protein-coding genes and lncRNAs, we performed differential analysis using the Cuffdiff program. By setting stringent criteria (FPKM > 3, fold-change >2, and p < 0.05), 5593 protein-coding genes and 6294 lncRNAs were found to be aberrantly expressed within the group comparison (Additional files 2, 3 and 4). Interestingly, the majority of the aberrantly expressed novel lncRNAs are downregulated in our dataset. The potential mechanism of dysregulation remains to be further explored.
miRNA-binding site prediction
To establish the lncRNA–miRNA–mRNA (ceRNA) network, Targetscan 6.2 and miRanda 3.3a were used for lncRNA targeted search. Using a high confidence score (see the Materials and Methods section), the predicted MREs showed that 54 miRNAs may interact with 83 lncRNAs (Additional file 5). Based on those 54 miRNAs, we used the experimentally validated microRNA–target interactions database (miRTarbase) to search for the miRNAs’ mRNA targets. The results showed that 51 miRNAs can interact with 2260 targets (Additional file 6). Most of these targets are cancer-associated genes such as PTEN, STAT3, VEGFA, KRAS, TP53, CCND1, CDK6, E2F1, FGFR1 and EGFR, with roles in cell proliferation, apoptosis, cell cycle, invasion and metastasis.
lncRNA–miRNA–mRNA ceRNA network construction
The ceRNA network was visualized by importing the above interactions into the Cytoscape software to assemble the regulation network. The genes within the above networks were further processed by gene function (GO; http://geneontology.org/) and KEGG (http://www.genome.jp/kegg/) pathway analysis.
Gene function and pathway analysis
Accumulated studies have revealed that ceRNAs could serve as post-transcriptional regulators of protein-coding gene expression by decoying miRNAs from other target transcripts, such as lncRNAs, mRNA, pseudogenes and circular RNAs (circRNAs) [41–43]. Some studies have confirmed that ceRNAs play an important role in the regulation of gene expression in cancers such as head and neck squamous cell carcinoma, prostate cancer, papillary thyroid carcinoma, pituitary gonadotrope tumors, ovarian cancer, and chronic lymphocytic leukemia [44–46].
In this study, based on the high-throughput RNA sequencing and microarray data, we constructed a putative ceRNA network by integrating lncRNA, miRNA and mRNA expression. In our ceRNA network, lncRNA UCA1 was reported to be an independent prognostic factor associated with tumor differentiation and location .
Subsequent studies reported that UCA1 could upregulate 3 target genes of miR-204-5p (BCL2, RAB22A and CREB1) though competitive sponging of miR-204-5p, and then promote proliferation and chemoresistance in CRC cells. . Yang et al. reported that the overexpressed UCA1 was correlated with metastasis of epithelial ovarian cancer (EOC) and functioned as a ceRNA to suppress the expression of matrix metallopeptidase 14 (MMP14) via competition for miR-485-5p. . Chen et al. reported on nuclear enriched abundant transcript 1 (NEAT1) lncRNA as a novel prognostic indicator for patients with ESCC, finding that it contributes to the malignant character of ESCC . NEAT1 was also reported to be a well sponge platform for many kinds of miRNA, such as miR-548 in the regulation of breast cancer cell apoptosis , miR-204 in regulation of epithelial-to-mesenchymal transition (EMT) and the radioresistance of NPC cells , and has-mir-98-5p in regulation of EGCG-induced CTR1 and cDDP sensitivity enhancement in NSCLC . These specific features could potentially be used to classify lncRNAs and identify those that participate in ceRNA networks.
We found that biological function was significantly enriched with signaling pathways, small molecule metabolic processes, apoptotic processes, small GTPase-mediated signal transduction, negative regulation of the apoptotic process, mitotic cell cycle, protein binding, ATP binding, DNA binding, protein kinase binding, protein kinase activity and GTP binding, among other things (Fig. 3; Additional file 8). These signaling pathways were often altered, especially in the cancer cells, resulting in phenotypes of uncontrolled growth and increased capability to invade the surrounding tissue. These crucial molecules involved in signaling pathways represent attractive targets for cancer therapy [54–57]. Agents targeting epidermal growth factor receptor (EGFR), PI3K, and mTOR have been developed for interfering with their signaling functions.
Previous studies demonstrated that the PTEN/PI3K/Akt pathway was essential to side population (SP) cells thanks to its involvement in the regulation of ABCG2 transporter function in primary ESCCs . The PTEN/PI3K/Akt pathway plays a modulating role in regulating stem-like cells in primary ESCCs and may provide essential clues for the development of novel therapeutic strategies and efficient drugs.
Consistent with these studies, we found the PI3K/Akt signaling pathway was the most significantly enriched pathway, based on KEGG pathway analysis (Fig. 5; Additional file 9). PTEN, a well-researched protein in cancer, encodes a plasma-membrane lipid phosphatase that functions as a negative regulator of the PI3K/Akt signaling pathway [58–60].
There is more evidence that both non-coding and protein-coding transcripts regulate PTEN levels via PTEN ceRNAs and then antagonize downstream PI3K/Akt signaling [61–64, 16]. Our analysis suggested that lncRNAs like NEAT1 and TCONS_00287673 may upregulate PTEN though competitive sponging of miR-26a-5p and miR-182-5p. The regulatory mechanism of these ceRNAs needs to be further evaluated.
Competing endogenous RNAs
Coding non-coding index
Coding potential calculator
Esophageal squamous cell carcinoma
Fragments per kilobase of exon per million fragments mapped
Long non-coding RNAs
We would like to acknowledge anonymous reviewers and our academic editor for their constructive suggestions on the manuscript. We would like to sincerely thank Dr. Jianning Liu of LC Biotech for his help with the bioinformatics analysis.
This research was financially supported by a National Natural Science Foundation of China (NSFC) grant funded by the government of China (U1304809).
Availability of data and materials
The microarray and RNAseq data sets were submitted to the Gene Expression Omnibus (GEO) database under accession numbers GSE72138 and GSE72273, respectively.
JCS and SGG contributed to design the experiment. JCS and JQY performed the experiments and drafted the manuscript. XZY, RNY and TYD contributed to the analysis and interpretation of the data. XSW and GQK contributed to analysis of the data. JCS, JQY and SGG contributed to revise the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin. 2011. doi:10.3322/caac.20107.Google Scholar
- Kano M, Seki N, Kikkawa N, Fujimura L, Hoshino I, Akutsu Y, Chiyomaru T, Enokida H, Nakagawa M, Matsubara H. miR-145, miR-133a and miR-133b: Tumor-suppressive miRNAs target FSCN1 in esophageal squamous cell carcinoma. Int J Cancer. 2010. doi:10.1002/ijc.25284.PubMedGoogle Scholar
- Kranzfelder M, Buchler P, Friess H. Surgery within multimodal therapy concepts for esophageal squamous cell carcinoma (ESCC): the MRI approach and review of the literature. Adv Med Sci. 2009. doi:10.2478/v10039-009-0044-1.PubMedGoogle Scholar
- Shen ZY, Cen S, CW, et al. Immortalization of human fetal esophageal epithelial cells induced by E6 and E7 genes of human papilloma virus 18. Chinese J Exp Clin Virol 1999;13:121–4
- Shen ZY, Xu LY, Chen XH, Cai WJ, Shen J, Chen JY, Huang TH, Zeng Y. The genetic events of HPV-immortalized esophageal epithelium cells. Int J Mol Med. 2001;8:537–42.PubMedGoogle Scholar
- Wang XY, Wu MH, Liu F, Li Y, Li N, Li GY, Shen SR. Differential miRNA expression and their target genes between NGX6-positive and negative colon cancer cells. Mol Cell Biochem. 2010. doi:10.1007/s11010-010-0582-7.PubMed CentralGoogle Scholar
- Feng XL, Wang ZM, Fillmore R, Xi YG. MiR-200, a new star miRNA in human cancer. Cancer Lett. 2014. doi:10.1016/j.canlet.2013.11.004.PubMed CentralGoogle Scholar
- Wang J, Wang H, Zhang Y, Zhen N, Zhang L, Qiao Y, Weng W, Liu X, Ma L, Xiao W, Yu W, Chu Q, Pan Q, Sun F. Mutual inhibition between YAP and SRSF1 maintains long non-coding RNA, Malat1-induced tumourigenesis in liver cancer. Cell Signal. 2014. doi:10.1016/j.cellsig.2014.01.022.Google Scholar
- Reddy KB. MicroRNA (miRNA) in cancer. Cancer Cell Int. 2015. doi:10.1186/S12935-015-0185-1.PubMedPubMed CentralGoogle Scholar
- Huang Y, Shen XJ, Zou QA, Wang SP, Tang SM, Zhang GZ. Biological functions of microRNAs: a review. J Physiol Biochem. 2011. doi:10.1007/s13105-010-0050-6.Google Scholar
- Huang Y, Shen XJ, Zou QA, Zhao QL. Biological functions of MicroRNAs. Russ J Bioorg Chem. 2010. doi:10.1134/S1068162010060026.Google Scholar
- Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009. doi:10.1038/Nrg2521.PubMedGoogle Scholar
- Chang THHY. Long noncoding RNA in genome regulation: prospects and mechanisms. RNA Biol. 2010;7:582–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Cao J. The functional role of long non-coding RNAs and epigenetics. Biological procedures online. 2014. doi:10.1186/1480-9222-16-11.PubMedPubMed CentralGoogle Scholar
- Halley P, Kadakkuzha BM, Faghihi MA, Magistri M, Zeier Z, Khorkova O, Coito C, Hsiao J, Lawrence M, Wahlestedt C. Regulation of the apolipoprotein gene cluster by a long noncoding RNA. Cell Rep. 2014. doi:10.1016/j.celrep.2013.12.015.PubMedPubMed CentralGoogle Scholar
- Sumazin P, Yang XR, Chiu HS, Chung WJ, Iyer A, Llobet-Navas D, Rajbhandari P, Bansal M, Guarnieri P, Silva J, Califano A. An extensive MicroRNA-mediated network of RNA-RNA interactions regulates established oncogenic pathways in glioblastoma. Cell. 2011. doi:10.1016/j.cell.2011.09.041.PubMedPubMed CentralGoogle Scholar
- Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, Tramontano A, Bozzoni I. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011. doi:10.1016/j.cell.2011.09.028.PubMedGoogle Scholar
- Wang Y, Xu Z, Jiang J, Xu C, Kang J, Xiao L, Wu M, Xiong J, Guo X, Liu H. Endogenous miRNA sponge lincRNA-RoR regulates Oct4, Nanog, and Sox2 in human embryonic stem cell self-renewal. Dev Cell. 2013. doi:10.1016/j.devcel.2013.03.002.Google Scholar
- Zhou X, Gao Q, Wang J, Zhang X, Liu K, Duan Z. Linc-RNA-RoR acts as a “sponge” against mediation of the differentiation of endometrial cancer stem cells by microRNA-145. Gynecol Oncol. 2014. doi:10.1016/j.ygyno.2014.02.033.Google Scholar
- Wang K, Liu F, Zhou LY, Long B, Yuan SM, Wang Y, Liu CY, Sun T, Zhang XJ, Li PF. The long noncoding RNA CHRF regulates cardiac hypertrophy by targeting miR-489. Circ Res. 2014. doi:10.1161/CIRCRESAHA.114.302476.Google Scholar
- Thomson DW, Dinger ME. Endogenous microRNA sponges: evidence and controversy. Nat Rev Genet. 2016. doi:10.1038/nrg.2016.20.PubMedGoogle Scholar
- Shi H, Gu YC, Yang J, Xu L, Mi WY, Yu WG. Lipocalin 2 promotes lung metastasis of murine breast cancer cells. J Exp Clin Canc Res. 2008. doi:10.1186/1756-9966-27-83.Google Scholar
- Gao XL, Gulari E, Zhou XC. In situ synthesis of oligonucleotide microarrays. Biopolymers. 2004. doi:10.1002/Bip.20005.PubMedGoogle Scholar
- Zhu Q, Hong A, Sheng N, Zhang X, Jun K-Y, Srivannavit O, Gulari E, Gao X, Zhou X. Microfluidic biochip for nucleic acid and protein analysis. Methods Mol Biol. 2007;382:287–312.View ArticlePubMedGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003. doi:10.1093/bioinformatics/19.2.185.PubMedGoogle Scholar
- Parkhomchuk D, Borodina T, Amstislavskiy V, Banaru M, Hallen L, Krobitsch S, Lehrach H, Soldatov A. Transcriptome analysis by strand-specific sequencing of complementary DNA. Nucleic Acids Res. 2009. doi:10.1093/nar/gkp596.PubMedPubMed CentralGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009. doi:10.1093/bioinformatics/btp120.PubMedPubMed CentralGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010. doi:10.1038/Nbt.1621.PubMedPubMed CentralGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, Pimentel H, Salzberg SL, Rinn JL, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012. doi:10.1038/nprot.2012.016.PubMedPubMed CentralGoogle Scholar
- Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008. doi:10.1038/Nmeth.1226.PubMedGoogle Scholar
- Harrow J, Frankish A, Gonzalez JM, Tapanari E, Diekhans M, Kokocinski F, Aken BL, Barrell D, Zadissa A, Searle S, Barnes I, Bignell A, Boychenko V, Hunt T, Kay M, Mukherjee G, Rajan J, Despacio-Reyes G, Saunders G, Steward C, Harte R, Lin M, Howald C, Tanzer A, Derrien T, Chrast J, Walters N, Balasubramanian S, Pei BK, Tress M, Rodriguez JM, Ezkurdia I, van Baren J, Brent M, Haussler D, Kellis M, Valencia A, Reymond A, Gerstein M, Guigo R, Hubbard TJ. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 2012. doi:10.1101/gr.135350.111.Google Scholar
- Roberts A, Pimentel H, Trapnell C, Pachter L. Identification of novel transcripts in annotated genomes using RNA-Seq. Bioinformatics. 2011. doi:10.1093/bioinformatics/btr355.Google Scholar
- Hebenstreit D, Fang MQ, Gu MX, Charoensawan V, van Oudenaarden A, Teichmann SA. RNA sequencing reveals two major classes of gene expression levels in metazoan cells. Mol Syst Biol. 2011. doi:10.1038/Msb.2011.28.PubMedPubMed CentralGoogle Scholar
- Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007. doi:10.1093/Nar/Gkm391.Google Scholar
- Sun L, Luo HT, Bu DC, Zhao GG, Yu KT, Zhang CH, Liu YN, Chen RS, Zhao Y. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013. doi:10.1093/nar/gkt646.Google Scholar
- Basu S, Muller F, Sanges R. Examples of sequence conservation analyses capture a subset of mouse long non-coding RNAs sharing homology with fish conserved genomic elements. BMC Bioinformatics. 2013. doi:10.1186/1471-2105-14-S7-S14.PubMedPubMed CentralGoogle Scholar
- Alvarez-Dominguez JR, Hu WQ, Yuan BB, Shi JH, Park SS, Gromatzky AA, van Oudenaarden A, Lodish HF. Global discovery of erythroid long noncoding RNAs reveals novel regulators of red cell maturation. Blood. 2014. doi:10.1182/blood-2013-10-530683.PubMedPubMed CentralGoogle Scholar
- Jiang W, Liu Y, Liu R, Zhang K, Zhang Y. The lncRNA DEANR1 facilitates human endoderm differentiation by activating FOXA2 expression. Cell Rep. 2015. doi:10.1016/j.celrep.2015.03.008.Google Scholar
- Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of Gene Ontology categories in Biological Networks. Bioinformatics. 2005. doi:10.1093/bioinformatics/bti551.PubMedGoogle Scholar
- Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011. doi:10.1016/j.cell.2011.07.014.Google Scholar
- de Giorgio A, Krell J, Harding V, Stebbing J, Castellano L. Emerging roles of competing endogenous RNAs in cancer: insights from the regulation of PTEN. Mol Cell Biol. 2013. doi:10.1128/MCB.00683-13.Google Scholar
- Karreth FA, Pandolfi PP. ceRNA cross-talk in cancer: when ce-bling rivalries go awry. Cancer Discovery. 2013. doi:10.1158/2159-8290.CD-13-0202.PubMedPubMed CentralGoogle Scholar
- Hansen TB, Kjems J, Damgaard CK. Circular RNA and miR-7 in cancer. Cancer Res. 2013. doi:10.1158/0008-5472.CAN-13-1568.PubMed CentralGoogle Scholar
- Qu S, Yang X, Li X, Wang J, Gao Y, Shang R, Sun W, Dou K, Li H. Circular RNA: a new star of noncoding RNAs. Cancer Lett. 2015. doi:10.1016/j.canlet.2015.06.003.Google Scholar
- Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000. doi:10.1016/S0092-8674(00)81683-9.PubMed CentralGoogle Scholar
- Bianco R, Melisi D, Ciardiello F, Tortora G. Key cancer cell signal transduction pathways as therapeutic targets. Eur J Cancer. 2006. doi:10.1016/j.ejca.2005.07.034.PubMedGoogle Scholar
- Li JY, Ma X, Zhang CB. Overexpression of long non-coding RNA UCA1 predicts a poor prognosis in patients with esophageal squamous cell carcinoma. Int J Clin Exp Pathol. 2014;7:7938–44.PubMedPubMed CentralGoogle Scholar
- Bian Z, Jin L, Zhang J, Yin Y, Quan C, Hu Y, Feng Y, Liu H, Fei B, Mao Y, Zhou L, Qi X, Huang S, Hua D, Xing C, Huang Z. LncRNA-UCA1 enhances cell proliferation and 5-fluorouracil resistance in colorectal cancer by inhibiting miR-204-5p. Sci Rep. 2016. doi:10.1038/srep23892.Google Scholar
- Yang Y, Jiang Y, Wan Y, Zhang L, Qiu J, Zhou S, Cheng W. UCA1 functions as a competing endogenous RNA to suppress epithelial ovarian cancer metastasis. Tumour Biol. 2016. doi:10.1007/s13277-016-4917-1.PubMed CentralGoogle Scholar
- Chen X, Kong J, Ma Z, Gao S, Feng X. Up regulation of the long non-coding RNA NEAT1 promotes esophageal squamous cell carcinoma cell progression and correlates with poor prognosis Am. J Cancer Res. 2015;5:2808–15.View ArticleGoogle Scholar
- Ke H, Zhao L, Feng X, Xu H, Zou L, Yang Q, Su X, Peng L, Jiao B. NEAT1 is required for survival of breast cancer cells through FUS and miR-548 gene. Regul Syst Bio. 2016. doi:10.4137/GRSB.S29414.Google Scholar
- Lu Y, Li T, Wei G, Liu L, Chen Q, Xu L, Zhang K, Zeng D, Liao R. The long non-coding RNA NEAT1 regulates epithelial to mesenchymal transition and radioresistance in through miR-204/ZEB1 axis in nasopharyngeal carcinoma. Tumour Biol. 2016. doi:10.1007/s13277-015-4773-4.Google Scholar
- Jiang P, Wu X, Wang X, Huang W, Feng Q. NEAT1 upregulates EGCG-induced CTR1 to enhance cisplatin sensitivity in lung cancer cells Oncotarget 2016. doi:10.18632/oncotarget.9712.
- Hainaut P, Plymoth A. Targeting the hallmarks of cancer: towards a rational approach to next-generation cancer therapy. Curr Opin Oncol. 2013. doi:10.1097/Cco.0b013e32835b651e.Google Scholar
- Ma X, Ma CX, Wang J. Endometrial carcinogenesis and molecular signaling pathways. American Journal of Molecular Biology. 2014. doi:10.4236/ajmb.2014.43015.Google Scholar
- Groen JN, Capraro D, Morris KV. The emerging role of pseudogene expressed non-coding RNAs in cellular functions. Int J Biochem Cell Biol. 2014. doi:10.1016/j.biocel.2014.05.008.PubMedPubMed CentralGoogle Scholar
- Hopkins BD, Hodakoski C, Barrows D, Mense SM, Parsons RE. PTEN function: the long and the short of it. Trends Biochem Sci. 2014. doi:10.1016/j.tibs.2014.02.006.PubMedPubMed CentralGoogle Scholar
- Salmena L, Carracedo A, Pandolfi PP. Tenets of PTEN tumor suppression. Cell. 2008. doi:10.1016/j.cell.2008.04.013.Google Scholar
- Maehama T, Dixon JE. The tumor suppressor, PTEN/MMAC1, dephosphorylates the lipid second messenger, phosphatidylinositol 3,4,5-trisphosphate. J Biol Chem. 1998;273:13375–8.View ArticlePubMedGoogle Scholar
- Song MS, Salmena L, Pandolfi PP. The functions and regulation of the PTEN tumour suppressor. Nat Rev Mol Cell Biol. 2012. doi:10.1038/nrm3330.PubMedGoogle Scholar
- Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010. doi:10.1038/nature09144.PubMedPubMed CentralGoogle Scholar
- Tay Y, Kats L, Salmena L, Weiss D, Tan SM, Ala U, Karreth F, Poliseno L, Provero P, Di Cunto F, Lieberman J, Rigoutsos I, Pandolfi PP. Coding-independent regulation of the tumor suppressor PTEN by competing endogenous mRNAs. Cell. 2011. doi:10.1016/j.cell.2011.09.029.Google Scholar
- Karreth FA, Tay Y, Perna D, Ala U, Tan SM, Rust AG, DeNicola G, Webster KA, Weiss D, Perez-Mancera PA, Krauthammer M, Halaban R, Provero P, Adams DJ, Tuveson DA, Pandolfi PP. In vivo identification of tumor- suppressive PTEN ceRNAs in an oncogenic BRAF-induced mouse model of melanoma. Cell. 2011. doi:10.1016/j.cell.2011.09.032.Google Scholar
- Lee DY, Jeyapalan Z, Fang L, Yang J, Zhang Y, Yee AY, Li M, Du WW, Shatseva T, Yang BB. Expression of versican 3′-untranslated region modulates endogenous microRNA functions. PLoS One. 2010. doi:10.1371/journal.pone.0013599.Google Scholar