The genomic landscape of 2,023 colorectal cancers.

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

Stratification to Neoadjuvant Radiotherapy in Rectal Cancer by Regimen and Transcriptional Signatures.

Response to neoadjuvant radiotherapy (RT) in rectal cancer has been associated with immune and stromal features that are captured by transcriptional signatures. However, how such associations perform across different chemoradiotherapy regimens and within individual consensus molecular subtypes (CMS) and how they affect survival remain unclear. In this study, gene expression and clinical data of pretreatment biopsies from nine cohorts of primary rectal tumors were combined (N = 826). Exploratory analyses were done with transcriptomic signatures for the endpoint of pathologic complete response (pCR), considering treatment regimen or CMS subtype. Relevant findings were tested for overall survival and recurrence-free survival. Immune and stromal signatures were strongly associated with pCR and lack of pCR, respectively, in RT and capecitabine (Cap)/5-fluorouracil (5FU)-treated patients (N = 387), in which the radiosensitivity signature (RSS) showed the strongest association. Upon addition of oxaliplatin (Ox; N = 123), stromal signatures switched direction and showed higher chances to achieve pCR than without Ox (p for interaction 0.02). Among Cap/5FU patients, most signatures performed similarly across CMS subtypes, except cytotoxic lymphocytes that were associated with pCR in CMS1 and CMS4 cases compared with other CMS subtypes (p for interaction 0.04). The only variables associated with survival were pCR and RSS. Although the frequency of pCR across different chemoradiation regimens is relatively similar, our data suggest that response rates may differ depending on the biological landscape of rectal cancer. Response to neoadjuvant RT in stroma-rich tumors may potentially be improved by the addition of Ox. RSS in preoperative biopsies provides predictive information for response specifically to neoadjuvant RT with 5FU. SIGNIFICANCE: Rectal cancers with stromal features may respond better to RT and 5FU/Cap with the addition of Ox. Within patients not treated with Ox, high levels of cytotoxic lymphocytes associate with response only in immune and stromal tumors. Our analyses provide biological insights about the outcome by different radiotherapy regimens in rectal cancer.

Identification and validation of a machine learning model of complete response to radiation in rectal cancer reveals immune infiltrate and TGFß as key predictors.

BACKGROUND: It is uncertain which biological features underpin the response of rectal cancer (RC) to radiotherapy. No biomarker is currently in clinical use to select patients for treatment modifications. METHODS: We identified two cohorts of patients (total N = 249) with RC treated with neoadjuvant radiotherapy (45Gy/25) plus fluoropyrimidine. This discovery set included 57 cases with pathological complete response (pCR) to chemoradiotherapy (23%). Pre-treatment cancer biopsies were assessed using transcriptome-wide mRNA expression and targeted DNA sequencing for copy number and driver mutations. Biological candidate and machine learning (ML) approaches were used to identify predictors of pCR to radiotherapy independent of tumour stage. Findings were assessed in 107 cases from an independent validation set (GSE87211). FINDINGS: Three gene expression sets showed significant independent associations with pCR: Fibroblast-TGFß Response Signature (F-TBRS) with radioresistance; and cytotoxic lymphocyte (CL) expression signature and consensus molecular subtype CMS1 with radiosensitivity. These associations were replicated in the validation cohort. In parallel, a gradient boosting machine model comprising the expression of 33 genes generated in the discovery cohort showed high performance in GSE87211 with 90% sensitivity, 86% specificity. Biological and ML signatures indicated similar mechanisms underlying radiation response, and showed better AUC and p-values than published transcriptomic signatures of radiation response in RC. INTERPRETATION: RCs responding completely to chemoradiotherapy (CRT) have biological characteristics of immune response and absence of immune inhibitory TGFß signalling. These tumours may be identified with a potential biomarker based on a 33 gene expression signature. This could help select patients likely to respond to treatment with a primary radiotherapy approach as for anal cancer. Conversely, those with predicted radioresistance may be candidates for clinical trials evaluating addition of immune-oncology agents and stromal TGFß signalling inhibition. FUNDING: The Stratification in Colorectal Cancer Consortium (S:CORT) was funded by the Medical Research Council and Cancer Research UK (MR/M016587/1).

Prognostic and Predictive Value of Immunoscore in Stage III Colorectal Cancer: Pooled Analysis of Cases From the SCOT and IDEA-HORG Studies.

PURPOSE: Immunoscore (IS) is prognostic in stage III colorectal cancer (CRC) and may predict benefit of duration (6 v 3 months) of adjuvant infusional fluorouracil, leucovorin, and oxaliplatin (FOLFOX) chemotherapy. We sought to determine IS prognostic and predictive value in stage-III CRC treated with adjuvant FOLFOX or oral capecitabine and infusional oxaliplatin (CAPOX) in the SCOT and IDEA-HORG trials. METHODS: Three thousand sixty-one cases had tumor samples, of which 2,643 (1,792 CAPOX) were eligible for IS testing. Predefined cutoffs (IS-Low and IS-High) were used to classify cases into two groups for analysis of disease-free survival (3-year DFS) and multivariable-adjusted hazard ratios (mvHRs) by Cox regression. RESULTS: IS was determined in 2,608 (99.5%) eligible cases, with 877 (33.7%) samples classified as IS-Low. IS-Low tumors were more commonly high-risk (T4 and/or N2; 52.9% IS-Low v 42.2% IS-High; P < .001) and in younger patients (P = .024). Patients with IS-Low tumors had significantly shorter DFS in the CAPOX, FOLFOX, and combined cohorts (mvHR, 1.52 [95% CI, 1.28 to 1.82]; mvHR, 1.58 [95% CI, 1.22 to 2.04]; and mvHR, 1.55 [95% CI, 1.34 to 1.79], respectively; P < .001 all comparisons), regardless of sex, BMI, clinical risk group, tumor location, treatment duration, or chemotherapy regimen. IS prognostic value was greater in younger (≤65 years) than older (>65 years) patients in the CAPOX cohort (mvHR, 1.92 [95% CI, 1.50 to 2.46] v 1.28 [95% CI, 1.01 to 1.63], PINTERACTION = .026), and in DNA mismatch repair proficient than deficient mismatch repair disease (mvHR, 1.68 [95% CI, 1.41 to 2.00] v 0.67 [95% CI, 0.30 to 1.49], PINTERACTION = .03), although these exploratory analyses were uncorrected for multiple testing. Adding IS to a model containing all clinical variables significantly improved prediction of DFS (likelihood ratio test, P < .001) regardless of MMR status. CONCLUSION: IS is prognostic in stage III CRC treated with FOLFOX or CAPOX, including within clinically relevant tumor subgroups. Possible variation in IS prognostic value by age and MMR status, and prediction of benefit from extended adjuvant therapy merit validation.

JAK/STAT3 represents a therapeutic target for colorectal cancer patients with stromal-rich tumors.

Colorectal cancer (CRC) is a heterogenous malignancy underpinned by dysregulation of cellular signaling pathways. Previous literature has implicated aberrant JAK/STAT3 signal transduction in the development and progression of solid tumors. In this study we investigate the effectiveness of inhibiting JAK/STAT3 in diverse CRC models, establish in which contexts high pathway expression is prognostic and perform in depth analysis underlying phenotypes. In this study we investigated the use of JAK inhibitors for anti-cancer activity in CRC cell lines, mouse model organoids and patient-derived organoids. Immunohistochemical staining of the TransSCOT clinical trial cohort, and 2 independent large retrospective CRC patient cohorts was performed to assess the prognostic value of JAK/STAT3 expression. We performed mutational profiling, bulk RNASeq and NanoString GeoMx® spatial transcriptomics to unravel the underlying biology of aberrant signaling. Inhibition of signal transduction with JAK1/2 but not JAK2/3 inhibitors reduced cell viability in CRC cell lines, mouse, and patient derived organoids (PDOs). In PDOs, reduced Ki67 expression was observed post-treatment. A highly significant association between high JAK/STAT3 expression within tumor cells and reduced cancer-specific survival in patients with high stromal invasion (TSPhigh) was identified across 3 independent CRC patient cohorts, including the TrasnSCOT clinical trial cohort. Patients with high phosphorylated STAT3 (pSTAT3) within the TSPhigh group had higher influx of CD66b + cells and higher tumoral expression of PDL1. Bulk RNAseq of full section tumors showed enrichment of NFκB signaling and hypoxia in these cases. Spatial deconvolution through GeoMx® demonstrated higher expression of checkpoint and hypoxia-associated genes in the tumor (pan-cytokeratin positive) regions, and reduced lymphocyte receptor signaling in the TME (pan-cytokeratin- and αSMA-) and αSMA (pan-cytokeratin- and αSMA +) areas. Non-classical fibroblast signatures were detected across αSMA + regions in cases with high pSTAT3. Therefore, in this study we have shown that inhibition of JAK/STAT3 represents a promising therapeutic strategy for patients with stromal-rich CRC tumors. High expression of JAK/STAT3 proteins within both tumor and stromal cells predicts poor outcomes in CRC, and aberrant signaling is associated with distinct spatially-dependant differential gene expression.

Multiplex analysis of intratumoural immune infiltrate and prognosis in patients with stage II-III colorectal cancer from the SCOT and QUASAR 2 trials: a retrospective analysis.

BACKGROUND: Tumour-infiltrating CD8+ cytotoxic T cells confer favourable prognosis in colorectal cancer. The added prognostic value of other infiltrating immune cells is unclear and so we sought to investigate their prognostic value in two large clinical trial cohorts. METHODS: We used multiplex immunofluorescent staining of tissue microarrays to assess the densities of CD8+, CD20+, FoxP3+, and CD68+ cells in the intraepithelial and intrastromal compartments from tumour samples of patients with stage II-III colorectal cancer from the SCOT trial (ISRCTN59757862), which examined 3 months versus 6 months of adjuvant oxaliplatin-based chemotherapy, and from the QUASAR 2 trial (ISRCTN45133151), which compared adjuvant capecitabine with or without bevacizumab. Both trials included patients aged 18 years or older with an Eastern Cooperative Oncology Group performance status of 0-1. Immune marker predictors were analysed by multiple regression, and the prognostic and predictive values of markers for colorectal cancer recurrence-free interval by Cox regression were assessed using the SCOT cohort for discovery and QUASAR 2 cohort for validation. FINDINGS: After exclusion of cases without tissue microarrays and with technical failures, and following quality control, we included 2340 cases from the SCOT trial and 1069 from the QUASAR 2 trial in our analysis. Univariable analysis of associations with recurrence-free interval in cases from the SCOT trial showed a strong prognostic value of intraepithelial CD8 (CD8IE) as a continuous variable (hazard ratio [HR] for 75th vs 25th percentile [75vs25] 0·73 [95% CI 0·68-0·79], p=2·5 × 10-16), and of intrastromal FoxP3 (FoxP3IS; 0·71 [0·64-0·78], p=1·5 × 10-13) but not as strongly in the epithelium (FoxP3IE; 0·89 [0·84-0·96], p=1·5 × 10-4). Associations of other markers with recurrence-free interval were moderate. CD8IE and FoxP3IS retained independent prognostic value in bivariable and multivariable analysis, and, compared with either marker alone, a composite marker including both markers (CD8IE-FoxP3IS) was superior when assessed as a continuous variable (adjusted [a]HR75 vs 25 0·70 [95% CI 0·63-0·78], p=5·1 × 10-11) and when categorised into low, intermediate, and high density groups using previously published cutpoints (aHR for intermediate vs high 1·68 [95% CI 1·29-2·20], p=1·3 × 10-4; low vs high 2·58 [1·91-3·49], p=7·9 × 10-10), with performance similar to the gold-standard Immunoscore. The prognostic value of CD8IE-FoxP3IS was confirmed in cases from the QUASAR 2 trial, both as a continuous variable (aHR75 vs 25 0·84 [95% CI 0·73-0·96], p=0·012) and as a categorical variable for low versus high density (aHR 1·80 [95% CI 1·17-2·75], p=0·0071) but not for intermediate versus high (1·30 [0·89-1·88], p=0·17). INTERPRETATION: Combined evaluation of CD8IE and FoxP3IS could help to refine risk stratification in colorectal cancer. Investigation of FoxP3IS cells as an immunotherapy target in colorectal cancer might be merited. FUNDING: Medical Research Council, National Institute for Health Research, Cancer Research UK, Swedish Cancer Society, Roche, and Promedica Foundation.

Intra-prostatic tumour evolution, steps in metastatic spread and histogenomic associations revealed by integration of multi-region whole-genome sequencing with histopathological features.

BACKGROUND: Extension of prostate cancer beyond the primary site by local invasion or nodal metastasis is associated with poor prognosis. Despite significant research on tumour evolution in prostate cancer metastasis, the emergence and evolution of cancer clones at this early stage of expansion and spread are poorly understood. We aimed to delineate the routes of evolution and cancer spread within the prostate and to seminal vesicles and lymph nodes, linking these to histological features that are used in diagnostic risk stratification. METHODS: We performed whole-genome sequencing on 42 prostate cancer samples from the prostate, seminal vesicles and lymph nodes of five treatment-naive patients with locally advanced disease. We spatially mapped the clonal composition of cancer across the prostate and the routes of spread of cancer cells within the prostate and to seminal vesicles and lymph nodes in each individual by analysing a total of > 19,000 copy number corrected single nucleotide variants. RESULTS: In each patient, we identified sample locations corresponding to the earliest part of the malignancy. In patient 10, we mapped the spread of cancer from the apex of the prostate to the seminal vesicles and identified specific genomic changes associated with the transformation of adenocarcinoma to amphicrine morphology during this spread. Furthermore, we show that the lymph node metastases in this patient arose from specific cancer clones found at the base of the prostate and the seminal vesicles. In patient 15, we observed increased mutational burden, altered mutational signatures and histological changes associated with whole genome duplication. In all patients in whom histological heterogeneity was observed (4/5), we found that the distinct morphologies were located on separate branches of their respective evolutionary trees. CONCLUSIONS: Our results link histological transformation with specific genomic alterations and phylogenetic branching. These findings have implications for diagnosis and risk stratification, in addition to providing a rationale for further studies to characterise the genetic changes causally linked to morphological transformation. Our study demonstrates the value of integrating multi-region sequencing with histopathological data to understand tumour evolution and identify mechanisms of prostate cancer spread.

Utilizing Human Genetics to Develop Chemoprevention for Cancer-Too Good an Opportunity to be Missed.

Large-scale genetic studies are reliably identifying many risk factors for disease in the general population. Several of these genetic risk factors encode potential drug targets, and genetics has already helped to introduce targeted agents for some diseases, an example being lipid-lowering drugs to reduce the incidence of cardiovascular disease. Multiple drugs have been developed to treat cancers based on somatic mutations and genomics, but in stark contrast, there seems to be a reluctance to use germline genetic data to develop drugs to prevent malignancy, despite the large numbers of people who could benefit, the potential for lowering cancer rates, and the widespread current use of non-pharmaceutical measures to reduce cancer risk factors such as tobacco, alcohol, and infectious diseases. We argue that concerted efforts for cancer prevention based on genetics, including genes influenced by common polymorphisms that modulate cancer risk, are urgently needed. There are enormous, yet underutilized, opportunities to develop novel targeted agents for chemoprevention of cancer based on human germline genetics. Such efforts are likely to require the support of a dedicated funding program by national and international agencies. See related commentary by Winham and Sherman, p. 13.

Adult stem cell activity in naked mole rats for long-term tissue maintenance.

The naked mole rat (NMR), Heterocephalus glaber, the longest-living rodent, provides a unique opportunity to explore how evolution has shaped adult stem cell (ASC) activity and tissue function with increasing lifespan. Using cumulative BrdU labelling and a quantitative imaging approach to track intestinal ASCs (Lgr5+) in their native in vivo state, we find an expanded pool of Lgr5+ cells in NMRs, and these cells specifically at the crypt base (Lgr5+CBC) exhibit slower division rates compared to those in short-lived mice but have a similar turnover as human LGR5+CBC cells. Instead of entering quiescence (G0), NMR Lgr5+CBC cells reduce their division rates by prolonging arrest in the G1 and/or G2 phases of the cell cycle. Moreover, we also observe a higher proportion of differentiated cells in NMRs that confer enhanced protection and function to the intestinal mucosa which is able to detect any chemical imbalance in the luminal environment efficiently, triggering a robust pro-apoptotic, anti-proliferative response within the stem/progenitor cell zone.

Estudio de la relación de la ancestría genética con el riesgo al cáncer colorrectal en Colombia / Study of the relationship between genetic ancestry and colorrectal cáncer risk in Colombia

Introducción. El cáncer colorrectal es una carga para la salud pública en Colombia y el mundo. Estudios de asociación genética han identificado regiones cromosómicas asociadas a esta enfermedad, mostrando riesgo variable entre poblaciones, debido a la historia demográfica y la ancestría genética. Objetivo. Estudiar el riesgo que aportan 20 marcadores al cáncer colorrectal en Colombia, empleando 955 casos y 972 controles del consorcio CHIBCHA, analizando conjuntamente el efecto de la ancestría genética global y local. Metodología. Las muestras se genotipificaron usando microarreglos Axyom Affymetrix LAT y CUSTOME, para obtener los genotipos genómicos globales, incluyendo 20 SNPs de riesgo. Los análisis estadísticos se realizaron en PLINK (asociaciones), ADMIXTURE (ancestría global), Elai (ancestría local) y R (modelos logísticos). Resultados. Once regiones cromosómicas resultaron asociadas presentando ORs entre 1.14 y 1.41 (p<0.05): 18q21.1, 19q13.11, 10p14, 14q.2.2, 20p12.3, 8q23.3, 6p21.2, 15q13.3 y 8q24.21. Una mayor ancestría europea se asoció con el riesgo a nivel global (OR=3.016, IC 95%:1.162-7.894, p=0.00325), y a nivel cromosómico local se detectaron las regiones 6q23.2 (ORajustado=1.378, IC95%: 1.202-1.580, Pajustado=4.2e-6) y 4p13 (ORajustado=1.301, IC95%:1.137-1.489; Pajustado=0.00013). Conclusiones. La ancestría podría considerarse un factor en la explicación de la susceptibilidad en Colombia, indicando que la mezcla genética de origen amerindio y europeo, influye en la estructura poblacional y explicaría las diferencias en la incidencia del CCR entre poblaciones latinas y europeas.

Introduction: Colorectal cancer is a public health burden in the world and Colombia. Recent genome wide association studies have identified chromosomal regions associated with the disease, depicting variable risk between populations, owing to the demographic history and genetic ancestry. Objective: We aimed to study the colorectal cancer risk in Colombia provided for 20 genetic markers, by using 955 cases and 972 controls from the CHIBCHA consortium, in the context of global and local genetic ancestry. Methodology: The samples were genotyped using Axyom Affymetrix LAT and CUSTOME array in order to obtain the global genome genotypes including 20 risk SNPs. Statistical analysis was performed in PLINK (associations), ADMIXTURE (global ancestry), Elai (local ancestry) and R language (logistic models). Results: Eleven chromosomal regions were associated with ORs ranging between 1.14-1.41 (p<0.05): 18q21.1, 19q13.11, 10p14, 14q.2.2, 20p12.3, 8q23.3, 6p21.2, 15q13.3 y 8q24.21. On average, a higher global European ancestry was associated with colorectal cancer risk (OR=3.016, IC 95%:1.162-7.894, p=0.00325). At the local chromosomal level two regions presented a significant increment of European ancestry 6q23.2 (OR adjusted=1.378, CI95%: 1.202-1.580, p adjusted =4.2e-6) and 4p13 (OR adjusted =1.301, CI95%:1.137-1.489; p adjusted =0.00013). Conclusions: Genetic ancestry can be considered as a relevant factor for the colorectal cancer susceptibility in Colombia. Both Native American and European ancestry are accounting for the most part of population structure in the sample we studied, which could explain the differences for the colorectal cancer incidence between Latin American and European populations.