Abstract
Most severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) diagnostic tests have relied on RNA extraction followed by reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays. Whereas automation improved logistics and different pooling strategies increased testing capacity, highly multiplexed next-generation sequencing (NGS) diagnostics remain a largely untapped resource. NGS tests have the potential to markedly increase throughput while providing crucial SARS-CoV-2 variant information. Current NGS-based detection and genotyping assays for SARS-CoV-2 are costly, mostly due to parallel sample processing through multiple steps. Here, we have established ApharSeq, in which samples are barcoded in the lysis buffer and pooled before reverse transcription. We validated this assay by applying ApharSeq to more than 500 clinical samples from the Clinical Virology Laboratory at Hadassah hospital in a robotic workflow. The assay was linear across five orders of magnitude, and the limit of detection was Ct 33 (~1000 copies/ml, 95% sensitivity) with >99.5% specificity. ApharSeq provided targeted high-confidence genotype information due to unique molecular identifiers incorporated into this method. Because of early pooling, we were able to estimate a 10- to 100-fold reduction in labor, automated liquid handling, and reagent requirements in high-throughput settings compared to current testing methods. The protocol can be tailored to assay other host or pathogen RNA targets simultaneously. These results suggest that ApharSeq can be a promising tool for current and future mass diagnostic challenges.