Harrison, Tabitha A.

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

Abstract Background: Colorectal cancer (CRC) incidence has increased in those 50 years or younger compared to older individuals in recent decades. While the intra-tumoral microbiome is a potent facilitator of tumorigenesis, the role of the microbiome in early onset (eo)CRC is unclear. Identification of the intra-tumoral microbiome may enhance prediction of adverse outcomes to improve survival of CRC overall age groups and eoCRC. Methods: We performed targeted sequencing of tumor and matched normal DNA samples from 3695 colorectal cancer cases within the Genetics and Colorectal Cancer Consortium (GECCO) and the Colon Cancer Family Registry (CCFR). Tumor DNA was extracted from formalin fixed paraffin embedded (FFPE) or fresh frozen tissue. Normal tissue DNA was extracted from blood or normal surrounding colorectal tissue. We designed probes for the bacterial pathogens enterotoxigenic Bacteroides fragilis (ETBF), polyketide synthase positive E. coli (pks+ EC), and Fusobacterium nucleatum (Fn). We used logistic regression to assess tumor-associated bacteria and clinical parameters, tumor molecular characteristics, and tumor stage adjusted for age, sex, and study and a Cox promotional hazard to assess association with CRC-specific survival, overall and in eoCRC, adjusted for age at diagnosis, sex, mutational burden, and microsatellite instability (MSI) and stratified baseline hazards by study. Survival data were available for 2357 and 308 patients, overall and in eoCRC, respectively. Results: The prevalence of different bacteria in tumors differed slightly overall vs. in eoCRC tumors: ETBF (5.2%, 6.3%), pks+ EC (12.9%, 14.8%), and Fn (10.4%, 8.2%). Overall, men had a lower probability of tumors with Fn (odds ratio (OR)=0.8, 95% confidence interval (CI)=0.64-0.99, p=0.04) but higher odds of tumors with pks+ EC (OR=1.38, CI=1.14-1.68, p=0.001), compared to women. The odds of prevalent Fn were higher for proximal vs distal tumors (OR=1.76, CI=1.38-2.26, p<0.001), but there was no association between tumor site and pks+ EC or ETBF. MSI-high tumors were more likely to carry Fn (OR=2.79, CI=2.08-3.75, p<0.001) and less likely to carry pks+ EC (OR=0.63, CI=0.44-0.88, p=0.008) when compared to patients with MSS tumors. In patients with non-hypermutated tumors, we observed that pks+ EC positive tumors were associated with better CRC-specific survival (HR=0.73, CI=0.57-0.93, p=0.01), whereas ETBF positive tumors (HR=1.99, CI=1,18-3.36,p=0.009) and Fn positive tumors (HR=1.33, CI=1.04-1.72, p=0.03) were associated with poorer survival. In eo-CRC patients, Fn positive status was associated with poorer survival (HR=2.4, CI= 1.056 - 5.492, p=0.026). Conclusions: We showed that the intra-tumoral microbiome differs in distribution across patient populations and patient survival over all age groups and, specifically, in eoCRC. Citation Format: Meredith A. Hullar, Keith R. Curtis, Tabitha Harrison, Yi Lin, Robert Steinfelder, Sonja I. Berndt, Daniel D. Buchanan, Andrew T. Chan, David A. Drew, Jane Figueiredo, Amy J. French, Tabitha A. Harrison, Mark A. Jenkins, Johannes Melaku, Victor Moreno, Tomotaka Ugai, Shuji Ogino, Conghui Qu, Chenxu Qu, Steven Thibodeau, Syed Zaidi, Amanda I. Phipps, Ulrike Peters, on behalf of GECCO and CCFR investigators. Evaluation of intra-tumoral pks+ E. coli, enterotoxigenic B. fragilis and Fusobacterium nucleatum, overall and in early onset disease, in colorectal cancer cases [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 3039.

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 Mendelian randomization (MR) studies are susceptible to metadata errors (e.g. incorrect specification of the effect allele column) and other analytical issues that can introduce substantial bias into analyses. We developed a quality control (QC) pipeline for the Fatty Acids in Cancer Mendelian Randomization Collaboration (FAMRC) that can be used to identify and correct for such errors. Methods We collated summary association statistics from fatty acid and cancer genome-wide association studies (GWAS) and subjected the collated data to a comprehensive QC pipeline. We identified metadata errors through comparison of study-specific statistics to external reference data sets (the National Human Genome Research Institute-European Bioinformatics Institute GWAS catalogue and 1000 genome super populations) and other analytical issues through comparison of reported to expected genetic effect sizes. Comparisons were based on three sets of genetic variants: (i) GWAS hits for fatty acids, (ii) GWAS hits for cancer and (iii) a 1000 genomes reference set. Results We collated summary data from 6 fatty acid and 54 cancer GWAS. Metadata errors and analytical issues with the potential to introduce substantial bias were identified in seven studies (11.6%). After resolving metadata errors and analytical issues, we created a data set of 219 842 genetic associations with 90 cancer types, generated in analyses of 566 665 cancer cases and 1 622 374 controls. Conclusions In this large MR collaboration, 11.6% of included studies were affected by a substantial metadata error or analytical issue. By increasing the integrity of collated summary data prior to their analysis, our protocol can be used to increase the reliability of downstream MR analyses. Our pipeline is available to other researchers via the CheckSumStats package (https://github.com/MRCIEU/CheckSumStats).

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 P trend <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%; P trend <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.

Abstract 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 &lt;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 &lt;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 &lt;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.

Abstract Background: The microbiome in colorectal cancer (CRC) tumors is now recognized as a potent facilitator of tumorigenesis. Identification of the intra-tumoral microbiome may provide an avenue to enhance prediction of adverse outcomes and ultimately lead to improved survival of CRC, the second leading cause of cancer death. Methods: We performed targeted sequencing of tumor and matched normal DNA samples for 3695 colorectal cancer cases within the Genetics and Colorectal Cancer Consortium (GECCO) and the Colon Cancer Family Registry (CCFR). Tumor DNA was extracted primarily from formalin fixed paraffin embedded (FFPE) or from fresh frozen tissue. FFPE slides were dissected to ensure &gt;70% tumor content. Normal DNA was primarily extracted from blood and otherwise normal surrounding colorectal tissue. We designed probes for bacterial pathogens including Bacteroides fragilis (ETBF), polyketide synthase positive E. coli (pks+ EC), enterotoxigenic and Fusobacterium nucleatum (Fn). Logistic regression was used to assess tumor-associated bacteria and various clinical parameters, tumor molecular characteristics, and tumor stage adjusted for age, sex, and study and a Cox promotional hazard to assess association with CRC-specific survival adjusted for age at diagnosis, sex, mutational burden and microsatellite instability (MSI) and stratified baseline hazards by study. Survival data were available for 2357 patients. Results: The prevalence of patients with bacteria in tumors varied for ETBF (5.2%), pks+ EC (12.9%), and Fn (10.4%). Men had a lower probability of tumors with Fn (odds ratio (OR)=0.8, 95% confidence interval (CI)=0.64-0.99,p=0.04) but higher odds of tumors with pks+ EC (OR=1.38, CI=1.14-1.68,p=0.001). Patients with proximal tumors had higher odds of carrying Fn (OR=1.76, CI=1.38-2.26, p&lt;0.001) as did patients with MSI high tumors (OR 2.79, CI=2.08-3.74) p&lt;0.001). Patients with MSI-high tumors had lower odds of carrying pks+ EC (OR=0.63, CI=0.44-0.88, p=0.008) when compared to patients with microsatellite stable tumors. Among non-hypermutated tumors, we observed that tumors positive for pks+ EC were associated with better CRC-specific survival (OR=0.73, CI=0.57-0.93,p=0.01 whereas in patients with ETBF positive tumors had shorter survival (OR=1.99,CI=1,18-3.36,p=0.009) as did patients with tumors positive for Fn (OR=1.33, CI=1.04-1.72, p=0.03). Conclusions: In this largest investigation of intra-tumoral bacteria in CRC, we demonstrate that inclusion of pathogenic intra-tumoral microbiome may enhance prediction of survival and can identify a subset of patients with CRC that may benefit from targeted therapies that modulate the tumor microbiome. Citation Format: Meredith A. J. Hullar, Keith R. Curtis, Tabitha A. Harrison, Yi Lin, Robert S. Steinfelder, Sonja I. Berndt, Daniel D. Buchanan, Andrew T. Chan, David Drew, Jane C. Figueiredo, Amy J. French, Ellen L. Goode, Mark A. Jenkins, Yohannes A. Melaku, Victor Moreno, Shuji Ogino, Stephen N. Thibideau, Amanda I. Phipps, Ulrike Peters. Evaluation of intra-tumoral pathogenic bacteria pks+ E. coli, enterotoxigenic B. fragilis and fusobacterium nucleatum in 3695 colorectal cancer 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 PR008.

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

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: 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 &lt; 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.