Error-Corrected Deep Targeted Sequencing of Circulating Cell-Free DNA from Colorectal Cancer Patients for Sensitive Detection of Circulating Tumor DNA.

Circulating tumor DNA (ctDNA) is a promising biomarker, reflecting the presence of tumor cells. Sequencing-based detection of ctDNA at low tumor fractions is challenging due to the crude error rate of sequencing. To mitigate this challenge, we developed ultra-deep mutation-integrated sequencing (UMIseq), a fixed-panel deep targeted sequencing approach, which is universally applicable to all colorectal cancer (CRC) patients. UMIseq features UMI-mediated error correction, the exclusion of mutations related to clonal hematopoiesis, a panel of normal samples for error modeling, and signal integration from single-nucleotide variations, insertions, deletions, and phased mutations. UMIseq was trained and independently validated on pre-operative (pre-OP) plasma from CRC patients (n = 364) and healthy individuals (n = 61). UMIseq displayed an area under the curve surpassing 0.95 for allele frequencies (AFs) down to 0.05%. In the training cohort, the pre-OP detection rate reached 80% at 95% specificity, while it was 70% in the validation cohort. UMIseq enabled the detection of AFs down to 0.004%. To assess the potential for detection of residual disease, 26 post-operative plasma samples from stage III CRC patients were analyzed. From this we found that the detection of ctDNA was associated with recurrence. In conclusion, UMIseq demonstrated robust performance with high sensitivity and specificity, enabling the detection of ctDNA at low allele frequencies.

Circulating tumor DNA for prognosis assessment and postoperative management after curative-intent resection of colorectal liver metastases.

The recurrence rate of colorectal liver metastases (CRLM) patients treated with curative intent is above 50%. Standard of care surveillance includes intensive computed tomographic (CT) imaging as well as carcinoembryonic antigen (CEA) measurements. Nonetheless, relapse detection often happens too late to resume curative treatment. This longitudinal cohort study enrolled 115 patients with plasma samples (N = 439) prospectively collected before surgery, postoperatively at day 30 and every third month for up to 3 years. Droplet digital PCR (ddPCR) was used to monitor serial plasma samples for somatic mutations. Assessment of ctDNA status either immediately after surgery, or serially during surveillance, stratified the patients into groups of high and low recurrence risk (hazard ratio [HR], 7.6; 95% CI, 3.0-19.7; P < .0001; and HR, 4.3; 95% CI, 2.3-8.1; P < .0001, respectively). The positive predictive value (PPV) of ctDNA was 100% in all postoperative analyses. In multivariable analyses, postoperative ctDNA status was the only consistently significant risk marker associated with relapse (P < .0001). Indeterminate CT findings were observed for 30.8% (21/68) of patients. All patients (9/21) that were ctDNA positive at the time of the indeterminate CT scan later relapsed, contrasting 42.6% (5/12) of those ctDNA negative (P = .0046). Recurrence diagnoses in patients with indeterminate CT findings were delayed (median 2.8 months, P < .0001). ctDNA status is strongly associated with detection of minimal residual disease and early detection of relapse. Furthermore, ctDNA status can potentially contribute to clinical decision-making in case of indeterminate CT findings, reducing time-to-intervention.

Validation of computational determination of microsatellite status using whole exome sequencing data from colorectal cancer patients.

BACKGROUND: Microsatellite instability (MSI), resulting from a defective mismatch repair system, occurs in approximately 15% of sporadic colorectal cancers (CRC). Since MSI is associated with a poor response to 5-fluorouracile based chemotherapy and is a positive predictive marker of immunotherapy, it is routine practice to evaluate the MSI status of resected tumors in CRC patients. MSIsensor is a novel computational tool for determining MSI status using Next Generation Sequencing. However, it is not widely used in the clinic and has not been independently validated in exome data from CRC. To facilitate clinical implementation of computational determination of MSI status, we compared MSIsensor to current gold standard methods for MSI testing. METHODS: MSI status was determined for 130 CRC patients (UICC stage I-IV) using immunohistochemistry, PCR based microsatellite stability testing and by applying MSIsensor to exome sequenced tumors and paired germline DNA. Furthermore, we investigated correlation between MSI status, mutational load and mutational signatures. RESULTS: Eighteen out of 130 (13.8%) patients were microsatellite instable. We found a 100% agreement between MSIsensor and gold standard methods for MSI testing. All MSI tumors were hypermutated. In addition, two microsatellite stable (MSS) tumors were hypermutated, which was explained by a dominant POLE signature and pathogenic POLE mutations (p.Pro286Arg and p.Ser459Phe). CONCLUSION: MSIsensor is a robust tool, which can be used to determine MSI status of tumor samples from exome sequenced CRC patients.

The potential role of Alu Y in the development of resistance to SN38 (Irinotecan) or oxaliplatin in colorectal cancer.

BACKGROUND: Irinotecan (SN38) and oxaliplatin are chemotherapeutic agents used in the treatment of colorectal cancer. However, the frequent development of resistance to these drugs represents a considerable challenge in the clinic. Alus as retrotransposons comprise 11% of the human genome. Genomic toxicity induced by carcinogens or drugs can reactivate Alus by altering DNA methylation. Whether or not reactivation of Alus occurs in SN38 and oxaliplatin resistance remains unknown. RESULTS: We applied reduced representation bisulfite sequencing (RRBS) to investigate the DNA methylome in SN38 or oxaliplatin resistant colorectal cancer cell line models. Moreover, we extended the RRBS analysis to tumor tissue from 14 patients with colorectal cancer who either did or did not benefit from capecitabine + oxaliplatin treatment. For the clinical samples, we applied a concept of 'DNA methylation entropy' to estimate the diversity of DNA methylation states of the identified resistance phenotype-associated methylation loci observed in the cell line models. We identified different loci being characteristic for the different resistant cell lines. Interestingly, 53% of the identified loci were Alu sequences- especially the Alu Y subfamily. Furthermore, we identified an enrichment of Alu Y sequences that likely results from increased integration of new copies of Alu Y sequence in the drug-resistant cell lines. In the clinical samples, SOX1 and other SOX gene family members were shown to display variable DNA methylation states in their gene regions. The Alu Y sequences showed remarkable variation in DNA methylation states across the clinical samples. CONCLUSION: Our findings imply a crucial role of Alu Y in colorectal cancer drug resistance. Our study underscores the complexity of colorectal cancer aggravated by mobility of Alu elements and stresses the importance of personalized strategies, using a systematic and dynamic view, for effective cancer therapy.

Whole-genome Mutational Analysis for Tumor-informed Detection of Circulating Tumor DNA in Patients with Urothelial Carcinoma.

BACKGROUND AND OBJECTIVE: Circulating tumor DNA (ctDNA) can be used for sensitive detection of minimal residual disease (MRD). However, the probability of detecting ctDNA in settings of low tumor burden is limited by the number of mutations analyzed and the plasma volume available. We used a whole-genome sequencing (WGS) approach for ctDNA detection in patients with urothelial carcinoma. METHODS: We used a tumor-informed WGS approach for ctDNA-based detection of MRD and evaluation of treatment responses. We analyzed 916 longitudinally collected plasma samples from 112 patients with localized muscle-invasive bladder cancer who received neoadjuvant chemotherapy (NAC) before radical cystectomy. Recurrence-free survival (primary endpoint), overall survival, and ctDNA dynamics during NAC were assessed. KEY FINDINGS AND LIMITATIONS: We found that WGS-based ctDNA detection is prognostic for patient outcomes with a median lead time of 131 d over radiographic imaging. WGS-based ctDNA assessment after radical cystectomy identified recurrence with sensitivity of 91% and specificity of 92%. In addition, genomic characterization of post-treatment plasma samples with a high ctDNA level revealed acquisition of platinum therapy-associated mutational signatures and copy number variations not present in the primary tumors. The sequencing depth is a limitation for studying tumor evolution. CONCLUSIONS AND CLINICAL IMPLICATIONS: Our results support the use of WGS for ultrasensitive ctDNA detection and highlight the possibility of plasma-based tracking of tumor evolution. WGS-based ctDNA detection represents a promising option for clinical use owing to the low volume of plasma needed and the ease of performing WGS, eliminating the need for personalized assay design. PATIENT SUMMARY: Detection of tumor DNA in blood samples from patients with cancer of the urinary tract is associated with poorer outcomes. Disease recurrence after surgery can be identified by the presence of tumor DNA in blood before it can be detected on radiography scans.

Ultrasensitive plasma-based monitoring of tumor burden using machine-learning-guided signal enrichment.

In solid tumor oncology, circulating tumor DNA (ctDNA) is poised to transform care through accurate assessment of minimal residual disease (MRD) and therapeutic response monitoring. To overcome the sparsity of ctDNA fragments in low tumor fraction (TF) settings and increase MRD sensitivity, we previously leveraged genome-wide mutational integration through plasma whole-genome sequencing (WGS). Here we now introduce MRD-EDGE, a machine-learning-guided WGS ctDNA single-nucleotide variant (SNV) and copy-number variant (CNV) detection platform designed to increase signal enrichment. MRD-EDGESNV uses deep learning and a ctDNA-specific feature space to increase SNV signal-to-noise enrichment in WGS by ~300× compared to previous WGS error suppression. MRD-EDGECNV also reduces the degree of aneuploidy needed for ultrasensitive CNV detection through WGS from 1 Gb to 200 Mb, vastly expanding its applicability within solid tumors. We harness the improved performance to identify MRD following surgery in multiple cancer types, track changes in TF in response to neoadjuvant immunotherapy in lung cancer and demonstrate ctDNA shedding in precancerous colorectal adenomas. Finally, the radical signal-to-noise enrichment in MRD-EDGESNV enables plasma-only (non-tumor-informed) disease monitoring in advanced melanoma and lung cancer, yielding clinically informative TF monitoring for patients on immune-checkpoint inhibition.

Analysis of circulating tumor DNA during checkpoint inhibition in metastatic melanoma using a tumor-agnostic panel.

Immunotherapy has revolutionized treatment of patients diagnosed with metastatic melanoma, where nearly half of patients receive clinical benefit. However, immunotherapy is also associated with immune-related adverse events, which may be severe and persistent. It is therefore important to identify patients that do not benefit from therapy early. Currently, regularly scheduled CT scans are used to investigate size changes in target lesions to evaluate progression and therapy response. This study aims to explore if panel-based analysis of circulating tumor DNA (ctDNA) taken at 3-week intervals may provide a window into the growing cancer, can be used to identify nonresponding patients early, and determine genomic alterations associated with acquired resistance to checkpoint immunotherapy without analysis of tumor tissue biopsies. We designed a gene panel for ctDNA analysis and sequenced 4-6 serial plasma samples from 24 patients with unresectable stage III or IV melanoma and treated with first-line checkpoint inhibitors enrolled at the Department of Oncology at Aarhus University Hospital, Denmark. TERT was the most mutated gene found in ctDNA and associated with a poor prognosis. We detected more ctDNA in patients with high metastatic load, which indicates that more aggressive tumors release more ctDNA into the bloodstream. Although we did not find evidence of specific mutations associated with acquired resistance, we did demonstrate in this limited cohort of 24 patients that untargeted, panel-based ctDNA analysis has the potential to be used as a minimally invasive tool in clinical practice to identify patients where the benefits of immunotherapy outweigh the drawbacks.

Impact of Whole Genome Doubling on Detection of Circulating Tumor DNA in Colorectal Cancer.

OBJECTIVE: Circulating tumor DNA (ctDNA) is a candidate biomarker of cancer with practice-changing potential in the detection of both early and residual disease. Disease stage and tumor size affect the probability of ctDNA detection, whereas little is known about the influence of other tumor characteristics on ctDNA detection. This study investigates the impact of tumor cell whole-genome doubling (WGD) on the detection of ctDNA in plasma collected preoperatively from newly diagnosed colorectal cancer (CRC) patients. METHODS: WGD was estimated from copy numbers derived from whole-exome sequencing (WES) data of matched tumor and normal DNA from 833 Danish CRC patients. To explore if tumor WGD status impacts ctDNA detection, we applied tumor-informed ctDNA analysis to preoperative plasma samples from all patients. RESULTS: Patients with WGD+ tumors had 53% increased odds of being ctDNA positive (OR = 1.53, 95%CI: 1.12-2.09). After stratification for UICC stage, the association persisted for Stage I (OR = 2.44, 95%CI: 1.22-5.03) and Stage II (OR = 1.76, 95%CI: 1.11-2.81) but not for Stage III (OR = 0.83, 95%CI: 0.44-1.53) patients. CONCLUSION: The presence of WGD significantly increases the probability of detecting ctDNA, particularly for early-stage disease. In patients with more advanced disease, the benefit of WGD on ctDNA detection is less pronounced, consistent with increased DNA shedding from these tumors, making ctDNA detection less dependent on the amount of ctDNA released per tumor cell.

Comparing single-target and multitarget approaches for postoperative circulating tumour DNA detection in stage II-III colorectal cancer patients.

Circulating tumour DNA (ctDNA) detection for postoperative risk stratification in cancer patients has great clinical potential. However, low ctDNA abundances complicates detection. Multitarget (MT) detection strategies have been developed to increase sensitivity. Yet, empirical evidence supporting performance gains of MT vs. single-target (ST) strategies in a postoperative setting is limited. We compared ctDNA detection in 379 paired plasma samples from 112 stage II-III colorectal cancer patients by ST digital PCR and MT sequencing of 16 patient-specific variants. The strategies exhibited good concordance (90%, Cohen's Kappa 0.79), with highly correlated ctDNA quantifications (Pearson r = 0.985). A difference was observed in ctDNA detection preoperatively (ST 72/92, MT 88/92). However, no difference was observed immediately after surgery in recurrence (ST 11/22, MT 10/22) or nonrecurrence (both 2/34) patients. In serial samples, detection was similar within recurrence (ST 13/16, MT 14/16) and nonrecurrence (ST 3/49, MT 1/49) patients. Both approaches yielded similar lead times to standard-of-care radiology (ST 4.0 months, MT 4.1 months). Our findings do not support significant performance gains of the MT strategy over the ST strategy for postoperative ctDNA detection.

Antisense transcription in gammaretroviruses as a mechanism of insertional activation of host genes.

Transcription of retroviruses is initiated at the U3-R region boundary in the integrated provirus and continues unidirectionally to produce genomic and mRNA products of positive polarity. Several studies have recently demonstrated the existence of naturally occurring protein-encoding transcripts of negative polarity in complex retroviruses. We report here on the identification of transcripts of negative polarity in simple murine leukemia virus (MLV). In T-cell and B-cell lymphomas induced by SL3-3 and Akv MLV, antisense transcripts initiated in the U3 region of the proviral 5' long terminal repeat (LTR) and continued into the cellular proto-oncogenes Jdp2 and Bach2 to create chimeric transcripts consisting of viral and host sequence. The phenomenon was validated in vivo using a knock-in mouse model homozygous for a single LTR at a position known to activate Nras in B-cell lymphomas. A 5' rapid amplification of cDNA ends (RACE) analysis indicated a broad spectrum of initiation sites within the U3 region of the 5' LTR. Our data show for the first time transcriptional activity of negative polarity initiating in the U3 region of simple retroviruses and suggest a novel mechanism of insertional activation of host genes. Elucidation of the nature and potential regulatory role of 5' LTR antisense transcription will be relevant to the design of therapeutic vectors and may contribute to the increasing recognition of pervasive eukaryotic transcription.

Activation of alternative Jdp2 promoters and functional protein isoforms in T-cell lymphomas by retroviral insertion mutagenesis.

Retroviral insertional mutagenesis has been instrumental for the identification of genes important in cancer development. The molecular mechanisms involved in retroviral-mediated activation of proto-oncogenes influence the distribution of insertions within specific regions during tumorigenesis and hence may point to novel gene structures. From a retroviral tagging screen on tumors of 1767 SL3-3 MLV-infected BALB/c mice, intron 2 of the AP-1 repressor Jdp2 locus was found frequently targeted by proviruses resulting in upregulation of non-canonical RNA subspecies. We identified several promoter regions within 1000 bp upstream of exon 3 that allowed for the production of Jdp2 protein isoforms lacking the histone acetylase inhibitory domain INHAT present in canonical Jdp2. The novel Jdp2 isoforms localized to the nucleus and over-expression in murine fibroblast cells induced cell death similar to canonic Jdp2. When expressed in the context of oncogenic NRAS both full length Jdp2 and the shorter isoforms increased anchorage-independent growth. Our results demonstrate a biological function of Jdp2 lacking the INHAT domain and suggest a post-genomic application for the use of retroviral tagging data in identifying new gene products with a potential role in tumorigenesis.

Correlation between early dynamics in circulating tumour DNA and outcome from FOLFIRI treatment in metastatic colorectal cancer.

Chemotherapy resistance remains a challenge in the clinical management of metastatic colorectal cancer (mCRC). Here, early changes in cell-free circulating tumour DNA (ctDNA) levels were explored as a marker of therapeutic efficacy. Twenty-four mCRC patients were enrolled and treated with FOLFIRI based first-line therapy. Blood samples collected pre-treatment, at day 7, 14, 21, 60 and at progression were analysed for cell-free DNA (cfDNA) and ctDNA levels using digital droplet PCR. A subset of samples were additionally analysed by targeted sequencing. Patients with high pre-treatment ctDNA or cfDNA levels (≥75th centile) had significantly shorter progression free survival (PFS) than patients with lower levels. Despite an overall decline in ctDNA levels from pre-treatment to first CT-scan, serial analysis identified seven patients with temporary increases in ctDNA consistent with growth of resistant cells. These patients had shorter PFS and shorter overall survival. Targeted sequencing analyses of cfDNA revealed dramatic changes in the clonal composition in response to treatment. Our study suggests that increasing ctDNA levels during the first cycles of first-line FOLFIRI treatment is a predictor of incipient progressive disease and poorer survival. Thus, we demonstrate the importance of monitoring ctDNA levels as early as one week after treatment onset to enable early detection of treatment failure.

Novel DNA methylation biomarkers show high sensitivity and specificity for blood-based detection of colorectal cancer-a clinical biomarker discovery and validation study.

BACKGROUND: Early detection plays an essential role to reduce colorectal cancer (CRC) mortality. While current screening methods suffer from poor compliance, liquid biopsy-based strategies for cancer detection is rapidly gaining promise. Here, we describe the development of TriMeth, a minimal-invasive blood-based test for detection of early-stage colorectal cancer. The test is based on assessment of three tumour-specific DNA methylation markers in circulating cell-free DNA. RESULTS: A thorough multi-step biomarker discovery study based on DNA methylation profiles of more than 5000 tumours and blood cell populations identified CRC-specific DNA methylation markers. The DNA methylation patterns of biomarker candidates were validated by bisulfite sequencing and methylation-specific droplet digital PCR in CRC tumour tissue and peripheral blood leucocytes. The three best performing markers were first applied to plasma from 113 primarily early-stage CRC patients and 87 age- and gender-matched colonoscopy-verified controls. Based on this, the test scoring algorithm was locked, and then TriMeth was validated in an independent cohort comprising 143 CRC patients and 91 controls. Three DNA methylation markers, C9orf50, KCNQ5, and CLIP4, were identified, each capable of discriminating plasma from colorectal cancer patients and healthy individuals (areas under the curve 0.86, 0.91, and 0.88). When combined in the TriMeth test, an average sensitivity of 85% (218/256) was observed (stage I: 80% (33/41), stage II: 85% (121/143), stage III: 89% (49/55), and stage IV: 88% (15/17)) at 99% (176/178) specificity in two independent plasma cohorts. CONCLUSION: TriMeth enables detection of early-stage colorectal cancer with high sensitivity and specificity. The reported results underline the potential utility of DNA methylation-based detection of circulating tumour DNA in the clinical management of colorectal cancer.

Molecular-Subtype-Specific Biomarkers Improve Prediction of Prognosis in Colorectal Cancer.

Colorectal cancer (CRC) is characterized by major inter-tumor diversity that complicates the prediction of disease and treatment outcomes. Recent efforts help resolve this by sub-classification of CRC into natural molecular subtypes; however, this strategy is not yet able to provide clinicians with improved tools for decision making. We here present an extended framework for CRC stratification that specifically aims to improve patient prognostication. Using transcriptional profiles from 1,100 CRCs, including >300 previously unpublished samples, we identify cancer cell and tumor archetypes and suggest the tumor microenvironment as a major prognostic determinant that can be influenced by the microbiome. Notably, our subtyping strategy allowed identification of archetype-specific prognostic biomarkers that provided information beyond and independent of UICC-TNM staging, MSI status, and consensus molecular subtyping. The results illustrate that our extended subtyping framework, combining subtyping and subtype-specific biomarkers, could contribute to improved patient prognostication and may form a strong basis for future studies.

miRNA profiling of circulating EpCAM(+) extracellular vesicles: promising biomarkers of colorectal cancer.

Cancer cells secrete small membranous extracellular vesicles (EVs) into their microenvironment and circulation. These contain biomolecules, including proteins and microRNAs (miRNAs). Both circulating EVs and miRNAs have received much attention as biomarker candidates for non-invasive diagnostics. Here we describe a sensitive analytical method for isolation and subsequent miRNA profiling of epithelial-derived EVs from blood samples of patients with colorectal cancer (CRC). The epithelial-derived EVs were isolated by immunoaffinity-capture using the epithelial cell adhesion molecule (EpCAM) as marker. This approach mitigates some of the specificity issues observed in earlier studies of circulating miRNAs, in particular the negative influence of miRNAs released by erythrocytes, platelets and non-epithelial cells. By applying this method to 2 small-scale patient cohorts, we showed that blood plasma isolated from CRC patients prior to surgery contained elevated levels of 13 EpCAM(+)-EV miRNAs compared with healthy individuals. Upon surgical tumour removal, the plasma levels of 8 of these were reduced (miR-16-5p, miR-23a-3p, miR-23b-3p, miR-27a-3p, miR-27b-3p, miR-30b-5p, miR-30c-5p and miR-222-3p). These findings indicate that the miRNAs are of tumour origin and may have potential as non-invasive biomarkers for detection of CRC. This work describes a non-invasive blood-based method for sensitive detection of cancer with potential for clinical use in relation to diagnosis and screening. We used the method to study CRC; however, it is not restricted to this disease. It may in principle be used to study any cancer that release epithelial-derived EVs into circulation.

miR-625-3p regulates oxaliplatin resistance by targeting MAP2K6-p38 signalling in human colorectal adenocarcinoma cells.

Oxaliplatin resistance in colorectal cancers (CRC) is a major medical problem, and predictive markers are urgently needed. Recently, miR-625-3p was reported as a promising predictive marker. Herein, we show that miR-625-3p functionally induces oxaliplatin resistance in CRC cells, and identify the signalling networks affected by miR-625-3p. We show that the p38 MAPK activator MAP2K6 is a direct target of miR-625-3p, and, accordingly, is downregulated in non-responder patients of oxaliplatin therapy. miR-625-3p-mediated resistance is reversed by anti-miR-625-3p treatment and ectopic expression of a miR-625-3p insensitive MAP2K6 variant. In addition, reduction of p38 signalling by using siRNAs, chemical inhibitors or expression of a dominant-negative MAP2K6 protein induces resistance to oxaliplatin. Transcriptome, proteome and phosphoproteome profiles confirm inactivation of MAP2K6-p38 signalling as one likely mechanism of oxaliplatin resistance. Our study shows that miR-625-3p induces oxaliplatin resistance by abrogating MAP2K6-p38-regulated apoptosis and cell cycle control networks, and corroborates the predictive power of miR-625-3p.

Activation of the brain-specific neurogranin gene in murine T-cell lymphomas by proviral insertional mutagenesis.

Neurogranin (Nrgn) is a highly expressed brain-specific protein, which sequesters calmodulin at low Ca(2+)-levels. We report here on retroviral activation of the Nrgn gene in tumors induced by the T-cell lymphomagenic SL3-3 murine leukemia virus. We have performed a systematic expression analysis of Nrgn in various mouse tissues and SL3-3 induced T-cell tumors. This demonstrated that insertional activation of Nrgn increased RNA and protein expression levels to that observed in brain. Furthermore, elevated Nrgn expression was also observed in some T-cell tumors with no detected provirus integrations into this genomic region. The presented data demonstrate that Nrgn can be produced at high levels outside the brain, and suggest a novel oncogenic role in T-cell lymphomas in mice.