FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA.

Analysis of DNA methylation in cell-free DNA reveals clinically relevant biomarkers but requires specialized protocols such as whole-genome bisulfite sequencing. Meanwhile, millions of cell-free DNA samples are being profiled by whole-genome sequencing. Here, we develop FinaleMe, a non-homogeneous Hidden Markov Model, to predict DNA methylation of cell-free DNA and, therefore, tissues-of-origin, directly from plasma whole-genome sequencing. We validate the performance with 80 pairs of deep and shallow-coverage whole-genome sequencing and whole-genome bisulfite sequencing data.

FinaleMe: Predicting DNA methylation by the fragmentation patterns of plasma cell-free DNA.

Analysis of DNA methylation in cell-free DNA (cfDNA) reveals clinically relevant biomarkers but requires specialized protocols and sufficient input material that limits its applicability. Millions of cfDNA samples have been profiled by genomic sequencing. To maximize the gene regulation information from the existing dataset, we developed FinaleMe, a non-homogeneous Hidden Markov Model (HMM), to predict DNA methylation of cfDNA and, therefore, tissues-of-origin directly from plasma whole-genome sequencing (WGS). We validated the performance with 80 pairs of deep and shallow-coverage WGS and whole-genome bisulfite sequencing (WGBS) data.

Single duplex DNA sequencing with CODEC detects mutations with high sensitivity.

Detecting mutations from single DNA molecules is crucial in many fields but challenging. Next-generation sequencing (NGS) affords tremendous throughput but cannot directly sequence double-stranded DNA molecules ('single duplexes') to discern the true mutations on both strands. Here we present Concatenating Original Duplex for Error Correction (CODEC), which confers single duplex resolution to NGS. CODEC affords 1,000-fold higher accuracy than NGS, using up to 100-fold fewer reads than duplex sequencing. CODEC revealed mutation frequencies of 2.72 × 10-8 in sperm of a 39-year-old individual, and somatic mutations acquired with age in blood cells. CODEC detected genome-wide, clonal hematopoiesis mutations from single DNA molecules, single mutated duplexes from tumor genomes and liquid biopsies, microsatellite instability with 10-fold greater sensitivity and mutational signatures, and specific tumor mutations with up to 100-fold fewer reads. CODEC enables more precise genetic testing and reveals biologically significant mutations, which are commonly obscured by NGS errors.

Massively parallel enrichment of low-frequency alleles enables duplex sequencing at low depth.

Assaying for large numbers of low-frequency mutations requires sequencing at extremely high depth and accuracy. Increasing sequencing depth aids the detection of low-frequency mutations yet limits the number of loci that can be simultaneously probed. Here we report a method for the accurate tracking of thousands of distinct mutations that requires substantially fewer reads per locus than conventional hybrid-capture duplex sequencing. The method, which we named MAESTRO (for minor-allele-enriched sequencing through recognition oligonucleotides), combines massively parallel mutation enrichment with duplex sequencing to track up to 10,000 low-frequency mutations, with up to 100-fold fewer reads per locus. We show that MAESTRO can be used to test for chimaerism by tracking donor-exclusive single-nucleotide polymorphisms in sheared genomic DNA from human cell lines, to validate whole-exome sequencing and whole-genome sequencing for the detection of mutations in breast-tumour samples from 16 patients, and to monitor the patients for minimal residual disease via the analysis of cell-free DNA from liquid biopsies. MAESTRO improves the breadth, depth, accuracy and efficiency of mutation testing by sequencing.

Liquid biopsy detection of genomic alterations in pediatric brain tumors from cell-free DNA in peripheral blood, CSF, and urine.

BACKGROUND: The ability to identify genetic alterations in cancers is essential for precision medicine; however, surgical approaches to obtain brain tumor tissue are invasive. Profiling circulating tumor DNA (ctDNA) in liquid biopsies has emerged as a promising approach to avoid invasive procedures. Here, we systematically evaluated the feasibility of profiling pediatric brain tumors using ctDNA obtained from plasma, cerebrospinal fluid (CSF), and urine. METHODS: We prospectively collected 564 specimens (257 blood, 240 urine, and 67 CSF samples) from 258 patients across all histopathologies. We performed ultra-low-pass whole-genome sequencing (ULP-WGS) to assess copy number variations and estimate tumor fraction and developed a pediatric CNS tumor hybrid capture panel for deep sequencing of specific mutations and fusions. RESULTS: ULP-WGS detected copy number alterations in 9/46 (20%) CSF, 3/230 (1.3%) plasma, and 0/153 urine samples. Sequencing detected alterations in 3/10 (30%) CSF, 2/74 (2.7%) plasma, and 0/2 urine samples. The only positive results were in high-grade tumors. However, most samples had insufficient somatic mutations (median 1, range 0-39) discoverable by the sequencing panel to provide sufficient power to detect tumor fractions of greater than 0.1%. CONCLUSIONS: Children with brain tumors harbor very low levels of ctDNA in blood, CSF, and urine, with CSF having the most DNA detectable. Molecular profiling is feasible in a small subset of high-grade tumors. The level of clonal aberrations per genome is low in most of the tumors, posing a challenge for detection using whole-genome or even targeted sequencing methods. Substantial challenges therefore remain to genetically characterize pediatric brain tumors from liquid biopsies.

Detection of Leptomeningeal Disease Using Cell-Free DNA From Cerebrospinal Fluid.

Importance: Leptomeningeal disease (LMD) is a devastating complication of cancer that is frequently underdiagnosed owing to the low sensitivity of cerebrospinal fluid (CSF) cytologic assessment, the current benchmark diagnostic method. Improving diagnostic sensitivity may lead to improved treatment decisions. Objective: To assess whether cell-free DNA (cfDNA) analysis of CSF may be used to diagnose LMD more accurately than cytologic analysis. Design, Setting, and Participants: This diagnostic study conducted in a neuro-oncology clinic at 2 large, tertiary medical centers assessed the use of genomic sequencing of CSF samples obtained from 30 patients with suspected or confirmed LMD from 2015 through 2018 to identify tumor-derived cfDNA. From the same CSF samples, cytologic analyses were conducted, and the results of the 2 tests were compared. This study consisted of 2 patient populations: 22 patients with cytologically confirmed LMD without parenchymal tumors abutting their CSF and 8 patients with parenchymal brain metastases with no evidence of LMD. Patients were considered positive for the presence of LMD if previous CSF cytologic analysis was positive for malignant cells. The analysis was conducted from 2015 to 2018. Main Outcomes and Measures: The primary outcome was the diagnostic accuracy of cfDNA analysis, defined as the number of tests that resulted in correct diagnoses out of the total number of tests assayed. Hypotheses were formed before data collection. Results: In total, 30 patients (23 women [77%]; median age, 51 years [range, 28-81 years]), primarily presenting with metastatic solid malignant neoplasms, participated in this study. For 48 follow-up samples from patients previously diagnosed via cytologic analysis as having LMD with no parenchymal tumor abutting CSF, cfDNA findings were accurate in the assessment of LMD in 45 samples (94%; 95% CI, 83%-99%), whereas cytologic analysis was accurate in 36 samples (75%; 95% CI, 60%-86%), a significant difference (P = .02). Of 43 LMD-positive samples, CSF cfDNA analysis was sensitive to LMD in 40 samples (93%; 95% CI, 81%-99%), and cytologic analysis was sensitive to LMD in 31 samples (72%; 95% CI, 56%-85%), a significant difference (P = .02). For 3 patients with parenchymal brain metastases abutting the CSF and no suspicion of LMD, cytologic findings were negative for LMD in all 3 patients, whereas cfDNA findings were positive in all 3 patients. Conclusions and Relevance: This diagnostic study found improved sensitivity and accuracy of cfDNA CSF testing vs cytologic assessment for diagnosing LMD with the exception of parenchymal tumors abutting CSF, suggesting improved ability to diagnosis LMD. Consideration of incorporating CSF cfDNA analysis into clinical care is warranted.

Circulating Cell-Free DNA as Biomarker of Taxane Resistance in Metastatic Castration-Resistant Prostate Cancer.

There are no biomarkers predictive of resistance to docetaxel or cabazitaxel validated for patients with metastatic castration-resistant prostate cancer (mCRPC). We assessed the association between ABCB1 amplification and primary resistance to docetaxel or cabazitaxel for patients with mCRPC, using circulating cell-free DNA (cfDNA). Patients with ≥1 plasma sample drawn within 12 months before starting docetaxel (cohort A) or cabazitaxel (cohort B) for mCRPC were identified from the Dana-Farber Cancer Institute IRB approved database. Sparse whole genome sequencing was performed on the selected cfDNA samples and tumor fractions were estimated using the computational tool ichorCNA. We evaluated the association between ABCB1 amplification or other copy number alterations and primary resistance to docetaxel or cabazitaxel. Of the selected 176 patients, 45 samples in cohort A and 21 samples in cohort B had sufficient tumor content. No significant association was found between ABCB1 amplification and primary resistance to docetaxel (p = 0.58; odds ratio (OR) = 1.49) or cabazitaxel (p = 0.97; OR = 1.06). No significant association was found between exploratory biomarkers and primary resistance to docetaxel or cabazitaxel. In this study, ABCB1 amplification did not predict primary resistance to docetaxel or cabazitaxel for mCRPC. Future studies including ABCB1 amplification in a suite of putative biomarkers and a larger cohort may aid in drawing definitive conclusions.

Sensitive Detection of Minimal Residual Disease in Patients Treated for Early-Stage Breast Cancer.

PURPOSE: Existing cell-free DNA (cfDNA) methods lack the sensitivity needed for detecting minimal residual disease (MRD) following therapy. We developed a test for tracking hundreds of patient-specific mutations to detect MRD with a 1,000-fold lower error rate than conventional sequencing. EXPERIMENTAL DESIGN: We compared the sensitivity of our approach to digital droplet PCR (ddPCR) in a dilution series, then retrospectively identified two cohorts of patients who had undergone prospective plasma sampling and clinical data collection: 16 patients with ER+/HER2- metastatic breast cancer (MBC) sampled within 6 months following metastatic diagnosis and 142 patients with stage 0 to III breast cancer who received curative-intent treatment with most sampled at surgery and 1 year postoperative. We performed whole-exome sequencing of tumors and designed individualized MRD tests, which we applied to serial cfDNA samples. RESULTS: Our approach was 100-fold more sensitive than ddPCR when tracking 488 mutations, but most patients had fewer identifiable tumor mutations to track in cfDNA (median = 57; range = 2-346). Clinical sensitivity was 81% (n = 13/16) in newly diagnosed MBC, 23% (n = 7/30) at postoperative and 19% (n = 6/32) at 1 year in early-stage disease, and highest in patients with the most tumor mutations available to track. MRD detection at 1 year was strongly associated with distant recurrence [HR = 20.8; 95% confidence interval, 7.3-58.9]. Median lead time from first positive sample to recurrence was 18.9 months (range = 3.4-39.2 months). CONCLUSIONS: Tracking large numbers of individualized tumor mutations in cfDNA can improve MRD detection, but its sensitivity is driven by the number of tumor mutations available to track.

A nuclease-polymerase chain reaction enables amplification of probes used for capture-based DNA target enrichment.

DNA target enrichment via hybridization capture is a commonly adopted approach which remains expensive due in-part to using biotinylated-probe panels. Here we provide a novel isothermal amplification reaction to amplify rapidly existing probe panels without knowledge of the sequences involved, thereby decreasing a major portion of the overall sample preparation cost. The reaction employs two thermostable enzymes, BST-polymerase and duplex-specific nuclease DSN. DSN initiates random 'nicks' on double-stranded-DNA which enable BST to polymerize DNA by displacing the nicked-strand. Displaced strands re-hybridize and the process leads to an exponential chain-reaction generating biotinylated DNA fragments within minutes. When starting from single-stranded-DNA, DNA is first converted to double-stranded-DNA via terminal-deoxynucleotidyl-transferase (TdT) prior to initiation of BST-DSN reaction. Biotinylated probes generated by TdT-BST-DSN (TBD) reactions using panels of 33, 190 or 7186 DNA targets are used for hybrid-capture-based target enrichment from amplified circulating-DNA, followed by targeted re-sequencing. Polymerase-nuclease isothermal-chain-reactions generate random amplified probes with no apparent sequence dependence. One round of target-capture using TBD probes generates a modest on-target sequencing ratio, while two successive rounds of capture generate >80% on-target reads with good sequencing uniformity. TBD-reactions generate enough capture-probes to increase by approximately two to three orders-of-magnitude the target-enrichment experiments possible from an initial set of probes.

Tumor fraction in cell-free DNA as a biomarker in prostate cancer.

BACKGROUND: Tumor content in circulating cell-free DNA (cfDNA) is a promising biomarker, but longitudinal dynamics of tumor-derived and non-tumor-derived cfDNA through multiple courses of therapy have not been well described. METHODS: CfDNA from 663 plasma samples from 140 patients with castration-resistant prostate cancer (CRPC) was subject to sparse whole genome sequencing. Tumor fraction (TFx) estimated using the computational tool ichorCNA was correlated with clinical features and responses to therapy. RESULTS: TFx associated with the number of bone metastases (median TFx = 0.014 with no bone metastases, 0.047 with 1-3 bone metastases, 0.190 for 4+ bone metastases; P < 0.0001) and with visceral metastases (P < 0.0001). In multivariable analysis, TFx remained associated with metastasis location (P = 0.042); TFx was positively correlated with alkaline phosphatase (P = 0.0227) and negatively correlated with hemoglobin (Hgb) (P < 0.001), but it was not correlated with prostate specific antigen (PSA) (P = 0.75). Tumor-derived and non-tumor-derived cfDNA track together and do not increase with generalized tissue damage from chemotherapy or radiation at the time scales examined. All new treatments that led to ≥30% PSA decline at 6 weeks were associated with TFx decline when baseline TFx was >7%; however, TFx in patients being subsequently maintained on secondary hormonal therapy was quite dynamic. CONCLUSION: TFx correlates with clinical features associated with overall survival in CRPC, and TFx decline is a promising biomarker for initial therapeutic response. TRIAL REGISTRATION: Dana-Farber/Harvard Cancer Center (DF/HCC) protocol no. 18-135. FUNDING: Wong Family Award in Translational Oncology, Dana Farber Cancer Institute Medical Oncology grant, Gerstner Family Foundation, Janssen Pharmaceuticals Inc., and Koch Institute Support (core) grant P30-CA14051 from the National Cancer Institute (NCI).

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

Association of Cell-Free DNA Tumor Fraction and Somatic Copy Number Alterations With Survival in Metastatic Triple-Negative Breast Cancer.

Purpose Cell-free DNA (cfDNA) offers the potential for minimally invasive genome-wide profiling of tumor alterations without tumor biopsy and may be associated with patient prognosis. Triple-negative breast cancer (TNBC) is characterized by few mutations but extensive somatic copy number alterations (SCNAs), yet little is known regarding SCNAs in metastatic TNBC. We sought to evaluate SCNAs in metastatic TNBC exclusively via cfDNA and determine if cfDNA tumor fraction is associated with overall survival in metastatic TNBC. Patients and Methods In this retrospective cohort study, we identified 164 patients with biopsy-proven metastatic TNBC at a single tertiary care institution who received prior chemotherapy in the (neo)adjuvant or metastatic setting. We performed low-coverage genome-wide sequencing of cfDNA from plasma. Results Without prior knowledge of tumor mutations, we determined tumor fraction of cfDNA for 96.3% of patients and SCNAs for 63.9% of patients. Copy number profiles and percent genome altered were remarkably similar between metastatic and primary TNBCs. Certain SCNAs were more frequent in metastatic TNBCs relative to paired primary tumors and primary TNBCs in publicly available data sets The Cancer Genome Atlas and METABRIC, including chromosomal gains in drivers NOTCH2, AKT2, and AKT3. Prespecified cfDNA tumor fraction threshold of ≥ 10% was associated with significantly worse metastatic survival (median, 6.4 v 15.9 months) and remained significant independent of clinicopathologic factors (hazard ratio, 2.14; 95% CI, 1.4 to 3.8; P < .001). Conclusion We present the largest genomic characterization of metastatic TNBC to our knowledge, exclusively from cfDNA. Evaluation of cfDNA tumor fraction was feasible for nearly all patients, and tumor fraction ≥ 10% is associated with significantly worse survival in this large metastatic TNBC cohort. Specific SCNAs are enriched and prognostic in metastatic TNBC, with implications for metastasis, resistance, and novel therapeutic approaches.

Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors.

Whole-exome sequencing of cell-free DNA (cfDNA) could enable comprehensive profiling of tumors from blood but the genome-wide concordance between cfDNA and tumor biopsies is uncertain. Here we report ichorCNA, software that quantifies tumor content in cfDNA from 0.1× coverage whole-genome sequencing data without prior knowledge of tumor mutations. We apply ichorCNA to 1439 blood samples from 520 patients with metastatic prostate or breast cancers. In the earliest tested sample for each patient, 34% of patients have ≥10% tumor-derived cfDNA, sufficient for standard coverage whole-exome sequencing. Using whole-exome sequencing, we validate the concordance of clonal somatic mutations (88%), copy number alterations (80%), mutational signatures, and neoantigens between cfDNA and matched tumor biopsies from 41 patients with ≥10% cfDNA tumor content. In summary, we provide methods to identify patients eligible for comprehensive cfDNA profiling, revealing its applicability to many patients, and demonstrate high concordance of cfDNA and metastatic tumor whole-exome sequencing.

Regulation of fibrin-mediated tumor cell adhesion to the endothelium using anti-thrombin aptamer.

Molecular intervention during transient stages of various metastatic pathways may lead to development of promising therapeutic technologies. One of such involves soluble fibrin (sFn) that has been implicated as a cross-linker between circulating blood or tumor cells and endothelial cell receptors, promoting cell arrest on the endothelium during circulation. sFn generation is a result of thrombin-mediated fibrinogen (Fg) cleavage due to either vascular injuries or a tumor microenvironment. For cancer therapy, thrombin-mediated conversions of Fg to sFn thus serve as potential intervention points to decrease circulating tumor cell adhesion to the endothelium and subsequent metastatic events. The purpose of this work was to investigate the function of an anti-thrombin oligonucleotide aptamer in reducing tumor cell arrest. Both molecular and cellular interactions were examined to demonstrate the binding and inhibitory effects of anti-thrombin aptamer. The results show that the aptamer is capable of inhibiting thrombin-mediated Fg conversion, thereby reducing sFn-mediated tumor cell adhesion in a concentration-dependent manner. Notably, the aptamer is able to bind thrombin under dynamic flow conditions and reduce tumor cell adhesive events at various physiological shear rates. This study further indicates that oligonucleotide aptamers hold great promise as therapeutic regulators of tumor cell adhesion, and consequently, metastatic activity.

Aptamer-based polyvalent ligands for regulated cell attachment on the hydrogel surface.

Natural biomolecules are often used to functionalize materials to achieve desired cell-material interactions. However, their applications can be limited owing to denaturation during the material functionalization process. Therefore, efforts have been made to develop synthetic ligands with polyvalence as alternatives to natural affinity biomolecules for the synthesis of functional materials and the control of cell-material interactions. This work was aimed at investigating the capability of a hydrogel functionalized with a novel polyvalent aptamer in inducing cell attachment in dynamic flow and releasing the attached cells in physiological conditions through a hybridization reaction. The results show that the polyvalent aptamer could induce cell attachment on the hydrogel in dynamic flow. Moreover, cell attachment on the hydrogel surface was significantly influenced by the value of shear stress. The cell density on the hydrogel was increased from 40 cells/mm(2) to nearly 700 cells/mm(2) when the shear stress was decreased from 0.05 to 0.005 Pa. After the attachment onto the hydrogel surface, approximately 95% of the cells could be triggered to detach within 20 min by using an oligonucleotide complementary sequence that displaced polyvalent aptamer strands from the hydrogel surface. While it was found that the cell activity was reduced, the live/dead staining results show that ≥98% of the detached cells were viable. Therefore, this work has suggested that the polyvalent aptamer is a promising synthetic ligand for the functionalization of materials for regulated cell attachment.