Widespread adaptation of liquid biopsy for the early detection of cancer has yet to reach clinical utility. Circulating tumor DNA is commonly detected though the presence of genetic alterations, but only a minor fraction of tumor-derived cell-free DNA (cfDNA) fragments exhibit mutations. The cellular processes occurring in cancer development mark the chromatin. These epigenetic marks are reflected by modifications in the cfDNA methylation, fragment size, and structure. In this review, we describe how going beyond DNA sequence information alone, by analyzing cfDNA epigenetic and immune signatures, boosts the potential of liquid biopsy for the early detection of cancer.

Glioma‐derived cell‐free DNA (cfDNA) is challenging to detect using liquid biopsy because quantities in body fluids are low. We determined the glioma‐derived DNA fraction in cerebrospinal fluid (CSF), plasma, and urine samples from patients using sequencing of personalized capture panels guided by analysis of matched tumor biopsies. By sequencing cfDNA across thousands of mutations, identified individually in each patient’s tumor, we detected tumor‐derived DNA in the majority of CSF (7/8), plasma (10/12), and urine samples (10/16), with a median tumor fraction of 6.4 × 10−3, 3.1 × 10−5, and 4.7 × 10−5, respectively. We identified a shift in the size distribution of tumor‐derived cfDNA fragments in these body fluids. We further analyzed cfDNA fragment sizes using whole‐genome sequencing, in urine samples from 35 glioma patients, 27 individuals with non‐malignant brain disorders, and 26 healthy individuals. cfDNA in urine of glioma patients was significantly more fragmented compared to urine from patients with non‐malignant brain disorders (P = 1.7 × 10−2) and healthy individuals (P = 5.2 × 10−9). Machine learning models integrating fragment length could differentiate urine samples from glioma patients (AUC = 0.80–0.91) suggesting possibilities for truly non‐invasive cancer detection.

Background Assays that account for the biological properties and fragmentation of cell-free DNA (cfDNA) can improve the performance of liquid biopsy. However, preanalytic and physiological differences between individuals on fragmentomic analysis are poorly defined. Methods We analyzed the impact of collection tube, plasma processing time, and physiology on the size distribution of cfDNA, their genome-wide representation, and sequence diversity at the cfDNA fragment ends using shallow whole-genome sequencing. Results Neither different stabilizing collection tubes nor processing times affected the cfDNA fragment sizes, but could impact the genome-wide fragmentation patterns and fragment-end sequences of cfDNA. In addition, beyond differences depending on the gender, the physiological conditions tested between 63 individuals (age, body mass index, use of medication, and chronic conditions) minimally influenced the outcome of fragmentomic methods. Conclusions Fragmentomic approaches have potential for implementation in the clinic, pending clear traceability of analytical and physiological factors.

The structure, fragmentation pattern, length and terminal sequence of cell-free DNA (cfDNA) is under the influence of nucleases present in the blood. We hypothesized that differences in the diversity of bases at the end of cfDNA fragments can be leveraged on a genome-wide scale to enhance the sensitivity for detecting the presence of tumor signals in plasma. We surveyed the cfDNA termini in 572 plasma samples from 319 patients with 18 different cancer types using low-coverage whole genome sequencing. The fragment-end sequence and diversity were altered in all cancer types in comparison to 76 healthy controls. We converted the fragment end sequences into a quantitative metric and observed that this correlates with circulating tumor DNA tumor fraction (R = 0.58, p < 0.001, Spearman). Using these metrics, we were able to classify cancer samples from control at a low tumor content (AUROC of 91% at 1% tumor fraction) and shallow sequencing coverage (mean AUROC = 0.99 at >1M fragments). Combining fragment-end sequences and diversity using machine learning, we classified cancer from healthy controls (mean AUROC = 0.99, SD = 0.01). Using unsupervised clustering we showed that early-stage lung cancer (n = 13) can be classified from control or later stages based on fragment-end sequences. We observed that fragment-end sequences can be used for prognostication (hazard ratio: 0.49) and residual disease detection in resectable esophageal adenocarcinoma patients, moving fragmentomics toward a greater clinical implementation.

Cell-free DNA (cfDNA) can be isolated from blood and/or urine of cancer patients and analyzed with sequencing. Unfortunately, most conventional short-read sequencing methods are technically challenging, labor intensive and time consuming, requiring several days but more typically weeks to obtain interpretable data which are limited by a bias for short cfDNA fragments. Here, we demonstrate that with Oxford Nanopore Technologies sequencing we can achieve economical and ultra-fast delivery of clinical data from liquid biopsies. Our ITSFASTR approach is able to deliver copy number aberrations, and cfDNA fragmentation profiles in less than 24 hours from sample collection. The tumor-derived cfDNA fraction calculated from lung cancer patient plasma and urine from bladder cancer patients was highly correlated (R=0.98) to the tumor fraction calculated from conventional short-read sequencing of the same samples. cfDNA size profile and fragmentation patterns in plasma and urine exhibited the typical cfDNA features yet with a significantly higher proportion of fragments that exceed 300bp, exhibiting similar tumor fraction than shorter size fragments. Notably, comprehensive fragment-end composition and nucleosome profiling near transcription start sites can be determined from the same data. We propose that ITSFASTR is the first point-of-care solution for obtaining genomic and fragmentomic results from liquid biopsies.

423 Background: Despite the advent of precision medicine, prediction of survival outcome of esophageal cancer patients remains a challenge. Here we aim to investigate the value of prediction models integrating multi-signal data including radiomics and circulating tumor DNA (ctDNA) data in addition to clinical data for the prediction of resectable esophageal adenocarcinoma (rEAC) related outcomes. Methods: In total n=111 rEAC patients treated with neoadjuvant chemoradiotherapy (nCRT; n=71) +/- anti-PD-L1 (n=40) were included. Baseline clinical variables (n=9) were based on the SOURCE survival prediction model (van den Boorn et al. JNCCN. 2021). The baseline ctDNA data from plasma was derived from fragmentomic and copy number aberrations (ichorCNA) from shallow whole genome sequencing (<5-fold coverage) and a custom next-generation sequencing panel (n=23 genes). Baseline radiomic original features were extracted by the pyradiomics package from CT-image delineated tumor volumes. An initial redundancy filtering was performed to remove correlating variables (r>0.6). We evaluated the predictive performance of baseline ctDNA and radiomics features on overall survival (OS), progression free survival (PFS), and time to progression (TTP), through fitting Cox-regression models. Four ctDNA features were included in the models: P20-150, ichorCNA, fragment end score and mutation detection. For the radiomics features we performed an additional back- and forward variable selection based on Akaike’s Information Criterion. Using the likelihood ratio test we tested if the model fit changed after adding ctDNA and radiomics features to a model with clinical variables. Results: The addition of radiomics to clinical variables improved model fit for OS (p=0.017). Baseline prediction of OS resulted in a C-index of 0.65 using clinical variables only, 0.65 with ctDNA, 0.68 with radiomics and 0.68 with ctDNA and radiomics combined. For PFS model fit improved after adding radiomics (p=0.020) and ctDNA and radiomics combined (p=0.017). Baseline prediction of PFS resulted in a C-index of 0.64 using clinical variables, 0.65 with ctDNA, 0.67 with radiomics, and 0.68 with ctDNA and radiomics combined. For TTP model fit improved after adding radiomics (p=0.008) and radiomics and ctDNA combined (p=0.002). Baseline prediction of TTP resulted in a C-index of 0.64 with clinical variables, 0.65 with ctDNA, 0.71 with radiomics, and 0.72 with ctDNA and radiomics combined. Based on the cox-regression models using clinical variables and radiomics, risk stratification by splitting the cohort in a high and low risk group was possible for OS, PFS and TTP (p<0.001). Conclusions: Combining clinical variables from SOURCE with radiomics data improved predictions of OS, PFS, and TTP among patients with rEAC. Multi-signal integration of clinical and radiomics variables could potentially be used to identify risk groups.

Assays that account for the biological properties and fragmentation of cell-free DNA (cfDNA) can improve the performance of liquid biopsy. However, pre-analytic and physiological differences between individuals on fragmentomic analysis are poorly defined. We analyzed the impact of collection tube, plasma processing time and physiology on the size distribution of cfDNA, their genome-wide representation and sequence diversity at the cfDNA fragment-ends using shallow Whole Genome Sequencing. We observed that using different stabilizing collection tubes, or processing times does not affect the cfDNA fragment sizes, but can impact the genome-wide fragmentation patterns and fragment-end sequences of cfDNA. In addition, beyond differences depending on the gender, the physiological conditions tested between 63 individuals (age, body mass index, use of medication and chronic conditions) minimally influenced the outcome of fragmentomic methods. Our results highlight that fragmentomic approaches have potential for implementation in the clinic, pending clear traceability of analytical and physiological factors.

Liquid biopsies contain multiple analytes that can be mined to improve the detection and management of cancer. Beyond cell-free DNA (cfDNA), mutations have been detected in DNA associated with extracellular vesicles (EV-DNA). The genome-wide composition and structure of EV-DNA are poorly characterized, and it remains undecided whether circulating EVs are enriched in tumor signal compared to unfractionated cfDNA. Here, using whole genome sequencing from selected lung cancer patients with a high cfDNA tumor content (>5%), we determined that the tumor fraction and heterogeneity are comparable between DNA associated with EVs and matched plasma cfDNA. DNA in EV fractions, obtained with standardized size-exclusion chromatography, are comprised of short ~150-180 bp fragments and long >1000 bp fragments that are poor in tumor signal. Other fractions only exhibit short fragments with similar tumor DNA content. The composition in bases at the end of EV-DNA fragments, as well as their fragmentation patterns are similar to plasma cfDNA. Mitochondrial DNA is relatively enriched in EV fractions. Our results highlight that cfDNA in plasma is of dual nature, either bound to proteins (including the nucleosome) but also associated to EV. cfDNA associated to small EV (including exosomes) is however not preferentially enriched in tumor signal.