MEDUSA: Maintaining Entire DNA Duplexes for Utmost Sequencing Accuracy

Highly accurate DNA sequencing is essential for many applications. Duplex sequencing-based methods achieve unparalleled accuracy by requiring reads from both strands of the original DNA duplex to match. Yet, methods to prepare dsDNA for sequencing may resynthesize portions of each DNA duplex and cause base damage errors on one strand to become indistinguishable from true mutations on both strands. Here, we report MEDUSA (Maintaining Entire DNA Duplexes for Utmost Sequencing Accuracy) which minimizes dsDNA resynthesis to maximize duplex sequencing accuracy and yield. MEDUSA carefully repairs and blunts fragmented dsDNA, then employs apyrase to digest residual dNTPs, followed by restricted dA-tailing to prevent resynthesis. MEDUSA affords full genome coverage in a simplified protocol that is broadly compatible with dsDNA fragmentation and library preparation kits. We benchmarked MEDUSA on sheared genomic DNA and cell-free DNA and found a residual SNV frequency within a median 1.23-fold (range 0.92-1.88; p < 0.001) of what was expected if resynthesis was mostly blocked, but with full genome coverage and duplex yields within a median 1.02-fold (range 0.33 - 1.46; p = 0.258) of traditional methods that do not limit resynthesis. In all, MEDUSA could enable high breadth or depth of duplex sequencing while limiting false discovery.