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 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.

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.

miRNA involvement in cell cycle regulation in colorectal cancer cases.

Uncontrolled cell replication is a key component of carcinogenesis. MicroRNAs (miRNAs) regulate genes involved in checkpoints, DNA repair, and genes encoding for key proteins regulating the cell cycle. We investigated how miRNAs and mRNAs in colorectal cancer subjects interact to regulate the cell cycle. Using RNA-Seq data from 217 individuals, we analyzed differential expression (carcinoma minus normal mucosa) of 123 genes within the cell cycle pathway with differential miRNA expression, adjusting for age and sex. Multiple comparison adjustments for gene/miRNA associations were made at the gene level using an FDR <0.05. Differentially expressed miRNAs and mRNAs were tested for associations with colorectal cancer survival. MRNA and miRNA sequences were compared to identify seed region matches to support biological interpretation of the observed associations. Sixty-seven mRNAs were dysregulated with a fold change (FC) <0.67 or >1.50. Thirty-two mRNAs were associated with 48 miRNAs; 102 of 290 total associations had identified seed matches; of these, ten had negative beta coefficients. Hsa-miR-15a-5p and hsa-miR-20b-5p were associated with colorectal cancer survival with an FDR <0.05 (HR 0.86 95% CI 0.79, 0.94; HR 0.83 95% CI 0.75, 0.91 respectively). Our findings suggest that miRNAs impact mRNA translation at multiple levels within the cell cycle.

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.

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.

Expression of Wnt-signaling pathway genes and their associations with miRNAs in colorectal cancer.

The Wnt-signaling pathway functions in regulating cell growth and thus is involved in the carcinogenic process of several cancers, including colorectal cancer. We tested the hypothesis that multiple genes in this signaling pathway are dysregulated and that miRNAs are associated with these dysregulated genes. We used data from 217 colorectal cancer (CRC) cases to evaluate differences in Wnt-signaling pathway gene expression between paired CRC 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 strongly associated with CRC (fold change (FC) of >1.5 or <0.67) and that were statistically significant after adjustment for multiple comparisons. Of the 138 Wnt-signaling pathway genes examined, 27 were significantly down-regulated (FC<0.67) and 32 genes were significantly up-regulated (FC>1.50) after adjusting for multiple comparisons. Thirteen of the 66 Wnt-signaling genes that were differentially expressed in CRC tumors were associated with differential expression of miRNAs. A total of 93 miRNA:mRNA associations were detected for these 13 genes. Of these 93 associations, 36 miRNA seed-region matches were observed, suggesting that miRNAs have both direct and indirect effects on Wnt-signaling pathway genes. In summary, our data supports the hypothesis that the Wnt-signaling pathway is dysregulated in CRC and suggest that miRNAs may importantly influence that dysregulation.

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

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

The NF-κB signalling pathway in colorectal cancer: associations between dysregulated gene and miRNA expression.

BACKGROUND: The nuclear factor-kappa B (NF-κB) signalling pathway is a regulator of immune response and inflammation that has been implicated in the carcinogenic process. We examined differentially expressed genes in this pathway and miRNAs to determine associations with colorectal cancer (CRC). METHODS: We used data from 217 CRC cases to evaluate differences in NF-κB signalling 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 analysed. We evaluated genes most strongly associated and differentially expressed (fold change (FC) of > 1.5 or < 0.67) that were statistically significant after adjustment for multiple comparisons. RESULTS: Of the 92 genes evaluated, 22 were significantly downregulated and nine genes were significantly upregulated in all tumours. Two additional genes (CD14 and CSNK2A1) were dysregulated in MSS tumours and two genes (CARD11 and VCAM1) were downregulated and six genes were upregulated (LYN, TICAM2, ICAM1, IL1B, CCL4 and PTGS2) in MSI tumours. Sixteen of the 21 dysregulated genes were associated with 40 miRNAs. There were 76 miRNA:mRNA associations of which 38 had seed-region matches. Genes were associated with multiple miRNAs, with TNFSRF11A (RANK) being associated with 15 miRNAs. Likewise several miRNAs were associated with multiple genes (miR-150-5p with eight genes, miR-195-5p with four genes, miR-203a with five genes, miR-20b-5p with four genes, miR-650 with six genes and miR-92a-3p with five genes). CONCLUSIONS: Focusing on the genes and their associated miRNAs within the entire signalling pathway provides a comprehensive understanding of this complex pathway as it relates to CRC and offers insight into potential therapeutic agents.

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.

MicroRNA-messenger RNA interactions involving JAK-STAT signaling genes in colorectal cancer.

JAK-STAT signaling influences many downstream processes that, unchecked, contribute to carcinogenesis and metastasis. MicroRNAs (miRNAs) are hypothesized as a mechanism to prevent uncontrolled growth from continuous JAK-STAT activation. We investigated differential expression between paired carcinoma and normal colorectal mucosa of messenger RNAs (mRNAs) and miRNAs using RNA-Seq and Agilent Human miRNA Microarray V19.0 data, respectively, using a negative binomial mixed effects model to test 122 JAK-STAT-signaling genes in 217 colorectal cancer (CRC) cases. Overall, 42 mRNAs were differentially expressed with a fold change of >1.50 or <0.67, remaining significant with a false discovery rate of < 0.05; four were dysregulated in microsatellite stable (MSS) tumors, eight were for microsatellite unstable (MSI)-specific tumors. Of these 54 mRNAs, 17 were associated with differential expression of 46 miRNAs, comprising 116 interactions: 16 were significant overall, one for MSS tumors only. Twenty of the 29 interactions with negative beta coefficients involved miRNA seed sequence matches with mRNAs, supporting miRNA-mediated mRNA repression; 17 of these mRNAs encode for receptor molecules. Receptor molecule degradation is an established JAK-STAT signaling control mechanism; our results suggest that miRNAs facilitate this process. Interactions involving positive beta coefficients may illustrate downstream effects of disrupted STAT activity, and subsequent miRNA upregulation.

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 influence survival after diagnosis with colorectal cancer.

Half of miRNAs expressed in colorectal tissue are expressed < 50% of the population. Many infrequently expressed miRNAs have low levels of expression. We hypothesize that less frequently expressed miRNAs, when expressed at higher levels, influence both disease stage and survival after diagnosis with colorectal cancer (CRC); low levels of expression may be background noise. We examine 304 infrequently expressed miRNAs in 1893 population-based cases of CRC with paired carcinoma and normal mucosa miRNA profiles. We evaluate miRNAs with disease stage and survival after adjusting for age, study center, sex, MSI status, and AJCC stage. These miRNAs were further evaluated with RNA-Seq data to identify miRNA::mRNA associations that may provide insight into the functionality of miRNAs. Eleven miRNAs were associated with advanced disease stage among colon cancer patients (Q value = 0.10). Eight infrequently expressed miRNAs influenced survival if highly expressed in overall CRC. Of these, five increased likelihood of dying if they were highly expressed, i.e. miR-124-3p, miR-143-5p, miR-145-3p, miR31-5p, and miR-99b-5p, while three were associated with better survival if highly expressed, i.e. miR-362-5p, miR-374a-5p, and miR-590-5p. Thirteen miRNAs infrequently expressed in colon-specific carcinoma tissue were associated with CRC survival if highly expressed. Evaluation of miRNAs::mRNA associations showed that mRNA expression influenced by infrequently expressed miRNA contributed to networks and pathways shown to influence disease progression and prognosis. Our large study enabled us to examine the implications of infrequently expressed miRNAs after removal of background noise. These results require replication in other studies. Confirmation of our findings in other studies could lead to important markers for prognosis.

An Assessment of Database-Validated microRNA Target Genes in Normal Colonic Mucosa: Implications for Pathway Analysis.

BACKGROUND: Determination of functional pathways regulated by microRNAs (miRNAs), while an essential step in developing therapeutics, is challenging. Some miRNAs have been studied extensively; others have limited information. In this study, we focus on 254 miRNAs previously identified as being associated with colorectal cancer and their database-identified validated target genes. METHODS: We use RNA-Seq data to evaluate messenger RNA (mRNA) expression for 157 subjects who also had miRNA expression data. In the replication phase of the study, we replicated associations between 254 miRNAs associated with colorectal cancer and mRNA expression of database-identified target genes in normal colonic mucosa. In the discovery phase of the study, we evaluated expression of 18 miR-NAs (those with 20 or fewer database-identified target genes along with miR-21-5p, miR-215-5p, and miR-124-3p which have more than 500 database-identified target genes) with expression of 17 434 mRNAs to identify new targets in colon tissue. Seed region matches between miRNA and newly identified targeted mRNA were used to help determine direct miRNA-mRNA associations. RESULTS: From the replication of the 121 miRNAs that had at least 1 database-identified target gene using mRNA expression methods, 97.9% were expressed in normal colonic mucosa. Of the 8622 target miRNA-mRNA associations identified in the database, 2658 (30.2%) were associated with gene expression in normal colonic mucosa after adjusting for multiple comparisons. Of the 133 miRNAs with database-identified target genes by non-mRNA expression methods, 97.2% were expressed in normal colonic mucosa. After adjustment for multiple comparisons, 2416 miRNA-mRNA associations remained significant (19.8%). Results from the discovery phase based on detailed examination of 18 miRNAs identified more than 80 000 miRNA-mRNA associations that had not previously linked to the miRNA. Of these miRNA-mRNA associations, 15.6% and 14.8% had seed matches for CRCh38 and CRCh37, respectively. CONCLUSIONS: Our data suggest that miRNA target gene databases are incomplete; pathways derived from these databases have similar deficiencies. Although we know a lot about several miRNAs, little is known about other miRNAs in terms of their targeted genes. We encourage others to use their data to continue to further identify and validate miRNA-targeted genes.

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.

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.

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.

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.

Diet and lifestyle factors associated with miRNA expression in colorectal tissue.

MicroRNAs (miRNAs) are small non-protein-coding RNA molecules that regulate gene expression. Diet and lifestyle factors have been hypothesized to be involved in the regulation of miRNA expression. In this study it was hypothesized that diet and lifestyle factors are associated with miRNA expression. Data from 1,447 cases of colorectal cancer to evaluate 34 diet and lifestyle variables using miRNA expression in normal colorectal mucosa as well as for differential expression between paired carcinoma and normal tissue were used. miRNA data were obtained using an Agilent platform. Multiple comparisons were adjusted for using the false discovery rate q-value. There were 250 miRNAs differentially expressed between carcinoma and normal colonic tissue by level of carbohydrate intake and 198 miRNAs differentially expressed by the level of sucrose intake. Of these miRNAs, 166 miRNAs were differentially expressed for both carbohydrate intake and sucrose intake. Ninety-nine miRNAs were differentially expressed by the level of whole grain intake in normal colonic mucosa. Level of oxidative balance score was associated with 137 differentially expressed miRNAs between carcinoma and paired normal rectal mucosa. Additionally, 135 miRNAs were differentially expressed in colon tissue based on recent NSAID use. Other dietary factors, body mass index, waist and hip circumference, and long-term physical activity levels did not alter miRNA expression after adjustment for multiple comparisons. These results suggest that diet and lifestyle factors regulate miRNA level. They provide additional support for the influence of carbohydrate, sucrose, whole grains, NSAIDs, and oxidative balance score on colorectal cancer risk.

Genetic variants in the TGFß-signaling pathway influence expression of miRNAs in colon and rectal normal mucosa and tumor tissue.

The TGF-ß signaling pathway is involved in regulation of cell growth, angiogenesis, and metastasis. We test the hypothesis that genetic variation in the TGF-ß signaling pathway alters miRNA expression.We use data from 1188 colorectal cancer cases to evaluate associations between 80 SNPs in 21 genes.Seven variants eIF4E rs12498533, NFκB1 rs230510, TGFB1 rs4803455, TGFBR1 rs1571590 and rs6478974, SMAD3 rs3743343, and RUNX1 rs8134179 were associated with expression level of miRNAs in normal colorectal mucosa. RUNX2 rs12333172 and BMPR1B rs13134042 were associated with miRNAs in normal colon mucosa; eIF4EBP3 rs250425, SMAD3 rs12904944, SMAD7 rs3736242, and PTEN rs532678 were associated with miRNA expression in normal rectal mucosa. Evaluation of the differential expression between carcinoma and normal mucosa showed that SMAD3 rs12708491 and rs2414937, NFκB1 rs230510 and rs3821958, and RUNX3 rs6672420 were associated with several miRNAs for colorectal carcinoma. Evaluation of site-specific differential miRNA expression showed that BMPR1B rs2120834, BMPR2 rs2228545, and eIF4EBP3 rs250425 were associated with differential miRNA expression in colon tissue and SMAD3 rs12901071, rs1498506, and rs2414937, BMPR2 rs2228545, and RUNX2 rs2819854, altered differential miRNA expression in rectal tissue.These data support the importance of the TGF-ß signaling pathway to the carcinogenic process, possibly through their influence on miRNA expression levels.