Pimlapas Leekitcharoenphon

The transcriptional landscape and small RNAs of Salmonella enterica serovar Typhimurium

More than 50 y of research have provided great insight into the physiology, metabolism, and molecular biology of Salmonella enterica serovar Typhimurium (S. Typhimurium), but important gaps in our knowledge remain. It is clear that a precise choreography of gene expression is required for Salmonella infection, but basic genetic information such as the global locations of transcription start sites (TSSs) has been lacking. We combined three RNA-sequencing techniques and two sequencing platforms to generate a robust picture of transcription in S. Typhimurium. Differential RNA sequencing identified 1,873 TSSs on the chromosome of S. Typhimurium SL1344 and 13% of these TSSs initiated antisense transcripts. Unique findings include the TSSs of the virulence regulators phoP, slyA, and invF. Chromatin immunoprecipitation revealed that RNA polymerase was bound to 70% of the TSSs, and two-thirds of these TSSs were associated with σ70 (including phoP, slyA, and invF) from which we identified the −10 and −35 motifs of σ70-dependent S. Typhimurium gene promoters. Overall, we corrected the location of important genes and discovered 18 times more promoters than identified previously. S. Typhimurium expresses 140 small regulatory RNAs (sRNAs) at early stationary phase, including 60 newly identified sRNAs. Almost half of the experimentally verified sRNAs were found to be unique to the Salmonella genus, and <20% were found throughout the Enterobacteriaceae. This description of the transcriptional map of SL1344 advances our understanding of S. Typhimurium, arguably the most important bacterial infection model.

Detection of mcr-1 encoding plasmid-mediated colistin-resistant Escherichia coli isolates from human bloodstream infection and imported chicken meat, Denmark 2015.

The plasmid-mediated colistin resistance gene, mcr-1, was detected in an Escherichia coli isolate from a Danish patient with bloodstream infection and in five E. coli isolates from imported chicken meat. One isolate from chicken meat belonged to the epidemic spreading sequence type ST131. In addition to IncI2, an incX4 replicon was found to be linked to mcr-1. This report follows a recent detection of mcr-1 in E. coli from animals, food and humans in China.

Solving the problem of comparing whole bacterial genomes across different sequencing platforms.

Whole genome sequencing (WGS) shows great potential for real-time monitoring and identification of infectious disease outbreaks. However, rapid and reliable comparison of data generated in multiple laboratories and using multiple technologies is essential. So far studies have focused on using one technology because each technology has a systematic bias making integration of data generated from different platforms difficult. We developed two different procedures for identifying variable sites and inferring phylogenies in WGS data across multiple platforms. The methods were evaluated on three bacterial data sets and sequenced on three different platforms (Illumina, 454, Ion Torrent). We show that the methods are able to overcome the systematic biases caused by the sequencers and infer the expected phylogenies. It is concluded that the cause of the success of these new procedures is due to a validation of all informative sites that are included in the analysis. The procedures are available as web tools.

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 reveling 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.

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

BackgroundThe 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.ResultsHere 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.The server was evaluated using four published bacterial WGS data sets (V. cholerae, S. aureus CC398, S. Typhimurium and M. tuberculosis). 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.ConclusionsThe 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/.

Genomic variation in Salmonella enterica core genes for epidemiological typing

BackgroundTechnological 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.ResultsWe identify a set of 2,882 core genes clusters based on 73 publicly available Salmonella enterica 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.ConclusionsGenomic 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.

Global Genomic Epidemiology of Salmonella enterica Serovar Typhimurium DT104

ABSTRACT It has been 30 years since the initial emergence and subsequent rapid global spread of multidrug-resistant Salmonella enterica serovar Typhimurium DT104 (MDR DT104). Nonetheless, its origin and transmission route have never been revealed. We used whole-genome sequencing (WGS) and temporally structured sequence analysis within a Bayesian framework to reconstruct temporal and spatial phylogenetic trees and estimate the rates of mutation and divergence times of 315 S. Typhimurium DT104 isolates sampled from 1969 to 2012 from 21 countries on six continents. DT104 was estimated to have emerged initially as antimicrobial susceptible in ∼1948 (95% credible interval [CI], 1934 to 1962) and later became MDR DT104 in ∼1972 (95% CI, 1972 to 1988) through horizontal transfer of the 13-kb Salmonella genomic island 1 (SGI1) MDR region into susceptible strains already containing SGI1. This was followed by multiple transmission events, initially from central Europe and later between several European countries. An independent transmission to the United States and another to Japan occurred, and from there MDR DT104 was probably transmitted to Taiwan and Canada. An independent acquisition of resistance genes took place in Thailand in ∼1975 (95% CI, 1975 to 1990). In Denmark, WGS analysis provided evidence for transmission of the organism between herds of animals. Interestingly, the demographic history of Danish MDR DT104 provided evidence for the success of the program to eradicate Salmonella from pig herds in Denmark from 1996 to 2000. The results from this study refute several hypotheses on the evolution of DT104 and suggest that WGS may be useful in monitoring emerging clones and devising strategies for prevention of Salmonella infections.

Genomic Signature of Multidrug-Resistant Salmonella enterica Serovar Typhi Isolates Related to a Massive Outbreak in Zambia between 2010 and 2012

ABSTRACT Retrospectively, we investigated the epidemiology of a massive Salmonella enterica serovar Typhi outbreak in Zambia during 2010 to 2012. Ninety-four isolates were susceptibility tested by MIC determinations. Whole-genome sequence typing (WGST) of 33 isolates and bioinformatic analysis identified the multilocus sequence type (MLST), haplotype, plasmid replicon, antimicrobial resistance genes, and genetic relatedness by single nucleotide polymorphism (SNP) analysis and genomic deletions. The outbreak affected 2,040 patients, with a fatality rate of 0.5%. Most (83.0%) isolates were multidrug resistant (MDR). The isolates belonged to MLST ST1 and a new variant of the haplotype, H58B. Most isolates contained a chromosomally translocated region containing seven antimicrobial resistance genes, catA1, bla TEM-1, dfrA7, sul1, sul2, strA, and strB, and fragments of the incompatibility group Q1 (IncQ1) plasmid replicon, the class 1 integron, and the mer operon. The genomic analysis revealed 415 SNP differences overall and 35 deletions among 33 of the isolates subjected to whole-genome sequencing. In comparison with other genomes of H58, the Zambian isolates separated from genomes from Central Africa and India by 34 and 52 SNPs, respectively. The phylogenetic analysis indicates that 32 of the 33 isolates sequenced belonged to a tight clonal group distinct from other H58 genomes included in the study. The small numbers of SNPs identified within this group are consistent with the short-term transmission that can be expected over a period of 2 years. The phylogenetic analysis and deletions suggest that a single MDR clone was responsible for the outbreak, during which occasional other S. Typhi lineages, including sensitive ones, continued to cocirculate. The common view is that the emerging global S. Typhi haplotype, H58B, containing the MDR IncHI1 plasmid is responsible for the majority of typhoid infections in Asia and sub-Saharan Africa; we found that a new variant of the haplotype harboring a chromosomally translocated region containing the MDR islands of IncHI1 plasmid has emerged in Zambia. This could change the perception of the term “classical MDR typhoid” currently being solely associated with the IncHI1 plasmid. It might be more common than presently thought that S. Typhi haplotype H58B harbors the IncHI1 plasmid or a chromosomally translocated MDR region or both.

The Role of the st313-td Gene in Virulence of Salmonella Typhimurium ST313

Multidrug-resistant Salmonella enterica serovar Typhimurium ST313 has emerged in sub-Saharan Africa causing severe infections in humans. Therefore, it has been speculated that this specific sequence type, ST313, carries factors associated with increased pathogenicity. We assessed the role in virulence of a gene with a yet unknown function, st313-td, detected in ST313 through comparative genomics. Additionally, the structure of the genomic island ST313-GI, harbouring the gene was determined. The gene st313-td was cloned into wild type S. Typhimurium 4/74 (4/74-C) as well as knocked out in S. Typhimurium ST313 02–03/002 (Δst313-td) followed by complementation (02-03/002-C). Δst313-td was less virulent in mice following i.p. challenge than the wild type and this phenotype could be partly complemented in trans, indicating that st313-td plays a role during systemic infection. The gene st313-td was shown not to affect invasion of cultured epithelial cells, while the absence of the gene significantly affects uptake and intracellular survival within macrophages. The gene st313-td was proven to be strongly associated to invasiveness, harboured by 92.5% of S. Typhimurium blood isolates (n = 82) and 100% of S. Dublin strains (n = 50) analysed. On the contrary, S. Typhimurium isolates of animal and food origin (n = 82) did not carry st313-td. Six human, non-blood isolates of S. Typhimurium from Belarus, China and Nepal harboured the gene and belonged to sequence types ST398 and ST19. Our data showed a global presence of the st313-td gene and in other sequence types than ST313. The gene st313-td was shown to be expressed during logarithmic phase of growth in 14 selected Salmonella strains carrying the gene. This study reveals that st313-td plays a role in S. Typhimurium ST313 pathogenesis and adds another chapter to understanding of the virulence of S. Typhimurium and in particular of the emerging sequence type ST313.

Extremely Drug-Resistant Salmonella enterica Serovar Senftenberg Infections in Patients in Zambia

ABSTRACT Two cases of extremely drug-resistant Salmonella enterica serovar Senftenberg isolated from patients in Zambia were investigated by utilizing MIC determinations and whole-genome sequencing. The isolates were resistant to, and harbored genes toward, nine drug classes, including fluoroquinolones and extended-spectrum cephalosporins, contained two plasmid replicons, and differed by 93 single-nucleotide polymorphisms.