Correlation assessment of SARS-CoV-2 variants and their subvariants present in clinical and wastewater samples in Oregon, USA (February 7, 2021 - February 26, 2022) using the Freyja bioinformatics approach

Background Wastewater surveillance is a valuable tool for monitoring SARS-CoV-2 at the community level. As the virus diversified into many variants and subvariants that share overlapping mutations, resolving them accurately from wastewater becomes a key bioinformatic challenge. Objectives and aims This study evaluated two distinct bioinformatic approaches, multilocus sequence typing (MLST) and Freyja, for identifying SARS-CoV-2 variants and subvariants in Oregon wastewater samples collected from February 2021 to February 2022. Methods The MLST approach identified SARS-CoV-2 variants using unique mutations curated from clinical samples. In contrast, the Freyja approach resolved variant and subvariant abundances using genome wide mutation profiles weighted by sequencing depth. In this study, the variant and subvariants relative abundances produced by both approaches were compared against those observed in clinical surveillance data. Results Both approaches identified SARS-CoV-2 variants at relative abundances that agreed closely with those observed in clinical surveillance data. However, only the Freyja approach identified over 200 Delta subvariants, divided into three clades (21A, 21I and 21J) and two levels (Level 1 and 2) based on Pango subvariants. Delta subvariants showed strong agreement at Level 1 subvariants (r s = 0.892–0.944), while agreement at Level 2 subvariants was inconsistent (r s = 0.324–0.903). Conclusions The Freyja approach provided enhanced resolution of SARS-CoV-2 variants and subvariants in wastewater, at abundances that agreed with clinical surveillance. This added resolution is a critical advantage for public health surveillance as SARS-CoV-2 continues to evolve and share mutations across variants and subvariants.