A Genomic Analysis Workflow for Colorectal Cancer Precision Oncology

Transcriptome-wide gene expression outlier analysis pinpoints therapeutic vulnerabilities in colorectal cancer

Multiple strategies are continuously being explored to expand the drug target repertoire in solid tumors. We devised a novel computational workflow for transcriptome-wide gene expression outlier analysis that allows the systematic identification of both overexpression and underexpression events in cancer cells. Here, it was applied to expression values obtained through RNA sequencing in 226 colorectal cancer (CRC) cell lines that were also characterized by whole-exome sequencing and microarray-based DNA methylation profiling. We found cell models displaying an abnormally high or low expression level for 3533 and 965 genes, respectively. Gene expression abnormalities that have been previously associated with clinically relevant features of CRC cell lines were confirmed. Moreover, by integrating multi-omics data, we identified both genetic and epigenetic alternations underlying outlier expression values. Importantly, our atlas of CRC gene expression outliers can guide the discovery of novel drug targets and biomarkers. As a proof of concept, we found that CRC cell lines lacking expression of the MTAP gene are sensitive to treatment with a PRMT5-MTA inhibitor (MRTX1719). Finally, other tumor types may also benefit from this approach.

Mutational signatures of colorectal cancers according to distinct computational workflows

Tumor mutational signatures have gained prominence in cancer research, yet the lack of standardized methods hinders reproducibility and robustness. Leveraging colorectal cancer (CRC) as a model, we explored the influence of computational parameters on mutational signature analyses across 230 CRC cell lines and 152 CRC patients. Results were validated in three independent datasets: 483 endometrial cancer patients stratified by mismatch repair (MMR) status, 35 lung cancer patients by smoking status and 12 patient-derived organoids (PDOs) annotated for colibactin exposure. Assessing various bioinformatic tools, reference datasets and input data sizes including whole genome sequencing, whole exome sequencing and a pan-cancer gene panel, we demonstrated significant variability in the results. We report that the use of distinct algorithms and references led to statistically different results, highlighting how arbitrary choices may induce variability in the mutational signature contributions. Furthermore, we found a differential contribution of mutational signatures between coding and intergenic regions and defined the minimum number of somatic variants required for reliable mutational signature assignment. To facilitate the identification of the most suitable workflows, we developed Comparative Mutational Signature analysis on Coding and Extragenic Regions (CoMSCER), a bioinformatic tool which allows researchers to easily perform comparative mutational signature analysis by coupling the results from several tools and public reference datasets and to assess mutational signature contributions in coding and non-coding genomic regions. In conclusion, our study provides a comparative framework to elucidate the impact of distinct computational workflows on mutational signatures.

Dendritic cell-targeted therapy expands CD8 T cell responses to bona-fide neoantigens in lung tumors

Cross-presentation by type 1 DCs (cDC1) is critical to induce and sustain antitumoral CD8 T cell responses to model antigens, in various tumor settings. However, the impact of cross-presenting cDC1 and the potential of DC-based therapies in tumors carrying varied levels of bona-fide neoantigens (neoAgs) remain unclear. Here we develop a hypermutated model of non-small cell lung cancer in female mice, encoding genuine MHC-I neoepitopes to study neoAgs-specific CD8 T cell responses in spontaneous settings and upon Flt3L + αCD40 (DC-therapy). We find that cDC1 are required to generate broad CD8 responses against a range of diverse neoAgs. DC-therapy promotes immunogenicity of weaker neoAgs and strongly inhibits the growth of high tumor-mutational burden (TMB) tumors. In contrast, low TMB tumors respond poorly to DC-therapy, generating mild CD8 T cell responses that are not sufficient to block progression. scRNA transcriptional analysis, immune profiling and functional assays unveil the changes induced by DC-therapy in lung tissues, which comprise accumulation of cDC1 with increased immunostimulatory properties and less exhausted effector CD8 T cells. We conclude that boosting cDC1 activity is critical to broaden the diversity of anti-tumoral CD8 T cell responses and to leverage neoAgs content for therapeutic advantage.

Non-canonical antigens are the largest fraction of peptides presented by MHC class I in mismatch repair deficient murine colorectal cancer

BACKGROUND

Immunotherapy based on checkpoint inhibitors is highly effective in mismatch repair deficient (MMRd) colorectal cancer (CRC). These tumors carry a high number of mutations, which are predicted to translate into a wide array of neoepitopes; however, a systematic classification of the neoantigen repertoire in MMRd CRC is lacking. Mass spectrometry peptidomics has demonstrated the existence of MHC class I associated peptides (MAPs) originating from non-coding DNA regions. Based on these premises we investigated DNA genomic regions responsible for generating MMRd-induced peptides.

METHODS

We exploited mouse CRC models in which the MMR gene Mlh1 was genetically inactivated. Isogenic cell lines CT26 Mlh1+/+ and Mlh1-/- were inoculated in immunocompromised and immunocompetent mice. Whole genome and RNA sequencing data were generated from samples obtained before and after injection in murine hosts. First, peptide databases were built from transcriptomes of isogenic cell lines. We then compiled a database of peptides lost after tumor cells injection in immunocompetent mice, likely due to immune editing. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) and matched next-generation sequencing databases were employed to identify the DNA regions from which the immune-targeted MAPs originated. Finally, we adopted in vitro T cell assays to verify whether MAP-specific T cells were part of the in vivo immune response against Mlh1-/- cells.

RESULTS

Whole genome sequencing analyses revealed an unbalanced distribution of immune edited alterations across the genome in Mlh1-/- cells grown in immunocompetent mice. Specifically, untranslated (UTR) and coding regions exhibited the largest fraction of mutations leading to highly immunogenic peptides. Moreover, the integrated computational and LC-MS/MS analyses revealed that MAPs originate mainly from atypical translational events in both Mlh1+/+ and Mlh1-/- tumor cells. In addition, mutated MAPs-derived from UTRs and out-of-frame translation of coding regions-were highly enriched in Mlh1-/- cells. The MAPs trigger T-cell activation in mice primed with Mlh1-/- cells.

CONCLUSIONS

Our results suggest that-in comparison to MMR proficient CRC-MMRd tumors generate a significantly higher number of non-canonical mutated peptides able to elicit T cell responses. These results reveal the importance of evaluating the diversity of neoepitope repertoire in MMRd tumors.

Coexisting cancer stem cells with heterogeneous gene amplifications, transcriptional profiles, and malignancy are isolated from single glioblastomas

Glioblastoma (GBM) is known as an intractable, highly heterogeneous tumor encompassing multiple subclones, each supported by a distinct glioblastoma stem cell (GSC). The contribution of GSC genetic and transcriptional heterogeneity to tumor subclonal properties is debated. In this study, we describe the systematic derivation, propagation, and characterization of multiple distinct GSCs from single, treatment-naive GBMs (GSC families). The tumorigenic potential of each GSC better correlates with its transcriptional profile than its genetic make-up, with classical GSCs being inherently more aggressive and mesenchymal more dependent on exogenous growth factors across multiple GBMs. These GSCs can segregate and recapitulate different histopathological aspects of the same GBM, as shown in a paradigmatic tumor with two histopathologically distinct components, including a conventional GBM and a more aggressive primitive neuronal component. This study provides a resource for investigating how GSCs with distinct genetic and/or phenotypic features contribute to individual GBM heterogeneity and malignant escalation.

Consensus clustering methodology to improve molecular stratification of non-small cell lung cancer

Recent advances in machine learning research, combined with the reduced sequencing costs enabled by modern next-generation sequencing, paved the way to the implementation of precision medicine through routine multi-omics molecular profiling of tumours. Thus, there is an emerging need of reliable models exploiting such data to retrieve clinically useful information. Here, we introduce an original consensus clustering approach, overcoming the intrinsic instability of common clustering methods based on molecular data. This approach is applied to the case of non-small cell lung cancer (NSCLC), integrating data of an ongoing clinical study (PROMOLE) with those made available by The Cancer Genome Atlas, to define a molecular-based stratification of the patients beyond, but still preserving, histological subtyping. The resulting subgroups are biologically characterized by well-defined mutational and gene-expression profiles and are significantly related to disease-free survival (DFS). Interestingly, it was observed that (1) cluster B, characterized by a short DFS, is enriched in KEAP1 and SKP2 mutations, that makes it an ideal candidate for further studies with inhibitors, and (2) over- and under-representation of inflammation and immune systems pathways in squamous-cell carcinomas subgroups could be potentially exploited to stratify patients treated with immunotherapy.

Molecular methods for colorectal cancer screening: Progress with next-generation sequencing evolution

Currently, colorectal cancer (CRC) represents the third most common malignancy and the second most deadly cancer worldwide, with a higher incidence in developed countries. Like other solid tumors, CRC is a heterogeneous genomic disease in which various alterations, such as point mutations, genomic rearrangements, gene fusions or chromosomal copy number alterations, can contribute to the disease development. However, because of its orderly natural history, easily accessible onset location and high lifetime incidence, CRC is ideally suited for preventive intervention, but the many screening efforts of the last decades have been compromised by performance limitations and low penetrance of the standard screening tools. The advent of next-generation sequencing (NGS) has both facilitated the identification of previously unrecognized CRC features such as its relationship with gut microbial pathogens and revolutionized the speed and throughput of cataloguing CRC-related genomic alterations. Hence, in this review, we summarized the several diagnostic tools used for CRC screening in the past and the present, focusing on recent NGS approaches and their revolutionary role in the identification of novel genomic CRC characteristics, the advancement of understanding the CRC carcinogenesis and the screening of clinically actionable targets for personalized medicine.

Targeting the DNA Damage Response Pathways and Replication Stress in Colorectal Cancer

PURPOSE

Genomic instability is a hallmark of cancer and targeting DNA damage response (DDR) is emerging as a promising therapeutic strategy in different solid tumors. The effectiveness of targeting DDR in colorectal cancer has not been extensively explored.

EXPERIMENTAL DESIGN

We challenged 112 cell models recapitulating the genomic landscape of metastatic colorectal cancer with ATM, ATR, CHK1, WEE1, and DNA-PK inhibitors, in parallel with chemotherapeutic agents. We focused then on ATR inhibitors (ATRi) and, to identify putative biomarkers of response and resistance, we analyzed at multiple levels colorectal cancer models highly sensitive or resistant to these drugs.

RESULTS

We found that around 30% of colorectal cancers, including those carrying KRAS and BRAF mutations and unresponsive to targeted agents, are sensitive to at least one DDR inhibitor. By investigating potential biomarkers of response to ATRi, we found that ATRi-sensitive cells displayed reduced phospho-RPA32 foci at basal level, while ATRi-resistant cells showed increased RAD51 foci formation in response to replication stress. Lack of ATM and RAD51C expression was associated with ATRi sensitivity. Analysis of mutational signatures and HRDetect score identified a subgroup of ATRi-sensitive models. Organoids derived from patients with metastatic colorectal cancer recapitulated findings obtained in cell lines.

CONCLUSIONS

In conclusion, a subset of colorectal cancers refractory to current therapies could benefit from inhibitors of DDR pathways and replication stress. A composite biomarker involving phospho-RPA32 and RAD51 foci, lack of ATM and RAD51C expression, as well as analysis of mutational signatures could be used to identify colorectal cancers likely to respond to ATRi.

DNA damage response and repair genes in advanced bone and soft tissue sarcomas: An 8-gene signature as a candidate predictive biomarker of response to trabectedin and olaparib combination

Background

Advanced and unresectable bone and soft tissue sarcomas (BSTS) still represent an unmet medical need. We demonstrated that the alkylating agent trabectedin and the PARP1-inhibitor olaparib display antitumor activity in BSTS preclinical models. Moreover, in a phase Ib clinical trial (NCT02398058), feasibility, tolerability and encouraging results have been observed and the treatment combination is currently under study in a phase II trial (NCT03838744).

Methods

Differential expression of genes involved in DNA Damage Response and Repair was evaluated by Nanostring^® technology, extracting RNA from pre-treatment tumor samples of 16 responder (≥6-month progression free survival) and 16 non-responder patients. Data validation was performed by quantitative real-time PCR, RNA in situ hybridization, and immunohistochemistry. The correlation between the identified candidate genes and both progression-free survival and overall survival was investigated in the publicly available dataset “Sarcoma (TCGA, The Cancer Genome Atlas)”.

Results

Differential RNA expression analysis revealed an 8-gene signature (CDKN2A, PIK3R1, SLFN11, ATM, APEX2, BLM, XRCC2, MAD2L2) defining patients with better outcome upon trabectedin+olaparib treatment. In responder vs. non-responder patients, a significant differential expression of these genes was further confirmed by RNA in situ hybridization and by qRT-PCR and immunohistochemistry in selected experiments. Correlation between survival outcomes and genetic alterations in the identified genes was shown in the TCGA sarcoma dataset.

Conclusions

This work identified an 8-gene expression signature to improve prediction of response to trabectedin+olaparib combination in BSTS. The predictive role of these potential biomarkers warrants further investigation.

Temozolomide Treatment Alters Mismatch Repair and Boosts Mutational Burden in Tumor and Blood of Colorectal Cancer Patients

The majority of metastatic colorectal cancers (mCRC) are mismatch repair (MMR) proficient and unresponsive to immunotherapy, whereas MMR-deficient (MMRd) tumors often respond to immune-checkpoint blockade. We previously reported that the treatment of colorectal cancer preclinical models with temozolomide (TMZ) leads to MMR deficiency, increased tumor mutational burden (TMB), and sensitization to immunotherapy. To clinically translate these findings, we designed the ARETHUSA clinical trial whereby O6-methylguanine-DNA-methyltransferase (MGMT)-deficient, MMR-proficient, RAS-mutant mCRC patients received priming therapy with TMZ. Analysis of tissue biopsies and circulating tumor DNA (ctDNA) revealed the emergence of a distinct mutational signature and increased TMB after TMZ treatment. Multiple alterations in the nucleotide context favored by the TMZ signature emerged in MMR genes, and the p.T1219I MSH6 variant was detected in ctDNA and tissue of 94% (16/17) of the cases. A subset of patients whose tumors displayed the MSH6 mutation, the TMZ mutational signature, and increased TMB achieved disease stabilization upon pembrolizumab treatment.

SIGNIFICANCE

MMR-proficient mCRCs are unresponsive to immunotherapy. We provide the proof of concept that inactivation of MMR genes can be achieved pharmacologically with TMZ and molecularly monitored in the tissue and blood of patients with mCRC. This strategy deserves additional evaluation in mCRC patients whose tumors are no longer responsive to standard-of-care treatments. See related commentary by Willis and Overman, p. 1612. This article is highlighted in the In This Issue feature, p. 1599.

Identification of four novel hub genes as monitoring biomarkers for colorectal cancer

Background

It must be admitted that the incidence of colorectal cancer (CRC) was on the rise all over the world, but the related treatment had not caught up. Further research on the underlying pathogenesis of CRC was conducive to improving the survival status of current CRC patients.

Methods

Differentially expressed genes (DEGs) screening were conducted based on “limma” and “RobustRankAggreg” package of R software. Weighted gene co-expression network analysis (WGCNA) was performed in the integrated DEGs that from The Cancer Genome Atlas (TCGA), and all samples of validation were from Gene Expression Omnlbus (GEO) dataset.

Results

The terms obtained in the functional annotation for primary DEGs indicated that they were associated with CRC. The MEyellow stand out whereby showed the significant correlation with clinical feature (disease), and 4 hub genes, including ABCC13, AMPD1, SCNN1B and TMIGD1, were identified in yellow module. Nine datasets from Gene Expression Omnibus database confirmed these four genes were significantly down-regulated and the survival estimates for the low-expression group of these genes were lower than for the high-expression group in Kaplan-Meier survival analysis section. MEXPRESS suggested that down-regulation of some top hub genes may be caused by hypermethylation. Receiver operating characteristic curves indicated that these genes had certain diagnostic efficacy. Moreover, tumor-infiltrating immune cells and gene set enrichment analysis for hub genes suggested that there were some associations between these genes and the pathogenesis of CRC.

Conclusion

This study identified modules that were significantly associated with CRC, four novel hub genes, and further analysis of these genes. This may provide a little new insights and directions into the potential pathogenesis of CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s41065-021-00216-7.

Circulating tumor DNA to guide rechallenge with panitumumab in metastatic colorectal cancer: the phase 2 CHRONOS trial

Anti-epidermal growth factor receptor (EGFR) monoclonal antibodies are approved for the treatment of RAS wild-type (WT) metastatic colorectal cancer (mCRC), but the emergence of resistance mutations restricts their efficacy. We previously showed that RAS, BRAF and EGFR mutant alleles, which appear in circulating tumor DNA (ctDNA) during EGFR blockade, decline upon therapy withdrawal. We hypothesized that monitoring resistance mutations in blood could rationally guide subsequent therapy with anti-EGFR antibodies. We report here the results of CHRONOS, an open-label, single-arm phase 2 clinical trial exploiting blood-based identification of RAS/BRAF/EGFR mutations levels to tailor a chemotherapy-free anti-EGFR rechallenge with panitumumab (ClinicalTrials.gov: NCT03227926 ; EudraCT 2016-002597-12). The primary endpoint was objective response rate. Secondary endpoints were progression-free survival, overall survival, safety and tolerability of this strategy. In CHRONOS, patients with tissue-RAS WT tumors after a previous treatment with anti-EGFR-based regimens underwent an interventional ctDNA-based screening. Of 52 patients, 16 (31%) carried at least one mutation conferring resistance to anti-EGFR therapy and were excluded. The primary endpoint of the trial was met; and, of 27 enrolled patients, eight (30%) achieved partial response and 17 (63%) disease control, including two unconfirmed responses. These clinical results favorably compare with standard third-line treatments and show that interventional liquid biopsies can be effectively and safely exploited in a timely manner to guide anti-EGFR rechallenge therapy with panitumumab in patients with mCRC. Further larger and randomized trials are warranted to formally compare panitumumab rechallenge with standard-of-care therapies in this patient setting.

PARP1 Inhibitor and Trabectedin Combination Does Not Increase Tumor Mutational Burden in Advanced Sarcomas-A Preclinical and Translational Study

Simple Summary

Immunotherapy has revolutionized cancer treatment, but not for all tumor types. Indeed, sarcomas are considered “immune-cold” tumors, which are relatively unresponsive to immunotherapy. One strategy to potentiate immunotherapy efficacy is to increase tumor immunogenicity, for instance by boosting the number of candidate targets (neoantigens) to be recognized by the immune system. Tumor mutational burden indicates the number of somatic mutations identified in the tumor and normalized per megabase. Tumor mutational burden is considered as an acceptable, measurable surrogate of tumor neoantigens. Here, we explored whether the combination of two DNA-damaging agents, trabectedin and olaparib, could increase tumor mutational burden in sarcomas, to prime subsequent immunotherapy. We found no variation in tumor mutational burden after trabectedin + olaparib in preclinical and clinical samples. Therefore, other aspects should be considered to increase sarcoma immunogenicity, by exploiting different pathways such as the potential modulation of the tumor microenvironment induced by trabectedin + olaparib.

Abstract

Drug-induced tumor mutational burden (TMB) may contribute to unleashing the immune response in relatively “immune-cold” tumors, such as sarcomas. We previously showed that PARP1 inhibition perpetuates the DNA damage induced by the chemotherapeutic agent trabectedin in both preclinical models and sarcoma patients. In the present work, we explored acquired genetic changes in DNA repair genes, mutational signatures, and TMB in a translational platform composed of cell lines, xenografts, and tumor samples from patients treated with trabectedin and olaparib combination, compared to cells treated with temozolomide, an alkylating agent that induces hypermutation. Whole-exome and targeted panel sequencing data analyses revealed that three cycles of trabectedin and olaparib combination neither affected the mutational profiles, DNA repair gene status, or copy number alterations, nor increased TMB both in homologous recombinant-defective and proficient cells or in xenografts. Moreover, TMB was not increased in tumor specimens derived from trabectedin- and olaparib-treated patients (5–6 cycles) when compared to pre-treatment biopsies. Conversely, repeated treatments with temozolomide induced a massive TMB increase in the SJSA-1 osteosarcoma model. In conclusion, a trabectedin and olaparib combination did not show mutagenic effects and is unlikely to prime subsequent immune-therapeutic interventions based on TMB increase. On the other hand, these findings are reassuring in the increasing warning of treatment-induced hematologic malignancies correlated to PARP1 inhibitor use.

Werner Helicase Is a Synthetic-Lethal Vulnerability in Mismatch Repair-Deficient Colorectal Cancer Refractory to Targeted Therapies, Chemotherapy, and Immunotherapy

Targeted therapies, chemotherapy, and immunotherapy are used to treat patients with mismatch repair-deficient (dMMR)/microsatellite instability-high (MSI-H) colorectal cancer. The clinical effectiveness of targeted therapy and chemotherapy is limited by resistance and drug toxicities, and about half of patients receiving immunotherapy have disease that is refractory to immune checkpoint inhibitors. Loss of Werner syndrome ATP-dependent helicase (WRN) is a synthetic lethality in dMMR/MSI-H cells. To inform the development of WRN as a therapeutic target, we performed WRN knockout or knockdown in 60 heterogeneous dMMR colorectal cancer preclinical models, demonstrating that WRN dependency is an almost universal feature and a robust marker for patient selection. Furthermore, models of resistance to clinically relevant targeted therapy, chemotherapy, and immunotherapy retain WRN dependency. These data show the potential of therapeutically targeting WRN in patients with dMMR/MSI-H colorectal cancer and support WRN as a therapeutic option for patients with dMMR/MSI-H cancers refractory to current treatment strategies. SIGNIFICANCE: We found that a large, diverse set of dMMR/MSI-H colorectal cancer preclinical models, including models of treatment-refractory disease, are WRN-dependent. Our results support WRN as a promising synthetic-lethal target in dMMR/MSI-H colorectal cancer tumors as a monotherapy or in combination with targeted agents, chemotherapy, or immunotherapy.This article is highlighted in the In This Issue feature, p. 1861.

CD4 T Cell-Dependent Rejection of Beta-2 Microglobulin Null Mismatch Repair-Deficient Tumors

Inactivation of beta-2 microglobulin (B2M) is considered a determinant of resistance to immune checkpoint inhibitors (ICPi) in melanoma and lung cancers. In contrast, B2M loss does not appear to affect response to ICPis in mismatch repair-deficient (MMRd) colorectal tumors where biallelic inactivation of B2M is frequently observed. We inactivated B2m in multiple murine MMRd cancer models. Although MMRd cells would not readily grow in immunocompetent mice, MMRd B2m null cells were tumorigenic and regressed when treated with anti-PD-1 and anti-CTLA4. The efficacy of ICPis against MMRd B2m null tumors did not require CD8⁺ T cells but relied on the presence of CD4⁺ T cells. Human tumors expressing low levels of B2M display increased intratumoral CD4⁺ T cells. We conclude that B2M inactivation does not blunt the efficacy of ICPi in MMRd tumors, and we identify a unique role for CD4⁺ T cells in tumor rejection. SIGNIFICANCE: B2M alterations, which impair antigen presentation, occur frequently in microsatellite-unstable colorectal cancers. Although in melanoma and lung cancers B2M loss is a mechanism of resistance to immune checkpoint blockade, we show that MMRd tumors respond to ICPis through CD4⁺ T-cell activation.This article is highlighted in the In This Issue feature, p. 1601.

Radiogenomics in Colorectal Cancer

Simple Summary

Colorectal carcinoma is characterized by intratumoral heterogeneity that can be assessed by radiogenomics. Radiomics, high-throughput quantitative data extracted from medical imaging, combined with molecular analysis, through genomic and transcriptomic data, is expected to lead to significant advances in personalized medicine. However, a radiogenomics approach in colorectal cancer is still in its early stages and many problems remain to be solved. Here we review the progress and challenges in this field at its current stage, as well as future developments.

Abstract

The steady improvement of high-throughput technologies greatly facilitates the implementation of personalized precision medicine. Characterization of tumor heterogeneity through image-derived features—radiomics and genetic profile modifications—genomics, is a rapidly evolving field known as radiogenomics. Various radiogenomics studies have been dedicated to colorectal cancer so far, highlighting the potential of these approaches to enhance clinical decision-making. In this review, a general outline of colorectal radiogenomics literature is provided, discussing the current limitations and suggested further developments.

Omics Integration Analysis Unravel the Landscape of Driving Mechanisms of Colorectal Cancer

Colorectal cancer (CRC) is one of the most malignant cancers and results in a substantial rate of morbidity and mortality. Diagnosis of this malignancy in early stages increases the chance of effective treatment. High-throughput data analyses reveal omics signatures and also provide the possibility of developing computational models for early detection of this disease. Such models would be able to use as complementary tools for early detection of different types of cancers including CRC. In this study, using gene expression data, the Flux balance analysis (FBA) applied to decode metabolic fluxes in cancer and normal cells. Moreover, transcriptome and genome analyses revealed driver agents of CRC in a biological network scheme. By applying comprehensive publicly available data from TCGA, different aspect of CRC regulome including the regulatory effect of gene expression, methylation, microRNA, copy number aberration and point mutation profile over protein levels investigated and the results provide a regulatory picture underlying CRC. Compiling omics profiles indicated snapshots of changes in different omics levels and flux rate of CRC. In conclusion, considering obtained CRC signatures and their role in biological operating systems of cells, the results suggest reliable driver regulatory modules that could potentially serve as biomarkers and therapeutic targets and furthermore expand our understanding of driving mechanisms of this disease.
.

Whole exome sequencing analysis of urine trans-renal tumour DNA in metastatic colorectal cancer patients

Background

The analysis of circulating free tumour DNA (ctDNA) in blood, commonly referred as liquid biopsy, is being used to characterise patients with solid cancers. Tumour-specific genetic variants can also be present in DNA isolated from other body fluids, such as urine. Unlike blood, urine sampling is non-invasive, can be self-performed, and allows recurrent longitudinal monitoring. The features of tumour DNA that clears from the glomerular filtration barrier, named trans-renal tumour DNA (trtDNA), are largely unexplored.

Patients and methods

Specimens were collected from 24 patients with KRAS or BRAF mutant metastatic colorectal cancer (mCRC). Driver mutations were assessed by droplet digital PCR (ddPCR) in ctDNA from plasma and trtDNA from urine. Whole exome sequencing (WES) was performed in DNA isolated from tissue, plasma and urine.

Results

Out of the 24 CRC cases, only four had sufficient DNA to allow WES analyses in urine and plasma. We found that tumour alterations primarily reside in low molecular weight fragments (less than 112 bp). In patients whose trtDNA was more than 2.69% of the urine derived DNA, cancer-specific molecular alterations, mutational signatures and copy number profiles identified in urine DNA are comparable with those detected in plasma ctDNA.

Conclusions

With current technologies, WES analysis of trtDNA is feasible in a small fraction of mCRC patients. Tumour-related genetic information is mainly present in low molecular weight DNA fragments. Although the limited amounts of trtDNA poses analytical challenges, enrichment of low molecular weight DNAs and optimised computational tools can improve the detection of tumour-specific genetic information in urine.

A Subset of Colorectal Cancers with Cross-Sensitivity to Olaparib and Oxaliplatin

PURPOSE

Defects in the homologous recombination (HR) repair pathway are of clinical interest due to sensitivity of HR-deficient cells to PARP inhibitors. We were interested in defining PARP vulnerability in patients with metastatic colorectal cancer (mCRC) carrying KRAS and BRAF mutations who display poor prognosis, have limited therapeutic options, and represent an unmet clinical need.

EXPERIMENTAL DESIGN

We tested colorectal cancer cell lines, patient-derived organoids (PDO), and patient-derived xenografts (PDX) enriched for KRAS and BRAF mutations for sensitivity to the PARP inhibitor olaparib, and the chemotherapeutic agents oxaliplatin and 5-fluorouracil (5-FU). Genomic profiles and DNA repair proficiency of colorectal cancer models were compared with pharmacologic response.

RESULTS

Thirteen of 99 (around 13%) colorectal cancer cell lines were highly sensitive to clinically active concentrations of olaparib and displayed functional deficiency in HR. Response to PARP blockade was positively correlated with sensitivity to oxaliplatin in colorectal cancer cell lines as well as patient-derived organoids. Treatment of PDXs with olaparib impaired tumor growth and maintenance therapy with PARP blockade after initial oxaliplatin response delayed disease progression in mice.

CONCLUSIONS

These results indicate that a colorectal cancer subset characterized by poor prognosis and limited therapeutic options is vulnerable to PARP inhibition and suggest that PDO-based drug-screening assays can be used to identify patients with colorectal cancer likely to benefit from olaparib. As patients with mCRC almost invariably receive therapies based on oxaliplatin, "maintenance" treatment with PARP inhibitors warrants further clinical investigation.

Adaptive mutability of colorectal cancers in response to targeted therapies

The emergence of drug resistance limits the efficacy of targeted therapies in human tumors. The prevalent view is that resistance is a fait accompli: when treatment is initiated, cancers already contain drug-resistant mutant cells. Bacteria exposed to antibiotics transiently increase their mutation rates (adaptive mutability), thus improving the likelihood of survival. We investigated whether human colorectal cancer (CRC) cells likewise exploit adaptive mutability to evade therapeutic pressure. We found that epidermal growth factor receptor (EGFR)/BRAF inhibition down-regulates mismatch repair (MMR) and homologous recombination DNA-repair genes and concomitantly up-regulates error-prone polymerases in drug-tolerant (persister) cells. MMR proteins were also down-regulated in patient-derived xenografts and tumor specimens during therapy. EGFR/BRAF inhibition induced DNA damage, increased mutability, and triggered microsatellite instability. Thus, like unicellular organisms, tumor cells evade therapeutic pressures by enhancing mutability.

Clinical Application of Next-Generation Sequencing as A Liquid Biopsy Technique in Advanced Colorectal Cancer: A Trick or A Treat?

Owing to its advantages over prior relevant technologies, massive parallel or next-generation sequencing (NGS) is rapidly evolving, with growing applications in a wide range of human diseases. The burst in actionable molecular alterations in many cancer types advocates for the practicality of using NGS in the clinical setting, as it permits the parallel characterization of multiple genes in a cost- and time-effective way, starting from low-input DNA. In advanced clinical practice, the oncological management of colorectal cancer requires prior knowledge of KRAS, NRAS, and BRAF status, for the design of appropriate therapeutic strategies, with more gene mutations still surfacing as potential biomarkers. Tumor heterogeneity, as well as the need for serial gene profiling due to tumor evolution and the emergence of novel genetic alterations, have promoted the use of liquid biopsies-especially in the form of circulating tumor DNA (ctDNA)-as a promising alternative to tissue molecular analysis. This review discusses recent studies that have used plasma NGS in advanced colorectal cancer and summarizes the clinical applications, as well as the technical challenges involved in adopting this technique in a clinically beneficial oncological practice.

How liquid biopsies can change clinical practice in oncology

Cell-free DNA fragments are shed into the bloodstream by tumor cells. The analysis of circulating tumor DNA (ctDNA), commonly known as liquid biopsy, can be exploited for a variety of clinical applications. ctDNA is being used to genotype solid cancers non-invasively, to track tumor dynamics and to detect the emergence of drug resistance. In a few settings, liquid biopsies have already entered clinical practice. For example, ctDNA is used to guide treatment in a subset of lung cancers. In this review, we discuss how recent improvements in the sensitivity and accuracy of ctDNA analyses have led to unprecedented advances in this research field. We further consider what is required for the routine deployment of liquid biopsies in the clinical diagnostic space. We pinpoint technical hurdles that liquid biopsies have yet to overcome, including preanalytical and analytical challenges. We foresee how liquid biopsies will transform clinical practice: by complementing (or replacing) imaging to monitor treatment response and by detecting minimal residual disease after surgery with curative intent.

Patient-Derived Xenografts and Matched Cell Lines Identify Pharmacogenomic Vulnerabilities in Colorectal Cancer

PURPOSE

Patient-derived xenograft (PDX) models accurately recapitulate the tumor of origin in terms of histopathology, genomic landscape, and therapeutic response, but some limitations due to costs associated with their maintenance and restricted amenability for large-scale screenings still exist. To overcome these issues, we established a platform of 2D cell lines (xeno-cell lines, XL), derived from PDXs of colorectal cancer with matched patient germline gDNA available.

EXPERIMENTAL DESIGN

Whole-exome and transcriptome sequencing analyses were performed. Biomarkers of response and resistance to anti-HER therapy were annotated. Dependency on the WRN helicase gene was assessed in MSS, MSI-H, and MSI-like XLs using a reverse genetics functional approach.

RESULTS

XLs recapitulated the entire spectrum of colorectal cancer transcriptional subtypes. Exome and RNA-seq analyses delineated several molecular biomarkers of response and resistance to EGFR and HER2 blockade. Genotype-driven responses observed in vitro in XLs were confirmed in vivo in the matched PDXs. MSI-H models were dependent upon WRN gene expression, while loss of WRN did not affect MSS XLs growth. Interestingly, one MSS XL with transcriptional MSI-like traits was sensitive to WRN depletion.

CONCLUSIONS

The XL platform represents a preclinical tool for functional gene validation and proof-of-concept studies to identify novel druggable vulnerabilities in colorectal cancer.

Evolving neoantigen profiles in colorectal cancers with DNA repair defects

Background

Neoantigens that arise as a consequence of tumor-specific mutations can be recognized by T lymphocytes leading to effective immune surveillance. In colorectal cancer (CRC) and other tumor types, a high number of neoantigens is associated with patient response to immune therapies. The molecular processes governing the generation of neoantigens and their turnover in cancer cells are poorly understood. We exploited CRC as a model system to understand how alterations in DNA repair pathways modulate neoantigen profiles over time.

Methods

We performed whole exome sequencing (WES) and RNA sequencing (RNAseq) in CRC cell lines, in vitro and in vivo, and in CRC patient-derived xenografts (PDXs) to track longitudinally genomic profiles, clonal evolution, mutational signatures, and predicted neoantigens.

Results

The majority of CRC models showed remarkably stable mutational and neoantigen profiles; however, those carrying defects in DNA repair genes continuously diversified. Rapidly evolving and evolutionary stable CRCs displayed characteristic genomic signatures and transcriptional profiles. Downregulation of molecules implicated in antigen presentation occurred selectively in highly mutated and rapidly evolving CRC.

Conclusions

These results indicate that CRCs carrying alterations in DNA repair pathways display dynamic neoantigen patterns that fluctuate over time. We define CRC subsets characterized by slow and fast evolvability and link this phenotype to downregulation of antigen-presenting cellular mechanisms. Longitudinal monitoring of the neoantigen landscape could be relevant in the context of precision medicine.

Electronic supplementary material

The online version of this article (10.1186/s13073-019-0654-6) contains supplementary material, which is available to authorized users.