scholarly journals Whole-genome epidemiology links phage-mediated acquisition of a virulence gene to the clonal expansion of a pandemic Salmonella enterica serovar Typhimurium clone

2020 ◽  
Vol 6 (11) ◽  
Author(s):  
Eleonora Tassinari ◽  
Matt Bawn ◽  
Gaetan Thilliez ◽  
Oliver Charity ◽  
Luke Acton ◽  
...  

Epidemic and pandemic clones of bacterial pathogens with distinct characteristics continually emerge, replacing those previously dominant through mechanisms that remain poorly characterized. Here, whole-genome-sequencing-powered epidemiology linked horizontal transfer of a virulence gene, sopE, to the emergence and clonal expansion of a new epidemic Salmonella enterica serovar Typhimurium (S. Typhimurium) clone. The sopE gene is sporadically distributed within the genus Salmonella and rare in S . enterica Typhimurium lineages, but was acquired multiple times during clonal expansion of the currently dominant pandemic monophasic S. Typhimurium sequence type (ST) 34 clone. Ancestral state reconstruction and time-scaled phylogenetic analysis indicated that sopE was not present in the common ancestor of the epidemic clade, but later acquisition resulted in increased clonal expansion of sopE-containing clones that was temporally associated with emergence of the epidemic, consistent with increased fitness. The sopE gene was mainly associated with a temperate bacteriophage mTmV, but recombination with other bacteriophage and apparent horizontal gene transfer of the sopE gene cassette resulted in distribution among at least four mobile genetic elements within the monophasic S . enterica Typhimurium ST34 epidemic clade. The mTmV prophage lysogenic transfer to other S. enterica serovars in vitro was limited, but included the common pig-associated S . enterica Derby (S. Derby). This may explain mTmV in S. Derby co-circulating on farms with monophasic S. Typhimurium ST34, highlighting the potential for further transfer of the sopE virulence gene in nature. We conclude that whole-genome epidemiology pinpoints potential drivers of evolutionary and epidemiological dynamics during pathogen emergence, and identifies targets for subsequent research in epidemiology and bacterial pathogenesis.

2021 ◽  
Vol 7 (12) ◽  
Author(s):  
David J. Edwards ◽  
Sebastián Duchene ◽  
Bernard Pope ◽  
Kathryn E. Holt

Homoplasic SNPs are considered important signatures of strong (positive) selective pressure, and hence of adaptive evolution for clinically relevant traits such as antibiotic resistance and virulence. Here we present a new tool, SNPPar, for efficient detection and analysis of homoplasic SNPs from large whole genome sequencing datasets (>1000 isolates and/or >100 000 SNPs). SNPPar takes as input an SNP alignment, tree and annotated reference genome, and uses a combination of simple monophyly tests and ancestral state reconstruction (ASR, via TreeTime) to assign mutation events to branches and identify homoplasies. Mutations are annotated at the level of codon and gene, to facilitate analysis of convergent evolution. Testing on simulated data (120 Mycobacterium tuberculosis alignments representing local and global samples) showed SNPPar can detect homoplasic SNPs with very high specificity (zero false-positives in all tests) and high sensitivity (zero false-negatives in 89 % of tests). SNPPar analysis of three empirically sampled datasets ( Elizabethkingia anophelis , Burkholderia dolosa and M. tuberculosis ) produced results that were in concordance with previous studies, in terms of both individual homoplasies and evidence of convergence at the codon and gene levels. SNPPar analysis of a simulated alignment of ~64 000 genome-wide SNPs from 2000 M. tuberculosis genomes took ~23 min and ~2.6 GB of RAM to generate complete annotated results on a laptop. This analysis required ASR be conducted for only 1.25 % of SNPs, and the ASR step took ~23 s and 0.4 GB of RAM. SNPPar automates the detection and annotation of homoplasic SNPs efficiently and accurately from large SNP alignments. As demonstrated by the examples included here, this information can be readily used to explore the role of homoplasy in parallel and/or convergent evolution at the level of nucleotide, codon and/or gene.


2021 ◽  
Vol 70 (12) ◽  
Author(s):  
Letícia T. Oliveira ◽  
Lívia A. Alves ◽  
Erika N. Harth-Chu ◽  
Ryota Nomura ◽  
Kazuhiko Nakano ◽  
...  

Introduction. Streptococcus mutans , a common species of the oral microbiome, expresses virulence genes promoting cariogenic dental biofilms, persistence in the bloodstream and cardiovascular infections. Gap statement. Virulence gene expression is variable among S. mutans strains and controlled by the transcription regulatory systems VicRK and CovR. Aim. This study investigates polymorphisms in the vicRK and covR loci in S. mutans strains isolated from the oral cavity or from the bloodstream, which were shown to differ in expression of covR, vicRK and downstream genes. Methodology. The transcriptional activities of covR, vicR and vicK were compared by RT-qPCR between blood and oral strains after exposure to human serum. PCR-amplified promoter and/or coding regions of covR and vicRK of 18 strains (11 oral and 7 blood) were sequenced and compared to the reference strain UA159. Results. Serum exposure significantly reduced covR and vicR/K transcript levels in most strains (P<0.05), but reductions were higher in oral than in blood strains. Single-nucleotide polymorphisms (SNPs) were detected in covR regulatory and coding regions, but SNPs affecting the CovR effector domain were only present in two blood strains. Although vicR was highly conserved, vicK showed several SNPs, and SNPs affecting VicK regions important for autokinase activity were found in three blood strains. Conclusions. This study reveals transcriptional and structural diversity in covR and vicR/K, and identifies polymorphisms of functional relevance in blood strains, indicating that covR and vicRK might be important loci for S. mutans adaptation to host selective pressures associated with virulence diversity.


2017 ◽  
Vol 5 (46) ◽  
Author(s):  
Najwa Syahirah Roslan ◽  
Shagufta Jabeen ◽  
Nurulfiza Mat Isa ◽  
Abdul Rahman Omar ◽  
Mohd Hair Bejo ◽  
...  

ABSTRACT Salmonella enterica subsp. enterica serovar Typhimurium is one of several well-categorized Salmonella serotypes recognized globally. Here, we report the whole-genome sequence of S. Typhimurium strain UPM 260, isolated from a broiler chicken.


Author(s):  
Hisami Kobayashi ◽  
Yasuhiro Tanizawa ◽  
Mitsuo Sakamoto ◽  
Moriya Ohkuma ◽  
Masanori Tohno

The taxonomic status of the species Clostridium methoxybenzovorans was assessed. The 16S rRNA gene sequence, whole-genome sequence and phenotypic characterizations suggested that the type strain deposited in the American Type Culture Collection ( C. methoxybenzovorans ATCC 700855T) is a member of the species Eubacterium callanderi . Hence, C. methoxybenzovorans ATCC 700855T cannot be used as a reference for taxonomic study. The type strain deposited in the German Collection of Microorganism and Cell Cultures GmbH (DSM 12182T) is no longer listed in its online catalogue. Also, both the 16S rRNA gene and the whole-genome sequences of the original strain SR3T showed high sequence identity with those of Lacrimispora indolis (recently reclassified from Clostridium indolis ) as the most closely related species. Analysis of the two genomes showed average nucleotide identity based on blast and digital DNA–DNA hybridization values of 98.3 and 87.9 %, respectively. Based on these results, C. methoxybenzovorans SR3T was considered to be a member of L. indolis .


2013 ◽  
Vol 63 (Pt_10) ◽  
pp. 3920-3926 ◽  
Author(s):  
Julia S. Bennett ◽  
Keith A. Jolley ◽  
Martin C. J. Maiden

Phylogenies generated from whole genome sequence (WGS) data provide definitive means of bacterial isolate characterization for typing and taxonomy. The species status of strains recently defined with conventional taxonomic approaches as representing Neisseria oralis was examined by the analysis of sequences derived from WGS data, specifically: (i) 53 Neisseria ribosomal protein subunit (rps) genes (ribosomal multi-locus sequence typing, rMLST); and (ii) 246 Neisseria core genes (core genome MLST, cgMLST). These data were compared with phylogenies derived from 16S and 23S rRNA gene sequences, demonstrating that the N. oralis strains were monophyletic with strains described previously as representing ‘ Neisseria mucosa var. heidelbergensis’ and that this group was of equivalent taxonomic status to other well-described species of the genus Neisseria . Phylogenetic analyses also indicated that Neisseria sicca and Neisseria macacae should be considered the same species as Neisseria mucosa and that Neisseria flavescens should be considered the same species as Neisseria subflava . Analyses using rMLST showed that some strains currently defined as belonging to the genus Neisseria were more closely related to species belonging to other genera within the family; however, whole genome analysis of a more comprehensive selection of strains from within the family Neisseriaceae would be necessary to confirm this. We suggest that strains previously identified as representing ‘ N. mucosa var. heidelbergensis’ and deposited in culture collections should be renamed N. oralis . Finally, one of the strains of N. oralis was able to ferment lactose, due to the presence of β-galactosidase and lactose permease genes, a characteristic previously thought to be unique to Neisseria lactamica , which therefore cannot be thought of as diagnostic for this species; however, the rMLST and cgMLST analyses confirm that N. oralis is most closely related to N. mucosa .


2021 ◽  
Vol 7 (7) ◽  
Author(s):  
Casper Jamin ◽  
Sien De Koster ◽  
Stefanie van Koeveringe ◽  
Dieter De Coninck ◽  
Klaas Mensaert ◽  
...  

Whole-genome sequencing (WGS) is becoming the de facto standard for bacterial typing and outbreak surveillance of resistant bacterial pathogens. However, interoperability for WGS of bacterial outbreaks is poorly understood. We hypothesized that harmonization of WGS for outbreak surveillance is achievable through the use of identical protocols for both data generation and data analysis. A set of 30 bacterial isolates, comprising of various species belonging to the Enterobacteriaceae family and Enterococcus genera, were selected and sequenced using the same protocol on the Illumina MiSeq platform in each individual centre. All generated sequencing data were analysed by one centre using BioNumerics (6.7.3) for (i) genotyping origin of replications and antimicrobial resistance genes, (ii) core-genome multi-locus sequence typing (cgMLST) for Escherichia coli and Klebsiella pneumoniae and whole-genome multi-locus sequencing typing (wgMLST) for all species. Additionally, a split k-mer analysis was performed to determine the number of SNPs between samples. A precision of 99.0% and an accuracy of 99.2% was achieved for genotyping. Based on cgMLST, a discrepant allele was called only in 2/27 and 3/15 comparisons between two genomes, for E. coli and K. pneumoniae, respectively. Based on wgMLST, the number of discrepant alleles ranged from 0 to 7 (average 1.6). For SNPs, this ranged from 0 to 11 SNPs (average 3.4). Furthermore, we demonstrate that using different de novo assemblers to analyse the same dataset introduces up to 150 SNPs, which surpasses most thresholds for bacterial outbreaks. This shows the importance of harmonization of data-processing surveillance of bacterial outbreaks. In summary, multi-centre WGS for bacterial surveillance is achievable, but only if protocols are harmonized.


Microbiology ◽  
2020 ◽  
Vol 166 (8) ◽  
pp. 785-793
Author(s):  
Shou Miura ◽  
Yukino Tamamura ◽  
Mariko Takayasu ◽  
Miwa Sasaki ◽  
Natsuko Nishimura ◽  
...  

Salmonella enterica subspecies enterica serovar Typhimurium (S. Typhimurium) definitive phage type 104 (DT104), S. enterica subspecies enterica serovar Worthington (S. Worthington) and S. bongori produce ArtA and ArtB (ArtAB) toxin homologues, which catalyse ADP-ribosylation of pertussis toxin-sensitive G protein. ArtAB gene (artAB) is encoded on prophage in DT104 and its expression is induced by mitomycin C (MTC) and hydrogen peroxide (H2O2) that trigger the bacterial SOS response. Although the genetic regulatory mechanism associated with artAB expression is not characterized, it is thought to be associated with prophage induction, which occurs when the RecA-mediated SOS response is triggered. Here we show that subinhibitory concentration of quinolone antibiotics that are SOS-inducing agents, also induce ArtAB production in these Salmonella strains. Both MTC and fluoroquinolone antibiotics such as enrofloxacin-induced artA and recA transcription and artAB-encoding prophage (ArtAB-prophage) in DT104 and S. Worthington. However, in S. bongori , which harbours artAB genes on incomplete prophage, artA transcription was induced by MTC and enrofloxacin, but prophage induction was not observed. Taken together, these results suggest that SOS response followed by induction of artAB transcription is essential for ArtAB production. H2O2-mediated induction of ArtAB prophage and efficient production of ArtAB was observed in DT104 but not in S. Worthington and S. bongori . Therefore, induction of artAB expression with H2O2 is strain-specific, and the mode of action of H2O2 as an SOS-inducing agent might be different from those of MTC and quinolone antibiotics.


2021 ◽  
Vol 7 (11) ◽  
Author(s):  
Isabelle Bernaquez ◽  
Christiane Gaudreau ◽  
Pierre A. Pilon ◽  
Sadjia Bekal

Many public health laboratories across the world have implemented whole-genome sequencing (WGS) for the surveillance and outbreak detection of foodborne pathogens. PulseNet-affiliated laboratories have determined that most single-strain foodborne outbreaks are contained within 0–10 multi-locus sequence typing (MLST)-based allele differences and/or core genome single-nucleotide variants (SNVs). In addition to being a food- and travel-associated outbreak pathogen, most Shigella spp. cases occur through continuous person-to-person transmission, predominantly involving men who have sex with men (MSM), leading to long-term and recurrent outbreaks. Continuous transmission patterns coupled to genetic evolution under antibiotic treatment pressure require an assessment of existing WGS-based subtyping methods and interpretation criteria for cluster inclusion/exclusion. An evaluation of 4 WGS-based subtyping methods [SNVPhyl, coreMLST, core genome MLST (cgMLST) and whole-genome MLST (wgMLST)] was performed on 9 foodborne-, travel- and MSM-related retrospective outbreaks from a collection of 91 Shigella flexneri and 232  Shigella sonnei isolates to determine the methods’ epidemiological concordance, discriminatory power, robustness and ability to generate stable interpretation criteria. The discriminatory powers were ranked as follows: coreMLST<SNVPhyl<cgMLST<wgMLST (range: 0.970–1.000). The genetic differences observed for non-MSM-related Shigella spp. outbreaks respect the standard 0–10 allele/SNV guideline; however, mobile genetic element (MGE)-encoded loci caused inflated genetic variation and discrepant phylogenies for prolonged MSM-related S. sonnei outbreaks via wgMLST. The S. sonnei correlation coefficients of wgMLST were also the lowest at 0.680, 0.703 and 0.712 for SNVPhyl, coreMLST and cgMLST, respectively. Plasmid maintenance, mobilization and conjugation-associated genes were found to be the main source of genetic distance inflation in addition to prophage-related genes. Duplicated alleles arising from the repeated nature of IS elements were also responsible for many false cg/wgMLST differences. The coreMLST approach was shown to be the most robust, followed by SNVPhyl and wgMLST for inter-laboratory comparability. Our results highlight the need for validating species-specific subtyping methods based on microbial genome plasticity and outbreak dynamics in addition to the importance of filtering confounding MGEs for cluster detection.


2021 ◽  
Vol 7 (9) ◽  
Author(s):  
Sébastien O. Leclercq ◽  
Maxime Branger ◽  
David G. E. Smith ◽  
Pierre Germon

Escherichia coli is a very versatile species for which diversity has been explored from various perspectives highlighting, for example, phylogenetic groupings and pathovars, as well as a wide range of O serotypes. The highly variable O-antigen, the most external part of the lipopolysaccharide (LPS) component of the outer membrane of E. coli , is linked to the innermost lipid A through the core region of LPS of which five different structures, denominated K-12, R1, R2, R3 and R4, have been characterized so far. The aim of the present study was to analyse the prevalence of these LPS core types in the E. coli species and explore their distribution in the different E. coli phylogenetic groups and in relationship with the virulence gene repertoire. Results indicated an uneven distribution of core types between the different phylogroups, with phylogroup A strains being the most diverse in terms of LPS core types, while phylogroups B1, D and E strains were dominated by the R3 type, and phylogroups B2 and C strains were dominated by the R1 type. Strains carrying the LEE virulence operon were mostly of the R3 type whatever the phylogroup while, within phylogroup B2, strains carrying a K-12 core all belonged to the complex STc131, one of the major clones of extraintestinal pathogenic E. coli (ExPEC) strains. The origin of this uneven distribution is discussed but remains to be fully explained, as well as the consequences of carrying a specific core type on the wider aspects of bacterial phenotype.


2019 ◽  
Vol 5 (7) ◽  
Author(s):  
Charles H. D. Williamson ◽  
Nathan E. Stone ◽  
Amalee E. Nunnally ◽  
Heidie M. Hornstra ◽  
David M. Wagner ◽  
...  

Clostridioides difficile is a ubiquitous, diarrhoeagenic pathogen often associated with healthcare-acquired infections that can cause a range of symptoms from mild, self-limiting disease to toxic megacolon and death. Since the early 2000s, a large proportion of C. difficile cases have been attributed to the ribotype 027 (RT027) lineage, which is associated with sequence type 1 (ST1) in the C. difficile multilocus sequence typing scheme. The spread of ST1 has been attributed, in part, to resistance to fluoroquinolones used to treat unrelated infections, which creates conditions ideal for C. difficile colonization and proliferation. In this study, we analysed 27 isolates from a healthcare network in northern Arizona, USA, and 1352 publicly available ST1 genomes to place locally sampled isolates into a global context. Whole genome, single nucleotide polymorphism analysis demonstrated that at least six separate introductions of ST1 were observed in healthcare facilities in northern Arizona over an 18-month sampling period. A reconstruction of transmission networks identified potential nosocomial transmission of isolates, which were only identified via whole genome sequence analysis. Antibiotic resistance heterogeneity was observed among ST1 genomes, including variability in resistance profiles among locally sampled ST1 isolates. To investigate why ST1 genomes are so common globally and in northern Arizona, we compared all high-quality C. difficile genomes and identified that ST1 genomes have gained and lost a number of genomic regions compared to all other C. difficile genomes; analyses of other toxigenic C. difficile sequence types demonstrate that this loss may be anomalous and could be related to niche specialization. These results suggest that a combination of antimicrobial resistance and gain and loss of specific genes may explain the prominent association of this sequence type with C. difficile infection cases worldwide. The degree of genetic variability in ST1 suggests that classifying all ST1 genomes into a quinolone-resistant hypervirulent clone category may not be appropriate. Whole genome sequencing of clinical C. difficile isolates provides a high-resolution surveillance strategy for monitoring persistence and transmission of C. difficile and for assessing the performance of infection prevention and control strategies.


Sign in / Sign up

Export Citation Format

Share Document