Unveiling the genomic landscape of Salmonella enterica serotypes Typhimurium, Newport, and Infantis in Latin American surface waters: a comparative analysis

Surface waters are considered ecological habitats where Salmonella enterica can persist and disseminate to fresh produce production systems. This study aimed to explore the genomic profiles of S. enterica serotypes Typhimurium, Newport, and Infantis from surface waters in Chile, Mexico, and Brazil collected between 2019 and 2022. We analyzed the whole genomes of 106 S. Typhimurium, 161 S. Newport, and 113 S. Infantis isolates. Our phylogenetic analysis exhibited distinct groupings of isolates by their respective countries except for a notable case involving a Chilean S. Newport isolate closely related to two Mexican isolates, showing 4 and 13 single nucleotide polymorphisms of difference, respectively. The patterns of the most frequently detected antimicrobial resistance genes varied across countries and serotypes. A strong correlation existed between integron carriage and genotypic multidrug resistance (MDR) across serotypes in Chile and Mexico (R > 0.90, P < 0.01), while integron(s) were not detected in any of the Brazilian isolates. By contrast, we did not identify any strong correlation between plasmid carriage and genotypic MDR across diverse countries and serotypes.IMPORTANCEUnveiling the genomic landscape of S. enterica in Latin American surface waters is pivotal for ensuring public health. This investigation sheds light on the intricate genomic diversity of S. enterica in surface waters across Chile, Mexico, and Brazil. Our research also addresses critical knowledge gaps, pioneering a comprehensive understanding of surface waters as a reservoir for multidrug-resistant S. enterica. By integrating our understanding of integron carriage as biomarkers into broader MDR control strategies, we can also work toward targeted interventions that mitigate the emergence and dissemination of MDR in S. enterica in surface waters. Given its potential implications for food safety, this study emphasizes the critical need for informed policies and collaborative initiatives to address the risks associated with S. enterica in surface waters.

Advancing metagenome-assembled genome-based pathogen identification: unraveling the power of long-read assembly algorithms in Oxford Nanopore sequencing

Oxford Nanopore sequencing is one of the high-throughput sequencing technologies that facilitates the reconstruction of metagenome-assembled genomes (MAGs). This study aimed to assess the potential of long-read assembly algorithms in Oxford Nanopore sequencing to enhance the MAG-based identification of bacterial pathogens using both simulated and mock communities. Simulated communities were generated to mimic those on fresh spinach and in surface water. Long reads were produced using R9.4.1+SQK-LSK109 and R10.4 + SQK-LSK112, with 0.5, 1, and 2 million reads. The simulated bacterial communities included multidrug-resistant Salmonella enterica serotypes Heidelberg, Montevideo, and Typhimurium in the fresh spinach community individually or in combination, as well as multidrug-resistant Pseudomonas aeruginosa in the surface water community. Real data sets of the ZymoBIOMICS HMW DNA Standard were also studied. A bioinformatic pipeline (MAGenie, freely available at https://github.com/jackchen129/MAGenie) that combines metagenome assembly, taxonomic classification, and sequence extraction was developed to reconstruct draft MAGs from metagenome assemblies. Five assemblers were evaluated based on a series of genomic analyses. Overall, Flye outperformed the other assemblers, followed by Shasta, Raven, and Unicycler, while Canu performed least effectively. In some instances, the extracted sequences resulted in draft MAGs and provided the locations and structures of antimicrobial resistance genes and mobile genetic elements. Our study showcases the viability of utilizing the extracted sequences for precise phylogenetic inference, as demonstrated by the consistent alignment of phylogenetic topology between the reference genome and the extracted sequences. R9.4.1+SQK-LSK109 was more effective in most cases than R10.4+SQK-LSK112, and greater sequencing depths generally led to more accurate results.IMPORTANCEBy examining diverse bacterial communities, particularly those housing multiple Salmonella enterica serotypes, this study holds significance in uncovering the potential of long-read assembly algorithms to improve metagenome-assembled genome (MAG)-based pathogen identification through Oxford Nanopore sequencing. Our research demonstrates that long-read assembly stands out as a promising avenue for boosting precision in MAG-based pathogen identification, thus advancing the development of more robust surveillance measures. The findings also support ongoing endeavors to fine-tune a bioinformatic pipeline for accurate pathogen identification within complex metagenomic samples.

Complete genome sequences of 17 Salmonella enterica serovar Schwarzengrund isolates carrying an IncFIB-IncFIC (FII) fusion plasmid

Seventeen Salmonella enterica serovar Schwarzengrund isolates from chicken (n = 9) and clinical samples including stool (n = 6), urine (n = 1), and gallbladder (n = 1) were sequenced and found to carry an IncFIB-IncFIC (FII) fusion plasmid of approximately 145 Kb. This information provides reference genomic data for comparative studies of S. Schwarzengrund pathogenicity and plasmid genetics.

Application of quasimetagenomics methods to define microbial diversity and subtype Listeria monocytogenes in dairy and seafood production facilities

IMPORTANCE

In developed countries, the human diet is predominated by food commodities, which have been manufactured, processed, and stored in a food production facility. Little is known about the application of metagenomic sequencing approaches for detecting foodborne pathogens, such as L. monocytogenes, and characterizing microbial diversity in food production ecosystems. In this work, we investigated the utility of 16S rRNA amplicon and quasimetagenomic sequencing for the taxonomic and phylogenetic classification of Listeria culture enrichments of environmental swabs collected from dairy and seafood production facilities. We demonstrated that single-nucleotide polymorphism (SNP) analyses of L. monocytogenes metagenome-assembled genomes (MAGs) from quasimetagenomic data sets can achieve similar resolution as culture isolate whole-genome sequencing. To further understand the impact of genome coverage on MAG SNP cluster resolution, an in silico downsampling approach was employed to reduce the percentage of target pathogen sequence reads, providing an initial estimate of required MAG coverage for subtyping resolution of L. monocytogenes.

bettercallsal: better calling of Salmonella serotypes from enrichment cultures using shotgun metagenomic profiling and its application in an outbreak setting

Most current Salmonella subtyping analyses rely on whole genome sequencing (WGS), which focuses on the high-resolution analysis of single genomes or multiple single genomes from the isolated colonies on microbiological agar plates. In this study, we introduce bioinformatics innovations for a metagenomic outbreak response workflow that accurately identifies multiple Salmonella serovars at the same time. bettercallsal is one of the first analysis tools to identify multiple Salmonella enterica serotypes from metagenomic or quasi-metagenomic datasets with high accuracy, allowing these isolate-independent methods to be incorporated into surveillance and root cause investigations. It was tested on an in silico benchmark dataset comprising 29 unique Salmonella serovars, 46 non-Salmonella bacterial genomes, and 10 viral genomes at varying read depths and on previously well-characterized and sequenced non-selective primary and selective enrichments of papaya and peach samples from separate outbreak investigations that resulted in the identification of multiple Salmonella serovars using traditional isolate culturing and WGS as well as nucleic acid assays. Analyses were also conducted on these datasets using a custom-built k-mer tool, SeqSero2, and Kallisto to compare serotype calling to bettercallsal. The in silico dataset analyzed with bettercallsal achieved the maximum precision, recall, and accuracy of 100, 83, and 94%, respectively. In the papaya outbreak samples, bettercallsal identified the presence of multiple serovars in agreement with the Luminex® xMAP assay results and also identified more serovars per sample, as evidenced by NCBI SNP clustering. In peach outbreak samples, bettercallsal identified two serovars in concordance with k-mer analysis and the Luminex xMAP assay. The genome hit reported by bettercallsal clustered with the chicken isolate genome, as reported by the FDA peach outbreak investigation from sequenced isolates (WGS). Overall, bettercallsal outperformed k-mer, Seqsero2, and Kallisto in identifying multiple serovars from enrichment cultures using shotgun metagenomic sequencing.

Epidemiological Cutoff Values for Standard Broth Microdilution and Disk Diffusion Susceptibility Testing of Aeromonas hydrophila Isolated from Fish

Aeromonas hydrophila and other closely related Aeromonas species cause motile aeromonad septicemia, a common fish disease. The disease affects many aquaculture sectors potentially requiring antimicrobial treatments. Therefore, researchers and laboratory diagnosticians need criteria called epidemiological cutoff values (ECVs) to determine whether a bacterial isolate has developed decreased susceptibility to an antimicrobial. To generate ECVs for this bacterium, we assembled a diverse collection of 245 isolates previously identified as A. hydrophila from fish. Using rpoD sequencing, we confirmed that 97 of the 245 isolates were A. hydrophila. We allocated the isolates among three laboratories and tested their susceptibility against eight antimicrobials using standard Clinical and Laboratory Standards Institute (CLSI) disk diffusion and broth microdilution methods. The resulting frequency distributions were statistically analyzed to determine wild-type cutoff estimates, which, along with scatterplots, were used to estimate potential ECVs. In collaboration with the CLSI, aquaculture working group, we proposed ECVs for six of the eight antimicrobials tested. Subsequently, the CLSI Subcommittee on Veterinary Antimicrobial Susceptibility Testing reviewed our data and approved the ECVs to be added to the 2020 edition of the VET04 performance standards for antimicrobial susceptibility testing of aquatic bacteria.

Microbiome Population Dynamics of Cold-Smoked Sockeye Salmon during Refrigerated Storage and after Culture Enrichment

ABSTRACT

Cold-smoked salmon is a ready-to-eat seafood product of high commercial importance. The processing and storage steps facilitate the introduction, growth, and persistence of foodborne pathogens and spoilage bacteria. The growth of commensal bacteria during storage and once the product is opened also influence the quality and safety of cold-smoked salmon. Here we investigated the microbial community through targeted 16S rRNA gene and shotgun metagenomic sequencing as means to better understand the interactions among bacteria in cold-smoked salmon. Cold-smoked salmon samples were tested over 30 days of aerobic storage at 4°C and cultured at each time point in a buffered Listeria enrichment broth (BLEB) commonly used to detect Listeria in foods. The microbiomes were composed of Firmicutes and Proteobacteria, namely, Carnobacterium, Brochothrix, Pseudomonas, Serratia, and Psychrobacter. Pseudomonas species were the most diverse species, with 181 taxa identified. In addition, we identified potential homologs to 10 classes of bacteriocins in microbiomes of cold-smoked salmon stored at 4°C and corresponding BLEB culture enrichments. The findings presented here contribute to our understanding of microbiome population dynamics in cold-smoked salmon, including changes in bacterial taxa during aerobic cold storage and after culture enrichment. This may facilitate improvements to pathogen detection and quality preservation of this food.

HIGHLIGHTS

Development and Validation of a Quantitative PCR Method for Species Verification and Serogroup Determination of Listeria monocytogenes Isolates

ABSTRACT

Listeria monocytogenes (Lm) is one of the leading causes of death because of foodborne illness, affecting the elderly, pregnant women, neonates, and people who are immunocompromised. Serologically, Lm can be classified into 13 serotypes, although only 4 are typically linked with food contamination and illness. Since 2000, a shift in serotypes involved in listeriosis outbreaks has been observed, suggesting that tracking of serotypes could help identify emerging trends. A PCR method developed in 2004 allowed detection of the four major serotypes as molecular serogroups, corresponding to broad phylogenetic groups. In this study, a novel quantitative PCR (qPCR) method was developed that uses two multiplex qPCRs, one to confirm the Listeria genus and Lm species and the second for Lm molecular serogrouping. This method was compared with the U.S. Food and Drug Administration Bacteriological Analytical Manual (BAM) method for Lm and the seroagglutination method, using a 208-strain panel. Comparison of the genus and species qPCR assay with the BAM methods found an equal or slightly higher accuracy for the qPCR method (>98%), compared with the BAM protocol (>96%), when evaluated against independent characterization data. Molecular serogrouping using the qPCR method (96.6%) was more accurate than the seroagglutination assay (75.6%). The qPCR method identified Lm 4bV strains, which could not be resolved using seroagglutination. The qPCR could not identify lineage III and IV serotype 4b strains but did correctly identify 16 of 18 lineage III and IV strains. The qPCR method performed genus identification for the Listeria species Lm, L. innocua, L. welshimeri, L. ivanovii, and L. seeligeri. In addition, the method performed species identification for Lm and classified Lm into six molecular serogroups: 2A, 2B, 2C, 4B, NT, and 4bV. This method provided a rapid and accurate confirmation of Lm and serogroup determinations; furthermore, it could help identify otherwise unlinked strains by enabling whole genome sequencing analysis based on broad phylogeny, independent of other information.

HIGHLIGHTS

Analysis of the Molecular Diversity Among Cronobacter Species Isolated From Filth Flies Using Targeted PCR, Pan Genomic DNA Microarray, and Whole Genome Sequencing Analyses

Cronobacter species are opportunistic pathogens capable of causing life-threatening infections in humans, with serious complications arising in neonates, infants, immuno-compromised individuals, and elderly adults. The genus is comprised of seven species: Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite a multiplicity of genomic data for the genus, little is known about likely transmission vectors. Using DNA microarray analysis, in parallel with whole genome sequencing, and targeted PCR analyses, the total gene content of two C. malonaticus, three C. turicensis, and 14 C. sakazaki isolated from various filth flies was assessed. Phylogenetic relatedness among these and other strains obtained during surveillance and outbreak investigations were comparatively assessed. Specifically, microarray analysis (MA) demonstrated its utility to cluster strains according to species-specific and sequence type (ST) phylogenetic relatedness, and that the fly strains clustered among strains obtained from clinical, food and environmental sources from United States, Europe, and Southeast Asia. This combinatorial approach was useful in data mining for virulence factor genes, and phage genes and gene clusters. In addition, results of plasmidotyping were in agreement with the species identity for each strain as determined by species-specific PCR assays, MA, and whole genome sequencing. Microarray and BLAST analyses of Cronobacter fly sequence datasets were corroborative and showed that the presence and absence of virulence factors followed species and ST evolutionary lines even though such genes were orthologous. Additionally, zebrafish infectivity studies showed that these pathotypes were as virulent to zebrafish embryos as other clinical strains. In summary, these findings support a striking phylogeny amongst fly, clinical, and surveillance strains isolated during 2010–2015, suggesting that flies are capable vectors for transmission of virulent Cronobacter spp.; they continue to circulate among United States and European populations, environments, and that this “pattern of circulation” has continued over decades.

Research Note: Effect of butyric acid glycerol esters on ileal and cecal mucosal and luminal microbiota in chickens challenged with Eimeria maxima

Coccidiosis is one of the most prevalent diseases seen in the poultry industry leading to excessive economic losses. The aim of this study was to investigate the effect of butyric acid glycerol esters (BE) on the ileal and cecal microbiota in birds challenged with Eimeria maxima (EM). Ross 708 male broilers were fed a diet supplemented with 0 (control) or 0.25% BE from day 1. On day 21, half of the birds were infected with 10³ EM oocysts. For determing microbiota, ileal and cecal contents and epithelial scrapings were collected at 7 and 10 D postinfection (PI). Alpha diversity of bacterial communities was mostly affected (P < 0.05) by time PI and EM infection. The richness of luminal bacterial populations in the ileum and ceca was affected (P < 0.05) by addition of BE and by time PI × EM × BE interaction, respectively. In the ileal and cecal luminal and mucosal bacterial communities, permutational multivariate analysis of variance (PERMANOVA, unweighted UniFrac) showed significant (P < 0.05) differences because of time PI and interaction between time PI, EM, and BE. Significant (P < 0.05) differences in taxonomic composition at the family level were observed in microbiota of luminal and mucosal populations of the ileum and ceca owing to time PI, EM, BE, and their interactions. The bacterial community present in the cecal lumen was characterized by the lowest number of differential bacteria, whereas the cecal mucosal community was characterized by the highest number of differentially abundant bacteria. In conclusion, our results show that EM infection and time PI has the biggest impact on microbial diversity in the chicken gut. The presence of BE in the diet had a limited effect on gut microbiota.

The Role of the Salmonella spvB IncF Plasmid and Its Resident Entry Exclusion Gene traS on Plasmid Exclusion

Salmonella enterica cause significant illnesses worldwide. There has been a marked increase in resistance to fluoroquinolones and β-lactams/cephalosporins, antibiotics commonly used to treat salmonellosis. However, S. enterica serovars vary in their resistance to these and other antibiotics. The systemic virulence of some Salmonella serovars is due to a low copy number, IncF plasmid (65–100 kb) that contains the ADP-ribosylating toxin, SpvB. This virulence plasmid is present in only nine Salmonella serovars. It is possible that the spvB-virulence plasmid excludes other plasmids and may explain why antibiotic resistance is slow to develop in certain Salmonella serovars such as S. Enteritidis. The distribution of plasmid entry exclusion genes traS/traT and traY/excA are variable in Salmonella IncF and IncI plasmids, respectively and may account for differences in emergent antimicrobial resistance for some Salmonella serovars. The goal of this study is to determine the contribution of the Salmonella spvB-virulence plasmid in F-plasmid exclusion. From conjugation experiments, S. Typhimurium exhibited lower conjugation frequency with incFI and incFII plasmids when the spvB-virulence plasmid is present. Furthermore, introduction of cloned incFI traS into a “plasmidless” S. Typhimurium LT2 strain and Escherichia coli DH5α excluded incFI plasmid. However, deletion of the virulence plasmid traS did not affect plasmid exclusion significantly compared to a spvB control deletion. In addition, differences in F plasmid conjugation in natural Salmonella isolates did not correlate with IncF or SpvB-virulence plasmid genotype. There appear to be other plasmid or chromosomal genes at play in plasmid exclusion that may be responsible for the slow development of antibiotic resistance in certain serovars.

Global transcriptomic analyses of Salmonella enterica in Iron-depleted and Iron-rich growth conditions

Background

Salmonella enterica possess several iron acquisition systems, encoded on the chromosome and plasmids. Recently, we demonstrated that incompatibility group (Inc) FIB plasmid-encoded iron acquisition systems (Sit and aerobactin) likely play an important role in persistence of Salmonella in human intestinal epithelial cells (Caco-2). In this study, we sought to determine global transcriptome analyses of S. enterica in iron-rich (IR) and iron-depleted (ID) growth conditions.

Results

The number of differentially-expressed genes were substantially higher for recipient (SE819) (n = 966) and transconjugant (TC) (n = 945) compared to the wild type (WT) (SE163A) (n = 110) strain in ID as compared to IR growth conditions. Several virulence-associated factors including T3SS, flagellin, cold-shock protein (cspE), and regulatory genes were upregulated in TC in ID compared to IR conditions. Whereas, IS1 and acrR/tetR transposases located on the IncFIB plasmid, ferritin and several regulatory genes were downregulated in TC in ID conditions. Enterobactin transporter (entS), iron ABC transporter (fepCD), colicin transporter, IncFIB-encoded enolase, cyclic di-GMP regulator (cdgR) and other regulatory genes of the WT strain were upregulated in ID compared to IR conditions. Conversely, ferritin, ferrous iron transport protein A (feoA), IncFIB-encoded IS1 and acrR/tetR transposases and ArtA toxin of WT were downregulated in ID conditions. SDS-PAGE coupled with LC-MS/MS analyses revealed that siderophore receptor proteins such as chromosomally-encoded IroN and, IncFIB-encoded IutA were upregulated in WT and TC in ID growth conditions. Both chromosome and IncFIB plasmid-encoded SitA was overexpressed in WT, but not in TC or recipient in ID conditions. Increased expression of flagellin was detected in recipient and TC, but not in WT in ID conditions.

Conclusion

Iron concentrations in growth media influenced differential gene expressions both at transcriptional and translational levels, including genes encoded on the IncFIB plasmid. Limited iron availability within the host may promote pathogenic Salmonella to differentially express subsets of genes encoded by chromosome and/or plasmids, facilitating establishment of successful infection.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5768-0) contains supplementary material, which is available to authorized users.

Microbiomes Associated With Foods From Plant and Animal Sources

Food microbiome composition impacts food safety and quality. The resident microbiota of many food products is influenced throughout the farm to fork continuum by farming practices, environmental factors, and food manufacturing and processing procedures. Currently, most food microbiology studies rely on culture-dependent methods to identify bacteria. However, advances in high-throughput DNA sequencing technologies have enabled the use of targeted 16S rRNA gene sequencing to profile complex microbial communities including non-culturable members. In this study we used 16S rRNA gene sequencing to assess the microbiome profiles of plant and animal derived foods collected at two points in the manufacturing process; post-harvest/pre-retail (cilantro) and retail (cilantro, masala spice mixes, cucumbers, mung bean sprouts, and smoked salmon). Our findings revealed microbiome profiles, unique to each food, that were influenced by the moisture content (dry spices, fresh produce), packaging methods, such as modified atmospheric packaging (mung bean sprouts and smoked salmon), and manufacturing stage (cilantro prior to retail and at retail). The masala spice mixes and cucumbers were comprised mainly of Proteobacteria, Firmicutes, and Actinobacteria. Cilantro microbiome profiles consisted mainly of Proteobacteria, followed by Bacteroidetes, and low levels of Firmicutes and Actinobacteria. The two brands of mung bean sprouts and the three smoked salmon samples differed from one another in their microbiome composition, each predominated by either by Firmicutes or Proteobacteria. These data demonstrate diverse and highly variable resident microbial communities across food products, which is informative in the context of food safety, and spoilage where indigenous bacteria could hamper pathogen detection, and limit shelf life.

A Comprehensive Evaluation of the Genetic Relatedness of Listeria monocytogenes Serotype 4b Variant Strains

Recently, we have identified a link between four listeriosis incidents/outbreaks to a variant of Listeria monocytogenes (Lm) serotype 4b strains, 4bV. Although 4bV strains have been reported from clinical specimens as well as from foods, listeriosis outbreaks occurring in 2014–2016 were the first reported outbreaks involving 4bV in the USA. Since traditional typing methods do not detect members of this group, we undertook a systematic and retrospective analysis of all Lm in the NCBI WGS Sequence Read Archive database to investigate the burden of 4bV strains among all listeriosis cases. This analysis identified the presence of isolates causing sporadic cases as well as those associated with the aforementioned outbreaks, as determined by WGS and traditional epidemiology. In total, approximately 350 Lm 4bV strains were identified from multiple parts of the USA as well as from Australia and Chile, dating back to 2001. The genomic relatedness of these strains was compared using the CFSAN SNP Pipeline and multi-virulence-locus sequence typing (MVLST). Using the CFSAN Pipeline tool, the 4bV strains were found to group into seven clusters that were separate from 4b strains. All seven clades appeared to contain isolates from both clinical and non-clinical sources. Conversely, the MVLST analysis revealed that practically all of the strains belonged to a single clade, suggesting that 4bV strains from disparate geographic regions and sources are under varied selective pressure, restricting diversity across these six virulence loci while allowing more variability across the genome as a whole. Further evaluation of these 4bV strains identified genes potentially acquired from a lineage II source external to the lmo0733–lmo0739 region, as well as highly conserved SNPs unique to the 4bV strains when compared to those from other lineages. Taken together, these data suggest that 4bV strains have undergone adaptive responses to selective pressures that may enhance survival in the environment while maintaining the pathogenic potential of serotype 4b strains.

High-Resolution Microbiome Profiling for Detection and Tracking of Salmonella enterica

16S rRNA community profiling continues to be a useful tool to study microbiome composition and dynamics, in part due to advances in next generation sequencing technology that translate into reductions in cost. Reliable taxonomic identification to the species-level, however, remains difficult, especially for short-read sequencing platforms, due to incomplete coverage of the 16S rRNA gene. This is especially true for Salmonella enterica, which is often found as a low abundant member of the microbial community, and is often found in combination with several other closely related enteric species. Here, we report on the evaluation and application of Resphera Insight, an ultra-high resolution taxonomic assignment algorithm for 16S rRNA sequences to the species level. The analytical pipeline achieved 99.7% sensitivity to correctly identify S. enterica from WGS datasets extracted from the FDA GenomeTrakr Bioproject, while demonstrating 99.9% specificity over other Enterobacteriaceae members. From low-diversity and low-complexity samples, namely ice cream, the algorithm achieved 100% specificity and sensitivity for Salmonella detection. As demonstrated using cilantro and chili powder, for highly complex and diverse samples, especially those that contain closely related species, the detection threshold will likely have to be adjusted higher to account for misidentifications. We also demonstrate the utility of this approach to detect Salmonella in the clinical setting, in this case, bloodborne infections.

Genotypic and phenotypic characterization of multidrug resistant Salmonella Typhimurium and Salmonella Kentucky strains recovered from chicken carcasses

Salmonella Typhimurium is the leading cause of human non-typhoidal gastroenteritis in the US. S. Kentucky is one the most commonly recovered serovars from commercially processed poultry carcasses. This study compared the genotypic and phenotypic properties of two Salmonella enterica strains Typhimurium (ST221_31B) and Kentucky (SK222_32B) recovered from commercially processed chicken carcasses using whole genome sequencing, phenotype characterizations and an intracellular killing assay. Illumina MiSeq platform was used for sequencing of two Salmonella genomes. Phylogenetic analysis employing homologous alignment of a 1,185 non-duplicated protein-coding gene in the Salmonella core genome demonstrated fully resolved bifurcating patterns with varying levels of diversity that separated ST221_31B and SK222_32B genomes into distinct monophyletic serovar clades. Single nucleotide polymorphism (SNP) analysis identified 2,432 (ST19) SNPs within 13 Typhimurium genomes including ST221_31B representing Sequence Type ST19 and 650 (ST152) SNPs were detected within 13 Kentucky genomes including SK222_32B representing Sequence Type ST152. In addition to serovar-specific conserved coding sequences, the genomes of ST221_31B and SK222_32B harbor several genomic regions with significant genetic differences. These included phage and phage-like elements, carbon utilization or transport operons, fimbriae operons, putative membrane associated protein-encoding genes, antibiotic resistance genes, siderophore operons, and numerous hypothetical protein-encoding genes. Phenotype microarray results demonstrated that ST221_31B is capable of utilizing certain carbon compounds more efficiently as compared to SK222_3B; namely, 1,2-propanediol, M-inositol, L-threonine, α-D-lactose, D-tagatose, adonitol, formic acid, acetoacetic acid, and L-tartaric acid. ST221_31B survived for 48 h in macrophages, while SK222_32B was mostly eliminated. Further, a 3-fold growth of ST221_31B was observed at 24 hours post-infection in chicken granulosa cells while SK222_32B was unable to replicate in these cells. These results suggest that Salmonella Typhimurium can survive host defenses better and could be more invasive than Salmonella Kentucky and provide some insights into the genomic determinants responsible for these differences.

A clade of Listeria monocytogenes serotype 4b variant strains linked to recent listeriosis outbreaks associated with produce from a defined geographic region in the US

Four listeriosis incidences/outbreaks, spanning 19 months, have been linked to Listeria monocytogenes serotype 4b variant (4bV) strains. Three of these incidents can be linked to a defined geographical region, while the fourth is likely to be linked. In this study, whole genome sequencing (WGS) of strains from these incidents was used for genomic comparisons using two approached. The first was JSpecies tetramer, which analyzed tetranucleotide frequency to assess relatedness. The second, the CFSAN SNP Pipeline, was used to perform WGS SNP analyses against three different reference genomes to evaluate relatedness by SNP distances. In each case, unrelated strains were included as controls. The analyses showed that strains from these incidents form a highly related clade with SNP differences of ≤101 within the clade and >9000 against other strains. Multi-Virulence-Locus Sequence Typing, a third standardized approach for evaluation relatedness, was used to assess the genetic drift in six conserved, known virulence loci and showed a different clustering pattern indicating possible differences in selection pressure experienced by these genes. These data suggest a high degree of relatedness among these 4bV strains linked to a defined geographic region and also highlight the possibility of alterations related to adaptation and virulence.

Use of a Pan-Genomic DNA Microarray in Determination of the Phylogenetic Relatedness among Cronobacter spp. and Its Use as a Data Mining Tool to Understand Cronobacter Biology

Cronobacter (previously known as Enterobacter sakazakii) is a genus of Gram-negative, facultatively anaerobic, oxidase-negative, catalase-positive, rod-shaped bacteria of the family Enterobacteriaceae. These organisms cause a variety of illnesses such as meningitis, necrotizing enterocolitis, and septicemia in neonates and infants, and urinary tract, wound, abscesses or surgical site infections, septicemia, and pneumonia in adults. The total gene content of 379 strains of Cronobacter spp. and taxonomically-related isolates was determined using a recently reported DNA microarray. The Cronobacter microarray as a genotyping tool gives the global food safety community a rapid method to identify and capture the total genomic content of outbreak isolates for food safety, environmental, and clinical surveillance purposes. It was able to differentiate the seven Cronobacter species from one another and from non-Cronobacter species. The microarray was also able to cluster strains within each species into well-defined subgroups. These results also support previous studies on the phylogenic separation of species members of the genus and clearly highlight the evolutionary sequence divergence among each species of the genus compared to phylogenetically-related species. This review extends these studies and illustrates how the microarray can also be used as an investigational tool to mine genomic data sets from strains. Three case studies describing the use of the microarray are shown and include: (1) the determination of allelic differences among Cronobacter sakazakii strains possessing the virulence plasmid pESA3; (2) mining of malonate and myo-inositol alleles among subspecies of Cronobacter dublinensis strains to determine subspecies identity; and (3) lastly using the microarray to demonstrate sequence divergence and phylogenetic relatedness trends for 13 outer-membrane protein alleles among 240 Cronobacter and phylogenetically-related strains. The goal of this review is to describe microarrays as a robust tool for genomics research of this assorted and important genus, a criterion toward the development of future preventative measures to eliminate this foodborne pathogen from the global food supply.

Early Recovery of Salmonella from Food Using a 6-Hour Non-selective Pre-enrichment and Reformulation of Tetrathionate Broth

Culture based methods are commonly employed to detect pathogens in food and environmental samples. These methods are time consuming and complex, requiring multiple non-selective and selective enrichment broths, and usually take at least 1 week to recover and identify pathogens. Improving pathogen detection in foods is a primary goal for regulatory agencies and industry. Salmonella detection in food relies on a series of culture steps in broth formulations optimized to resuscitate Salmonella and reduce the abundance of competitive bacteria. Examples of non-selective pre-enrichment broths used to isolate Salmonella from food include Lactose, Universal Pre-enrichment, BPW, and Trypticase Soy broths. Tetrathionate (TT) and Rappaport-Vassiliadis (RV) broths are employed after a 24-h non-selective enrichment to select for Salmonella and hamper the growth of competitive bacteria. In this study, we tested a new formulation of TT broth that lacks brilliant green dye and has lower levels of TT . We employed this TT broth formulation in conjunction with a 6-h non-selective pre-enrichment period and determined that Salmonella recovery was possible one day earlier than standard food culture methods. We tested the shortened culture method in different non-selective enrichment broths, enumerated Salmonella in the non-selective enrichments, and used 16S rRNA gene sequencing to determine the proportional abundances of Salmonella in the TT and RV selective enrichments. Together these data revealed that a 6-h non-selective pre-enrichment reduces the levels of competitive bacteria inoculated into the selective TT and RV broths, enabling the recovery of Salmonella 1 day earlier than standard culture enrichment methods.

Whole Genome Sequence Analysis Using JSpecies Tool Establishes Clonal Relationships between Listeria monocytogenes Strains from Epidemiologically Unrelated Listeriosis Outbreaks

In an effort to build a comprehensive genomic approach to food safety challenges, the FDA has implemented a whole genome sequencing effort, GenomeTrakr, which involves the sequencing and analysis of genomes of foodborne pathogens. As a part of this effort, we routinely sequence whole genomes of Listeria monocytogenes (Lm) isolates associated with human listeriosis outbreaks, as well as those isolated through other sources. To rapidly establish genetic relatedness of these genomes, we evaluated tetranucleotide frequency analysis via the JSpecies program to provide a cursory analysis of strain relatedness. The JSpecies tetranucleotide (tetra) analysis plots standardized (z-score) tetramer word frequencies of two strains against each other and uses linear regression analysis to determine similarity (r²). This tool was able to validate the close relationships between outbreak related strains from four different outbreaks. Included in this study was the analysis of Lm strains isolated during the recent caramel apple outbreak and stone fruit incident in 2014. We identified that many of the isolates from these two outbreaks shared a common 4b variant (4bV) serotype, also designated as IVb-v1, using a qPCR protocol developed in our laboratory. The 4bV serotype is characterized by the presence of a 6.3 Kb DNA segment normally found in serotype 1/2a, 3a, 1/2c and 3c strains but not in serotype 4b or 1/2b strains. We decided to compare these strains at a genomic level using the JSpecies Tetra tool. Specifically, we compared several 4bV and 4b isolates and identified a high level of similarity between the stone fruit and apple 4bV strains, but not the 4b strains co-identified in the caramel apple outbreak or other 4b or 4bV strains in our collection. This finding was further substantiated by a SNP-based analysis. Additionally, we were able to identify close relatedness between isolates from clinical cases from 1993–1994 and a single case from 2011 as well as links between two isolates from over 30 years ago. The identification of these potential links shows that JSpecies Tetra analysis can be a useful tool in rapidly assessing genetic relatedness of Lm isolates during outbreak investigations and for comparing historical isolates. Our analyses led to the identification of a highly related clonal group involved in two separate outbreaks, stone fruit and caramel apple, and suggests the possibility of a new genotype that may be better adapted for certain foods and/or environment.

Occurrence and Diversity of Clinically Important Vibrio Species in the Aquatic Environment of Georgia

Among the more than 70 different Vibrio species inhabiting marine, estuarine, and freshwater ecosystems, 12 are recognized as human pathogens. The warm subtropical climate of the Black Sea coastal area and inland regions of Georgia likely provides a favorable environment for various Vibrio species. From 2006 to 2009, the abundance, ecology, and diversity of clinically important Vibrio species were studied in different locations in Georgia and across seasons. Over a 33-month period, 1,595 presumptive Vibrio isolates were collected from the Black Sea (n = 657) and freshwater lakes around Tbilisi (n = 938). Screening of a subset of 440 concentrated and enriched water samples by PCR-electrospray ionization/mass spectrometry (PCR-ESI/MS) detected the presence of DNA from eight clinically important Vibrio species: V. cholerae, V. parahaemolyticus, V. vulnificus, V. mimicus, V. alginolyticus, V. harveyi, V. metschnikovii, and V. cincinnatiensis. Almost 90% of PCR/ESI-MS samples positive for Vibrio species were collected from June through November. Three important human-pathogenic Vibrio species (V. cholerae, V. parahaemolyticus, and V. vulnificus) were detected in 62.8, 37.8, and 21.4% of samples testing positive for Vibrios, respectively. The results of these activities suggest that natural reservoirs for human-pathogenic Vibrios exist in Georgian aquatic environments. Water temperature at all sampling sites was positively correlated with the abundance of clinically important Vibrio spp. (except V. metschnikovii), and salinity was correlated with species composition at particular Black Sea sites as well as inland reservoirs.

Cilantro microbiome before and after nonselective pre-enrichment for Salmonella using 16S rRNA and metagenomic sequencing

BACKGROUND

Salmonella enterica is a common cause of foodborne gastroenteritis in the United States and is associated with outbreaks in fresh produce such as cilantro. Salmonella culture-based detection methods are complex and time consuming, and improvments to increase detection sensitivity will benefit consumers. In this study, we used 16S rRNA sequencing to determine the microbiome of cilantro. We also investigated changes to the microbial community prior to and after a 24-hour nonselective pre-enrichment culture step commonly used by laboratory analysts to resuscitate microorganisms in foods suspected of contamination with pathogens. Cilantro samples were processed for Salmonella detection according to the method in the United States Food and Drug Administration Bacteriological Analytical Manual. Genomic DNA was extracted from culture supernatants prior to and after a 24-hour nonselective pre-enrichment step and 454 pyrosequencing was performed on 16S rRNA amplicon libraries. A database of Enterobacteriaceae 16S rRNA sequences was created, and used to screen the libraries for Salmonella, as some samples were known to be culture positive. Additionally, culture positive cilantro samples were examined for the presence of Salmonella using shotgun metagenomics on the Illumina MiSeq.

RESULTS

Time zero uncultured samples had an abundance of Proteobacteria while the 24-hour enriched samples were composed mostly of Gram-positive Firmicutes. Shotgun metagenomic sequencing of Salmonella culture positive cilantro samples revealed variable degrees of Salmonella contamination among the sequenced samples.

CONCLUSIONS

Our cilantro study demonstrates the use of high-throughput sequencing to reveal the microbiome of cilantro, and how the microbiome changes during the culture-based protocols employed by food safety laboratories to detect foodborne pathogens. Finding that culturing the cilantro shifts the microbiome to a predominance of Firmicutes suggests that changing our culture-based methods will improve detection sensitivity for foodborne enteric pathogens.

Development of a Custom-Designed, Pan Genomic DNA Microarray to Characterize Strain-Level Diversity among Cronobacter spp

Cronobacter species cause infections in all age groups; however neonates are at highest risk and remain the most susceptible age group for life-threatening invasive disease. The genus contains seven species:Cronobacter sakazakii, Cronobacter malonaticus, Cronobacter turicensis, Cronobacter muytjensii, Cronobacter dublinensis, Cronobacter universalis, and Cronobacter condimenti. Despite an abundance of published genomes of these species, genomics-based epidemiology of the genus is not well established. The gene content of a diverse group of 126 unique Cronobacter and taxonomically related isolates was determined using a pan genomic-based DNA microarray as a genotyping tool and as a means to identify outbreak isolates for food safety, environmental, and clinical surveillance purposes. The microarray constitutes 19,287 independent genes representing 15 Cronobacter genomes and 18 plasmids and 2,371 virulence factor genes of phylogenetically related Gram-negative bacteria. The Cronobacter microarray was able to distinguish the seven Cronobacter species from one another and from non-Cronobacter species; and within each species, strains grouped into distinct clusters based on their genomic diversity. These results also support the phylogenic divergence of the genus and clearly highlight the genomic diversity among each member of the genus. The current study establishes a powerful platform for further genomics research of this diverse genus, an important prerequisite toward the development of future countermeasures against this foodborne pathogen in the food safety and clinical arenas.

Re-examination of the taxonomic status of Enterobacter helveticus, Enterobacter pulveris and Enterobacter turicensis as members of the genus Cronobacter and their reclassification in the genera Franconibacter gen. nov. and Siccibacter gen. nov. as Franconibacter helveticus comb. nov., Franconibacter pulveris comb. nov. and Siccibacter turicensis comb. nov., respectively

Recently, a taxonomical re-evaluation of the genus Enterobacter, based on multi-locus sequence typing (MLST) analysis, has led to the proposal that the species Enterobacter pulveris, Enterobacter helveticus and Enterobacter turicensis should be reclassified as novel species of the genus Cronobacter. In the present work, new genome-scale analyses, including average nucleotide identity, genome-scale phylogeny and k-mer analysis, coupled with previously reported DNA–DNA hybridization values and biochemical characterization strongly indicate that these three species of the genus Enterobacter are not members of the genus Cronobacter, nor do they belong to the re-evaluated genus Enterobacter. Furthermore, data from this polyphasic study indicated that all three species constitute two new genera. We propose reclassifying Enterobacter pulveris and Enterobacter helveticus in the genus Franconibacter gen. nov. as Franconibacter pulveris comb. nov. (type strain 601/05^T = LMG 24057^T = DSM 19144^T) and Franconibacter helveticus comb. nov. (type strain 513/05^T = LMG 23732^T = DSM 18396^T), respectively, and Enterobacter turicensis in the genus Siccibacter gen. nov. as Siccibacter turicensis comb. nov. (type strain 508/05^T = LMG 23730^T = DSM 18397^T).

Genomic and phenotypic characterization of Vibrio cholerae non-O1 isolates from a US Gulf Coast cholera outbreak

Between November 2010, and May 2011, eleven cases of cholera, unrelated to a concurrent outbreak on the island of Hispaniola, were recorded, and the causative agent, Vibrio cholerae serogroup O75, was traced to oysters harvested from Apalachicola Bay, Florida. From the 11 diagnosed cases, eight isolates of V. cholerae were isolated and their genomes were sequenced. Genomic analysis demonstrated the presence of a suite of mobile elements previously shown to be involved in the disease process of cholera (ctxAB, VPI-1 and -2, and a VSP-II like variant) and a phylogenomic analysis showed the isolates to be sister taxa to toxigenic V. cholerae V51 serogroup O141, a clinical strain isolated 23 years earlier. Toxigenic V. cholerae O75 has been repeatedly isolated from clinical cases in the southeastern United States and toxigenic V. cholerae O141 isolates have been isolated globally from clinical cases over several decades. Comparative genomics, phenotypic analyses, and a Caenorhabditis elegans model of infection for the isolates were conducted. This analysis coupled with isolation data of V. cholerae O75 and O141 suggests these strains may represent an underappreciated clade of cholera-causing strains responsible for significant disease burden globally.

Molecular diversity and predictability of Vibrio parahaemolyticus along the Georgian coastal zone of the Black Sea

Vibrio parahaemolyticus is a leading cause of seafood-related gastroenteritis and is also an autochthonous member of marine and estuarine environments worldwide. One-hundred seventy strains of V. parahaemolyticus were isolated from water and plankton samples collected along the Georgian coast of the Black Sea during 28 months of sample collection. All isolated strains were tested for presence of tlh, trh, and tdh. A subset of strains were serotyped and tested for additional factors and markers of pandemicity. Twenty-six serotypes, five of which are clinically relevant, were identified. Although all 170 isolates were negative for tdh, trh, and the Kanagawa Phenomenon, 7 possessed the GS-PCR sequence and 27 the 850 bp sequence of V. parahaemolyticus pandemic strains. The V. parahaemolyticus population in the Black Sea was estimated to be genomically heterogeneous by rep-PCR and the serodiversity observed did not correlate with rep-PCR genomic diversity. Statistical modeling was used to predict presence of V. parahaemolyticus as a function of water temperature, with strongest concordance observed for Green Cape site samples (Percent of total variance = 70, P < 0.001). Results demonstrate a diverse population of V. parahaemolyticus in the Black Sea, some of which carry pandemic markers, with increased water temperature correlated to an increase in abundance of V. parahaemolyticus.

Cronobacter: an emergent pathogen causing meningitis to neonates through their feeds

The recognition of Cronobacter as a public health concern was raised when powdered infant formula (PIF) was linked to several neonatal meningitis outbreaks. It is an opportunistic pathogen that causes necrotising enterocolitis, infantile septicaemia, and meningitis which carries a high mortality rate among neonates. It has been also linked with cases of infection in adults and elderly. Over the past decade, much focus has been made on developing sensitive and specific characterisation, detection, and isolation methods to ascertain the quality of foods, notably contamination of PIF with Cronobacter and to understand its ability to cause disease. Whole genome sequencing has unveiled several putative virulence factors, yet the full capacity of the pathogenesis of Cronobacter has not yet been elucidated.

Identification of virulence properties in Salmonella Typhimurium DT104 using Caenorhabditis elegans

Salmonella enterica serover Typhimurium definitive phage type DT104, resistant to multiple antibiotics, is one of the most widespread Salmonella species in human infection worldwide. Although several cohort studies indicate that DT104 carrying the multidrug resistance (MDR) locus on salmonella genomic island 1 is a possible hyper-virulent strain compared to DT104 strains without MDR, or other Salmonella enterica serotypes, existing experimental evidence regarding virulence properties associated with the MDR region is controversial. To address this question, we constructed an isogenic MDR deletion (∆MDR) mutant strain of DT104, SNS12, by allelic exchange and used Caenorhabditis elegans as a host model to assess differences in virulence between these two strains. SNS12 exhibited decreased virulence in C. elegans, and we observed increased colonization and proliferation of the intestine of C. elegans by DT104. The immune response against MDR-carrying DT104 appears to function through a non-canonical Unfolded Protein Response (UPR) pathway, namely prion-like-(QN-rich)-domain-bearing protein pathway (PQN), in a ced-1 dependent manner in C. elegans. Further, we also demonstrate that genes of the PQN pathway and antimicrobial peptide gene abf-2, are expressed at higher transcriptional levels in worms immediately following exposure to DT104, in comparison with worms exposed to SNS12. Altogether, our results suggest that the MDR region of Salmonella Typhimurium DT104 has a direct role in virulence against Caenorhabditis elegans.

Complete genome sequence and phenotype microarray analysis of Cronobacter sakazakii SP291: a persistent isolate cultured from a powdered infant formula production facility

Outbreaks of human infection linked to the powdered infant formula (PIF) food chain and associated with the bacterium Cronobacter, are of concern to public health. These bacteria are regarded as opportunistic pathogens linked to life-threatening infections predominantly in neonates, with an under developed immune system. Monitoring the microbiological ecology of PIF production sites is an important step in attempting to limit the risk of contamination in the finished food product. Cronobacter species, like other microorganisms can adapt to the production environment. These organisms are known for their desiccation tolerance, a phenotype that can aid their survival in the production site and PIF itself. In evaluating the genome data currently available for Cronobacter species, no sequence information has been published describing a Cronobacter sakazakii isolate found to persist in a PIF production facility. Here we report on the complete genome sequence of one such isolate, Cronobacter sakazakii SP291 along with its phenotypic characteristics. The genome of C. sakazakii SP291 consists of a 4.3-Mb chromosome (56.9% GC) and three plasmids, denoted as pSP291-1, [118.1-kb (57.2% GC)], pSP291-2, [52.1-kb (49.2% GC)], and pSP291-3, [4.4-kb (54.0% GC)]. When C. sakazakii SP291 was compared to the reference C. sakazakii ATCC BAA-894, which is also of PIF origin, the annotated genome data identified two interesting functional categories, comprising of genes related to the bacterial stress response and resistance to antimicrobial and toxic compounds. Using a phenotypic microarray (PM), we provided a full metabolic profile comparing C. sakazakii SP291 and the previously sequenced C. sakazakii ATCC BAA-894. These data extend our understanding of the genome of this important neonatal pathogen and provides further insights into the genotypes associated with features that can contribute to its persistence in the PIF environment.

Pan-genome analysis of the emerging foodborne pathogen Cronobacter spp. suggests a species-level bidirectional divergence driven by niche adaptation

Background

Members of the genus Cronobacter are causes of rare but severe illness in neonates and preterm infants following the ingestion of contaminated infant formula. Seven species have been described and two of the species genomes were subsequently published. In this study, we performed comparative genomics on eight strains of Cronobacter, including six that we sequenced (representing six of the seven species) and two previously published, closed genomes.

Results

We identified and characterized the features associated with the core and pan genome of the genus Cronobacter in an attempt to understand the evolution of these bacteria and the genetic content of each species. We identified 84 genomic regions that are present in two or more Cronobacter genomes, along with 45 unique genomic regions. Many potentially horizontally transferred genes, such as lysogenic prophages, were also identified. Most notable among these were several type six secretion system gene clusters, transposons that carried tellurium, copper and/or silver resistance genes, and a novel integrative conjugative element.

Conclusions

Cronobacter have diverged into two clusters, one consisting of C. dublinensis and C. muytjensii (Cdub-Cmuy) and the other comprised of C. sakazakii, C. malonaticus, C. universalis, and C. turicensis, (Csak-Cmal-Cuni-Ctur) from the most recent common ancestral species. While several genetic determinants for plant-association and human virulence could be found in the core genome of Cronobacter, the four Cdub-Cmuy clade genomes contained several accessory genomic regions important for survival in a plant-associated environmental niche, while the Csak-Cmal-Cuni-Ctur clade genomes harbored numerous virulence-related genetic traits.

Role of GbpA protein, an important virulence-related colonization factor, for Vibrio cholerae's survival in the aquatic environment

Vibrio cholerae N-acetyl glucosamine-binding protein A (GbpA) is a chitin binding protein and a virulence factor involved in the colonization of human intestine. We investigated the distribution and genetic variations of gbpA in 488 V. cholerae strains of environmental and clinical origin, belonging to different serogroups and biotypes. We found that the gene is consistently present and highly conserved including an environmental V. cholerae-related strain of ancestral origin. The gene was also consistently expressed in a number of representative V. cholerae strains cultured in laboratory aquatic microcosms under conditions simulating those found in temperate marine environments. Functional analysis carried out on V. cholerae O1 El Tor N16961 showed that GbpA is not involved in adhesion to inorganic surfaces but promotes interaction with environmental biotic substrates (plankton and bivalve hepatopancreas cells) representing known marine reservoir or host for the bacterium. It is suggested that the ability of GbpA to colonize human intestinal cells most probably originated from its primary function in the aquatic environment.

Vibrio cholerae in an Historically Cholera-Free Country

We report the autochthonous existence of Vibrio cholerae in coastal waters of Iceland, a geothermally active country where cholera is absent and has never been reported. Seawater, mussel, and macroalgae samples were collected close to and distant from sites where geothermal activity causes a significant increase in water temperature during low tides. V. cholerae was detected only at geothermal-influenced sites during low-tides. None of the V. cholerae isolates encoded cholera toxin (ctxAB) and all were non-O1/non-O139 serogroups. However, all isolates encoded other virulence factors that are associated with cholera as well as extra-intestinal V. cholerae infections. The virulence factors were functional at temperatures of coastal waters of Iceland, suggesting an ecological role. It is noteworthy that V. cholerae was isolated from samples collected at sites distant from anthropogenic influence, supporting the conclusion that V. cholerae is autochthonous to the aquatic environment of Iceland.

Interaction of Vibrio cholerae non-O1/non-O139 with copepods, cladocerans and competing bacteria in the large alkaline lake Neusiedler See, Austria

Vibrio cholerae is a human pathogen and natural inhabitant of aquatic environments. Serogroups O1/O139 have been associated with epidemic cholera, while non-O1/non-O139 serogroups usually cause human disease other than classical cholera. V. cholerae non-O1/non-O139 from the Neusiedler See, a large Central European lake, have caused ear and wound infections, including one case of fatal septicaemia. Recent investigations demonstrated rapid planktonic growth of V. cholerae non-O1/non-O139 and correlation with zooplankton biomass. The aim of this study was to elucidate the interaction of autochthonous V. cholerae with two dominant crustacean zooplankton species in the lake and investigate the influence of the natural bacterial community on this interaction. An existing data set was evaluated for statistical relationships between zooplankton species and V. cholerae and co-culture experiments were performed in the laboratory. A new fluorescence in situ hybridisation protocol was applied for quantification of V. cholerae non-O1/non-O139 cells, which significantly reduced analysis time. The experiments clearly demonstrated a significant relationship of autochthonous V. cholerae non-O1/non-O139 with cladocerans by promoting growth of V. cholerae non-O1/non-O139 in the water and on the surfaces of the cladocerans. In contrast, copepods had a negative effect on the growth of V. cholerae non-O1/non-O139 via competing bacteria from their surfaces. Thus, beside other known factors, biofilm formation by V. cholerae on crustacean zooplankton appears to be zooplankton taxon specific and may be controlled by the natural bacterial community.

Occurrence of the Vibrio cholerae seventh pandemic VSP-I island and a new variant

Using comparative genomics, we identified a new variant of the Vibrio Seventh Pandemic Island-I (VSP-I). Results of polymerase chain reaction (PCR) screening for both known VSP-I variants indicate that the novel variant is present only in non-O1/non-O139 strains of V. cholerae and Vibrio mimicus. Comparative genomics revealed little sequence divergence in the seventh pandemic VSP-I; however, a second insertion site located on the smaller chromosome was identified. Although the seventh pandemic VSP-I genomic island was detected in all seventh pandemic V. cholerae serogroup O1 and O139 isolates examined in this study, unique genes of the island cannot be used alone as an identifying target, because the seventh pandemic VSP-I was also present in three non-seventh pandemic strains of V. cholerae isolated from Chesapeake Bay. As an alternative, a PCR assay targeting the VC2346 gene was found to be confirmatory for seventh pandemic isolates of V. cholerae.

Comparative genomic analysis reveals evidence of two novel Vibrio species closely related to V. cholerae

Background

In recent years genome sequencing has been used to characterize new bacterial species, a method of analysis available as a result of improved methodology and reduced cost. Included in a constantly expanding list of Vibrio species are several that have been reclassified as novel members of the Vibrionaceae. The description of two putative new Vibrio species, Vibrio sp. RC341 and Vibrio sp. RC586 for which we propose the names V. metecus and V. parilis, respectively, previously characterized as non-toxigenic environmental variants of V. cholerae is presented in this study.

Results

Based on results of whole-genome average nucleotide identity (ANI), average amino acid identity (AAI), rpoB similarity, MLSA, and phylogenetic analysis, the new species are concluded to be phylogenetically closely related to V. cholerae and V. mimicus. Vibrio sp. RC341 and Vibrio sp. RC586 demonstrate features characteristic of V. cholerae and V. mimicus, respectively, on differential and selective media, but their genomes show a 12 to 15% divergence (88 to 85% ANI and 92 to 91% AAI) compared to the sequences of V. cholerae and V. mimicus genomes (ANI <95% and AAI <96% indicative of separate species). Vibrio sp. RC341 and Vibrio sp. RC586 share 2104 ORFs (59%) and 2058 ORFs (56%) with the published core genome of V. cholerae and 2956 (82%) and 3048 ORFs (84%) with V. mimicus MB-451, respectively. The novel species share 2926 ORFs with each other (81% Vibrio sp. RC341 and 81% Vibrio sp. RC586). Virulence-associated factors and genomic islands of V. cholerae and V. mimicus, including VSP-I and II, were found in these environmental Vibrio spp.

Conclusions

Results of this analysis demonstrate these two environmental vibrios, previously characterized as variant V. cholerae strains, are new species which have evolved from ancestral lineages of the V. cholerae and V. mimicus clade. The presence of conserved integration loci for genomic islands as well as evidence of horizontal gene transfer between these two new species, V. cholerae, and V. mimicus suggests genomic islands and virulence factors are transferred between these species.

Discovery of novel Vibrio cholerae VSP-II genomic islands using comparative genomic analysis

This report describes Vibrio seventh pandemic island II (VSP-II) and three novel variants revealed by comparative genomics of 23 Vibrio cholerae strains and their presence among a large and diverse collection of V. cholerae isolates. Three VSP-II variants were reported previously and our results demonstrate the presence of three novel VSP-II in clinical and environmental V. cholerae marked by major deletions and genetic rearrangements. A new VSP-II cluster was found in the seventh pandemic V. cholerae O1 El Tor strain CIRS101, which is dominant (95%) among the recent (2004-2007) seven pandemic V. cholerae O1 El Tor isolates from two endemic sites, but was not found in older strains from the same region. Two other variants were found in V. cholerae TMA21 and RC385, two environmental strains from coastal Brazil and the Chesapeake Bay, respectively, the latter being prevalent among environmental V. cholerae non-O1/non-O139 and Vibrio mimicus. The results of this study indicate that the VSP-II island has undergone significant rearrangement through a complex evolutionary pathway in V. cholerae. Interestingly, one of the new VSP-II revealed the presence of 'old' and 'new' V. cholerae O1 El Tor pandemic clones circulating in some of the areas where cholera is endemic.

Detection of toxigenic Vibrio cholerae O1 in freshwater lakes of the former Soviet Republic of Georgia

Three freshwater lakes, Lisi Lake, Kumisi Lake and Tbilisi Sea, near Tbilisi, Georgia, were studied from January 2006 to December 2007 to determine the presence of Vibrio cholerae employing both bacteriological culture method and direct detection methods, namely PCR and direct fluorescent antibody (DFA). For PCR, DNA extracted from water samples was tested for presence of V. cholerae and genes coding for selected virulence factors. Vibrio cholerae non-O1/non-O139 was routinely isolated by culture from all three lakes; whereas V. cholerae O1 and O139 were not. Water samples collected during the summer months from Lisi Lake and Kumisi Lake were positive for both V. cholerae and V. cholerae ctxA, tcpA, zot, ompU and toxR by PCR. Water samples collected during the same period from both Lisi and Kumisi Lake were also positive for V. cholerae serogroup O1 by DFA. All of the samples were negative for V. cholerae serotype O139. The results of this study provide evidence for an environmental presence of toxigenic V. cholerae O1, which may represent a potential source of illness as these lakes serve as recreational water in Tbilisi, Georgia.

The pre-seventh pandemic Vibrio cholerae BX 330286 El Tor genome: evidence for the environment as a genome reservoir

Vibrio cholerae O1 El Tor BX 330286 was isolated from a water sample in Australia in 1986, 9 years after an indigenous outbreak of cholera occurred in that region. This environmental strain encodes virulence factors highly similar to those of clinical strains, suggesting an ability to cause disease in humans. We demonstrate its high similarity in gene content and genome-wide nucleotide sequence to clinical V. cholerae strains, notably to pre-seventh pandemic O1 El Tor strains isolated in 1910 (V. cholerae NCTC 8457) and 1937 (V. cholerae MAK 757), as well as seventh pandemic strains isolated after 1960 globally. Here we demonstrate that this strain represents a transitory clone with shared characteristics between pre-seventh and seventh pandemic strains of V. cholerae. Interestingly, this strain was isolated 25 years after the beginning of the seventh pandemic, suggesting the environment as a genome reservoir in areas where cholera does not occur in sporadic, endemic or epidemic form.

Genomic analysis of a novel integrative conjugative element in Vibrio cholerae

Integrative conjugative elements (ICEs) are a class of self-transmissible mobile elements that mediate horizontal gene transfer in bacteria, and play an important role in bacterial evolution. Since 1992, ICEs of the SXT/R391 family have been found to be widely distributed among Vibrio cholerae strains isolated in Asian countries. Here we describe ICEVchB33, an ICE found in the genomes of two V. cholerae O1 Eltor strains, one isolated in India, 1994, and the other from Mozambique, 2004. ICEVchB33 revealed a new genetic organization, different from other ICEs of the SXT/R391 family, demonstrating the genomic plasticity of these elements.

Biofilms in water, its role and impact in human disease transmission

Understanding the mechanism of biofilm formation is the first step in determining its function and, thereby, its impact and role in the environment. Extensive studies accomplished during the past few years have elucidated the genetics and biochemistry of biofilm formation. Cell-to-cell communication, that is, quorum sensing, is a key factor in the initiation of biofilm. Occurrence of viable but nonculturable bacteria, including Vibrio cholerae in biofilms has been reported and most likely such cells were overlooked previously because appropriate methods of detection were not employed. For this reason discovery and investigation of this important bacterial ecological niche in the environment were impeded.