Evaluation of five commercial DNA extraction kits using Salmonella as a model for implementation of rapid Nanopore sequencing in routine diagnostic laboratories.

Oxford Nanopore long-read sequencing offers advantages over Illumina short reads for the identification and characterization of bacterial pathogens for outbreak detection and surveillance activities within a diagnostic public health laboratory context. Compared to Illumina, Nanopore is more cost-effective for small batches, has a lower capital cost and has a faster turnaround time, in addition to the ability to assemble complete bacterial genomes. The quantity and quality of DNA required for Nanopore sequencing are greater than for Illumina, and the DNA extraction methods recommended for obtaining high-molecular-weight DNA are different from those typically used in diagnostic laboratories. Using a Salmonella isolate with a previously closed PacBio genome as a model Enterobacteriaceae organism, we evaluated the quantity, quality and fragmentation of five commercial DNA extraction kits. Nanopore sequencing performance was evaluated for the top three methods: Qiagen EZ1 DNA Tissue, Qiagen DNeasy Blood and Tissue, and a modified, in-house version of the MasterPure Complete DNA and RNA purification. To evaluate the effect of post-extraction DNA purification methods, we subjected extracted DNA from the three selected extraction methods to purification by AMPure beads or ethanol precipitation and compared these outputs with untreated DNA as a control. All methods are suitable for routine whole-genome sequencing (WGS), since all 60 replicates had very high genome recovery rates, with ≥98 % of the reference genome covered by mapped Nanopore reads. For 85 % of the replicates, assembly was able to produce a complete, circular chromosome using either Flye or Canu. In most cases, it is recommended to move directly from extraction to sequencing, as untreated DNA had the highest rates of genome closure regardless of extraction method. Using our evaluation criteria, the Qiagen DNeasy Blood and Tissue kit was found to be the best overall method due to its low cost, ability to scale from single tubes to 96-well plates, and high consistency in yield and sequencing performance.

Exploring the mobilome and resistome of Enterococcus faecium in a One Health context across two continents.

Enterococcus faecium is a ubiquitous opportunistic pathogen that is exhibiting increasing levels of antimicrobial resistance (AMR). Many of the genes that confer resistance and pathogenic functions are localized on mobile genetic elements (MGEs), which facilitate their transfer between lineages. Here, features including resistance determinants, virulence factors and MGEs were profiled in a set of 1273 E. faecium genomes from two disparate geographic locations (in the UK and Canada) from a range of agricultural, clinical and associated habitats. Neither lineages of E. faecium, type A and B, nor MGEs are constrained by geographic proximity, but our results show evidence of a strong association of many profiled genes and MGEs with habitat. Many features were associated with a group of clinical and municipal wastewater genomes that are likely forming a new human-associated ecotype within type A. The evolutionary dynamics of E. faecium make it a highly versatile emerging pathogen, and its ability to acquire, transmit and lose features presents a high risk for the emergence of new pathogenic variants and novel resistance combinations. This study provides a workflow for MGE-centric surveillance of AMR in Enterococcus that can be adapted to other pathogens.

Combining analytical epidemiology and genomic surveillance to identify risk factors associated with the spread of antimicrobial resistance in Salmonella enterica subsp. enterica serovar Heidelberg.

Antimicrobial resistance (AMR) has become a critical threat to public health worldwide. The use of antimicrobials in food and livestock agriculture, including the production of poultry, is thought to contribute to the dissemination of antibiotic resistant bacteria (ARB) and the genes and plasmids that confer the resistant phenotype (ARG). However, the relative contribution of each of these processes to the emergence of resistant pathogens in poultry production and their potential role in the transmission of resistant pathogens in human infections, requires a deeper understanding of the dynamics of ARB and ARG in food production and the factors involved in the increased risk of transmission.

ECTyper: in silico Escherichia coli serotype and species prediction from raw and assembled whole-genome sequence data.

Escherichia coli is a priority foodborne pathogen of public health concern and phenotypic serotyping provides critical information for surveillance and outbreak detection activities. Public health and food safety laboratories are increasingly adopting whole-genome sequencing (WGS) for characterizing pathogens, but it is imperative to maintain serotype designations in order to minimize disruptions to existing public health workflows. Multiple in silico tools have been developed for predicting serotypes from WGS data, including SRST2, SerotypeFinder and EToKi EBEis, but these tools were not designed with the specific requirements of diagnostic laboratories, which include: speciation, input data flexibility (fasta/fastq), quality control information and easily interpretable results. To address these specific requirements, we developed ECTyper (https://github.com/phac-nml/ecoli_serotyping) for performing both speciation within Escherichia and Shigella, and in silico serotype prediction. We compared the serotype prediction performance of each tool on a newly sequenced panel of 185 isolates with confirmed phenotypic serotype information. We found that all tools were highly concordant, with 92-97 % for O-antigens and 98-100 % for H-antigens, and ECTyper having the highest rate of concordance. We extended the benchmarking to a large panel of 6954 publicly available E. coli genomes to assess the performance of the tools on a more diverse dataset. On the public data, there was a considerable drop in concordance, with 75-91 % for O-antigens and 62-90 % for H-antigens, and ECTyper and SerotypeFinder being the most concordant. This study highlights that in silico predictions show high concordance with phenotypic serotyping results, but there are notable differences in tool performance. ECTyper provides highly accurate and sensitive in silico serotype predictions, in addition to speciation, and is designed to be easily incorporated into bioinformatic workflows.

Mobility of ß-lactam resistance under ampicillin treatment in gut microbiota suffering from pre-disturbance.

Ingestion of food- or waterborne antibiotic-resistant bacteria may lead to dissemination of antibiotic resistance genes (ARGs) in the gut microbiota. The gut microbiota often suffers from various disturbances. It is not clear whether and how disturbed microbiota may affect ARG mobility under antibiotic treatments. For proof of concept, in the presence or absence of streptomycin pre-treatment, mice were inoculated orally with a ß-lactam-susceptible Salmonella enterica serovar Heidelberg clinical isolate (recipient) and a ß-lactam resistant Escherichia coli O80:H26 isolate (donor) carrying a blaCMY-2 gene on an IncI2 plasmid. Immediately following inoculation, mice were treated with or without ampicillin in drinking water for 7 days. Faeces were sampled, donor, recipient and transconjugant were enumerated, blaCMY-2 abundance was determined by quantitative PCR, faecal microbial community composition was determined by 16S rRNA amplicon sequencing and cecal samples were observed histologically for evidence of inflammation. In faeces of mice that received streptomycin pre-treatment, the donor abundance remained high, and the abundance of S. Heidelberg transconjugant and the relative abundance of Enterobacteriaceae increased significantly during the ampicillin treatment. Co-blooming of the donor, transconjugant and commensal Enterobacteriaceae in the inflamed intestine promoted significantly (P<0.05) higher and possibly wider dissemination of the blaCMY-2 gene in the gut microbiota of mice that received the combination of streptomycin pre-treatment and ampicillin treatment (Str-Amp) compared to the other mice. Following cessation of the ampicillin treatment, faecal shedding of S. Heidelberg transconjugant persisted much longer from mice in the Str-Amp group compared to the other mice. In addition, only mice in the Str-Amp group shed a commensal E. coli O2:H6 transconjugant, which carries three copies of the blaCMY-2 gene, one on the IncI2 plasmid and two on the chromosome. The findings highlight the significance of pre-existing gut microbiota for ARG dissemination and persistence during and following antibiotic treatments of infectious diseases.

Rapid and accurate SNP genotyping of clonal bacterial pathogens with BioHansel.

Hierarchical genotyping approaches can provide insights into the source, geography and temporal distribution of bacterial pathogens. Multiple hierarchical SNP genotyping schemes have previously been developed so that new isolates can rapidly be placed within pre-computed population structures, without the need to rebuild phylogenetic trees for the entire dataset. This classification approach has, however, seen limited uptake in routine public health settings due to analytical complexity and the lack of standardized tools that provide clear and easy ways to interpret results. The BioHansel tool was developed to provide an organism-agnostic tool for hierarchical SNP-based genotyping. The tool identifies split k-mers that distinguish predefined lineages in whole genome sequencing (WGS) data using SNP-based genotyping schemes. BioHansel uses the Aho-Corasick algorithm to type isolates from assembled genomes or raw read sequence data in a matter of seconds, with limited computational resources. This makes BioHansel ideal for use by public health agencies that rely on WGS methods for surveillance of bacterial pathogens. Genotyping results are evaluated using a quality assurance module which identifies problematic samples, such as low-quality or contaminated datasets. Using existing hierarchical SNP schemes for Mycobacterium tuberculosis and Salmonella Typhi, we compare the genotyping results obtained with the k-mer-based tools BioHansel and SKA, with those of the organism-specific tools TBProfiler and genotyphi, which use gold-standard reference-mapping approaches. We show that the genotyping results are fully concordant across these different methods, and that the k-mer-based tools are significantly faster. We also test the ability of the BioHansel quality assurance module to detect intra-lineage contamination and demonstrate that it is effective, even in populations with low genetic diversity. We demonstrate the scalability of the tool using a dataset of ~8100 S. Typhi public genomes and provide the aggregated results of geographical distributions as part of the tool's output. BioHansel is an open source Python 3 application available on PyPI and Conda repositories and as a Galaxy tool from the public Galaxy Toolshed. In a public health context, BioHansel enables rapid and high-resolution classification of bacterial pathogens with low genetic diversity.

On-Farm Anaerobic Digestion of Dairy Manure Reduces the Abundance of Antibiotic Resistance-Associated Gene Targets and the Potential for Plasmid Transfer.

The present study investigated the impact of on-farm anaerobic digestion on the abundance of enteric bacteria, antibiotic resistance-associated gene targets, and the horizontal transfer potential of extended-spectrum ß-lactamase (ESBL) genes. Samples of raw and digested manure were obtained from six commercial dairy farms in Ontario, Canada. Digestion significantly abated populations of viable coliforms in all six farms. Conjugative transfer of plasmids carrying ß-lactamase genes from manure bacteria enriched overnight with buffered peptone containing 4 mg/liter cefotaxime into a ß-lactam-sensitive green fluorescent protein (GFP)-labeled Escherichia coli recipient strain was evaluated in patch matings. Digestion significantly decreased the frequency of the horizontal transfer of ESBL genes. Twenty-five transconjugants were sequenced, revealing six distinct plasmids, ranging in size from 40 to 180 kb. A variety of ESBL genes were identified: blaCTX-M-1, blaCTX-M-14, blaCTX-M-15, blaCTX-M-27, blaCTX-M-55, and blaPER-1. blaCTX-M-15 was the most prevalent ESBL gene detected on plasmids harbored by transconjugants. Various mobile genetic elements were found located proximal to resistance genes. Ten gene targets, including sul1, str(A), str(B), erm(B), erm(F), intI1, aadA, incW, blaPSE, and blaOXA-20, were quantified by quantitative PCR on a subset of 18 raw and 18 digested samples. Most targets were significantly more abundant in raw manure; however, erm(B) and erm(F) targets were more abundant in digested samples. Overall, on-farm digestion of dairy manure abated coliform bacteria, a number of antibiotic resistance-associated gene targets, and the potential for in vitro conjugation of plasmids conferring resistance to extended-spectrum ß-lactams and other classes of antibiotics into E. coli CV601. IMPORTANCE Using livestock manure for fertilization can entrain antibiotic-resistant bacteria into soil. Manure on some dairy farms is anaerobically digested before being land applied. Recommending the widespread implementation of the practice should be founded on understanding the impact of this treatment on various endpoints of human health concern. Although lab-scale anaerobic treatments have shown potential for reducing the abundance of antibiotic resistance genes, there are very few data from commercial farms. Anaerobic digestion of manure on six dairy farms efficiently abated coliform bacteria, E. coli, and a majority of antibiotic resistance-associated gene targets. In addition, the conjugation potential of plasmids carrying ESBL genes into introduced E. coli strain CV601 was reduced. Overall, anaerobic digestion abated coliform bacteria, the genes that they carry, and the potential for ESBL-carrying plasmid transfer.

Complete Genome Sequences for 36 Canadian Salmonella enterica Serovar Typhimurium and I 1,4,[5],12:i:- Isolates from Clinical and Animal Sources.

Here, we report the complete genome sequences for 36 Canadian isolates of Salmonella enterica subsp. enterica serovar Typhimurium and its monophasic variant I 1,4,[5]:12:i:- from both clinical and animal sources. These genome sequences will provide useful references for understanding the genetic variation within this prominent serotype.

Universal whole-sequence-based plasmid typing and its utility to prediction of host range and epidemiological surveillance.

Bacterial plasmids play a large role in allowing bacteria to adapt to changing environments and can pose a significant risk to human health if they confer virulence and antimicrobial resistance (AMR). Plasmids differ significantly in the taxonomic breadth of host bacteria in which they can successfully replicate, this is commonly referred to as 'host range' and is usually described in qualitative terms of 'narrow' or 'broad'. Understanding the host range potential of plasmids is of great interest due to their ability to disseminate traits such as AMR through bacterial populations and into human pathogens. We developed the MOB-suite to facilitate characterization of plasmids and introduced a whole-sequence-based classification system based on clustering complete plasmid sequences using Mash distances (https://github.com/phac-nml/mob-suite). We updated the MOB-suite database from 12 091 to 23 671 complete sequences, representing 17 779 unique plasmids. With advances in new algorithms for rapidly calculating average nucleotide identity (ANI), we compared clustering characteristics using two different distance measures - Mash and ANI - and three clustering algorithms on the unique set of plasmids. The plasmid nomenclature is designed to group highly similar plasmids together that are unlikely to have multiple representatives within a single cell. Based on our results, we determined that clusters generated using Mash and complete-linkage clustering at a Mash distance of 0.06 resulted in highly homogeneous clusters while maintaining cluster size. The taxonomic distribution of plasmid biomarker sequences for replication and relaxase typing, in combination with MOB-suite whole-sequence-based clusters have been examined in detail for all high-quality publicly available plasmid sequences. We have incorporated prediction of plasmid replication host range into the MOB-suite based on observed distributions of these sequence features in combination with known plasmid hosts from the literature. Host range is reported as the highest taxonomic rank that covers all of the plasmids which share replicon or relaxase biomarkers or belong to the same MOB-suite cluster code. Reporting host range based on these criteria allows for comparisons of host range between studies and provides information for plasmid surveillance.

Mobility of ß-Lactam Resistance Under Bacterial Co-infection and Ampicillin Treatment in a Mouse Model.

Ingestion of food- or waterborne antibiotic-resistant bacteria may lead to the dissemination of antibiotic-resistance genes in the gut microbiota and the development of antibiotic-resistant bacterial infection, a significant threat to animal and public health. Food or water may be contaminated with multiple resistant bacteria, but animal models on gene transfer were mainly based on single-strain infections. In this study, we investigated the mobility of ß-lactam resistance following infection with single- versus multi-strain of resistant bacteria under ampicillin treatment. We characterized three bacterial strains isolated from food-animal production systems, Escherichia coli O80:H26 and Salmonella enterica serovars Bredeney and Heidelberg. Each strain carries at least one conjugative plasmid that encodes a ß-lactamase. We orally infected mice with each or all three bacterial strain(s) in the presence or absence of ampicillin treatment. We assessed plasmid transfer from the three donor bacteria to an introduced E. coli CV601gfp recipient in the mouse gut, and evaluated the impacts of the bacterial infection on gut microbiota and gut health. In the absence of ampicillin treatment, none of the donor or recipient bacteria established in the normal gut microbiota and plasmid transfer was not detected. In contrast, the ampicillin treatment disrupted the gut microbiota and enabled S. Bredeney and Heidelberg to colonize and transfer their plasmids to the E. coli CV601gfp recipient. E. coli O80:H26 on its own failed to colonize the mouse gut. However, during co-infection with the two Salmonella strains, E. coli O80:H26 colonized and transferred its plasmid to the E. coli CV601gfp recipient and a residential E. coli O2:H6 strain. The co-infection significantly increased plasmid transfer frequency, enhanced Proteobacteria expansion and resulted in inflammation in the mouse gut. Our findings suggest that single-strain infection models for evaluating in vivo gene transfer may underrepresent the consequences of multi-strain infections following the consumption of heavily contaminated food or water.

Development of a multi-locus typing scheme for an Enterobacteriaceae linear plasmid that mediates inter-species transfer of flagella.

Due to the public health importance of flagellar genes for typing, it is important to understand mechanisms that could alter their expression or presence. Phenotypic novelty in flagellar genes arise predominately through accumulation of mutations but horizontal transfer is known to occur. A linear plasmid termed pBSSB1 previously identified in Salmonella Typhi, was found to encode a flagellar operon that can mediate phase variation, which results in the rare z66 flagella phenotype. The identification and tracking of homologs of pBSSB1 is limited because it falls outside the normal replicon typing schemes for plasmids. Here we report the generation of nine new pBSSB1-family sequences using Illumina and Nanopore sequence data. Homologs of pBSSB1 were identified in 154 genomes representing 25 distinct serotypes from 67,758 Salmonella public genomes. Pangenome analysis of pBSSB1-family contigs was performed using roary and we identified three core genes amenable to a minimal pMLST scheme. Population structure analysis based on the newly developed pMLST scheme identified three major lineages representing 35 sequence types, and the distribution of these sequence types was found to span multiple serovars across the globe. This in silico pMLST scheme has shown utility in tracking and subtyping pBSSB1-family plasmids and it has been incorporated into the plasmid MLST database under the name "pBSSB1-family".

Targeting discriminatory SNPs in Salmonella enterica serovar Heidelberg genomes using RNase H2-dependent PCR.

We report a novel RNase H2-dependent PCR (rhPCR) genotyping assay for a small number of discriminatory single-nucleotide polymorphisms (SNPs) that identify lineages and sub-lineages of the highly clonal pathogen Salmonella Heidelberg (SH). Standard PCR primers targeting numerous SNP locations were initially designed in silico, modified to be RNase H2-compatible, and then optimized by laboratory testing. Optimization often required repeated cycling through variations in primer design, assay conditions, reagent concentrations and selection of alternative SNP targets. The final rhPCR assay uses 28 independent rhPCR reactions to target 14 DNA bases that can distinguish 15 possible lineages and sub-lineages of SH. On evaluation, the assay correctly identified the 12 lineages and sub-lineages represented in a panel of 75 diverse SH strains. Non-specific amplicons were observed in 160 (15.2%) of the 1050 reactions, but due to their low intensity did not compromise assay performance. Furthermore, in silico analysis of 500 closed genomes from 103 Salmonella serovars and laboratory rhPCR testing of five prevalent Salmonella serovars including SH indicated the assay can identify Salmonella isolates as SH, since only SH isolates generated amplicons from all 14 target SNPs. The genotyping results can be fully correlated with whole genome sequencing (WGS) data in silico. This fast and economical assay, which can identify SH isolates and classify them into related or unrelated lineages and sub-lineages, has potential applications in outbreak identification, source attribution and microbial source tracking.

An improved DNA array-based classification method for the identification of Salmonella serotypes shows high concordance between traditional and genotypic testing.

Previously we developed and tested the Salmonella GenoSerotyping Array (SGSA), which utilized oligonucleotide probes for O- and H- antigen biomarkers to perform accurate molecular serotyping of 57 Salmonella serotypes. Here we describe the development and validation of the ISO 17025 accredited second version of the SGSA (SGSA v. 2) with reliable and unambiguous molecular serotyping results for 112 serotypes of Salmonella which were verified both in silico and in vitro. Improvements included an expansion of the probe sets along with a new classifier tool for prediction of individual antigens and overall serotype from the array probe intensity results. The array classifier and probe sequences were validated in silico to high concordance using 36,153 draft genomes of diverse Salmonella serotypes assembled from public repositories. We obtained correct and unambiguous serotype assignments for 31,924 (88.30%) of the tested samples and a further 3,916 (10.83%) had fully concordant antigen predictions but could not be assigned to a single serotype. The SGSA v. 2 can directly use bacterial colonies with a limit of detection of 860 CFU/mL or purified DNA template at a concentration of 1.0 x 10-1 ng/µl. The SGSA v. 2 was also validated in the wet laboratory and certified using panel of 406 samples representing 185 different serotypes with correct antigen and serotype determinations for 60.89% of the panel and 18.31% correctly identified but an ambiguous overall serotype determination.

Complete Genome and Plasmid Sequences of 32 Salmonella enterica Strains from 30 Serovars.

We report here 32 completed closed genome sequences of strains representing 30 serotypes of Salmonella. These genome sequences will provide useful references for understanding the genetic variation within Salmonella enterica serotypes, particularly as references to aid in comparative genomics studies, as well as providing information for improving in silico serotyping accuracy.

Draft Genome Assembly of a Potentially Zoonotic Cryptosporidium parvum Isolate, UKP1.

Cryptosporidium parvum is a zoonotic protozoan parasite that causes food and waterborne gastrointestinal disease and whose major animal reservoirs are cattle and small ruminants. We report here on a draft whole-genome sequence of a zoonotic isolate of C. parvum isolated from a person with cryptosporidiosis.

MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies.

Large-scale bacterial population genetics studies are now routine due to cost-effective Illumina short-read sequencing. However, analysing plasmid content remains difficult due to incomplete assembly of plasmids. Bacterial isolates can contain any number of plasmids and assembly remains complicated due to the presence of repetitive elements. Numerous tools have been developed to analyse plasmids but the performance and functionality of the tools are variable. The MOB-suite was developed as a set of modular tools for reconstruction and typing of plasmids from draft assembly data to facilitate characterization of plasmids. Using a set of closed genomes with publicly available Illumina data, the MOB-suite identified contigs of plasmid origin with both high sensitivity and specificity (95 and 88 %, respectively). In comparison, plasmidfinder demonstrated high specificity (99 %) but limited sensitivity (50 %). Using the same dataset of 377 known plasmids, MOB-recon accurately reconstructed 207 plasmids so that they were assigned to a single grouping without other plasmid or chromosomal sequences, whereas plasmidSPAdes was only able to accurately reconstruct 102 plasmids. In general, plasmidSPAdes has a tendency to merge different plasmids together, with 208 plasmids undergoing merge events. The MOB-suite reduces the number of errors but produces more hybrid plasmids, with 84 plasmids undergoing both splits and merges. The MOB-suite also provides replicon typing similar to plasmidfinder but with the inclusion of relaxase typing and prediction of conjugation potential. The MOB-suite is written in Python 3 and is available from https://github.com/phac-nml/mob-suite.

The Sequence of Two Bacteriophages with Hypermodified Bases Reveals Novel Phage-Host Interactions.

Bacteriophages SP-15 and ΦW-14 are members of the Myoviridae infecting Bacillus subtilis and Delftia (formerly Pseudomonas) acidovorans, respectively. What links them is that in both cases, approximately 50% of the thymine residues are replaced by hypermodified bases. The consequence of this is that the physico-chemical properties of the DNA are radically altered (melting temperature (Tm), buoyant density and susceptibility to restriction endonucleases). Using 454 pyrosequencing technology, we sequenced the genomes of both viruses. Phage ΦW-14 possesses a 157-kb genome (56.3% GC) specifying 236 proteins, while SP-15 is larger at 222 kb (38.6 mol % G + C) and encodes 318 proteins. In both cases, the phages can be considered genomic singletons since they do not possess BLASTn homologs. While no obvious genes were identified as being responsible for the modified base in ΦW-14, SP-15 contains a cluster of genes obviously involved in carbohydrate metabolism.

Completed Genome Sequences of Strains from 36 Serotypes of Salmonella.

We report here the completed closed genome sequences of strains representing 36 serotypes of Salmonella These genome sequences will provide useful references for understanding the genetic variation between serotypes, particularly as references for mapping of raw reads or to create assemblies of higher quality, as well as to aid in studies of comparative genomics of Salmonella.