Standardized droplet preamplification method for downstream circulating cell-free DNA analysis

preprint OA: closed
📄 Open PDF Full text JSON View at publisher

Abstract

Circulating cell-free DNA (ccfDNA) can be found in blood and other biofluids and is a minimally invasive biomarker for several pathological processes. As tumors become more invasive, an increasing amount of circulating tumor DNA (ctDNA) is also shed into the peripheral circulation. Combined analysis of ccfDNA and ctDNA has demonstrated prognostic and predictive value in metastatic disease. However, localized tumors shed significantly less ccfDNA/ctDNA and accurate detection remains a technical challenge. To overcome this barrier, droplet preamplification has been used to perform robust multiplexed analysis of low input samples. To reduce false positives, it is essential to use a high-fidelity polymerase with 3’-5’ exonuclease activity. However, attempts to combine high-fidelity polymerases with commercial droplet digital chemistries have had limited success. There is also no standardized method for efficient amplicon recovery from droplets. In this work, we present a method to reliably stabilize emulsions and recover preamplified template. We systematically compared our protocol with different destabilization methods and found an average 41% improvement in recovery efficiency. We anticipate this standardized method will increase the consistency and reproducibility of ccfDNA/ctDNA analyses. This technique could be readily translated to other low-input or low-biomass samples, such as urine, saliva, or archived biopsy specimens.
Full text 38,923 characters · extracted from preprint-html · click to expand
Standardized droplet preamplification method for downstream circulating cell-free DNA analysis | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Standardized droplet preamplification method for downstream circulating cell-free DNA analysis Colin Skeen , View ORCID Profile Erica D. Pratt doi: https://doi.org/10.1101/2025.02.09.637351 Colin Skeen a Department of Biomedical Engineering, Boston University , Boston, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site Erica D. Pratt a Department of Biomedical Engineering, Boston University , Boston, USA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Erica D. Pratt For correspondence: epratt{at}bu.edu Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Circulating cell-free DNA (ccfDNA) can be found in blood and other biofluids and is a minimally invasive biomarker for several pathological processes. As tumors become more invasive, an increasing amount of circulating tumor DNA (ctDNA) is also shed into the peripheral circulation. Combined analysis of ccfDNA and ctDNA has demonstrated prognostic and predictive value in metastatic disease. However, localized tumors shed significantly less ccfDNA/ctDNA and accurate detection remains a technical challenge. To overcome this barrier, droplet preamplification has been used to perform robust multiplexed analysis of low input samples. To reduce false positives, it is essential to use a high-fidelity polymerase with 3’-5’ exonuclease activity. However, attempts to combine high-fidelity polymerases with commercial droplet digital chemistries have had limited success. There is also no standardized method for efficient amplicon recovery from droplets. In this work, we present a method to reliably stabilize emulsions and recover preamplified template. We systematically compared our protocol with different destabilization methods and found an average 41% improvement in recovery efficiency. We anticipate this standardized method will increase the consistency and reproducibility of ccfDNA/ctDNA analyses. This technique could be readily translated to other low-input or low-biomass samples, such as urine, saliva, or archived biopsy specimens. Introduction Circulating cell-free DNA (ccfDNA) is heavily fragmented double-stranded DNA present in low concentrations in the blood, urine, and other body fluids. In the bloodstream, ccfDNA is primarily released by hematopoietic cells undergoing apoptosis or necrosis. Since the half-life of ccfDNA is under three hours, persistence at elevated levels can be an indicator of ongoing cell death. Total ccfDNA concentration increases in several pathological processes, such as surgical trauma[ 1 ], infection[ 2 ], myocardial infarction[ 3 ], and cancer[ 4 – 9 ]. In cancer patients, a fraction of ccfDNA is tumor-derived and generally increases with disease stage[ 10 ]. This circulating tumor DNA (ctDNA) reflects the genetic heterogeneity of the primary tumor and any metastatic lesions. A growing body of literature has demonstrated the prognostic and predictive value of ctDNA in many high-incidence and high-mortality metastatic cancers, such as breast, colorectal, gastroesophageal, and pancreatic cancer. ccfDNA and ctDNA can serve as minimally invasive ‘liquid’ biopsies that provide continuous information at every stage of patient treatment. However, amplification and detection of circulating DNA in early-stage disease remains challenging. Compared to metastatic disease, both overall ccfDNA concentration and ctDNA copy number are significantly reduced[ 11 , 12 ]. Additionally, certain high-mortality cancers (e.g. glioma, pancreatic) shed much less ctDNA into the circulation than others[ 13 ]. Multiple studies have confirmed that next-generation sequencing (NGS), digital droplet PCR (ddPCR), and other PCR-reliant platforms have reduced sensitivity and increased intra-assay variability under low input conditions[ 14 – 17 ]. Stochastic amplification at low concentrations can skew the distribution of wild-type and tumor-derived sequences[ 18 , 19 ]. Sequence-dependent polymerase bias and errors can also prevent detection of rare mutations. The combination of amplification bias and PCR-induced errors can obscure rare, but clinically meaningful, mutations. One example is the CXGG ‘poison’ motif present in KRAS (codon 12) and TP53 (codon 248 and 282) mutations. These motifs are hypothesized to reduce polymerase processivity and impede amplification of mutant sequences[ 20 , 21 ]. Improving low-input assay performance could enable earlier and more effective interventions. Several pre-PCR enrichment strategies have been developed to improve ctDNA detection from low input samples. Qualitative methods selectively amplify mutant fragments[ 22 , 23 ], block wild-type ccfDNA amplification [ 24 , 25 ], or remove wild-type ccfDNA entirely using nucleases[ 26 , 27 ], CRISPR/Cas9[ 28 , 29 ], or other site-selective cleavage proteins[ 30 , 31 ]. However, these approaches cannot capture quantitative changes in ccfDNA that have been associated with progression and response to therapy. Targeted preamplification of specific loci increases total DNA input while preserving both wild-type and tumor-derived sequences[ 32 – 34 ]. We and others have developed specialized protocols to perform targeted enrichment in bulk PCR reactions[ 33 ] and in droplets[ 35 – 39 ]. Multiple studies have reported that using a high-fidelity polymerase with 3’-5’ exonuclease activity reduces misincorporation errors and improves assay performance[ 33 , 36 ]. Okada et al. recently demonstrated that droplet preamplification significantly improved ctDNA detection and more accurately predicted the presence of KRAS mutant tumors (AUC 0.855 – 0.903) compared to standard ddPCR (AUC 0.430 – 0.528) [ 40 ]. Droplet preamplification has been used for ultrasensitive profiling of early-stage cancer and precancerous lesions using plasma, fine needle aspirates, and formalin-fixed paraffin-embedded tissues[ 35 – 38 , 41 ]. However, there is no standardized protocol for ensuring specific and robust pre-PCR enrichment, or efficient amplicon recovery from droplets. The lack of standardization makes it difficult to interpret results from different studies and laboratories. An ideal pre-PCR enrichment protocol would enable both 1) unbiased amplification of multiple sequences while minimizing PCR-induced errors and 2) efficient template recovery for downstream analysis using NGS or ddPCR. Based on these criteria, the primary objectives of this study were to identify a method to stably emulsify high-fidelity PCR master mixes using Bio-Rad™ droplet digital chemistry and to reliably extract preamplified template for downstream analysis. Four different droplet destabilization strategies were evaluated based on DNA yield, DNA integrity, and protocol repeatability. Protocol performance was evaluated using contrived samples and quantified using ddPCR. Materials & methods Reagents and equipment The dMIQE checklist and experimental details can be found in Supplementary Materials. Primers and probes were used as previously described[ 36 ]. Reference materials 125 bp double-stranded gBlock gene fragments were used to simulate fragmented ccfDNA [ 42 ]. The positive control gBlock contained a KRAS G12D mutation (GAT) while the negative control contained the wild-type sequence (GGT). All gBlocks were synthesized by IDT (IA, USA). KRAS wild-type and KRAS G12D stock solutions were mixed at a 50/50 ratio to create a duplex reference standard (KRAS WT/G12D ). Preparation details are available in Supplementary Materials. Droplet preamplification Duplex preamplification reactions were prepared using the KRAS WT/G12D . Droplets were stabilized by adding 1 – 7 µL of 2x ddPCR Buffer Control for Probes (Bio-Rad, CA, USA) to a total reaction volume of 22 µL. Droplets were generated directly into thin-walled PCR tube strips using the Bio-Rad AutoDG Droplet Digital PCR System. Dropletized reactions were thermal cycled with the following PCR conditions: 3 min at 98°C, 9 cycles at 98°C for 20s, 66.5°C for 20s and 72°C for 3 min, followed by 72°C for 3 min. All steps were performed using a Bio-Rad T100 thermal cycler. Master mix details are available in Table S1.3. Droplet destabilization Digital preamplification reactions were set up as described above, but the polymerase was replaced with an equivalent volume of water. To accurately calculate the initial number of copies in each sample, an equivalent volume of KRAS WT/G12D standard was analyzed directly using ddPCR (n = 3). Thermal cycled droplets underwent three different methods of droplet disruption (chemical, physical, and thermal). For chemical and thermal protocols, the aqueous phase containing amplified template was processed using a MinElute PCR Purification kit (Qiagen, Hilden, Germany) following a modified version of the vendor’s instructions. pH indicator was not added to Buffer PB and an additional 5-minute incubation at 35°C was added prior to elution. The final elution volume was 10 µL. The kwpower function in R package MultNonParam was used to calculate sample size for each protocol (n = 24). Based on our analysis, we had 80% power to detect differences in destabilization efficiency across groups. Chemical The chloroform extraction protocol outlined in the Bio-Rad ddPCR application guide was used as a positive control [ 43 ]. Briefly, excess droplet generation oil was removed, and samples were diluted using 20 µL of TE Buffer. Next, 70 µL chloroform was added to all samples. After vortexing at max speed for 1 minute, samples were centrifuged at 15,500 x g for 10 min. The aqueous phase containing template was collected for purification. Physical Physical droplet disruption was performed using the GeneJet Purification kit (Fisher Scientific, MA, USA) as previously described[ 39 ]. 40 µL Binding Buffer and 40 µL isopropanol were added to each sample. Samples were vortexed at top speed for 1 minute, transferred to silica columns, and centrifuged for 1 minute at 16,000 x g. 700 µL Wash Buffer was added, the sample was centrifuged, and the flow-through was discarded. Columns were centrifuged again for 1 minute. 20 µL Elution Buffer was added, and columns were centrifuged for 1 minute. The physical disruption protocol simultaneously purified the template and was analyzed directly. Thermal (Thermal a ) PCR tube strips were submerged in liquid nitrogen for 1 minute, then thawed to room temperature as previously described [ 44 ]. The aqueous phase was collected, and the total volume was recorded prior to purification. Modified Thermal (Thermal b ) Based on observed results using previously reported methods, we created a revised thermal protocol. Droplets were diluted using 40 µL of TE Buffer and vortexed at max speed for 10 min. Samples were centrifuged at 2,000 x g for 1 min, snap frozen using liquid nitrogen for 1 min, and thawed at RT. This process was repeated and then the aqueous phase containing template was collected for purification. DNA integrity following thermal disruption was analyzed using TapeStation cfDNA ScreenTape (Agilent Technologies, CA, USA). ddPCR workflow All reagents were thawed to room temperature for 30 minutes before preparing ddPCR reactions (Table S1.4). Droplets were generated in semi-skirted plates and transferred to a C1000 Touch™ thermal cycler. Reactions were thermal cycled with the following PCR conditions: 10 minutes at 95°C; 40 cycles of 30s at 94°C, 60s at 61°C; 10 min. at 98°C, and a final hold at 12°C. The ramp rate was 2°C/s. Once the thermal cycler lid volume was below 50°C, plates were removed and analyzed using the QX600 Droplet Reader. Samples with more than 10,000 intact droplets were analyzed using QX Manager 2.0 (Bio-Rad). Gating thresholds were set manually based on the droplet amplitude of non-template controls (NTCs) and applied uniformly to all experimental samples. Performance metrics DNA preamplification efficiency was defined as follows: Where No. Copiesinitial was defined as the average measured concentration of the unamplified sample using ddPCR (n = 3). The dilution factor (DF) was defined as the ratio of the volume of template analyzed and total template volume. The percentage of DNA recovered was defined as follows: Statistical analysis Unless otherwise noted, data is reported as median and standard deviation (SD) using the form median (SD). To compare statistical differences between multiple samples, the non-parametric Kruskal-Wallis test was used. Results significant by Kruskal-Wallis were followed by a Dunn’s post-test with Benjamini-Hochberg correction for multiple comparisons. Comparisons of two groups were made using a non-parametric two-tailed Wilcoxon-Mann-Whitney test. Homogeneity of variances was tested using the non-parametric Fligner-Killeen test. For all hypothesis testing, the minimum level of statistical significance (p-value) was 0.05. R version 4.1.3 (The R Foundation for Statistical Computing, Vienna, Austria) was used for data analysis and figure construction. Results & discussion Stabilization of droplets containing high-fidelity polymerase reagents Prior reports using non-proprietary PCR mixes or additives with Bio-Rad™ droplet digital chemistries have reported stability issues including polydispersity, coalescence, or failure to generate droplets entirely[ 38 , 45 , 46 ]. Adding a small volume of ddPCR Supermix was reported to stabilize LAMP reactions[ 46 ]. However, the polymerase present in ddPCR Supermix has a high false-positive rate (∼0.10%) and is unsuitable for preamplification of low input samples[ 38 ]. ddPCR Buffer Control has a viscosity similar to ddPCR Supermix but does not contain a polymerase. It is sold as a blank control and is fully compatible with the Bio-Rad Automated Droplet Generator (AutoDG). Based on these findings, we tested whether Buffer Control could be used to stabilize emulsions containing high-fidelity PCR reagents. All experiments used Q5 High-Fidelity polymerase, which has been previously shown to enable efficient and accurate preamplification of low input ccfDNA samples[ 33 , 36 , 37 ]. To evaluate the impact of Buffer Control on droplet stability, a titration series experiment was conducted (1 µL – 7 µL) and droplets were thermal cycled and evaluated by visual inspection. Without Buffer Control, droplets were successfully generated but would coalesce after thermal cycling ( Fig 1A ). A discrete droplet layer was observed starting at 2 µL and droplet homogeneity appeared to increase with Buffer Control added. To identify the minimum required volume of Buffer Control, we measured the preamplification efficiency across a narrower range (1.5 – 4.5 µL). Each reaction contained 2,850 copies of the KRAS WT/G12D reference standard. Amplification efficiency increased with Buffer Control added until 3.5 µL, after which there was no statistically significant difference ( Fig 1B ). Two-dimensional cluster plots of the duplex reaction showed clear separation of the two molecular targets ( Fig 1C ). Non-template and wild-type only controls showed no observable increase in false positives compared to conventional ddPCR (Fig S3, Fig S4). We selected 4 µL of Buffer Control as the optimized volume to ensure small variations in sample quality would not impact assay performance. The median amplification ratio for the 4 µL group was 60.15 (14.27%) which matches our prior results[ 36 ]. The preamplification efficiency of KRAS WT and G12D were similar (p = 0.297; Wilcoxon-Mann-Whitney), suggesting no significant amplification bias ( Fig 1D ). Based on these results, we used the optimized preamplification protocol for subsequent experiments. Download figure Open in new tab Figure 1. Optimization of droplet preamplification reaction. A) Representative droplets generated using Q5 high-fidelity polymerase and varying amounts of Bio-Rad Buffer Control for Probes before (top) and after (bottom) thermal cycling. B) Efficiency of PCR preamplification as a function of Buffer Control added to the reaction mix (Kruskal-Wallis). C) Representative 2D plot for pessimal (left) and optimal (right) reaction mixes. D) Comparison of total number of amplified copies for KRAS WT and KRAS G12D in a duplex reaction. ****: p<0.0001; ***: p<0.001; **: p<0.01; *: p<0.05. Comparison of droplet destabilization methods For accurate detection of low abundance ctDNA in early-stage cancer, efficient recovery of amplicons from droplets is critical. Bio-Rad recommends a chloroform destabilization protocol to extract DNA from droplets [ 43 ]. More recently, physical [ 39 ] and thermal [ 44 ] destabilization protocols have reported similar DNA recoveries with fewer processing steps. Chloroform displaces surfactants from the oil-water droplet interface, inducing coalescence. However, chemical destabilizers are partially soluble in the recovered aqueous layer and can reduce the efficiency of downstream PCR reactions [ 47 ]. Thermal destabilization methods use liquid nitrogen freeze-thaws to expand the inner aqueous layer, causing the interfacing surfactant to become sparser and decreasing droplet stability. Physical destabilization requires passing the droplets through a DNA spin column with centrifugation steps to physically disrupt the droplets. Both methods reduce opportunities for cross-contamination and sample loss. To determine which method worked best with our optimized protocol, we compared the number of copies present in a reference standard of known concentration before and after dropletization and destabilization ( Fig 2A ). Preamplification reactions were set up as previously described, but the polymerase was replaced with an equivalent volume of water. This yielded stabilized DNA-containing droplets that could be thermal cycled, but no amplification would take place. In our hands, the Bio-Rad recommended chemical protocol outperformed both physical and thermal protocols ( Table 1 ). Our results are in line with the previously reported recoveries for chemical (45%) and physical (36%) protocols [ 44 ]. However, our measured recovery for thermal destabilization (Thermal a ) was significantly lower than previously reported (27.4% versus 68%). A recent study using droplet preamplification to improve NGS-based ccfDNA analysis also reported lower recovery when using thermal destabilization[ 39 ]. Our master mix modifications to emulsify high-fidelity reagents could increase droplet stability, which might explain this difference in results. View this table: View inline View popup Download powerpoint Table 1. Efficiency and variability of four different droplet destabilization methods. 1 standard deviation 2 coefficient of variation Download figure Open in new tab Figure 2. Evaluation of droplet destabilization methods. A) Experimental design for evaluating each protocol. B) The proportion of input template recovered using chemical, physical, and thermal protocols (Kruskal Wallis; Dunn’s post-hoc). ****: p<0.0001; ***: p<0.001; **: p<0.01; *: p<0.05. We then created a hybrid protocol containing elements of the chemical and thermal protocols, but only using TE buffer as a diluent. We found that the new thermal protocol (Thermal b ) had the highest percent recovery compared to prior methods (p <0.00001; Kruskal-Wallis) ( Fig 2B ). The introduction of a dilution and vortexing step did not affect assay repeatability, as variance in DNA recovery was not significantly different across the four methods (p < 0.05; Fligner-Killeen). We also tested variables such as centrifugation speed, and altered the order of the dilution, vortexing, and freeze-thaw steps. None of these modifications improved DNA recovery (Fig S5). Evaluation of DNA integrity A concern with thermal methods is that ccfDNA can be damaged by multiple freeze-thaw cycles [ 48 , 49 ]. To evaluate whether the two liquid nitrogen steps in the Thermal b protocol impacted DNA quality, we analyzed gBlock samples before and after destabilization using capillary electrophoresis. The original 125 bp fragment and 95 bp amplicon were clearly distinguishable with no evidence of additional fragmentation or degradation ( Fig 3 , Fig S6). Download figure Open in new tab Figure 3: Representative Tapestation electropherograms of gBlocks pre- and post-droplet preamplification. D1: 125 bp KRAS gBlock before preamplification. H1 & E1: Samples taken after thermal destabilization (protocol Thermal b ). Only two peaks are visible, the original unamplified gBlock and the expected 95 bp amplicon. Lower is the lower marker used to align ladder data with samples of unknown size. Conclusion The goal of this study was to develop a unified protocol for ccfDNA preamplification using Bio-Rad droplet digital chemistry. To stabilize high-fidelity polymerase reagents, we leveraged existing Bio-Rad reagents. This approach is fully compatible with automated droplet-generator systems and minimizes any chance of contamination or line fouling. To recover preamplified template, we developed a destabilization protocol with an efficiency of 68.6%. Using gBlocks as contrived ccfDNA samples allowed for standardized evaluation of our method and previously published techniques. Our approach outperformed chemical, physical, and thermal methods by 33%, 53%, and 150%, respectively. Future Perspective Droplet preamplification can resolve subsampling effects that are a key barrier preventing accurate quantification and analysis of low input plasma samples. This allows for repeated testing of each plasma sample without exhausting the original specimen, improving assay robustness and enabling analyses using orthogonal methods. Future studies could implement synthetic spike-in controls[ 49 ] to quantify amplification efficiency and DNA recovery from droplets with even higher accuracy. While this study leverages existing commercial ddPCR reagents, this approach could be readily adapted to other droplet microfluidics platforms. This pre-PCR enrichment method maximizes the amount of information that can be obtained in a single blood draw. Acknowledgements We thank Dr. Susan Tsai, Dr. Poojitha Sitaram, and Julia Sexton for helpful discussion, suggestions, and early technical assistance. Footnotes Funding: This work was supported in part by ACS IRG 22-153-42 from the American Cancer Society. Minor in-text edits for clarity and readability. References [1]. ↵ Konishi S , Narita T , Hatakeyama S , et al. Utility of total cell-free DNA levels for surgical damage evaluation in patients with urological surgeries . Sci Rep . 2021 ; 11 ( 1 ): 22103 . OpenUrl CrossRef PubMed [2]. ↵ Hu Y , Zhao Y , Zhang Y , et al. Cell-free DNA: a promising biomarker in infectious diseases . Trends Microbiol . 2024 ;S0966-842X(24)00168-9. [3]. ↵ Chang CP-Y , Chia R-H , Wu T-L , et al. Elevated cell-free serum DNA detected in patients with myocardial infarction . Clin Chim Acta Int J Clin Chem . 2003 ; 327 ( 1–2 ): 95 – 101 . OpenUrl CrossRef [4]. ↵ Diehl F , Schmidt K , Choti MA , et al. Circulating mutant DNA to assess tumor dynamics . Nat Med . 2008 ; 14 ( 9 ): 985 – 990 . OpenUrl CrossRef PubMed Web of Science [5]. Dawson S-J , Tsui DWY , Murtaza M , et al. Analysis of Circulating Tumor DNA to Monitor Metastatic Breast Cancer . N Engl J Med . 2013 ; 368 ( 13 ): 1199 – 1209 . OpenUrl CrossRef PubMed Web of Science [6]. Thompson JC , Yee SS , Troxel AB , et al. Detection of therapeutically targetable driver and resistance mutations in lung cancer patients by next generation sequencing of cell-free circulating tumor DNA . Clin Cancer Res Off J Am Assoc Cancer Res . 2016 ; 22 ( 23 ): 5772 – 5782 . OpenUrl CrossRef [7]. Rozenblum AB , Ilouze M , Dudnik E , et al. Clinical Impact of Hybrid Capture-Based Next-Generation Sequencing on Changes in Treatment Decisions in Lung Cancer . J Thorac Oncol Off Publ Int Assoc Study Lung Cancer . 2017 ; 12 ( 2 ): 258 – 268 . OpenUrl [8]. Kim ST , Banks KC , Lee S-H , et al. Prospective Feasibility Study for Using Cell-Free Circulating Tumor DNA–Guided Therapy in Refractory Metastatic Solid Cancers: An Interim Analysis . JCO Precis Oncol . 2017 ;( 1 ): 1 – 15 . [9]. ↵ Schwaederlé MC , Patel SP , Husain H , et al. Utility of genomic assessment of blood-derived circulating tumor DNA (ctDNA) in patients with advanced lung adenocarcinoma . Clin Cancer Res Off J Am Assoc Cancer Res . 2017 ; 23 ( 17 ): 5101 – 5111 . OpenUrl CrossRef [10]. ↵ Bettegowda C , Sausen M , Leary RJ , et al. Detection of Circulating Tumor DNA in Early- and Late-Stage Human Malignancies . Sci Transl Med . 2014 ; 6 ( 224 ): 224ra24 – 224ra24 . OpenUrl Abstract / FREE Full Text [11]. ↵ Kamat AA , Baldwin M , Urbauer D , et al. Plasma cell-free DNA in ovarian cancer . Cancer . 2010 ; 116 ( 8 ): 1918 – 1925 . OpenUrl CrossRef PubMed Web of Science [12]. ↵ Bedin C , Enzo MV , Del Bianco P , et al. Diagnostic and prognostic role of cell-free DNA testing for colorectal cancer patients . Int J Cancer . 2017 ; 140 ( 8 ): 1888 – 1898 . OpenUrl CrossRef PubMed [13]. ↵ Husain H , Pavlick DC , Fendler BJ , et al. Tumor Fraction Correlates With Detection of Actionable Variants Across > 23,000 Circulating Tumor DNA Samples . JCO Precis Oncol . 2022 ;( 6 ): e2200261 . [14]. ↵ Brychta N , Krahn T , von Ahsen O. Detection of KRAS Mutations in Circulating Tumor DNA by Digital PCR in Early Stages of Pancreatic Cancer . Clin Chem . 2016 ; 62 ( 11 ): 1482 – 1491 . OpenUrl Abstract / FREE Full Text [15]. McNulty SN , Mann PR , Robinson JA , et al. Impact of Reducing DNA Input on Next-Generation Sequencing Library Complexity and Variant Detection . J Mol Diagn . 2020 ; 22 ( 5 ): 720 – 727 . OpenUrl CrossRef PubMed [16]. Yu L , Lopez G , Rassa J , et al. Direct comparison of circulating tumor DNA sequencing assays with targeted large gene panels . PLOS ONE . 2022 ; 17 ( 4 ): e0266889 . OpenUrl CrossRef PubMed [17]. ↵ Zhang Y , Xu Y , Zhong W , et al. Total DNA input is a crucial determinant of the sensitivity of plasma cell-free DNA EGFR mutation detection using droplet digital PCR . Oncotarget . 2016 ; 8 ( 4 ): 5861 – 5873 . OpenUrl [18]. ↵ Akbari M , Doré Hansen M , Halgunset J , et al. Low Copy Number DNA Template Can Render Polymerase Chain Reaction Error Prone in a Sequence-Dependent Manner . J Mol Diagn JMD . 2005 ; 7 ( 1 ): 36 – 39 . OpenUrl CrossRef PubMed [19]. ↵ Kebschull JM , Zador AM . Sources of PCR-induced distortions in high-throughput sequencing data sets . Nucleic Acids Res . 2015 ; 43 ( 21 ): e143 . OpenUrl CrossRef PubMed [20]. ↵ Barnard R , Futo V , Pecheniuk N , et al. PCR Bias Toward the Wild-Type k-ras and p53 Sequences: Implications for PCR Detection of Mutations and Cancer Diagnosis . BioTechniques . 1998 ; 25 ( 4 ): 684 – 691 . OpenUrl CrossRef PubMed Web of Science [21]. ↵ Liu Q , Thorland EC , Sommer SS . Inhibition of PCR Amplification by a Point Mutation Downstream of a Primer . BioTechniques . 1997 ; 22 ( 2 ): 292 – 300 . OpenUrl CrossRef PubMed Web of Science [22]. ↵ Li J , Wang L , Mamon H , et al. Replacing PCR with COLD-PCR enriches variant DNA sequences and redefines the sensitivity of genetic testing . Nat Med . 2008 ; 14 ( 5 ): 579 – 584 . OpenUrl CrossRef PubMed Web of Science [23]. ↵ How Kit A , Mazaleyrat N , Daunay A , et al. Sensitive detection of KRAS mutations using enhanced-ice-COLD-PCR mutation enrichment and direct sequence identification . Hum Mutat . 2013 ; 34 ( 11 ): 1568 – 1580 . OpenUrl CrossRef PubMed [24]. ↵ Dominguez PL , Kolodney MS . Wild-type blocking polymerase chain reaction for detection of single nucleotide minority mutations from clinical specimens . Oncogene . 2005 ; 24 ( 45 ): 6830 – 6834 . OpenUrl CrossRef PubMed Web of Science [25]. ↵ Oldenburg RP , Liu MS , Kolodney MS . Selective amplification of rare mutations using locked nucleic acid oligonucleotides that competitively inhibit primer binding to wild-type DNA . J Invest Dermatol . 2008 ; 128 ( 2 ): 398 – 402 . OpenUrl CrossRef PubMed Web of Science [26]. ↵ Zhulidov PA , Bogdanova EA , Shcheglov AS , et al. Simple cDNA normalization using kamchatka crab duplex-specific nuclease . Nucleic Acids Res . 2004 ; 32 ( 3 ): e37 . OpenUrl CrossRef PubMed [27]. ↵ Song C , Liu Y , Fontana R , et al. Elimination of unaltered DNA in mixed clinical samples via nuclease-assisted minor-allele enrichment . Nucleic Acids Res . 2016 ; 44 ( 19 ): e146 . OpenUrl CrossRef PubMed [28]. ↵ Gu W , Crawford ED , O’Donovan BD , et al. Depletion of Abundant Sequences by Hybridization (DASH): using Cas9 to remove unwanted high-abundance species in sequencing libraries and molecular counting applications . Genome Biol . 2016 ; 17 : 41 . OpenUrl CrossRef PubMed [29]. ↵ Phallen J , Sausen M , Adleff V , et al. Direct detection of early-stage cancers using circulating tumor DNA . Sci Transl Med . 2017 ; 9 ( 403 ): eaan2415 . OpenUrl Abstract / FREE Full Text [30]. ↵ Song J , Hegge JW , Mauk MG , et al. Highly specific enrichment of rare nucleic acid fractions using Thermus thermophilus argonaute with applications in cancer diagnostics . Nucleic Acids Res . 2020 ; 48 ( 4 ): e19 . OpenUrl CrossRef PubMed [31]. ↵ Liu Q , Guo X , Xun G , et al. Argonaute integrated single-tube PCR system enables supersensitive detection of rare mutations . Nucleic Acids Res . 2021 ; 49 ( 13 ): e75 . OpenUrl CrossRef PubMed [32]. ↵ Livak KJ , Wills QF , Tipping AJ , et al. Methods for qPCR gene expression profiling applied to 1440 lymphoblastoid single cells . Methods San Diego Calif . 2013 ; 59 ( 1 ): 71 – 79 . OpenUrl CrossRef PubMed Web of Science [33]. ↵ Jackson JB , Choi DS , Luketich JD , et al. Multiplex Preamplification of Serum DNA to Facilitate Reliable Detection of Extremely Rare Cancer Mutations in Circulating DNA by Digital PCR . J Mol Diagn . 2016 ; 18 ( 2 ): 235 – 243 . OpenUrl CrossRef PubMed [34]. ↵ Guan Y , Mayba O , Sandmann T , et al. High-Throughput and Sensitive Quantification of Circulating Tumor DNA by Microfluidic-Based Multiplex PCR and Next-Generation Sequencing . J Mol Diagn JMD . 2017 ; 19 ( 6 ): 921 – 932 . OpenUrl CrossRef PubMed [35]. ↵ Ono Y , Sugitani A , Karasaki H , et al. An improved digital polymerase chain reaction protocol to capture low-copy KRAS mutations in plasma cell-free DNA by resolving ‘subsampling’ issues . Mol Oncol . 2017 ; 11 ( 10 ): 1448 – 1458 . OpenUrl CrossRef PubMed [36]. ↵ Pratt ED , Cowan RW , Manning SL , et al. Multiplex Enrichment and Detection of Rare KRAS Mutations in Liquid Biopsy Samples using Digital Droplet Pre-Amplification . Anal Chem . 2019 ; 91 ( 12 ): 7516 – 7523 . OpenUrl CrossRef [37]. ↵ Pratt ED , Zhen DB , Cowan RW , et al. Ultra-low Input Circulating Tumor DNA Detection by MED-Amp in Early-Stage Pancreatic Cancer [Internet] . bioRxiv ; 2021 [cited 2022 Sep 26]. p. 2021.03.28.437388. Available from: https://www.biorxiv.org/content/10.1101/2021.03.28.437388v1 . [38]. ↵ Maeda C , Ono Y , Hayashi A , et al. Multiplex Digital PCR Assay to Detect Multiple KRAS and GNAS Mutations Associated with Pancreatic Carcinogenesis from Minimal Specimen Amounts . J Mol Diagn JMD . 2023 ; 25 ( 6 ): 367 – 377 . OpenUrl CrossRef PubMed [39]. ↵ Kinoo A , Caye-Eude A , Vial Y , et al. Partitioning for Easy Multiplexing . J Mol Diagn . 2023 ;S152515782300123X. [40]. ↵ Okada T , Mizukami Y , Ono Y , et al. Digital PCR-based plasma cell-free DNA mutation analysis for early-stage pancreatic tumor diagnosis and surveillance . J Gastroenterol . 2020 ; 55 ( 12 ): 1183 – 1193 . OpenUrl CrossRef PubMed [41]. ↵ Kopp JL , Dubois CL , Schaeffer DF , et al. Loss of Pten and Activation of Kras Synergistically Induce Formation of Intraductal Papillary Mucinous Neoplasia From Pancreatic Ductal Cells in Mice . Gastroenterology . 2018 ; 154 ( 5 ): 1509 - 1523.e5 . OpenUrl CrossRef PubMed [42]. ↵ Underhill HR , Kitzman JO , Hellwig S , et al. Fragment Length of Circulating Tumor DNA . PLOS Genet . 2016 ; 12 ( 7 ): e1006162 . OpenUrl CrossRef PubMed [43]. ↵ Bio-Rad. Droplet DigitalTM PCR Droplet DigitalTM PCR Applications Guide. [Internet] . 2018 . Available from: https://www.bio-rad.com/webroot/web/pdf/lsr/literature/Bulletin_6407.pdf . [44]. ↵ Dutra L , Franz O , Puupponen V-M , et al. DNA recovery from Droplet DigitalTM PCR emulsions using liquid nitrogen . BioTechniques . 2020 ; 69 ( 6 ): 450 – 454 . OpenUrl PubMed [45]. ↵ Corless BC , Chang GA , Cooper S , et al. Development of Novel Mutation-Specific Droplet Digital PCR Assays Detecting TERT Promoter Mutations in Tumor and Plasma Samples . J Mol Diagn . 2019 ; 21 ( 2 ): 274 – 285 . OpenUrl CrossRef PubMed [46]. ↵ Hussain M. Isothermal droplet digital PCR method for quantification of CHO residual DNA . J Pharm Biomed Anal . 2022 ; 211 : 114564 . OpenUrl CrossRef PubMed [47]. ↵ Karbaschi M , Shahi P , Abate AR . Rapid, chemical-free breaking of microfluidic emulsions with a hand-held antistatic gun . Biomicrofluidics . 2017 ; 11 ( 4 ): 044107 . OpenUrl CrossRef PubMed [48]. ↵ Chan KCA . Effects of Preanalytical Factors on the Molecular Size of Cell-Free DNA in Blood | Clinical Chemistry | Oxford Academic [Internet] . [cited 2023 Sep 13]. Available from: https://academic.oup.com/clinchem/article/51/4/781/5629599?login=false . [49]. ↵ El Messaoudi S , Rolet F , Mouliere F , et al. Circulating cell free DNA: Preanalytical considerations . Clin Chim Acta . 2013 ; 424 : 222 – 230 . OpenUrl CrossRef PubMed Web of Science View the discussion thread. Back to top Previous Next Posted February 11, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Standardized droplet preamplification method for downstream circulating cell-free DNA analysis Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Standardized droplet preamplification method for downstream circulating cell-free DNA analysis Colin Skeen , Erica D. Pratt bioRxiv 2025.02.09.637351; doi: https://doi.org/10.1101/2025.02.09.637351 Share This Article: Copy Citation Tools Standardized droplet preamplification method for downstream circulating cell-free DNA analysis Colin Skeen , Erica D. Pratt bioRxiv 2025.02.09.637351; doi: https://doi.org/10.1101/2025.02.09.637351 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Bioengineering Subject Areas All Articles Animal Behavior and Cognition (7637) Biochemistry (17705) Bioengineering (13899) Bioinformatics (41968) Biophysics (21460) Cancer Biology (18603) Cell Biology (25526) Clinical Trials (138) Developmental Biology (13385) Ecology (19910) Epidemiology (2067) Evolutionary Biology (24328) Genetics (15614) Genomics (22513) Immunology (17741) Microbiology (40423) Molecular Biology (17193) Neuroscience (88646) Paleontology (667) Pathology (2835) Pharmacology and Toxicology (4827) Physiology (7647) Plant Biology (15160) Scientific Communication and Education (2046) Synthetic Biology (4302) Systems Biology (9825) Zoology (2271)

Text is read by the "Ask this paper" AI Q&A widget below. Extraction quality varies by source — PMC NXML preserves structure cleanly, OA-HTML may include some navigation residue, and OA-PDF can have broken hyphenation. The publisher copy (via DOI) is the canonical version.

My notes (saved in your browser only)

Ask this paper AI returns verbatim quotes from the full text · source: preprint-html

Answers must be backed by verbatim quotes from this paper's full text. Hallucinated quotes are dropped automatically; if no verbatim passage answers the question, we say so. How this works

Citation neighborhood (no data yet)

We don't have any in-corpus citations linked to this paper yet. This is a recent paper (2025) — citers typically take a year or two to land, and the OpenAlex reference graph may still be filling in.

Source provenance

europepmc
last seen: 2026-05-20T01:45:00.602351+00:00