Huang, Wen‐Yi

Background: Considering the pathophysiological importance of the intestine-passing microbiota, we hypothesized that the mutational and microbial profiles of colorectal cancer might gradually change along detailed anatomical locations and that these trends might differ by the age of diagnosis and the microsatellite instability (MSI) status. Methods: We assessed somatic mutations and selected bacterial species in colorectal carcinomas along the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum, utilizing targeted sequencing data ( N = 5685) from the Genetics and Epidemiology of Colorectal Cancer Consortium and tissue microbial data ( N = 1389) from the Health Professionals Follow-up Study and Nurses’ Health Study. Results: There were decreasing trends from the cecum to rectum in the prevalence of gene mutations in the IGF/PI3K, receptor tyrosine kinase (RTK)/RAS, TGFB, and WNT signaling pathways (all P trend < 0.0001), as well as high levels of Bacteroides fragilis ( P trend < 0.0001) and high levels of Fusobacterium nucleatum ( P trend = 0.0029) (but not enterotoxigenic Bacteroides fragilis or pks+ Escherichia coli ). Conversely, TP53 pathway mutations showed an increasing cecum-to-rectum trend ( P trend < 0.0001). Similar to later-onset cases diagnosed at ages 50 years and older, in early-onset cases diagnosed under the age of 50 years, we found generally decreasing cecum-to-rectum trends of IGF/PI3K, RTK/RAS, TGFB, and WNT pathway mutations and an increasing cecum-to-rectum trend in TP53 pathway mutations. By contrast, the cecum-to-rectum prevalence-changing trends in IGF/PI3K, RTK/RAS, TGFB, and TP53 pathway mutations were present in non–MSI-high tumors but not in MSI-high tumors (all P interaction < 0.0052). Conclusion: Non–MSI-high colorectal cancers exhibited decreasing cecum-to-rectum trends of IGF/PI3K, RTK/RAS, and TGFB pathway mutations (and an increasing trend of TP53 pathway mutations), in contrast to MSI-high tumors that demonstrated no such trends. Our study highlights the importance of detailed colorectal subsite information combined with tumor molecular and microbial characteristics in colorectal cancer research.

Abstract Background Physical activity is associated with lower colorectal cancer (CRC) risk, but its association with molecular subtypes defined by genetic and epigenetic alterations of the disease is unclear. Such information may enhance the understanding of the mechanisms related to the benefits of physical activity. Methods Pooled observational (cases: n = 5386; controls: n = 6798; studies n = 5) and genome-wide association data (cases: n = 8178; controls: n = 10 472; studies n = 5) were used. We used multivariable logistic regression models and Mendelian randomization to assess the association between physical activity and the risk of CRC subtypes defined by individual tumor markers (and marker combinations), namely microsatellite instability status, CpG island methylator phenotype status, and BRAF and KRAS mutations. We used case-only analysis to test for differences between molecular subtypes. We applied Bonferroni correction to account for multiple tests. Results In the pooled observational analysis, higher levels of physical activity were associated with lower CRC risk (Obs-per 1SD, odds ratio [OR] = 0.94, 95% confidence interval [CI] = 0.90 to 0.97), with an association that was stronger in males (Obs-per 1SD, OR = 0.91, 95% CI = 0.87 to 0.96) than in females (Obs-per 1SD, OR = 0.97, 95% CI = 0.91 to 1.03; Pinteraction = .04). Higher physical activity was associated with a lower risk of CRC across all molecular subtypes, especially in males. There was no difference in the associations by subtypes by pooled observational or Mendelian randomization analyses. The findings did not differ by study design, anatomical site, and early or late age onset of CRC. Conclusions Our findings suggest that physical activity is not differentially associated with the 4 major molecular subtypes involved in colorectal carcinogenesis, indicating that its benefits extend broadly across colorectal cancer pathogenesis.

Abstract Background: Waist circumference (WC) and its allometric counterpart, “a body shape index” (ABSI), are risk factors for colorectal cancer; however, it is uncertain whether associations with these body measurements are limited to specific molecular subtypes of the disease. Methods: Data from 2,772 colorectal cancer cases and 3,521 controls were pooled from four cohort studies within the Genetics and Epidemiology of Colorectal Cancer Consortium. Four molecular markers (BRAF mutation, KRAS mutation, CpG island methylator phenotype, and microsatellite instability) were analyzed individually and in combination (Jass types). Multivariable logistic and multinomial logistic models were used to assess the associations of WC and ABSI with overall colorectal cancer risk and, in case-only analyses, to evaluate heterogeneity by molecular subtype, respectively. Results: Higher WC (ORper 5 cm = 1.06, 95% confidence interval, 1.04–1.09) and ABSI (ORper 1-SD = 1.07, 95% confidence interval, 1.00–1.14) were associated with elevated colorectal cancer risk. There was no evidence of heterogeneity between the molecular subtypes. No difference was observed regarding the influence of WC and ABSI on the four major molecular markers in proximal colon, distal colon, and rectal cancers, as well as in early- and late-onset colorectal cancers. Associations did not differ in the Jass-type analysis. Conclusions: Higher WC and ABSI were associated with elevated colorectal cancer risk; however, they do not differentially influence all four major molecular mutations involved in colorectal carcinogenesis but underscore the importance of maintaining a healthy body weight in colorectal cancer prevention. Impact: The proposed results have potential utility in colorectal cancer prevention.

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

Abstract Background: The genotoxin colibactin, produced by Escherichia coli containing the pks island, is associated with a tumor single base substitution (SBS) mutational signature SBS88. It is unknown whether lifestyle and environmental factors related to gut dysbiosis impact colorectal cancer (CRC) risk and survival differently by colibactin signature status. Methods: Utilizing the resources of the Genetic Epidemiology of Colorectal cancer Consortium and Colon Cancer Family Registry (GECCO-CCFR), a consortium including both case-control and cohort study designs, SBS88 tumor mutational signature was measured in 4308 CRC cases who were microsatellite stable (MSS) and had at least five somatic single nucleotide variants for signature decomposition. Associations of 13 lifestyle and environmental factors with CRC risk by colibactin signature status were assessed using logistic and multinomial regression; associations with CRC survival were assessed using Cox proportional hazards regression. Results: 392 (9%) of the 4308 MSS tumors were positive for SBS88. Having a BMI ≥30 kg/m2 (versus a BMI 18.5-24.9 kg/m2) was strongly associated with worse CRC-specific survival among those with the SBS88 signature, but not among those with SBS88 negative CRC [further adjusted hazard ratio (HR) = 3.40 95% confidence interval (CI) 1.47, 7.84, and HR = 0.97 95% CI 0.78, 1.21, respectively, P for heterogeneity = 0.066]. Associations of fruit intake with CRC risk may be differential by SBS88 status. In minimally adjusted models, those in the highest quartile of fruit intake (versus the first quartile) had 0.53 (95% CI 0.37, 0.76) times the risk of SBS88 positive CRC, as opposed to 0.75 (95% CI 0.66, 0.85) times the risk of SBS88 negative CRC (P for heterogeneity = 0.047). Among cohort studies specifically, associations between CRC risk and BMI, alcohol, and fruit intake also differed by SBS88 status (P for heterogeneity ranging 0.02-0.05). Conclusions: Higher BMI may be associated with worse CRC-specific survival among those with SBS88 positive tumors. BMI, alcohol, and dietary factors may be differentially associated with CRC risk based on SBS88 status. Citation Format: Claire Elizabeth Thomas, Peter Georgeson, Conghui Qu, Robert S. Steinfelder, Daniel D. Buchanan, Mark A. Jenkins, Andrea Gsur, Jane C. Figueiredo, Polly A. Newcomb, Peter T. Campbell, Christina Newton, Caroline Y. Um, Lisa Boardman, Victor Moreno, Jenny Chang-Claude, Michael Hoffmeister, Marc J. Gunter, Andrew T. Chan, Shuji Ogino, Ellen L. Goode, Brigid Lynch, Yohannes A. Melaku, Steven Gallinger, Sonja I. Berndt, Wen-Yi Huang, Ulrike Peters, Amanda I. Phipps. Lifestyle and environmental factors in relation to colorectal cancer risk and survival by colibactin tumor mutational signature status [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6453.

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

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

Abstract Background 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 (GxE) interactions and may provide information on the underlying biological pathway. Objective We investigated interactions between folate consumption and predicted gene expression on colorectal cancer risk across the genome. Methods We used variant weights from the PrediXcan models of colon tissue-specific gene expression as a priori variant information for a set-based GxE 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 4,839 genes with available genetically predicted expression. We meta-analyzed results across 23 studies for a total of 13,498 cases with colorectal tumors and 13,918 controls of European ancestry. Results We found suggestive evidence of interaction with folate intake 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). 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. Conclusion We identified three genes involved in preventing or repairing DNA damage that may interact with folate consumption to alter CRC risk.

Abstract Colorectal cancer (CRC) is a heterogenous disease that develops through somatic mutations in driver genes, leading to activation of diverse neoplastic pathways. To systematically examine if somatically mutated genes and pathways impact survival, we sequenced tumor and normal DNA samples for 4,512 CRC cases using a targeted panel. We performed Cox regression to estimate hazard ratios (HR) and 95% confidence intervals (CIs) for associations of disease-specific (DS)-survival with somatically mutated genes and pathways, adjusting for age and sex, and stratifying baseline hazards by study population. We assessed statistical significance using Bonferroni p-value thresholds to account for multiple testing of 214 genes (2.34x10-4) and 6 key CRC pathways (8.3x10-3). We limited analyses to non-silent mutations. We observed that DS-survival was significantly more favorable among individuals with hypermutated (HM) tumors, primarily consisting of microsatellite unstable and POLE-mutated tumors, as compared to those with non-hypermutated (NHM) tumors (HR=0.4, 95% CI: 0.3-0.5, p=1.2x10-17). BRAF V600E mutations were associated with poorer survival (HR=2.0, 95% CI: 1.6-2.5, p=1.6x10-10). This association was more pronounced among NHM tumors (HR=2.3, 95% CI: 1.8-2.9, p=4.3x10-12). We identified suggestive associations (p-values &lt; 5.0x10-3 in overall or stratified analyses) between DS-survival and mutations in B2M, TP53, and SMAD4. Mutations in B2M may provide more favorable prognosis for survival (HR=0.5, 95% CI: 0.3-0.8, p=4.4x10-3), with similar effect sizes in HM and NHM tumors. Poorer survival may be associated with mutations in TP53 (HR=1.2, 95% CI:1.2-1.4, p=9.3x10-4) and SMAD4 (HR=1.3, 95% CI:1.1-1.6, p=3.2x10-3) in NHM tumors. We further observed statistically significant associations between survival and three mutated pathways: TP53/ATM (HR=1.2, 95% CI:1.1-1.4, p=9.0x10-4), receptor tyrosine kinases (RTK) and RAS (HR=1.3, 95% CI:1.1-1.5, p=4.5x10-5), and TGF-beta (HR=1.3, 95% CI:1.1-1.5, p=2.9x10-4). Findings for TP53/ATM and RTK/RAS were primarily due to mutations in one gene within each pathway (TP53 and BRAF, respectively). However, the TGF-beta pathway finding, which was more pronounced in NHM tumors, included two genes with p-values less than 0.05 (SMAD4 p=3.2x10-3 and TGFBR2 p=0.01). Despite our large sample size, relatively few somatically mutated genes were significantly associated with DS-survival. It may be that pathway-level testing affords more power for analyses and that larger sample sizes are needed. It is of interest that non-silent mutations in B2M, whose product is a component of the class I major histocompatibility complex, was associated with improved DS-survival, though our finding was not statistically significant after multiple testing correction so further exploration and replication is needed. Citation Format: Tabitha A. Harrison, Syed H. Zaidi, Conghui Qu, Amanda I. Phipps, Robert S. Steinfelder, Quang M. Trinh, Sonja I. Berndt, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Kim F. Doheny, David A. Drew, Jane C. Figueiredo, Steven J. Gallinger, Marios Giannakis, Andrea Gsur, Marc J. Gunter, Michael Hoffmeister, Wen-Yi Huang, Paul J. Limburg, Victor Moreno, Polly A. Newcomb, Shuji Ogino, Ross L. Prentice, Tameka Shelford, Wei Sun, Stephen N. Thibodeau, Li Hsu, Ulrike Peters. Associations of somatically mutated genes and pathways with colorectal cancer specific survival in 4,500 colorectal cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB090.

Objective An understanding of the etiologic heterogeneity of colorectal cancer (CRC) is critical for improving precision prevention, including individualized screening recommendations and the discovery of novel drug targets and repurposable drug candidates for chemoprevention. Known differences in molecular characteristics and environmental risk factors among tumors arising in different locations of the colorectum suggest partly distinct mechanisms of carcinogenesis. The extent to which the contribution of inherited genetic risk factors for CRC differs by anatomical subsite of the primary tumor has not been examined. Design To identify new anatomical subsite-specific risk loci, we performed genome-wide association study (GWAS) meta-analyses including data of 48 214 CRC cases and 64 159 controls of European ancestry. We characterised effect heterogeneity at CRC risk loci using multinomial modelling. Results We identified 13 loci that reached genome-wide significance (p<5×10 −8 ) and that were not reported by previous GWASs for overall CRC risk. Multiple lines of evidence support candidate genes at several of these loci. We detected substantial heterogeneity between anatomical subsites. Just over half (61) of 109 known and new risk variants showed no evidence for heterogeneity. In contrast, 22 variants showed association with distal CRC (including rectal cancer), but no evidence for association or an attenuated association with proximal CRC. For two loci, there was strong evidence for effects confined to proximal colon cancer. Conclusion Genetic architectures of proximal and distal CRC are partly distinct. Studies of risk factors and mechanisms of carcinogenesis, and precision prevention strategies should take into consideration the anatomical subsite of the tumour.