Public health genomics capacity assessment: readiness for large-scale pathogen genomic surveillance in Canada's public health laboratories.

BACKGROUND: Along with rapid diagnostic testing, contact tracing, and public health measures, an effective pandemic response incorporates genomics-based surveillance. Large-scale SARS-CoV-2 genome sequencing is a crucial component of the global response to COVID-19. Characterizing the state of genomics readiness among Canada's public health laboratories was necessary to inform strategic planning and deployment of capacity-building resources in the early stages of the pandemic. METHODS: We used a qualitative study design and focus group discussions, encompassing both technical and leadership perspectives, to perform an in-depth evaluation of the state of pathogen genomics readiness in Canada. RESULTS: We found substantial diversity in the state of readiness for SARS-CoV-2 genomic surveillance across Canada. Despite this variability, we identified common barriers and needs in the areas of specimen access, data flow and sharing, computing infrastructure, and access to highly qualified bioinformatics personnel. CONCLUSIONS: These findings enable the strategic prioritization and deployment of resources to increase Canada's ability to perform effective public health genomic surveillance for COVID-19 and prepare for future emerging infectious diseases. They also provide a unique qualitative research model for use in capacity building.

A primer on microbial bioinformatics for nonbioinformaticians.

BACKGROUND: Presently, the bottleneck in the deployment of high-throughput sequencing technology is the ability to analyse the increasing amount of data produced in a fit-for-purpose manner. The field of microbial bioinformatics is thriving and quickly adapting to technological changes, which creates difficulties for nonbioinformaticians in following the complexity and increasingly obscure jargon of this field. AIMS: This review is directed towards nonbioinformaticians who wish to gain understanding of the overall microbial bioinformatic processes, from raw data obtained from sequencers to final outputs. SOURCES: The software and analytical strategies reviewed are based on the personal experience of the authors. CONTENT: The bioinformatic processes of transforming raw reads to actionable information in a clinical and epidemiologic context is explained. We review the advantages and limitations of two major strategies currently applied: read mapping, which is the comparison with a predefined reference genome, and de novo assembly, which is the unguided assembly of the raw data. Finally, we discuss the main analytical methodologies and the most frequently used freely available software and its application in the context of bacterial infectious disease management. IMPLICATIONS: High-throughput sequencing technologies are overhauling outbreak investigation and epidemiologic surveillance while creating new challenges due to the amount and complexity of data generated. The continuously evolving field of microbial bioinformatics is required for stakeholders to fully harness the power of these new technologies.

Neisseria gonorrhoeae Sequence Typing for Antimicrobial Resistance, a Novel Antimicrobial Resistance Multilocus Typing Scheme for Tracking Global Dissemination of N. gonorrhoeae Strains.

A curated Web-based user-friendly sequence typing tool based on antimicrobial resistance determinants in Neisseria gonorrhoeae was developed and is publicly accessible (https://ngstar.canada.ca). The N. gonorrhoeae Sequence Typing for Antimicrobial Resistance (NG-STAR) molecular typing scheme uses the DNA sequences of 7 genes (penA, mtrR, porB, ponA, gyrA, parC, and 23S rRNA) associated with resistance to ß-lactam antimicrobials, macrolides, or fluoroquinolones. NG-STAR uses the entire penA sequence, combining the historical nomenclature for penA types I to XXXVIII with novel nucleotide sequence designations; the full mtrR sequence and a portion of its promoter region; portions of ponA, porB, gyrA, and parC; and 23S rRNA sequences. NG-STAR grouped 768 isolates into 139 sequence types (STs) (n = 660) consisting of 29 clonal complexes (CCs) having a maximum of a single-locus variation, and 76 NG-STAR STs (n = 109) were identified as unrelated singletons. NG-STAR had a high Simpson's diversity index value of 96.5% (95% confidence interval [CI] = 0.959 to 0.969). The most common STs were NG-STAR ST-90 (n = 100; 13.0%), ST-42 and ST-91 (n = 45; 5.9%), ST-64 (n = 44; 5.72%), and ST-139 (n = 42; 5.5%). Decreased susceptibility to azithromycin was associated with NG-STAR ST-58, ST-61, ST-64, ST-79, ST-91, and ST-139 (n = 156; 92.3%); decreased susceptibility to cephalosporins was associated with NG-STAR ST-90, ST-91, and ST-97 (n = 162; 94.2%); and ciprofloxacin resistance was associated with NG-STAR ST-26, ST-90, ST-91, ST-97, ST-150, and ST-158 (n = 196; 98.0%). All isolates of NG-STAR ST-42, ST-43, ST-63, ST-81, and ST-160 (n = 106) were susceptible to all four antimicrobials. The standardization of nomenclature associated with antimicrobial resistance determinants through an internationally available database will facilitate the monitoring of the global dissemination of antimicrobial-resistant N. gonorrhoeae strains.

Usefulness of High-Quality Core Genome Single-Nucleotide Variant Analysis for Subtyping the Highly Clonal and the Most Prevalent Salmonella enterica Serovar Heidelberg Clone in the Context of Outbreak Investigations.

Salmonella enterica serovar Heidelberg is the second most frequently occurring serovar in Quebec and the third-most prevalent in Canada. Given that conventional pulsed-field gel electrophoresis (PFGE) subtyping for common Salmonella serovars, such as S. Heidelberg, yields identical subtypes for the majority of isolates recovered, public health laboratories are desperate for new subtyping tools to resolve highly clonal S. Heidelberg strains involved in outbreak events. As PFGE was unable to discriminate isolates from three epidemiologically distinct outbreaks in Quebec, this study was conducted to evaluate whole-genome sequencing (WGS) and phylogenetic analysis as an alternative to conventional subtyping tools. Genomes of 46 isolates from 3 Quebec outbreaks (2012, 2013, and 2014) supported by strong epidemiological evidence were sequenced and analyzed using a high-quality core genome single-nucleotide variant (hqSNV) bioinformatics approach (SNV phylogenomics [SNVphyl] pipeline). Outbreaks were indistinguishable by conventional PFGE subtyping, exhibiting the same PFGE pattern (SHEXAI.0001/SHEBNI.0001). Phylogenetic analysis based on hqSNVs extracted from WGS separated the outbreak isolates into three distinct groups, 100% concordant with the epidemiological data. The minimum and maximum number of hqSNVs between isolates from the same outbreak was 0 and 4, respectively, while >59 hqSNVs were measured between 2 previously indistinguishable outbreaks having the same PFGE and phage type, thus corroborating their distinction as separate unrelated outbreaks. This study demonstrates that despite the previously reported high clonality of this serovar, the WGS-based hqSNV approach is a superior typing method, capable of resolving events that were previously indistinguishable using classic subtyping tools.

Clostridium botulinum Group II Isolate Phylogenomic Profiling Using Whole-Genome Sequence Data.

Clostridium botulinum group II isolates (n = 163) from different geographic regions, outbreaks, and neurotoxin types and subtypes were characterized in silico using whole-genome sequence data. Two clusters representing a variety of botulinum neurotoxin (BoNT) types and subtypes were identified by multilocus sequence typing (MLST) and core single nucleotide polymorphism (SNP) analysis. While one cluster included BoNT/B4/F6/E9 and nontoxigenic members, the other comprised a wide variety of different BoNT/E subtype isolates and a nontoxigenic strain. In silico MLST and core SNP methods were consistent in terms of clade-level isolate classification; however, core SNP analysis showed higher resolution capability. Furthermore, core SNP analysis correctly distinguished isolates by outbreak and location. This study illustrated the utility of next-generation sequence-based typing approaches for isolate characterization and source attribution and identified discrete SNP loci and MLST alleles for isolate comparison.

Two novel toxin variants revealed by whole-genome sequencing of 175 Clostridium botulinum type E strains.

We sequenced 175 Clostridium botulinum type E strains isolated from food, clinical, and environmental sources from northern Canada and analyzed their botulinum neurotoxin (bont) coding sequences (CDSs). In addition to bont/E1 and bont/E3 variant types, neurotoxin sequence analysis identified two novel BoNT type E variants termed E10 and E11. Strains producing type E10 were found along the eastern coastlines of Hudson Bay and the shores of Ungava Bay, while strains producing type E11 were only found in the Koksoak River region of Nunavik. Strains producing BoNT/E3 were widespread throughout northern Canada, with the exception of the coast of eastern Hudson Bay.

Public health genomics and the new molecular epidemiology of bacterial pathogens.

Laboratory methods that can unambiguously fingerprint pathogenic microbes are needed to investigate the transmission of human infectious diseases from diverse sources, such as from the community, from the environment, within hospitals, or from contaminated food or water sources. Public health investigations currently rely on laboratory subtyping methods that ultimately provide only a fraction of the total genetic information of a pathogen, and although there is widespread success using existing subtyping methods, they do not always provide sufficient evidence to link disease cases together into outbreaks or to link these human cases to the culprit source. Alternatively, whole-genome sequencing of bacterial pathogens provides an unabridged examination of the genetic content of individual pathogen isolates, enabling public health laboratories to benefit from comparative analyses of total genetic content. In this context, whole-genome sequencing represents the ultimate epidemiological typing method - a universally applicable, highly detailed typing platform capable of providing the entire genetic blueprint of a pathogen and distinguishing strains to the single nucleotide level. These new genomic methods, if implemented within existing public health laboratory response programs, promise to revolutionize the ability of the laboratory to provide information and evidence on the evolution, transmission and virulence for bacterial pathogens - and this revolution is launching the new field of 'genomicepidemiology'.

PepTool and GeneTool: platform-independent tools for biological sequence analysis.

Although we are unable to discuss all of the functionality available in PepTool and GeneTool, it should be evident from this brief review that both packages offer a great deal in terms of functionality and ease-of-use. Furthermore, a number of useful innovations including platform-independent GUI design, networked parallelism, direct internet connectivity, database compression, and a variety of enhanced or improved algorithms should make these two programs particularly useful in the rapidly changing world of biological sequence analysis. More complete descriptions of the programs, algorithms and operation of PepTool and GeneTool are available on the BioTools web site (www.biotools.com), in the associated program user manuals and in the on-line Help pages.