Search

Duplex sequencing (DuplexSeq) is an error-corrected next-generation sequencing (ecNGS) method in which molecular barcodes informatically link PCR-copies back to their source DNA strands, enabling computational removal of errors by comparing grouped strand sequencing reads. The resulting background of less than one artifactual mutation per 107nucleotides allows for direct detection of somatic mutations. TwinStrand Biosciences, Inc. has developed a DuplexSeq-based mutagenesis assay to sample the rat genome, which can be applied to genetic toxicity testing. To evaluate this assay for early detection of mutagenesis, a time-course study was conducted using male Hsd:Sprague Dawley SD rats (3 per group) administered a single dose of 40 mg/kg N-ethyl-N-nitrosourea (ENU) via gavage, with mutation frequency (MF) and spectrum analyzed in stomach, bone marrow, blood, and liver tissues at 3 h, 24 h, 7 d, and 28 d post-exposure. Significant increases in MF were observed in ENU-exposed rats as early as 24 h for stomach (site of contact) and bone marrow (a highly proliferative tissue) and at 7 d for liver and blood. The canonical, mutational signature of ENU was established by 7 d post-exposure in all four tissues. Interlaboratory analysis of a subset of samples from different tissues and time points demonstrated remarkable reproducibility for both MF and spectrum. These results demonstrate that MF and spectrum can be evaluated successfully by directly sequencing targeted regions of DNA obtained from various tissues, a considerable advancement compared to currently usedin vivogene mutation assays.HIGHLIGHTSDuplexSeq is an ultra-accurate NGS technology that directly quantifies mutationsENU-dependent mutagenesis was detected 24 h post-exposure in proliferative tissuesMultiple tissues exhibited the canonical ENU mutation spectrum 7 d after exposureResults obtained with DuplexSeq were highly concordant between laboratoriesThe Rat-50 Mutagenesis Assay is promising for applications in genetic toxicology

Mutagenicity testing is an essential component of health safety assessment. Duplex Sequencing (DS), an emerging high-accuracy DNA sequencing technology, may provide substantial advantages over conventional mutagenicity assays. DS could be used to eliminate reliance on standalone reporter assays and provide mechanistic information alongside mutation frequency (MF) data. However, the performance of DS must be thoroughly assessed before it can be routinely implemented for standard testing. We used DS to study spontaneous and procarbazine (PRC)-induced mutations in the bone marrow (BM) of MutaMouse males across a panel of 20 diverse genomic targets. Mice were exposed to 0, 6.25, 12.5, or 25 mg/kg-bw/day for 28 days by oral gavage and BM sampled 42 days post-exposure. Results were compared with those obtained using the conventionallacZviral plaque assay on the same samples. DS detected significant increases in mutation frequencies and changes to mutation spectra at all PRC doses. Low intra-group variability within DS samples allowed for detection of increases at lower doses than thelacZassay. While thelacZassay initially yielded a higher fold-change in mutant frequency than DS, inclusion of clonal mutations in DS mutation frequencies reduced this discrepancy. Power analyses suggested that three animals per dose group and 500 million duplex base pairs per sample is sufficient to detect a 1.5-fold increase in mutations with >80% power. Overall, we demonstrate several advantages of DS over classical mutagenicity assays and provide data to support efforts to identify optimal study designs for the application of DS as a regulatory test.

Error-corrected Duplex Sequencing (DuplexSeq) enables direct quantification of low-frequency mutations and offers tremendous potential for chemical mutagenicity assessment. We investigated the utility of DuplexSeq to quantify induced mutation frequency (MF) and spectrum in human lymphoblastoid TK6 cells exposed to a prototypical DNA alkylating agent,N-ethyl-N-nitrosourea (ENU). Furthermore, we explored appropriate experimental parameters for this application, and assessed inter-laboratory reproducibility. In two independent experiments in two laboratories, TK6 cells were exposed to ENU (25-200 µM) and DNA was sequenced 48, 72, and 96 h post-exposure. A DuplexSeq mutagenicity panel targeting twenty 2.4-kb regions distributed across the genome was used to sample diverse, genome-representative sequence contexts. A robust increase in MF that was unaffected by time was observed in both laboratories. Concentration-response in the MF from the two laboratories was strongly positively correlated (R2=0.95). C:G>T:A, T:A>C:G, T:A>A:T, and T:A>G:C mutations increased in consistent, concentration-dependent manners in both laboratories, with high proportions of C:G>T:A at all time points. The target sites responded similarly between the two laboratories and revealed a higher average MF in intergenic regions. These results, demonstrating remarkable reproducibility across time and laboratory for both MF and spectrum, support the high value of DuplexSeq for characterizing chemical mutagenicity in both research and regulatory evaluation.

The organotypic human air‐liquid‐interface (ALI) airway tissue model has been used as an in vitro cell culture system for evaluating the toxicity of inhaled substances. ALI airway cultures are highly differentiated, which has made it challenging to evaluate genetic toxicology endpoints. In the current study, we assayed DNA damage with the high‐throughput CometChip assay and quantified mutagenesis with Duplex Sequencing, an error‐corrected next‐generation sequencing method capable of detecting a single mutation per 107 base pairs. Fully differentiated human ALI airway cultures were treated from the basolateral side with 6.25 to 100 μg/mL ethyl methanesulfonate (EMS) over a period of 28 days. CometChip assays were conducted after 3 and 28 days of treatment, and Duplex Sequencing after 28 days of treatment. Treating the airway cultures with EMS resulted in time‐ and concentration‐dependent increases in DNA damage and a concentration‐dependent increase in mutant frequency. The mutations observed in the EMS‐treated cultures were predominantly C → T transitions and exhibited a unique trinucleotide signature relative to the negative control. Measurement of physiological endpoints indicated that the EMS treatments had no effect on anti‐p63‐positive basal cell frequency, but produced concentration‐responsive increases in cytotoxicity and perturbations in cell morphology, along with concentration‐responsive decreases in culture viability, goblet cell and anti‐Ki67‐positive proliferating cell frequency, cilia beating frequency, and mucin secretion. The results indicate that a unified 28‐day study can be used to measure several important safety endpoints in physiologically relevant human in vitro ALI airway cultures, including DNA damage, mutagenicity, and tissue‐specific general toxicity.

Error-corrected sequencing of genomic targets enriched by probe-based capture has become a standard approach for detecting single-nucleotide variants (SNVs) and small insertion/deletions (indels) present at very low variant allele frequencies. Less attention has been given to strategies for comparable detection of rare structural variant (SV) junctions, where different error mechanisms must be addressed. Working from cell samples with known SV properties, we demonstrate that Duplex Sequencing (DuplexSeq), which demands confirmation of variants on both strands of a source DNA molecule, eliminates false SV junctions arising from chimeric PCR. DuplexSeq could not address frequent intermolecular ligation artifacts that arise during Y-adapter addition prior to strand denaturation without requiring multiple source molecules. In contrast, tagmentation libraries coupled with data filtering based on strand family size greatly reduced both artifact classes and enabled efficient and specific detection of even single-molecule SV junctions. The throughput of SV capture sequencing (svCapture) and the high base-level accuracy of DuplexSeq provided detailed views of the microhomology profile and limited occurrence of de novo SNVs near the junctions of hundreds of sub-clonal and newly created SVs, suggesting end joining as a predominant formation mechanism. The open source svCapture pipeline enables rare SV detection as a routine addition to SNVs/indels in properly prepared capture sequencing libraries.

Sensitive and specific detection of low frequency mutations in acute myeloid leukemia (AML) is critical in research of minimal residual disease (MRD). At variant allele frequencies (VAF) below approximately 1%, PCR and sequencing errors result in prohibitive signal-to-noise ratios with next generation sequencing (NGS). Duplex Sequencing (DS) relates the original top and bottom DNA strands to make double stranded consensus sequences to greatly reduce errors. In two experiments we show further improvement of DS of AML-related genes by the use of enzymatic fragmentation (EF) and an updated gene panel that incorporates 2022 European LeukemiaNet (ELN) recommendations. In the first DS study, 29 AML-related genes (59 kb) were targeted in hybrid capture. Mutant cell line DNA was mixed into DNA from a healthy young donor to simulate MRD. Expected VAFs in mixtures were 1.0-0.003%. Samples comprised a mix of 4 insertions and deletions (indel), a mix of 15 single nucleotide variants (SNV), and 4 serial dilutions of a FLT3 ITD plus an NPM1 insertion. Pure diluent DNA was used as a negative control. Each sample was prepared in quadruplicate with mechanical fragmentation (MF) vs EF. DNA input mass was 1,500 ng per replicate, except for 50-250 ng for the FLT3/NPM1 mixes with expected VAF of 1% and 0.1%, respectively. Input masses were set to ensure >95% probability of detection of all mutations in the combined data. All targeted variants down to 0.003% VAF were detected in spike-in mixtures, with expected vs observed VAF highly correlated whether using MF or EF (R2 >0.98 for all mixtures). Duplex molecular depth was 1.2-2.0x higher with EF vs MF across input masses. Panel-wide mean duplex depth per 1,500 ng replicate was 30,376x for MF and 48,036x for EF, for combined mean depths of 121,508x for MF and 192,144x for EF. In the negative control, background mutational calls at spike-in positions were 4/2,993,429 duplex bases using MF and 0/4,780,491 duplex bases with EF, reflecting increased specificity with EF. Next, a revised panel targeting 36 AML-related genes (80 kb) was used for DS with EF. Here differing mutant cell line DNA was mixed into the same diluent from above (the negative control). This mixture harbors 27 variants with expected VAFs of 0.125-0.006% (4 indels, 21 SNVs, a FLT3 ITD and an NPM1 insertion), 17 of which (2 indels, 14 SNVs and the NPM1 insertion) span VAFs of 0.011-0.009% to establish a limit of detection of 0.01%. With 2,000 ng (+/- 10%) input, for 14 samples, panel-wide mean duplex depth was 56,528x. For variants with expected VAF ≥ 0.009%, the assay had >98% sensitivity, >96% specificity, and >98% accuracy at genomic positions of known true positive variants. In summary, DS with EF yields more data per nanogram of DNA than MF, maximizing the use of precious samples. This method exhibits extremely low background mutation signal and a low limit of detection across targets that are valuable in AML MRD research. Citation Format: Elizabeth Schmidt, Devon M. Fitzgerald, Camila Zanette, Gavin D. Meredith, Mark A. McElwain, Raul Burciaga, Kevin C. Vavra, Thomas H. Smith, Jesse J. Salk, Jake Higgins. Accurate detection of low frequency AML-associated mutations in vitro using Duplex Sequencing with enzymatic fragmentation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 2 (Clinical Trials and Late-Breaking Research); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(8_Suppl):Abstract nr LB254.

Mutations have a profound effect on human health, particularly through an increased risk of carcinogenesis and genetic disease. The strong correlation between mutagenesis and carcinogenesis has been a driving force behind genotoxicity research for more than 50 years. The stochastic and infrequent nature of mutagenesis makes it challenging to observe and to study. Indeed, decades have been spent developing increasingly sophisticated assays and methods to study these low-frequency genetic errors, in hopes of better predicting which chemicals may be carcinogens, understanding their mode of action, and informing guidelines to prevent undue human exposure. While effective, widely used genetic selection-based technologies have a number of limitations that have hampered major advancements in the field of genotoxicity. Emerging new tools, in the form of enhanced next-generation sequencing platforms and methods, are changing this paradigm. In this review, we discuss rapidly evolving sequencing tools and technologies, such as error-corrected sequencing and single cell analysis, which we anticipate will fundamentally reshape the field. In addition, we consider a variety emerging applications for these new technologies, including the detection of DNA adducts, inference of mutational processes based on genomic site and local sequence contexts, and evaluation of genome engineering fidelity, as well as other cutting-edge challenges for the next 50 years of environmental and molecular mutagenesis research. Environ. Mol. Mutagen. 61:135-151, 2020. © 2019 The Authors. Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Polyguanine tracts (PolyGs) are short guanine homopolymer repeats that are prone to accumulating mutations when cells divide. This feature makes them especially suitable for cell lineage tracing, which has been exploited to detect and characterize precancerous and cancerous somatic evolution. PolyG genotyping, however, is challenging because of the inherent biochemical difficulties in amplifying and sequencing repetitive regions. To overcome this limitation, we developed PolyG-DS, a next-generation sequencing (NGS) method that combines the error-correction capabilities of duplex sequencing (DS) with enrichment of PolyG loci using CRISPR-Cas9-targeted genomic fragmentation. PolyG-DS markedly reduces technical artifacts by comparing the sequences derived from the complementary strands of each original DNA molecule. We demonstrate that PolyG-DS genotyping is accurate, reproducible, and highly sensitive, enabling the detection of low-frequency alleles (

Clonal hematopoiesis (CH aka CHIP) is the process whereby otherwise healthy individuals accumulate clones bearing low-frequency somatic mutations in hematologic malignancy (HM) driver genes. Historically thought to be a phenomenon of the elderly that is associated with increased risk of HM and cardiovascular disease, CH has been observed in increasingly younger cohorts as sequencing technology has advanced. CH in hematopoietic stem cell transplant (HSCT) donors has also been associated with adverse recipient outcomes. In this retrospective analysis, we asked whether CH mutations preferentially expanded in recipients vs. donors in paired samples collected a median of 34 years post-HSCT (range 6-45). 16 donor-recipient pairs were analyzed, 11 of which had a donor sample collected pre-HSCT. All individuals were healthy at the time of contemporary sampling. Median donor age was 26 years (range 12-84) at time of HSCT. We took advantage of the ultra-low error rate ( The mean DS error-corrected molecular depth across all samples was 26,731x for AML-associated genes and 32,812x in the neutral regions. CH was detected in 100% of donors and recipients, with a trend toward more CH clones in older vs. younger donors. We defined CH as mutations observed in at least 2 unique DNA molecules, rather than by a threshold variant allele frequency (VAF). Among the 11 donors with matched historical and contemporary DNA samples, CH mutations in the latter had a higher average VAF and comprised more unique mutations than in the former, reflecting clonal expansions over time. The magnitude of CH expansions was greater for the AML-associated genes than for neutral loci, consistent with positive selection. Post-HSCT the overall number of unique mutations was similar between donors and recipients. Most donor CH variants did not preferentially expand in the recipient. A minority of recipients had a modestly increased overall average VAF relative to paired donors in AML genes. While the VAF of some likely pathogenic variants increased in the recipient vs the donor, this was not the case for most CH mutations in most recipients and the reciprocal was sometimes observed. DS represents a powerful tool for the high-resolution study of CH evolution. In a unique historical cohort comprising longitudinal samples from some of the earliest patients and donors to undergo HSCT, we observed limited instances of accelerated CH in recipients relative to donors. While this cohort is biased due to enriched inclusion of long-term survivors, our data exemplifies how healthy hematopoiesis, indistinguishable between donors and recipients, can be maintained almost half a century following transplant. Citation Format: Masumi Ueda, Jake Higgins, Thomas Smith, Elizabeth Schmidt, Fang Y. Lo, Charles C. Valentine, Jesse J. Salk, Rainer Storb, Jerald P. Radich. Duplex Sequencing reveals ubiquitous clonal hematopoiesis and complex donor-recipient clonal dynamics following HSCT [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr LB063.

Significance Error-corrected next-generation sequencing (ecNGS) can be used to rapidly detect and quantify the in vivo mutagenic impact of environmental exposures or endogenous processes in any tissue, from any species, at any genomic location. The greater speed, higher scalability, richer data outputs, and cross-species and cross-locus applicability of ecNGS compared to existing methods make it a powerful new tool for mutational research, regulatory safety testing, and emerging clinical applications., The ability to accurately measure mutations is critical for basic research and identifying potential drug and chemical carcinogens. Current methods for in vivo quantification of mutagenesis are limited because they rely on transgenic rodent systems that are low-throughput, expensive, prolonged, and do not fully represent other species such as humans. Next-generation sequencing (NGS) is a conceptually attractive alternative for detecting mutations in the DNA of any organism; however, the limit of resolution for standard NGS is poor. Technical error rates (∼1 × 10−3) of NGS obscure the true abundance of somatic mutations, which can exist at per-nucleotide frequencies ≤1 × 10−7. Using duplex sequencing, an extremely accurate error-corrected NGS (ecNGS) technology, we were able to detect mutations induced by three carcinogens in five tissues of two strains of mice within 31 d following exposure. We observed a strong correlation between mutation induction measured by duplex sequencing and the gold-standard transgenic rodent mutation assay. We identified exposure-specific mutation spectra of each compound through trinucleotide patterns of base substitution. We observed variation in mutation susceptibility by genomic region, as well as by DNA strand. We also identified a primordial marker of carcinogenesis in a cancer-predisposed strain of mice, as evidenced by clonal expansions of cells carrying an activated oncogene, less than a month after carcinogen exposure. These findings demonstrate that ecNGS is a powerful method for sensitively detecting and characterizing mutagenesis and the early clonal evolutionary hallmarks of carcinogenesis. Duplex sequencing can be broadly applied to basic mutational research, regulatory safety testing, and emerging clinical applications.

The ability to accurately measure mutations is critical for basic research and identification of potential drug and chemical carcinogens. Current methods for in vivo quantification of mutagenesis are limited because they rely on transgenic rodent systems that are low-throughput, expensive, prolonged, and don’t fully represent other species such as humans. Next generation sequencing (NGS) is a conceptually attractive alternative for mutation detection in the DNA of any organism, however, the limit of resolution for standard NGS is poor. Technical error rates (~1×10−3) of NGS obscure the true abundance of somatic mutations, which can exist at pernucleotide frequencies ≤1×10−7. Using Duplex Sequencing, an extremely accurate error-corrected NGS (ecNGS) technology, we were able to detect mutations induced by 3 carcinogens in 5 tissues of 2 strains of mice within 31 days following exposure. We observed a strong correlation between mutation induction measured by Duplex Sequencing and the gold-standard transgenic rodent mutation assay. We identified exposure-specific mutation spectra of each compound through trinucleotide patterns of base substitution. We observed variation in mutation susceptibility by genomic region, as well as by DNA strand. We also identified the primordial signs of carcinogenesis in a cancer-predisposed strain of mice, as evidenced by clonal expansions of cells carrying an activated oncogene, less than a month after carcinogen exposure. These findings demonstrate that ecNGS is a powerful method for sensitively detecting and characterizing mutagenesis and the early clonal evolutionary hallmarks of carcinogenesis. Duplex Sequencing can be broadly applied to chemical safety testing, basic mutational research, and related clinical uses.SIGNIFICANCE STATEMENTError-corrected next generation sequencing (ecNGS) can be used to rapidly detect and quantify the in vivo mutagenic impact of environmental exposures or endogenous processes in any tissue, from any species, at any genomic location. The greater speed, higher scalability, richer data outputs, as well as cross-species and cross-locus applicability of ecNGS compared to existing methods make it a powerful new tool for mutational research, regulatory safety testing, and emerging clinical applications.

Mutations of ABL1 are the dominant mechanism of relapse in Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph + ALL). We performed highly accurate Duplex Sequencing of exons 4–10 of ABL1 on bone marrow or peripheral blood samples from 63 adult patients with previously untreated Ph + ALL who received induction with intensive chemotherapy plus a BCR-ABL1 TKI. We identified ABL1 mutations prior to BCR-ABL1 TKI exposure in 78% of patients. However, these mutations were generally present at extremely low levels (median variant allelic frequency 0.008% [range, 0.004%–3.71%] and did not clonally expand and lead to relapse in any patient, even when the pretreatment mutation was known to confer resistance to the TKI received. In relapse samples harboring a TKI-resistant ABL1 mutation, the corresponding mutation could not be detected pretreatment, despite validated sequencing sensitivity of Duplex Sequencing down to 0.005%. In samples under the selective pressure of ongoing TKI therapy, we detected low-level, emerging resistance mutations up to 5 months prior to relapse. These findings suggest that pretreatment ABL1 mutation assessment should not guide upfront TKI selection in Ph + ALL, although serial testing while on TKI therapy may allow for early detection of clinically actionable resistant clones.

Rationale: Accurate and sensitive detection of measurable residual disease (MRD) is critical for clinical management and the development of new therapies for pediatric acute myeloid leukemia (AML). The presence of MRD is a well-established predictor of disease relapse and mortality. Ideally, MRD detection would be both broadly applicable across patients and highly specific, thus enabling clinicians to make timely treatment decisions to improve outcomes. Background: In pediatric AML, unlike adult AML, there are few cross-patient recurrent somatic driver mutations. This makes it impractical to use a fixed gene panel for MRD monitoring. One solution is to perform disease tracking with patient-specific mutations that are present in individual AMLs at time of diagnosis (TOD). However, following modern therapies, should residual disease remain, the telltale mutations will be present in only a very small subpopulation of cells–far below the ~1% limit-of-detection of conventional NGS methods. Duplex sequencing (DS) is an ultra-accurate NGS error-correction method which increases the accuracy of base calls by more than 10,000-fold. Here we demonstrate the use of tumor-informed DS to detect patient-specific MRD in pediatric AML samples with unprecedented sensitivity and specificity. Methods: Marrow aspirates were collected from 13 pediatric AML patients at TOD, during treatment, remission, and at relapse. Among the cohort, 9/13 patients ultimately relapsed. Personalized hybrid capture panels were designed to target hundreds of somatic mutations that were present in each patient’s TOD sample, as identified by whole genome sequencing of tumor-normal pairs. DS was used for longitudinal monitoring of these mutations in blood or marrow at up to 5 time points following initial induction. Results: We detected MRD in one or more samples from all 9 patients who went on to relapse. However, we also observed high-confidence TOD variants in the last collected sample from all 5 patients who did not relapse during at least 24 months of clinical follow-up. As a technical control, we tested a subset of patient panels on a blood sample from a healthy subject and found no disease-specific mutations. Interestingly, we observed a marked degree of heterogeneity of response-to-therapy in the cohort. In some cases, all mutations logarithmically decreased uniformly, whereas in others there was clear dispersion of subclones over time, with some declining 100-times more quickly than others. Conclusions: While DS in conjunction with personalized panels represents a powerful approach for MRD tracking, simply the existence of very low frequency TOD mutations in pediatric AML is not specific enough for predicting relapse without additional interpretation. Lingering mutations in non-relapsers might represent preneoplastic field defects, pre-existing somatic mosaicism, or perhaps even true residual disease under immune suppression. Our study highlights an unappreciated degree of complexity that arises when using ultra-accurate mutation detection for cancer monitoring. Citation Format: Thomas H. Smith, Jake E. Higgins, Fang Y. Lo, Michael J. Hipp, Charles C. Valentine, Amanda R. Leonti, Rhonda E. Ries, Soheil Meshinchi, Benjamin J. Huang, Jesse J. Salk. Longitudinal monitoring of pediatric acute myeloid leukemia using duplex sequencing of patient-specific panels reveals ultra-low frequency MRD that marks persistent carcinogenesis and complex clonal evolution [abstract]. In: Proceedings of the AACR Special Conference on the Evolutionary Dynamics in Carcinogenesis and Response to Therapy; 2022 Mar 14-17. Philadelphia (PA): AACR; Cancer Res 2022;82(10 Suppl):Abstract nr A009.

Next-generation sequencing methods suffer from low recovery, uneven coverage, and false mutations. DNA fragmentation by sonication is a major contributor to these problems because it produces randomly sized fragments, PCR amplification bias, and end artifacts. In addition, oligonucleotide-based hybridization capture, a common target enrichment method, has limited efficiency for small genomic regions, contributing to low recovery. This becomes a critical problem in clinical applications, which value cost-effective approaches focused on the sequencing of small gene panels. To address these issues, we developed a targeted genome fragmentation approach based on CRISPR/Cas9 digestion that produces DNA fragments of similar length. These fragments can be enriched by a simple size selection, resulting in targeted enrichment of up to approximately 49,000-fold. Additionally, homogenous length fragments significantly reduce PCR amplification bias and maximize read usability. We combined this novel target enrichment approach with Duplex Sequencing, which uses double-strand molecular tagging to correct for sequencing errors. The approach, termed CRISPR-DS, enables efficient target enrichment of small genomic regions, even coverage, ultra-accurate sequencing, and reduced DNA input. As proof of principle, we applied CRISPR-DS to the sequencing of the exonic regions of TP53 and performed side-by-side comparisons with standard Duplex Sequencing. CRISPR-DS detected previously reported pathogenic TP53 mutations present as low as 0.1% in peritoneal fluid of women with ovarian cancer, while using 10- to 100-fold less DNA than standard Duplex Sequencing. Whether used as standalone enrichment or coupled with high-accuracy sequencing methods, CRISPR-based fragmentation offers a simple solution for fast and efficient small target enrichment.

Background: The presence of measurable residual disease (MRD) is strongly associated with treatment outcomes in acute myeloid leukemia (AML). However, > 20% of MRD negative cases (as assessed by flow cytometry) subsequently relapse. We sought to determine if ultrasensitive Duplex Sequencing (DS), which relies on double-stranded consensus-making to achieve an error rate below one-in-ten-million, yields better prognostic performance via molecular MRD detection. Methods: Retrospective targeted DNA sequencing of 29 genes recurrently mutated in adult AML was performed on paired diagnostic and remission bone marrow samples from patients enrolled on the SWOG trial S0106 (randomized 7+3 versus 7+3 + gemtuzumab ozogamicin (GO)). Patients were selected if they had remission samples with flow cytometry results (n=67). Non-error corrected sequencing was performed on diagnostic samples (average depth 279x) and DS was performed on remission samples (average duplex molecular depth 27,002x). For each patient, potential germline variants were identified and excluded from the analysis if the variant allele fraction (VAF) was ≥ 35% at both diagnosis and remission, or ≥ 40% at either time point and a gnomAD allele frequency ≥ 0.05. Somatic variants present at diagnosis were classified as potentially deleterious if computationally predicted as such and with a VAF ≥ 5% (≥ 1% for FLT3-ITD/NPM1 insertions). For analysis of residual disease in remission, we evaluated the following outcomes (events): overall survival (OS; death), relapse-free survival (RFS) and time to relapse (TTR; relapse with death a competing event). All outcomes were measured from date of morphologic remission to date of event, with patients without event censored at date of last contact. Associations between residual disease and outcomes were assessed using Cox regression models (cause-specific model for TTR). Results: The median age was 48 years (range 8-60). 32 patients were randomized to 7+3 and 30 to 7+3+GO. A total of 172 potentially deleterious variants were identified in the diagnostic samples. Variants had an average VAF of 31% (range 1.4-91.5%) at diagnosis and were detected in 23 of the 29 genes, with FLT3 being the most frequently mutated. Of the 67 patients analyzed, 93% (n=62) had at least one variant detected at diagnosis (median 2, range 0-9) and 68% (n=42) had at least one residual diagnostic variant also found in the remission sample. We defined the presence of DS MRD as non-DTA (DNMT3A, TET2, ASXL1) time-of-diagnosis mutations identified at the remission time point with a VAF > 0.1% and/or an NPM1 VAF > 0.01% (PMID:31860405). DS MRD was strongly associated with all outcomes, with hazard ratios (and 95% CI) for TTR: 7.1 (2.7-18.9); RFS: 4.9 (2.2-10.9) and OS: 5.1 (2.1-12.3). As a comparator, we correlated treatment outcomes with the results of a flow cytometry MRD assay previously carried out on the same samples during the S0106 trial. The prognostic association of flow MRD with TTR, RFS and OS was less strong (i.e., a smaller hazard ratio) than DS MRD, with TTR: 2.5 (0.9-6.7); RFS: 2.2 (0.9-5.4) and OS: 2.4 (1.0-6.1). RFS and TTR for DS MRD and flow cytometry are plotted below for the 62 patients with a variant detected at diagnosis. Comparing DS MRD with flow cytometry, discordance was found in 20 cases: 15 cases where DS was positive and flow negative, and 5 cases in the opposite direction. Among the 15 discordant cases with DS positive and flow negative, 9 relapsed and 2 died without relapse and 4 were alive without relapse at last contact. Among the 5 discordant cases with DS negative and flow positive, 1 relapsed, 1 died without relapse, and 3 were alive without relapse at last contact. Conclusions: Among the 67 patients evaluated in this prospectively collected study, the presence of MRD defined by DS was strongly associated with adverse disease outcomes. The vast majority of patients had at least one time-of-diagnosis mutation that could be tracked as a measure of MRD. When compared with flow cytometry, DS exhibited superior negative and positive predictive values for foretelling future relapse, although this could potentially reflect the historical version of the flow cytometry assay used. These findings suggest that DS is a powerful tool that could be used in patient management and for early treatment assessment in clinical trials. Figure 1 Figure 1. Disclosures Hourigan: Sellas: Research Funding. Radich: Amgen: Membership on an entity's Board of Directors or advisory committees; Novartis: Membership on an entity's Board of Directors or advisory committees; Genentech: Membership on an entity's Board of Directors or advisory committees; BMS: Membership on an entity's Board of Directors or advisory committees.

Sensitive and specific detection of measurable residual disease (MRD) after treatment in pediatric acute myeloid leukemia (AML) is prognostic of relapse and is important for clinical decision making. Mutation-based methods are increasingly being used, but are hampered by the limited number of common driver gene mutations to target as clone markers. Additional targets would greatly increase MRD detection power. However, even in cases with many AML-defining mutations, it is the limited accuracy of current molecular methods which establishes the lower bounds of sensitivity. Here we describe an ultrasensitive approach for disease monitoring with personalized hybrid capture panels targeting hundreds of somatic mutations identified by whole genome sequencing (WGS), and using extremely accurate Duplex Sequencing (DS) in longitudinal samples. In a pilot cohort of 13 patients we demonstrate detection sensitivities several orders of magnitude beyond currently available single locus testing or less accurate sequencing. With multi-target panels, overall power for MRD detection is cumulative across sites. For example, if a patient has MRD at a true frequency of 1/30,000, sequencing a single mutant site to 10,000x molecular depth would be unlikely to detect MRD. However, sequencing 10 sites each to 10,000x would effectively total 100,000x informative site depth, increasing power to >95%. However, standard sequencing assays are insufficiently accurate to achieve this theoretical limit of detection (LOD). DS enables accurate detection of individual variants to Marrow aspirates were collected from 13 uniformly treated pediatric AML patients at time of diagnosis (TOD), during treatment (end of induction, EOI), in remission (end of therapy, EOT), and at relapse. 9/13 patients relapsed. DNA from TOD was analyzed by WGS. Germline variants were excluded and somatic single nucleotide variants (SNVs) were targeted by a custom probe panel designed for each patient. An average of 170 SNVs were targeted per patient (range 53-200). More than 90% of the SNVs were noncoding. Longitudinal samples were then analyzed with DS, which compares sequences from both strands of each DNA molecule to eliminate technical noise and reveal biological mutation signal with extreme accuracy and sensitivity. A median of 82% of WGS SNVs were validated by DS in the TOD DNA, and the vast majority of those were also present at relapse. Relapsers had more SNVs at diagnosis than non-relapsers. EOT samples were sequenced to an average Duplex molecular depth of 29,400x, with a maximum of 61,283x. The figure below shows time course plots tracking SNVs at diagnosis, EOT and relapse for 2 patients. Among mutations validated in TOD samples, a median of only 8 (5%) were detected per EOT sample (range 0-66 mutations). MRD was detected in 8/9 relapsers. Targeting 1 or even 10 SNVs would therefore have missed MRD in the majority of these patients. Among relapsers, median EOT SNV VAF was 0.069%. The lowest single VAF detected per EOT sample ranged from 0.036% to 0.002%. The presence of an SNV at diagnosis and relapse implies that it must truly be present at EOT, whether or not it is detected. Therefore, if a small minority of leukemic mutations are detected at EOT, the true overall MRD frequency is much lower than the LOD at any single site. In the only relapser where MRD was not detected, targeted SNVs were present at diagnosis and relapse, so additional sequencing depth at EOT would eventually reveal ultra-low frequency mutations. Among patients that did not relapse by the end of the study, median VAF at EOI (the latest time point DNA available) was 0.0258%. Therefore, non-relapsers have a lower median VAF at EOI than relapsers do even later at EOT, potentially indicating very early on that treatment is more successful. This study shows excellent performance of DS-based assays for detecting MRD with patient-specific panels. We have demonstrated that among large panels of validated somatic SNVs present at time of diagnosis, a median of 5% are identified at EOT in eventual relapsers. DS detected MRD in 8/9 patients, and at a median VAF well below the limit of detection of any other sequencing technology. Comprehensive personalized hybrid selection panels coupled with DS represents a powerful option for MRD monitoring in pediatric AML and potentially other cancers. Figure Disclosures Higgins: TwinStrand Biosciences: Current Employment. Lo:TwinStrand Biosciences: Current Employment. Hipp:TwinStrand Biosciences: Current Employment. Valentine:TwinStrand Biosciences: Current Employment. Williams:TwinStrand Biosciences: Current Employment. Radich:TwinStrand Biosciences: Research Funding. Salk:TwinStrand Biosciences: Current Employment.

Significance Cancers evolve many mutations. Clonal driver mutations are selected early. Subsequent evolution occurs in a branching fashion, possibly without selection (“neutral evolution”). Rarer mutations occur later on smaller branches of the evolutionary tree. Using a DNA-sequencing method, duplex sequencing, with unprecedented accuracy and sensitivity, we quantified rare unique subclonal mutations in diagnostic specimens from 5 human colorectal cancers. Rarer subclones probe later evolutionary time points than previously possible. We confirm neutral evolution at later times and find many more subclonal mutations than expected. A theoretical method allowed us to extrapolate further forward in time to diagnosis. At diagnosis, every DNA base is mutated in at least one cancer cell. In particular, any therapy resistance mutation would be present., Human colorectal cancers (CRCs) contain both clonal and subclonal mutations. Clonal driver mutations are positively selected, present in most cells, and drive malignant progression. Subclonal mutations are randomly dispersed throughout the genome, providing a vast reservoir of mutant cells that can expand, repopulate the tumor, and result in the rapid emergence of resistance, as well as being a major contributor to tumor heterogeneity. Here, we apply duplex sequencing (DS) methodology to quantify subclonal mutations in CRC tumor with unprecedented depth (104) and accuracy (

Introduction/Background Poor survival of high-grade serous pelvic cancer (HGSC) is caused by a lack of effective screening measures. Uterine and tubal lavage (UtL) is the sampling technique for HGSC detection with the highest sensitivity ever reported. We combined UtL with the highest accuracy sequencing technology for TP53 mutations, Duplex Sequencing (DS) in order to develop a test potentially feasible for HGSC screening. Methodology In study I, the clinical feasibility of UtL using a new catheter was examined in 93 patients undergoing surgery for suspected HGSC. In study II, pain at UtL under local anesthesia was evaluated in 22 healthy women. In Study III, uterine lavages from 10 HGSC patients and 11 controls were subjected to DS to identify TP53 mutations. Results In study I, UtL was performed successfully in 92 (98.9%) of 93 cases by 16 gynecologists. In 84.8% of patients no cervical dilatation was required. For benign conditions, the need for dilatation was associated with menopause status (OR, 4.929; 95% CI, 1.439–16.884; P=0.016). In study II, pain during UtL was rated with a median VAS score of 1.6. During 4 weeks after UtL, no serious complications occurred. UtL took 6.5 minutes on average. The amount of extracted DNA was above the lower limit for sensitive, deep-sequencing mutation analysis in all cases. In study III, in 80% of all HGSC tumor-derived TP53 mutations in UtL samples were detected by DS. However, we detected low frequency TP53 mutations in nearly all samples tested. Mutations increased with age and shared the selection traits of clonal TP53 mutations commonly found in human tumors. However, tumor-derived TP53 mutations were present at frequencies at least 10-fold higher than background mutations. Conclusion The UtL concept fulfils all prerequisites to be used in a potential screening setting. Interpretation of DS results needs to consider TP53 biologic background mutations. Disclosure Philipp Harter COI: Honoraria: Astra Zeneca, Roche, Sotio, Tesaro, Stryker, ASCO, Zai Lab, MSD; Advisory Board: Astra Zeneca, Roche, Tesaro, Lilly, Clovis, Immunogen, MSD/Merck; Research funding (Inst): Astra Zeneca, Roche, GSK, Boehringer Ingelheim, Medac, DFG, European Union, DKH, Tesaro, Genmab Rosana Risques COI: RAR shares equity in NanoString Technologies Inc. and is the principal investigator on an NIH SBIR subcontract research agreement with TwinStrand Biosciences Inc. Florian Heitz COI: Honoraria: AstraZeneca, Roche, Tesaro, Clovis; Advisory Board: Astra Zeneca, Roche, Tesaro, Clovis Christoph Grimm COI: Consultant: AstraZeneca, Celgene, MSD, PharmaMar, Roche, Tesaro, Vifor Pharma Speaker: Amgen, AstraZeneca, MSD, PharmaMar, Roche, Tesaro Direct research funding: Meda Pharma, Roche Diagnostics Ignace Vergote COI: CONSULTING = ADVISORY BOARD: Advaxis, Inc., Eisai Inc., MSD Belgium, Roche NV, Genmab, F. Hoffmann-La Roche Ltd, Pharmamar, Millennium Pharmaceuticals, Clovis Oncology Inc., AstraZeneca NV, Tesaro, Oncoinvent AS, Immunogen Inc, Sotio CONTRACTED RESEARCH (via KULeuven): Oncoinvent AS, Genmab GRANT = CORPORATE SPONSORED RESEARCH: Amgen, Roche, Stichting tegen Kanker Accomodations, travel Fabian Trillsch COI:Consultant: AstraZeneca, Roche, Tesaro Speaker/Travel Expenses: AstraZeneca, Medac, PharmaMar, Roche, Tesaro Paul Speisr: Share Holder OVARTEC GmbH.

Human colorectal cancers (CRC) contain both clonal and subclonal mutations. Clonal mutations are positively selected, present in most cells and drive malignant progression. Subclonal mutations are randomly dispersed throughout the genome; they provide a vast reservoir of mutant cells that can expand, repopulate the tumor and result in the rapid emergence of resistance, as well as being a major contributor to tumor heterogeneity. Here, we apply Duplex Sequencing (DS) methodology to quantify subclonal mutations in CRC tumor with unprecedented depth (104) and accuracy (−7). We measured mutation frequencies in genes encoding replicative DNA polymerases and in genes frequently mutated in CRC, and found an unexpectedly high effective mutation rate, 7.1 × 10−7. The curve of subclonal mutation accumulation as a function of sequencing depth, using DNA obtained from five different tumors, is in accord with a neutral model of tumor evolution. We present a new theoretical approach to model neutral evolution independent of the infinite sites assumption (which states that a particular mutation arises only in one tumor cell at any given time). Our analysis indicates that the infinite sites assumption is not applicable once the number of tumor cells exceeds the reciprocal of the mutation rate, a circumstance relevant to even the smallest clinically diagnosable tumor. Our methods allow accurate estimation of the total mutation burden in clinical cancers. Our results indicate that no DNA locus is wild type in every malignant cell within a tumor at the time of diagnosis (probability of all cells wild type = 10−308).Significance StatementCancers evolve many mutations. Clonal mutations are selected early. Subsequent evolution occurs in a branching fashion, possibly without selection (“neutral evolution”). Rarer mutations occur later on smaller branches of the evolutionary tree. Using a DNA sequencing method, duplex sequencing, with unprecedented accuracy and sensitivity, we quantified very rare unique subclonal mutations in diagnostic specimens from five human colorectal cancers. Rarer subclones probe later evolutionary timepoints than previously possible. We confirm neutral evolution at later times and find many more subclonal mutations than expected. A novel theoretical method allowed us to extrapolate further forward in time to diagnosis. At diagnosis, every base in DNA is mutated in at least one cancer cell. In particular any therapy resistance mutation would be present.

Aristolochic acids (AA) are a potent class of mutagen that occur in plants of the genera Aristolochia and Asarum. Exposures to AA are caused both by intentional ingestion of traditional herbal medicines containing Aristolochia and through inadvertent consumption from contaminated crops. The carcinogenic potential and epidemiologic connection of AA exposure with upper-tract urothelial cancers (UTUC), liver cancer and kidney failure have been known for years; however, the extent of AA exposure and its influence on cancer risk and mortality have yet to be fully explored in many world populations. Using whole exome sequencing of tumor samples, a recent study determined that 78% of HCCs from Taiwan exhibit distinct mutation signatures of AA exposure (Ng, Poon, Huang, et al., 2017). Exposure to AA is geographically widespread and there is a major opportunity for primary and secondary prevention of AA-associated diseases if it were possible to detect the exposure through non-invasively sampled body fluids such as blood or urine sediment. Although AA adducts can be measured in body fluids, this measurement only reflects recent ingestion, not the life-integrated exposure that better approximates disease risk. A mutation burden-based assay could circumvent this problem; however, a major challenge is that AA-mediated mutations are extremely rare in bulk tissue or fluid samples from non-tumor sources. Duplex Sequencing (DS) is a tag-based NGS error-correction method that enables the detection of ultra-rare mutations that exist at the mutation frequency of normal tissues (~10e-8). DS can be applied to DNA from any tissue of any organism without the need for tumor formation or ex vivo cloning. We developed a DS-based assay that measures mutations across a 50 kb genome-representative panel and applied it to the tissues and liquid biopsies of cancer patients of known AA-exposure status. We sampled the blood, urine, kidney, normal-adjacent urothelium and ureteral tumors from 12 UTUC patients where 6 of the patients were exposed to AA, and 6 were not. Under blinded conditions, we positively identified AA mutation signatures in the normal adjacent urothelium tissues of 5 AA-exposed patients. Remarkably, we also detected AA signatures in all available blood and urine samples from the same patients. The proportion of mutations attributed to AA-exposure ranged from 2.5% in blood to 67.8% in kidney. These mutations, on the order of one-in-a-million, were exposure-related, not tumor derived. We did not detect AA mutation signatures in any samples from the non-AA-exposed patients. To our knowledge, this is the first time an AA mutation signature has been detected in non-invasively sampled body fluids. DS is a promising approach for the detection of mutagenic signatures caused by environmental carcinogens, and we foresee it being a powerful metric that can be used to measure life-integrated carcinogenic processes and cancer risk. Citation Format: Fang Yin Lo, Charles Valentine, Elizabeth Schmidt, Lindsey Williams, Arnoud Boot, Steve Rozen, Jesse Salk. Non-invasive detection of aristolochic acid exposure using ultra-sensitive duplex sequencing [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3147.

Hematologic malignancies (HM) are driven by clonal evolution—from preneoplastic clones through acquired drug resistance. Clonal hematopoiesis of indeterminate potential (CHIP) is a process whereby healthy individuals accumulate low-frequency clones with driver gene mutations. CHIP is associated with increased risk of HM but is challenging to detect. Measurable residual disease (MRD) after therapy in acute myeloid leukemia (AML) and other HM is a strong predictor of relapse; however, current MRD assays either have limited sensitivity or are narrowly applicable. Tyrosine kinase inhibitors (TKI) target BCR-ABL1 in chronic myeloid leukemia (CML) and Ph+ acute lymphocytic leukemia (ALL), yet TKI-resistance (TKI-R) mutations can drive relapse. Early detection of MRD and TKI-R mutations can inform important clinical interventions. Next-generation sequencing (NGS) can detect HM subclones, but not below a frequency of approximately 1%. Common error correction methods have difficulty detecting mutations below 0.1%. Duplex sequencing (DS) is a sophisticated NGS error correction assay capable of detecting individual mutations down to ultralow frequencies of at least 0.001% with a background error frequency below 1 in 10 million bases. We used DS for three clinically relevant HM applications: (1) CHIP: We performed DS on DNA from 10 hematopoietic cell transplant (HCT) donor/recipient pairs 7-46 years post-HCT. We targeted 29 genes recurrently mutated in AML, many of which have been described in CHIP. DS error-corrected depth up to 45,402x was generated. We identified dozens of clonal exonic variants (CEV) in every donor and recipient, most nonsynonymous, down to variant allele frequency (VAF) of 5.7E-5. The majority of the 154 shared donor/recipient CEVs were present at VAF This abstract is also being presented as Poster A48. Citation Format: Jake Higgins, Lindsey N. Williams, Charles C. Valentine, Gabriel Pratt, Rashmi Kanagal-Shamanna, Nicholas J. Short, Masumi Ueda, Rainer Storb, Jerald P. Radich, Jesse J. Salk. Clonal landscapes of hematologic malignancies redefined by ultrasensitive duplex sequencing [abstract]. In: Proceedings of the AACR Special Conference on Advances in Liquid Biopsies; Jan 13-16, 2020; Miami, FL. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(11_Suppl):Abstract nr PR04.

Summary: High-accuracy next-generation DNA sequencing promises a paradigm shift in early cancer detection by enabling the identification of mutant cancer molecules in minimally invasive body fluid samples. We demonstrate 80% sensitivity for ovarian cancer detection using ultra-accurate Duplex Sequencing to identify TP53 mutations in uterine lavage. However, in addition to tumor DNA, we also detect low-frequency TP53 mutations in nearly all lavages from women with and without cancer. These mutations increase with age and share the selection traits of clonal TP53 mutations commonly found in human tumors. We show that low-frequency TP53 mutations exist in multiple healthy tissues, from newborn to centenarian, and progressively increase in abundance and pathogenicity with older age across tissue types. Our results illustrate that subclonal cancer evolutionary processes are a ubiquitous part of normal human aging, and great care must be taken to distinguish tumor-derived from age-associated mutations in high-sensitivity clinical cancer diagnostics. : Salk et al. demonstrate that ultra-sensitive DNA sequencing to identify TP53 mutations among cells shed into uterine fluid shows potential for minimally invasive ovarian cancer detection. Yet they also reveal ubiquitous age-related accumulations of cancer-like TP53 mutations in the normal tissues of healthy women. This highlights an important challenge of using tumor driver mutations for cancer screening. Keywords: Duplex Sequencing, next-generation sequencing, TP53, clonal evolution, early detection, somatic mutations, high-grade serous ovarian cancer, uterine lavage, aging, gynecologic oncology

High accuracy next-generation DNA sequencing promises a paradigm shift in early cancer detection by enabling the identification of mutant cancer molecules in minimally-invasive body fluid samples. We demonstrate 80% sensitivity for ovarian cancer detection using ultra-accurate Duplex Sequencing to identify TP53 mutations in uterine lavage. However, in addition to tumor DNA, we also detect low frequency TP53 mutations in nearly all lavages from women with and without cancer. These mutations increase with age and share the selection traits of clonal TP53 mutations commonly found in human tumors. We show that low frequency TP53 mutations exist in multiple healthy tissues, from newborn to centenarian, and progressively increase in abundance and pathogenicity with older age across tissue types. Our results illustrate that subclonal cancer evolutionary processes are a ubiquitous part of normal human aging and great care must be taken to distinguish tumor-derived, from age-associated mutations in high sensitivity clinical cancer diagnostics.

Despite the remarkable throughput of next-generation sequencing technologies, standard techniques are limited by the difficulty in distinguishing sequencing errors from genuine low-frequency DNA variants within heterogeneous cellular or molecular populations. This Review discusses sequencing methodologies and bioinformatic strategies that have been devised for the reliable detection of rare mutations and describes various important applications in diverse fields including cancer, ageing and metagenomics. Mutations, the fuel of evolution, are first manifested as rare DNA changes within a population of cells. Although next-generation sequencing (NGS) technologies have revolutionized the study of genomic variation between species and individual organisms, most have limited ability to accurately detect and quantify rare variants among the different genome copies in heterogeneous mixtures of cells or molecules. We describe the technical challenges in characterizing subclonal variants using conventional NGS protocols and the recent development of error correction strategies, both computational and experimental, including consensus sequencing of single DNA molecules. We also highlight major applications for low-frequency mutation detection in science and medicine, describe emerging methodologies and provide our vision for the future of DNA sequencing.

Cumulative evidence indicates that a significant proportion of cancer evolution may occur before the development of histological abnormalities. While recent improvements in DNA sequencing technology have begun to reveal the presence of these early preneoplastic clones, the concept of 'premalignant field' was already introduced by Slaughter more than half a century ago. Also referred to as 'field effect', 'field defect' or 'field cancerization', these terms describe the phenomenon by which molecular alterations develop in normal-appearing tissue and expand to form premalignant patches with the potential to progress to dysplasia and cancer. Field effects have been well-characterized in ulcerative colitis, an inflammatory bowel disease that increases the risk of colorectal cancer. The study of the molecular alterations that define these fields is informative of mechanisms of tumor initiation and progression and has provided potential targets for early cancer detection. Herein, we summarize the current knowledge about the molecular alterations that comprise the field effect in ulcerative colitis and the clinical utility of these fields for cancer screening and prevention.

Hematopoietic cell transplantation (HCT) is used in the treatment of hematologic malignancies (HM). Clonal hematopoiesis of indeterminate potential (CHIP) is the age-associated process whereby healthy individuals accumulate low-frequency clones driven by mutations in genes recurrently mutated in myeloid cancers; CHIP is associated with an increased risk of HM and cardiovascular disease. It has been hypothesized that CHIP clones present in an HCT donor might disproportionately proliferate in the recipient due to a relative growth advantage during immune reconstitution. Here we use Duplex Sequencing (DS) to investigate, years after HCT, whether the clonal makeup of CHIP in prior donors differs from that of the donor-derived hematopoietic system long after allogeneic transfer to a recipient. DS is an ultra-sensitive NGS error-correction technology that relies on sequencing of the individual paired strands of original source DNA molecules to dramatically reduce technical background compared to other NGS methodologies. Here we performed DS on 10 pairs of donor/recipient peripheral blood DNA samples collected 7-46 years post-transplant, as well as DNA from a healthy young control. We used a sequencing panel comprising 29 genes recurrently mutated in acute myeloid leukemia (AML) and generated mean DS molecular depth of 27,580x per sample, with maximum DS depth reaching 45,402x. We detected an average of 356 and 283 total variants per donor and recipient, respectively, and 462 total variants in the control. These included 291 single-count exonic variants (SEV) in the control and an average of 215 SEVs in donors (Fig. A). To enrich for CHIP (clonal) variants, we filtered out single-count (non-clonal) events, as well as intronic (non-splice site) variants and germline SNPs. Whereas the young control had only 5 clonal exonic variants (CEV), meaning nearly all of the variants in the control were SEVs, SNPs or intronic (Fig. A,B), HCT donors and recipients respectively averaged 44 and 38 CEVs. 91% of CEVs were non-synonymous (Fig. B), which is higher than expected by chance. The average age of the donors was 36 years old at the time of HCT (range 12-84), and the oldest donor (M1) had the most CEVs. 24/29 genes in the panel exhibited CEVs in the control or HCT donors (Fig. C). DNMT3A, TET2, BRINP3 (FAM5C) and ASXL1 had the greatest number of CEVs and this imbalance was not explained purely by gene size. The young control sample had an overall somatic mutation frequency (MF) of 4.7x10-7. HCT donors and recipients had average MFs of 1.0x10-6 and 1.7x10-6, respectively. 154 CEVs were shared in at least 1 donor-recipient pair, with the majority present at variant allele frequencies (VAFs) 10-fold, all in DNMT3A, TET2 and RAD21. Only 1 CEV decreased >10-fold: TP53 R273H (4.9x10-4 to 4.2x10-5). Analysis of the same 10 donor-recipient pairs using a commercial sequencing service found that, with a VAF cutoff of 2% and filtering for SNPs, only 7 shared CEVs could be detected that increased or decreased in any recipient. Our ultra-high accuracy DS data indicates that the structure of CHIP is more polyclonal than has been previously reported. In a handful of older healthy donors we detected dozens of concurrent CHIP mutations, with a relative gene distribution nearly identical to that previously reported in studies on thousands of patients using less sensitive method. We observed a similar age-dependent effect. Furthermore, our findings indicate that the clonal dynamics of CHIP clones between HCT donors and recipients are substantially more complex than previously known. Nearly all somatic variants shared between HCT donors and recipients were observed at Figure Disclosures Higgins: TwinStrand Biosciences: Employment. Pratt:TwinStrand Biosciences: Employment. Valentine:TwinStrand Biosciences: Employment. Williams:TwinStrand Biosciences: Employment. Radich:Novartis: Other: RNA Sequencing; TwinStrand Biosciences: Research Funding. Salk:TwinStrand Biosciences: Employment.

Background: Some previous studies in Ph+ ALL suggested that T315I mutations of the ABL1 KD are present in many pts prior to tyrosine kinase inhibitor (TKI) treatment. However, these reports used RT-PCR to generate complementary DNA (cDNA) prior to analysis, which may introduce random errors and lead to false positive results. We hypothesized that ultrasensitive DS of genomic DNA (gDNA) would more accurately detect low-level pretreatment ABL1 mutations in Ph+ ALL because it does not rely upon error prone RT-PCR. Methods: DS compares the nucleotide sequences of each strand of double-stranded molecules, improving the accuracy of conventional next-generation sequencing by more than 10,000-fold. DS of exons 4-10 of ABL1 was performed to an average molecular depth of >10,000x. A mixture of known ABL1 variants at varying low variant allelic frequencies (VAFs) versus a negative control were used to determine sensitivity and specificity of DS. DS of ABL1 was performed on 64 pts with newly diagnosed Ph+ ALL prior to receiving frontline hyper-CVAD plus a TKI (imatinib, n=5; dasatinib, n=38; ponatinib, n=21). DS was also performed on cDNA samples from RT-PCR to assess the RT-induced error rate. On RNA from relapse samples, the KD (codons 221 through 500) of BCR-ABL1 was sequenced by the Sanger method, using a nested PCR approach, with a detection limit of 10-20%. Results: A total of 115 pretreatment ABL1 mutations were detected by DS, with ≥1 mutation detected in 47/64 pts (73%). These mutations were generally present at very low levels (median VAF 0.008% [range, 0.004%-0.649%]). Eleven TKI resistance (TKI-R) mutations (10% of all detected mutations) were identified by DS in 7 pts (Table 1). Eighteen pts (28%) relapsed; however, none of the pts with TKI-R mutations relapsed, despite 5 pts receiving a TKI at least intermediately resistant to the detected mutation. Using Sanger sequencing, TKI-R mutations were detected at relapse in 9 pts (T315I, n=6; F317I, V229L and V338G, n=1 each). None of these TKI-R mutations was detected by DS in pretreatment samples. Together, these results suggest that low-level pretreatment ABL1 KD mutations are not clinically meaningful in Ph+ ALL. To validate the DS methodology, we examined its sensitivity and specificity using 6 gDNA samples containing ABL1 KD mutations (9 different SNVs, including 6 TKI-R) mixed with control gDNA from a healthy young individual for predicted VAFs of 1/250 to 1/25,000. Mutations were identified by DS with 100% sensitivity down to VAF In the control samples, the background mutation frequency of ABL1 was ~5x10-7. Using a null model of random distribution of mutations in normal samples at this frequency, there was no significant difference in the proportion of TKI-R mutations detected with DS versus the expected background frequency, suggesting that there is no biased detection of TKI-R mutations by DS. To assess whether cDNA may be susceptible to artefactual mutations induced by RT steps, DS was performed on cDNA from 6 pts in whom no pretreatment ABL1 KD mutation was detected by DS. cDNA generated by RT had a 50-100x increased mutation frequency relative to the normal gDNA samples (Figure 2). False TKI-R mutations were observed in all cDNA samples (ranging from 11-16 TKI-R variants per sample). Overall, the fraction of TKI-R variants relative to total variants detected was significantly higher in all 6 cDNA samples compared to the negative control. Conclusion: DS had 100% sensitivity and specificity for detection of TKI-R mutations in ABL1. In contrast, RT-generated cDNA increased mutation frequency by 50- to 100-fold, suggesting that RT-based methods may falsely assess the presence of TKI-R mutations. Prior to TKI treatment, DS detected ABL1 mutations at a rate similar to healthy controls, and even in cases when a low-level TKI-R mutation was detected by DS, there was no association with relapse. Together, these findings suggest that pretreatment ABL1 mutations in Ph+ ALL are due to random age-related mutagenesis and do not represent clinically meaningful subclones that contribute to relapse. The superior accuracy of DS compared to RT-PCR has broad implications in the detection of low-level mutations across cancers, which may allow for early interventions to target emerging resistant subclones. Disclosures Short: Takeda Oncology: Consultancy, Research Funding; AstraZeneca: Consultancy; Amgen: Honoraria. Kantarjian:Ariad: Research Funding; AbbVie: Honoraria, Research Funding; Takeda: Honoraria; Agios: Honoraria, Research Funding; Amgen: Honoraria, Research Funding; Cyclacel: Research Funding; Daiichi-Sankyo: Research Funding; Astex: Research Funding; Jazz Pharma: Research Funding; Pfizer: Honoraria, Research Funding; Immunogen: Research Funding; BMS: Research Funding; Actinium: Honoraria, Membership on an entity's Board of Directors or advisory committees; Novartis: Research Funding. Sasaki:Otsuka: Honoraria; Pfizer: Consultancy. Ravandi:Macrogenix: Consultancy, Research Funding; Cyclacel LTD: Research Funding; Menarini Ricerche: Research Funding; Selvita: Research Funding; Xencor: Consultancy, Research Funding; Amgen: Honoraria, Membership on an entity's Board of Directors or advisory committees, Research Funding. Cortes:Astellas Pharma: Consultancy, Honoraria, Research Funding; Immunogen: Consultancy, Honoraria, Research Funding; Sun Pharma: Research Funding; Merus: Consultancy, Honoraria, Research Funding; BiolineRx: Consultancy; Bristol-Myers Squibb: Consultancy, Research Funding; Biopath Holdings: Consultancy, Honoraria; Takeda: Consultancy, Research Funding; Pfizer: Consultancy, Honoraria, Research Funding; Jazz Pharmaceuticals: Consultancy, Research Funding; Forma Therapeutics: Consultancy, Honoraria, Research Funding; Novartis: Consultancy, Honoraria, Research Funding; Daiichi Sankyo: Consultancy, Honoraria, Research Funding. Konopleva:Reata Pharmaceuticals: Equity Ownership, Patents & Royalties; Astra Zeneca: Research Funding; Agios: Research Funding; Genentech: Honoraria, Research Funding; Ascentage: Research Funding; Kisoji: Consultancy, Honoraria; Calithera: Research Funding; Stemline Therapeutics: Consultancy, Honoraria, Research Funding; Forty-Seven: Consultancy, Honoraria; Eli Lilly: Research Funding; AbbVie: Consultancy, Honoraria, Research Funding; Cellectis: Research Funding; Amgen: Consultancy, Honoraria; Ablynx: Research Funding; F. Hoffman La-Roche: Consultancy, Honoraria, Research Funding. Garcia-Manero:Novartis: Research Funding; Amphivena: Consultancy, Research Funding; Helsinn: Research Funding; AbbVie: Research Funding; Celgene: Consultancy, Research Funding; Astex: Consultancy, Research Funding; Onconova: Research Funding; H3 Biomedicine: Research Funding; Merck: Research Funding. Pratt:TwinStrand Biosciences: Employment. Williams:TwinStrand Biosciences: Employment. Valentine:TwinStrand Biosciences: Employment. Salk:TwinStrand Biosciences: Employment. Radich:TwinStrand Biosciences: Research Funding; Novartis: Other: RNA Sequencing. Jabbour:Pfizer: Consultancy, Research Funding; AbbVie: Consultancy, Research Funding; Takeda: Consultancy, Research Funding; Amgen: Consultancy, Research Funding; Adaptive: Consultancy, Research Funding; BMS: Consultancy, Research Funding; Cyclacel LTD: Research Funding.

Cancer is a disease of somatic evolution, characterized by the natural selection of genomic mutations that facilitate enhanced cell survival and proliferation. Although many thousands of tumor genomes have now been sequenced, our ability to identify early genetic patterns of clonal selection in both humans and model organisms have been hampered by inadequately sensitive methodologies for identifying mutations during the long period between their occurrence and the final outgrowth of a clinically apparent tumor. Sensitive molecular tests are capable of measuring minute levels of cancer cells in tissue samples, however no existing method is satisfactory at scaling to high-throughput detection at the rate of one cancer-associated variant per more than a million normal cells. NGS is used to study the genetic variation in cell populations, although the accuracy of standard NGS methods limit our ability to detect sub-clones below ~1%. Duplex Sequencing (DS) is a sensitive NGS error-correction method which increases the accuracy of base calls by more than 100,000-fold. DS enables precise reconstruction of the original double-stranded source molecule while overcoming both chemical and technical artifacts that arise during library preparation and sequencing. Here we present the use of Duplex Sequencing, in both human and mouse tissues, for measuring sub-clones at allelic fractions less than 1×10-4 for an early glimpse into pre-neoplastic evolution. We illustrate how, merely weeks after mutagen exposure, we observe emerging clones of cells bearing cancer-driving mutations in histologically normal mouse tissue. Furthermore, we illustrate how similar patterns of clonal selection can be seen in multiple otherwise healthy tissues of humans as part of normal aging. We discuss how variations in the pattern and rate of accumulation of rare cancer-associated mutations offers a new preclinical tool for evaluating carcinogenicity of chemicals, as well as a potential clinical tool for assessing life-integrated carcinogenic processes and cancer risk in humans. Citation Format: Charles C. Valentine, Mark Fielden, Robert Young, Jake Higgins, Lindsey Williams, Tan Li, Rohan Kulkarni, Sheroy Minocherhomji, Rosana Risques, Patrick Danaher, Jesse Salk. Duplex Sequencing detects rare subclonal variants that mark early carcinogenesis and preneoplastic clonal evolution [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4649.

The purpose of this study is to develop a universally applicable Duplex Sequencing (DS) assay for extremely sensitive detection of minimal residual disease (MRD) in acute myeloid leukemia (AML). The presence of MRD plays a critical role in AML recurrence and mortality, and is an important prognostic marker to guide clinical care and clinical trial risk stratification. Unfortunately, current MRD detection methods, including light microscopy, cytogenetics, qPCR, and flow cytometry have shortcomings in terms of sensitivity, general applicability, and/or inter-laboratory reproducibility. Duplex Sequencing eliminates technical artefacts observed with other sequencing technologies and delivers an unprecedented level of accuracy and sensitivity, on widely available NGS platforms. Here, we report validation of our DS MRD mutation detection assay, with a background mutation level more than 100-fold below single-stranded tag-based error correction (SSCS). Our complete DS MRD probe panel targets mutation hotspots or full coding sequences in 29 genes recurrently mutated in AML, representing mutations found in approximately 90% of patients. Supplementing with a targeted RNA-seq panel for recurrent gene fusions would approach 100%. To simulate low frequency mutations in MRD, we generated 2 DNA dilution series: either AML patient DNA or cell line DNA with known AML mutations that were serially diluted into normal human DNA. Replicate libraries were prepared with custom targeted probe panels and sequenced to over 50,000x or over 1,000,000x total Duplex depth (original Duplex DNA molecules sequenced per site), respectively. Expected mutations were detected down to MAFs of 1.9x10-5 and 6.0x10-6. Background mutation levels of 5.7x10-7 and 5.8x10-7, were respectively seen in each series, and likely reflect the true age-associated biological background, rather than technical noise. SSCS artefacts in these samples were approximately 100x higher than the frequency of DS mutations, and obscured expected mutations in the lower dilutions. We are further validating our method in both retrospective and prospective studies. In retrospective samples from a completed clinical trial, we are comparing the sensitivity and accuracy of DS vs. multicolor flow cytometry for MRD with known clinical outcomes. The relative prognostic value of DNA from blood vs. bone marrow will be quantified with DS as well. We will also perform DS on prospective samples to examine the prognostic value of specific genes and MAFs to eliminate uninformative variants, which will also allow accurate quantification of low level mutations in non-leukemic clones that accumulate with aging and can confound detection of MRD. Elimination of sequencing artefacts is critical for the application of NGS assays to MRD detection, especially when expected mutations are unknown, and DS is an important new tool in high sensitivity cancer diagnostics. Citation Format: Jacob E. Higgins, Lindsey N. Williams, Christopher S. Hourigan, Jerald P. Radich, Jesse J. Salk. Duplex sequencing for MRD detection in acute myeloid leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 422.

Background: Patients with ulcerative colitis (UC) are at risk of developing colorectal cancer. We have previously reported that cancer progression is associated with the presence of clonal expansions and shorter telomeres in nondysplastic mucosa. We sought to validate these findings in an independent case–control study. Methods: This study included 33 patients with UC: 14 progressors (patients with high-grade dysplasia or cancer) and 19 nonprogressors. For each patient, a mean of 5 nondysplastic biopsies from proximal, mid, and distal colon were assessed for clonal expansions, as determined by clonal length altering mutations in polyguanine tracts, and telomere length, as measured by quantitative PCR. Both parameters were compared with individual clinicopathological characteristics. Results: Clonal expansions and shorter telomeres were more frequent in nondysplastic biopsies from UC progressors than nonprogressors, but only for patients with early-onset of UC (diagnosis at younger than 50 years of age). Late-onset progressor patients had very few or no clonal expansions and longer telomeres. A few nonprogressors exhibited clonal expansions, which were associated with older age and shorter telomeres. In progressors, clonal expansions were associated with proximity to dysplasia. The mean percentage of clonally expanded mutations distinguished early-onset progressors from nonprogressors with 100% sensitivity and 80% specificity. Conclusions: Early-onset progressors develop cancer in a field of clonally expanded epithelium with shorter telomeres. The detection of these clones in a few random nondysplastic colon biopsies is a promising cancer biomarker in early-onset UC. Curiously, patients with late-onset UC seem to develop cancer without the involvement of such fields. (Inflamm Bowel Dis 2013;19:2593–2602)

Background: Kinase domain (KD) mutations in ABL1 are the dominant mechanism of relapse in patients (pts) with Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph+ ALL). Previous studies in Ph+ ALL have suggested that some pts harbor subclonal T315I mutations prior to tyrosine kinase inhibitor (TKI) treatment. However, nearly all reports have relied upon inherently error prone RT-PCR to generate template cDNA prior to mutation analysis. We hypothesized that conventional assays might over-estimate the incidence of pre-existing resistance mutations, and that improved sequencing accuracy might yield important information for risk stratification and TKI selection. Methods: Duplex sequencing (DS) is a molecular tagging method that improves the accuracy of conventional next-generation sequencing by more than 10,000-fold, by comparing the nucleotide sequences of each strand of double-stranded molecules. ABL1 DS was performed on genomic DNA from 64 pts with newly diagnosed Ph+ ALL treated with hyper-CVAD plus a TKI. DS of exons 4-10 of ABL1 was performed to an average molecular depth of >10,000x. Among pts who relapsed, using RNA extracted from relapse samples, the KD (codons 221 through 500) of the BCR-ABL1 fusion transcript was sequenced by the Sanger method using a nested PCR approach, with a detection limit of 10-20%. Results: The median age of the cohort was 54 years (range, 20-80 years). The TKI used was imatinib in 5 pts, dasatinib in 38 pts, and ponatinib in 21 pts. All pts achieved complete remission, and 12 pts (19%) underwent allogeneic stem cell transplantation in first remission. A total of 115 pretreatment ABL1 KD mutations were detected among 47 pts (73%). The median number of pretreatment ABL1 KD mutations was 2 (range, 0-6 mutations). The median variant allelic frequency (VAF) of the detected somatic mutations was 0.008% (range, 0.004%-0.649%). Five mutations (4%) and 40 mutations (35%) were present at a VAF ≥0.1% and ≥0.01%, respectively. Eleven mutations known to confer resistance to at least one TKI were detected in 7 pts (11%), and included: F317L in 4 pts, E225K in 2 pts, and E225V, L384M, M244V, Q252H and T315I in 1 pt each. Five mutations were detected in ≥1 pretreatment sample (F317L in 4 pts, M244V in 3 pts, and E255K, E459K and V355V in 2 pts each). Of these 7 pts with pretreatment resistance mutations, 5 pts received a TKI known not to be sensitive to the mutation(s); 2 pts who received ponatinib had mutations at least intermediately sensitive to ponatinib (1 pt with F317L and 1 pt with both E255K and M244V mutations). With a median duration of follow-up of 54 months (range, 1-124 months), 18 pts have relapsed. None of the 7 pts with known pretreatment resistance mutations relapsed. There was no difference in the number of pretreatment mutations between pts who relapsed and those who did not (median mutations: 1 [range, 0-4] and 2 [range, 0-6], respectively; P=0.26). Of the 18 pts who relapsed, 14 underwent Sanger sequencing for ABL1 KD mutations at the time of relapse. Clonal resistance mutations were detected at relapse in 9 pts (64% of sequenced samples): T315I in 6 pts, and F317I, V229L and V338G in 1 pt each (all in pts treated with dasatinib, except V338G in pt on imatinib). Relapse mutations were not observed in pretreatment samples in any of the pts. Conclusions: DS identified very low level pretreatment ABL1 KD mutations in a majority of pts with newly diagnosed Ph+ ALL but these appear inconsequential. Ninety percent of mutations identified have not been described as resistance mutations, suggesting that they may be either synonymous or functionally neutral amino acid changes resulting from normal aging. There was no apparent association of these mutations and risk of relapse, even in the minority of cases in which known resistance mutations were detected at baseline. These data suggest that pretreatment testing for ABL1 KD mutations in Ph+ ALL is unlikely to affect treatment decisions. Disclosures Short: Takeda Oncology: Consultancy. Sasaki:Otsuka Pharmaceutical: Honoraria. Ravandi:Jazz: Honoraria; Orsenix: Honoraria; Orsenix: Honoraria; Seattle Genetics: Research Funding; Sunesis: Honoraria; Xencor: Research Funding; Bristol-Myers Squibb: Research Funding; Amgen: Honoraria, Research Funding, Speakers Bureau; Abbvie: Research Funding; Astellas Pharmaceuticals: Consultancy, Honoraria; Seattle Genetics: Research Funding; Bristol-Myers Squibb: Research Funding; Astellas Pharmaceuticals: Consultancy, Honoraria; Amgen: Honoraria, Research Funding, Speakers Bureau; Macrogenix: Honoraria, Research Funding; Jazz: Honoraria; Abbvie: Research Funding; Xencor: Research Funding; Macrogenix: Honoraria, Research Funding; Sunesis: Honoraria. Cortes:Pfizer: Consultancy, Research Funding; Astellas Pharma: Consultancy, Research Funding; Novartis: Consultancy, Research Funding; Daiichi Sankyo: Consultancy, Research Funding; Arog: Research Funding. Konopleva:Stemline Therapeutics: Research Funding. Radich:TwinStrand Biosciences: Research Funding. Jabbour:Takeda: Consultancy, Research Funding; Bristol-Myers Squibb: Consultancy, Research Funding; Novartis: Research Funding; Pfizer: Consultancy, Research Funding; Abbvie: Research Funding.

Acute myeloid leukemia (AML) is a challenging disease to treat: most patients achieve remission after induction chemotherapy, but the majority eventually relapse. Minimal residual disease (MRD) after initial treatment is the best predictor of relapse and is thus a critical metric around which to develop new treatments. However, conventional MRD diagnostics, including cytology and flow cytometry, are of variable sensitivity and often only perform well in specialty centers, and there is no gold standard. Molecular tests developed to measure trace MRD in other hematological malignancies (i.e. CML) are high-resolution but assay a single, universally-present mutation, while many different genetic drivers exist in AML and these are spread among dozens of genes. Virtually every AML patient harbors a unique combination of mutations, making it difficult to design an effective universal assay. As such, most reported molecular AML MRD assays are either sensitive for mutations that are only found in a narrow subset of patients, or can screen many potential sites of mutation, but with low sensitivity. Here we present a broadly applicable Duplex Sequencing-based AML MRD assay that can readily detect mutant allele frequencies (MAF) below 1/10,000 across a large panel of genes, and below 1/100,000 in a focused panel. Conventional next generation sequencing (NGS) introduces errors during amplification and sequencing, creating a background of artifactual noise that obscures true mutations present below ~1%. Duplex Sequencing improves accuracy >100,000-fold through a molecular tagging approach whereby both strands of each original DNA duplex are ligated with a unique molecular barcode and amplified such that the reads generated from each strand can be related back to their unique original duplex. Reads can also be distinguished from those of their mate strand, thus the two strands of each DNA duplex can be compared and any discrepant nucleotide positions are discounted as errors. Our complete AML panel targets 151 exons or hotspot codons in 29 genes with a 59 kilobase (kb) hybrid-capture footprint. This region comprises loci containing single-nucleotide and short indel mutations found in approximately 90% of adult AMLs. A mean Duplex error-corrected sequence depth of 10,837 and a maximum Duplex depth of 14,967 was obtained across these targets from a single library preparation using 250 ng of sheared leukocyte DNA (Fig. 1). Duplex depth can be readily increased by preparing additional Duplex libraries from the same source DNA to achieve proportionally higher sensitivity for rarer variants. This stands in contrast to conventional NGS where, beyond a modest level, an increase in depth simply increases the number of background errors identified (Fig. 2A). We simulated low-level residual disease by mixing control DNA from a healthy young blood donor with DNA from 9 human cell lines harboring known AML mutations at dilutions from 1:100 to 1:100,000 (Table 2). The genomic loci of these 9 mutations in NRAS, KRAS and TP53 were captured with a small 1 kb probe panel. This mixture was sequenced to a mean Duplex depth of >1,000,000-fold, with the highest and lowest MAFs shown in Fig. 2B. All were close to expected frequencies (r2=0.96) with MAF as low as 6x10-6 (Fig. 3). As proof of specificity, we examined all coding nucleotide positions (excluding the 9 expected variants) and identified only 241 background variant counts out of 414,452,402 total Duplex BP, for an aggregate mutation frequency of 5.8x10-7, consistent with the estimated background of normal human aging. Our Duplex Sequencing-based AML MRD assay is flexible, broadly applicable and extremely sensitive. The assay is easily implemented using standard NGS equipment and automated cloud-based analysis software. The ~90% of AML patients served by this SNV-focused panel can be expanded to nearly 100% with complementary indel detection via targeted NGS RT-PCR. Optionally, when a patient's mutation profile from time-of-diagnosis is known, MRD testing can focus exclusively on those targets using a subset of pre-validated probes to reduce sequencing cost. Improved MRD testing will facilitate accurate prognostication, better selection among treatment options, and could serve as a surrogate endpoint in clinical trials to bring new treatments to patients faster. We are currently evaluating Duplex Sequencing MRD tests in both retrospective and prospective clinical trials. Disclosures Higgins: TwinStrand Biosciences: Employment. Williams:TwinStrand Biosciences: Employment. Buckley:CTI Biopharma: Employment; TwinStrand Biosciences: Consultancy. Radich:TwinStrand Biosciences: Research Funding. Salk:TwinStrand Biosciences: Employment, Equity Ownership.

Cancer recapitulates Darwinian evolution. Mutations acquired during life that provide cells with a growth or survival advantage will preferentially multiply to form a tumor. As a result of The Cancer Genome Atlas Project, we have gathered detailed information on the nucleotide sequence changes in a number of human cancers. The sources of mutations in cancer are diverse, and the complexity of those found to be clonally present in tumors has increasingly made it difficult to identify key rate-limiting genes for tumor growth that could serve as potential targets for directed therapies. The impact of DNA sequencing on future cancer research and personalized therapy is likely to be profound and merits critical evaluation.

Chronic inflammation predisposes to a variety of human cancers. Affected tissues slowly accumulate mutations, some of which affect growth regulation and drive successive waves of clonal evolution, whereas a far greater number are functionally neutral and serve only to passively mark expanding clones. Ulcerative colitis (UC) is an inflammatory bowel disease, in which up to 10% of patients eventually develop colon cancer. Here we have mapped mutations in hypermutable intergenic and intronic polyguanine tracts in patients with UC to delineate the extent of clonal expansions associated with carcinogenesis. We genotyped colon biopsies for length altering mutations at 28 different polyguanine markers. In eight patients without neoplasia, we detected only two mutations in a single individual from among 37 total biopsies. In contrast, for 11 UC patients with neoplasia elsewhere in the colon, we identified 63 mutations in 51 nondysplastic biopsies, and every patient possessed at least one mutant clone. A subset of clones were large and extended over many square centimeters of colon. Of these, some occurred as isolated populations in nondysplastic tissue, considerably distant from neoplastic lesions. Other large clones included regions of cancer, suggesting that the tumor arose within a preexisting clonal field. Our results demonstrate that neutral mutations in polyguanine tracts serve as a unique tool for identifying fields of clonal expansions, which may prove clinically useful for distinguishing a subset of UC patients who are at risk for developing cancer.

Clinical oncology is being revolutionized by the increasing use of molecularly targeted therapies. This paradigm holds great promise for improving cancer treatment; however, allocating specific therapies to the patients who are most likely to derive a durable benefit continues to represent a considerable challenge. It is becoming increasingly clear that cancers are characterized by extensive intratumour genetic heterogeneity, and that patients being considered for treatment with a targeted agent might, therefore, already possess resistance to the drug in a minority of cells. Indeed, multiple examples of pre-existing subclonal resistance mutations to various molecularly targeted agents have been described, which we review herein. Early detection of pre-existing or emerging drug resistance could enable more personalized use of targeted cancer therapy, as patients could be stratified to receive the therapies that are most likely to be effective. We consider how monitoring of drug resistance could be incorporated into clinical practice to optimize the use of targeted therapies in individual patients.

High-throughput DNA sequencing has fundamentally changed practice standards in clinical oncology. We can now systematically examine the dysfunctional genes and pathways of individual tumors and can often pair these with targeted therapies. Moreover, NGS has improved our ability to resolve molecular heterogeneity, such as results from evolving subclones within a tumor or admixtures of cancer-derived and normal DNA in liquid biopsies. However, this remains imperfect. With standard NGS, sequencing errors create a background that obscures variants present at less than 1-2%. Methods of tag-based error correction improve accuracy, but most rely on hemimolecular tagging and are thus still limited by certain types of errors caused by DNA damage. Our group pioneered Duplex Sequencing, an error correction technology that independently genotypes and compares both single-stranded halves of individual DNA duplexes. This approach improves sequencing accuracy by >100,000-fold and allows detection of a single mutation among ten million unmutated nucleotides. We carried out Duplex Sequencing on uterine lavages from women with and without high-grade serous ovarian cancer to determine the sensitivity for identifying tumor-derived TP53 mutations as a screening test. Under blinded conditions we detected tumor-derived mutations in 80% of cancer cases. However, in addition, we found that nearly all women, including those without cancer, carried very low frequency TP53 mutations that did not originate from a tumor. These could not be explained by technical errors, as they exhibited a strong signature of positive selection, indicated by higher-than-random non-synonymous fraction, substantial clustering in cancer-associated hotspots and a spectrum, exonic distribution and predicted pathogenicity closely mirroring TP53 mutations in cancer databases. Furthermore, abundance of these “biological background” mutations increased with age. To more broadly investigate the extent of clonal evolution in normal aging, we used Duplex Sequencing to examine tissue specimens from two individuals without cancer, differing in age by 101 years. We assessed relative mutation load in 22 cancer-associated genes and one non-cancer associated gene in 19 different tissues. Initial results indicate that the extent of age-associated mutation accumulation varies by organ and the pattern of apparent clonal selection is influenced by both tissue type and genetic context. Our findings illustrate the near-universal presence of cancer-like mutations arising from age-associated processes in multiple non-neoplastic tissues. Somatic clonal evolution appears to be a ubiquitous part of human life. Beyond many biological implications, this highlights an important specificity challenge that must be accounted for when using ultra-accurate sequencing technologies for cancer screening and other sensitive diagnostics. Citation Format: Jesse J. Salk, Elisabeth Maritschnegg, Paul Speiser, Robert Zeillinger, Rosana Risques, Lawrence A. Loeb. Duplex Sequencing detects cancer-associated mutations arising during normal aging: Clonal evolution over a century of human lifetime [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3041. doi:10.1158/1538-7445.AM2017-3041

Duplex Sequencing (DS) is a next-generation sequencing methodology capable of detecting a single mutation among >1 × 10(7) wild-type nucleotides, thereby enabling the study of heterogeneous populations and very-low-frequency genetic alterations. DS can be applied to any double-stranded DNA sample, but it is ideal for small genomic regions of

Intratumor genetic heterogeneity reflects the evolutionary history of a cancer and is thought to influence treatment outcomes. Here we report that a simple PCR-based assay interrogating somatic variation in hypermutable polyguanine (poly-G) repeats can provide a rapid and reliable assessment of mitotic history and clonal architecture in human cancer. We use poly-G repeat genotyping to study the evolution of colon carcinoma. In a cohort of 22 patients, we detect poly-G variants in 91% of tumors. Patient age is positively correlated with somatic mutation frequency, suggesting that some poly-G variants accumulate before the onset of carcinogenesis during normal division in colonic stem cells. Poorly differentiated tumors have fewer mutations than well-differentiated tumors, possibly indicating a shorter mitotic history of the founder cell in these cancers. We generate poly-G mutation profiles of spatially separated samples from primary carcinomas and matched metastases to build well-supported phylogenetic trees that illuminate individual patients' path of metastatic progression. Our results show varying degrees of intratumor heterogeneity among patients. Finally, we show that poly-G mutations can be found in other cancers than colon carcinoma. Our approach can generate reliable maps of intratumor heterogeneity in large numbers of patients with minimal time and cost expenditure.

Lou et al. (1) report a technique for increased sensitivity of DNA sequencing that they call “circle sequencing.” The authors’ method involves circularizing single-stranded DNA and performing rolling-circle amplification (RCA). Multiple duplicates are thereby generated, and errors can be removed by comparing the sequences of the duplicates. Lou et al. report accuracy similar to that of single-stranded tagging techniques (2, 3), but less than that achieved with double-stranded tagging (4). The authors propose their method as a high-efficiency tool for deep-sequencing of heterogeneous DNA samples.

The extent of heterogeneity in DNA sequence within human cancers has been difficult to fully assess, due to the absence of methods with sufficient sensitivity to detect minority variants. For example, next-generation sequencing could in principle detect sub-clonal mutations at any genomic position; however due to amplification errors and inaccuracies in the sequencing platform itself, such approaches are limited to detection of mutations present in >1% of cells. We have overcome this limitation by developing an approach, Duplex Sequencing, which improves the accuracy of next-generation sequencing by >100,000 fold. Duplex Sequencing is based upon separately tagging and sequencing the two strands of single molecules of duplex DNA. True mutations are present at the same position in both DNA strands and are complementary, whereas artifacts arising from amplification or sequencing errors are seen in only one strand. The calculated error rate of our approach is less than one artifactual mutation per billion nucleotides sequenced (ref. 1). We have also developed an enrichment approach based on sequential rounds of hybridization to biotinylated probes which enables efficient sequencing of targeted regions of the genome (ref. 2). Duplex Sequencing thereby enables exceptionally sensitive detection of sequence heterogeneity within any set of genes. We have now applied Duplex Sequencing to study of sequence heterogeneity in both normal and malignant samples. In acute myeloid leukemia (AML), we find multiple sub-clonal mutations which fall below the detection limit of conventional approaches, some of which encode mutations that can drive chemotherapy resistance and cancer progression. Our results indicate that prior studies have under-estimated the burden of sub-clonal mutations in AML by more than 1,000-fold. We have additionally studied chronic myeloid leukemia (CML); treatment of CML with inhibitors directed against the Abl kinase represents the prototypical targeted cancer therapy. With Duplex Sequencing, we find that Abl mutations are extremely uncommon at the time of CML diagnosis. In contrast, refractory patients who fail therapy frequently possess multiple sub-clonal mutations within the Abl kinase, most of which are below the resolution of conventional approaches but are readily detected with our method. In ongoing work, we have found multiple sub-clonal mutations within human prostate cancer and colon cancer, suggestive of widespread intratumor heterogeneity that may limit the efficacy of single-agent therapy in these diseases. References: 1. Schmitt MW, Kennedy SR, Salk JJ, Fox EJ, Hiatt JB, Loeb LA (2012). Detection of ultra-rare mutations by next-generation sequencing. Proceedings of the National Academy of Sciences USA. 109(36):14508-13. 2. Schmitt MW, Fox EJ, Prindle MJ, Reid-Bayliss KS, True LD, Radich JP, Loeb LA (2015). Sequencing small genomic targets with high efficiency and extreme accuracy. Nature Methods. 12(5):423-5. Citation Format: Michael W. Schmitt, Jesse J. Salk, Edward J. Fox, Justin R. Pritchard, J Graeme Hodgson, Victor M. Rivera, Pamela S. Becker, Jerald P. Radich, Lawrence A. Loeb. Direct assessment of sequence heterogeneity in human cancers by Duplex Sequencing. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 2395.

Next-generation DNA sequencing promises to revolutionize clinical medicine and basic research. However, while this technology has the capacity to generate hundreds of billions of nucleotides of DNA sequence in a single experiment, the error rate of ∼1% results in hundreds of millions of sequencing mistakes. These scattered errors can be tolerated in some applications but become extremely problematic when “deep sequencing” genetically heterogeneous mixtures, such as tumors or mixed microbial populations. To overcome limitations in sequencing accuracy, we have developed a method termed Duplex Sequencing. This approach greatly reduces errors by independently tagging and sequencing each of the two strands of a DNA duplex. As the two strands are complementary, true mutations are found at the same position in both strands. In contrast, PCR or sequencing errors result in mutations in only one strand and can thus be discounted as technical error. We determine that Duplex Sequencing has a theoretical background error rate of less than one artifactual mutation per billion nucleotides sequenced. In addition, we establish that detection of mutations present in only one of the two strands of duplex DNA can be used to identify sites of DNA damage. We apply the method to directly assess the frequency and pattern of random mutations in mitochondrial DNA from human cells.

Genome instability is regarded as a hallmark of cancer. Human tumors frequently carry clonally expanded mutations in their mitochondrial DNA (mtDNA), some of which may drive cancer progression and metastasis. The high prevalence of clonal mutations in tumor mtDNA has commonly led to the assumption that the mitochondrial genome in cancer is genetically unstable, yet this hypothesis has not been experimentally tested. In this study, we directly measured the frequency of non-clonal (random) de novo single base substitutions in the mtDNA of human colorectal cancers. Remarkably, tumor tissue exhibited a decreased prevalence of these mutations relative to adjacent non-tumor tissue. The difference in mutation burden was attributable to a reduction in C∶G to T∶A transitions, which are associated with oxidative damage. We demonstrate that the lower random mutation frequency in tumor tissue was also coupled with a shift in glucose metabolism from oxidative phosphorylation to anaerobic glycolysis, as compared to non-neoplastic colon. Together these findings raise the intriguing possibility that fidelity of mitochondrial genome is, in fact, increased in cancer as a result of a decrease in reactive oxygen species-mediated mtDNA damage., Author Summary Mitochondria are the intracellular organelles responsible for energy production in eukaryotic cells. They are unique in that they contain their own DNA (mtDNA), which encodes genes important for mitochondrial function and is the cell's only genetic material stored outside the nucleus. Mutations in both nuclear and mtDNA are believed to have a role in tumor growth and metastasis. While it is well established that nuclear DNA has an increased overall burden of mutations in human cancers, this has not been investigated in mtDNA. Here we use a highly sensitive assay to determine mutational load in mtDNA isolated from patient-matched normal and cancerous colonic tissues. Surprisingly, we discovered that the frequency of mutations in mtDNA is actually decreased in tumors relative to healthy tissues, suggesting that, unlike in nuclear DNA, accelerated mutagenesis in mtDNA does not facilitate a cancer's development and may even hinder it. These findings raise the possibility that mitochondrial-targeted cancer therapeutics focused on directly increasing mtDNA damage and mutation might effectively suppress malignant growth.