Draft Genome Sequences of Two Clostridium botulinum Group II Strains Carrying Phage-Like Plasmids.

Clostridium botulinum is responsible for botulism, a potentially lethal foodborne intoxication. Here, we report the draft genome sequences of C. botulinum group II strains 202F (serotype F) and Hazen (serotype E). The genomes share many similarities, including multiple mobile genetic elements.

Salmonella enterica subsp. enterica virulence potential can be linked to higher survival within a dynamic in vitro human gastrointestinal model.

Salmonella enterica subsp. enterica is one of the leading causes of human foodborne infections and several outbreaks are now associated with the consumption of fresh fruit and vegetables. This study aims at evaluating whether Salmonella virulence can be linked to an enhanced ability to survive successive digestive environments. Thirteen S. enterica strains were selected according to high and low virulence phenotypes. Lettuce inoculated separately with each S. enterica strain was used as food matrix in the TNO gastrointestinal model (TIM-1) of the human upper gastrointestinal tract. During the passage in the stomach, counts determined using PMA-qPCR were 2-5 logs higher than the cultivable counts for all strains indicating the presence of viable but non-cultivable cells. Bacterial growth was observed in the duodenum compartment after 180 min for all but one strain and growth continued into the ileal compartment. After passage through the simulated gastrointestinal tract, both virulent and avirulent S. enterica strains survived but high virulence strains had a significantly (p = 0.004) better average survival rate (1003 %-3753 %) than low virulence strains (from 25 % to 3730%). The survival rates of S. enterica strains could be linked to the presence of genes associated with acid and bile resistance and their predicted products. The presence of single nucleotide polymorphisms may also impact the function of virulence associated genes and play a role in the resulting phenotype. These data provide an understanding of the relationship between measured virulence potential and survival of S. enterica during dynamic simulated gastrointestinal transit.

Genomic and phenotypic analysis of SspH1 identifies a new Salmonella effector, SspH3.

Salmonella is a major foodborne pathogen and is responsible for a range of diseases. Not all Salmonella contributes to severe health outcomes as there is a large degree of genetic heterogeneity among the 2,600 serovars within the genus. This variability across Salmonella serovars is linked to numerous genetic elements that dictate virulence. While several genetic elements encode virulence factors with well-documented contributions to pathogenesis, many genetic elements implicated in Salmonella virulence remain uncharacterized. Many pathogens encode a family of E3 ubiquitin ligases that are delivered into the cells that they infect using a Type 3 Secretion System (T3SS). These effectors, known as NEL-domain E3s, were first characterized in Salmonella. Most Salmonella encodes the NEL-effectors sspH2 and slrP, whereas only a subset of Salmonella encodes sspH1. SspH1 has been shown to ubiquitinate the mammalian protein kinase PKN1, which has been reported to negatively regulate the pro-survival program Akt. We discovered that SspH1 mediates the degradation of PKN1 during infection of a macrophage cell line but that this degradation does not impact Akt signaling. Genomic analysis of a large collection of Salmonella genomes identified a putative new gene, sspH3, with homology to sspH1. SspH3 is a novel NEL-domain effector.

Bacteriophage-Induced Lipopolysaccharide Mutations in Escherichia coli Lead to Hypersensitivity to Food Grade Surfactant Sodium Dodecyl Sulfate.

Bacteriophages (phages) are considered as one of the most promising antibiotic alternatives in combatting bacterial infectious diseases. However, one concern of employing phage application is the emergence of bacteriophage-insensitive mutants (BIMs). Here, we isolated six BIMs from E. coli B in the presence of phage T4 and characterized them using genomic and phenotypic methods. Of all six BIMs, a six-amino acid deletion in glucosyltransferase WaaG likely conferred phage resistance by deactivating the addition of T4 receptor glucose to the lipopolysaccharide (LPS). This finding was further supported by the impaired phage adsorption to BIMs and glycosyl composition analysis which quantitatively confirmed the absence of glucose in the LPS of BIMs. Since LPSs actively maintain outer membrane (OM) permeability, phage-induced truncations of LPSs destabilized the OM and sensitized BIMs to various substrates, especially to the food-grade surfactant sodium dodecyl sulfate (SDS). This hypersensitivity to SDS was exploited to design a T4-SDS combination which successfully prevented the generation of BIMs and eliminated the inoculated bacteria. Collectively, phage-driven modifications of LPSs immunized BIMs from T4 predation but increased their susceptibilities as a fitness cost. The findings of this study suggest a novel strategy to enhance the effectiveness of phage-based food safety interventions.

The Salmonella enterica Plasmidome as a Reservoir of Antibiotic Resistance.

The emergence of multidrug-resistant bacterial strains worldwide has become a serious problem for public health over recent decades. The increase in antimicrobial resistance has been expanding via plasmids as mobile genetic elements encoding antimicrobial resistance (AMR) genes that are transferred vertically and horizontally. This study focuses on Salmonella enterica, one of the leading foodborne pathogens in industrialized countries. S. enterica is known to carry several plasmids involved not only in virulence but also in AMR. In the current paper, we present an integrated strategy to detect plasmid scaffolds in whole genome sequencing (WGS) assemblies. We developed a two-step procedure to predict plasmids based on i) the presence of essential elements for plasmid replication and mobility, as well as ii) sequence similarity to a reference plasmid. Next, to confirm the accuracy of the prediction in 1750 S. enterica short-read sequencing data, we combined Oxford Nanopore MinION long-read sequencing with Illumina MiSeq short-read sequencing in hybrid assemblies for 84 isolates to evaluate the proportion of plasmid that has been detected. At least one scaffold with an origin of replication (ORI) was predicted in 61.3% of the Salmonella isolates tested. The results indicated that IncFII and IncI1 ORIs were distributed in many S. enterica serotypes and were the most prevalent AMR genes carrier, whereas IncHI2A/IncHI2 and IncA/C2 were more serotype restricted but bore several AMR genes. Comparison between hybrid and short-read assemblies revealed that 81.1% of plasmids were found in the short-read sequencing using our pipeline. Through this process, we established that plasmids are prevalent in S. enterica and we also substantially expand the AMR genes in the resistome of this species.

Combining Whole-Genome Sequencing and Multimodel Phenotyping To Identify Genetic Predictors of Salmonella Virulence.

Salmonella comprises more than 2,600 serovars. Very few environmental and uncommon serovars have been characterized for their potential role in virulence and human infections. A complementary in vitro and in vivo systematic high-throughput analysis of virulence was used to elucidate the association between genetic and phenotypic variations across Salmonella isolates. The goal was to develop a strategy for the classification of isolates as a benchmark and predict virulence levels of isolates. Thirty-five phylogenetically distant strains of unknown virulence were selected from the Salmonella Foodborne Syst-OMICS (SalFoS) collection, representing 34 different serovars isolated from various sources. Isolates were evaluated for virulence in 4 complementary models of infection to compare virulence traits with the genomics data, including interactions with human intestinal epithelial cells, human macrophages, and amoeba. In vivo testing was conducted using the mouse model of Salmonella systemic infection. Significant correlations were identified between the different models. We identified a collection of novel hypothetical and conserved proteins associated with isolates that generate a high burden. We also showed that blind prediction of virulence of 33 additional strains based on the pan-genome was high in the mouse model of systemic infection (82% agreement) and in the human epithelial cell model (74% agreement). These complementary approaches enabled us to define virulence potential in different isolates and present a novel strategy for risk assessment of specific strains and for better monitoring and source tracking during outbreaks.IMPORTANCESalmonella species are bacteria that are a major source of foodborne disease through contamination of a diversity of foods, including meat, eggs, fruits, nuts, and vegetables. More than 2,600 different Salmonella enterica serovars have been identified, and only a few of them are associated with illness in humans. Despite the fact that they are genetically closely related, there is enormous variation in the virulence of different isolates of Salmonella enterica Identification of foodborne pathogens is a lengthy process based on microbiological, biochemical, and immunological methods. Here, we worked toward new ways of integrating whole-genome sequencing (WGS) approaches into food safety practices. We used WGS to build associations between virulence and genetic diversity within 83 Salmonella isolates representing 77 different Salmonella serovars. Our work demonstrates the potential of combining a genomics approach and virulence tests to improve the diagnostics and assess risk of human illness associated with specific Salmonella isolates.

Phenomic and genomic approaches to studying the inhibition of multiresistant Salmonella enterica by microcin J25.

In livestock production, antibiotics are used to promote animal growth, control infections and thereby increase profitability. This practice has led to the emergence of multiresistant bacteria such as Salmonella, of which some serovars are disseminated in the environment. The objective of this study is to evaluate microcin J25 as an inhibitor of Salmonella enterica serovars of various origins including human, livestock and food. Among the 116 isolates tested, 37 (31.8%) were found resistant to at least one antibiotic, and 28 were multiresistant with 19 expressing the penta-resistant phenotype ACSSuT. Microcin J25 inhibited all isolates, with minimal inhibitory concentration values ranging from 0.06 µg/ml (28.4 nM) to 400 µg/ml (189 µM). Interestingly, no cross-resistance was found between microcin J25 and antibiotics. Multiple sequence alignments of genes encoding for the different proteins involved in the recognition and transport of microcin J25 showed that only ferric-hydroxamate uptake is an essential determinant for susceptibility of S. enterica to microcin J25. Examination of Salmonella strains exposed to microcin J25 by transmission electronic microscopy showed for the first-time involvement of a pore formation mechanism. Microcin J25 was a strong inhibitor of several multiresistant isolates of Salmonella and may have a great potential as an alternative to antibiotics.

A megaplasmid family driving dissemination of multidrug resistance in Pseudomonas.

Multidrug resistance (MDR) represents a global threat to health. Here, we used whole genome sequencing to characterise Pseudomonas aeruginosa MDR clinical isolates from a hospital in Thailand. Using long-read sequence data we obtained complete sequences of two closely related megaplasmids (>420 kb) carrying large arrays of antibiotic resistance genes located in discrete, complex and dynamic resistance regions, and revealing evidence of extensive duplication and recombination events. A comprehensive pangenomic and phylogenomic analysis indicates that: 1) these large plasmids comprise an emerging family present in different members of the Pseudomonas genus, and associated with multiple sources (geographical, clinical or environmental); 2) the megaplasmids encode diverse niche-adaptive accessory traits, including multidrug resistance; 3) the accessory genome of the megaplasmid family is highly flexible and diverse. The history of the megaplasmid family, inferred from our analysis of the available database, suggests that members carrying multiple resistance genes date back to at least the 1970s.

Bacteriophage-Insensitive Mutants of Antimicrobial-Resistant Salmonella Enterica are Altered in their Tetracycline Resistance and Virulence in Caco-2 Intestinal Cells.

Bacteriophages have shown promise as therapeutic alternatives to antibiotics for the control of infectious bacteria, including the human pathogen Salmonella. However, the development of effective phage-based applications requires the elucidation of key interactions between phages and target hosts, particularly since host resistance to phage is inevitable. Little is known about the alteration of host phenotypes following the development of resistance to phage. The aim of this study is to evaluate the antibiotic susceptibility and virulence of a Salmonella isolate following the development of resistance to bacteriophage SI1. We observed enhanced susceptibility to tetracycline and decreased invasion capacity in a differentiated Caco-2 intestinal cell line. Whole genome sequence analysis revealed an array of mutations, most notably, truncations in vgrG1_2, a core gene involved in Type VI secretion and mutations in the lipopolysaccharide, thereby indicating the plausible attachment site of phage SI1. These findings shed light on understanding the underlying mechanism for phage immunity within the host. Importantly, we reveal an associated genetic cost to the bacterial host with developing resistance to phages. Taken together, these results will aid in advancing strategies to delay or eliminate the development of host resistance when designing informed phage-based antimicrobials.

Major Release of 161 Whole-Genome Sequences from the International Pseudomonas Consortium Database.

Pseudomonas aeruginosa is an environmental bacterium and opportunistic pathogen. Here, we present draft genome sequences for 161 isolates from diverse clinical and environmental sources. This set of genome sequences complements other major public data releases from the International Pseudomonas Consortium Database.

Genomics of antibiotic-resistance prediction in Pseudomonas aeruginosa.

Antibiotic resistance is a worldwide health issue spreading quickly among human and animal pathogens, as well as environmental bacteria. Misuse of antibiotics has an impact on the selection of resistant bacteria, thus contributing to an increase in the occurrence of resistant genotypes that emerge via spontaneous mutation or are acquired by horizontal gene transfer. There is a specific and urgent need not only to detect antimicrobial resistance but also to predict antibiotic resistance in silico. We now have the capability to sequence hundreds of bacterial genomes per week, including assembly and annotation. Novel and forthcoming bioinformatics tools can predict the resistome and the mobilome with a level of sophistication not previously possible. Coupled with bacterial strain collections and databases containing strain metadata, prediction of antibiotic resistance and the potential for virulence are moving rapidly toward a novel approach in molecular epidemiology. Here, we present a model system in antibiotic-resistance prediction, along with its promises and limitations. As it is commonly multidrug resistant, Pseudomonas aeruginosa causes infections that are often difficult to eradicate. We review novel approaches for genotype prediction of antibiotic resistance. We discuss the generation of microbial sequence data for real-time patient management and the prediction of antimicrobial resistance.

The Pseudomonas aeruginosa Pan-Genome Provides New Insights on Its Population Structure, Horizontal Gene Transfer, and Pathogenicity.

The huge increase in the availability of bacterial genomes led us to a point in which we can investigate and query pan-genomes, for example, the full set of genes of a given bacterial species or clade. Here, we used a data set of 1,311 high-quality genomes from the human pathogen Pseudomonas aeruginosa, 619 of which were newly sequenced, to show that a pan-genomic approach can greatly refine the population structure of bacterial species, provide new insights to define species boundaries, and generate hypotheses on the evolution of pathogenicity. The 665-gene P. aeruginosa core genome presented here, which constitutes only 1% of the entire pan-genome, is the first to be in the same order of magnitude as the minimal bacterial genome and represents a conservative estimate of the actual core genome. Moreover, the phylogeny based on this core genome provides strong evidence for a five-group population structure that includes two previously undescribed groups of isolates. Comparative genomics focusing on antimicrobial resistance and virulence genes showed that variation among isolates was partly linked to this population structure. Finally, we hypothesized that horizontal gene transfer had an important role in this respect, and found a total of 3,010 putative complete and fragmented plasmids, 5% and 12% of which contained resistance or virulence genes, respectively. This work provides data and strategies to study the evolutionary trajectories of resistance and virulence in P. aeruginosa.

Genomic characterisation of an international Pseudomonas aeruginosa reference panel indicates that the two major groups draw upon distinct mobile gene pools.

Pseudomonas aeruginosa is an important opportunistic pathogen, especially in the context of infections of cystic fibrosis (CF). In order to facilitate coordinated study of this pathogen, an international reference panel of P. aeruginosa isolates was assembled. Here we report the genome sequencing and analysis of 33 of these isolates and 7 reference genomes to further characterise this panel. Core genome single nucleotide variant phylogeny demonstrated that the panel strains are widely distributed amongst the P. aeruginosa population. Common loss-of-function mutations reported as adaptive during CF (such as in mucA and mexA) were identified amongst isolates from chronic respiratory infections. From the 40 strains analysed, 37 unique resistomes were predicted, based on the Resistance Gene Identifier method using the Comprehensive Antibiotic Resistance Database. Notably, hierarchical clustering and phylogenetic reconstructions based on the presence/absence of genomic islands (GIs), prophages and other regions of genome plasticity (RGPs) supported the subdivision of P. aeruginosa into two main groups. This is the largest, most diverse analysis of GIs and associated RGPs to date, and the results suggest that, at least at the largest clade grouping level (group 1 vs group 2), each group may be drawing upon distinct mobile gene pools.

Salmonella enterica Prophage Sequence Profiles Reflect Genome Diversity and Can Be Used for High Discrimination Subtyping.

Non-typhoidal Salmonella is a leading cause of foodborne illness worldwide. Prompt and accurate identification of the sources of Salmonella responsible for disease outbreaks is crucial to minimize infections and eliminate ongoing sources of contamination. Current subtyping tools including single nucleotide polymorphism (SNP) typing may be inadequate, in some instances, to provide the required discrimination among epidemiologically unrelated Salmonella strains. Prophage genes represent the majority of the accessory genes in bacteria genomes and have potential to be used as high discrimination markers in Salmonella. In this study, the prophage sequence diversity in different Salmonella serovars and genetically related strains was investigated. Using whole genome sequences of 1,760 isolates of S. enterica representing 151 Salmonella serovars and 66 closely related bacteria, prophage sequences were identified from assembled contigs using PHASTER. We detected 154 different prophages in S. enterica genomes. Prophage sequences were highly variable among S. enterica serovars with a median ± interquartile range (IQR) of 5 ± 3 prophage regions per genome. While some prophage sequences were highly conserved among the strains of specific serovars, few regions were lineage specific. Therefore, strains belonging to each serovar could be clustered separately based on their prophage content. Analysis of S. Enteritidis isolates from seven outbreaks generated distinct prophage profiles for each outbreak. Taken altogether, the diversity of the prophage sequences correlates with genome diversity. Prophage repertoires provide an additional marker for differentiating S. enterica subtypes during foodborne outbreaks.

Whole-Genome Sequencing of Lactobacillus Species from Two Commercial Probiotic Products.

Eight Lactobacillus strains, each intrinsically resistant to an antibiotic, were isolated from two commercial probiotic products. Whole-genome sequencing identified two efflux transporters, a multidrug and extrusion protein (MATE) efflux transporter, and LmrCD, which may contribute to their intrinsic antibiotic resistance and may therefore facilitate their survival in the intestinal microbiota following antibiotic therapy.

Genomic characterisation of environmental Pseudomonas aeruginosa isolated from dental unit waterlines revealed the insertion sequence ISPa11 as a chaotropic element.

The bacterium Pseudomonas aeruginosa is well known to have a remarkable adaptive capacity allowing it to colonise many environments. A recent study on environmental isolates from dental unit waterlines (DUWLs) hinted at a genetic clustering into two groups. Isolates from one of these groups, named cluster III, were shown to have unusual phenotypes for environmental isolates, such as an increased biofilm production. To have a better ecological view, more specifically on isolates from cluster III, the complete genomes of 39 isolates including 16 from DUWLs were sequenced. In addition to an investigation of antibiotic resistance and secretion system gene content, a molecular phylogeny allowed confirmation of the split of the 16 environmental isolates in two groups and also sheds light on a correlation between the phylogenetic positions and the serotypes of the isolates. Isolates from cluster III harboured insertion sequences ISPa11 inserted into the O-specific antigen biosynthesis clusters and the gene lasR, encoding for a master regulator of the quorum sensing. Investigation of key regulators revealed another truncated gene, gacS. Alteration in lasR and gacS genes was consistent with phenotypic assays confirming their inactivation. These results bring new perspectives regarding the ecological adaptive potential of P. aeruginosa.

Prophage Integrase Typing Is a Useful Indicator of Genomic Diversity in Salmonella enterica.

Salmonella enterica is a bacterial species that is a major cause of illness in humans and food-producing animals. S. enterica exhibits considerable inter-serovar diversity, as evidenced by the large number of host adapted serovars that have been identified. The development of methods to assess genome diversity in S. enterica will help to further define the limits of diversity in this foodborne pathogen. Thus, we evaluated a PCR assay, which targets prophage integrase genes, as a rapid method to investigate S. enterica genome diversity. To evaluate the PCR prophage integrase assay, 49 isolates of S. enterica were selected, including 19 clinical isolates from clonal serovars (Enteritidis and Heidelberg) that commonly cause human illness, and 30 isolates from food-associated Salmonella serovars that rarely cause human illness. The number of integrase genes identified by the PCR assay was compared to the number of integrase genes within intact prophages identified by whole genome sequencing and phage finding program PHASTER. The PCR assay identified a total of 147 prophage integrase genes within the 49 S. enterica genomes (79 integrase genes in the food-associated Salmonella isolates, 50 integrase genes in S. Enteritidis, and 18 integrase genes in S. Heidelberg). In comparison, whole genome sequencing and PHASTER identified a total of 75 prophage integrase genes within 102 intact prophages in the 49 S. enterica genomes (44 integrase genes in the food-associated Salmonella isolates, 21 integrase genes in S. Enteritidis, and 9 integrase genes in S. Heidelberg). Collectively, both the PCR assay and PHASTER identified the presence of a large diversity of prophage integrase genes in the food-associated isolates compared to the clinical isolates, thus indicating a high degree of diversity in the food-associated isolates, and confirming the clonal nature of S. Enteritidis and S. Heidelberg. Moreover, PHASTER revealed a diversity of 29 different types of prophages and 23 different integrase genes within the food-associated isolates, but only identified four different phages and integrase genes within clonal isolates of S. Enteritidis and S. Heidelberg. These results demonstrate the potential usefulness of PCR based detection of prophage integrase genes as a rapid indicator of genome diversity in S. enterica.

Plasmid composition in Aeromonas salmonicida subsp. salmonicida 01-B526 unravels unsuspected type three secretion system loss patterns.

BACKGROUND: Aeromonas salmonicida subsp. salmonicida is a ubiquitous psychrophilic waterborne bacterium and a fish pathogen. The numerous mobile elements, especially insertion sequences (IS), in its genome promote rearrangements that impact its phenotype. One of the main virulence factors of this bacterium, its type three secretion system (TTSS), is affected by these rearrangements. In Aeromonas salmonicida subsp. salmonicida most of the TTSS genes are encoded in a single locus on a large plasmid called pAsa5, and may be lost when the bacterium is cultivated at a higher temperature (25 °C), producing non-virulent mutants. In a previous study, pAsa5-rearranged strains that lacked the TTSS locus on pAsa5 were produced using parental strains, including 01-B526. Some of the generated deletions were explained by homologous recombination between ISs found on pAsa5, whereas the others remained unresolved. To investigate those rearrangements, short- and long-read high-throughput sequencing technologies were used on the A. salmonicida subsp. salmonicida 01-B526 whole genome. RESULTS: Whole genome sequencing of the 01-B526 strain revealed that its pAsa5 has an additional IS copy, an ISAS5, compared to the reference strain (A449) sequence, which allowed for a previously unknown rearrangement to occur. It also appeared that 01-B526 bears a second large plasmid, named pAsa9, which shares 40 kbp of highly similar sequences with pAsa5. Following these discoveries, previously unexplained deletions were elucidated by genotyping. Furthermore, in one of the derived strains a fusion of pAsa5 and pAsa9, involving the newly discovered ISAS5 copy, was observed. CONCLUSION: The loss of TTSS and hence virulence is explained by one consistent mechanism: IS-driven homologous recombination. The similarities between pAsa9 and pAsa5 also provide another example of genetic diversity driven by ISs.

A Syst-OMICS Approach to Ensuring Food Safety and Reducing the Economic Burden of Salmonellosis.

The Salmonella Syst-OMICS consortium is sequencing 4,500 Salmonella genomes and building an analysis pipeline for the study of Salmonella genome evolution, antibiotic resistance and virulence genes. Metadata, including phenotypic as well as genomic data, for isolates of the collection are provided through the Salmonella Foodborne Syst-OMICS database (SalFoS), at https://salfos.ibis.ulaval.ca/. Here, we present our strategy and the analysis of the first 3,377 genomes. Our data will be used to draw potential links between strains found in fresh produce, humans, animals and the environment. The ultimate goals are to understand how Salmonella evolves over time, improve the accuracy of diagnostic methods, develop control methods in the field, and identify prognostic markers for evidence-based decisions in epidemiology and surveillance.