Correlation between Phenotypic and In Silico Detection of Antimicrobial Resistance in Salmonella enterica in Canada Using Staramr.

Whole genome sequencing (WGS) of Salmonella supports both molecular typing and detection of antimicrobial resistance (AMR). Here, we evaluated the correlation between phenotypic antimicrobial susceptibility testing (AST) and in silico prediction of AMR from WGS in Salmonella enterica (n = 1321) isolated from human infections in Canada. Phenotypic AMR results from broth microdilution testing were used as the gold standard. To facilitate high-throughput prediction of AMR from genome assemblies, we created a tool called Staramr, which incorporates the ResFinder and PointFinder databases and a custom gene-drug key for antibiogram prediction. Overall, there was 99% concordance between phenotypic and genotypic detection of categorical resistance for 14 antimicrobials in 1321 isolates (18,305 of 18,494 results in agreement). We observed an average sensitivity of 91.2% (range 80.5-100%), a specificity of 99.7% (98.6-100%), a positive predictive value of 95.4% (68.2-100%), and a negative predictive value of 99.1% (95.6-100%). The positive predictive value of gentamicin was 68%, due to seven isolates that carried aac(3)-IVa, which conferred MICs just below the breakpoint of resistance. Genetic mechanisms of resistance in these 1321 isolates included 64 unique acquired alleles and mutations in three chromosomal genes. In general, in silico prediction of AMR in Salmonella was reliable compared to the gold standard of broth microdilution. WGS can provide higher-resolution data on the epidemiology of resistance mechanisms and the emergence of new resistance alleles.

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.

Shared genome analyses of notable listeriosis outbreaks, highlighting the critical importance of epidemiological evidence, input datasets and interpretation criteria.

The persuasiveness of genomic evidence has pressured scientific agencies to supplement or replace well-established methodologies to inform public health and food safety decision-making. This study of 52 epidemiologically defined Listeria monocytogenes isolates, collected between 1981 and 2011, including nine outbreaks, was undertaken (1) to characterize their phylogenetic relationship at finished genome-level resolution, (2) to elucidate the underlying genetic diversity within an endemic subtype, CC8, and (3) to re-evaluate the genetic relationship and epidemiology of a CC8-delimited outbreak in Canada in 2008. Genomes representing Canadian Listeria outbreaks between 1981 and 2010 were closed and manually annotated. Single nucleotide variants (SNVs) and horizontally acquired traits were used to generate phylogenomic models. Phylogenomic relationships were congruent with classical subtyping and epidemiology, except for CC8 outbreaks, wherein the distribution of SNV and prophages revealed multiple co-evolving lineages. Chronophyletic reconstruction of CC8 evolution indicates that prophage-related genetic changes among CC8 strains manifest as PFGE subtype reversions, obscuring the relationship between CC8 isolates, and complicating the public health interpretation of subtyping data, even at maximum genome resolution. The size of the shared genome interrogated did not change the genetic relationship measured between highly related isolates near the tips of the phylogenetic tree, illustrating the robustness of these approaches for routine public health applications where the focus is recent ancestry. The possibility exists for temporally and epidemiologically distinct events to appear related even at maximum genome resolution, highlighting the continued importance of epidemiological evidence.

Neptune: a bioinformatics tool for rapid discovery of genomic variation in bacterial populations.

The ready availability of vast amounts of genomic sequence data has created the need to rethink comparative genomics algorithms using 'big data' approaches. Neptune is an efficient system for rapidly locating differentially abundant genomic content in bacterial populations using an exact k-mer matching strategy, while accommodating k-mer mismatches. Neptune's loci discovery process identifies sequences that are sufficiently common to a group of target sequences and sufficiently absent from non-targets using probabilistic models. Neptune uses parallel computing to efficiently identify and extract these loci from draft genome assemblies without requiring multiple sequence alignments or other computationally expensive comparative sequence analyses. Tests on simulated and real datasets showed that Neptune rapidly identifies regions that are both sensitive and specific. We demonstrate that this system can identify trait-specific loci from different bacterial lineages. Neptune is broadly applicable for comparative bacterial analyses, yet will particularly benefit pathogenomic applications, owing to efficient and sensitive discovery of differentially abundant genomic loci. The software is available for download at: http://github.com/phac-nml/neptune.

SNVPhyl: a single nucleotide variant phylogenomics pipeline for microbial genomic epidemiology.

The recent widespread application of whole-genome sequencing (WGS) for microbial disease investigations has spurred the development of new bioinformatics tools, including a notable proliferation of phylogenomics pipelines designed for infectious disease surveillance and outbreak investigation. Transitioning the use of WGS data out of the research laboratory and into the front lines of surveillance and outbreak response requires user-friendly, reproducible and scalable pipelines that have been well validated. Single Nucleotide Variant Phylogenomics (SNVPhyl) is a bioinformatics pipeline for identifying high-quality single-nucleotide variants (SNVs) and constructing a whole-genome phylogeny from a collection of WGS reads and a reference genome. Individual pipeline components are integrated into the Galaxy bioinformatics framework, enabling data analysis in a user-friendly, reproducible and scalable environment. We show that SNVPhyl can detect SNVs with high sensitivity and specificity, and identify and remove regions of high SNV density (indicative of recombination). SNVPhyl is able to correctly distinguish outbreak from non-outbreak isolates across a range of variant-calling settings, sequencing-coverage thresholds or in the presence of contamination. SNVPhyl is available as a Galaxy workflow, Docker and virtual machine images, and a Unix-based command-line application. SNVPhyl is released under the Apache 2.0 license and available at http://snvphyl.readthedocs.io/ or at https://github.com/phac-nml/snvphyl-galaxy.

Bacillus cereus Biovar Anthracis Causing Anthrax in Sub-Saharan Africa-Chromosomal Monophyly and Broad Geographic Distribution.

Through full genome analyses of four atypical Bacillus cereus isolates, designated B. cereus biovar anthracis, we describe a distinct clade within the B. cereus group that presents with anthrax-like disease, carrying virulence plasmids similar to those of classic Bacillus anthracis. We have isolated members of this clade from different mammals (wild chimpanzees, gorillas, an elephant and goats) in West and Central Africa (Côte d'Ivoire, Cameroon, Central African Republic and Democratic Republic of Congo). The isolates shared several phenotypic features of both B. anthracis and B. cereus, but differed amongst each other in motility and their resistance or sensitivity to penicillin. They all possessed the same mutation in the regulator gene plcR, different from the one found in B. anthracis, and in addition, carry genes which enable them to produce a second capsule composed of hyaluronic acid. Our findings show the existence of a discrete clade of the B. cereus group capable of causing anthrax-like disease, found in areas of high biodiversity, which are possibly also the origin of the worldwide distributed B. anthracis. Establishing the impact of these pathogenic bacteria on threatened wildlife species will require systematic investigation. Furthermore, the consumption of wildlife found dead by the local population and presence in a domestic animal reveal potential sources of exposure to humans.

A Study of the Infant Nasal Microbiome Development over the First Year of Life and in Relation to Their Primary Adult Caregivers Using cpn60 Universal Target (UT) as a Phylogenetic Marker.

Whereas the infant gut microbiome is the subject of intense study, relatively little is known regarding the nares microbiome in newborns and during early life. This study aimed to survey the typical composition and diversity of human anterior nare microflora for developing infants over time, and to explore how these correlate to their primary caregivers. Single nare swabs were collected at five time points over a one-year period for each subject from infant-caregiver pairs. Our study comprised of 50 infants (recruited at 2 weeks, post delivery) and their 50 primary caregivers. Applying the chaperonin-60 (cpn60) universal target (UT) amplicon as our molecular barcoding marker to census survey the microbial communities, we longitudinally surveyed infant nares microbiota at 5 time points over the course of the first year of life. The inter- and intra-subject diversity was catalogued and compared, both longitudinally and relative to their adult primary caregivers. Although within-subject variability over time and inter-subject variability were both observed, the assessment detected only one or two predominant genera for individual infant samples, belonging mainly to phyla Actinobacteria, Firmicutes, and Proteobacteria. Consistent with previously observed microbial population dynamics in other body sites, the diversity of nares microflora increased over the first year of life and infants showed differential operational taxonomic units (OTUs) relative to their matched primary caregiver. The collected evidence also support that both temporal and seasonal changes occur with respect to carriage of potentially pathogenic bacteria (PPBs), which may influence host predisposition to infection. This pilot study surveying paired infant/caregiver nare microbiomes provides novel longitudinal diversity information that is pertinent to better understanding nare microbiome development in infants.

Comparison of Sample Preparation Methods Used for the Next-Generation Sequencing of Mycobacterium tuberculosis.

The advent and widespread application of next-generation sequencing (NGS) technologies to the study of microbial genomes has led to a substantial increase in the number of studies in which whole genome sequencing (WGS) is applied to the analysis of microbial genomic epidemiology. However, microorganisms such as Mycobacterium tuberculosis (MTB) present unique problems for sequencing and downstream analysis based on their unique physiology and the composition of their genomes. In this study, we compare the quality of sequence data generated using the Nextera and TruSeq isolate preparation kits for library construction prior to Illumina sequencing-by-synthesis. Our results confirm that MTB NGS data quality is highly dependent on the purity of the DNA sample submitted for sequencing and its guanine-cytosine content (or GC-content). Our data additionally demonstrate that the choice of library preparation method plays an important role in mitigating downstream sequencing quality issues. Importantly for MTB, the Illumina TruSeq library preparation kit produces more uniform data quality than the Nextera XT method, regardless of the quality of the input DNA. Furthermore, specific genomic sequence motifs are commonly missed by the Nextera XT method, as are regions of especially high GC-content relative to the rest of the MTB genome. As coverage bias is highly undesirable, this study illustrates the importance of appropriate protocol selection when performing NGS studies in order to ensure that sound inferences can be made regarding mycobacterial genomes.

Tatumella saanichensis sp. nov., isolated from a cystic fibrosis patient.

Polyphasic taxonomic analysis was performed on a clinical isolate (NML 06-3099T) from a cystic fibrosis patient, including whole-genome sequencing, proteomics, phenotypic testing, electron microscopy, chemotaxonomy and a clinical investigation. Comparative whole-genome sequence analysis and multilocus sequence analysis (MLSA) between Tatumella ptyseos ATCC 33301T and clinical isolate NML 06-3099T suggested that the clinical isolate was closely related to, but distinct from, the species T. ptyseos. By 16S rRNA gene sequencing, the clinical isolate shared 98.7 % sequence identity with T. ptyseos ATCC 33301T. A concatenate of six MLSA loci (totalling 4500 bp) revealed < 93.9 % identity between T. ptyseos ATCC 33301T, other members of the genus and the clinical isolate. A whole-genome sequence comparison between NML 06-3099T and ATCC 33301T determined that the average nucleotide identity was 76.24 %. The overall DNA G+C content of NML 06-3099T was 51.27 %, consistent with members of the genus Tatumella. By matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) MS analysis, NML 06-3099T had a genus-level match, but not a species-level match, to T. ptyseos. By shotgun proteomics, T. ptyseos ATCC 33301T and NML 06-3099T were found to have unique proteomes. The two strains had similar morphologies and multiple fimbriae, as observed by transmission electron microscopy, but were distinguishable by phenotypic testing. Cellular fatty acids found were typical for members of the Enterobacteriaceae. NML 06-3099T was susceptible to commonly used antibiotics. Based on these data, NML 06-3099T represents a novel species in the genus Tatumella, for which the name Tatumella saanichensis sp. nov. is proposed (type strain NML 06-3099T = CCUG 55408T = DSM 19846T).

Sylvatic plague in a Canadian black-tailed prairie dog (Cynomys ludovicianus).

In 2010, a black-tailed prairie dog (Cynomys ludovicianus) was found dead in Grasslands National Park, Saskatchewan, Canada. Postmortem gross and histologic findings indicated bacterial septicemia, likely due to Yersinia pestis, which was confirmed by molecular analysis. This is the first report of Y. pestis in the prairie dog population within Canada.

DNA sequence heterogeneity of Campylobacter jejuni CJIE4 prophages and expression of prophage genes.

Campylobacter jejuni carry temperate bacteriophages that can affect the biology or virulence of the host bacterium. Known effects include genomic rearrangements and resistance to DNA transformation. C. jejuni prophage CJIE1 shows sequence variability and variability in the content of morons. Homologs of the CJIE1 prophage enhance both adherence and invasion to cells in culture and increase the expression of a specific subset of bacterial genes. Other C. jejuni temperate phages have so far not been well characterized. In this study we describe investigations into the DNA sequence variability and protein expression in a second prophage, CJIE4. CJIE4 sequences were obtained de novo from DNA sequencing of five C. jejuni isolates, as well as from whole genome sequences submitted to GenBank by other research groups. These CJIE4 DNA sequences were heterogenous, with several different insertions/deletions (indels) in different parts of the prophage genome. Two variants of a 3-4 kb region inserted within CJIE4 had different gene content that distinguished two major conserved CJIE4 prophage families. Additional indels were detected throughout the prophage. Detection of proteins in the five isolates characterized in our laboratory in isobaric Tags for Relative and Absolute Quantitation (iTRAQ) experiments indicated that prophage proteins within each of the two large indel variants were expressed during growth of the bacteria on Mueller Hinton agar plates. These proteins included the extracellular DNase associated with resistance to DNA transformation and prophage repressor proteins. Other proteins associated with known or suspected roles in prophage biology were also expressed from CJIE4, including capsid protein, the phage integrase, and MazF, a type II toxin-antitoxin system protein. Together with the results previously obtained for the CJIE1 prophage these results demonstrate that sequence variability and expression of moron genes are both general properties of temperate bacteriophages in C. jejuni.

Phylogenetic relationship and virulence inference of Streptococcus Anginosus Group: curated annotation and whole-genome comparative analysis support distinct species designation.

BACKGROUND: The Streptococcus Anginosus Group (SAG) represents three closely related species of the viridans group streptococci recognized as commensal bacteria of the oral, gastrointestinal and urogenital tracts. The SAG also cause severe invasive infections, and are pathogens during cystic fibrosis (CF) pulmonary exacerbation. Little genomic information or description of virulence mechanisms is currently available for SAG. We conducted intra and inter species whole-genome comparative analyses with 59 publically available Streptococcus genomes and seven in-house closed high quality finished SAG genomes; S. constellatus (3), S. intermedius (2), and S. anginosus (2). For each SAG species, we sequenced at least one numerically dominant strain from CF airways recovered during acute exacerbation and an invasive, non-lung isolate. We also evaluated microevolution that occurred within two isolates that were cultured from one individual one year apart. RESULTS: The SAG genomes were most closely related to S. gordonii and S. sanguinis, based on shared orthologs and harbor a similar number of proteins within each COG category as other Streptococcus species. Numerous characterized streptococcus virulence factor homologs were identified within the SAG genomes including; adherence, invasion, spreading factors, LPxTG cell wall proteins, and two component histidine kinases known to be involved in virulence gene regulation. Mobile elements, primarily integrative conjugative elements and bacteriophage, account for greater than 10% of the SAG genomes. S. anginosus was the most variable species sequenced in this study, yielding both the smallest and the largest SAG genomes containing multiple genomic rearrangements, insertions and deletions. In contrast, within the S. constellatus and S. intermedius species, there was extensive continuous synteny, with only slight differences in genome size between strains. Within S. constellatus we were able to determine important SNPs and changes in VNTR numbers that occurred over the course of one year. CONCLUSIONS: The comparative genomic analysis of the SAG clarifies the phylogenetics of these bacteria and supports the distinct species classification. Numerous potential virulence determinants were identified and provide a foundation for further studies into SAG pathogenesis. Furthermore, the data may be used to enable the development of rapid diagnostic assays and therapeutics for these pathogens.

whole-genome sequence of livestock-associated st398 methicillin-resistant staphylococcus aureus Isolated from Humans in Canada.

Despite reports of high colonization rates of ST398 livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) among pigs and pig farmers, the incidence of LA-MRSA infection in the general population in Canada appears to be rare in comparison to that in some European countries. In this study, the complete genome sequence of a Canadian representative LA-MRSA isolate (08BA02176) from a human postoperative surgical site infection was acquired and compared to the sequenced genome of an LA-MRSA isolate (S0385) from Europe to identify genetic traits that may explain differences in the success of these particular strains in some locales.

Comparative genomics of Vibrio cholerae from Haiti, Asia, and Africa.

Cholera was absent from the island of Hispaniola at least a century before an outbreak that began in Haiti in the fall of 2010. Pulsed-field gel electrophoresis (PFGE) analysis of clinical isolates from the Haiti outbreak and recent global travelers returning to the United States showed indistinguishable PFGE fingerprints. To better explore the genetic ancestry of the Haiti outbreak strain, we acquired 23 whole-genome Vibrio cholerae sequences: 9 isolates obtained in Haiti or the Dominican Republic, 12 PFGE pattern-matched isolates linked to Asia or Africa, and 2 nonmatched outliers from the Western Hemisphere. Phylogenies for whole-genome sequences and core genome single-nucleotide polymorphisms showed that the Haiti outbreak strain is genetically related to strains originating in India and Cameroon. However, because no identical genetic match was found among sequenced contemporary isolates, a definitive genetic origin for the outbreak in Haiti remains speculative.