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.

Prevalence and Correlates of Post-Diagnosis Alcohol Use among Cancer Survivors

BACKGROUND

Alcohol is a risk factor for cancer and may pose unique risks for cancer survivors. Population-based studies of confirmed cancer cases are needed to estimate the extent of drinking among cancer survivors and to understand which survivors are most at risk of alcohol-related health problems.

METHODS

Cancer survivors who resided in the Puget Sound Surveillance, Epidemiology, and End Results (SEER) region, were ages 21 to 74 years at diagnosis, and were 6 to 17 months post-diagnosis at the start of the recruitment period (April 2020-December 2020) were sent a survey that included demographics, substance use, mental health, and cancer-related items. Data from returned surveys (n = 1,488) were weighted to represent the characteristics of the Puget Sound SEER region. We estimated the prevalence of post-diagnosis alcohol use as well as demographic, behavioral, and clinical correlates of three levels of drinking: any drinking, drinking exceeding cancer prevention guidelines, and hazardous drinking.

RESULTS

The weighted prevalence of any drinking, drinking exceeding cancer prevention guidelines, and hazardous drinking was 71%, 46.2%, and 31.6%, respectively. Higher income and cannabis use were associated with increased odds of all three drinking levels. Lower physical health quality of life, having non-colorectal gastrointestinal cancer, and receiving chemotherapy within the last month were associated with decreased odds of all three drinking levels.

CONCLUSIONS

The prevalence of any drinking and at-risk drinking was higher than in previous studies and differed based on sociodemographic, substance use, and cancer-related factors.

IMPACT

Findings highlight the importance of identifying and addressing risky alcohol use in cancer care settings.

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.

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.

Cannabis use after a cancer diagnosis in a population-based sample of cancer survivors

PURPOSE

This study aimed to characterize the prevalence and correlates of cannabis use and the methods and reasons for use among recently diagnosed cancer survivors in a population sample within Washington state.

METHODS

We identified individuals diagnosed with invasive cancers in the prior 6 to 17 months from April 2020 to December 2020 using the Seattle-Puget Sound Surveillance, Epidemiology, and End Results (SEER) cancer registry. Participants (n = 1,515) completed a questionnaire, including demographics, medical history, cannabis use, and other substance use. Cancer characteristics and date of diagnosis were obtained from SEER registry data. We calculated weighted prevalence estimates and logistic regression models to evaluate correlates of cannabis use.

RESULTS

Overall, 41.3% of survivors reported cannabis use at any time after diagnosis, most commonly via edibles (60.5%) and smoking (43.8%). The most frequently reported reasons for use were sleep (54.5%), mood, stress, anxiety, and depression (44.3%), pain (42.3%), and recreation (42.3%). Cannabis use was associated with younger age, race (White vs. Asian), less education, former or current smoking, consuming more than 2 alcohol-containing drinks per day, having late-stage cancer, and cancer site.

CONCLUSION

In this first evaluation of cannabis use in a registry-linked, population-based sample of survivors of all cancer types, based in a state where recreational and medical cannabis have been legal for a decade, approximately 2 in 5 survivors reported post-diagnosis use. Given how common cannabis use is among cancer survivors, there is a great need to understand its impact on cancer treatment outcomes and the overall health of cancer survivors.

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.

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.

Dietary Factors Differ Between Young-Onset and Older-Onset Colorectal Cancer Patients

We aimed to evaluate differences in dietary factors between young-onset (diagnosed at ages <50) and older-onset colorectal cancer (CRC). CRC patients diagnosed from 1998 to 2018 reported to the Puget Sound Surveillance, Epidemiology, and End Results registry were recruited using mail and telephone. Consented patients completed questionnaires assessing demographics, medical history, and CRC risk factors, including dietary factors. We used multi-variable logistic regression to calculate adjusted odds ratios (ORs) and 95% confidence intervals (CIs) comparing dietary intake in young-onset vs. older-onset CRC. Analyses included 1,087 young- and 2,554 older-onset CRC patients. Compared to older-onset CRC, young-onset CRC patients had lower intake of vegetables (OR for highest intake vs. lowest = 0.59 CI: 0.55, 0.64) and fruit (OR for highest intake vs. lowest = 0.94 CI: 0.88, 0.99) and higher intake of processed meat (OR for highest intake vs. lowest = 1.82 CI: 1.11, 2.99) and spicy food (OR for highest intake vs. lowest = 1.69 CI: 1.09, 2.61). There was no statistically significant difference between young- and older-onset CRC patients for red meat consumption. Dietary patterns differed between young- and older-onset CRC; young-onset CRC patients had lower intake of vegetables and fruit and higher intakes of processed meat and spicy food.

Colorectal cancer pre-diagnostic symptoms are associated with anatomic cancer site

BACKGROUND

Signs and red flag symptoms in colorectal cancer (CRC) patients who are below the recommended screening age are often overlooked, leading to delayed diagnosis and worse prognosis. This study investigates how patient pre-diagnostic symptoms are associated with anatomic site of their cancer and whether the association varies by age at CRC diagnosis.

METHODS

We ascertained CRC patients' experienced symptoms and screening through medical abstractions from an ongoing population-based study of CRC patients identified through a SEER cancer registry (N = 626). We used logistic regression to estimate odds ratios and 95% confidence intervals for the association between symptoms and CRC anatomic site. Additional analyses were stratified by age at diagnosis. Early-onset was defined as less than 50 years of age at CRC diagnosis.

RESULTS

Participants who experienced blood in stool were more likely (odds ratio (95% confidence interval)) to have rectal (vs. colon) cancer (4.37 (3.02, 6.33)), as were patients who experienced changes to stool (1.78 (1.21, 2.60)). Patients diagnosed with colon cancer were more likely to present with abdominal pain (0.30 (0.19, 0.47)), anemia (0.40 (0.21, 0.75)), other symptoms (0.33 (0.19, 0.55)) and no symptoms (0.68 (0.44, 1.04)). When stratifying by age at diagnosis, we found that the association between blood in stool and rectal tumor location was particularly pronounced for patients with early-onset CRC (6.48 (2.73, 15.41)).

CONCLUSIONS

Common pre-diagnostic red flag symptoms are associated with CRC anatomic site. These findings can inform best practices for gastroenterologist triage of care and early evaluation of CRC and are of key importance given the rise of early-onset (pre-screening age) CRC.

TRIAL REGISTRATION

Not applicable to this study and analysis.

Impact of weight loss on cancer-related proteins in serum: results from a cluster randomised controlled trial of individuals with type 2 diabetes

BACKGROUND

Type 2 diabetes is associated with higher risk of several cancer types. However, the biological intermediates driving this relationship are not fully understood. As novel interventions for treating and managing type 2 diabetes become increasingly available, whether they also disrupt the pathways leading to increased cancer risk is currently unknown. We investigated the effect of a type 2 diabetes intervention, in the form of intentional weight loss, on circulating proteins associated with cancer risk to gain insight into potential mechanisms linking type 2 diabetes and adiposity with cancer development.

METHODS

Fasting serum samples from participants with diabetes enrolled in the Diabetes Remission Clinical Trial (DiRECT) receiving the Counterweight-Plus weight-loss programme (intervention, N = 117, mean weight-loss 10 kg, 46% diabetes remission) or best-practice care by guidelines (control, N = 143, mean weight-loss 1 kg, 4% diabetes remission) were subject to proteomic analysis using the Olink Oncology-II platform (48% of participants were female; 52% male). To identify proteins which may be altered by the weight-loss intervention, the difference in protein levels between groups at baseline and 1 year was examined using linear regression. Mendelian randomization (MR) was performed to extend these results to evaluate cancer risk and elucidate possible biological mechanisms linking type 2 diabetes and cancer development. MR analyses were conducted using independent datasets, including large cancer meta-analyses, UK Biobank, and FinnGen, to estimate potential causal relationships between proteins modified during intentional weight loss and the risk of colorectal, breast, endometrial, gallbladder, liver, and pancreatic cancers.

FINDINGS

Nine proteins were modified by the intervention: glycoprotein Nmb; furin; Wnt inhibitory factor 1; toll-like receptor 3; pancreatic prohormone; erb-b2 receptor tyrosine kinase 2; hepatocyte growth factor; endothelial cell specific molecule 1 and Ret proto-oncogene (Holm corrected P-value <0.05). Mendelian randomization analyses indicated a causal relationship between predicted circulating furin and glycoprotein Nmb on breast cancer risk (odds ratio (OR) = 0.81, 95% confidence interval (CI) = 0.67-0.99, P-value = 0.03; and OR = 0.88, 95% CI = 0.78-0.99, P-value = 0.04 respectively), though these results were not supported in sensitivity analyses examining violations of MR assumptions.

INTERPRETATION

Intentional weight loss among individuals with recently diagnosed diabetes may modify levels of cancer-related proteins in serum. Further evaluation of the proteins identified in this analysis could reveal molecular pathways that mediate the effect of adiposity and type 2 diabetes on cancer risk.

FUNDING

The main sources of funding for this work were Diabetes UK, Cancer Research UK, World Cancer Research Fund, and Wellcome.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

BACKGROUND

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

METHODS

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,294 cancer cases and up to 1,238,345 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 × 10-8) cis-acting SNPs (i.e., in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") <0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) >70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes. Findings were replicated in the FinnGen study and then pooled using meta-analysis.

FINDINGS

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR: 1.19, 95% CI: 1.10-1.29, q-value = 0.033, PPH4 = 84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR: 1.42, 95% CI: 1.20-1.69, q-value = 0.055, PPH4 = 73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR: 0.66, 95% CI: 0.53-0.81, q-value = 0.067, PPH4 = 81.8%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR: 0.92, 95% CI: 0.88-0.97, q-value = 0.15, PPH4 = 85.6%). These findings were replicated in pooled analyses with the FinnGen study. Though suggestive evidence was found to support an association of macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR: 2.46, 95% CI: 1.48-4.10, q-value = 0.072, PPH4 = 76.1%), this finding was not replicated when pooled with the FinnGen study. For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

INTERPRETATION

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 4 circulating inflammatory markers in risk of 4 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

FUNDING

Cancer Research UK (C68933/A28534, C18281/A29019, PPRCPJT∖100005), World Cancer Research Fund (IIG_FULL_2020_022), National Institute for Health Research (NIHR202411, BRC-1215-20011), Medical Research Council (MC_UU_00011/1, MC_UU_00011/3, MC_UU_00011/6, and MC_UU_00011/4), Academy of Finland Project 326291, European Union's Horizon 2020 grant agreement no. 848158 (EarlyCause), French National Cancer Institute (INCa SHSESP20, 2020-076), Versus Arthritis (21173, 21754, 21755), National Institutes of Health (U19 CA203654), National Cancer Institute (U19CA203654).

Observational and genetic associations between cardiorespiratory fitness and cancer: a UK Biobank and international consortia study

BACKGROUND

The association of fitness with cancer risk is not clear.

METHODS

We used Cox proportional hazards models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for risk of lung, colorectal, endometrial, breast, and prostate cancer in a subset of UK Biobank participants who completed a submaximal fitness test in 2009-12 (N = 72,572). We also investigated relationships using two-sample Mendelian randomisation (MR), odds ratios (ORs) were estimated using the inverse-variance weighted method.

RESULTS

After a median of 11 years of follow-up, 4290 cancers of interest were diagnosed. A 3.5 ml O2⋅min-1⋅kg-1 total-body mass increase in fitness (equivalent to 1 metabolic equivalent of task (MET), approximately 0.5 standard deviation (SD)) was associated with lower risks of endometrial (HR = 0.81, 95% CI: 0.73-0.89), colorectal (0.94, 0.90-0.99), and breast cancer (0.96, 0.92-0.99). In MR analyses, a 0.5 SD increase in genetically predicted O2⋅min-1⋅kg-1 fat-free mass was associated with a lower risk of breast cancer (OR = 0.92, 95% CI: 0.86-0.98). After adjusting for adiposity, both the observational and genetic associations were attenuated.

DISCUSSION

Higher fitness levels may reduce risks of endometrial, colorectal, and breast cancer, though relationships with adiposity are complex and may mediate these relationships. Increasing fitness, including via changes in body composition, may be an effective strategy for cancer prevention.

Association of financial worry with substance use, mental health, and quality of life in cancer patients

PURPOSE

Financial worry is an under-appreciated negative effect of cancer. The relationship of financial worry and health behaviors in cancer is poorly characterized and has important clinical implications. This study examined the association of financial worry with substance misuse, mood, and quality of life.

METHODS

People with cancer (n = 1473; 6 to 20 months after diagnosis) were recruited from a SEER cancer registry in the Pacific Northwest. Participants completed an online survey assessing financial worry; misuse of cannabis, alcohol, and prescription drugs; tobacco smoking status; quality of life (physical and mental dimensions); anxiety; and depression. Multivariable regressions tested the association of financial worry to each health indicator and outcome.

RESULTS

In adjusted analyses, financial worry was associated with being a current vs. never smoker (odds ratio (OR) = 1.91, 95% confidence interval (CI): 1.01, 3.60), and a positive screen for an anxiety (OR = 3.01, 95% CI: 1.93, 4.68) and depressive (OR = 3.08, 95% CI: 1.89, 5.00) disorder. Financial worry was not associated with cannabis, alcohol, or prescription drug misuse (all ps > 0.05), but was associated with a decrease in physical (β =  - 2.97, 95% CI: - 4.15, - 1.79) and mental (β =  - 5.27, 95% CI: - 6.59, - 3.96) quality of life.

CONCLUSION

Financial worry among cancer survivors is associated with anxiety, depression, and worse quality of life. Of the evaluated substances, there was only an increased odds of current tobacco use with financial worry. Future longitudinal studies should inform the relationships between these factors.

IMPLICATIONS FOR CANCER SURVIVORS

Financial worry and material hardship may both need to be addressed in cancer survivorship.

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.

Combining Asian and European genome-wide association studies of colorectal cancer improves risk prediction across racial and ethnic populations

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expand PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS are 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1681-3651 cases and 8696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They are significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values < 0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

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.

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.

Association between circulating inflammatory markers and adult cancer risk: a Mendelian randomization analysis

Background

Tumour-promoting inflammation is a "hallmark" of cancer and conventional epidemiological studies have reported links between various inflammatory markers and cancer risk. The causal nature of these relationships and, thus, the suitability of these markers as intervention targets for cancer prevention is unclear.

Methods

We meta-analysed 6 genome-wide association studies of circulating inflammatory markers comprising 59,969 participants of European ancestry. We then used combined cis-Mendelian randomization and colocalisation analysis to evaluate the causal role of 66 circulating inflammatory markers in risk of 30 adult cancers in 338,162 cancer cases and up to 824,556 controls. Genetic instruments for inflammatory markers were constructed using genome-wide significant (P < 5.0 x 10-8) cis-acting SNPs (i.e. in or ±250 kb from the gene encoding the relevant protein) in weak linkage disequilibrium (LD, r2 < 0.10). Effect estimates were generated using inverse-variance weighted random-effects models and standard errors were inflated to account for weak LD between variants with reference to the 1000 Genomes Phase 3 CEU panel. A false discovery rate (FDR)-corrected P-value ("q-value") < 0.05 was used as a threshold to define "strong evidence" to support associations and 0.05 ≤ q-value < 0.20 to define "suggestive evidence". A colocalisation posterior probability (PPH4) > 70% was employed to indicate support for shared causal variants across inflammatory markers and cancer outcomes.

Results

We found strong evidence to support an association of genetically-proxied circulating pro-adrenomedullin concentrations with increased breast cancer risk (OR 1.19, 95% CI 1.10-1.29, q-value=0.033, PPH4=84.3%) and suggestive evidence to support associations of interleukin-23 receptor concentrations with increased pancreatic cancer risk (OR 1.42, 95% CI 1.20-1.69, q-value=0.055, PPH4=73.9%), prothrombin concentrations with decreased basal cell carcinoma risk (OR 0.66, 95% CI 0.53-0.81, q-value=0.067, PPH4=81.8%), macrophage migration inhibitory factor concentrations with increased bladder cancer risk (OR 1.14, 95% CI 1.05-1.23, q-value=0.072, PPH4=76.1%), and interleukin-1 receptor-like 1 concentrations with decreased triple-negative breast cancer risk (OR 0.92, 95% CI 0.88-0.97, q-value=0.15), PPH4=85.6%). For 22 of 30 cancer outcomes examined, there was little evidence (q-value ≥ 0.20) that any of the 66 circulating inflammatory markers examined were associated with cancer risk.

Conclusion

Our comprehensive joint Mendelian randomization and colocalisation analysis of the role of circulating inflammatory markers in cancer risk identified potential roles for 5 circulating inflammatory markers in risk of 5 site-specific cancers. Contrary to reports from some prior conventional epidemiological studies, we found little evidence of association of circulating inflammatory markers with the majority of site-specific cancers evaluated.

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.

Mortality by age, gene and gender in carriers of pathogenic mismatch repair gene variants receiving surveillance for early cancer diagnosis and treatment: a report from the prospective Lynch syndrome database

Background

The Prospective Lynch Syndrome Database (PLSD) collates information on carriers of pathogenic or likely pathogenic MMR variants (path_MMR) who are receiving medical follow-up, including colonoscopy surveillance, which aims to the achieve early diagnosis and treatment of cancers. Here we use the most recent PLSD cohort that is larger and has wider geographical representation than previous versions, allowing us to present mortality as an outcome, and median ages at cancer diagnoses for the first time.

Methods

The PLSD is a prospective observational study without a control group that was designed in 2012 and updated up to October 2022. Data for 8500 carriers of path_MMR variants from 25 countries were included, providing 71,713 years of follow up. Cumulative cancer incidences at 65 years of age were combined with 10-year crude survival following cancer, to derive estimates of mortality up to 75 years of age by organ, gene, and gender.

Findings

Gynaecological cancers were more frequent than colorectal cancers in path_MSH2, path_MSH6 and path_PMS2 carriers [cumulative incidence: 53.3%, 49.6% and 23.3% at 75 years, respectively]. Endometrial, colon and ovarian cancer had low mortality [8%, 13% and 15%, respectively] and prostate cancers were frequent in male path_MSH2 carriers [cumulative incidence: 39.7% at 75 years]. Pancreatic, brain, biliary tract and ureter and kidney and urinary bladder cancers were associated with high mortality [83%, 66%, 58%, 27%, and 29%, respectively]. Among path_MMR carriers undergoing colonoscopy surveillance, particularly path_MSH2 carriers, more deaths followed non-colorectal Lynch syndrome cancers than colorectal cancers.

Interpretation

In path_MMR carriers undergoing colonoscopy surveillance, non-colorectal Lynch syndrome cancers were associated with more deaths than were colorectal cancers. Reducing deaths from non-colorectal cancers presents a key challenge in contemporary medical care in Lynch syndrome.

Funding

We acknowledge funding from the Norwegian Cancer Society, contract 194751-2017.

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.

Genome-wide Interaction Study with Smoking for Colorectal Cancer Risk Identifies Novel Genetic Loci Related to Tumor Suppression, Inflammation, and Immune Response

BACKGROUND

Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer.

METHODS

A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia.

RESULTS

Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10-8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20-1.30] compared with the other genotypes (OR <1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10-8) and 8q24.23 (rs7005722, P = 2.88 × 10-8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09-1.16) compared with the other genotypes (OR <1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07-1.28) compared with the other genotypes (OR <1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33).

CONCLUSIONS

Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response.

IMPACT

These findings can guide potential prevention treatments.

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.

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.

Combining Asian-European Genome-Wide Association Studies of Colorectal Cancer Improves Risk Prediction Across Race and Ethnicity

Polygenic risk scores (PRS) have great potential to guide precision colorectal cancer (CRC) prevention by identifying those at higher risk to undertake targeted screening. However, current PRS using European ancestry data have sub-optimal performance in non-European ancestry populations, limiting their utility among these populations. Towards addressing this deficiency, we expanded PRS development for CRC by incorporating Asian ancestry data (21,731 cases; 47,444 controls) into European ancestry training datasets (78,473 cases; 107,143 controls). The AUC estimates (95% CI) of PRS were 0.63(0.62-0.64), 0.59(0.57-0.61), 0.62(0.60-0.63), and 0.65(0.63-0.66) in independent datasets including 1,681-3,651 cases and 8,696-115,105 controls of Asian, Black/African American, Latinx/Hispanic, and non-Hispanic White, respectively. They were significantly better than the European-centric PRS in all four major US racial and ethnic groups (p-values<0.05). Further inclusion of non-European ancestry populations, especially Black/African American and Latinx/Hispanic, is needed to improve the risk prediction and enhance equity in applying PRS in clinical practice.

DNA methylation profile in CpG-depleted regions uncovers a high-risk subtype of early-stage colorectal cancer

BACKGROUND

The current risk stratification system defined by clinicopathological features does not identify the risk of recurrence in early-stage (stage I-II) colorectal cancer (CRC) with sufficient accuracy. We aimed to investigate whether DNA methylation could serve as a novel biomarker for predicting prognosis in early-stage CRC patients.

METHODS

We analyzed the genome-wide methylation status of CpG loci using Infinium MethylationEPIC array run on primary tumor tissues and normal mucosa of early-stage CRC patients to identify potential methylation markers for prognosis. The machine-learning approach was applied to construct a DNA methylation-based prognostic classifier for early-stage CRC (MePEC) using the 4 gene methylation markers FAT3, KAZN, TLE4, and DUSP3. The prognostic value of the classifier was evaluated in 2 independent cohorts (n = 438 and 359, respectively).

RESULTS

The comprehensive analysis identified an epigenetic subtype with high risk of recurrence based on a group of CpG loci in the CpG-depleted region. In multivariable analysis, the MePEC classifier was independently and statistically significantly associated with time to recurrence in validation cohort 1 (hazard ratio = 2.35, 95% confidence interval = 1.47 to 3.76, P < .001) and cohort 2 (hazard ratio = 3.20, 95% confidence interval = 1.92 to 5.33, P < .001). All results were further confirmed after each cohort was stratified by clinicopathological variables and molecular subtypes.

CONCLUSIONS

We demonstrated the prognostic statistical significance of a DNA methylation profile in the CpG-depleted region, which may serve as a valuable source for tumor biomarkers. MePEC could identify an epigenetic subtype with high risk of recurrence and improve the prognostic accuracy of current clinical variables in early-stage CRC.

Separating the effects of early and later life adiposity on colorectal cancer risk: a Mendelian randomization study

BACKGROUND

Observational studies have linked childhood obesity with elevated risk of colorectal cancer; however, it is unclear if this association is causal or independent from the effects of obesity in adulthood on colorectal cancer risk.

METHODS

We conducted Mendelian randomization (MR) analyses to investigate potential causal relationships between self-perceived body size (thinner, plumper, or about average) in early life (age 10) and measured body mass index in adulthood (mean age 56.5) with risk of colorectal cancer. The total and independent effects of body size exposures were estimated using univariable and multivariable MR, respectively. Summary data were obtained from a genome-wide association study of 453,169 participants in UK Biobank for body size and from a genome-wide association study meta-analysis of three colorectal cancer consortia of 125,478 participants.

RESULTS

Genetically predicted early life body size was estimated to increase odds of colorectal cancer (odds ratio [OR] per category change: 1.12, 95% confidence interval [CI]: 0.98-1.27), with stronger results for colon cancer (OR: 1.16, 95% CI: 1.00-1.35), and distal colon cancer (OR: 1.25, 95% CI: 1.04-1.51). After accounting for adult body size using multivariable MR, effect estimates for early life body size were attenuated towards the null for colorectal cancer (OR: 0.97, 95% CI: 0.77-1.22) and colon cancer (OR: 0.97, 95% CI: 0.76-1.25), while the estimate for distal colon cancer was of similar magnitude but more imprecise (OR: 1.27, 95% CI: 0.90-1.77). Genetically predicted adult life body size was estimated to increase odds of colorectal (OR: 1.27, 95% CI: 1.03, 1.57), colon (OR: 1.32, 95% CI: 1.05, 1.67), and proximal colon (OR: 1.57, 95% CI: 1.21, 2.05).

CONCLUSIONS

Our findings suggest that the positive association between early life body size and colorectal cancer risk is likely due to large body size retainment into adulthood.

Quantitative Pathologic Analysis of Digitized Images of Colorectal Carcinoma Improves Prediction of Recurrence-Free Survival

BACKGROUND & AIMS

To examine whether quantitative pathologic analysis of digitized hematoxylin and eosin slides of colorectal carcinoma (CRC) correlates with clinicopathologic features, molecular alterations, and prognosis.

METHODS

A quantitative segmentation algorithm (QuantCRC) was applied to 6468 digitized hematoxylin and eosin slides of CRCs. Fifteen parameters were recorded from each image and tested for associations with clinicopathologic features and molecular alterations. A prognostic model was developed to predict recurrence-free survival using data from the internal cohort (n = 1928) and validated on an internal test (n = 483) and external cohort (n = 938).

RESULTS

There were significant differences in QuantCRC according to stage, histologic subtype, grade, venous/lymphatic/perineural invasion, tumor budding, CD8 immunohistochemistry, mismatch repair status, KRAS mutation, BRAF mutation, and CpG methylation. A prognostic model incorporating stage, mismatch repair, and QuantCRC resulted in a Harrell's concordance (c)-index of 0.714 (95% confidence interval [CI], 0.702-0.724) in the internal test and 0.744 (95% CI, 0.741-0.754) in the external cohort. Removing QuantCRC from the model reduced the c-index to 0.679 (95% CI, 0.673-0.694) in the external cohort. Prognostic risk groups were identified, which provided a hazard ratio of 2.24 (95% CI, 1.33-3.87, P = .004) for low vs high-risk stage III CRCs and 2.36 (95% CI, 1.07-5.20, P = .03) for low vs high-risk stage II CRCs, in the external cohort after adjusting for established risk factors. The predicted median 36-month recurrence rate for high-risk stage III CRCs was 32.7% vs 13.4% for low-risk stage III and 15.8% for high-risk stage II vs 5.4% for low-risk stage II CRCs.

CONCLUSIONS

QuantCRC provides a powerful adjunct to routine pathologic reporting of CRC. A prognostic model using QuantCRC improves prediction of recurrence-free survival.

Abstract A010: Molecular characteristics and prognosis of early-onset colorectal cancer according to detailed anatomical locations: Comparison to later-onset cases

Background: Evidence suggests that, compared to later-onset colorectal cancer (CRC), early-onset CRC (diagnosed before 50 years of age) occur more frequently in rectal location and less frequently in proximal colon location. Studies also indicate heterogeneity of molecular characteristics within early-onset CRCs and between early-onset and later-onset CRCs. Understanding how tumor molecular features and biological aggressiveness of early-onset CRC differ by detailed tumor location is important in personalized patient management. Methods: Using 14,004 CRC cases including 3,089 early-onset cases from The Cancer Genome Atlas and the Genetics and Epidemiology of CRC Consortium (GECCO), we examined microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF mutations in cancers of the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum. We further conducted CRC-specific survival analyses of detailed tumor location stratified by molecular characteristics using multivariable Cox proportional hazards models. Results: The proportions of MSI-high, CIMP-high, and BRAF-mutated early-onset tumors gradually increased from the rectum (8.8%, 3.4%, and 3.5%, respectively) to the ascending colon (46% MSI-high; 15% CIMP-high) or transverse colon (8.6% BRAF-mutated), followed by declines in the cecum for MSI-high (36%) and BRAF-mutated (4.6%) tumors. Notably, later-onset MSI-high tumors showed a continuous decrease in KRAS mutation prevalence from the rectum (36%) to ascending colon (9%; Ptrend &lt;0.001) followed by an increase in the cecum (14%), while early-onset MSI-high cancer showed no such trend. In survival analyses, there was a trend for better CRC-specific survival from the cecum to sigmoid colon (Ptrend&lt;0.001) (but not to rectum), which did not significantly differ by age of diagnosis (Pinteraction=0.61). Multivariable-adjusted hazard ratios (HRs) for sigmoid vs. cecum were 0.86 (95% CI, 0.64-1.14) for early-onset cases and 0.81 (95% CI, 0.70-0.93) for later-onset cases. Notably, the prognostic association of tumor location differed by MSI status (Pinteraction=0.001). Non-MSI-high tumors exhibited the cecum-to-sigmoid trend for better CRC-specific survival [Ptrend &lt;0.001; multivariable-adjusted HR for the sigmoid colon (vs. cecum), 0.80; 95% CI, 0.70-0.92], whereas MSI-high tumors demonstrated an opposite cecum-to-sigmoid trend for worse survival (Ptrend = 0.020; the corresponding HR, 2.13; 95% CI, 1.15-3.92). Conclusions: Our findings support biogeographical and pathogenic heterogeneity of CRC in different age groups and colorectal subsites. Both detailed colorectal location and tumor molecular features need to be accounted for prognostication of early-onset and later-onset CRCs to advance precision medicine.

Citation Format: Tomotaka Ugai, Yin Cao, Koichiro Haruki, Tabitha A. Harrison, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Steven Gallinger, Marc J. Gunter, Michael Hoffmeister, Mark A. Jenkins, Roger L. Milne, Polly A. Newcomb, Rish K. Pai, Andrew J. Pellatt, Robert E. Schoen, Bethany Van Guelpen, Michael O. Woods, Amanda I. Phipps, Ulrike Peters, Shuji Ogino. Molecular characteristics and prognosis of early-onset colorectal cancer according to detailed anatomical locations: Comparison to later-onset cases [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer; 2022 Oct 1-4; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_1):Abstract nr A010.

T cell-inflamed gene expression profile is associated with favorable disease-specific survival in non-hypermutated microsatellite-stable colorectal cancer patients

BACKGROUND

The anti-tumor immune response plays a key role in colorectal cancer (CRC) progression and survival. The T cell-inflamed gene expression profile (GEP) is a biomarker predicting response to checkpoint inhibitor immunotherapy across immunogenic cancer types, but the prognostic value in CRC is unknown. We evaluated associations with disease-specific survival, somatic mutations, and examined its differentially expressed genes and pathways among 84 sporadic CRC patients from the Seattle Colon Cancer Family Registry.

METHODS

Gene expression profiling was performed using Nanostring's nCounter PanCancer IO 360 panel. Somatic mutations were identified by a targeted DNA sequencing panel.

RESULTS

The T cell-inflamed GEP was positively associated with tumor mutation burden and microsatellite instability high (MSI-H). Higher T cell-inflamed GEP had favorable CRC-specific survival (hazard ratio [HR] per standard deviation unit = 0.50, p = 0.004) regardless of hypermutation or MSI status. Analysis of recurrently mutated genes having at least 10 mutation carriers, suggested that the T cell-inflamed GEP is positively associated with RYR1, and negatively associated with APC. However, these associations were attenuated after adjusting for hypermutation or MSI status. We also found that expression of genes RPL23, EPCAM, AREG and ITGA6, and the Wnt signaling pathway was negatively associated with the T cell-inflamed GEP, which might indicate immune-inhibitory mechanisms.

CONCLUSIONS

Our results show that the T cell-inflamed GEP is a prognostic biomarker in non-hypermutated microsatellite-stable CRC. This also suggests that patient stratification for immunotherapy within this CRC subgroup should be explored further. Moreover, reported immune-inhibitory gene expression signals may suggest targets for therapeutic combination with immunotherapy.

Interactions between folate intake and genetic predictors of gene expression levels associated with colorectal cancer risk

Observational studies have shown higher folate consumption to be associated with lower risk of colorectal cancer (CRC). Understanding whether and how genetic risk factors interact with folate could further elucidate the underlying mechanism. Aggregating functionally relevant genetic variants in set-based variant testing has higher power to detect gene-environment (G × E) interactions and may provide information on the underlying biological pathway. We investigated interactions between folate consumption and predicted gene expression on colorectal cancer risk across the genome. We used variant weights from the PrediXcan models of colon tissue-specific gene expression as a priori variant information for a set-based G × E approach. We harmonized total folate intake (mcg/day) based on dietary intake and supplemental use across cohort and case-control studies and calculated sex and study specific quantiles. Analyses were performed using a mixed effects score tests for interactions between folate and genetically predicted expression of 4839 genes with available genetically predicted expression. We pooled results across 23 studies for a total of 13,498 cases with colorectal tumors and 13,918 controls of European ancestry. We used a false discovery rate of 0.2 to identify genes with suggestive evidence of an interaction. We found suggestive evidence of interaction with folate intake on CRC risk for genes including glutathione S-Transferase Alpha 1 (GSTA1; p = 4.3E-4), Tonsuko Like, DNA Repair Protein (TONSL; p = 4.3E-4), and Aspartylglucosaminidase (AGA: p = 4.5E-4). We identified three genes involved in preventing or repairing DNA damage that may interact with folate consumption to alter CRC risk. Glutathione is an antioxidant, preventing cellular damage and is a downstream metabolite of homocysteine and metabolized by GSTA1. TONSL is part of a complex that functions in the recovery of double strand breaks and AGA plays a role in lysosomal breakdown of glycoprotein.

Predictive Polygenic Score for Outcome after First-Line Oxaliplatin-Based Chemotherapy in Colorectal Cancer Patients Using Supervised Principal Component Analysis

BACKGROUND

Associations between candidate germline genetic variants and treatment outcome of oxaliplatin, a drug commonly used for patients with colorectal cancer, have been reported but not robustly established. This study aimed to construct polygenic hazard scores (PHSs) as predictive markers for oxaliplatin treatment outcome by using a supervised principal component approach (PCA).

METHODS

Genome-wide association analysis for overall survival, including interaction terms (SNP*treatment type) was carried out using two phase III trials, 3,098 resected stage III colon cancer (rCC) patients of NCCTG N0147 and 506 metastatic colorectal cancer (mCRC) patients of NCCTG N9741, separately. SNPs showing interaction with genome-wide significance (P < 5 × 10-8) were selected for PCA to derive a PHS. PHS interaction with treatment was included in Cox regression models to predict outcome. Replication of prediction models was performed in an independent cohort, DACHS.

RESULTS

The two PHSs based on the first two principal components of selected SNPs (15SNPs for rCC and 13SNPs for mCRC) were used to construct interaction terms with treatment type and included in models adjusted for clinical covariables. However, in the DACHS study, the addition of the two PHS terms to clinical models did not improve the prediction error in either patients with rCC or mCRC. PHS interaction was also not replicated.

CONCLUSIONS

The PHSs derived using principal components efficiently combined multiple predictive SNPs for estimating likelihood of benefit from oxaliplatin versus other treatment but could not be replicated.

IMPACT

These results highlight the potential but also challenges in generating evidence for a predictive polygenic score for oxaliplatin efficacy.

Colorectal cancer incidences in Lynch syndrome: a comparison of results from the prospective lynch syndrome database and the international mismatch repair consortium

OBJECTIVE

To compare colorectal cancer (CRC) incidences in carriers of pathogenic variants of the MMR genes in the PLSD and IMRC cohorts, of which only the former included mandatory colonoscopy surveillance for all participants.

METHODS

CRC incidences were calculated in an intervention group comprising a cohort of confirmed carriers of pathogenic or likely pathogenic variants in mismatch repair genes (path_MMR) followed prospectively by the Prospective Lynch Syndrome Database (PLSD). All had colonoscopy surveillance, with polypectomy when polyps were identified. Comparison was made with a retrospective cohort reported by the International Mismatch Repair Consortium (IMRC). This comprised confirmed and inferred path_MMR carriers who were first- or second-degree relatives of Lynch syndrome probands.

RESULTS

In the PLSD, 8,153 subjects had follow-up colonoscopy surveillance for a total of 67,604 years and 578 carriers had CRC diagnosed. Average cumulative incidences of CRC in path_MLH1 carriers at 70 years of age were 52% in males and 41% in females; for path_MSH2 50% and 39%; for path_MSH6 13% and 17% and for path_PMS2 11% and 8%. In contrast, in the IMRC cohort, corresponding cumulative incidences were 40% and 27%; 34% and 23%; 16% and 8% and 7% and 6%. Comparing just the European carriers in the two series gave similar findings. Numbers in the PLSD series did not allow comparisons of carriers from other continents separately. Cumulative incidences at 25 years were < 1% in all retrospective groups.

CONCLUSIONS

Prospectively observed CRC incidences (PLSD) in path_MLH1 and path_MSH2 carriers undergoing colonoscopy surveillance and polypectomy were higher than in the retrospective (IMRC) series, and were not reduced in path_MSH6 carriers. These findings were the opposite to those expected. CRC point incidence before 50 years of age was reduced in path_PMS2 carriers subjected to colonoscopy, but not significantly so.

Genome-Wide Interaction Analysis of Genetic Variants With Menopausal Hormone Therapy for Colorectal Cancer Risk

BACKGROUND

The use of menopausal hormone therapy (MHT) may interact with genetic variants to influence colorectal cancer (CRC) risk.

METHODS

We conducted a genome-wide, gene-environment interaction between single nucleotide polymorphisms and the use of any MHT, estrogen only, and combined estrogen-progestogen therapy with CRC risk, among 28 486 postmenopausal women (11 519 CRC patients and 16 967 participants without CRC) from 38 studies, using logistic regression, 2-step method, and 2- or 3-degree-of-freedom joint test. A set-based score test was applied for rare genetic variants.

RESULTS

The use of any MHT, estrogen only and estrogen-progestogen were associated with a reduced CRC risk (odds ratio [OR] = 0.71, 95% confidence interval [CI] = 0.64 to 0.78; OR = 0.65, 95% CI = 0.53 to 0.79; and OR = 0.73, 95% CI = 0.59 to 0.90, respectively). The 2-step method identified a statistically significant interaction between a GRIN2B variant rs117868593 and MHT use, whereby MHT-associated CRC risk was statistically significantly reduced in women with the GG genotype (OR = 0.68, 95% CI = 0.64 to 0.72) but not within strata of GC or CC genotypes. A statistically significant interaction between a DCBLD1 intronic variant at 6q22.1 (rs10782186) and MHT use was identified by the 2-degree-of-freedom joint test. The MHT-associated CRC risk was reduced with increasing number of rs10782186-C alleles, showing odds ratios of 0.78 (95% CI = 0.70 to 0.87) for TT, 0.68 (95% CI = 0.63 to 0.73) for TC, and 0.66 (95% CI = 0.60 to 0.74) for CC genotypes. In addition, 5 genes in rare variant analysis showed suggestive interactions with MHT (2-sided P < 1.2 × 10-4).

CONCLUSION

Genetic variants that modify the association between MHT and CRC risk were identified, offering new insights into pathways of CRC carcinogenesis and potential mechanisms involved.

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.

Type-2 Diabetes Mellitus and Risk of Colorectal Polyps: A Colonoscopy-Based Study Using Natural Language Processing

PURPOSE

Although type-2 diabetes (T2D) has been associated with colorectal cancer in previous studies, the association of T2D with colorectal polyps is unknown.

METHODS

Using pathology reports from the University of Utah (UU) Enterprise Data Warehouse (EDW), we developed a rule-based natural language processing (NLP) pipeline to extract colorectal polyp diagnoses and features (site, shape, number, size) on 15,679 patients who underwent a colonoscopy at the UU Gastroenterology clinic from 2013-2016. The NLP pipeline was validated by manual abstraction of 350 pathology reports, and demonstrated excellent performance (accuracy 91%). Patient characteristics, including age, sex, race, diabetes status, smoking, BMI, and medication use, were abstracted from the EDW. Odds ratios (OR) and 95% confidence limits (95% CI) adjusted for abstracted variables were calculated using multivariable polytomous logistic regression.

RESULTS

Participants were on average 56 years old, 85% White, 50% male, with a mean BMI of 29 kg/m2. About 27% of the participants reported history of T2D; 71% of whom used anti-diabetes medication. Participants were classified as having adenomas (30%), serrated polyps (16%), synchronous adenomas and serrated polyps (19%) or as polyp-free controls (35%). T2D was associated with a statistically significant lower risk of colorectal polyps [0.83(0.73,0.92)]. When evaluated by polyp subtype, T2D was marginally associated with reduced adenoma risk [0.90(0.80,1.02)], and inversely associated with risk of serrated polyps [0.80(0.67,0.93)]. The associations did not vary by lesion severity within polyp subtypes. There was a statistically significant deceased risk for polyps among anti-diabetes medication users [0.84(0.69,0.99)].

CONCLUSIONS

Overall, T2D was associated with a statistically significant reduced risk of colorectal polyps; this reduced risk was consistent for both adenomas and serrated polyps. As T2D has previously been shown to increase colorectal cancer risk, this differential association with colorectal polyps may possibly be due to a variable effect of anti-diabetes medication use. Further studies are needed to better understand the mechanisms through which diabetes and its treatment may be differentially associated with colorectal polyps.

Association between germline variants and somatic mutations in colorectal cancer

Colorectal cancer (CRC) is a heterogeneous disease with evidence of distinct tumor types that develop through different somatically altered pathways. To better understand the impact of the host genome on somatically mutated genes and pathways, we assessed associations of germline variations with somatic events via two complementary approaches. We first analyzed the association between individual germline genetic variants and the presence of non-silent somatic mutations in genes in 1375 CRC cases with genome-wide SNPs data and a tumor sequencing panel targeting 205 genes. In the second analysis, we tested if germline variants located within previously identified regions of somatic allelic imbalance were associated with overall CRC risk using summary statistics from a recent large scale GWAS (n≃125 k CRC cases and controls). The first analysis revealed that a variant (rs78963230) located within a CNA region associated with TLR3 was also associated with a non-silent mutation within gene FBXW7. In the secondary analysis, the variant rs2302274 located in CDX1/PDGFRB frequently gained/lost in colorectal tumors was associated with overall CRC risk (OR = 0.96, p = 7.50e-7). In summary, we demonstrate that an integrative analysis of somatic and germline variation can lead to new insights about CRC.

Large-scale Integrated Analysis of Genetics and Metabolomic Data Reveals Potential Links Between Lipids and Colorectal Cancer Risk

BACKGROUND

The etiology of colorectal cancer is not fully understood.

METHODS

Using genetic variants and metabolomics data including 217 metabolites from the Framingham Heart Study (n = 1,357), we built genetic prediction models for circulating metabolites. Models with prediction R2 > 0.01 (Nmetabolite = 58) were applied to predict levels of metabolites in two large consortia with a combined sample size of approximately 46,300 cases and 59,200 controls of European and approximately 21,700 cases and 47,400 controls of East Asian (EA) descent. Genetically predicted levels of metabolites were evaluated for their associations with colorectal cancer risk in logistic regressions within each racial group, after which the results were combined by meta-analysis.

RESULTS

Of the 58 metabolites tested, 24 metabolites were significantly associated with colorectal cancer risk [Benjamini-Hochberg FDR (BH-FDR) < 0.05] in the European population (ORs ranged from 0.91 to 1.06; P values ranged from 0.02 to 6.4 × 10-8). Twenty one of the 24 associations were replicated in the EA population (ORs ranged from 0.26 to 1.69, BH-FDR < 0.05). In addition, the genetically predicted levels of C16:0 cholesteryl ester was significantly associated with colorectal cancer risk in the EA population only (OREA: 1.94, 95% CI, 1.60-2.36, P = 2.6 × 10-11; OREUR: 1.01, 95% CI, 0.99-1.04, P = 0.3). Nineteen of the 25 metabolites were glycerophospholipids and triacylglycerols (TAG). Eighteen associations exhibited significant heterogeneity between the two racial groups (PEUR-EA-Het < 0.005), which were more strongly associated in the EA population. This integrative study suggested a potential role of lipids, especially certain glycerophospholipids and TAGs, in the etiology of colorectal cancer.

CONCLUSIONS

This study identified potential novel risk biomarkers for colorectal cancer by integrating genetics and circulating metabolomics data.

IMPACT

The identified metabolites could be developed into new tools for risk assessment of colorectal cancer in both European and EA populations.

Associations Between Glycemic Traits and Colorectal Cancer: A Mendelian Randomization Analysis

BACKGROUND

Glycemic traits-such as hyperinsulinemia, hyperglycemia, and type 2 diabetes-have been associated with higher colorectal cancer risk in observational studies; however, causality of these associations is uncertain. We used Mendelian randomization (MR) to estimate the causal effects of fasting insulin, 2-hour glucose, fasting glucose, glycated hemoglobin (HbA1c), and type 2 diabetes with colorectal cancer.

METHODS

Genome-wide association study summary data were used to identify genetic variants associated with circulating levels of fasting insulin (n = 34), 2-hour glucose (n = 13), fasting glucose (n = 70), HbA1c (n = 221), and type 2 diabetes (n = 268). Using 2-sample MR, we examined these variants in relation to colorectal cancer risk (48 214 case patient and 64 159 control patients).

RESULTS

In inverse-variance models, higher fasting insulin levels increased colorectal cancer risk (odds ratio [OR] per 1-SD = 1.65, 95% confidence interval [CI] = 1.15 to 2.36). We found no evidence of any effect of 2-hour glucose (OR per 1-SD = 1.02, 95% CI = 0.86 to 1.21) or fasting glucose (OR per 1-SD = 1.04, 95% CI = 0.88 to 1.23) concentrations on colorectal cancer risk. Genetic liability to type 2 diabetes (OR per 1-unit increase in log odds = 1.04, 95% CI = 1.01 to 1.07) and higher HbA1c levels (OR per 1-SD = 1.09, 95% CI = 1.00 to 1.19) increased colorectal cancer risk, although these findings may have been biased by pleiotropy. Higher HbA1c concentrations increased rectal cancer risk in men (OR per 1-SD = 1.21, 95% CI = 1.05 to 1.40), but not in women.

CONCLUSIONS

Our results support a causal effect of higher fasting insulin, but not glucose traits or type 2 diabetes, on increased colorectal cancer risk. This suggests that pharmacological or lifestyle interventions that lower circulating insulin levels may be beneficial in preventing colorectal tumorigenesis.

Beyond GWAS of Colorectal Cancer: Evidence of Interaction with Alcohol Consumption and Putative Causal Variant for the 10q24.2 Region

BACKGROUND

Currently known associations between common genetic variants and colorectal cancer explain less than half of its heritability of 25%. As alcohol consumption has a J-shape association with colorectal cancer risk, nondrinking and heavy drinking are both risk factors for colorectal cancer.

METHODS

Individual-level data was pooled from the Colon Cancer Family Registry, Colorectal Transdisciplinary Study, and Genetics and Epidemiology of Colorectal Cancer Consortium to compare nondrinkers (≤1 g/day) and heavy drinkers (>28 g/day) with light-to-moderate drinkers (1-28 g/day) in GxE analyses. To improve power, we implemented joint 2df and 3df tests and a novel two-step method that modifies the weighted hypothesis testing framework. We prioritized putative causal variants by predicting allelic effects using support vector machine models.

RESULTS

For nondrinking as compared with light-to-moderate drinking, the hybrid two-step approach identified 13 significant SNPs with pairwise r2 > 0.9 in the 10q24.2/COX15 region. When stratified by alcohol intake, the A allele of lead SNP rs2300985 has a dose-response increase in risk of colorectal cancer as compared with the G allele in light-to-moderate drinkers [OR for GA genotype = 1.11; 95% confidence interval (CI), 1.06-1.17; OR for AA genotype = 1.22; 95% CI, 1.14-1.31], but not in nondrinkers or heavy drinkers. Among the correlated candidate SNPs in the 10q24.2/COX15 region, rs1318920 was predicted to disrupt an HNF4 transcription factor binding motif.

CONCLUSIONS

Our study suggests that the association with colorectal cancer in 10q24.2/COX15 observed in genome-wide association study is strongest in nondrinkers. We also identified rs1318920 as the putative causal regulatory variant for the region.

IMPACT

The study identifies multifaceted evidence of a possible functional effect for rs1318920.

Genetic Regulation of DNA Methylation Yields Novel Discoveries in GWAS of Colorectal Cancer

BACKGROUND

Colorectal cancer has a strong epigenetic component that is accompanied by frequent DNA methylation (DNAm) alterations in addition to heritable genetic risk. It is of interest to understand the interrelationship of germline genetics, DNAm, and colorectal cancer risk.

METHODS

We performed a genome-wide methylation quantitative trait locus (meQTL) analysis in 1,355 people, assessing the pairwise associations between genetic variants and lymphocytes methylation data. In addition, we used penalized regression with cis-genetic variants ± 1 Mb of methylation to identify genome-wide heritable DNAm. We evaluated the association of genetically predicted methylation with colorectal cancer risk based on genome-wide association studies (GWAS) of over 125,000 cases and controls using the multivariate sMiST as well as univariately via examination of marginal association with colorectal cancer risk.

RESULTS

Of the 142 known colorectal cancer GWAS loci, 47 were identified as meQTLs. We identified four novel colorectal cancer-associated loci (NID2, ATXN10, KLHDC10, and CEP41) that reside over 1 Mb outside of known colorectal cancer loci and 10 secondary signals within 1 Mb of known loci.

CONCLUSIONS

Leveraging information of DNAm regulation into genetic association of colorectal cancer risk reveals novel pathways in colorectal cancer tumorigenesis. Our summary statistics-based framework sMiST provides a powerful approach by combining information from the effect through methylation and residual direct effects of the meQTLs on disease risk. Further validation and functional follow-up of these novel pathways are needed.

IMPACT

Using genotype, DNAm, and GWAS, we identified four new colorectal cancer risk loci. We studied the landscape of genetic regulation of DNAm via single-SNP and multi-SNP meQTL analyses.

Genetic variants associated with circulating C-reactive protein levels and colorectal cancer survival: Sex-specific and lifestyle factors specific associations

Elevated blood levels of C-reactive protein (CRP) have been linked to colorectal cancer (CRC) survival. We evaluated genetic variants associated with CRP levels and their interactions with sex and lifestyle factors in association with CRC-specific mortality. Our study included 16 142 CRC cases from the International Survival Analysis in Colorectal Cancer Consortium. We identified 618 common single nucleotide polymorphisms (SNPs) associated with CRP levels from the NHGRI-EBI GWAS Catalog. Cox proportional hazards regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for associations between SNPs and CRC-specific mortality adjusting for age, sex, genotyping platform/study and principal components. We investigated their interactions with sex and lifestyle factors using likelihood ratio tests. Of 5472 (33.9%) deaths accrued over up to 10 years of follow-up, 3547 (64.8%) were due to CRC. No variants were associated with CRC-specific mortality after multiple comparison correction. We observed strong evidence of interaction between variant rs1933736 at FRK gene and sex in relation to CRC-specific mortality (corrected Pinteraction  = .0004); women had higher CRC-specific mortality associated with the minor allele (HR = 1.11, 95% CI = 1.04-1.19) whereas an inverse association was observed for men (HR = 0.88, 95% CI = 0.82-0.94). There was no evidence of interactions between CRP-associated SNPs and alcohol, obesity or smoking. Our study observed a significant interaction between sex and a CRP-associated variant in relation to CRC-specific mortality. Future replication of this association and functional annotation of the variant are needed.

Risk Stratification for Early-Onset Colorectal Cancer Using a Combination of Genetic and Environmental Risk Scores: An International Multi-Center Study

BACKGROUND

The incidence of colorectal cancer (CRC) among individuals aged younger than 50 years has been increasing. As screening guidelines lower the recommended age of screening initiation, concerns including the burden on screening capacity and costs have been recognized, suggesting that an individualized approach may be warranted. We developed risk prediction models for early-onset CRC that incorporate an environmental risk score (ERS), including 16 lifestyle and environmental factors, and a polygenic risk score (PRS) of 141 variants.

METHODS

Relying on risk score weights for ERS and PRS derived from studies of CRC at all ages, we evaluated risks for early-onset CRC in 3486 cases and 3890 controls aged younger than 50 years. Relative and absolute risks for early-onset CRC were assessed according to values of the ERS and PRS. The discriminatory performance of these scores was estimated using the covariate-adjusted area under the receiver operating characteristic curve.

RESULTS

Increasing values of ERS and PRS were associated with increasing relative risks for early-onset CRC (odds ratio per SD of ERS = 1.14, 95% confidence interval [CI] = 1.08 to 1.20; odds ratio per SD of PRS = 1.59, 95% CI = 1.51 to 1.68), both contributing to case-control discrimination (area under the curve = 0.631, 95% CI = 0.615 to 0.647). Based on absolute risks, we can expect 26 excess cases per 10 000 men and 21 per 10 000 women among those scoring at the 90th percentile for both risk scores.

CONCLUSIONS

Personal risk scores have the potential to identify individuals at differential relative and absolute risk for early-onset CRC. Improved discrimination may aid in targeted CRC screening of younger, high-risk individuals, potentially improving outcomes.

Mortality risk and physical activity across the lifespan in endometrial cancer survivors

PURPOSE

Physical activity (pre- and post-diagnosis) has been studied in prevention and survivorship contexts for endometrial cancer. However, the association of physical activity (PA) across the lifespan on mortality risk among endometrial cancer survivors is understudied. The study's objective was to identify the association of lifetime PA on mortality risk in endometrial cancer survivors.

METHODS

Seven hundred forty-five endometrial cancer survivors drawn from a population-based cancer registry (diagnosed between 1991 and 1994) reported the frequency (sessions/week) of moderate- and vigorous intensity physical activity (MVPA) at age 12, age 20, and 5 years pre-interview (post-diagnosis). Cox proportional hazards were used to estimate hazard ratios (HR) and 95% confidence intervals for the association between PA, all-cause, and cardiovascular disease mortality as assessed in 2016. MVPA was modeled using natural cubic splines.

RESULTS

Diagnosis age, body mass index, and smoking (pack-years) were each positively associated with increased all-cause mortality risk. Those who did one session of MVPA 5 years pre-interview had a lower mortality risk (HR 0.61; 95% CI 0.41-0.92) compared to those with no MVPA. Those reporting one session of MVPA was similarly observed at age 12 (HR 0.95; 95% CI 0.86-1.06) and at age 20 (HR 0.87; 95% CI 0.65-1.16).

CONCLUSION

Those who participated in PA, compared to those who did not, in the 5 years before diagnosis had a lower mortality risk. While PA was not independently protective against mortality risk at ages 12 or 20, PA is still important for endometrial cancer survivors for other non-mortality outcomes.

Genome-wide association study identifies tumor anatomical site-specific risk variants for colorectal cancer survival

Identification of new genetic markers may improve the prediction of colorectal cancer prognosis. Our objective was to examine genome-wide associations of germline genetic variants with disease-specific survival in an analysis of 16,964 cases of colorectal cancer. We analyzed genotype and colorectal cancer-specific survival data from a consortium of 15 studies. Approximately 7.5 million SNPs were examined under the log-additive model using Cox proportional hazards models, adjusting for clinical factors and principal components. Additionally, we ran secondary analyses stratifying by tumor site and disease stage. We used a genome-wide p-value threshold of 5 × 10-8 to assess statistical significance. No variants were statistically significantly associated with disease-specific survival in the full case analysis or in the stage-stratified analyses. Three SNPs were statistically significantly associated with disease-specific survival for cases with tumors located in the distal colon (rs698022, HR = 1.48, CI 1.30-1.69, p = 8.47 × 10-9) and the proximal colon (rs189655236, HR = 2.14, 95% CI 1.65-2.77, p = 9.19 × 10-9 and rs144717887, HR = 2.01, 95% CI 1.57-2.58, p = 3.14 × 10-8), whereas no associations were detected for rectal tumors. Findings from this large genome-wide association study highlight the potential for anatomical-site-stratified genome-wide studies to identify germline genetic risk variants associated with colorectal cancer-specific survival. Larger sample sizes and further replication efforts are needed to more fully interpret these findings.

Molecular and Pathology Features of Colorectal Tumors and Patient Outcomes Are Associated with Fusobacterium nucleatum and Its Subspecies animalis

BACKGROUND

Fusobacterium nucleatum (F. nucleatum) activates oncogenic signaling pathways and induces inflammation to promote colorectal carcinogenesis.

METHODS

We characterized F. nucleatum and its subspecies in colorectal tumors and examined associations with tumor characteristics and colorectal cancer-specific survival. We conducted deep sequencing of nusA, nusG, and bacterial 16s rRNA genes in tumors from 1,994 patients with colorectal cancer and assessed associations between F. nucleatum presence and clinical characteristics, colorectal cancer-specific mortality, and somatic mutations.

RESULTS

F. nucleatum, which was present in 10.3% of tumors, was detected in a higher proportion of right-sided and advanced-stage tumors, particularly subspecies animalis. Presence of F. nucleatum was associated with higher colorectal cancer-specific mortality (HR, 1.97; P = 0.0004). This association was restricted to nonhypermutated, microsatellite-stable tumors (HR, 2.13; P = 0.0002) and those who received chemotherapy [HR, 1.92; confidence interval (CI), 1.07-3.45; P = 0.029). Only F. nucleatum subspecies animalis, the main subspecies detected (65.8%), was associated with colorectal cancer-specific mortality (HR, 2.16; P = 0.0016), subspecies vincentii and nucleatum were not (HR, 1.07; P = 0.86). Additional adjustment for tumor stage suggests that the effect of F. nucleatum on mortality is partly driven by a stage shift. Presence of F. nucleatum was associated with microsatellite instable tumors, tumors with POLE exonuclease domain mutations, and ERBB3 mutations, and suggestively associated with TP53 mutations.

CONCLUSIONS

F. nucleatum, and particularly subspecies animalis, was associated with a higher colorectal cancer-specific mortality and specific somatic mutated genes.

IMPACT

Our findings identify the F. nucleatum subspecies animalis as negatively impacting colorectal cancer mortality, which may occur through a stage shift and its effect on chemoresistance.

A genome-wide search for determinants of survival in 1926 patients with advanced colorectal cancer with follow-up in over 22,000 patients

BACKGROUND

While genome-wide association studies (GWAS) have identified germline variants influencing the risk of developing colorectal cancer (CRC), there has been limited examination of the possible role of inherited variation as a determinant of patient outcome.

PATIENTS AND METHODS

We performed a GWAS for overall survival (OS) in 1926 patients with advanced CRC from the COIN and COIN-B clinical trials. For single nucleotide polymorphisms (SNPs) showing an association with OS (P < 1.0 × 10-5), we conducted sensitivity analyses based on the time from diagnosis to death and sought independent replications in 5675 patients from the Study of Colorectal Cancer in Scotland (SOCCS) and 16,964 patients from the International Survival Analysis in Colorectal cancer Consortium (ISACC). We analysed the Human Protein Atlas to determine if ERBB4 expression was associated with survival in 438 patients with colon adenocarcinomas.

RESULTS

The most significant SNP associated with OS was rs79612564 in ERBB4 (hazard ratio [HR] = 1.24, 95% confidence interval [CI] = 1.16-1.32, P = 1.9 × 10-7). SNPs at 17 loci had suggestive associations for OS and all had similar effects on the time from diagnosis to death. No lead SNPs were independently replicated in the meta-analysis of all patients from SOCCS and ISACC. However, rs79612564 was significant in stage-IV patients from SOCCS (P = 2.1 × 10-2) but not ISACC (P = 0.89) and SOCCS combined with COIN and COIN-B attained genome-wide significance (P = 1.7 × 10-8). Patients with high ERBB4 expression in their colon adenocarcinomas had worse survival (HR = 1.50, 95% CI = 1.1-1.9, P = 4.6 × 10-2).

CONCLUSIONS

Genetic and expression data support a potential role for rs79612564 in the receptor tyrosine kinase ERBB4 as a predictive biomarker of survival.

A Statistical Method for Association Analysis of Cell Type Compositions

Gene expression data are often collected from tissue samples that are composed of multiple cell types. Studies of cell type composition based on gene expression data from tissue samples have recently attracted increasing research interest and led to new method development for cell type composition estimation. This new information on cell type composition can be associated with individual characteristics (e.g., genetic variants) or clinical outcomes (e.g., survival time). Such association analysis can be conducted for each cell type separately followed by multiple testing correction. An alternative approach is to evaluate this association using the composition of all the cell types, thus aggregating association signals across cell types. A key challenge of this approach is to account for the dependence across cell types. We propose a new method to quantify the distances between cell types while accounting for their dependencies, and use this information for association analysis. We demonstrate our method in two applied examples: to assess the association between immune cell type composition in tumor samples of colorectal cancer patients versus survival time and SNP genotypes. We found immune cell composition has prognostic value, and our distance metric leads to more accurate survival time prediction than other distance metrics that ignore cell type dependencies. In addition, survival time-associated SNPs are enriched among the SNPs associated with immune cell composition.

Salicylic Acid and Risk of Colorectal Cancer: A Two-Sample Mendelian Randomization Study

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes' coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84-1.27 and OR: 1.08, 95% CI: 0.86-1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.