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.
It has been proposed that miRNAs are involved in the control of telomeres. We test that hypothesis by examining the association between miRNAs and telomere length (TL). Additionally, we evaluate if ...genetic variation in telomerase reverse transcriptase (TERT) is associated with miRNA expression levels. We use data from a population-based study of colorectal cancer (CRC), where we have previously shown associations between TL and TERT and CRC, to test associations between TL and miRNA expression and TERT and miRNA expression. To gain insight into functions of miRNAs associated with TERT we tested linear associations between miRNAs and their targeted gene mRNAs. An Agilent platform that contained information on over 2000 miRNAs was used. TL was measured using a multiplexed quantitative PCR (qPCR). RNAseq was used to assess gene expression. Our sample consisted of 1152 individuals with SNP data and miRNA expression data; 363 individuals with both TL and miRNA; and 148 individuals with miRNA and mRNA data. Thirty-three miRNAs were directly associated with TL after adjusting for age and sex (false discovery rate (FDR) of 0.05). TERT rs2736118 was associated with differences in miRNA expression between carcinoma and normal colonic mucosa for 75 miRNAs (FDR <0.05). Genes regulated by these miRNAs, as indicated by mRNA/miRNA associations, were associated with major signaling pathways beyond their TL-related functions, including PTEN, and PI3K/AKT signaling. Our data support a direct association between miRNAs and TL; differences in miRNA expression levels by TERT genotype were observed. Based on miRNA and targeted mRNA associations our data suggest that TERT is involved in non-TL-related functions by acting through altered miRNA expression.
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.
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.
Missing data can arise in bioinformatics applications for a variety of reasons, and imputation methods are frequently applied to such data. We are motivated by a colorectal cancer study where miRNA ...expression was measured in paired tumor-normal samples of hundreds of patients, but data for many normal samples were missing due to lack of tissue availability. We compare the precision and power performance of several imputation methods, and draw attention to the statistical dependence induced by K-Nearest Neighbors (KNN) imputation. This imputation-induced dependence has not previously been addressed in the literature. We demonstrate how to account for this dependence, and show through simulation how the choice to ignore or account for this dependence affects both power and type I error rate control.
Physical activity has been shown to reduce risk of colon cancer. Some studies have shown site-specific associations while others have not. The inverse association between physical activity and colon ...cancer is consistent although only 7 of 13 studies that have collected both colon and rectal cancer data in the same manner report reduced risk for rectal cancer; four of these studies detected statistically significant inverse associations. The frequency, duration and intensity of activity are important components of a public health message to reduce risk of colon cancer through performance of physical activity. However, difficulties in estimating the exact amount of activity needed and frequency and intensity of activity result in only crude estimates of dose needed for a protective effect. Much of the literature suggest that more intense activity is needed to reduce colon cancer risk and that somewhere between 3.5 and 4 hours of vigorous activity per week may be needed to optimise protection. Several biological mechanisms have been proposed to explain the association between physical activity and colon cancer; many of these mechanisms also support the observation that intense activities are most protective. Biological mechanisms include: physical activity increasing gut motility; enhancing the immune system; decreasing insulin and insulin-like growth factor levels; decreasing obesity; enhancing free radical scavenger systems; and influencing prostaglandin levels. The evidence taken together provides strong support for lack of physical activity being causally related to colon cancer. It has been estimated that 12-14% of colon cancer could be attributed to lack of frequent involvement in vigorous physical activity.
Known genetic loci explain only a small proportion of the familial relative risk of colorectal cancer (CRC). We conducted a genome-wide association study of CRC in East Asians with 14,963 cases and ...31,945 controls and identified 6 new loci associated with CRC risk (P = 3.42 × 10(-8) to 9.22 × 10(-21)) at 10q22.3, 10q25.2, 11q12.2, 12p13.31, 17p13.3 and 19q13.2. Two of these loci map to genes (TCF7L2 and TGFB1) with established roles in colorectal tumorigenesis. Four other loci are located in or near genes involved in transcriptional regulation (ZMIZ1), genome maintenance (FEN1), fatty acid metabolism (FADS1 and FADS2), cancer cell motility and metastasis (CD9), and cell growth and differentiation (NXN). We also found suggestive evidence for three additional loci associated with CRC risk near genome-wide significance at 8q24.11, 10q21.1 and 10q24.2. Furthermore, we replicated 22 previously reported CRC-associated loci. Our study provides insights into the genetic basis of CRC and suggests the involvement of new biological pathways.
MicroRNAs (miRNAs) are small, non-protein-coding RNA molecules that are commonly dysregulated in colorectal tumors. The objective of this study was to identify smaller subsets of highly predictive ...miRNAs.
Data come from population-based studies of colorectal cancer conducted in Utah and the Kaiser Permanente Medical Care Program. Tissue samples were available for 1,953 individuals, of which 1,894 had carcinoma tissue and 1,599 had normal mucosa available for statistical analysis. Agilent Human miRNA Microarray V.19.0 was used to generate miRNA expression profiles; validation of expression levels was carried out using quantitative PCR. We used random forest analysis and verified findings with logistic modeling in separate data sets. Important microRNAs are identified and bioinformatics tools are used to identify target genes and related biological pathways.
We identified 16 miRNAs for colon and 17 miRNAs for rectal carcinoma that appear to differentiate between carcinoma and normal mucosa; of these, 12 were important for both colon and rectal cancer, hsa-miR-663b, hsa-miR-4539, hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-21-5p, hsa-miR-4506, hsa-miR-92a-3p, hsa-miR-93-5p, hsa-miR-145-5p, hsa-miR-3651, hsa-miR-378a-3p, and hsa-miR-378i. Estimated misclassification rates were low at 4.83% and 2.5% among colon and rectal observations, respectively. Among independent observations, logistic modeling reinforced the importance of these miRNAs, finding the primary principal components of their variation statistically significant (P<0.001 among both colon and rectal observations) and again producing low misclassification rates. Repeating our analysis without those miRNAs initially identified as important identified other important miRNAs; however, misclassification rates increased and distinctions between remaining miRNAs in terms of classification importance were reduced.
Our data support the hypothesis that while many miRNAs are dysregulated between carcinoma and normal mucosa, smaller subsets of these miRNAs are useful and informative in discriminating between these tissues.
Colorectal cancer (CRC) accounts for about 8% of all new cancer cases diagnosed in the US. We used whole exome sequence data from triplet samples (colon carcinoma, colon adenoma, and normal tissue) ...from 18 individuals to assess gene mutation rates. Of the 2 204 genes that were mutated, APC, TTN, TP53, KRAS, OBSCN, SOX9, PCDH17, SIGLEC10, MYH6, and BRD9 were consistent with genes being an early driver of carcinogenesis, in that they were mutated in multiple adenomas and multiple carcinomas. Fifty‐two genes were mutated in ≥12.5% of microsatellite stable (MSS) carcinomas but not in any of the adenomas, in line with the profile of a late driver event involved in tumor progression. Thirty‐eight genes were sequenced in a larger independent set of 148 carcinoma/normal tissue pairs to obtain more precise mutation frequencies. Eight of the genes, APC, TP53, ATM, CSMD3, LRP1B, RYR2, BIRC6, and MUC17, contained mutations in >20% of the carcinomas. Interestingly, mutations in four genes in addition to APC that are associated with dysregulation of Wnt signaling, were all classified as early driver events. Most of the genes that are commonly associated with colon cancer, including APC, TP53, and KRAS, were all classified as being early driver genes being mutated in both adenomas and carcinomas. Classifying genes as potential early and late driver events points to candidate genes that may help dissect pathways involved in both tumor initiation and progression.
There is no model to estimate absolute invasive breast cancer risk for Hispanic women.
The San Francisco Bay Area Breast Cancer Study (SFBCS) provided data on Hispanic breast cancer case patients ...(533 US-born, 553 foreign-born) and control participants (464 US-born, 947 foreign-born). These data yielded estimates of relative risk (RR) and attributable risk (AR) separately for US-born and foreign-born women. Nativity-specific absolute risks were estimated by combining RR and AR information with nativity-specific invasive breast cancer incidence and competing mortality rates from the California Cancer Registry and Surveillance, Epidemiology, and End Results program to develop the Hispanic risk model (HRM). In independent data, we assessed model calibration through observed/expected (O/E) ratios, and we estimated discriminatory accuracy with the area under the receiver operating characteristic curve (AUC) statistic.
The US-born HRM included age at first full-term pregnancy, biopsy for benign breast disease, and family history of breast cancer; the foreign-born HRM also included age at menarche. The HRM estimated lower risks than the National Cancer Institute's Breast Cancer Risk Assessment Tool (BCRAT) for US-born Hispanic women, but higher risks in foreign-born women. In independent data from the Women's Health Initiative, the HRM was well calibrated for US-born women (observed/expected O/E ratio = 1.07, 95% confidence interval CI = 0.81 to 1.40), but seemed to overestimate risk in foreign-born women (O/E ratio = 0.66, 95% CI = 0.41 to 1.07). The AUC was 0.564 (95% CI = 0.485 to 0.644) for US-born and 0.625 (95% CI = 0.487 to 0.764) for foreign-born women.
The HRM is the first absolute risk model that is based entirely on data specific to Hispanic women by nativity. Further studies in Hispanic women are warranted to evaluate its validity.