Abstract:The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

Rate variation and recurrent sequence errors in pandemic-scale phylogenetics

Variation and Evolution Analysis of SARS-CoV-2 Using Self-Game Sequence Optimization

Deciphering the Subtype Differentiation History of SARS-CoV-2 Based on a New Breadth-First Searching Optimized Alignment Method Over a Global Data Set of 24,768 Sequences

Robust expansion of phylogeny for fast-growing genome sequence data

Topological data analysis identifies emerging adaptive mutations in SARS-CoV-2

Multi-scale phylodynamic modelling of rapid punctuated pathogen evolution

Addressing pandemic-wide systematic errors in the SARS-CoV-2 phylogeny

Refining SARS-CoV-2 Intra-host Variation by Leveraging Large-scale Sequencing Data

TEMPO: A Transformer-Based Mutation Prediction Framework for SARS-CoV-2 Evolution

Rapid SARS-CoV-2 surveillance using clinical, pooled, or wastewater sequence as a sensor for population change

Evolutionary analysis and lineage designation of SARS-CoV-2 genomes

Emergence of genomic diversity and recurrent mutations in SARS-CoV-2

Information Clustering and Pathogen Evolution

An accelerating, decreasing phylogenetic trend in SARS-CoV-2 genome compositional heterogeneity during the pandemic

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance

Accounting for reporting delays in real-time phylodynamic analyses with preferential sampling

Uncovering heterogeneous inter-community disease transmission from neutral allele frequency time series

Intrahost dynamics, together with genetic and phenotypic effects predict the success of viral mutations

Phylogenetic Analysis of SARS-CoV-2 Data Is Difficult

Worldwide tracing of mutations and the evolutionary dynamics of SARS CoV 2

Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows