Buchanan, Daniel D.

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 Recent research suggests higher circulating lymphocyte counts may protect against colorectal cancer (CRC). However, the role of specific lymphocyte subtypes and activation states remain unclear. CD4+ T cells—a highly dynamic lymphocyte subtype—undergo gene expression changes upon activation that are critical to their effector function. Previous studies using bulk tissue have limited our understanding of their role in CRC risk to static associations. We applied Mendelian randomization (MR) and genetic colocalisation to evaluate causal relationships of gene expression on CRC risk across multiple CD4+ T cell subtypes and activation states. Genetic proxies were obtained from single-cell transcriptomic data, allowing us to investigate the causal effect of expression of 1,805 genes across CD4+ T cell activation states on CRC risk (78,473 cases; 107,143 controls). Analyses were stratified by CRC anatomical subsites and sex, with sensitivity analyses assessing whether the observed effect estimates were likely to be CD4+ T cell-specific. We identified 6 genes—FADS2, FHL3, HLA-DRB1, HLA-DRB5, RPL28, and TMEM258—with strong evidence for a causal role in CRC development (FDR-P &lt; 0.05; colocalisation H4 &gt; 0.8). Causal estimates varied by CD4+ T cell subtype, activation state, CRC subsite and sex. However, many of genetic proxies used to instrument gene expression in CD4+ T cells also act as eQTLs in other tissues, highlighting the challenges of using genetic proxies to instrument tissue-specific expression changes. We demonstrate the importance of capturing the dynamic nature of CD4+ T cells in understanding CRC risk, and prioritize genes for further investigation in cancer prevention.

Abstract Introduction There is convincing evidence that overall adiposity increases the risks of several cancers. Whether the distribution of adiposity plays a similar role is unclear. Methods We used 2-sample Mendelian randomization (MR) to examine causal relationships of 5 adiposity distribution traits (abdominal subcutaneous adipose tissue (ASAT); visceral adipose tissue (VAT); gluteofemoral adipose tissue (GFAT); liver fat; and pancreas fat) with the risks of 12 obesity-related cancers (endometrial, ovarian, breast, colorectal, pancreas, multiple myeloma, liver, kidney (renal cell), thyroid, gallbladder, esophageal adenocarcinoma, and meningioma). Results Sample size across all genome-wide association studies (GWAS) ranged from 8407 to 728 896 (median: 57 249). We found evidence that higher genetically predicted ASAT increased the risks of endometrial cancer, liver cancer, and esophageal adenocarcinoma (odds ratios (OR) and 95% confidence intervals (CI) per standard deviation (SD) higher ASAT = 1.79 (1.18 to 2.71), 3.83 (1.39 to 10.53), and 2.34 (1.15 to 4.78), respectively). Conversely, we found evidence that higher genetically predicted GFAT decreased the risks of breast cancer and meningioma (ORs and 95% CIs per SD higher genetically predicted GFAT = 0.77 (0.62 to 0.97) and 0.53 (0.32 to 0.90), respectively). We also found evidence for an effect of higher genetically predicted VAT and liver fat on increased liver cancer risk (ORs and 95% CIs per SD higher genetically predicted adiposity trait = 4.29 (1.41 to 13.07) and 4.09 (2.29 to 7.28), respectively). Discussion Our analyses provide novel insights into the relationship between adiposity distribution and cancer risk. These insights highlight the potential importance of adipose tissue distribution alongside maintaining a healthy weight for cancer prevention.

Abstract Background Physical activity (PA) is an established protective factor for colorectal cancer (CRC), but it is unclear if genetic variants modify this effect. To investigate this possibility, we conducted a genome-wide gene–PA interaction analysis. Methods Using logistic regression and two-step and joint tests, we analyzed interactions between common genetic variants across the genome and PA in relation to CRC risk. Self-reported PA levels were categorized as active (≥ 8.75 MET-h/wk) vs. inactive (< 8.75 MET-h/wk) and as study- and sex-specific quartiles of activity. Results PA had an overall protective effect on CRC (OR [active vs. inactive] = 0.85; 95%CI = 0.81–0.90). The two-step GxE method identified an interaction between rs4779584, an intergenic variant near the GREM1 and SCG5 genes, and PA for CRC risk (p-interaction = 2.6×10 − 8 ). Stratification by genotype at this locus showed a significant reduction in CRC risk by 20% in active vs. inactive participants with the CC genotype (OR = 0.80; 95%CI = 0.75–0.85), but no significant PA–CRC association among CT or TT carriers. When PA was modeled as quartiles, the 1-d.f. GxE test identified that rs56906466, an intergenic variant near the KCNG1 gene, modified the association between PA and CRC (p-interaction = 3.5×10 − 8 ). Stratification at this locus showed that increase in PA (highest vs. lowest quartile) was associated with a lower CRC risk solely among TT carriers (OR = 0.77; 95%CI = 0.72–0.82). Conclusions In summary, we identified two genetic variants that modified the association between PA and CRC risk. One of them, related to GREM1 and SCG5 , suggests that the bone morphogenetic protein (BMP)-related, inflammatory, and/or insulin signaling pathways may be associated with the protective influence of PA on colorectal carcinogenesis.

Abstract Background Being able to estimate the risk of metachronous disease in a patient with colorectal cancer (CRC) could enable risk-appropriate surveillance. The aim of this study was to develop a risk-prediction model to estimate individual 10-year risk of metachronous disease following a CRC diagnosis. Methods A population-based cohort of patients with CRC was recruited soon after diagnosis between 1997 and 2012 from the United States, Canada, and Australia. Cox regression with the least absolute shrinkage and selection operator penalization was used to identify factors that predicted the risk of a new primary CRC diagnosed at least 1 year after the initial CRC diagnosis. Potential predictors included demography, anthropometry, lifestyle factors, comorbidities, personal and family cancer history, medication use, age at diagnosis, and pathological features of the first CRC. Internal validation through bootstrapping was used to evaluate the discrimination and calibration. Results We included 6085 CRC cases; 138 (2.3%) of these cases were diagnosed with metachronous disease over a median of 12 years (IQR = 5-17 years). Metachronous CRC risk was predicted by body mass index; smoking status; level of physical activity; family history of cancer and synchronous CRC; stage, grade, histological type, and DNA mismatch repair status; and age at diagnosis of the first CRC. The model was valid with a C statistic of 0.65 (95% CI = 0.63 to 0.68) and a calibration slope of 0.873 (SD = 0.087). Conclusions Metachronous CRC can be predicted with reasonable accuracy using a prediction model that consists of clinical variables collected as part of routine practice.

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 Menopausal users of hormone replacement therapy (HRT) are at increased breast cancer risk and decreased colorectal cancer (CRC) risk compared with individuals who have never used HRT, but these opposing associations may differ by familial risk of breast cancer and CRC. We harmonized data from 3 cohorts and generated separate breast cancer and CRC familial risk scores based on cancer family history. We defined moderate or strong family history as a risk score of 0.4 or higher, where 0.4 was equivalent to a 50-year-old woman with 1 parent diagnosed with either breast cancer or CRC at 55 years of age. Of 24 486 women assessed, 1243 and 405 were diagnosed with incident breast cancer and CRC, respectively. For breast cancer, menopausal HRT ever use versus never use hazard ratios were 1.27 (95% CI = 1.11 to 1.45) for a breast cancer familial risk score below 0.4 and 1.01 (95% CI = 0.82 to 1.25) for a breast cancer familial risk score of 0.4 or higher (Pdifference = .08). For CRC, menopausal HRT hazard ratios were 0.63 (95% CI = 0.50 to 0.78) for a CRC familial risk score below 0.4 and 1.21 (95% CI = 0.73 to 2.00) for a CRC familial risk score of 0.4 or higher (Pdifference = .03). Associations with menopausal HRT use that apply to the general population may not hold for women at moderate or strong familial risk of these cancers.

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

Abstract Deep learning may detect biologically important signals embedded in tumor morphologic features that confer distinct prognoses. Tumor morphologic features were quantified to enhance patient risk stratification within DNA mismatch repair (MMR) groups using deep learning. Using a quantitative segmentation algorithm (QuantCRC) that identifies 15 distinct morphologic features, we analyzed 402 resected stage III colon carcinomas [191 deficient (d)-MMR; 189 proficient (p)-MMR] from participants in a phase III trial of FOLFOX-based adjuvant chemotherapy. Results were validated in an independent cohort (176 d-MMR; 1,094 p-MMR). Association of morphologic features with clinicopathologic variables, MMR, KRAS, BRAFV600E, and time-to-recurrence (TTR) was determined. Multivariable Cox proportional hazards models were developed to predict TTR. Tumor morphologic features differed significantly by MMR status. Cancers with p-MMR had more immature desmoplastic stroma. Tumors with d-MMR had increased inflammatory stroma, epithelial tumor-infiltrating lymphocytes (TIL), high-grade histology, mucin, and signet ring cells. Stromal subtype did not differ by BRAFV600E or KRAS status. In p-MMR tumors, multivariable analysis identified tumor-stroma ratio (TSR) as the strongest feature associated with TTR [HRadj 2.02; 95% confidence interval (CI), 1.14–3.57; P = 0.018; 3-year recurrence: 40.2% vs. 20.4%; Q1 vs. Q2–4]. Among d-MMR tumors, extent of inflammatory stroma (continuous HRadj 0.98; 95% CI, 0.96–0.99; P = 0.028; 3-year recurrence: 13.3% vs. 33.4%, Q4 vs. Q1) and N stage were the most robust prognostically. Association of TSR with TTR was independently validated. In conclusion, QuantCRC can quantify morphologic differences within MMR groups in routine tumor sections to determine their relative contributions to patient prognosis, and may elucidate relevant pathophysiologic mechanisms driving prognosis. Significance: A deep learning algorithm can quantify tumor morphologic features that may reflect underlying mechanisms driving prognosis within MMR groups. TSR was the most robust morphologic feature associated with TTR in p-MMR colon cancers. Extent of inflammatory stroma and N stage were the strongest prognostic features in d-MMR tumors. TIL density was not independently prognostic in either MMR group.