To date 142 species have been described in the Vibrionaceae family of bacteria, classified into seven genera; Aliivibrio, Echinimonas, Enterovibrio, Grimontia, Photobacterium, Salinivibrio and Vibrio. As vibrios are widespread in marine environments and show versatile metabolisms and ecologies, these bacteria are recognized as one of the most diverse and important marine heterotrophic bacterial groups for elucidating the correlation between genome evolution and ecological adaptation. However, on the basis of 16S rRNA gene phylogeny, we could not find any robust monophyletic lineages in any of the known genera. We needed further attempts to reconstruct their evolutionary history based on multilocus sequence analysis (MLSA) and/or genome wide taxonomy of all the recognized species groups. In our previous report in 2007, we conducted the first broad multilocus sequence analysis (MLSA) to infer the evolutionary history of vibrios using nine housekeeping genes (the 16S rRNA gene, gapA, gyrB, ftsZ, mreB, pyrH, recA, rpoA, and topA), and we proposed 14 distinct clades in 58 species of Vibrionaceae. Due to the difficulty of designing universal primers that can amplify the genes for MLSA in every Vibrionaceae species, some clades had yet to be defined. In this study, we present a better picture of an updated molecular phylogeny for 86 described vibrio species and 10 genome sequenced Vibrionaceae strains, using 8 housekeeping gene sequences. This new study places special emphasis on (1) eight newly identified clades (Damselae, Mediterranei, Pectenicida, Phosphoreum, Profundum, Porteresiae, Rosenbergii, and Rumoiensis); (2) clades amended since the 2007 proposal with recently described new species; (3) orphan clades of genomospecies F6 and F10; (4) phylogenetic positions defined in 3 genome-sequenced strains (N418, EX25, and EJY3); and (5) description of V. tritonius sp. nov., which is a member of the “Porteresiae” clade.

In bacterial genome and metagenome sequencing, Illumina sequencers are most frequently used due to their high throughput capacity, and multiple library preparation kits have been developed for Illumina platforms. Here, we systematically analysed and compared the sequencing bias generated by currently available library preparation kits for Illumina sequencing. Our analyses revealed that a strong sequencing bias is introduced in low-GC regions by the Nextera XT kit. The level of bias introduced is dependent on the level of GC content; stronger bias is generated as the GC content decreases. Other analysed kits did not introduce this strong sequencing bias. The GC content-associated sequencing bias introduced by Nextera XT was more remarkable in metagenome sequencing of a mock bacterial community and seriously affected estimation of the relative abundance of low-GC species. The results of our analyses highlight the importance of selecting proper library preparation kits according to the purposes and targets of sequencing, particularly in metagenome sequencing, where a wide range of microbial species with various degrees of GC content is present. Our data also indicate that special attention should be paid to which library preparation kit was used when analysing and interpreting publicly available metagenomic data.

Enterohemorrhagic E. coli (EHEC) causes diarrhea and hemorrhagic colitis with life-threatening complications, such as hemolytic uremic syndrome. Their major virulence factor is Shiga toxin (Stx), which is encoded by bacteriophages. Of the two types of Stx, the production of Stx2, particularly that of Stx2a (a subtype of Stx2), is a major risk factor for severe EHEC infections, but the Stx2 production level is highly variable between strains. Here, we define four major and two minor subtypes of Stx2a-encoding phages according to their replication proteins. The subtypes are correlated with Stx2a titers produced by the host O157 strains, suggesting a critical role of the phage subtype in determining the Stx2a production level. We further show that one of the two subclades in the clade 8, a proposed hyper-virulent lineage of O157, carries the Stx2 phage subtype that confers the highest Stx2 production to the host strain. The presence of this subclade may explain the proposed high virulence potential of clade 8. These results provide novel insights into the variation in virulence among O157 strains and highlight the role of phage variation in determining the production level of the virulence factors that phages encode.

Microorganisms often live in symbiosis with their hosts, and some are considered mutualists, where all species involved benefit from the interaction. How free-living microorganisms have evolved to become mutualists is unclear. Here we report an experimental system in which non-symbiotic Escherichia coli evolves into an insect mutualist. The stinkbug Plautia stali is typically associated with its essential gut symbiont, Pantoea sp., which colonizes a specialized symbiotic organ. When sterilized newborn nymphs were infected with E. coli rather than Pantoea sp., only a few insects survived, in which E. coli exhibited specific localization to the symbiotic organ and vertical transmission to the offspring. Through transgenerational maintenance with P. stali, several hypermutating E. coli lines independently evolved to support the host’s high adult emergence and improved body colour; these were called ‘mutualistic’ E. coli. These mutants exhibited slower bacterial growth, smaller size, loss of flagellar motility and lack of an extracellular matrix. Transcriptomic and genomic analyses of ‘mutualistic’ E. coli lines revealed independent mutations that disrupted the carbon catabolite repression global transcriptional regulator system. Each mutation reproduced the mutualistic phenotypes when introduced into wild-type E. coli, confirming that single carbon catabolite repression mutations can make E. coli an insect mutualist. These findings provide an experimental system for future work on host–microbe symbioses and may explain why microbial mutualisms are omnipresent in nature.

Escherichia albertii is a recently recognized close relative of Escherichia coli . This emerging enteropathogen possesses a type III secretion system (T3SS) encoded by the locus of enterocyte effacement, similar to enteropathogenic and enterohemorrhagic E. coli (EPEC and EHEC). Shiga toxin-producing strains have also been identified. The genomic features of E. albertii , particularly differences from other Escherichia species, have not yet been well clarified. Here, we sequenced the genome of 29 E. albertii strains (3 complete and 26 draft sequences) isolated from multiple sources and performed intraspecies and intragenus genomic comparisons. The sizes of the E. albertii genomes range from 4.5 to 5.1 Mb, smaller than those of E. coli strains. Intraspecies genomic comparisons identified five phylogroups of E. albertii . Intragenus genomic comparison revealed that the possible core genome of E. albertii comprises 3,250 genes, whereas that of the genus Escherichia comprises 1,345 genes. Our analysis further revealed several unique or notable genetic features of E. albertii , including those responsible for known biochemical features and virulence factors and a possibly active second T3SS known as ETT2 ( E. coli T3SS 2) that is inactivated in E. coli . Although this organism has been observed to be nonmotile in vitro, genes for flagellar biosynthesis are fully conserved; chemotaxis-related genes have been selectively deleted. Based on these results, we have developed a nested polymerase chain reaction system to directly detect E. albertii . Our data define the genomic features of E. albertii and provide a valuable basis for future studies of this important emerging enteropathogen.

Serratia marcescens is an important nosocomial pathogen causing various opportunistic infections, such as urinary tract infections, bacteremia and sometimes even hospital outbreaks. The recent emergence and spread of multidrug-resistant (MDR) strains further pose serious threats to global public health. This bacterium is also ubiquitously found in natural environments, but the genomic differences between clinical and environmental isolates are not clear, including those between S. marcescens and its close relatives. In this study, we performed a large-scale genome analysis of S. marcescens and closely related species (referred to as the ‘ S. marcescens complex’), including more than 200 clinical and environmental strains newly sequenced here. Our analysis revealed their phylogenetic relationships and complex global population structure, comprising 14 clades, which were defined based on whole-genome average nucleotide identity. Clades 10, 11, 12 and 13 corresponded to S. nematodiphila , S. marcescens sensu stricto, S. ureilytica and S. surfactantfaciens, respectively. Several clades exhibited distinct genome sizes and GC contents and a negative correlation of these genomic parameters was observed in each clade, which was associated with the acquisition of mobile genetic elements (MGEs), but different types of MGEs, plasmids or prophages (and other integrative elements), were found to contribute to the generation of these genomic variations. Importantly, clades 1 and 2 mostly comprised clinical or hospital environment isolates and accumulated a wide range of antimicrobial resistance genes, including various extended-spectrum β-lactamase and carbapenemase genes, and fluoroquinolone target site mutations, leading to a high proportion of MDR strains. This finding suggests that clades 1 and 2 represent hospital-adapted lineages in the S. marcescens complex although their potential virulence is currently unknown. These data provide an important genomic basis for reconsidering the classification of this group of bacteria and reveal novel insights into their evolution, biology and differential importance in clinical settings.