Genome-wide interaction study of dietary intake of fibre, fruits, and vegetables with risk of colorectal cancer.

BACKGROUND: Consumption of fibre, fruits and vegetables have been linked with lower colorectal cancer (CRC) risk. A genome-wide gene-environment (G × E) analysis was performed to test whether genetic variants modify these associations. METHODS: A pooled sample of 45 studies including up to 69,734 participants (cases: 29,896; controls: 39,838) of European ancestry were included. To identify G × E interactions, we used the traditional 1--degree-of-freedom (DF) G × E test and to improve power a 2-step procedure and a 3DF joint test that investigates the association between a genetic variant and dietary exposure, CRC risk and G × E interaction simultaneously. FINDINGS: The 3-DF joint test revealed two significant loci with p-value <5 × 10-8. Rs4730274 close to the SLC26A3 gene showed an association with fibre (p-value: 2.4 × 10-3) and G × fibre interaction with CRC (OR per quartile of fibre increase = 0.87, 0.80, and 0.75 for CC, TC, and TT genotype, respectively; G × E p-value: 1.8 × 10-7). Rs1620977 in the NEGR1 gene showed an association with fruit intake (p-value: 1.0 × 10-8) and G × fruit interaction with CRC (OR per quartile of fruit increase = 0.75, 0.65, and 0.56 for AA, AG, and GG genotype, respectively; G × E -p-value: 0.029). INTERPRETATION: We identified 2 loci associated with fibre and fruit intake that also modify the association of these dietary factors with CRC risk. Potential mechanisms include chronic inflammatory intestinal disorders, and gut function. However, further studies are needed for mechanistic validation and replication of findings. FUNDING: National Institutes of Health, National Cancer Institute. Full funding details for the individual consortia are provided in acknowledgments.

Two genome-wide interaction loci modify the association of nonsteroidal anti-inflammatory drugs with colorectal cancer.

Regular, long-term aspirin use may act synergistically with genetic variants, particularly those in mechanistically relevant pathways, to confer a protective effect on colorectal cancer (CRC) risk. We leveraged pooled data from 52 clinical trial, cohort, and case-control studies that included 30,806 CRC cases and 41,861 controls of European ancestry to conduct a genome-wide interaction scan between regular aspirin/nonsteroidal anti-inflammatory drug (NSAID) use and imputed genetic variants. After adjusting for multiple comparisons, we identified statistically significant interactions between regular aspirin/NSAID use and variants in 6q24.1 (top hit rs72833769), which has evidence of influencing expression of TBC1D7 (a subunit of the TSC1-TSC2 complex, a key regulator of MTOR activity), and variants in 5p13.1 (top hit rs350047), which is associated with expression of PTGER4 (codes a cell surface receptor directly involved in the mode of action of aspirin). Genetic variants with functional impact may modulate the chemopreventive effect of regular aspirin use, and our study identifies putative previously unidentified targets for additional mechanistic interrogation.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Disentangling the relationships of body mass index and circulating sex hormone concentrations in mammographic density using Mendelian randomization.

PURPOSE: Mammographic density phenotypes, adjusted for age and body mass index (BMI), are strong predictors of breast cancer risk. BMI is associated with mammographic density measures, but the role of circulating sex hormone concentrations is less clear. We investigated the relationship between BMI, circulating sex hormone concentrations, and mammographic density phenotypes using Mendelian randomization (MR). METHODS: We applied two-sample MR approaches to assess the association between genetically predicted circulating concentrations of sex hormones [estradiol, testosterone, sex hormone-binding globulin (SHBG)], BMI, and mammographic density phenotypes (dense and non-dense area). We created instrumental variables from large European ancestry-based genome-wide association studies and applied estimates to mammographic density phenotypes in up to 14,000 women of European ancestry. We performed analyses overall and by menopausal status. RESULTS: Genetically predicted BMI was positively associated with non-dense area (IVW: ß = 1.79; 95% CI = 1.58, 2.00; p = 9.57 × 10-63) and inversely associated with dense area (IVW: ß = - 0.37; 95% CI = - 0.51,- 0.23; p = 4.7 × 10-7). We observed weak evidence for an association of circulating sex hormone concentrations with mammographic density phenotypes, specifically inverse associations between genetically predicted testosterone concentration and dense area (ß = - 0.22; 95% CI = - 0.38, - 0.053; p = 0.009) and between genetically predicted estradiol concentration and non-dense area (ß = - 3.32; 95% CI = - 5.83, - 0.82; p = 0.009), although results were not consistent across a range of MR approaches. CONCLUSION: Our findings support a positive causal association between BMI and mammographic non-dense area and an inverse association between BMI and dense area. Evidence was weaker and inconsistent for a causal effect of circulating sex hormone concentrations on mammographic density phenotypes. Based on our findings, associations between circulating sex hormone concentrations and mammographic density phenotypes are weak at best.

Genetic risk impacts the association of menopausal hormone therapy with colorectal cancer risk.

BACKGROUND: Menopausal hormone therapy (MHT), a common treatment to relieve symptoms of menopause, is associated with a lower risk of colorectal cancer (CRC). To inform CRC risk prediction and MHT risk-benefit assessment, we aimed to evaluate the joint association of a polygenic risk score (PRS) for CRC and MHT on CRC risk. METHODS: We used data from 28,486 postmenopausal women (11,519 cases and 16,967 controls) of European descent. A PRS based on 141 CRC-associated genetic variants was modeled as a categorical variable in quartiles. Multiplicative interaction between PRS and MHT use was evaluated using logistic regression. Additive interaction was measured using the relative excess risk due to interaction (RERI). 30-year cumulative risks of CRC for 50-year-old women according to MHT use and PRS were calculated. RESULTS: The reduction in odds ratios by MHT use was larger in women within the highest quartile of PRS compared to that in women within the lowest quartile of PRS (p-value = 2.7 × 10-8). At the highest quartile of PRS, the 30-year CRC risk was statistically significantly lower for women taking any MHT than for women not taking any MHT, 3.7% (3.3%-4.0%) vs 6.1% (5.7%-6.5%) (difference 2.4%, P-value = 1.83 × 10-14); these differences were also statistically significant but smaller in magnitude in the lowest PRS quartile, 1.6% (1.4%-1.8%) vs 2.2% (1.9%-2.4%) (difference 0.6%, P-value = 1.01 × 10-3), indicating 4 times greater reduction in absolute risk associated with any MHT use in the highest compared to the lowest quartile of genetic CRC risk. CONCLUSIONS: MHT use has a greater impact on the reduction of CRC risk for women at higher genetic risk. These findings have implications for the development of risk prediction models for CRC and potentially for the consideration of genetic information in the risk-benefit assessment of MHT use.

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.

Robust best linear weighted estimator with missing covariates in survival analysis.

Missing data in covariates can result in biased estimates and loss of power to detect associations. We consider Cox regression in which some covariates are subject to missing. The inverse probability weighted approach is often applied to regression analysis with missing covariates. Inverse probability weighted estimators typically are less efficient than likelihood-based estimators, but in general are more robust against model misspecification. In this article, we propose a robust best linear weighted estimator for Cox regression with missing covariates. Our proposed estimator is the projection of the simple inverse probability weighted estimator onto the orthogonal complement of the score space based on a working regression model of the observed data. The efficiency gain is from the use of the association between the survival outcome variable and the available covariates, which is the working regression model. The asymptotic distribution is derived, and the finite sample performance of the proposed estimator is examined via extensive simulation studies. The methods are applied to a colorectal cancer study to assess the association of the microsatellite instability status with colorectal cancer-specific mortality.

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.

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.

Genome-Wide Gene-Environment Interaction Analyses to Understand the Relationship between Red Meat and Processed Meat Intake and Colorectal Cancer Risk.

BACKGROUND: High red meat and/or processed meat consumption are established colorectal cancer risk factors. We conducted a genome-wide gene-environment (GxE) interaction analysis to identify genetic variants that may modify these associations. METHODS: A pooled sample of 29,842 colorectal cancer cases and 39,635 controls of European ancestry from 27 studies were included. Quantiles for red meat and processed meat intake were constructed from harmonized questionnaire data. Genotyping arrays were imputed to the Haplotype Reference Consortium. Two-step EDGE and joint tests of GxE interaction were utilized in our genome-wide scan. RESULTS: Meta-analyses confirmed positive associations between increased consumption of red meat and processed meat with colorectal cancer risk [per quartile red meat OR = 1.30; 95% confidence interval (CI) = 1.21-1.41; processed meat OR = 1.40; 95% CI = 1.20-1.63]. Two significant genome-wide GxE interactions for red meat consumption were found. Joint GxE tests revealed the rs4871179 SNP in chromosome 8 (downstream of HAS2); greater than median of consumption ORs = 1.38 (95% CI = 1.29-1.46), 1.20 (95% CI = 1.12-1.27), and 1.07 (95% CI = 0.95-1.19) for CC, CG, and GG, respectively. The two-step EDGE method identified the rs35352860 SNP in chromosome 18 (SMAD7 intron); greater than median of consumption ORs = 1.18 (95% CI = 1.11-1.24), 1.35 (95% CI = 1.26-1.44), and 1.46 (95% CI = 1.26-1.69) for CC, CT, and TT, respectively. CONCLUSIONS: We propose two novel biomarkers that support the role of meat consumption with an increased risk of colorectal cancer. IMPACT: The reported GxE interactions may explain the increased risk of colorectal cancer in certain population subgroups.

Genome-wide interaction analysis of folate for colorectal cancer risk.

BACKGROUND: Epidemiological and experimental evidence suggests that higher folate intake is associated with decreased colorectal cancer (CRC) risk; however, the mechanisms underlying this relationship are not fully understood. Genetic variation that may have a direct or indirect impact on folate metabolism can provide insights into folate's role in CRC. OBJECTIVES: Our aim was to perform a genome-wide interaction analysis to identify genetic variants that may modify the association of folate on CRC risk. METHODS: We applied traditional case-control logistic regression, joint 3-degree of freedom, and a 2-step weighted hypothesis approach to test the interactions of common variants (allele frequency >1%) across the genome and dietary folate, folic acid supplement use, and total folate in relation to risk of CRC in 30,550 cases and 42,336 controls from 51 studies from 3 genetic consortia (CCFR, CORECT, GECCO). RESULTS: Inverse associations of dietary, total folate, and folic acid supplement with CRC were found (odds ratio [OR]: 0.93; 95% confidence interval [CI]: 0.90, 0.96; and 0.91; 95% CI: 0.89, 0.94 per quartile higher intake, and 0.82 (95% CI: 0.78, 0.88) for users compared with nonusers, respectively). Interactions (P-interaction < 5×10-8) of folic acid supplement and variants in the 3p25.2 locus (in the region of Synapsin II [SYN2]/tissue inhibitor of metalloproteinase 4 [TIMP4]) were found using traditional interaction analysis, with variant rs150924902 (located upstream to SYN2) showing the strongest interaction. In stratified analyses by rs150924902 genotypes, folate supplementation was associated with decreased CRC risk among those carrying the TT genotype (OR: 0.82; 95% CI: 0.79, 0.86) but increased CRC risk among those carrying the TA genotype (OR: 1.63; 95% CI: 1.29, 2.05), suggesting a qualitative interaction (P-interaction = 1.4×10-8). No interactions were observed for dietary and total folate. CONCLUSIONS: Variation in 3p25.2 locus may modify the association of folate supplement with CRC risk. Experimental studies and studies incorporating other relevant omics data are warranted to validate this finding.

Genome-wide study of genetic polymorphisms predictive for outcome from first-line oxaliplatin-based chemotherapy in colorectal cancer patients.

We conducted the first large genome-wide association study to identify novel genetic variants that predict better (or poorer) prognosis in colorectal cancer patients receiving standard first-line oxaliplatin-based chemotherapy vs chemotherapy without oxaliplatin. We used data from two phase III trials, NCCTG N0147 and NCCTG N9741 and a population-based patient cohort, DACHS. Multivariable Cox proportional hazards models were employed, including an interaction term between each SNP and type of treatment for overall survival (OS) and progression-free survival. The analysis was performed for studies individually, and the results were combined using fixed-effect meta-analyses separately for resected stage III colon cancer (3098 patients from NCCTG N0147 and 549 patients from DACHS) and mCRC (505 patients from NCCTG N9741 and 437 patients from DACHS). We further performed gene-based analysis as well as in silico bioinformatics analysis for CRC-relevant functional genomic annotation of identified loci. In stage III colon cancer patients, a locus on chr22 (rs11912167) was associated with significantly poorer OS after oxaliplatin-based chemotherapy vs chemotherapy without oxaliplatin (Pinteraction < 5 × 10-8 ). For mCRC patients, three loci on chr1 (rs1234556), chr12 (rs11052270) and chr15 (rs11858406) were found to be associated with differential OS (P < 5 × 10-7 ). The locus on chr1 located in the intronic region of RCSD1 was replicated in an independent cohort of 586 mCRC patients from ALGB/SWOG 80405 (Pinteraction = .04). The GWA gene-based analysis yielded for RCSD1 the most significant association with differential OS in mCRC (P = 6.6 × 10-6 ). With further investigation into its biological mechanisms, this finding could potentially be used to individualize first-line treatment and improve clinical outcomes.

Probing the diabetes and colorectal cancer relationship using gene - environment interaction analyses.

BACKGROUND: Diabetes is an established risk factor for colorectal cancer. However, the mechanisms underlying this relationship still require investigation and it is not known if the association is modified by genetic variants. To address these questions, we undertook a genome-wide gene-environment interaction analysis. METHODS: We used data from 3 genetic consortia (CCFR, CORECT, GECCO; 31,318 colorectal cancer cases/41,499 controls) and undertook genome-wide gene-environment interaction analyses with colorectal cancer risk, including interaction tests of genetics(G)xdiabetes (1-degree of freedom; d.f.) and joint testing of Gxdiabetes, G-colorectal cancer association (2-d.f. joint test) and G-diabetes correlation (3-d.f. joint test). RESULTS: Based on the joint tests, we found that the association of diabetes with colorectal cancer risk is modified by loci on chromosomes 8q24.11 (rs3802177, SLC30A8 - ORAA: 1.62, 95% CI: 1.34-1.96; ORAG: 1.41, 95% CI: 1.30-1.54; ORGG: 1.22, 95% CI: 1.13-1.31; p-value3-d.f.: 5.46 × 10-11) and 13q14.13 (rs9526201, LRCH1 - ORGG: 2.11, 95% CI: 1.56-2.83; ORGA: 1.52, 95% CI: 1.38-1.68; ORAA: 1.13, 95% CI: 1.06-1.21; p-value2-d.f.: 7.84 × 10-09). DISCUSSION: These results suggest that variation in genes related to insulin signaling (SLC30A8) and immune function (LRCH1) may modify the association of diabetes with colorectal cancer risk and provide novel insights into the biology underlying the diabetes and colorectal cancer relationship.

Using DEPendency of Association on the Number of Top Hits (DEPTH) as a Complementary Tool to Identify Novel Colorectal Cancer Susceptibility Loci.

BACKGROUND: DEPendency of association on the number of Top Hits (DEPTH) is an approach to identify candidate susceptibility regions by considering the risk signals from overlapping groups of sequential variants across the genome. METHODS: We applied a DEPTH analysis using a sliding window of 200 SNPs to colorectal cancer data from the Colon Cancer Family Registry (CCFR; 5,735 cases and 3,688 controls), and Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO; 8,865 cases and 10,285 controls) studies. A DEPTH score > 1 was used to identify candidate susceptibility regions common to both analyses. We compared DEPTH results against those from conventional genome-wide association study (GWAS) analyses of these two studies as well as against 132 published susceptibility regions. RESULTS: Initial DEPTH analysis revealed 2,622 (CCFR) and 3,686 (GECCO) candidate susceptibility regions, of which 569 were common to both studies. Bootstrapping revealed 40 and 49 candidate susceptibility regions in the CCFR and GECCO data sets, respectively. Notably, DEPTH identified at least 82 regions that would not be detected using conventional GWAS methods, nor had they been identified by previous colorectal cancer GWASs. We found four reproducible candidate susceptibility regions (2q22.2, 2q33.1, 6p21.32, 13q14.3). The highest DEPTH scores were in the human leukocyte antigen locus at 6p21 where the strongest associated SNPs were rs762216297, rs149490268, rs114741460, and rs199707618 for the CCFR data, and rs9270761 for the GECCO data. CONCLUSIONS: DEPTH can identify candidate susceptibility regions for colorectal cancer not identified using conventional analyses of larger datasets. IMPACT: DEPTH has potential as a powerful complementary tool to conventional GWAS analyses for discovering susceptibility regions within the genome.

A Genetic Locus within the FMN1/GREM1 Gene Region Interacts with Body Mass Index in Colorectal Cancer Risk.

Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer. SIGNIFICANCE: This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions.

BACKGROUND: The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. METHODS: We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. RESULTS: We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. CONCLUSIONS: Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Validation of a Genetic-Enhanced Risk Prediction Model for Colorectal Cancer in a Large Community-Based Cohort.

BACKGROUND: Polygenic risk scores (PRS) which summarize individuals' genetic risk profile may enhance targeted colorectal cancer screening. A critical step towards clinical implementation is rigorous external validations in large community-based cohorts. This study externally validated a PRS-enhanced colorectal cancer risk model comprising 140 known colorectal cancer loci to provide a comprehensive assessment on prediction performance. METHODS: The model was developed using 20,338 individuals and externally validated in a community-based cohort (n = 85,221). We validated predicted 5-year absolute colorectal cancer risk, including calibration using expected-to-observed case ratios (E/O) and calibration plots, and discriminatory accuracy using time-dependent AUC. The PRS-related improvement in AUC, sensitivity and specificity were assessed in individuals of age 45 to 74 years (screening-eligible age group) and 40 to 49 years with no endoscopy history (younger-age group). RESULTS: In European-ancestral individuals, the predicted 5-year risk calibrated well [E/O = 1.01; 95% confidence interval (CI), 0.91-1.13] and had high discriminatory accuracy (AUC = 0.73; 95% CI, 0.71-0.76). Adding the PRS to a model with age, sex, family and endoscopy history improved the 5-year AUC by 0.06 (P < 0.001) and 0.14 (P = 0.05) in the screening-eligible age and younger-age groups, respectively. Using a risk-threshold of 5-year SEER colorectal cancer incidence rate at age 50 years, adding the PRS had a similar sensitivity but improved the specificity by 11% (P < 0.001) in the screening-eligible age group. In the younger-age group it improved the sensitivity by 27% (P = 0.04) with similar specificity. CONCLUSIONS: The proposed PRS-enhanced model provides a well-calibrated 5-year colorectal cancer risk prediction and improves discriminatory accuracy in the external cohort. IMPACT: The proposed model has potential utility in risk-stratified colorectal cancer prevention.

Molecular Characteristics of Early-Onset Colorectal Cancer According to Detailed Anatomical Locations: Comparison With Later-Onset Cases.

INTRODUCTION: Early-onset colorectal cancer diagnosed before the age of 50 years has been increasing. Likely reflecting the pathogenic role of the intestinal microbiome, which gradually changes across the entire colorectal length, the prevalence of certain tumor molecular characteristics gradually changes along colorectal subsites. Understanding how colorectal tumor molecular features differ by age and tumor location is important in personalized patient management. METHODS: Using 14,004 cases with colorectal cancer including 3,089 early-onset cases, we examined microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF mutations in carcinomas of the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum and compared early-onset cases with later-onset cases. RESULTS: The proportions of MSI-high, CIMP-high, and BRAF -mutated early-onset tumors were lowest in the rectum (8.8%, 3.4%, and 3.5%, respectively) and highest in the ascending colon (46% MSI-high; 15% CIMP-high) or transverse colon (8.6% BRAF -mutated) (all Ptrend <0.001 across the rectum to ascending colon). Compared with later-onset tumors, early-onset tumors showed a higher prevalence of MSI-high status and a lower prevalence of CIMP-high status and BRAF mutations in most subsites. KRAS mutation prevalence was higher in the cecum compared with that in the other subsites in both early-onset and later-onset tumors ( P < 0.001). Notably, later-onset MSI-high tumors showed a continuous decrease in KRAS mutation prevalence from the rectum (36%) to ascending colon (9%; Ptrend <0.001), followed by an increase in the cecum (14%), while early-onset MSI-high cancers showed no such trend. DISCUSSION: Our findings support biogeographical and pathogenic heterogeneity of colorectal carcinomas in different colorectal subsites and age groups.

Prognostic role of detailed colorectal location and tumor molecular features: analyses of 13,101 colorectal cancer patients including 2994 early-onset cases.

BACKGROUND: The pathogenic effect of colorectal tumor molecular features may be influenced by several factors, including those related to microbiota, inflammation, metabolism, and epigenetics, which may change along colorectal segments. We hypothesized that the prognostic association of colon cancer location might differ by tumor molecular characteristics. METHODS: Utilizing a consortium dataset of 13,101 colorectal cancer cases, including 2994 early-onset cases, we conducted survival analyses of detailed tumor location stratified by statuses of microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF oncogenic mutation. RESULTS: There was a statistically significant trend for better colon cancer-specific survival in relation to tumor location from the cecum to sigmoid colon (Ptrend = 0.002), excluding the rectum. The prognostic association of colon location differed by MSI status (Pinteraction = 0.001). Non-MSI-high tumors exhibited the cecum-to-sigmoid trend for better colon cancer-specific survival [Ptrend < 0.001; multivariable hazard ratio (HR) for the sigmoid colon (vs. cecum), 0.80; 95% confidence interval (CI) 0.70-0.92], whereas MSI-high tumors demonstrated a suggestive cecum-to-sigmoid trend for worse survival (Ptrend = 0.020; the corresponding HR, 2.13; 95% CI 1.15-3.92). The prognostic association of colon tumor location also differed by CIMP status (Pinteraction = 0.003) but not significantly by age, stage, or other features. Furthermore, MSI-high status was a favorable prognostic indicator in all stages. CONCLUSIONS: Both detailed colonic location and tumor molecular features need to be accounted for colon cancer prognostication to advance precision medicine. Our study indicates the important role of large-scale studies to robustly examine detailed colonic subsites in molecular oncology research.