Epidemiology of Danish Aeromonas salmonicida subsp salmonicida in Fish Farms Using Whole Genome Sequencing

Furunculosis, a serious infection caused by the bacterium Aeromonas salmonicida subsp. salmonicida is common in sea-reared rainbow trout production in Denmark. Developing an effective control strategy requires knowledge of the epidemiology, as well as the genomic and virulent variability of the Danish A. salmonicida subsp. salmonicida isolates. To obtain this, the genomes of 101 A. salmonicida subsp. salmonicida, including 99 Danish isolates, one Scottish strain and the type strain NCIMB 1102, were sequenced using the Illumina HiSeq platform. Isolates were de novo assembled, examined for presence of plasmids, virulence and iron acquisition proteins, genomic islands, and antibiotic resistance genes. Single Nucleotide Polymorphisms were aligned and subjected to Bayesian temporal phylogenetic and maximum likelihood tree reconstruction using the published genome of A. salmonicida subsp. salmonicida A449 as reference. Bayesian temporal phylogenetic reconstruction suggests that four major introductions of A. salmonicida subsp. salmonicida into Denmark have occurred. The introductions correlate with the freshwater and subsequent seawater expansion of rainbow trout production. Initial transmission of the bacterium could have been from seawater to freshwater or vice versa, and most minor clades include a mixture of strains from different fresh- and seawater farms. Genomic variation of A. salmonicida subsp. salmonicida mostly appeared to be associated with their plasmids and plasmid encoded virulence factors. Nine A. salmonicida subsp. salmonicida isolates harbored worldwide known antibiotic resistance genes against several antibiotics and there is an indication that 33% of the isolates contained the genomic island AsaGEI1b. These findings not only support the usefulness of whole genome sequencing for genetic studies of homogeneous bacteria in general, but provide novel information about the Danish A. salmonicida subsp. salmonicida population, with implications for vaccine development in efforts to better protect Danish rainbow trout in the future

Genomic variation in Salmonella enterica core genes for epidemiological typing

<p>Abstract</p> <p>Background</p> <p>Technological advances in high throughput genome sequencing are making whole genome sequencing (WGS) available as a routine tool for bacterial typing. Standardized procedures for identification of relevant genes and of variation are needed to enable comparison between studies and over time. The core genes--the genes that are conserved in all (or most) members of a genus or species--are potentially good candidates for investigating genomic variation in phylogeny and epidemiology.</p> <p>Results</p> <p>We identify a set of 2,882 core genes clusters based on 73 publicly available <it>Salmonella enterica </it>genomes and evaluate their value as typing targets, comparing whole genome typing and traditional methods such as 16S and MLST. A consensus tree based on variation of core genes gives much better resolution than 16S and MLST; the pan-genome family tree is similar to the consensus tree, but with higher confidence. The core genes can be divided into two categories: a few highly variable genes and a larger set of conserved core genes, with low variance. For the most variable core genes, the variance in amino acid sequences is higher than for the corresponding nucleotide sequences, suggesting that there is a positive selection towards mutations leading to amino acid changes.</p> <p>Conclusions</p> <p>Genomic variation within the core genome is useful for investigating molecular evolution and providing candidate genes for bacterial genome typing. Identification of genes with different degrees of variation is important especially in trend analysis.</p

snpTree - a web-server to identify and construct SNP trees from whole genome sequence data.

<p>Abstract</p> <p>Background</p> <p>The advances and decreasing economical cost of whole genome sequencing (WGS), will soon make this technology available for routine infectious disease epidemiology. In epidemiological studies, outbreak isolates have very little diversity and require extensive genomic analysis to differentiate and classify isolates. One of the successfully and broadly used methods is analysis of single nucletide polymorphisms (SNPs). Currently, there are different tools and methods to identify SNPs including various options and cut-off values. Furthermore, all current methods require bioinformatic skills. Thus, we lack a standard and simple automatic tool to determine SNPs and construct phylogenetic tree from WGS data.</p> <p>Results</p> <p>Here we introduce snpTree, a server for online-automatic SNPs analysis. This tool is composed of different SNPs analysis suites, perl and python scripts. snpTree can identify SNPs and construct phylogenetic trees from WGS as well as from assembled genomes or contigs. WGS data in fastq format are aligned to reference genomes by BWA while contigs in fasta format are processed by Nucmer. SNPs are concatenated based on position on reference genome and a tree is constructed from concatenated SNPs using FastTree and a perl script. The online server was implemented by HTML, Java and python script.</p> <p>The server was evaluated using four published bacterial WGS data sets (<it>V. cholerae</it>, <it>S. aureus </it>CC398, <it>S</it>. Typhimurium and <it>M. tuberculosis</it>). The evalution results for the first three cases was consistent and concordant for both raw reads and assembled genomes. In the latter case the original publication involved extensive filtering of SNPs, which could not be repeated using snpTree.</p> <p>Conclusions</p> <p>The snpTree server is an easy to use option for rapid standardised and automatic SNP analysis in epidemiological studies also for users with limited bioinformatic experience. The web server is freely accessible at http://www.cbs.dtu.dk/services/snpTree-1.0/.</p

Evaluation of whole genome sequencing for outbreak detection of Salmonella enterica

Salmonella enterica is a common cause of minor and large food borne outbreaks. To achieve successful and nearly 'real-time' monitoring and identification of outbreaks, reliable sub-typing is essential. Whole genome sequencing (WGS) shows great promises for using as a routine epidemiological typing tool. Here we evaluate WGS for typing of S. Typhimurium including different approaches for analyzing and comparing the data. A collection of 34 S. Typhimurium isolates was sequenced. This consisted of 18 isolates from six outbreaks and 16 epidemiologically unrelated background strains. In addition, 8 S. Enteritidis and 5 S. Derby were also sequenced and used for comparison. A number of different bioinformatics approaches were applied on the data; including pan-genome tree, k-mer tree, nucleotide difference tree and SNP tree. The outcome of each approach was evaluated in relation to the association of the isolates to specific outbreaks. The pan-genome tree clustered 65% of the S. Typhimurium isolates according to the pre-defined epidemiology, the k-mer tree 88%, the nucleotide difference tree 100% and the SNP tree 100% of the strains within S. Typhimurium. The resulting outcome of the four phylogenetic analyses were also compared to PFGE revealing that WGS typing achieved the greater performance than the traditional method. In conclusion, for S. Typhimurium, SNP analysis and nucleotide difference approach of WGS data seem to be the superior methods for epidemiological typing compared to other phylogenetic analytic approaches that may be used on WGS. These approaches were also superior to the more classical typing method, PFGE. Our study also indicates that WGS alone is insufficient to determine whether strains are related or un-related to outbreaks. This still requires the combination of epidemiological data and whole genome sequencing results

Characterization and genetic variation of vibrio cholerae isolated from clinical and environmental sources in Thailand

Cholera is still an important public health problem in several countries, including Thailand. In this study, a collection of clinical and environmental V. cholerae serogroup O1, O139, and non-O1/non-O139 strains originating from Thailand (1983 to 2013) was characterized to determine phenotypic and genotypic traits and to investigate the genetic relatedness. Using a combination of conventional methods and whole genome sequencing (WGS), 78 V. cholerae strains were identified. WGS was used to determine the serogroup, biotype, virulence, mobile genetic elements, and antimicrobial resistance genes using online bioinformatics tools. In addition, phenotypic antimicrobial resistance was determined by the minimal inhibitory concentration (MIC) test. The 78 V. cholerae strains belonged to the following serogroups O1: (n = 44), O139 (n = 16) and non-O1/non-O139 (n = 18). Interestingly, we found that the typical El Tor O1 strains were the major cause of clinical cholera during 1983-2000 with two Classical O1 strains detected in 2000. In 2004-2010, the El Tor variant strains revealed genotypes of the Classical biotype possessing either only ctxB or both ctxB and rstR while they harbored tcpA of the El Tor biotype. Thirty O1 and eleven O139 clinical strains carried CTXϕ (Cholera toxin) and tcpA as well four different pathogenic islands (PAIs). Beside non-O1/non-O139, the O1 environmental strains also presented chxA and Type Three Secretion System (TTSS). The in silico MultiLocus Sequence Typing (MLST) discriminated the O1 and O139 clinical strains from other serogroups and environmental strains. ST69 was dominant in the clinical strains belonging to the 7th pandemic clone. Non-O1/non-O139 and environmental strains showed various novel STs indicating genetic variation. Multidrug-resistant (MDR) strains were observed and conferred resistance to ampicillin, azithromycin, nalidixic acid, sulfamethoxazole, tetracycline, and trimethoprim and harboured variants of the SXT elements. For the first time since 1986, the presence of V. cholerae O1 Classical was reported causing cholera outbreaks in Thailand. In addition, we found that V. cholerae O1 El Tor variant and O139 were pre-dominating the pathogenic strains in Thailand. Using WGS and bioinformatic tools to analyze both historical and contemporary V. cholerae circulating in Thailand provided a more detailed understanding of the V. cholerae epidemiology, which ultimately could be applied for control measures and management of cholera in Thailand