The p53-signaling pathway and colorectal cancer: Interactions between downstream p53 target genes and miRNAs.

INTRODUCTION: We examined expression of genes in the p53-signaling pathway. We determine if genes that have significantly different expression in carcinoma tissue compared to normal mucosa also have significantly differentially expressed miRNAs. We utilize a sample of 217 CRC cases. METHODS: We focused on fold change (FC) > 1.50 or <0.67 for genes and miRNAs, that were statistically significant after adjustment for multiple comparisons. We evaluated the linear association between the differential expression of miRNA and mRNA. miRNA:mRNA seed-region matches also were determined. RESULTS: Eleven dysregulated genes were associated with 37 dysregulated miRNAs; all were down-stream from the TP53 gene. MiR-150-5p (HR = 0.82) and miR-196b-5p (HR 0.73) significantly reduced the likelihood of dying from CRC when miRNA expression increased in rectal tumors. CONCLUSIONS: Our data suggest that activation of p53 from cellular stress, could target downstream genes that in turn could influence cell cycle arrest, apoptosis, and angiogenesis through mRNA:miRNA interactions.

MicroRNA-transcription factor interactions and their combined effect on target gene expression in colon cancer cases.

Transcription factors (TFs) and microRNAs (miRNAs) regulate gene expression: TFs by influencing messenger RNA (mRNA) transcription and miRNAs by influencing mRNA translation and transcript degradation. Additionally, miRNAs and TFs alter each other's expression, making it difficult to ascertain the effect either one has on target gene (TG) expression. In this investigation, we use a two-way interaction model with the TF and miRNA as independent variables to investigate whether miRNAs and TFs work together to influence TG expression levels in colon cancer subjects. We used known TF binding sites and validated miRNA targets to determine potential miRNA-TF-TG interactions, restricting interactions to those with a TF previously associated with altered risk of colorectal cancer death. We analyzed interactions using normal colonic mucosa expression as well as differential expression, which is measured as colonic carcinoma expression minus normal colonic mucosa expression. We analyzed 3518 miRNA-TF-TG triplets using normal mucosa expression and 617 triplets using differential expression. Normal colonic RNA-Seq data were available for 168 individuals; of these, 159 also had carcinoma RNA-Seq data. Thirteen unique miRNA-TF-TG interactions, comprising six miRNAs, four TFs, and 11 TGs, were statistically significant after adjustment for multiple comparisons in normal colonic mucosa, and 14 unique miRNA-TF-TG interactions, comprising two miRNAs, two TFs, and 13 TGs, were found for carcinoma-normal differential expression. Our results show that TG expression is influenced by both miRNAs as well as TFs, and the influence of one regulator impacts the effect of the other on the shared TG expression.

Dysregulated genes and miRNAs in the apoptosis pathway in colorectal cancer patients.

Apoptosis is genetically regulated and involves intrinsic and extrinsic pathways. We examined 133 genes within these pathways to identify whether they are expressed differently in colorectal carcinoma (CRC) and normal tissue (N = 217) and if they are associated with similar differential miRNA expression. Gene expression data (RNA-Seq) and miRNA expression data (Agilent Human miRNA Microarray V19.0) were generated. We focused on dysregulated genes with a fold change (FC) of > 1.50 or < 0.67, that were significant after adjustment for multiple comparisons. miRNA:mRNA seed-region matches were determined. Twenty-three genes were significantly downregulated (FC < 0.67) and 18 were significantly upregulated (FC > 1.50). Of these 41 genes, 11 were significantly associated with miRNA differential expression. BIRC5 had the greatest number of miRNA associations (14) and the most miRNAs with a seed-region match (10). Four of these matches, miR-145-5p, miR-150-5p, miR-195-5p, and miR-650, had a negative beta coefficient. CSF2RB was associated with ten total miRNAs (five with a seed-region match, and one miRNA, miR-92a-3p, with a negative beta coefficient). Of the three miRNAs associated with CTSS, miR-20b-5p, and miR-501-3p, had a seed-region match and a negative beta coefficient between miRNA:mRNA pairs. Several miRNAs that were associated with dysregulated gene expression, seed-region matches, and negative beta coefficients also were associated with CRC-specific survival. Our data suggest that miRNAs could influence several apoptosis-related genes. BIRC5, CTSS, and CSF2R all had seed-region matches with miRNAs that would favor apoptosis. Our study identifies several miRNA associated with apoptosis-related genes, that if validated, could be important therapeutic targets.

The PI3K/AKT signaling pathway: Associations of miRNAs with dysregulated gene expression in colorectal cancer.

The PI3K/AKT-signaling pathway is one of the most frequently activated signal-transduction pathways in cancer. We examined how dysregulated gene expression is associated with miRNA expression in this pathway in colorectal cancer (CRC). We used data from 217 CRC cases to evaluate differential pathway gene expression between paired carcinoma and normal mucosa and identify miRNAs that are associated with these genes. Gene expression data from RNA-Seq and miRNA expression data from Agilent Human miRNA Microarray V19.0 were analyzed. We focused on genes most associated with CRC (fold change (FC) of >1.5 or <0.67) that were statistically significant after adjustment for multiple comparisons. Of the 304 genes evaluated, 76 had a FC of <0.67, and 57 had a FC of >1.50; 47 of these genes were associated with miRNA differential expression. There were 145 mRNA:miRNA seed-region matches of which 26 were inversely associated suggesting a greater likelihood of a direct association. Most miRNA:mRNA associations were with factors that stimulated the pathway. For instance, both IL6R and PDGFRA had inverse seed-region matches with seven miRNAs, suggesting that these miRNAs have a direct effect on these genes and may be key elements in activation of the pathway. Other miRNA:mRNA associations with similar impact on the pathway were miR-203a with ITGA4, miR-6071 with ITGAV, and miR-375 with THBS2, all genes involved in extracellular matrix function that activate PI3Ks. Gene expression in the PI3K/Akt-signaling pathway is dysregulated in CRC. MiRNAs were associated with many of these dysregulated genes either directly or in an indirect manner.

Infrequently expressed miRNAs in colorectal cancer tissue and tumor molecular phenotype.

This corrects the article DOI: 10.1038/modpathol.2017.38.

The TGFß-signaling pathway and colorectal cancer: associations between dysregulated genes and miRNAs.

BACKGROUND: The TGFß-signaling pathway plays an important role in the pathogenesis of colorectal cancer (CRC). Loss of function of several genes within this pathway, such as bone morphogenetic proteins (BMPs) have been seen as key events in CRC progression. METHODS: In this study we comprehensively evaluate differential gene expression (RNASeq) of 81 genes in the TGFß-signaling pathway and evaluate how dysregulated genes are associated with miRNA expression (Agilent Human miRNA Microarray V19.0). We utilize paired carcinoma and normal tissue from 217 CRC cases. We evaluate the associations between differentially expressed genes and miRNAs and sex, age, disease stage, and survival months. RESULTS: Thirteen genes were significantly downregulated and 14 were significantly upregulated after considering fold change (FC) of > 1.50 or < 0.67 and multiple comparison adjustment. Bone morphogenetic protein genes BMP5, BMP6, and BMP2 and growth differentiation factor GDF7 were downregulated. BMP4, BMP7, INHBA (Inhibin beta A), TGFBR1, TGFB2, TGIF1, TGIF2, and TFDP1 were upregulated. In general, genes with the greatest dysregulation, such as BMP5 (FC 0.17, BMP6 (FC 0.25), BMP2 (FC 0.32), CDKN2B (FC 0.32), MYC (FC 3.70), BMP7 (FC 4.17), and INHBA (FC 9.34) showed dysregulation in the majority of the population (84.3, 77.4, 81.1, 80.2, 82.0, 51.2, and 75.1% respectively). Four genes, TGFBR2, ID4, ID1, and PITX2, were un-associated or slightly upregulated in microsatellite-stable (MSS) tumors while downregulated in microsatellite-unstable (MSI) tumors. Eight dysregulated genes were associated with miRNA differential expression. E2F5 and THBS1 were associated with one or two miRNAs; RBL1, TGFBR1, TGIF2, and INHBA were associated with seven or more miRNAs with multiple seed-region matches. Evaluation of the joint effects of mRNA:miRNA identified interactions that were stronger in more advanced disease stages and varied by survival months. CONCLUSION: These data support an interaction between miRNAs and genes in the TGFß-signaling pathway in association with CRC risk. These interactions are associated with unique clinical characteristics that may provide targets for further investigations.

Single nucleotide polymorphisms within MicroRNAs, MicroRNA targets, and MicroRNA biogenesis genes and their impact on colorectal cancer survival.

We have shown that single nucleotide polymorphisms (SNPs) in microRNA (miRNA) genes, miRNA target genes, and miRNA biogenesis genes minimally contribute to colon cancer risk. It is possible that these SNPs alter survival. We analyzed 565 SNPs in or adjacent to microRNAs, target genes, or biogenesis genes, using 1,115 cases and 1,173 controls; 837 cases had survival information. We tested SNPs for associations with colorectal cancer (CRC) survival using a Cox proportional hazard model adjusting for age, study center, gender, AJCC disease stage, and MSI tumor status. Multiple comparison adjustments were made using the step-down Bonferroni correction. SNPs associated with survival (Praw  < 0.05) also were assessed with messenger RNA (mRNA). Seven of the 565 SNPs analyzed were associated significantly with CRC survival after adjustment for multiple comparisons. Six of these increased risk of dying, and one, rs12140 (ADAMTS1) decreased risk of dying from CRC (HRR = 0.44, 95% CI (0.24, 0.83; PHolm  = 0.011). Six SNPs altered colon cancer risk and five were associated with altered mRNA expression across genotypes. One SNP, rs2059691 (PRKRA), was associated with increased mRNA expression and worse survival, and one SNP, rs6598964 (LIN28A), decreased risk of developing colon cancer [OR = 0.77 95% CI (0.61, 0.98)] and increased risk of dying from CRC (HRR = 2.26 95% CI (1.52, 3.36). PHolm  = 0.003). The few SNPs associated with CRC survival, colon cancer risk, or with mRNA expression, resided in genes that influence metastasis and angiogenesis. © 2017 Wiley Periodicals, Inc.

The co-regulatory networks of tumor suppressor genes, oncogenes, and miRNAs in colorectal cancer.

Tumor suppressor genes (TSGs) and oncogenes (OG) are involved in carcinogenesis. MiRNAs also contribute to cellular pathways leading to cancer. We use data from 217 colorectal cancer (CRC) cases to evaluate differences in TSGs and OGs expression between paired CRC and normal mucosa and evaluate how TSGs and OGs are associated with miRNAs. Gene expression data from RNA-Seq and miRNA expression data from Agilent Human miRNA Microarray V19.0 were used. We focus on genes most strongly associated with CRC (fold change (FC) of ≥1.5 or ≤0.67) that were statistically significant after adjustment for multiple comparisons. Of the 74 TSGs evaluated, 22 were associated with carcinoma/normal mucosa differential expression. Ten TSGs were up-regulated (FAM123B, RB1, TP53, RUNX1, MSH2, BRCA1, BRCA2, SOX9, NPM1, and RNF43); six TSGs were down-regulated (PAX5, IZKF1, GATA3, PRDM1, TET2, and CYLD); four were associated with MSI tumors (MLH1, PTCH1, and CEBPA down-regulated and MSH6 up-regulated); and two were associated with MSS tumors (PHF6 and ASXL1 up-regulated). Thirteen of these TSGs were associated with 44 miRNAs. Twenty-seven of the 59 OGs evaluated were dysregulated: 14 down-regulated (KLF4, BCL2, SSETBP1, FGFR2, TSHR, MPL, KIT, PDGFRA, GNA11, GATA2, FGFR3, AR, CSF1R, and JAK3), seven up-regulated (DNMT1, EZH2, PTPN11, SKP2, CCND1, MET, and MYC); three down-regulated for MSI (FLT3, CARD11, and ALK); two up-regulated for MSI (IDH2 and HRAS); and one up-regulated with MSS tumors (CTNNB1). These findings suggest possible co-regulatory function between TSGs, OGs, and miRNAs, involving both direct and indirect associations that operate through feedback and feedforward loops.

The miRNA landscape of colorectal polyps.

The genomic landscape of adenomas and polyps may help define disease pathways. Expression of miRNAs in adenomas and polyps may importantly contribute to these pathways. We evaluated miRNA expression in 293 polyp-normal colorectal mucosa pairs. Polyps were classified as either adenomatous polyp (AD), hyperplastic polyp (HP), or sessile serrated polyp (SSP). We compared these miRNA expression profiles in polyps to miRNA expression in microsatellite unstable (MSI) and stable (MSS) tumors. A False Discovery Rate of 0.05 based on Benjamini and Hochberg was used to adjust for multiple comparisons. There were 70 miRNAs with differential expression by polyp type with a fold change <0.75 or >1.34 after adjusting for multiple comparisons. The major differences in miRNA expression were observed between AD and SSP and AD and HP, with few differences in expression noted for SSP and HP. AD polyps were more likely to be upregulated from normal colonic mucosa, while SSP and HP were more likely to be downregulated from normal colonic mucosa. MiRNA expression in the SSP and HP tumors almost uniformly go in opposite directions from the MSS tumor miRNA expression and was mixed with MSI tumors. We conclude that different types of polyps have unique miRNA expression profiles. © 2017 Wiley Periodicals, Inc.

Transcription factor-microRNA associations and their impact on colorectal cancer survival.

MicroRNAs (miRNAs) and Transcription Factors (TFs) both influence messenger RNA (mRNA) expression, disrupting biological pathways involved in carcinogenesis and prognosis. As many miRNAs target multiple mRNAs, thus influencing a multitude of biological pathways, deciphering which miRNAs are important for cancer development and survival is difficult. In this study, we (i) determine associations between TF and survival (N = 168 colon cancer cases); (ii) identify miRNAs associated with TFs related to survival; and (iii) determine if factors derived from TF-specific miRNA principal component analysis (PCA) influence survival. Cox Proportional hazard models were run for each PCA factor to determine Hazard Ratios (HR) and 95% Confidence Intervals (CI) adjusting for age, center, and AJCC stage. Thirty TFs improved survival when differential expression increased; 27 of these were associated significantly with normal colonic mucosa expression of 65 unique miRNAs when an FDR q-value of <0.05 was applied. Five factors, comprising 21 miRNAs, altered survival in rectal cancer subjects; four of these five factors improved survival and one factor reduced survival. One factor comprising four miRNAs reduced survival in colon cancer subjects. In summary, our data suggest that expression of TFs and their related miRNAs influence survival after diagnosis with colorectal cancer.

Alterations in microRNA expression associated with alcohol consumption in rectal cancer subjects.

PURPOSE: Alcohol consumption has been purported to influence many diseases. MicroRNAs (miRNAs) may be influenced by compounds found in alcohol. In this investigation, we test the hypothesis that total alcohol, beer, wine, and hard liquor influence miRNA expression. METHODS: We studied 1447 colorectal (CR) cancer cases with normal CR mucosa and carcinoma miRNA expression data along with alcohol consumption data. We analyzed long-term and long-term and current (LTC) alcohol use for beer, liquor, and wine with miRNA expression between paired carcinoma and normal colon and rectal tissues, adjusting for multiple comparisons using the positive false discovery rate q-value. MiRNAs associated significantly with alcohol were examined with all-cause mortality (ACM). MiRNAs associated significantly with ACM were examined with RNA-Seq data. RESULTS: Expression of 84 miRNAs was associated significantly with LTC wine use in normal rectal mucosa. Higher expression of two of these miRNAs significantly worsened ACM: hsa-miR-210 (Hazard Ratio [HR] 1.12, 95% CI (1.03, 1.21), p-value = 0.004), and hsa-miR-92a-1-5p (HR 1.20, 95% CI (1.04, 1.38), p-value = 0.013). These miRNAs were downregulated across levels of LTC wine consumption. CONCLUSIONS: Our results suggest that wine influences miRNA expression in rectal cancer, supporting the hypothesis that components in alcohol influence miRNA expression.

Infrequently expressed miRNAs in colorectal cancer tissue and tumor molecular phenotype.

We have previously shown that commonly expressed miRNAs influenced tumor molecular phenotype in colorectal cancer. We hypothesize that infrequently expressed miRNAs, when showing higher levels of expression, help to define tumor molecular phenotype. In this study, we examine 304 miRNAs expressed in at least 30 individuals, but in <50% of the population and with a mean level of expression above 1.0 relative florescent unit. We examine associations in 1893 individuals who have the tumor molecular phenotype data as well as miRNA expression levels for both carcinoma and normal colorectal tissue. We compare miRNAs uniquely associated with tumor molecular phenotype to the RNAseq data to identify genes associated with these miRNAs. This information is used to further identify unique pathways associated with tumor molecular phenotypes of TP53-mutated, KRAS-mutated, CpG island methylator phenotype and microsatellite instability tumors. Thirty-seven miRNAs were uniquely associated with TP53-mutated tumors; 30 of these miRNAs had higher level of expression in TP53-mutated tumors, while seven had lower levels of expression. Of the 34 miRNAs associated with CpG island methylator phenotype-high tumors, 16 were more likely to have a CpG island methylator phenotype-high tumor and 19 were less likely to be CpG island methylator phenotype-high. For microsatellite instability, 13 of the 22 infrequently expressed miRNAs were significantly less likely to be expressed in microsatellite unstable tumors. KRAS-mutated tumors were not associated with any miRNAs after adjustment for multiple comparisons. Of the dysregulated miRNAs, 17 were more likely to be TP53-mutated tumors while simultaneously being less likely to be CpG island methylator phenotype-high and/or microsatellite instability tumors. Genes regulated by these miRNAs were involved in numerous functions and pathways that influence cancer risk and progression. In summary, some infrequently expressed miRNAs, when expressed at higher levels, appear to have significant biological meaning in terms of tumor molecular phenotype and gene expression profiles.

MicroRNA profiles in colorectal carcinomas, adenomas and normal colonic mucosa: variations in miRNA expression and disease progression.

MiRNAs are small, non-protein-coding RNA molecules that regulate gene expression either by post-transcriptionally suppressing mRNA translation or by mRNA degradation. We examine differentially expressed miRNAs in colorectal carcinomas, adenomas and normal colonic mucosa. Data come from population-based studies of colorectal cancer conducted in Utah and the Kaiser Permanente Medical Care Program. A total of 1893 carcinoma/normal-paired samples and 290 adenoma tissue samples were run on the Agilent Human miRNA Microarray V19.0 which contained 2006 miRNAs. We tested for significant differences in miRNA expression between paired carcinoma/adenoma/normal colonic tissue samples. Fewer than 600 miRNAs were expressed in >80% of people for colonic tissue; of these 86.5% were statistically differentially expressed between carcinoma and normal colonic mucosa using a false discovery rate of 0.05. Roughly half of these differentially expressed miRNAs showed a progression in levels of expression from normal to adenoma to carcinoma tissue. Other miRNAs appeared to be altered at the normal to adenoma stage, while others were only altered at the adenoma to carcinoma stage or only at the normal to carcinoma stage. Evaluation of the Agilent platform showed a high degree of repeatability (r = 0.98) and reasonable agreement with the NanoString platform. Our data suggest that miRNAs are highly dysregulated in colorectal tissue among individuals with colorectal cancer; the pattern of disruption varies by miRNA as tissue progresses from normal to adenoma to carcinoma.

Dietary intake alters gene expression in colon tissue: possible underlying mechanism for the influence of diet on disease.

BACKGROUND: Although the association between diet and disease is well documented, the biologic mechanisms involved have not been entirely elucidated. In this study, we evaluate how dietary intake influences gene expression to better understand the underlying mechanisms through which diet operates. METHODS: We used data from 144 individuals who had comprehensive dietary intake and gene expression data from RNAseq using normal colonic mucosa. Using the DESeq2 statistical package, we identified genes that showed statistically significant differences in expression between individuals in high-intake and low-intake categories for several dietary variables of interest adjusting for age and sex. We examined total calories, total fats, vegetable protein, animal protein, carbohydrates, trans-fatty acids, mutagen index, red meat, processed meat, whole grains, vegetables, fruits, fiber, folate, dairy products, calcium, and prudent and western dietary patterns. RESULTS: Using a false discovery rate of less than 0.1, meat-related foods were statistically associated with 68 dysregulated genes, calcium with three dysregulated genes, folate with four dysregulated genes, and nonmeat-related foods with 65 dysregulated genes. With a more stringent false discovery rate of less than 0.05, there were nine meat-related dysregulated genes and 23 nonmeat-related genes. Ingenuity pathway analysis identified three major networks among genes identified as dysregulated with respect to meat-related dietary variables and three networks among genes identified as dysregulated with respect to nonmeat-related variables. The top networks (Ingenuity Pathway Analysis network score >30) associated with meat-related genes were (i) cancer, organismal injury, and abnormalities, tumor morphology, and (ii) cellular function and maintenance, cellular movement, cell death, and survival. Among genes related to nonmeat consumption variables, the top networks were (i) hematological system development and function, nervous system development and function, tissue morphology and (ii) connective tissue disorders, organismal injury, and abnormalities. CONCLUSION: Several dietary factors were associated with gene expression in our data. These findings provide insight into the possible mechanisms by which diet may influence disease processes.

Colorectal tumor molecular phenotype and miRNA: expression profiles and prognosis.

MiRNAs regulate gene expression by post-transcriptionally suppressing mRNA translation or by causing mRNA degradation. It has been proposed that unique miRNAs influence specific tumor molecular phenotype. In this paper, we test the hypotheses that miRNA expression differs by tumor molecular phenotype and that those differences may influence prognosis. Data come from population-based studies of colorectal cancer conducted in Utah and the Northern California Kaiser Permanente Medical Care Program. A total of 1893 carcinoma samples were run on the Agilent Human miRNA Microarray V19.0 containing 2006 miRNAs. We assessed differences in miRNA expression between TP53-mutated and non-mutated, KRAS-mutated and non-mutated, BRAF-mutated and non-mutated, CpG island methylator phenotype (CIMP) high and CIMP low, and microsatellite instability (MSI) and microsatellite stable (MSS) colon and rectal tumors. Using a Cox proportional hazard model we evaluated if those miRNAs differentially expressed by tumor phenotype influenced survival after adjusting for age, sex, and AJCC stage. There were 22 differentially expressed miRNAs for TP53-mutated colon tumors and 5 for TP53-mutated rectal tumors with a fold change of >1.49 (or <0.67). Additionally, 13 miRNAS were differentially expressed for KRAS-mutated rectal tumors, 8 differentially expressed miRNAs for colon CIMP high tumors, and 2 differentially expressed miRNAs for BRAF-mutated colon tumors. The majority of differentially expressed miRNAS were observed between MSI and MSS tumors (94 differentially expressed miRNAs for colon; 41 differentially expressed miRNAs for rectal tumors). Of these miRNAs differentially expressed between MSI and MSS tumors, the majority were downregulated. Ten of the differentially expressed miRNAs were associated with survival; after adjustment for MSI status, five miRNAS, miR-196b-5p, miR-31-5p, miR-99b-5p, miR-636, and miR-192-3p, were significantly associated with survival. In summary, it appears that the majority of miRNAs that are differentially expressed by tumor molecular phenotype are MSI tumors. However, these miRNAs appear to have minimal effect on prognosis.

Gene expression in colon cancer: A focus on tumor site and molecular phenotype.

Hundreds to thousands of genes are differentially expressed in tumors when compared to nontumor colonic tissue samples. We evaluated gene expression patterns to better understand differences in colon cancer by tumor site and tumor molecular phenotype. We analyzed RNA-seq data from tumor/normal paired samples from 175 colon cancer patients. We implemented a cross validation strategy with nonparametric tests to identify genes which displayed varying expression characteristics related to paired tumor/nontumor tissue across proximal and distal colon sites and by tumor molecular phenotypes, that is, TP53, KRAS, CpG Island Methylator Phenotype (CIMP), and microsatellite instability (MSI). We used Ingenuity Pathway Analysis (IPA) to determine networks associated with deregulated genes in our data. Genes showed significant differences in expression characteristics at the 0.01 level in both validation groups between tumor subsite (116 genes), CIMP high versus CIMP low (79 genes), MSI versus microsatellite stable (MSS) (49 genes), TP53-mutated versus not mutated (17genes), and KRAS-mutated versus not mutated (1 gene). Deregulated genes for CIMP high and MSI tumors were often down-regulated. In contrast to CIMP high and MSI tumors, genes that were deregulated in TP53 were likely to be up-regulated. ERK1, WNT, growth factors and inflammation-related factors were focal points of both CIMP and MSI IPA networks. The MUC family of genes was up-regulated MSI networks. Numerous genes showed differences in expression between proximal and distal tumors, nontumor proximal and distal tissue, and tumor molecular phenotype. Deregulated mucin genes appear to play an important role in MSI tumors.

An evaluation and replication of miRNAs with disease stage and colorectal cancer-specific mortality.

MicroRNAs (miRNAs) have been implicated in colorectal cancer (CRC) development and associated with prognostic indicators such as disease stage and survival. Prognostic associations are often based on few individuals and imprecise. In this study, we utilize population-based data from 1,141 CRC cases to replicate previously reported associations between 121 miRNAs and disease stage and survival. The Agilent Human miRNA Microarray V19.0 was used to generate miRNA data following a stringent quality control protocol. Assessment of survival was done using Cox Proportional Hazard models adjusting for age, disease stage and tumor molecular phenotype. Five miRNAs were associated with more advanced disease stage; hsa-miR-145-5p and hsa-miR-31-5p showed increased expression with more advanced tumor stage, while hsa-miR-200b-3p, hsa-miR-215 and hsa-miR-451a had decreased expression with more advanced tumors. Thirteen miRNAs were associated with CRC mortality among individuals diagnosed with colon cancer while 14 were associated with CRC mortality after a diagnosis with rectal cancer. Strongest associations were observed for those miRNAs that were expressed in a small subset of tumors. Most notable associations were for hsa-miR-145-3p [hazard ratio (HR) 2.94, 95% confidence interval (CI) 1.54, 5.61], and hsa-miR-9-3p (HR 10.28, 95% CI 1.31, 80.84) with colon cancer and hsa-miR-335-5p (HR 0.17, 95% CI 0.05, 0.54) for rectal cancer. hsa-miR-374a-5p, hsa-miR-570-3p and hsa-miR-18a-5p significantly reduced the hazard of dying for all cases, regardless of tumor site. Our findings illustrate the need for a large sample to evaluate the association of miRNAs with survival and disease stage in order to determine associations by tumor site.

Improved survival among colon cancer patients with increased differentially expressed pathways.

BACKGROUND: Studies of colorectal cancer (CRC) have shown that hundreds to thousands of genes are differentially expressed in tumors when compared to normal tissue samples. In this study, we evaluate how genes that are differentially expressed in colon versus normal tissue influence survival. METHODS: We performed RNA-seq on tumor/normal paired samples from 175 colon cancer patients. We implemented a cross validation strategy to determine genes that were significantly differentially expressed between tumor and normal samples. Differentially expressed genes were evaluated with Ingenuity Pathway Analysis to identify key pathways that were de-regulated. A summary differential pathway expression score (DPES) was developed to summarize hazard of dying while adjusting for age, American Joint Committee on Cancer (AJCC) stage, sex, and tumor molecular phenotype, i.e., MSI, TP53, KRAS, and CIMP. RESULTS: A total of 1,138 genes were up-regulated and 695 were down-regulated. These de-regulated genes were enriched for 19 Ingenuity Canonical Pathways, with the most significant pathways involving cell signaling and growth. Of the enriched pathways, 16 were significantly associated with CRC-specific mortality, including 1 metabolic pathway and 15 signaling pathways. In all instances, having a higher DPES (i.e., more de-regulated genes) was associated with better survival. Further assessment showed that individuals diagnosed at AJCC Stage 1 had more de-regulated genes than individuals diagnosed at AJCC Stage 4. CONCLUSIONS: Our data suggest that having more de-regulated pathways is associated with a good prognosis and may be a reaction to key events that are disabling to tumor progression. Please see related article: http://dx.doi.org/10.1186/s12916-015-0307-6 .

The functional role of miRNAs in colorectal cancer: insights from a large population-based study.

Identification of causal microRNAs (miRNAs) in colorectal cancer (CRC) is elusive, due to our lack of understanding of how specific miRNAs affect biological pathways and outcomes. An miRNA can regulate many mRNAs and an mRNA can be associated with many miRNAs; appreciation of these complex networks in which miRNAs operate is necessary to transition from identifying dysregulated miRNAs to identifying individual miRNAs or groups of miRNAs that are suitable for therapeutic purposes. The aim of the paper is to compile results from a population-based study (n = 1,954 cases with matched carcinoma/normal tissue) of miRNAs in CRC. The information gained allows for cohesive and comprehensive insight into miRNAs and CRC in terms of function and impact. Comparison of miRNA expression with mRNA expression from nine signaling pathways in carcinogenic processes allowed us to identify miRNA targets within a biological context. MiRNAs that directly influence mRNA expression may be effective biomarkers or therapeutic targets.

The MAPK-Signaling Pathway in Colorectal Cancer: Dysregulated Genes and Their Association With MicroRNAs.

Mitogen-activated protein kinase (MAPK) pathways regulate many cellular functions including cell proliferation and apoptosis. We examined associations of differential gene and microRNA (miRNA) expression between carcinoma and paired normal mucosa for 241 genes in the KEGG-identified MAPK-signaling pathway among 217 colorectal cancer (CRC) cases. Gene expression data (RNA-Seq) and miRNA expression data (Agilent Human miRNA Microarray V19.0; Agilent Technologies Inc., Santa Clara, CA, USA) were analyzed. We first identified genes most strongly associated with CRC using a fold change (FC) of >1.50 or <0.67) that were statistically significant after adjustment for multiple comparisons. We then determined miRNAs associated with dysregulated genes and through miRNA:mRNA (messenger RNA) seed region matches discerned genes with a greater likelihood of having a direct biological association. Ninety-nine genes had a meaningful FC for all CRC, microsatellite unstable-specific tumors, or microsatellite stable-specific tumors. Thirteen dysregulated genes were associated with miRNAs, totaling 68 miRNA:mRNA associations. Thirteen of the miRNA:mRNA associations had seed region matches where the differential expression between the miRNA and mRNA was inversely related suggesting a direct association as a result of their binding. Several direct associations, upstream of ERK1/ERK2, JNK, and p38, were found for PDGFRA with 7 miRNAs; RASGRP3 and PRKCB with miR-203a; and TGFBR1 with miR-6071 and miR-2117. Other associations between miRNAs and mRNAs are most likely indirect, resulting from feedback and feed forward loops. Our results suggest that miRNAs may alter MAPK signaling through direct binding with key genes in this pathway. We encourage others to validate results in targeted CRC experiments that can help solidify important therapeutic targets.