- Method
- Open access
- Published:
Oncogene Concatenated Enriched Amplicon Nanopore Sequencing for rapid, accurate, and affordable somatic mutation detection
Genome Biology volume 22, Article number: 227 (2021)
Abstract
We develop the Oncogene Concatenated Enriched Amplicon Nanopore Sequencing (OCEANS) method, in which variants with low variant allele frequency (VAFs) are amplified and subsequently concatenated for Nanopore Sequencing. OCEANS allows accurate detection of somatic mutations with VAF limits of detection between 0.05 and 1%. We construct 4 distinct multi-gene OCEANS panels targeting recurrent mutations in acute myeloid leukemia, melanoma, non-small- cell lung cancer, and hepatocellular carcinoma and validate them on clinical samples. By demonstrating detection of low VAF single nucleotide variant mutations using Nanopore Sequencing, OCEANS is poised to enable same-day clinical sequencing panels.
Background
High-throughput DNA sequencing is becoming a standard part of oncology care, with many laboratory-developed tests informing patient prognosis [1], therapy selection[2, 3], and minimal residual disease [4]. DNA sequencing is furthermore being explored as a method for enabling early screening of cancers in asymptomatic populations [5, 6]. The dominant platforms used for clinical high-throughput sequencing today are based on sequencing-by-synthesis (next-generation sequencing (NGS), e.g., the Illumina and Ion Torrent platforms).
Although NGS has many advantages including high throughput, high sensitivity, and high reproducibility/reliability, NGS also has three notable limitations: First, NGS read lengths are limited to roughly 300 nt, rendering it less suitable/sensitive to larger scale DNA alterations such as copy number variations and chromosomal translocations. Second, NGS is time-consuming, with multiple days needed for sequencing-by-synthesis chemistry, in addition to time- and labor-intensive library preparation and bioinformatic interpretation. Third, NGS requires a high capital investment on the order of <DOLLAR/>1,000,000 for a high-throughput instrument (e.g., Novaseq); more affordable NGS platforms such as the MiniSeq result in roughly 10-fold higher per-read sequencing costs.
Nanopore Sequencing overcomes all three of the above limitations of NGS, with reads as long as the fragments that are loaded, sequencing times as short as 15 min, and Oxford Nanopore MinION instrument plus starter pack costing approximately <DOLLAR/>1000 and size of a USB drive. On the other hand, while Nanopore Sequencing has historically had a higher error rate than NGS, there has been intense recent research effort in reducing Nanopore Sequencing error rates on both the hardware [7] and software [8] sides, which has brought the raw read error rate up to a competitive level [7, 9]. Detection of somatic single-base mutations at low variant allele frequencies (VAF) has not been demonstrated in peer reviewed publications [10–12]. New technologies are starting to be developed to enable low VAF detection using Nanopore Sequencing [35]. However, rapid mutation profiling with turn around times less than 1 day has not been demonstrated.
The importance of structural variations in cancer are beginning to be understood using long-read technologies. However, a large fraction, if not a majority, of actionable cancer DNA alterations reported to date were done so using short-read technologies and are therefore are single-base mutations [13], and in tumor tissue samples, there may be significant cancer heterogeneity and/or low tumor fraction. Consequently, reliable detection of single-base mutations at 5% VAF in formalin-fixed, paraffin-embedded (FFPE) tissues is currently considered a standard requirement for clinical NGS assays [14, 28]. Cell-free DNA in peripheral blood plasma is an emerging biospecimen for noninvasive cancer monitoring [15] and has even more stringent requirements on VAF limit of detection (LoD), with typical commercial NGS kits and services claiming LoDs of between 0.1 and 0.5% VAF [29, 30].
Here, we present a new method, Oncogene Concatenated Enriched Amplicon Nanopore Sequencing (OCEANS), that achieves ≤1% VAF LoD on single-base mutations in FFPE samples by Nanopore Sequencing (Fig. 1). First, we describe a amplicon concatenation method called Stochastic Amplicon Ligation (SAL) to utilize the long read capabilities of Nanopore Sequencing for sequencing short amplicons. SAL approach would be suitable for any amplicon-based short fragment library from samples like FFPE and cfDNA. The OCEANS method integrates the blocker displacement amplification (BDA) allele enrichment method [16, 17] with SAL, improving the VAF LoD by roughly 100-fold and throughput by roughly 10-fold using SQK-LSK109 sample prep chemistry on R9.4.1 flow cells. The entire OCEANS method takes less than 10 h from DNA to called variants, and the average sequencing cost per FFPE sample is roughly <DOLLAR/>7.5 for a 7-gene, 15-amplicon panel on the Oxford Nanopore MinION flow cells.
Results
Stochastic Amplicon Ligation. DNA samples for oncology sequencing are typically extracted from FFPE tissues and can have average lengths of less than 500 nt due to accumulated chemical damage [18]. We developed the Stochastic Amplicon Ligation (SAL) method to enzymatically concatenate many short DNA molecules together to utilize the long-read capability of Nanopore Sequencing and increase the effective throughput.
SAL is based on the Golden Gate assembly method used in synthetic biology to concatenate short oligos into synthetic genes [19]. In SAL (Fig. 2 b), amplicons are appended with engineered adapter sequences that possess a Type IIS restriction enzyme recognition site. After Type IIS cleavage, a 4 nt sticky end is left on the 5 ′ ends of both strands of the amplicons; these sticky ends allow the amplicons to transiently bind to each other, which then enzymatically ligate to form concatemers. Multiple cycles of enzymatic restriction and ligation are performed to increase the lengths of the concatemers, and we perform a SPRI (solid phase reversible immobilization) size selection afterwards to both enrich long concatemers and remove short recognition site oligos cleaved from the amplicons. The multiple temperature cycles between 37 and 16 ∘C improve the mean lengths of the concatemer assemblies by keeping the concentrations of activated monomer amplicons low. Prior literature [19] suggests that direct ligation of amplicons with 5 ′ sticky ends results in a much larger population of shorter concatemers.
SAL differs from traditional Golden Gate assembly in having universal sticky end sequences to allow stochastic incorporation of any amplicon with the appropriate adapters. This, in principle, allows unlimited growth of concatemers to longer lengths given sufficient monomer concentration and enough temperature cycles and reduces the possibility of unintended reactions due to nonspecific binding between non-cognate sticky ends. Experimentally, capillary electrophoresis indicated that SAL concatenated a mean of roughly 12 to 15 monomers per concatemer (Fig. 2 c). In addition to improving the throughput of Nanopore Sequencing, we also found that the SAL improved the quality of Nanopore Sequencing reads (Fig. 2 d). The Nanopore Sequencing results of a 340-nt amplicon had a mean phred quality score of 9.87, corresponding to an error rate of 10.3%. The concatemer, in contrast, had a mean phred score of 11.55, corresponding to an error rate of 7.0%. The lower quality score of shorter reads is due to lack of sufficient current signal information for proper normalization prior to basecalling by MinKNOW (Personal communication by ONT technical support).
The improved throughput and accuracy of Nanopore Sequencing of SAL concatemers allow calling somatic single-base mutations at 5% VAF when matched normal sample are available (Fig. 2 e). We first applied Nanopore Sequencing to amplicons from a 95%:5% mixture of the NA18537 and NA18562 human cell line genomic DNA. NA18537 and NA18562 were homozygous for different alleles at the rs3789806 and rs9648696 single nucleotide polymorphism (SNP) loci, so the mixture was 5% VAF in the NA18537 SNP alleles. In our analysis, it was difficult to call somatic mutations at 5% VAF, since there was a large number of loci on the amplicon with high variant read frequencies. When we subtracted the variant read fraction (VRF) from a 100% NA18537 sample, then the 5% VAF somatic mutations became more visible in the ΔVRF. In direct amplicon nanopore sequencing, the two 5% VAF SNP alleles were detected in the ΔVRF figure with +6.5 σ and +4.4 σ, respectively. The confidence of calling these 5% VAF variants were increased to +17.4 σ and +13.3 σ for SAL concatemers. For SAL concatemers, the long nanopore sequencing reads were bioinformatically deconcatenated using a custom python code [33].
The lower throughput and quality score of shorter DNA libraries are specific to SQK-LSK109 sample prep chemistry on R9.4.1 flow cells using MinKNOW 19.12.5 for basecalling. Sequencing throughput and quality of short libraries could be on par with long DNA libraries using ONT’s latest and upcoming improvements to sequencing chemistry and basecalling algorithms [36].
Rolling circle amplification (RCA) is an alternative method for generating long DNA from shorter DNA molecules, which has been used in the context of Nanopore Sequencing for improving accuracy [20]. RCA has several significant limitations compared to SAL, most notably that RCA requires an initial circularization of DNA which is known to have low efficiency, reducing the clinical sensitivity due to low conversion yield of biological DNA molecules in sequencing library. Additionally, RCA produces concatemers in which the segment sequences all reflect the sequence of the original molecule, rather than a uniform sampling of all DNA molecules on the loci of interest. Finally, RCA generates single-stranded DNA products rather than double-stranded DNA products, which should be converted into double-stranded DNA for efficient Nanopore Sequencing.
Integrating BDA allele enrichment. Frequently, matched normal FFPE tissue samples will not be available, so using SAL alone for Nanopore Sequencing detection of low VAF somatic mutations is unlikely to be impactful clinically. The OCEANS method employs blocker displacement amplification (BDA) [16, 17] to allow more robust detection of low VAF somatic mutations without requiring a matched normal sample. In brief, BDA includes a wildtype-binding blocker oligonucleotide that competes with a PCR primer in hybridizing to DNA templates of interest. The blocker binds more strongly than the primer to wildtype DNA sequences, preventing efficient PCR amplification. On DNA templates with sequence variants, the primer outcompetes the blocker and PCR proceeds as usual. Through the course of many PCR cycles (20–25 cycles), the VAF of sequence variants (including single nucleotide mutations) can be enriched by over 1000-fold.
In OCEANS, the DNA biospecimen is first mixed with multiple primers and blockers to undergo variant-selective PCR amplification (Fig. 3 a). The amplicons will over-represent sequence variants in genetic loci of interest, though some amplicons with wildtype sequences will still exist. The amplicons are subsequently appended with SAL adapters and concatenated into concatemers, and size-selected to remove short assemblies, primers, etc. The concatemers are then ligated to the standard Oxford Nanopore Sequencing adapters with attached motor proteins and loaded into the nanopore sequencing flow cell. The entire workflow takes about 10 h, including post-sequencing bioinformatics. On the same SNP alleles as in Fig. 2, the OCEANS results showed VRFs that were dramatically higher than the sample variant VAFs (Fig. 3 b), with 0.1% VAF SNP allele enriched to over 70% VRF. Thus, OCEANS allows robust variant calls of somatic mutations without the need for a matched normal DNA sample, which was not possible previously on the Nanopore Sequencing platform.
We next constructed two multiplexed OCEANS panels: a 7-amplicon panel covering recurrent mutations observed in acute myeloid leukemia (AML) and a 15-amplicon panel covering recurrent mutations observed in melanoma. The AML panel covers roughly 254 mutations in the COSMIC database across 7 genes (Fig. 4 a), and the melanoma panel covers roughly 370 mutations across 8 genes (Fig. 4 d). We first characterized the limit of detection for mutations covered by these OCEANS panels using synthetic spike-in reference samples, with VAFs ranging from 0.05 and 1%.
Variant calls were made using two different approaches: (1) based on the variant read frequency exceeding a threshold of 20% and (2) based on a Clair [21] score of above 180 (Personal communication from ONT). We found that both approaches were imperfect: considering VRF alone ignores the fact that Nanopore Sequencing has different error rates for certain sequences, e.g., homopolymers. Clair scores, on the other hand, are not monotonic with VAF and have been observed to be less accurate for indel calls [10, 25–27]. To ensure minimal false positives in variant calls, we require that a variant must be independently called by both methods in order to be reported. On our internal reference samples (Fig. 4 a, b, d), we observed VAF limits of detection between 0.05 and 1%. We did not observe any effect on enrichment of variants by excluding SAL from OCEANS workflow (Additional file 1: Figure S11). But SAL significantly increases the throughput for short amplicons by harnessing the long read capabilities of Nanopore Sequencing. We next applied our OCEANS panels to third-party reference samples from Horizon Discovery, with mutation VAFs at 5% (for AML) and 0.5% (for melanoma). The expected mutations were all called.
Interestingly, we also made a number of unexpected variant calls in the melanoma OCEANS panel (Fig. 4 e). The variant in the MAP2K2-207 amplicon was confirmed to be a non-pathogenic SNP (rs10250). The variants in the PIK3CA-542 amplicon were found to be aligned to the PIK3CA pseudogene (LOC100422375) [22], which were also preferentially enriched by BDA. The “variants” associated with the PIK3CA pseudogene were bioinformatically excluded from variant calls in subsequent experiments. Finally, we made confident variant calls for PIK3CA c.3140A >G and KRAS c.38G >A. We contacted Horizon Discovery customer support regarding these putative mutations, and the latter confirmed that these mutations are also present at low VAFs in the HD238 sample.
Validating OCEANS on clinical tissue samples. We next applied the melanoma OCEANS panel to clinical melanoma tissue samples, including both fresh/frozen (FF) and FFPE tissue (Fig. 5 a, b). As in the calibration experiments, we called somatic mutations only when the VRF was observed to be greater than 20%, and the Clair score was above 180. In total, DNA from 7 FF and 18 FFPE tissue samples were sequenced using both OCEANS and NGS. The melanoma OCEANS panels cover a total of 384 loci, corresponding to a total of 9600 total loci analyzed across the 25 samples.
Figure 5 c shows the comparison between OCEANS and NGS. All 16 somatic mutants called by NGS at above 5% VAF were also called by OCEANS, corresponding to a 100% OCEANS sensitivity relative to NGS. Of the 9584 NGS-negative loci, OCEANs called an additional 97 variants (Fig. 5 c); thus, relative to NGS, OCEANS had a 99.0% specificity.
We calculated the original sample VAF from OCEANS VRF using fold enrichment calculated for mutations in calibration experiments (Additional file 1: Section S5) [17]. The sequencing error rates combined with the saturation of VRFs near 100% after BDA enrichment means that our quantitation dynamic range is relatively small. However, estimation of sample VAF enables identification of high VAF (>5%) mutations to aid in making treatment decisions based on clinical diagnosis. OCEANS identified all mutations with NGS VRF >5% as high VAF mutations (Fig. 5 d). The OCEANS VRF and NGS VRF cannot be directly compared, since OCEANS VRF is not the true sample VAF. Therefore, precision-recall values were calculated from OCEANS calculated VAF and NGS VRF comparison in Fig. 5 d. By varying the OCEANS calculated VAF cutoff threshold, we can change the number of high VAF mutation calls that are verified by NGS as >5%, generating a set of precision/recall tradeoffs for detecting mutations with NGS VRF >5%, which can be plotted as a precision-recall curve (Fig. 5 e). Importantly, we believe that many of the 97 discordant called variants that were below the 5% NGS VRF cutoff could be real mutations, based on our calibration experiments. To confirm our discordant OCEANS mutation calls, we further performed digital PCR on 6 FFPE samples at 4 mutation loci (BRAF p. V600, KRAS p. G13D, KRAS p. E62K, and MAP2K1 p. P124L) and one fresh frozen sample for BRAF p. V600 loci (Supplementary excel table). Of these 25, 12 mutations were called positive by OCEANS and 13 were called negative by OCEANS. OCEANS was concordant with ddPCR for 11 positive samples and 11 negative samples (Table 1, Additional file 1: Section S6). Identification of such low VAF mutations makes OCEANS suitable for applications like Minimum Residual Disease (MRD) detection.
It is important to note that concordant positives between OCEANS and ddPCR indicate the existence of a DNA variant in the sample, which may not necessarily reflect a mutation in the patient. Cytosine deamination is a well-documented type of DNA damage frequently observed in DNA extracted from FFPE. We applied an FFPE damage repair kit to the FFPE DNA before performing OCEANS library preparation, but do not necessarily expect that all cytosine deaminations are repaired or excised. In particular, any repair kit based on cleaving/repairing uracils formed through the deamination of standard cytosine would not be able to rectify deamination of methylcytosines into thymines.
Because each ddPCR mutation requires a separate reaction, the ddPCR results required 4 times more input DNA than OCEANS just to cover these 4 mutations. For analysis of clinical biopsy samples, OCEANS would have significantly higher clinical sensitivity due to being able to analyze all mutations in the panel from a single sample.
Next, we wished to characterize the reproducibility and robustness of the OCEANS panel on different types of Nanopore Sequencing instruments and flow cells. The Oxford Nanopore Flongle flow cell, in particular, is relatively inexpensive at <DOLLAR/>70 and can further reduce turnaround time relative to MinION by reducing the need for sample batching before sequencing. We performed the OCEANS panel on all 25 melanoma samples on the Flongle, and observed highly quantitatively similar VRFs as our results on the MinION (Fig. 5 f).
NSCLC and HCC OCEANS panels. We next constructed two additional OCEANS panels: a 28-amplicon panel for non-small cell lung cancer (NSCLC) and an 11-amplicon panel for hepatocarcinoma (HCC) to show the generality of our approach. The NSCLC OCEANS panel covers roughly 1121 mutations in the COSMIC database across 13 genes (AKT1, ALK, BRAF, DDR2, EGFR, KRAS, NRAS, MAP2K1, MET, PIK3CA, PTEN, ROS1, and TP53, see Additional file 1: Section S4). DNA from 5 FF and 18 FFPE NSCLC tissue samples were sequenced using both OCEANS and NGS. Figure 6 a, c show the comparison between OCEANS and NGS. Nine out of 11 somatic mutants called by NGS at above 5% VAF were also called by OCEANS. The two mutations that had a Clair score less than 180 were indel mutations, for which Clair has been observed to be less accurate [10, 25].
The HCC OCEANS panel covers roughly 680 mutations across 7 genes (CTNNB1, ARID1A, AXIN, TERT, JAK1, PTEN, and TP53, see Additional file 1: Section S4). DNA from 5 FF and 16 FFPE HCC tissue samples were sequenced using both OCEANS and NGS. Figure 6 b, d show the comparison between OCEANS and NGS. Fourteen out of 17 somatic mutants called by NGS at above 5% VAF were also called by OCEANS. The 3 mutations not called by Clair were in the TERT amplicon within a homopolymer region (Additional filer~efMOESM1: Section S5). Higher Nanopore Sequencing error rates in homopolymer regions could be the reason for lower Clair scores despite the OCEANS VRF being >70% for these mutations. We observed 23 loci with OCEANS VRF greater than or equal to 20% for which the corresponding NGS VRF were 0%. We analyzed the NGS read depth for all loci with NGS VRF equal to 0% and their corresponding OCEANS VRF (Fig. 6 g). All 23 loci with OCEANS VRF greater than or equal to 20% had NGS read depth of less than 300. Overall, both OCEANS panels had high concordance between OCEANS and NGS. The area under the precision-recall curve was 92.34% for the NSCLC panel and 90.74% for the HCC panel based on OCEANS calculate VAF (Fig. 6 e, f).
Discussion
The rapid turnaround time and low instrument/consumables cost of Nanopore Sequencing render Nanopore Sequencing attractive for time-sensitive therapy selection and recurrence monitoring applications in oncology. Although multi-nucleotide indels, copy number variations, and gene fusions can be detected at low VAFs by Nanopore Sequencing, previously the rapid detection of single nucleotide somatic mutations with low VAFs has not been demonstrated on the Nanopore Sequencing platform. The OCEANS method we present detects single-base mutations from FFPE tissue-derived DNA with LoD of less than 1% with a 10-h workflow and positions Nanopore Sequencing to transition into clinical sequencing panels (Table 2). Within OCEANS, BDA enhances the confidence and sensitivity of somatic mutation calls, and SAL reduces the turnaround time by enhancing the throughput of Nanopore Sequencing. However, the latest and upcoming improvements to ONT’s sequencing chemistry and basecalling algorithms [36] may remove the need for SAL to get higher throughput with short DNA libraries, enabling BDA to be applied directly for Nanopore Sequencing.
For liquid tumors and for solid tumors where fresh/frozen tissue biopsy samples are available, the long read length of Nanopore Sequencing further enables accurate detection of large-scale DNA structural variations. Integrating OCEANS with existing Nanopore Sequencing experimental and bioinformatics methods for structural variant profiling would allow for the development of comprehensive oncology panels covering a broad range of structural alterations observed in leukemias and lymphomas. In contrast, the DNA from FFPE tissue samples and cell-free DNA from plasma are short and do not physically allow for direct long-read sequencing to identify large-scale alterations. We used Clair variant caller that is based on DeepVariant variant calling methodology, since it has been shown to have the highest accuracy for calling mutations on Nanopore Sequencing data [21, 32]. Improvements in variant calling software for Nanopore Sequencing will improve the specificity and sensitivity of OCEANS especially for in-del mutations. Here, we used SQK-LSK109 sample prep chemistry on R9.4.1 flow cells and MinKNOW 19.12.5 (fast basecalling model) for basecalling. The latest R10.3 flow cells in combination with ONT’s upcoming Q20+ sequencing chemistry with improved accuracy and throughput may further enhance OCEANS performance and turnaround time for low VAF detection.
In this work, our bioinformatics pipelines make qualitative variant calls based on a combination of observed VRF and Clair score. We used the fold-enrichment calculated for each BDA design and inferred the sample VAF from VRF for OCEANS (Additional file 1: Section S4) [17]. This enabled us to identify high VAF (≥5%) and low VAF mutations. However, the sequencing error rates combined with the saturation of VRFs after BDA enrichment mean that our quantitation dynamic range is relatively small. We underestimated the sample VAF if the VRF saturates or if the sequencing error rate at that particular mutation position was high. We overestimated the sample VAF for mutations like deletions and insertions of multiple nucleotides that resulted in a higher enrichment fold than that of the single nucleotide mutations used in calibration experiments. Integration of OCEANS with unique molecular identifier (UMI) barcodes [23, 24] will likely be necessary for accurate VAF quantitation on the Nanopore Sequencing platform.
Conclusions
We developed a method (OCEANS) to selectively amplify variants with low VAFs and subsequently concatenate the amplicons for Nanopore Sequencing. OCEANS allows accurate detection of somatic mutations with VAF limits of detection between 0.05 and 1% by Nanopore Sequencing. On clinical FFPE tumor samples OCEANS showed 99.79 to 99.99% area under the receiver-operator curve in comparison experiments against Illumina NGS. Comparison against digital PCR on 11 putative mutations at ≤1% VAF showed 10 concordant positive calls with VAFs between 0.02 and 0.66%. The rapid turnaround time and low instrument/consumables cost of Nanopore Sequencing is designed to enable same-day clinical sequencing panels. Integrating OCEANS with long-read and base modification detection capabilities of Nanopore Sequencing can enable development of comprehensive oncology panels.
Methods
Oligonucleotides. Primers and blockers (Sequences are provided in Additional file 3) were purchased from Integrated DNA Technologies, as well as gBlocks that serve as positive controls. All DNA oligos were purchased standard desalted and dissolved in 1x TE buffer (10 mM Tris-HCl, 1 mM EDTA, pH 8). Dilutions of gBlock DNA were done in 1x TE buffer with 0.2% Tween 20 (Sigma) and 100 ng/ μl Carrier RNA(poly A)(Qiagen, Catalog No.1017647).
Repository DNA samples. Human gDNA samples (NA18537 and NA18562) were purchased from Coriell Biorepository. BRAF V600E reference standard (HD238) and Myeloid DNA reference standard (HD829) were purchased from Horizon Discovery. Seven melanoma fresh frozen patient tissue samples were purchased from OriGene. Twenty-one melanoma FFPE patient tissue samples were obtained deidentified from MD Anderson Cancer Center. Five HCC and five NSCLC fresh frozen patient tissue samples were purchased from OriGene. Fourteen HCC FFPE patient tissue samples were purchased from US Biolab and two HCC FFPE patient tissue samples were purchased from OriGene. Eighteen NSCLC FFPE patient tissue samples were purchased from OriGene. DNA from fresh frozen samples were extracted using QIAamp DNA Mini Kit (Catalog No. 51304). DNA from FFPE samples were extracted using Qiagen GeneRead DNA FFPE kit (Catalog No.180134) and repaired using NEBNext FFPE DNA repair mix (NEB #M6630L).
Stochastic Amplicon Ligation (SAL). BsaI restriction sites and complementary overhangs for assembly were appended to DNA by PCR using SAL adapter primers and Phusion Hot Start Flex DNA polymerase (NEB M0535S). For single-plex reactions, the following thermocycling protocol was used: 98 ∘C, 30s (98 ∘C, 20s; 63 ∘C, 30s; 72 ∘C, 30s) × 3 (98 ∘C, 20s; 72 ∘C, 1 min) ×n; 72 ∘C, 5 min. For multiplex panels, the following thermocycling protocol was used: 98 ∘C, 30s (98 ∘C, 20s; 63 ∘C, 2 min; 72 ∘C, 2 min) × 3 (98 ∘C, 20s; 72 ∘C, 3 min) ×n; 72 ∘C, 5 min. The number of PCR cycles (n) was empirically determined based on the amount and nature of input DNA. Amplicons from PCR were purified by column purification with Monarch PCR & DNA Cleanup Kit (NEB #T1030). DNA was quantitated by Invitrogen Qubit dsDNA HS (High Sensitivity) Assay Kit (Catalog number: Q32851). NEB BsaI-HFv2 kit (E1601S) was used for amplicon assembly. A 20 μl reaction was set up containing 400–500 ng of DNA and 2 μl golden gate enzyme mix in 1x T4 DNA ligase reaction buffer. Reaction was cycled for 30 times, with each cycle containing an incubation at 37 ∘C for 2 min and at 16 ∘C for 2 min. The assembled DNA was size selected twice using 0.4x Agencourt AMPure XP beads (Beckman Coulter A63881).
Oncogene Concatenated Enriched Amplicon Nanopore Sequencing (OCEANS). Ten to 50 ng of human genomic DNA sample was mixed with primers and blockers and subject to PCR using Phusion Hot Start Flex DNA polymerase. For single-plex tests, 400 nM of the forward and reverse primer and 4 μM blocker were used. For the 7-plex AML panel, 75 nM of each forward and reverse primer and 750 nM of each blocker were used. For the 15-plex melanoma panel and 11-plex HCC panel, 50 nM of each forward and reverse primer and 500 nM of each blocker were used. For the 28-plex NSCLC panel, 15 nM of each forward and reverse primer and 150 nM of each blocker were used. For single-plex, the following thermocycling protocol was used: 98 ∘C, 30s (98 ∘C, 10s; 63 ∘C, 30s; 72 ∘C, 30s) × 23; 72 ∘C, 5 min. For AML and Melanoma multi-plex panels, the following thermocycling protocol was used: 98 ∘C, 30s (98 ∘C, 20s; 63 ∘C, 2 min; 72 ∘C, 2 min) × 23; 72 ∘C, 5 min. For HCC and NSCLC multi-plex panels, the following thermocycling protocol was used: 98 ∘C, 30s (98 ∘C, 10s; 63 ∘C, 5 min; 72 ∘C, 30s) × 23; 72 ∘C, 5 min.
The amplicons were then purified by column purification and used as input for PCR using SAL adapter primers and subject to SAL as described above. Assembled DNA was used for library preparation using Ligation sequencing kit (SQK-LSK109) following the protocol provided by ONT. Briefly, 50–200 ng of assembled DNA was end-repaired and dA tailed using NEBNext Ultra TM II End Repair/dA-Tailing Module (E7546) and purified using 1x AMPure XP beads. Native barcoding kit (EXP-NBD104 or EXP-NBD114) was used for barcoding. Forty to 100 fmol pooled barcoded DNA was used for adapter ligation. The ligation reaction was purified using 0.5x AMPure XP beads. The beads were washed using S Fragment Buffer (SFB) and eluted in Elution Buffer (EB). The eluate was quantified by Qubit and then loaded on to Minion R9.4.1 flow cells and sequenced for 40 min. For Flongle flow cells, 20–50 fmol pooled barcoded DNA was used for adapter ligation, and sequencing was run for 8–10 h.
Next-generation sequencing (NGS) verification of clinical samples on Illumina. Ten to 50 ng of genomic DNA extracted from FFPE samples was used as input, for multiplex PCR amplicon sequencing. The following thermocycling protocol was used: 98 ∘C, 30s (98 ∘C, 20s; 63 ∘C, 2 min; 72 ∘C, 2 min) × 15; 72 ∘C, 5 min. The amplicons were purified by column purification. NEBNext ultra II DNA library prep kit for Illumina (NEB # E7645S) was used for library preparation following the kit protocol. NEBNext Multiplex Oligos for Illumina (Dual Index Primers Set 1) (NEB #E7600S) were used for index PCR. Sequencing was done on Illumina MiSeq using the V2 kit. Each sample was sequenced to at least 10,000x coverage. Libraries were checked by running on Bioanalyzer before sequencing.
ddPCR quantitation protocol. ddPCR assays were performed on QX200 Droplet Digital PCR system (Bio-Rad) using BRAF V600E (Bio-Rad dHsaMDV2010027), BRAF V600K (Bio-Rad dHsaMDV2010035), KRAS p.E62K (Bio-Rad dHsaMDS453364969), KRAS p.G13D c.38G >A (Bio-Rad dHsaMDV2510598), and MAP2K1 p.P124L (Bio-Rad dHsaMDS897199134) mutation detection kits according to kit protocols.
Nanopore Sequencing Bioinformatic analysis. Nanopore Sequencing reads were basecalled using MinKNOW 19.12.5 (fast basecalling model). Barcoded Nanopore Sequencing fastq reads were demultiplexed using EPI2ME software (Oxford Nanopore Technologies). Typically, 10,000 Nanopore Sequencing reads were used per sample for analysis. The SAL reads were deconcatenated using a custom python script [33] and the deconcatenated reads were aligned to human reference genome (GRCh38) using minimap2 aligner (https://doi.org/https://github.com/lh3/minimap2) to generate a bam file. IGVtools (https://doi.org/http://software.broadinstitute.org/software/igv/download) was used to extract the number of A, C, G, & T basecalled nucleotides, insertions, and deletions at each position of the amplicon using the basecount command in IGV-command-line tools. The VRF at each nucleotide position was calculated as the highest frequency single-base change at that position.
Bioinformatic workflow for variant calling using Clair variant caller is summarized in Figure S9. Clair v1 (https://doi.org/https://github.com/aquaskyline/clair) was used for variant calling. First, the bam file generated using minimap2 was down-sampled to <150x coverage for each amplicon using samtools view command (Subsampled read names are provided in Additional file 5). The subsampling was necessary because Clair was trained on whole genome sequencing data with low coverage. Therefore, using high coverage results in a higher false positivity rate [21]. The subsampled bam file was then used to call variants using callVarBam submodule and ont r94-flipflop model. Allele frequency threshold of 0.2 and minimum coverage of 50x were used. Variant calls were filtered for score >180.
NGS Bioinformatic analysis. NGS reads were trimmed of adapter sequences and aligned to reference amplicon sequence using short read (sr) function of minimap2 aligner. IGVtools (http://software.broadinstitute.org/software/igv/download) was used to extract the number of A, C, G, & T basecalled nucleotides, insertions, and deletions at each position of the amplicon using the basecount command in IGV-command-line tools [31]. The VRF at each nucleotide position in the enrichment region (nucleotide positions covered by the blocker in BDA) was calculated as the highest frequency single-base change at that position. Mutations in the enrichment region with VRF >5% and minimum coverage of 80x were considered true positives.
Availability of data and materials
Custom bioinformatics software was written in python for Nanopore Sequencing reads processing and is available here: https://doi.org/https://github.com/aqueous87/snippyunder opensource license Mozilla Public License 2.0 and from Zenodo [33]. Nanopore and NGS sequencing data was deposited in NCBI Sequence Read Archive (SRA) under Bioproject accession PRJNA727742 [34]. Detailed list of SRA accessions is available in Additional file 4.
References
Dowsett M, Sestak I, Lopez-Knowles E, Sidhu K, Dunbier A, Cowens J, Cuzick J. Comparison of PAM50 risk of recurrence score with onco type DX and IHC4 for predicting risk of distant recurrence after endocrine therapy. J Clin Pathol. 2013; 31(22):2783–90.
Simon R, Roychowdhury S. Implementing personalized cancer genomics in clinical trials. Nat Rev Drug Discov. 2013; 12(5):358–69.
Steuer C, Ramalingam S. Tumor mutation burden: leading immunotherapy to the era of precision medicine. J Clin Oncol. 2018; 36(7):631–2.
Tie J, et al. Circulating tumor DNA analysis detects minimal residual disease and predicts recurrence in patients with stage II colon cancer. Sci Transl Med. 2016; 8(346):346ra92.
Chabon J, et al. Integrating genomic features for non-invasive early lung cancer detection. Nature. 2020; 580(7802):245–51.
Cohen J, et al. Detection and localization of surgically resectable cancers with a multi-analyte blood test. Science. 2018; 359(6378):926–30.
Laura C, Héctor R, Carlos F. Nanopore sequencing and its application to the study of microbial communities. Comput Struct Biotechnol J. 2021; 19:1497–511.
de Lannoy C, de Ridder D, Risse J. The long reads ahead: de novo genome assembly using the MinION. F1000Research. 2017; 6:1083.
Shafin K, Pesout T, Lorig-Roach R, Haukness M, Olsen H, Bosworth C, Paten B. Nanopore sequencing and the Shasta toolkit enable efficient de novo assembly of eleven human genomes. Nat Biotechnol. 2020; 38:1044–53.
Magi A, Semeraro R, Mingrino A, Giusti B, D’Aurizio R. Nanopore sequencing data analysis: state of the art, applications and challenges. Brief Bioinform. 2018; 19:1256–72.
Orsini P, et al. Design and MinION testing of a nanopore targeted gene sequencing panel for chronic lymphocytic leukemia. Sci Rep. 2018; 8:11798.
Suzuki A, et al. Sequencing and phasing cancer mutations in lung cancers using a long-read portable sequencer. DNA Res. 2017; 24:585–96.
Alexandrov L, Kim J, Haradhvala N, Huang M, Ng A, Wu Y, Islam S. The repertoire of mutational signatures in human cancer. Nature. 2020; 578(7793):94–101.
Blomquist TM, Crawford E, Willey J, Xu J. Impact of formalin time on targeted NGS performance in FFPE tissue. Cancer Res. 2019; 79(13 Supplement):3539.
Bettegowda C, Sausen M, Leary R, Kinde I, Wang Y, Agrawal N, Antonarakis E. Detection of circulating tumor DNA in early-and late-stage human malignancies. Sci Transl Med. 2014; 6(224):224ra24.
Wu L, Chen S, Wu Y, Patel A, Zhang D. Multiplexed enrichment of rare DNA variants via sequence-selective and temperature-robust amplification. Nat Biomed Eng. 2017; 1:714–23.
Song P, Chen SX, Yan YH, et al.Selective multiplexed enrichment for the detection and quantitation of low-fraction DNA variants via low-depth sequencing. Nat Biomed Eng. 2021; 5:690–701.
Stiller M, Sucker A, Griewank K, Aust D, Baretton G, Schadendorf D, Horn S. Single-strand DNA library preparation improves sequencing of formalin-fixed and paraffin-embedded (FFPE) cancer DNA. Oncotarget. 2016; 7(37):59115.
Engler C, Kandzia R, Marillonnet S. A one pot, one step, precision cloning method with high throughput capability. PloS one. 2008; 3(11):e3647.
Volden R, Palmer T, Byrne A, Cole C, Schmitz R, Green R, Vollmers C. Improving nanopore read accuracy with the R2C2 method enables the sequencing of highly multiplexed full-length single-cell cDNA. Proc Natl Acad Sci. 2018; 115(39):9726–31.
Luo R, Wong CL, Wong YS, Tang CI, Liu CM, Leung H, Lam T. Exploring the limit of using a deep neural network on pileup data for germline variant calling. Nat Mach Intell. 2020; 2:220–27.
Miyake T, Yoshino K, Enomoto T, Takata T, Ugaki H, Kim A, Kimura T. PIK3CA gene mutations and amplifications in uterine cancers, identified by methods that avoid confounding by PIK3CA pseudogene sequences. Cancer Lett. 2008; 261(1):120–6.
Kinde I, Wu J, Papadopoulos N, Kinzler K, Vogelstein B. Detection and quantification of rare mutations with massively parallel sequencing. Proc Natl Acad Sci. 2011; 108(23):9530–5.
Karst SM, Ziels RM, Kirkegaard RH, Sorensen EA, McDonald D, Zhu Q, Knight R, Albertsen M. High-accuracy long-read amplicon sequences using unique molecular identifiers with Nanopore or PacBio sequencing. Nat Methods. 2021; 18(2):165–69.
Jain M, et al. Improved data analysis for the MinION nanopore sequencer. Nat Methods. 2015; 12:351–6.
Wilson B, Eisenstein M, Soh H. High-Fidelity Nanopore Sequencing of Ultra-Short DNA Targets. Anal. Chem. 2019; 91:6783–9.
Luo R, Sedlazeck F, Lam T-W, Schatz M. A multi-task convolutional deep neural network for variant calling in single molecule sequencing. Nat Commun. 2019; 10:1–11.
Lih C, Harrington R, Sims D, et al. Analytical validation of the next-generation sequencing assay for a nationwide signal-finding clinical trial: molecular analysis for therapy choice clinical trial. J Mol Diagn. 2017; 19(2):313–27.
Johansson G, Andersson D, Filges S, Li J, Muth A, Godfrey T, Ståhlberg A. Considerations and quality controls when analyzing cell-free tumor DNA. Biomol Detect Quantif. 2019; 17:100078.
Vendrell J, Mau-Them F, Béganton B, Godreuil S, Coopman P, Solassol J. Circulating cell free tumor DNA detection as a routine tool forLung cancer patient management. International journal of molecular sciences. 2017; 18(2):264.
Robinson J, Thorvaldsdóttir H, Winckler W, Guttman M, Lander E, Getz G, Mesirov J. Integrative genomics viewer. Nat Biotechnol. 2011; 29(1):24–6.
Poplin R, Chang P, Alexander D, et al. A universal SNP and small-indel variant caller using deep neural networks. Nat Biotechnol. 2018; 36:983–7.
Thirunavukarasu D, Y Cheng L, Song P, X Chen S, J Borad M, Kwong L, James P, J Turner D, Yu Zhang D. SNIPPY to split concatenated Nanopore reads. 2021. https://github.com/aqueous87/snippy.DOI10.5281/zenodo.5034686.
https://www.ncbi.nlm.nih.gov/sra/PRJNA727742.
Alessio M, Myrthe J, Martin E, et.al. Accurate detection of circulating tumor DNA using nanopore consensus sequencing. bioRxiv (2020.07.14.202010). 2020. https://doi.org//10.1101/2020.07.14.202010. Accessed 7 July 2021.
Oxford Nanopore Technologies. Nanopore sequencing accuracy. 2021. https://doi.org/https://nanoporetech.com/accuracy. Accessed 8 July 2021.
Acknowledgements
The authors thank Jianyi Nie for editorial assistance. The authors thank Gang Bao for use of his group’s Bio-Rad ddPCR instrument. The authors thank Oxford Nanopore Technologies for providing MinION and Flongle flow cells. The authors thank Christina Hao for bioinformatics assistance.
Peer review information
Anahita Bishop was the primary editor of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.
Review history
The review history is available as Additional file 6.
Funding
This work was supported by NIH grant R01CA203964 to DYZ.
Author information
Authors and Affiliations
Contributions
DT and DYZ conceived the project. DT designed and performed SAL feasibility and characterization experiments. DT and LYC designed the OCEANS panels, with assistance from PS for BDA primer/blocker design. DT and SXC wrote bioinformatics software for inferring VRF from OCEANS and NGS FASTQ data. LK and MJB provided clinical samples. PJ and DJT performed independent validation of OCEANS panels at Oxford Nanopore Technologies and wrote bioinformatic deconcatenation software. DT and DYZ wrote the manuscript with input from all authors. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
IRB and/or ethics committee approvals have been obtained from Rice University- IRB-FY2017-248.
Consent for publication
The authors declare that they have no competing interests.
Competing interests
Correspondence may be addressed to DYZ (dyz1@rice.edu). DT, LYC, PS, SXC, and DYZ declare competing interests in the form of consulting for NuProbe. DYZ declares a competing interest in the form of consulting for Torus Biosystems and Pana Bio. PJ and DJT are employees of Oxford Nanopore Technologies.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
Supplementary information
Additional file 2
Mutation summary - Summary of Nanopore and NGS mutation calls in clinical samples
Additional file 3
OCEANS oligo list - Oligo sequences used
Additional file 4
SRA accessions - Detailed list of SRA accessions
Additional file 5
Subsamples read names - Subsampled read names used for variant calling with Clair
Additional file 6
Review history
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Thirunavukarasu, D., Cheng, L.Y., Song, P. et al. Oncogene Concatenated Enriched Amplicon Nanopore Sequencing for rapid, accurate, and affordable somatic mutation detection. Genome Biol 22, 227 (2021). https://doi.org/10.1186/s13059-021-02449-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13059-021-02449-1