Multi-omics machine learning to study host-microbiome interactions in early-onset colorectal cancer.

The incidence of early-onset colorectal cancer (eoCRC) is rising, and its pathogenesis is not completely understood. We hypothesized that machine learning utilizing paired tissue microbiome and plasma metabolome features could uncover distinct host-microbiome associations between eoCRC and average-onset CRC (aoCRC). Individuals with stages I-IV CRC (n = 64) were categorized as eoCRC (age ≤ 50, n = 20) or aoCRC (age ≥ 60, n = 44). Untargeted plasma metabolomics and 16S rRNA amplicon sequencing (microbiome analysis) of tumor tissue were performed. We fit DIABLO (Data Integration Analysis for Biomarker Discovery using Latent variable approaches for Omics studies) to construct a supervised machine-learning classifier using paired multi-omics (microbiome and metabolomics) data and identify associations unique to eoCRC. A differential association network analysis was also performed. Distinct clustering patterns emerged in multi-omic dimension reduction analysis. The metabolomics classifier achieved an AUC of 0.98, compared to AUC 0.61 for microbiome-based classifier. Circular correlation technique highlighted several key associations. Metabolites glycerol and pseudouridine (higher abundance in individuals with aoCRC) had negative correlations with Parasutterella, and Ruminococcaceae (higher abundance in individuals with eoCRC). Cholesterol and xylitol correlated negatively with Erysipelatoclostridium and Eubacterium, and showed a positive correlation with Acidovorax with higher abundance in individuals with eoCRC. Network analysis revealed different clustering patterns and associations for several metabolites e.g.: urea cycle metabolites and microbes such as Akkermansia. We show that multi-omics analysis can be utilized to study host-microbiome correlations in eoCRC and demonstrates promising biomarker potential of a metabolomics classifier. The distinct host-microbiome correlations for urea cycle in eoCRC may offer opportunities for therapeutic interventions.

Differences in Tumor-Associated T cell receptor repertoires between Early-Onset and Average-Onset colorectal cancer.

The incidence of colorectal cancer (CRC) among individuals younger than age 50 (early onset CRC; EOCRC) has substantially increased, yet the etiology and molecular mechanisms underlying this alarming rise remain unclear. We compared tumor-associated T cell repertoires between EOCRC and average-onset CRC (AOCRC) to uncover potentially unique immune microenvironment-related features by age of onset. Our discovery cohort included 242 patients who underwent surgical resection at Cleveland Clinic from 2000 to 2020. EOCRC was defined as age < 50 years at diagnosis (N = 126), and AOCRC as age ≥ 60 years (N = 116). T cell receptor (TCR) abundance and clonality were measured by immunosequencing of tumors. Logistic regression models were used to evaluate the associations between TCR repertoire features and age of onset, adjusting for sex, race, tumor location, and stage. Findings were replicated in 152 EOCRC and 1,984 AOCRC cases from the Molecular Epidemiology of Colorectal Cancer Study. EOCRC tumors had significantly higher TCR diversity compared to AOCRC tumors in the discovery cohort (Odds Ratio (OR):0.44, 95% Confidence Interval (CI):0.32-0.61, p < .0001). This association was also observed in the replication cohort (OR : 0.74, 95% CI : 0.62-0.89, p = .0013). No significant differences in TCR abundance were observed between EOCRC and AOCRC in either cohort. Higher TCR diversity, suggesting a more diverse intratumoral T cell response, is more frequently observed in EOCRC than AOCRC. Further studies are warranted to investigate the role of T cell diversity and the adaptive immune response more broadly in the etiology and outcomes of EOCRC.

Identification of potential mediators of the relationship between body mass index and colorectal cancer: a Mendelian randomization analysis.

BACKGROUND: Colorectal cancer (CRC) is the third-most-common cancer worldwide and its rates are increasing. Elevated body mass index (BMI) is an established risk factor for CRC, although the molecular mechanisms behind this association remain unclear. Using the Mendelian randomization (MR) framework, we aimed to investigate the mediating effects of putative biomarkers and other CRC risk factors in the association between BMI and CRC. METHODS: We selected as mediators biomarkers of established cancer-related mechanisms and other CRC risk factors for which a plausible association with obesity exists, such as inflammatory biomarkers, glucose homeostasis traits, lipids, adipokines, insulin-like growth factor 1 (IGF1), sex hormones, 25-hydroxy-vitamin D, smoking, physical activity (PA) and alcohol consumption. We used inverse-variance weighted MR in the main univariable analyses and performed sensitivity analyses (weighted-median, MR-Egger, Contamination Mixture). We used multivariable MR for the mediation analyses. RESULTS: Genetically predicted BMI was positively associated with CRC risk [odds ratio per SD (5 kg/m2) = 1.17, 95% CI: 1.08-1.24, P-value = 1.4 × 10-5] and robustly associated with nearly all potential mediators. Genetically predicted IGF1, fasting insulin, low-density lipoprotein cholesterol, smoking, PA and alcohol were associated with CRC risk. Evidence for attenuation was found for IGF1 [explained 7% (95% CI: 2-13%) of the association], smoking (31%, 4-57%) and PA (7%, 2-11%). There was little evidence for pleiotropy, although smoking was bidirectionally associated with BMI and instruments were weak for PA. CONCLUSIONS: The effect of BMI on CRC risk is possibly partly mediated through plasma IGF1, whereas the attenuation of the BMI-CRC association by smoking and PA may reflect confounding and shared underlying mechanisms rather than mediation.

Germline genetic regulation of the colorectal tumor immune microenvironment.

OBJECTIVE: To evaluate the contribution of germline genetics to regulating the briskness and diversity of T cell responses in CRC, we conducted a genome-wide association study to examine the associations between germline genetic variation and quantitative measures of T cell landscapes in 2,876 colorectal tumors from participants in the Molecular Epidemiology of Colorectal Cancer Study (MECC). METHODS: Germline DNA samples were genotyped and imputed using genome-wide arrays. Tumor DNA samples were extracted from paraffin blocks, and T cell receptor clonality and abundance were quantified by immunoSEQ (Adaptive Biotechnologies, Seattle, WA). Tumor infiltrating lymphocytes per high powered field (TILs/hpf) were scored by a gastrointestinal pathologist. Regression models were used to evaluate the associations between each variant and the three T-cell features, adjusting for sex, age, genotyping platform, and global ancestry. Three independent datasets were used for replication. RESULTS: We identified a SNP (rs4918567) near RBM20 associated with clonality at a genome-wide significant threshold of 5 × 10- 8, with a consistent direction of association in both discovery and replication datasets. Expression quantitative trait (eQTL) analyses and in silico functional annotation for these loci provided insights into potential functional roles, including a statistically significant eQTL between the T allele at rs4918567 and higher expression of ADRA2A (P = 0.012) in healthy colon mucosa. CONCLUSIONS: Our study suggests that germline genetic variation is associated with the quantity and diversity of adaptive immune responses in CRC. Further studies are warranted to replicate these findings in additional samples and to investigate functional genomic mechanisms.

Tissue-specific genetic variation suggests distinct molecular pathways between body shape phenotypes and colorectal cancer.

It remains unknown whether adiposity subtypes are differentially associated with colorectal cancer (CRC). To move beyond single-trait anthropometric indicators, we derived four multi-trait body shape phenotypes reflecting adiposity subtypes from principal components analysis on body mass index, height, weight, waist-to-hip ratio, and waist and hip circumference. A generally obese (PC1) and a tall, centrally obese (PC3) body shape were both positively associated with CRC risk in observational analyses in 329,828 UK Biobank participants (3728 cases). In genome-wide association studies in 460,198 UK Biobank participants, we identified 3414 genetic variants across four body shapes and Mendelian randomization analyses confirmed positive associations of PC1 and PC3 with CRC risk (52,775 cases/45,940 controls from GECCO/CORECT/CCFR). Brain tissue-specific genetic instruments, mapped to PC1 through enrichment analysis, were responsible for the relationship between PC1 and CRC, while the relationship between PC3 and CRC was predominantly driven by adipose tissue-specific genetic instruments. This study suggests distinct putative causal pathways between adiposity subtypes and CRC.

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.

Spectrum of somatic mutational features of colorectal tumors in ancestrally diverse populations.

Ancestrally diverse and admixed populations, including the Hispanic/Latino/a/x/e community, are underrepresented in cancer genetic and genomic studies. Leveraging the Latino Colorectal Cancer Consortium, we analyzed whole exome sequencing data on tumor/normal pairs from 718 individuals with colorectal cancer (128 Latino, 469 non-Latino) to map somatic mutational features by ethnicity and genetic ancestry. Global proportions of African, East Asian, European, and Native American ancestries were estimated using ADMIXTURE. Associations between global genetic ancestry and somatic mutational features across genes were examined using logistic regression. TP53 , APC , and KRAS were the most recurrently mutated genes. Compared to non-Latino individuals, tumors from Latino individuals had fewer KRAS (OR=0.64, 95%CI=0.41-0.97, p=0.037) and PIK3CA mutations (OR=0.55, 95%CI=0.31-0.98, p=0.043). Genetic ancestry was associated with presence of somatic mutations in 39 genes (FDR-adjusted LRT p<0.05). Among these genes, a 10% increase in African ancestry was associated with significantly higher odds of mutation in KNCN (OR=1.34, 95%CI=1.09-1.66, p=5.74×10 -3 ) and TMEM184B (OR=1.53, 95%CI=1.10-2.12, p=0.011). Among RMGs, we found evidence of association between genetic ancestry and mutation status in CDC27 (LRT p=0.0084) and between SMAD2 mutation status and AFR ancestry (OR=1.14, 95%CI=1.00-1.30, p=0.046). Ancestry was not associated with tumor mutational burden. Individuals with above-average Native American ancestry had a lower frequency of microsatellite instable (MSI-H) vs microsatellite stable tumors (OR=0.45, 95%CI=0.21-0.99, p=0.048). Our findings provide new knowledge about the relationship between ancestral haplotypes and somatic mutational profiles that may be useful in developing precision medicine approaches and provide additional insight into genomic contributions to cancer disparities. Significance: Our data in ancestrally diverse populations adds essential information to characterize mutational features in the colorectal cancer genome. These results will help enhance equity in the development of precision medicine strategies.

Distinct intratumoral microbiome of young-onset and average-onset colorectal cancer.

BACKGROUND: The unexplained rise of young-onset CRC (yoCRC, age <50 years) is of concern. Evidence suggests that microbial dysbiosis may be a contributing factor, but the tumor microbial profile of yoCRC in comparison to average-onset CRC (aoCRC, age >60) has not been fully investigated. METHODS: 16S rRNA amplicon sequencing was performed in tumor and paired adjacent non-malignant fresh frozen tissue specimens prospectively collected from 136 yoCRC and 140 aoCRC patients. Phyloseq, microbiomeSeq, metagenomeSeq, and NetComi were utilized for bioinformatics analysis. Statistical tests included Fisher's exact test, ANOVA, PERMANOVA with Bonferroni correction, linear regression, and Wilcoxon test. p-value <0.05 was considered statistically significant. FINDINGS: yoCRC patients were more likely to have left-sided (72.8 vs. 54.3%), rectal (36.7% vs. 25%), and stage IV (28% vs. 15%) tumors. yoCRC tumors had significantly higher microbial alpha diversity (p = 1.5 × 10-5) and varied beta diversity (R2 = 0.31, p = 0.013) than aoCRC tumors. yoCRC tumors were enriched with Akkermansia and Bacteroides, whereas aoCRC tumors showed greater relative abundances of Bacillus, Staphylococcus, Listeria, Enterococcus, Pseudomonas, Fusobacterium, and Escherichia/Shigella. Akkermansia had a predominantly negative correlation with the microbial communities in yoCRC tumors. yoCRC and aoCRC tumors had distinct microbial profiles associated with tumor location, sidedness, stage, and obesity. Fusobacterium (R2 = -0.23, p = 0.001) and Akkermansia (R2 = 0.05, p = 0.001) abundance correlated with overall survival in yoCRC. INTERPRETATION: Our study provides a comprehensive understanding of the microbial perturbations in yoCRC tumors. We identify microbial candidates that may highlight a distinct pathogenesis of yoCRC and serve as preventive, diagnostic, and therapeutic targets. FUNDING: Sondra and Stephen Hardis Chair in Oncology Research (A.A.K.).

Plasma metabolomic differences in early-onset compared to average-onset colorectal cancer.

Deleterious effects of environmental exposures may contribute to the rising incidence of early-onset colorectal cancer (eoCRC). We assessed the metabolomic differences between patients with eoCRC, average-onset CRC (aoCRC), and non-CRC controls, to understand pathogenic mechanisms. Patients with stage I-IV CRC and non-CRC controls were categorized based on age ≤ 50 years (eoCRC or young non-CRC controls) or  ≥ 60 years (aoCRC or older non-CRC controls). Differential metabolite abundance and metabolic pathway analyses were performed on plasma samples. Multivariate Cox proportional hazards modeling was used for survival analyses. All P values were adjusted for multiple testing (false discovery rate, FDR P < 0.15 considered significant). The study population comprised 170 patients with CRC (66 eoCRC and 104 aoCRC) and 49 non-CRC controls (34 young and 15 older). Citrate was differentially abundant in aoCRC vs. eoCRC in adjusted analysis (Odds Ratio = 21.8, FDR P = 0.04). Metabolic pathways altered in patients with aoCRC versus eoCRC included arginine biosynthesis, FDR P = 0.02; glyoxylate and dicarboxylate metabolism, FDR P = 0.005; citrate cycle, FDR P = 0.04; alanine, aspartate, and glutamate metabolism, FDR P = 0.01; glycine, serine, and threonine metabolism, FDR P = 0.14; and amino-acid t-RNA biosynthesis, FDR P = 0.01. 4-hydroxyhippuric acid was significantly associated with overall survival in all patients with CRC (Hazards ratio, HR = 0.4, 95% CI 0.3-0.7, FDR P = 0.05). We identified several unique metabolic alterations, particularly the significant differential abundance of citrate in aoCRC versus eoCRC. Arginine biosynthesis was the most enriched by the differentially altered metabolites. The findings hold promise in developing strategies for early detection and novel therapies.

Heterozygote advantage at HLA class I and II loci and reduced risk of colorectal cancer.

Objective: Reduced diversity at Human Leukocyte Antigen (HLA) loci may adversely affect the host's ability to recognize tumor neoantigens and subsequently increase disease burden. We hypothesized that increased heterozygosity at HLA loci is associated with a reduced risk of developing colorectal cancer (CRC). Methods: We imputed HLA class I and II four-digit alleles using genotype data from a population-based study of 5,406 cases and 4,635 controls from the Molecular Epidemiology of Colorectal Cancer Study (MECC). Heterozygosity at each HLA locus and the number of heterozygous genotypes at HLA class -I (A, B, and C) and HLA class -II loci (DQB1, DRB1, and DPB1) were quantified. Logistic regression analysis was used to estimate the risk of CRC associated with HLA heterozygosity. Individuals with homozygous genotypes for all loci served as the reference category, and the analyses were adjusted for sex, age, genotyping platform, and ancestry. Further, we investigated associations between HLA diversity and tumor-associated T cell repertoire features, as measured by tumor infiltrating lymphocytes (TILs; N=2,839) and immunosequencing (N=2,357). Results: Individuals with all heterozygous genotypes at all three class I genes had a reduced odds of CRC (OR: 0.74; 95% CI: 0.56-0.97, p= 0.031). A similar association was observed for class II loci, with an OR of 0.75 (95% CI: 0.60-0.95, p= 0.016). For class-I and class-II combined, individuals with all heterozygous genotypes had significantly lower odds of developing CRC (OR: 0.66, 95% CI: 0.49-0.87, p= 0.004) than those with 0 or one heterozygous genotype. HLA class I and/or II diversity was associated with higher T cell receptor (TCR) abundance and lower TCR clonality, but results were not statistically significant. Conclusion: Our findings support a heterozygote advantage for the HLA class-I and -II loci, indicating an important role for HLA genetic variability in the etiology of CRC.

Referral, Uptake, and Outcome of Genetic Counseling and Testing in Patients With Early-Onset Colorectal Cancer.

BACKGROUND: The incidence of early-onset colorectal cancer (EOCRC) is rapidly increasing. Pathogenic germline variants (PGVs) are detected in 16% to 20% of patients who have EOCRC, highlighting a need for genetic counseling (GC) and multigene panel testing in these patients. We aimed to determine the rate of referral to GC and uptake and outcomes of germline testing in patients with EOCRC. METHODS: We conducted a retrospective cohort study of patients aged <50 years diagnosed with colorectal cancer (CRC) from 2010 to 2019 at Cleveland Clinic. Demographic data were extracted, including age, sex, self-reported race, and family history of CRC. The proportions of patients with GC referral and completion of GC and genetic testing were investigated, and genetic testing results were analyzed. Multivariable logistic regression analysis was conducted to identify factors independently associated with GC referral and uptake. RESULTS: A total of 791 patients with EOCRC (57% male and 43% female) were included; 62% were referred for GC, and of those who were referred, 79% completed a GC appointment and 77% underwent genetic testing. Of those who underwent testing, 21% had a PGV detected; 82% were in known CRC-associated genes, with those associated with Lynch syndrome and familial adenomatous polyposis the most common, and 11% were in other actionable genes. Referral to GC was positively associated with family history of CRC (odds ratio [OR], 2.11; 95% CI, 1.51-2.96) and more recent year of diagnosis (2010-2013 vs 2017-2019; OR, 5.36; 95% CI, 3.59-8.01) but negatively associated with older age at diagnosis (OR, 0.89; 95% CI, 0.86-0.92). CONCLUSIONS: Referral to GC for patients with EOCRC is increasing over time; however, even in recent years, almost 25% of patients were not referred for GC. We found that 1 in 5 patients with EOCRC carry actionable PGVs, highlighting the need for health systems to implement care pathways to optimize GC referral and testing in all patients with EOCRC.

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.

Genetically determined circulating resistin concentrations and risk of colorectal cancer: a two-sample Mendelian randomization study.

PURPOSE: Resistin, a novel pro-inflammatory protein implicated in inflammatory processes, has been suggested to play a role in colorectal development. However, evidence from observational studies has been inconsistent. Mendelian randomization may be a complementary method to examine this association. METHODS: We conducted a two-sample Mendelian randomization to estimate the association between genetically determined circulating resistin concentrations and risk of colorectal cancer (CRC). Protein quantitative trait loci (pQTLs) from the SCALLOP consortium were used as instrumental variables (IVs) for resistin. CRC genetic summary data was obtained from GECCO/CORECT/CCFR (the Genetics and Epidemiology of Colorectal Cancer Consortium, Colorectal Cancer Transdisciplinary Study, and Colon Cancer Family Registry), and FinnGen (Finland Biobank). The inverse variance weighted method (IVW) was applied in the main analysis, and other robust methods were used as sensitivity analyses. Estimates for the association from the two data sources were then pooled using a meta-analysis approach. RESULTS: Thirteen pQTLs were identified as IVs explaining together 7.80% of interindividual variation in circulating resistin concentrations. Based on MR analyses, genetically determined circulating resistin concentrations were not associated with incident CRC (pooled-IVW-OR per standard deviation of resistin, 1.01; 95% CI 0.96, 1.06; p = 0.67. Restricting the analyses to using IVs within or proximal to the resistin-encoding gene (cis-IVs), or to IVs located elsewhere in the genome (trans-IVs) provided similar results. The association was not altered when stratified by sex or CRC subsites. CONCLUSIONS: We found no evidence of a relationship between genetically determined circulating resistin concentrations and risk of CRC.

Efforts to Grow Genomic Research in Ancestrally Diverse and Admixed Populations.

Recent initiatives by the research community to characterize the genomic and molecular landscapes of tumors in ancestrally diverse and admixed populations, including the publication by Ding and colleagues in this issue of Cancer Research, represent important efforts to improve our understanding of the entire spectrum of cancer genomic variation with potential clinical consequences. Ding and colleagues confirmed a similar prevalence of mutations in established breast cancer driver genes including PIK3CA, TP53, GATA3, MAP3K1, CDH1, CBFB, PTEN, and RUNX1 and recurrent amplifications in breast cancer drivers including MYC, FGFR1, CCND1, and ERBB2 in tumors from Hispanic/Latina women as compared with non-Hispanic White women. Importantly, they also identified Catalogue of Somatic Mutations in Cancer (COSMIC) signature 16 in a significant fraction of tumors from Hispanic/Latina women and a novel recurrent amplification on 17q11.2. This study highlights the potential for inclusion of participants from diverse populations to accelerate discoveries and advance equity in genomic medicine, as well as the need for even larger collaborative initiatives. See related article by Ding et al., p. 2600.

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.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals.

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Tumor immune infiltration estimated from gene expression profiles predicts colorectal cancer relapse.

A substantial fraction of patients with stage I-III colorectal adenocarcinoma (CRC) experience disease relapse after surgery with curative intent. However, biomarkers for predicting the likelihood of CRC relapse have not been fully explored. Therefore, we assessed the association between tumor infiltration by a broad array of innate and adaptive immune cell types and CRC relapse risk. We implemented a discovery-validation design including a discovery dataset from Moffitt Cancer Center (MCC; Tampa, FL) and three independent validation datasets: (1) GSE41258 (2) the Molecular Epidemiology of Colorectal Cancer (MECC) study, and (3) GSE39582. Infiltration by 22 immune cell types was inferred from tumor gene expression data, and the association between immune infiltration by each cell type and relapse-free survival was assessed using Cox proportional hazards regression. Within each of the four independent cohorts, CD4+ memory activated T cell (HR: 0.93, 95% CI: 0.90-0.96; FDR = 0.0001) infiltration was associated with longer time to disease relapse, independent of stage, microsatellite instability, and adjuvant therapy. Based on our meta-analysis across the four datasets, 10 innate and adaptive immune cell types associated with disease relapse of which 2 were internally validated using multiplex immunofluorescence. Moreover, immune cell type infiltration was a better predictors of disease relapse than Consensus Molecular Subtype (CMS) and other expression-based biomarkers (Immune-AICMCC:238.1-238.9; CMS-AICMCC: 241.0). These data suggest that transcriptome-derived immune profiles are prognostic indicators of CRC relapse and quantification of both innate and adaptive immune cell types may serve as candidate biomarkers for predicting prognosis and guiding frequency and modality of disease surveillance.

Rare Variants in the DNA Repair Pathway and the Risk of Colorectal Cancer.

BACKGROUND: Inherited susceptibility is an important contributor to colorectal cancer risk, and rare variants in key genes or pathways could account in part for the missing proportion of colorectal cancer heritability. METHODS: We conducted an exome-wide association study including 2,327 cases and 2,966 controls of European ancestry from three large epidemiologic studies. Single variant associations were tested using logistic regression models, adjusting for appropriate study-specific covariates. In addition, we examined the aggregate effects of rare coding variation at the gene and pathway levels using Bayesian model uncertainty techniques. RESULTS: In an exome-wide gene-level analysis, we identified ST6GALNAC2 as the top associated gene based on the Bayesian risk index (BRI) method [summary Bayes factor (BF)BRI = 2604.23]. A rare coding variant in this gene, rs139401613, was the top associated variant (P = 1.01 × 10-6) in an exome-wide single variant analysis. Pathway-level association analyses based on the integrative BRI (iBRI) method found extreme evidence of association with the DNA repair pathway (BFiBRI = 17852.4), specifically with the nonhomologous end joining (BFiBRI = 437.95) and nucleotide excision repair (BFiBRI = 36.96) subpathways. The iBRI method also identified RPA2, PRKDC, ERCC5, and ERCC8 as the top associated DNA repair genes (summary BFiBRI ≥ 10), with rs28988897, rs8178232, rs141369732, and rs201642761 being the most likely associated variants in these genes, respectively. CONCLUSIONS: We identified novel variants and genes associated with colorectal cancer risk and provided additional evidence for a role of DNA repair in colorectal cancer tumorigenesis. IMPACT: This study provides new insights into the genetic predisposition to colorectal cancer, which has potential for translation into improved risk prediction.

Genetic architectures of proximal and distal colorectal cancer are partly distinct.

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.