Comprehensive analysis of differential expression profiles of mRNAs and lncRNAs and identification of a 14-lncRNA prognostic signature for patients with colon adenocarcinoma.

The objective of this study was to identify potentially significant genes and long non-coding RNAs (lncRNAs) in colon cancer for a panel of lncRNA signatures that could be used as prognostic markers for colon adenocarcinoma (COAD) based on the data from The Cancer Genome Atlas (TCGA). RNA-seq V2 exon data of COAD were downloaded from the TCGA data portal for 285 tumor samples and 41 normal tissue samples adjacent to tumors. Differentially expressed mRNAs and lncRNAs were identified. A functional enrichment analysis of differentially expressed mRNAs was performed, followed by protein-protein interaction (PPI) network construction and significant module selection. Additionally, the regulatory relationships in differentially expressed mRNAs and lncRNAs were assessed, and an lncRNA-lncRNA co-regulation and functional synergistic analysis were performed. Furthermore, the risk score model and Cox regression analysis based on the expression levels of lncRNAs were used to develop a prognostic lncRNA signature. A total of 976 differentially expressed mRNAs and 169 differentially expressed lncRNAs were identified. MDFI and MEOX2 were the PPI network hubs. We found these lncRNAs to be mainly involved in vascular smooth muscle contraction and the cGMP-PKG signaling pathway. Several lncRNA-lncRNA pairs had co-regulatory relationships or functional synergistic effects, including BVES-AS1/MYLK-AS1, ADAMTS9-AS1/MYLK-AS1 and FENDRR/MYLK-AS1. The differential expression profile analysis of four candidate lncRNAs (MYLK-AS1, BVES-AS1, ADAMTS9-AS1, and FENDRR) in COAD tumors were confirmed by reverse transcription-quantitative PCR. Moreover, this study identified a 14-lncRNA signature that could predict the survival for COAD patients.

[1]  S. Hou,et al.  Differentially expressed long non-coding RNAs and the prognostic potential in colorectal cancer. , 2016, Neoplasma.

[2]  Michael Kohl,et al.  Cytoscape: software for visualization and analysis of biological networks. , 2011, Methods in molecular biology.

[3]  Gary D. Bader,et al.  An automated method for finding molecular complexes in large protein interaction networks , 2003, BMC Bioinformatics.

[4]  Xi Chen,et al.  Sulindac Selectively Inhibits Colon Tumor Cell Growth by Activating the cGMP/PKG Pathway to Suppress Wnt/β-Catenin Signaling , 2013, Molecular Cancer Therapeutics.

[5]  R. Wolff,et al.  MicroRNA Seed Region Length Impact on Target Messenger RNA Expression and Survival in Colorectal Cancer , 2016, PloS one.

[6]  David Bryant,et al.  DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists , 2007, Nucleic Acids Res..

[7]  P. Pauwels,et al.  Prognostic Role of Vasculogenic Mimicry in Colorectal Cancer , 2009, Diseases of the colon and rectum.

[8]  Jordan Anaya,et al.  The lncRNA DRAIC/PCAT29 Locus Constitutes a Tumor-Suppressive Nexus , 2015, Molecular Cancer Research.

[9]  M. Esteller,et al.  lncRNAs and microRNAs with a role in cancer development. , 2016, Biochimica et biophysica acta.

[10]  A. Rabson,et al.  MEOX2 regulates nuclear factor-kappaB activity in vascular endothelial cells through interactions with p65 and IkappaBbeta. , 2010, Cardiovascular research.

[11]  J. Tainer,et al.  Noncoding RNA joins Ku and DNA-PKcs for DNA-break resistance in breast cancer , 2016, Nature Structural &Molecular Biology.

[12]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[13]  Fudong Yu,et al.  LncRNA‐ATB mediated E‐cadherin repression promotes the progression of colon cancer and predicts poor prognosis , 2016, Journal of gastroenterology and hepatology.

[14]  Chyke A Doubeni,et al.  Adenoma detection rate and risk of colorectal cancer and death. , 2014, The New England journal of medicine.

[15]  Jianzhi Zhang,et al.  Why Do Hubs Tend to Be Essential in Protein Networks? , 2006, PLoS genetics.

[16]  Didier Y. R. Stainier,et al.  Molecular control of endothelial cell behaviour during blood vessel morphogenesis , 2011, Nature Reviews Molecular Cell Biology.

[17]  J. Kim,et al.  Polypyrimidine tract-binding protein 1-mediated down-regulation of ATG10 facilitates metastasis of colorectal cancer cells. , 2017, Cancer letters.

[18]  H. Ishwaran.  Variable importance in binary regression trees and forests , 2007, 0711.2434.

[19]  A. Hochberg,et al.  The H19 Long non-coding RNA in cancer initiation, progression and metastasis – a proposed unifying theory , 2015, Molecular Cancer.

[20]  Sandrine Imbeaud,et al.  DNA methylation‐based prognosis and epidrivers in hepatocellular carcinoma , 2015, Hepatology.

[21]  Jan P Stegemann,et al.  Biomechanics and Mechanotransduction in Cells and Tissues Mechanical , biochemical , and extracellular matrix effects on vascular smooth muscle cell phenotype , 2005 .

[22]  Bronwen L. Aken,et al.  GENCODE: The reference human genome annotation for The ENCODE Project , 2012, Genome research.

[23]  Ming Sun,et al.  Decreased expression of the long non-coding RNA FENDRR is associated with poor prognosis in gastric cancer and FENDRR regulates gastric cancer cell metastasis by affecting fibronectin1 expression , 2014, Journal of Hematology & Oncology.

[24]  Zhonghua Ma,et al.  A novel lncRNA, LL22NC03-N64E9.1, represses KLF2 transcription through binding with EZH2 in colorectal cancer , 2017, Oncotarget.

[25]  Chaofeng Liu,et al.  Adjusted Kaplan–Meier estimator and log‐rank test with inverse probability of treatment weighting for survival data , 2005, Statistics in medicine.

[26]  Howard Y. Chang,et al.  Unique features of long non-coding RNA biogenesis and function , 2015, Nature Reviews Genetics.

[27]  Limei Ma,et al.  C-Myc-activated long noncoding RNA CCAT1 promotes colon cancer cell proliferation and invasion , 2014, Tumor Biology.

[28]  Nathalia Amado,et al.  Isoquercitrin Suppresses Colon Cancer Cell Growth in Vitro by Targeting the Wnt/β-Catenin Signaling Pathway* , 2014, The Journal of Biological Chemistry.

[29]  W. Leung,et al.  Increasing incidence of colorectal cancer in Asia: implications for screening. , 2005, The Lancet. Oncology.

[30]  Charlotte Soneson,et al.  A comparison of methods for differential expression analysis of RNA-seq data , 2013, BMC Bioinformatics.

[31]  Hemant Ishwaran,et al.  Random Survival Forests , 2008, Wiley StatsRef: Statistics Reference Online.

[32]  A. Jemal,et al.  Cancer statistics, 2016 , 2016, CA: a cancer journal for clinicians.

[33]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[34]  G. Bellone,et al.  Low-doses of sequential-kinetic-activated interferon-γ enhance the ex vivo cytotoxicity of peripheral blood natural killer cells from patients with early-stage colorectal cancer. A preliminary study. , 2014, International immunopharmacology.

[35]  M. Kretz,et al.  Long non-coding RNA ZFAS1 interacts with CDK1 and is involved in p53-dependent cell cycle control and apoptosis in colorectal cancer , 2015, OncoTarget.

[36]  S. Babykutty,et al.  Insidious role of nitric oxide in migration/invasion of colon cancer cells by upregulating MMP-2/9 via activation of cGMP-PKG-ERK signaling pathways , 2012, Clinical & Experimental Metastasis.

[37]  Charity W. Law,et al.  voom: precision weights unlock linear model analysis tools for RNA-seq read counts , 2014, Genome Biology.

[38]  Maoni Guo,et al.  A potential prognostic long non-coding RNA signature to predict metastasis-free survival of breast cancer patients , 2015, Scientific Reports.

[39]  Matthew E. Ritchie,et al.  limma powers differential expression analyses for RNA-sequencing and microarray studies , 2015, Nucleic acids research.

[40]  Fang Zhou,et al.  LncRNA profile study reveals a three-lncRNA signature associated with the survival of patients with oesophageal squamous cell carcinoma , 2014, Gut.

[41]  Susumu Goto,et al.  The KEGG resource for deciphering the genome , 2004, Nucleic Acids Res..

[42]  Yi Pan,et al.  CytoNCA: A cytoscape plugin for centrality analysis and evaluation of protein interaction networks , 2015, Biosyst..

[43]  S. Tapscott,et al.  MyoD and the transcriptional control of myogenesis. , 2005, Seminars in cell & developmental biology.

[44]  Jaw-Yuan Wang,et al.  Nesfatin-1/Nucleobindin-2 enhances cell migration, invasion, and epithelial-mesenchymal transition via LKB1/AMPK/TORC1/ZEB1 pathways in colon cancer , 2016, Oncotarget.

[45]  A Riedl,et al.  Stromal-derived IGF2 promotes colon cancer progression via paracrine and autocrine mechanisms , 2017, Oncogene.

[46]  A. Hargrove,et al.  Biochemical Methods To Investigate lncRNA and the Influence of lncRNA:Protein Complexes on Chromatin. , 2016, Biochemistry.

[47]  N. Raab-Traub,et al.  I-mfa Domain Proteins Interact with Axin and Affect Its Regulation of the Wnt and c-Jun N-Terminal Kinase Signaling Pathways , 2002, Molecular and Cellular Biology.

[48]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[49]  S. Duan,et al.  DNA methylation of CMTM3, SSTR2, and MDFI genes in colorectal cancer. , 2017, Gene.

[50]  Yanying Sun,et al.  Comprehensive analysis of lncRNA expression profiles reveals a novel lncRNA signature to discriminate nonequivalent outcomes in patients with ovarian cancer , 2016, Oncotarget.

[51]  Guangchuang Yu,et al.  clusterProfiler: an R package for comparing biological themes among gene clusters. , 2012, Omics : a journal of integrative biology.

[52]  K. Mimori,et al.  A long noncoding RNA, lncRNA-ATB, is involved in the progression and prognosis of colorectal cancer. , 2015, Anticancer research.