Diabetes mellitus in relation to colorectal tumor molecular subtypes: A pooled analysis of more than 9000 cases.

Diabetes is an established risk factor for colorectal cancer. However, colorectal cancer is a heterogeneous disease and it is not well understood whether diabetes is more strongly associated with some tumor molecular subtypes than others. A better understanding of the association between diabetes and colorectal cancer according to molecular subtypes could provide important insights into the biology of this association. We used data on lifestyle and clinical characteristics from the Colorectal Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), including 9756 colorectal cancer cases (with tumor marker data) and 9985 controls, to evaluate associations between reported diabetes and risk of colorectal cancer according to molecular subtypes. Tumor markers included BRAF and KRAS mutations, microsatellite instability and CpG island methylator phenotype. In the multinomial logistic regression model, comparing colorectal cancer cases to cancer-free controls, diabetes was positively associated with colorectal cancer regardless of subtype. The highest OR estimate was found for BRAF-mutated colorectal cancer, n = 1086 (ORfully adj : 1.67, 95% confidence intervals [CI]: 1.36-2.05), with an attenuated association observed between diabetes and colorectal cancer without BRAF-mutations, n = 7959 (ORfully adj : 1.33, 95% CI: 1.19-1.48). In the case only analysis, BRAF-mutation was differentially associated with diabetes (Pdifference  = .03). For the other markers, associations with diabetes were similar across tumor subtypes. In conclusion, our study confirms the established association between diabetes and colorectal cancer risk, and suggests that it particularly increases the risk of BRAF-mutated tumors.

Comprehensively evaluating cis-regulatory variation in the human prostate transcriptome by using gene-level allele-specific expression.

The identification of cis-acting regulatory variation in primary tissues has the potential to elucidate the genetic basis of complex traits and further our understanding of transcriptomic diversity across cell types. Expression quantitative trait locus (eQTL) association analysis using RNA sequencing (RNA-seq) data can improve upon the detection of cis-acting regulatory variation by leveraging allele-specific expression (ASE) patterns in association analysis. Here, we present a comprehensive evaluation of cis-acting eQTLs by analyzing RNA-seq gene-expression data and genome-wide high-density genotypes from 471 samples of normal primary prostate tissue. Using statistical models that integrate ASE information, we identified extensive cis-eQTLs across the prostate transcriptome and found that approximately 70% of expressed genes corresponded to a significant eQTL at a gene-level false-discovery rate of 0.05. Overall, cis-eQTLs were heavily concentrated near the transcription start and stop sites of affected genes, and effects were negatively correlated with distance. We identified multiple instances of cis-acting co-regulation by using phased genotype data and discovered 233 SNPs as the most strongly associated eQTLs for more than one gene. We also noted significant enrichment (25/50, p = 2E-5) of previously reported prostate cancer risk SNPs in prostate eQTLs. Our results illustrate the benefit of assessing ASE data in cis-eQTL analyses by showing better reproducibility of prior eQTL findings than of eQTL mapping based on total expression alone. Altogether, our analysis provides extensive functional context of thousands of SNPs in prostate tissue, and these results will be of critical value in guiding studies examining disease of the human prostate.

DNA mismatch repair protein deficient non-neoplastic colonic crypts: a novel indicator of Lynch syndrome.

Lynch syndrome is the most common form of hereditary colorectal carcinoma. However, establishing the diagnosis of Lynch syndrome is challenging, and ancillary studies that distinguish between sporadic DNA mismatch repair (MMR) protein deficiency and Lynch syndrome are needed, particularly when germline mutation studies are inconclusive. The aim of this study was to determine if MMR protein-deficient non-neoplastic intestinal crypts can help distinguish between patients with and without Lynch syndrome. We evaluated the expression of MMR proteins in non-neoplastic intestinal mucosa obtained from colorectal surgical resection specimens from patients with Lynch syndrome-associated colorectal carcinoma (n = 52) and patients with colorectal carcinoma without evidence of Lynch syndrome (n = 70), including sporadic MMR protein-deficient colorectal carcinoma (n = 30), MMR protein proficient colorectal carcinoma (n = 30), and "Lynch-like" syndrome (n = 10). MMR protein-deficient non-neoplastic colonic crypts were identified in 19 of 122 (16%) patients. MMR protein-deficient colonic crypts were identified in 18 of 52 (35%) patients with Lynch syndrome compared to only 1 of 70 (1%) patients without Lynch syndrome (p < 0.001). This one patient had "Lynch-like" syndrome and harbored two MSH2-deficient non-neoplastic colonic crypts. MMR protein-deficient non-neoplastic colonic crypts were not identified in patients with sporadic MMR protein-deficient or MMR protein proficient colorectal carcinoma. Our findings suggest that MMR protein-deficient colonic crypts are a novel indicator of Lynch syndrome, and evaluation for MMR protein-deficient crypts may be a helpful addition to Lynch syndrome diagnostics.

FIRE: functional inference of genetic variants that regulate gene expression.

MOTIVATION: Interpreting genetic variation in noncoding regions of the genome is an important challenge for personal genome analysis. One mechanism by which noncoding single nucleotide variants (SNVs) influence downstream phenotypes is through the regulation of gene expression. Methods to predict whether or not individual SNVs are likely to regulate gene expression would aid interpretation of variants of unknown significance identified in whole-genome sequencing studies. RESULTS: We developed FIRE (Functional Inference of Regulators of Expression), a tool to score both noncoding and coding SNVs based on their potential to regulate the expression levels of nearby genes. FIRE consists of 23 random forests trained to recognize SNVs in cis-expression quantitative trait loci (cis-eQTLs) using a set of 92 genomic annotations as predictive features. FIRE scores discriminate cis-eQTL SNVs from non-eQTL SNVs in the training set with a cross-validated area under the receiver operating characteristic curve (AUC) of 0.807, and discriminate cis-eQTL SNVs shared across six populations of different ancestry from non-eQTL SNVs with an AUC of 0.939. FIRE scores are also predictive of cis-eQTL SNVs across a variety of tissue types. AVAILABILITY AND IMPLEMENTATION: FIRE scores for genome-wide SNVs in hg19/GRCh37 are available for download at https://sites.google.com/site/fireregulatoryvariation/. CONTACT: nilah@stanford.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sequential expression of miR-182 and miR-503 cooperatively targets FBXW7, contributing to the malignant transformation of colon adenoma to adenocarcinoma.

Genetic changes in colon cancer are known to parallel the tissue abnormalities associated with the disease, namely adenoma and adenocarcinoma. The role of microRNA dysregulation in dysplastic progression, however, is not well understood. Here, we show that miR-182 and miR-503 undergo sequential up-regulation and drive the progression of colon adenoma to adenocarcinoma by cooperatively down-regulating the tumour suppressor FBXW7. We identified that increased expression of miR-182 is a feature of adenomas. A subsequent increase in miR-503 expression works cooperatively with miR-182 to induce transformation of an adenoma to adenocarcinoma. We show that introducing miR-503 into AAC1 cells, which are derived from a benign adenoma, confers tumourigenic potential. We also demonstrated that blocking both miR-182 and miR-503 in HCT116 colon cancer cells resulted in increased FBXW7 expression and significantly reduced tumour size in xenograft models. We confirmed relevance of these results in patients by examining the expression levels of miR-182 and miR-503 in over 200 colon cancer patients with 12 year survival outcome data. Decreased patient survival was correlated with elevated expression of both miRNAs, suggesting that elevated levels of both miR-182 and miR-503 define a novel prognostic biomarker for colon cancer patients. In conclusion, we show that a sequential expression of miR-182 and miR-503 in benign adenoma cooperatively regulates the tumour suppressor FBXW7, contributing to the malignant transformation of colon adenoma to adenocarcinoma and miR-182 and miR-503 may prove to be novel therapeutic targets. Array data are available at: http://www.oncomir.umn.edu/

CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer.

Aberrant DNA methylation of CpG islands has been widely observed in human colorectal tumors and is associated with gene silencing when it occurs in promoter areas. A subset of colorectal tumors has an exceptionally high frequency of methylation of some CpG islands, leading to the suggestion of a distinct trait referred to as 'CpG island methylator phenotype', or 'CIMP'. However, the existence of CIMP has been challenged. To resolve this continuing controversy, we conducted a systematic, stepwise screen of 195 CpG island methylation markers using MethyLight technology, involving 295 primary human colorectal tumors and 16,785 separate quantitative analyses. We found that CIMP-positive (CIMP+) tumors convincingly represent a distinct subset, encompassing almost all cases of tumors with BRAF mutation (odds ratio = 203). Sporadic cases of mismatch repair deficiency occur almost exclusively as a consequence of CIMP-associated methylation of MLH1 . We propose a robust new marker panel to classify CIMP+ tumors.

Associations between colorectal cancer molecular markers and pathways with clinicopathologic features in older women.

BACKGROUND & AIMS: Colorectal tumors have a large degree of molecular heterogeneity. Three integrated pathways of carcinogenesis (ie, traditional, alternate, and serrated) have been proposed, based on specific combinations of microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and mutations in BRAF and KRAS. We used resources from the population-based Iowa Women's Health Study (n = 41,836) to associate markers of colorectal tumors, integrated pathways, and clinical and pathology characteristics, including survival times. METHODS: We assessed archived specimens from 732 incident colorectal tumors and characterized them as microsatellite stable (MSS), MSI high or MSI low, CIMP high or CIMP low, CIMP negative, and positive or negative for BRAF and/or KRAS mutations. Informative marker data were collected from 563 tumors (77%), which were assigned to the following integrated pathways: traditional (MSS, CIMP negative, BRAF mutation negative, and KRAS mutation negative; n = 170), alternate (MSS, CIMP low, BRAF mutation negative, and KRAS mutation positive; n = 58), serrated (any MSI, CIMP high, BRAF mutation positive, and KRAS mutation negative; n = 142), or unassigned (n = 193). Multivariable-adjusted Cox proportional hazards regression models were used to assess the associations of interest. RESULTS: Patients' mean age (P = .03) and tumors' anatomic subsite (P = .0001) and grade (P = .0001) were significantly associated with integrated pathway assignment. Colorectal cancer (CRC) mortality was not associated with the traditional, alternate, or serrated pathways, but was associated with a subset of pathway-unassigned tumors (MSS or MSI low, CIMP negative, BRAF mutation negative, and KRAS mutation positive) (n = 96 cases; relative risk = 1.76; 95% confidence interval, 1.07-2.89, compared with the traditional pathway). CONCLUSIONS: We identified clinical and pathology features associated with molecularly defined CRC subtypes. However, additional studies are needed to determine how these features might influence prognosis.

Genotoxic colibactin mutational signature in colorectal cancer is associated with clinicopathological features, specific genomic alterations and better survival.

Background and Aims: The microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain of Escherichia coli harboring the pks island that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown. Methods: SBS88-positive CRCs were identified from targeted sequencing data from 5,292 CRCs from 17 studies and tested for their association with clinico-pathological features, oncogenic pathways, genomic characteristics and survival. Results: In total, 7.5% (398/5,292) of the CRCs were SBS88-positive, of which 98.7% (392/398) were microsatellite stable/microsatellite instability low (MSS/MSI-L), compared with 80% (3916/4894) of SBS88 negative tumors (p=1.5x10-28). Analysis of MSS/MSI-L CRCs demonstrated that SBS88 positive CRCs were associated with the distal colon (OR=1.84, 95% CI=1.40-2.42, p=1x10-5) and rectum (OR=1.90, 95% CI=1.44-2.51, p=6x10-6) tumor sites compared with the proximal colon. The top seven recurrent somatic mutations associated with SBS88-positive CRCs demonstrated mutational contexts associated with colibactin-induced DNA damage, the strongest of which was the APC:c.835-8A>G mutation (OR=65.5, 95%CI=39.0-110.0, p=3x10-80). Large copy number alterations (CNAs) including CNA loss on 14q and gains on 13q, 16q and 20p were significantly enriched in SBS88-positive CRCs. SBS88-positive CRCs were associated with better CRC-specific survival (p=0.007; hazard ratio of 0.69, 95% CI=0.52-0.90) when stratified by age, sex, study, and by stage. Conclusion: SBS88-positivity, a biomarker of colibactin-induced DNA damage, can identify a novel subtype of CRC characterized by recurrent somatic mutations, copy number alterations and better survival. These findings provide new insights for treatment and prevention strategies for this subtype of CRC.

Body size and risk of colorectal cancer molecular defined subtypes and pathways: Mendelian randomization analyses.

BACKGROUND: Obesity has been positively associated with most molecular subtypes of colorectal cancer (CRC); however, the magnitude and the causality of these associations is uncertain. METHODS: We used Mendelian randomization (MR) to examine potential causal relationships between body size traits (body mass index [BMI], waist circumference, and body fat percentage) with risks of Jass classification types and individual subtypes of CRC (microsatellite instability [MSI] status, CpG island methylator phenotype [CIMP] status, BRAF and KRAS mutations). Summary data on tumour markers were obtained from two genetic consortia (CCFR, GECCO). FINDINGS: A 1-standard deviation (SD:5.1 kg/m2) increment in BMI levels was found to increase risks of Jass type 1MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype (odds ratio [OR]: 2.14, 95% confidence interval [CI]: 1.46, 3.13; p-value = 9 × 10-5) and Jass type 2non-MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype CRC (OR: 2.20, 95% CI: 1.26, 3.86; p-value = 0.005). The magnitude of these associations was stronger compared with Jass type 4non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-wildtype CRC (p-differences: 0.03 and 0.04, respectively). A 1-SD (SD:13.4 cm) increment in waist circumference increased risk of Jass type 3non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-mutated (OR 1.73, 95% CI: 1.34, 2.25; p-value = 9 × 10-5) that was stronger compared with Jass type 4 CRC (p-difference: 0.03). A higher body fat percentage (SD:8.5%) increased risk of Jass type 1 CRC (OR: 2.59, 95% CI: 1.49, 4.48; p-value = 0.001), which was greater than Jass type 4 CRC (p-difference: 0.03). INTERPRETATION: Body size was more strongly linked to the serrated (Jass types 1 and 2) and alternate (Jass type 3) pathways of colorectal carcinogenesis in comparison to the traditional pathway (Jass type 4). FUNDING: Cancer Research UK, National Institute for Health Research, Medical Research Council, National Institutes of Health, National Cancer Institute, American Institute for Cancer Research, Brigham and Women's Hospital, Prevent Cancer Foundation, Victorian Cancer Agency, Swedish Research Council, Swedish Cancer Society, Region Västerbotten, Knut and Alice Wallenberg Foundation, Lion's Cancer Research Foundation, Insamlingsstiftelsen, Umeå University. Full funding details are provided in acknowledgements.

Epidemiologic Factors in Relation to Colorectal Cancer Risk and Survival by Genotoxic Colibactin Mutational Signature.

BACKGROUND: The genotoxin colibactin causes a tumor single-base substitution (SBS) mutational signature, SBS88. It is unknown whether epidemiologic factors' association with colorectal cancer risk and survival differs by SBS88. METHODS: Within the Genetic Epidemiology of Colorectal Cancer Consortium and Colon Cancer Family Registry, we measured SBS88 in 4,308 microsatellite stable/microsatellite instability low tumors. Associations of epidemiologic factors with colorectal cancer risk by SBS88 were assessed using multinomial regression (N = 4,308 cases, 14,192 controls; cohort-only cases N = 1,911), and with colorectal cancer-specific survival using Cox proportional hazards regression (N = 3,465 cases). RESULTS: 392 (9%) tumors were SBS88 positive. Among all cases, the highest quartile of fruit intake was associated with lower risk of SBS88-positive colorectal cancer than SBS88-negative colorectal cancer [odds ratio (OR) = 0.53, 95% confidence interval (CI) 0.37-0.76; OR = 0.75, 95% CI 0.66-0.85, respectively, Pheterogeneity = 0.047]. Among cohort studies, associations of body mass index (BMI), alcohol, and fruit intake with colorectal cancer risk differed by SBS88. BMI ≥30 kg/m2 was associated with worse colorectal cancer-specific survival among those SBS88-positive [hazard ratio (HR) = 3.40, 95% CI 1.47-7.84], but not among those SBS88-negative (HR = 0.97, 95% CI 0.78-1.21, Pheterogeneity = 0.066). CONCLUSIONS: Most epidemiologic factors did not differ by SBS88 for colorectal cancer risk or survival. Higher BMI may be associated with worse colorectal cancer-specific survival among those SBS88-positive; however, validation is needed in samples with whole-genome or whole-exome sequencing available. IMPACT: This study highlights the importance of identification of tumor phenotypes related to colorectal cancer and understanding potential heterogeneity for risk and survival.

Folate intake and colorectal cancer risk according to genetic subtypes defined by targeted tumor sequencing.

BACKGROUND: Folate is involved in multiple genetic, epigenetic, and metabolic processes, and inadequate folate intake has been associated with an increased risk of cancer. OBJECTIVE: We examined whether folate intake is differentially associated with colorectal cancer (CRC) risk according to somatic mutations in genes linked to CRC using targeted sequencing. DESIGN: Participants within 2 large CRC consortia with available information on dietary folate, supplemental folic acid, and total folate intake were included. Colorectal tumor samples from cases were sequenced for the presence of nonsilent mutations in 105 genes and 6 signaling pathways (IGF2/PI3K, MMR, RTK/RAS, TGF-ß, WNT, and TP53/ATM). Multinomial logistic regression models were analyzed comparing mutated/nonmutated CRC cases to controls to compute multivariable-adjusted odds ratios (ORs) with 95% confidence interval (CI). Heterogeneity of associations of mutated compared with nonmutated CRC cases was tested in case-only analyses using logistic regression. Analyses were performed separately in hypermutated and nonhypermutated tumors, because they exhibit different clinical behaviors. RESULTS: We included 4339 CRC cases (702 hypermutated tumors, 16.2%) and 11,767 controls. Total folate intake was inversely associated with CRC risk (OR = 0.93; 95% CI: 0.90, 0.96). Among hypermutated tumors, 12 genes (AXIN2, B2M, BCOR, CHD1, DOCK3, FBLN2, MAP3K21, POLD1, RYR1, TET2, UTP20, and ZNF521) showed nominal statistical significance (P < 0.05) for heterogeneity by mutation status, but none remained significant after multiple testing correction. Among these genetic subtypes, the associations between folate variables and CRC were mostly inverse or toward the null, except for tumors mutated for DOCK3 (supplemental folic acid), CHD1 (total folate), and ZNF521 (dietary folate) that showed positive associations. We did not observe differential associations in analyses among nonhypermutated tumors, or according to the signaling pathways. CONCLUSIONS: Folate intake was not differentially associated with CRC risk according to mutations in the genes explored. The nominally significant differential mutation effects observed in a few genes warrants further investigation.

Postmenopausal hormone therapy and colorectal cancer risk by molecularly defined subtypes among older women.

BACKGROUND: Postmenopausal hormone (PMH) therapy may reduce colorectal cancer (CRC) risk, but existing data are inconclusive. OBJECTIVES: To evaluate associations between PMH therapy and incident CRC, overall and by molecularly defined subtypes, in the population-based Iowa Women's Health Study of older women. METHODS: Exposure data were collected from Iowa Women's Health Study participants (55-69 years) at baseline (1986). Archived, paraffin-embedded tissue specimens for 553 CRC cases were collected and analysed to determine microsatellite instability (MSI-L/MSS or MSI-H), CpG island methylator phenotype (CIMP-negative or CIMP-positive) and BRAF mutation (BRAF-wildtype or BRAF-mutated) status. Multivariable Cox regression models were fit to estimate RRs and 95% CIs. RESULTS: PMH therapy (ever vs never use) was inversely associated with incident CRC overall (RR=0.82; 95% CI 0.72 to 0.93), with a significantly lower risk for MSI-L/MSS tumours (RR=0.75; 95% CI 0.60 to 0.94), and borderline significantly lower risks for CIMP-negative (RR=0.79; 95% CI 0.63 to 1.01) and BRAF-wildtype (RR=0.83; 95% CI 0.66 to 1.04) tumours. For PMH therapy >5 years, the subtype-specific risk estimates for MSI-L/MSS, CIMP-negative and BRAF-wildtype tumours were: RR=0.60, 95% CI 0.40 to 0.91; RR=0.68, 95% CI 0.45 to 1.03; and RR=0.70, 95% CI 0.47 to 1.05, respectively. PMH therapy was not significantly associated with the MSI-H, CIMP-positive or BRAF-mutated CRC subtypes. CONCLUSIONS: In this prospective cohort study, PMH therapy was inversely associated with distinct molecularly defined CRC subtypes, which may be related to differential effects from oestrogen and/or progestin on heterogeneous pathways of colorectal carcinogenesis.

A High-Throughput Workflow for FFPE Tissue Proteomics.

Laser capture microdissection (LCM) has become an indispensable tool for mass spectrometry-based proteomic analysis of specific regions obtained from formalin-fixed paraffin-embedded (FFPE) tissue samples in both clinical and research settings. Low protein yields from LCM samples along with laborious sample processing steps present challenges for proteomic analysis without sacrificing protein and peptide recovery. Automation of sample preparation workflows is still under development, especially for samples such as laser-capture microdissected tissues. Here, we present a simplified and rapid workflow using adaptive focused acoustics (AFA) technology for sample processing for high-throughput FFPE-based proteomics. We evaluated three different workflows: standard extraction method followed by overnight trypsin digestion, AFA-assisted extraction and overnight trypsin digestion, and AFA-assisted extraction simultaneously performed with trypsin digestion. The use of AFA-based ultrasonication enables automated sample processing for high-throughput proteomic analysis of LCM-FFPE tissues in 96-well and 384-well formats. Further, accelerated trypsin digestion combined with AFA dramatically reduced the overall processing times. LC-MS/MS analysis revealed a slightly higher number of protein and peptide identifications in AFA accelerated workflows compared to standard and AFA overnight workflows. Further, we did not observe any difference in the proportion of peptides identified with missed cleavages or deamidated peptides across the three different workflows. Overall, our results demonstrate that the workflow described in this study enables rapid and high-throughput sample processing with greatly reduced sample handling, which is amenable to automation.

Body mass index and molecular subtypes of colorectal cancer.

BACKGROUND: Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease. METHODS: We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables. RESULTS: Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control). CONCLUSIONS: In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

Cigarette smoking and colorectal cancer risk by KRAS mutation status among older women.

OBJECTIVES: Existing data support a modest association between cigarette smoking and incident colorectal cancer (CRC) overall. In this study, we evaluated associations between cigarette smoking and CRC risk stratified by KRAS mutation status, using data and tissue resources from the Iowa Women's Health Study (IWHS). METHODS: The IWHS is a population-based cohort study of cancer incidence among 41,836 randomly selected Iowa women, ages 55-69 years of age at enrollment (1986). Exposure data, including cigarette smoking, were obtained by self-report at baseline. Incident CRCs (n=1,233) were ascertained by annual linkage with the Iowa Cancer Registry. Archived tissue specimens from CRC cases recorded through 2002 were recently requested for molecular epidemiology investigations. Tumor KRAS mutation status was determined by direct sequencing of exon 2, with informative results in 507/555 (91%) available CRC cases (342 mutation negative and 165 mutation positive). Multivariate Cox regression models were fit to estimate relative risks (RRs) and 95% confidence intervals (CIs) for associations between cigarette smoking variables and KRAS-defined CRC subtypes. RESULTS: Multiple smoking variables were associated with increased risk for KRAS mutation-negative tumors, including age at initiation (P=0.02), average number of cigarettes per day (P=0.01), cumulative pack-years (P=0.05), and induction period (P=0.04), with the highest point estimate observed for women who smoked ≥40 cigarettes per day on average (RR=2.38; 95% CI=1.25-4.51; compared with never smokers). Further consideration of CRC subsite suggested that cigarette smoking may be a stronger risk factor for KRAS mutation-negative tumors located in the proximal colon than in the distal colorectum. None of the smoking variables were significantly associated with KRAS mutation-positive CRCs (overall or stratified by anatomic subsite). CONCLUSIONS: Data from this prospective study of older women demonstrate differential associations between cigarette smoking and CRC subtypes defined by KRAS mutation status, and are consistent with the hypothesis that smoking adversely affects the serrated pathway of colorectal carcinogenesis.

Associations between intake of folate and related micronutrients with molecularly defined colorectal cancer risks in the Iowa Women's Health Study.

Folate and related micronturients may affect colorectal cancer (CRC) risk, but the molecular mechanism(s) remain incompletely defined. We analyzed associations between dietary folate, vitamin B6, vitamin B12, and methionine with incident CRC, overall and by microsatellite instability (MSS/MSI-L or MSI-H), CpG island methylator phenotype (CIMP-negative or CIMP-positive), BRAF mutation (negative or positive), and KRAS mutation (negative or positive) status in the prospective, population-based Iowa Women's Health Study (IWHS; 55-69 years at baseline; n = 41,836). Intake estimates were obtained from baseline, self-reported food frequency questionnaires. Molecular marker data were obtained for 514 incident CRC cases. Folate intake was inversely associated with overall CRC risk in age-adjusted Cox regression models, whereas methionine intake was inversely associated with overall CRC risk in multivariable-adjusted models [relative risk (RR) = 0.81; 95% CI = 0.69-0.95; P trend = 0.001 and RR = 0.72; 95% CI = 0.54-0.96; P trend = 0.03 for highest vs. lowest quartiles, respectively]. None of the dietary exposures were associated with MSI, CIMP, BRAF, or KRAS defined CRC subtypes. These data provide minimal support for major effects from the examined micronutrients on overall or molecularly defined CRC risks in the IWHS cohort.

A microRNA Transcriptome-wide Association Study of Prostate Cancer Risk.

Large genome-wide association studies have identified hundreds of single-nucleotide polymorphisms associated with increased risk of prostate cancer (PrCa), and many of these risk loci is presumed to confer regulatory effects on gene expression. While eQTL studies of long RNAs has yielded many potential risk genes, the relationship between PrCa risk genetics and microRNA expression dysregulation is understudied. We performed an microRNA transcriptome-wide association study of PrCa risk using small RNA sequencing and genome-wide genotyping data from N = 441 normal prostate epithelium tissue samples along with N = 411 prostate adenocarcinoma tumor samples from the Cancer Genome Atlas (TCGA). Genetically regulated expression prediction models were trained for all expressed microRNAs using the FUSION TWAS software. TWAS for PrCa risk was performed with both sets of models using single-SNP summary statistics from the recent PRACTICAL consortium PrCa case-control OncoArray GWAS meta-analysis. A total of 613 and 571 distinct expressed microRNAs were identified in the normal and tumor tissue datasets, respectively (overlap: 480). Among these, 79 (13%) normal tissue microRNAs demonstrated significant cis-heritability (median cis-h2 = 0.15, range: 0.03-0.79) for model training. Similar results were obtained from TCGA tumor samples, with 48 (9%) microRNA expression models successfully trained (median cis-h2 = 0.14, range: 0.06-0.60). Using normal tissue models, we identified two significant TWAS microRNA associations with PrCa risk: over-expression of mir-941 family microRNAs (PTWAS = 2.9E-04) and reduced expression of miR-3617-5p (PTWAS = 1.0E-03). The TCGA tumor TWAS also identified a significant association with miR-941 overexpression (PTWAS = 9.7E-04). Subsequent finemapping of the TWAS results using a multi-tissue database indicated limited evidence of causal status for each microRNA with PrCa risk (posterior inclusion probabilities <0.05). Future work will examine downstream regulatory effects of microRNA dysregulation as well as microRNA-mediated risk mechanisms via competing endogenous RNA relationships.

Identifying colorectal cancer caused by biallelic MUTYH pathogenic variants using tumor mutational signatures.

Carriers of germline biallelic pathogenic variants in the MUTYH gene have a high risk of colorectal cancer. We test 5649 colorectal cancers to evaluate the discriminatory potential of a tumor mutational signature specific to MUTYH for identifying biallelic carriers and classifying variants of uncertain clinical significance (VUS). Using a tumor and matched germline targeted multi-gene panel approach, our classifier identifies all biallelic MUTYH carriers and all known non-carriers in an independent test set of 3019 colorectal cancers (accuracy = 100% (95% confidence interval 99.87-100%)). All monoallelic MUTYH carriers are classified with the non-MUTYH carriers. The classifier provides evidence for a pathogenic classification for two VUS and a benign classification for five VUS. Somatic hotspot mutations KRAS p.G12C and PIK3CA p.Q546K are associated with colorectal cancers from biallelic MUTYH carriers compared with non-carriers (p = 2 × 10-23 and p = 6 × 10-11, respectively). Here, we demonstrate the potential application of mutational signatures to tumor sequencing workflows to improve the identification of biallelic MUTYH carriers.

Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival.

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR = 0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR = 1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features.

Human colon cancer profiles show differential microRNA expression depending on mismatch repair status and are characteristic of undifferentiated proliferative states.

BACKGROUND: Colon cancer arises from the accumulation of multiple genetic and epigenetic alterations to normal colonic tissue. microRNAs (miRNAs) are small, non-coding regulatory RNAs that post-transcriptionally regulate gene expression. Differential miRNA expression in cancer versus normal tissue is a common event and may be pivotal for tumor onset and progression. METHODS: To identify miRNAs that are differentially expressed in tumors and tumor subtypes, we carried out highly sensitive expression profiling of 735 miRNAs on samples obtained from a statistically powerful set of tumors (n = 80) and normal colon tissue (n = 28) and validated a subset of this data by qRT-PCR. RESULTS: Tumor specimens showed highly significant and large fold change differential expression of the levels of 39 miRNAs including miR-135b, miR-96, miR-182, miR-183, miR-1, and miR-133a, relative to normal colon tissue. Significant differences were also seen in 6 miRNAs including miR-31 and miR-592, in the direct comparison of tumors that were deficient or proficient for mismatch repair. Examination of the genomic regions containing differentially expressed miRNAs revealed that they were also differentially methylated in colon cancer at a far greater rate than would be expected by chance. A network of interactions between these miRNAs and genes associated with colon cancer provided evidence for the role of these miRNAs as oncogenes by attenuation of tumor suppressor genes. CONCLUSION: Colon tumors show differential expression of miRNAs depending on mismatch repair status. miRNA expression in colon tumors has an epigenetic component and altered expression that may reflect a reversion to regulatory programs characteristic of undifferentiated proliferative developmental states.