Search

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

Abstract Aberrant stem cell-like activity and impaired differentiation are central to the development of colorectal cancer (CRC). To identify functional mediators of these key cellular programs, we engineer a dual endogenous reporter system by genome-editing the SOX9 and KRT20 loci of human CRC cell lines to express fluorescent reporters, broadcasting aberrant stem cell-like and differentiation activity, respectively. By applying a CRISPR screen targeting 78 epigenetic regulators with 542 sgRNAs to this platform, we identify factors that contribute to stem cell-like activity and differentiation in CRC. Perturbation single cell RNA sequencing (Perturb-seq) of validated hits nominate SMARCB1 of the BAF complex (also known as SWI/SNF) as a negative regulator of differentiation across an array of neoplastic colon models. SMARCB1 is a dependency and required for in vivo growth of human CRC models. These studies highlight the utility of biologically designed endogenous reporter platforms to uncover regulators with therapeutic potential.

Summary: Background: Obesity has been positively associated with most molecular subtypes of colorectal cancer (CRC); however, the magnitude and the causality of these associations is uncertain. Methods: We used Mendelian randomization (MR) to examine potential causal relationships between body size traits (body mass index [BMI], waist circumference, and body fat percentage) with risks of Jass classification types and individual subtypes of CRC (microsatellite instability [MSI] status, CpG island methylator phenotype [CIMP] status, BRAF and KRAS mutations). Summary data on tumour markers were obtained from two genetic consortia (CCFR, GECCO). Findings: A 1-standard deviation (SD:5.1 kg/m2) increment in BMI levels was found to increase risks of Jass type 1MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype (odds ratio [OR]: 2.14, 95% confidence interval [CI]: 1.46, 3.13; p-value = 9 × 10−5) and Jass type 2non-MSI-high,CIMP-high,BRAF-mutated,KRAS-wildtype CRC (OR: 2.20, 95% CI: 1.26, 3.86; p-value = 0.005). The magnitude of these associations was stronger compared with Jass type 4non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-wildtype CRC (p-differences: 0.03 and 0.04, respectively). A 1-SD (SD:13.4 cm) increment in waist circumference increased risk of Jass type 3non-MSI-high,CIMP-low/negative,BRAF-wildtype,KRAS-mutated (OR 1.73, 95% CI: 1.34, 2.25; p-value = 9 × 10−5) that was stronger compared with Jass type 4 CRC (p-difference: 0.03). A higher body fat percentage (SD:8.5%) increased risk of Jass type 1 CRC (OR: 2.59, 95% CI: 1.49, 4.48; p-value = 0.001), which was greater than Jass type 4 CRC (p-difference: 0.03). Interpretation: Body size was more strongly linked to the serrated (Jass types 1 and 2) and alternate (Jass type 3) pathways of colorectal carcinogenesis in comparison to the traditional pathway (Jass type 4). Funding: Cancer Research UK, National Institute for Health Research, Medical Research Council, National Institutes of Health, National Cancer Institute, American Institute for Cancer Research, Brigham and Women's Hospital, Prevent Cancer Foundation, Victorian Cancer Agency, Swedish Research Council, Swedish Cancer Society, Region Västerbotten, Knut and Alice Wallenberg Foundation, Lion’s Cancer Research Foundation, Insamlingsstiftelsen, Umeå University. Full funding details are provided in acknowledgements.

Abstract Routine tumor-node-metastasis (TNM) staging of colorectal cancer is imperfect in predicting survival due to tumor pathobiological heterogeneity and imprecise assessment of tumor spread. We leveraged Bayesian additive regression trees (BART), a statistical learning technique, to comprehensively analyze patient-specific tumor characteristics for the improvement of prognostic prediction. Of 75 clinicopathologic, immune, microbial, and genomic variables in 815 stage II–III patients within two U.S.-wide prospective cohort studies, the BART risk model identified seven stable survival predictors. Risk stratifications (low risk, intermediate risk, and high risk) based on model-predicted survival were statistically significant (hazard ratios 0.19–0.45, vs. higher risk; P

Abstract Histopathologic assessment is indispensable for diagnosing colorectal cancer (CRC). However, manual evaluation of the diseased tissues under the microscope cannot reliably inform patient prognosis or genomic variations crucial for treatment selections. To address these challenges, we develop the Multi-omics Multi-cohort Assessment (MOMA) platform, an explainable machine learning approach, to systematically identify and interpret the relationship between patients’ histologic patterns, multi-omics, and clinical profiles in three large patient cohorts (n = 1888). MOMA successfully predicts the overall survival, disease-free survival (log-rank test P-value

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

The RNF43 G659fs mutation occurs frequently in colorectal cancer, but its function remains poorly understood. In this study, the authors show that RNF43 G659fs is an oncogenic colorectal cancer mutation and sensitizes tumor cells to PI3K/mTOR inhibition.

Germline biallelic pathogenic MUTYH variants predispose patients to colorectal cancer (CRC); however, approaches to identify MUTYH variant carriers are lacking. Here, the authors evaluated mutational signatures that could distinguish MUTYH carriers in large CRC cohorts, and found MUTYH-associated somatic mutations.

Abstract Objectives The CD274 (programmed cell death 1 ligand 1, PD‐L1)/PDCD1 (programmed cell death 1, PD‐1) immune checkpoint axis is known to regulate the antitumor immune response. Evidence also supports an immunosuppressive effect of Fusobacterium nucleatum. We hypothesised that tumor CD274 overexpression might be inversely associated with abundance of F. nucleatum in colorectal carcinoma. Methods We assessed tumor CD274 expression by immunohistochemistry and F. nucleatum DNA within tumor tissue by quantitative PCR in 812 cases among 4465 incident rectal and colon cancer cases that had occurred in two prospective cohort studies. Multivariable logistic regression analyses with inverse probability weighting were used to adjust for selection bias because of tissue data availability and potential confounders including microsatellite instability status, CpG island methylator phenotype, LINE‐1 methylation level and KRAS, BRAF and PIK3CA mutations. Results Fusobacterium nucleatum DNA was detected in tumor tissue in 109 (13%) cases. Tumor CD274 expression level was inversely associated with the amount of F. nucleatum in colorectal cancer tissue (P = 0.0077). For one category‐unit increase in three ordinal F. nucleatum categories (negative vs. low vs. high), multivariable‐adjusted odds ratios (with 95% confidence interval) of the low, intermediate and high CD274 categories (vs. negative) were 0.78 (0.41–1.51), 0.64 (0.32–1.28) and 0.50 (0.25–0.99), respectively (Ptrend = 0.032). Conclusions Tumor CD274 expression level was inversely associated with the amount of F. nucleatum in colorectal cancer tissue, suggesting that different immunosuppressive mechanisms (i.e. PDCD1 immune checkpoint activation and tumor F. nucleatum enrichment) tend to be used by different tumor subgroups.

Abstract Background Certain dietary patterns can elicit systemic and intestinal inflammatory responses, which may influence adaptive anti‐tumor immune responses and tumor behavior. We hypothesized that pro‐inflammatory diets might be associated with higher colorectal cancer mortality and that the association might be stronger for tumors with lower immune responses. Methods We calculated an empirical dietary inflammatory pattern (EDIP) score in 2829 patients among 3988 incident rectal and colon carcinoma cases in the Nurses’ Health Study and Health Professionals Follow‐up Study. Using Cox proportional hazards regression analyses, we examined the prognostic association of EDIP scores and whether it might be modified by histopathologic immune reaction (in 1192 patients with available data). Results Higher EDIP scores after colorectal cancer diagnosis were associated with worse survival, with multivariable‐adjusted hazard ratios (HRs) for the highest versus lowest tertile of 1.41 (95% confidence interval [CI]: 1.13–1.77; Ptrend = 0.003) for 5‐year colorectal cancer‐specific mortality and 1.44 (95% CI, 1.19‐1.74; Ptrend = 0.0004) for 5‐year all‐cause mortality. The association of post‐diagnosis EDIP scores with 5‐year colorectal cancer‐specific mortality differed by degrees of tumor‐infiltrating lymphocytes (TIL; Pinteraction = .002) but not by three other lymphocytic reaction patterns. The multivariable‐adjusted, 5‐year colorectal cancer‐specific mortality HRs for the highest versus lowest EDIP tertile were 1.59 (95% CI: 1.01–2.53) in TIL‐absent/low cases and 0.48 (95% CI: 0.16–1.48) in TIL‐intermediate/high cases. Conclusions Pro‐inflammatory diets after colorectal cancer diagnosis were associated with increased mortality, particularly in patients with absent or low TIL.

Abstract Background Plant‐based foods have been recommended for health. However, not all plant foods are healthy, and little is known about the association between plant‐based diets and specific molecular subtypes of colorectal cancer (CRC). We examined the associations of healthy and unhealthy plant‐based diets with the incidence of CRC and its molecular subtypes. Methods While 123 773 participants of the Nurses’ Health Study and the Health Professionals Follow‐up Study had been followed up (3 143 158 person‐years), 3077 of them had developed CRC. Healthy and unhealthy plant‐based diet indices (hPDI and uPDI, respectively) were calculated using repeated food frequency questionnaire data. We determined the tumoural status of microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and BRAF and KRAS mutations. Results Higher hPDI was associated with lower CRC incidence (multivariable hazard ratio [HR] comparing extreme quartiles, 0.86, 95% confidence interval [CI]: 0.77, 0.96; P‐trend = .04), whereas higher uPDI was associated with higher CRC incidence (multivariable HR comparing extreme quartiles, 1.16, 95% CI: 1.04, 1.29; P‐trend = .005). The association of hPDI significantly differed by KRAS status (P‐heterogeneity = .003) but not by other tumour markers. The hPDI was associated with lower incidence of KRAS‐wildtype CRC (multivariable HR comparing extreme quartiles, 0.74, 95% CI: 0.57, 0.96; P‐trend = .004) but not KRAS‐mutant CRC (P‐trend = .22). Conclusions While unhealthy plant‐based diet enriched with refined grains and sugar is associated with higher CRC incidence, healthy plant‐based diet rich in whole grains, fruits and vegetables is associated with lower incidence of CRC, especially KRAS‐wildtype CRC.

BackgroundThe relationships between tumor stromal features (such as desmoplastic reaction, myxoid stroma, and keloid-like collagen bundles) and immune cells in the colorectal carcinoma microenvironment have not yet been fully characterized.MethodsIn 908 tumors with available tissue among 4,465 incident colorectal adenocarcinoma cases in two prospective cohort studies, we examined desmoplastic reaction, myxoid stroma, and keloid-like collagen bundles. We conducted multiplex immunofluorescence for T cells [CD3, CD4, CD8, CD45RO (PTPRC), and FOXP3] and for macrophages [CD68, CD86, IRF5, MAF, and MRC1 (CD206)]. We used the inverse probability weighting method and the 4,465 incident cancer cases to adjust for selection bias.ResultsImmature desmoplastic reaction was associated with lower densities of intraepithelial CD3+CD8+CD45RO+ cells [multivariable odds ratio (OR) for the highest (vs. lowest) density category, 0.43; 95% confidence interval (CI), 0.29–0.62; Ptrend

Large scale sequencing study is of paramount importance to unravel the heterogeneity of colorectal cancer. Here, the authors sequenced 205 cancer genes in more than 2000 tumours and identified additional mutated driver genes, determined that mutational burden and specific mutations in TP53 are associated with survival odds.

Allelic losses occurring in cancer cells have been suggested as potential targets for therapy. Here, the authors show how recurring loss of heterozygosity of a drug metabolic gene in colorectal cancers can be exploited using a low molecular weight compound.

While advancements in genome sequencing have identified millions of somatic mutations in cancer, their functional impact is poorly understood. We previously developed the expression-based variant impact phenotyping (eVIP) method to use gene expression data to characterize the function of gene variants. The eVIP method uses a decision tree-based algorithm to predict the functional impact of somatic variants by comparing gene expression signatures induced by introduction of wild-type (WT) versus mutant cDNAs in cell lines. The method distinguishes between variants that are gain-of-function, loss-of-function, change-of-function, or neutral. We present eVIP2, software that allows for pathway analysis (eVIP Pathways) and usage with RNA-seq data. To demonstrate the eVIP2 software and approach, we characterized two recurrent frameshift variants in RNF43, a negative regulator of Wnt signaling, frequently mutated in colorectal, gastric, and endometrial cancer. RNF43 WT, RNF43 R117fs, RNF43 G659fs, or GFP control cDNA were overexpressed in HEK293T cells. Analysis with eVIP2 predicted that the frameshift at position 117 was a loss-of-function mutation, as expected. The second frameshift at position 659 has been previously described as a passenger mutation that maintains the RNF43 WT function as a negative regulator of Wnt. Surprisingly, eVIP2 predicted G659fs to be a change-of-function mutation. Additional eVIP Pathways analysis of RNF43 G659fs predicted 10 pathways to be significantly altered, including TNF-α via NFκB signaling, KRAS signaling, and hypoxia, highlighting the benefit of a more comprehensive approach when determining the impact of gene variant function. To validate these predictions, we performed reporter assays and found that each pathway activated by expression of RNF43 G659fs, but not expression of RNF43 WT, was identified as impacted by eVIP2, supporting that RNF43 G659fs is a change-of-function mutation and its effect on the identified pathways. Pathway activation was further validated by Western blot analysis. Lastly, we show primary colon adenocarcinoma patient samples with R117fs and G659fs variants have transcriptional profiles similar to BRAF missense mutations with activated RAS/MAPK signaling, consistent with KRAS signaling pathways being GOF in both variants. The eVIP2 method is an important step towards overcoming the current challenge of variant interpretation in the implementation of precision medicine. eVIP2 is available at https://github.com/BrooksLabUCSC/eVIP2.

Background Myeloid cells represent an abundant yet heterogeneous cell population in the colorectal cancer microenvironment, and their roles remain poorly understood.Methods We used multiplexed immunofluorescence combined with digital image analysis to identify CD14+ monocytic and CD15+ granulocytic cells and to evaluate their maturity (HLA-DR and CD33), immunosuppressive potential (ARG1) and proximity to cytokeratin (KRT)-positive tumor cells in 913 colorectal carcinomas. Using covariate data of 4465 incident colorectal cancers in two prospective cohort studies, the inverse probability weighting method was used with multivariable-adjusted Cox proportional hazards models to assess cancer-specific mortality according to ordinal quartiles (Q1–Q4) of myeloid cell densities. Immune cell–tumor cell proximity was measured with the nearest neighbor method and the G-cross function, which determines the likelihood of any tumor cell having at least one immune cell of the specified type within a certain radius.Results Higher intraepithelial (Ptrend=0.0002; HR for Q4 (vs Q1), 0.48, 95% CI 0.31 to 0.76) and stromal (Ptrend

Immune-checkpoint inhibition therapy has achieved success in a subset of patients. Here the authors profiled about 200 relevant metabolites in patient serum samples from three independent immunotherapy trials and found the serum kynurenine/tryptophan ratio increases to be associated with worse overall survival.

Immunotherapy targeting the CD274 (PD-L1)/PDCD1 (PD-1) immune checkpoint axis has emerged as a promising treatment strategy for various cancers. Experimental evidence suggests that phosphatidylinositol-4,5-bisphosphonate 3-kinase (PI3K) signaling may upregulate CD274 expression. Thus, we hypothesized that PIK3CA mutation, PTEN loss, or their combined status might be associated with CD274 overexpression in colorectal carcinoma. We assessed tumor CD274 and PTEN expression by immunohistochemistry and assessed PIK3CA mutation by pyrosequencing in 753 patients among 4,465 incident rectal and colon cancer cases that had occurred in two U.S.-wide prospective cohort studies. To adjust for potential confounders and selection bias due to tissue availability, inverse probability weighted multivariable ordinal logistic regression analyses used the 4,465 cases and tumoral data including microsatellite instability, CpG island methylator phenotype, KRAS and BRAF mutations. PIK3CA mutation and loss of PTEN expression were detected in 111 of 753 cases (15%) and 342 of 585 cases (58%), respectively. Tumor CD274 expression was negative in 306 (41%), low in 195 (26%), and high in 252 (33%) of 753 cases. PTEN loss was associated with CD274 overexpression [multivariable odds ratio (OR) 1.83; 95% confidence interval (CI), 1.22–2.75; P = .004]. PIK3CA mutation was statistically-insignificantly (P = .036 with the stringent alpha level of 0.005) associated with CD274 overexpression (multivariable OR, 1.54; 95% CI, 1.03–2.31). PIK3CA-mutated PTEN-lost tumors (n = 33) showed higher prevalence of CD274-positivity (82%) than PIK3CA-wild-type PTEN-lost tumors (n = 204; 70% CD274-positivity) and PTEN-expressed tumors (n = 147; 50% CD274-positivity) (P = .003). Our findings support the role of PI3K signaling in the CD274/PDCD1 pathway.

Background CD73–adenosine signaling in the tumor microenvironment is immunosuppressive and may be associated with aggressive renal cell carcinoma (RCC). We investigated the prognostic significance of CD73 protein expression in RCC leveraging nephrectomy samples. We also performed a complementary analysis using The Cancer Genome Atlas (TCGA) dataset to evaluate the correlation of CD73 (ecto-5′-nucleotidase (NT5E), CD39 (ectonucleoside triphosphate diphosphohydrolase 1 (ENTPD1)) and A2 adenosine receptor (A2AR; ADORA2A) transcript levels with markers of angiogenesis and antitumor immune response.Methods Patients with RCC with available archived nephrectomy samples were eligible for inclusion. Tumor CD73 protein expression was assessed by immunohistochemistry and quantified using a combined score (CS: % positive cells×intensity). Samples were categorized as CD73negative (CS=0), CD73low or CD73high (< and ≥median CS, respectively). Multivariable Cox regression analysis compared disease-free survival (DFS) and overall survival (OS) between CD73 expression groups. In the TCGA dataset, samples were categorized as low, intermediate and high NT5E, ENTPD1 and ADORA2A gene expression groups. Gene expression signatures for infiltrating immune cells, angiogenesis, myeloid inflammation, and effector T-cell response were compared between NT5E, ENTPD1 and ADORA2A expression groups.Results Among the 138 patients eligible for inclusion, ‘any’ CD73 expression was observed in 30% of primary tumor samples. High CD73 expression was more frequent in patients with M1 RCC (29% vs 12% M0), grade 4 tumors (27% vs 13% grade 3 vs 15% grades 1 and 2), advanced T-stage (≥T3: 22% vs T2: 19% vs T1: 12%) and tumors with sarcomatoid histology (50% vs 12%). In the M0 cohort (n=107), patients with CD73high tumor expression had significantly worse 5-year DFS (42%) and 10-year OS (22%) compared with those in the CD73negative group (DFS: 75%, adjusted HR: 2.7, 95% CI 1.3 to 5.9, p=0.01; OS: 64%, adjusted HR: 2.6, 95% CI 1.2 to 5.8, p=0.02) independent of tumor stage and grade. In the TCGA analysis, high NT5E expression was associated with significantly worse 5-year OS (p=0.008). NT5E and ENTPD1 expression correlated with higher regulatory T cell (Treg) signature, while ADORA2A expression was associated with increased Treg and angiogenesis signatures.Conclusions High CD73 expression portends significantly worse survival outcomes independent of stage and grade. Our findings provide compelling support for targeting the immunosuppressive and proangiogenic CD73–adenosine pathway in RCC.

Background: Tumour budding and poorly differentiated clusters (PDC) represent forms of tumour invasion. We hypothesised that T-cell densities (reflecting adaptive anti-tumour immunity) might be inversely associated with tumour budding and PDC in colorectal carcinoma. Methods: Utilising 915 colon and rectal carcinomas in two U.S.-wide prospective cohort studies, and multiplex immunofluorescence combined with machine learning algorithms, we assessed CD3, CD4, CD8, CD45RO (PTPRC), and FOXP3 co-expression patterns in lymphocytes. Tumour budding and PDC at invasive fronts were quantified by digital pathology and image analysis using the International tumour Budding Consensus Conference criteria. Using covariate data of 4,420 incident colorectal cancer cases, inverse probability weighting (IPW) was integrated with multivariable logistic regression analysis that assessed the association of T-cell subset densities with tumour budding and PDC while adjusting for selection bias due to tissue availability and potential confounders, including microsatellite instability status. Findings: Tumour budding counts were inversely associated with density of CD3+CD8+ [lowest vs. highest: multivariable odds ratio (OR), 0.50; 95% confidence interval (CI), 0.35–0.70; Ptrend < 0.001] and CD3+CD8+CD45RO+ cells (lowest vs. highest: multivariable OR, 0.44; 95% CI, 0.31–0.63; Ptrend < 0.001) in tumour epithelial region. Tumour budding levels were associated with higher colorectal cancer-specific mortality (multivariable hazard ratio, 2.13; 95% CI, 1.57–2.89; Ptrend < 0.001) in Cox regression analysis. There were no significant associations of PDC with T-cell subsets. Interpretation: Tumour epithelial naïve and memory cytotoxic T cell densities are inversely associated with tumour budding at invasive fronts, suggesting that cytotoxic anti-tumour immunity suppresses tumour microinvasion.

Evidence indicates the pathogenic role of epigenetic alterations in early-onset colorectal cancers diagnosed before age 50. However, features of colorectal cancers diagnosed at age 50–54 (hereafter referred to as “intermediate-onset”) remain less known. We hypothesized that tumor long interspersed nucleotide element-1 (LINE-1) hypomethylation might be increasingly more common with decreasing age of colorectal cancer diagnosis. In 1356 colorectal cancers, including 28 early-onset and 66 intermediate-onset cases, the tumor LINE-1 methylation level measured by bisulfite-PCR-pyrosequencing (scaled 0 to 100) showed a mean of 63.6 (standard deviation (SD) 10.1). The mean tumor LINE-1 methylation level decreased with decreasing age (mean 64.7 (SD 10.4) in age ≥70, 62.8 (SD 9.4) in age 55–69, 61.0 (SD 10.2) in age 50–54, and 58.9 (SD 12.0) in age p < 0.0001). In linear regression analysis, the multivariable-adjusted β coefficient (95% confidence interval (CI)) (vs. age ≥70) was −1.38 (−2.47 to −0.30) for age 55–69, −2.82 (−5.29 to −0.34) for age 50–54, and −4.54 (−8.24 to −0.85) for age Ptrend = 0.0003). Multivariable-adjusted hazard ratios (95% CI) for LINE-1 methylation levels of ≤45, 45–55, and 55–65 (vs. >65) were 2.33 (1.49–3.64), 1.39 (1.05–1.85), and 1.29 (1.02–1.63), respectively (Ptrend = 0.0005). In conclusion, tumor LINE-1 hypomethylation is increasingly more common with decreasing age of colorectal cancer diagnosis, suggesting a role of global DNA hypomethylation in colorectal cancer arising in younger adults.

Large-scale genomic characterization of tumors from prospective cohort studies may yield new insights into cancer pathogenesis. We performed whole-exome sequencing of 619 incident colorectal cancers (CRCs) and integrated the results with tumor immunity, pathology, and survival data. We identified recurrently mutated genes in CRC, such as BCL9L, RBM10, CTCF, and KLF5, that were not previously appreciated in this disease. Furthermore, we investigated the genomic correlates of immune-cell infiltration and found that higher neoantigen load was positively associated with overall lymphocytic infiltration, tumor-infiltrating lymphocytes (TILs), memory T cells, and CRC-specific survival. The association with TILs was evident even within microsatellite-stable tumors. We also found positive selection of mutations in HLA genes and other components of the antigen-processing machinery in TIL-rich tumors. These results may inform immunotherapeutic approaches in CRC. More generally, this study demonstrates a framework for future integrative molecular epidemiology research in colorectal and other malignancies.

Inhibitors targeting the PDCD1 (programmed cell death 1, PD-1) immune checkpoint pathway have revolutionized cancer treatment strategies. The TIME (Tumor Immunity in the MicroEnvironment) classification based on tumor CD274 (PDCD1 ligand 1, PD-L1) expression and tumor-infiltrating lymphocytes (TIL) has been proposed to predict response to immunotherapy. It remains to be determined clinical, pathological, and molecular features of TIME subtypes of colorectal cancer. Using 812 colon and rectal carcinoma cases from the Nurses' Health Study and Health Professionals Follow-up Study, we examined the association of tumor characteristics and survival outcomes with four TIME subtypes (TIME 1, CD274low/TILabsent; TIME 2, CD274high/TILpresent; TIME 3, CD274low/TILpresent; and TIME 4, CD274high/TILabsent). In survival analyses, Cox proportional hazards models were adjusted for potential confounders, including microsatellite instability (MSI) status, CpG island methylator phenotype (CIMP) status, LINE-1 methylation level, and KRAS, BRAF, and PIK3CA mutation status. TIME subtypes 1, 2, 3 and 4 had 218 (27%), 117 (14%), 103 (13%), and 374 (46%) colorectal cancer cases, respectively. Compared with TIL-absent subtypes (TIME 1 and 4), TIL-present subtypes (TIME 2 and 3) were associated with high-level MSI, high-degree CIMP, BRAF mutation, and higher amounts of neoantigens (p < 0.001). TIME subtypes were not significantly associated with colorectal cancer-specific or overall survival. In conclusion, TIL-present TIME subtypes of colorectal cancer are associated with high levels of MSI and neoantigen load, supporting better responsiveness to cancer immunotherapy. Further studies examining tumor molecular alterations and additional factors in the tumor microenvironment may inform development of immunoprevention and immunotherapy strategies.