Diversity and evolution of surface polysaccharide synthesis loci in Enterobacteriales

Two novel Enterobacter species, Enterobacter chinensis sp. nov. and Enterobacter rongchengensis sp. nov., recovered from clinical samples carrying multiple virulence factors

Two Enterobacter strains 170198T and 170250T were isolated from clinical blood samples from distinct patients in a hospital in Chengdu, China, in 2022. These isolates were subjected to whole-genome sequencing. A phylogenomic tree based on 2,096 concatenated core genes showed that the two strains were clustered within the genus Enterobacter. The average nucleotide identity (ANI) and in silico DNA-DNA hybridization (isDDH) values between each of the two strains and type strains of all currently known Enterobacter species were determined. The two strains belonged to two novel species as the highest ANI and isDDH values with type strains of all currently known Enterobacter species below the cutoff for species demarcation (96% for ANI and 70% for isDDH). Then the physiological and biochemical studies demonstrated that biochemical features and the profile of whole fatty acids of strains 170198T and 170250T were largely consistent with those known Enterobacter species. Nevertheless, the two novel species can be differentiated from all other Enterobacter species by certain biochemical characteristics. In conclusion, 170198T and 170250T represent two novel species of the genus Enterobacter, for which we propose Enterobacter chinensis sp. nov. and Enterobacter rongchengensis sp. nov., as the species names. The type strains of Enterobacter chinensis sp. nov., and Enterobacter rongchengensis sp. nov. are 170198T (=GDMCC 1.3549T=JCM 35826T) and 170250T (=GDMCC 1.3670T=JCM 36189T), respectively. The two novel species have clinical significance with the ability to cause bloodstream infections.IMPORTANCEEnterobacter is a group of bacteria comprising several common opportunistic pathogens and has a complicated taxonomy. Here, we reported two novel Enterobacter species. We demonstrated that the two novel species can be differentiated from other Enterobacter species by certain phenotypic characteristics and therefore provide information for designing tests for identification. We also showed that strains of the two novel species are able to cause human bloodstream infections and carry multiple virulence factors and therefore are of clinical significance. We highlight that the virulence of Enterobacter is less studied and warrants further exploration. We believe that the findings here are valuable for enhancing the appreciation toward Enterobacter, an important pathogen.

Genome sequence-based species classification of Enterobacter cloacae complex: a study among clinical isolates

Accurate species-level identification of Enterobacter cloacae complex (ECC) is crucial for related research. The classification of ECC is based on strain-to-strain phylogenetic congruence, as well as genomic features including average nucleotide identity (ANI) and digitalized DNA-DNA hybridization (dDDH). ANI and dDDH derived from whole-genome sequencing have emerged as a reliable metric for assessing genetic relatedness between genomes and are increasingly recognized as a standard for species delimitation. Up to now, there are two different classification methods for ECC. The first one categorizes E. hormaechei, a species within ECC, into five subspecies (E. hormaechei subsp. steigerwaltii, subsp. oharae, subsp. xiangfangensis, subsp. hoffmannii, and subsp. hormaechei). The second classifies E. hormaechei as three species: E. hormaechei, "E. xiangfangensis," "E. hoffmanii." While the former is well-accepted in the academic area, the latter may have a greater ability to distinguish different species of ECC. To assess the suitability of these identification criteria for clinical ECC isolates, we conducted a comprehensive analysis involving phylogenetic analysis, ANI and dDDH value alignment, virulence gene identification, and capsule typing on 256 clinical ECC strains isolated from the bloodstream. Our findings indicated that the method of categorizing E. hormaechei into five subspecies has better correlation and consistency with the molecular characteristics of clinical ECC isolates, as evidenced by phylogenetic analysis, virulence genes, and capsule typing. Therefore, the subspecies-based classification method appears more suitable for taxonomic assignments of clinical ECC isolates.

IMPORTANCE

Standardizing taxonomy of the Enterobacter cloacae complex (ECC) is necessary for data integration across diverse studies. The study utilized whole-genome data to accurately identify 256 clinical ECC isolated from bloodstream infections using average nucleotide identity (ANI), digitalized DNA-DNA hybridization (dDDH), and phylogenetic analysis. Through comprehensive assessments including phylogenetic analysis, ANI and dDDH comparisons, virulence gene, and capsule typing of the 256 clinical isolates, it was concluded that the classification method based on subspecies exhibited better correlation and consistency with the molecular characteristics of clinical ECC isolates. In summary, this research contributes to the precise identification of clinical ECC at the species level and expands our understanding of ECC.

Bacterial capsules: Occurrence, mechanism, and function

In environments characterized by extended multi-stress conditions, pathogens develop a variety of immune escape mechanisms to enhance their ability to infect the host. The capsules, polymers that bacteria secrete near their cell wall, participates in numerous bacterial life processes and plays a crucial role in resisting host immune attacks and adapting to their niche. Here, we discuss the relationship between capsules and bacterial virulence, summarizing the molecular mechanisms of capsular regulation and pathogenesis to provide new insights into the research on the pathogenesis of pathogenic bacteria.

Large-scale genomic survey with deep learning-based method reveals strain-level phage specificity determinants

BACKGROUND

Phage therapy, reemerging as a promising approach to counter antimicrobial-resistant infections, relies on a comprehensive understanding of the specificity of individual phages. Yet the significant diversity within phage populations presents a considerable challenge. Currently, there is a notable lack of tools designed for large-scale characterization of phage receptor-binding proteins, which are crucial in determining the phage host range.

RESULTS

In this study, we present SpikeHunter, a deep learning method based on the ESM-2 protein language model. With SpikeHunter, we identified 231,965 diverse phage-encoded tailspike proteins, a crucial determinant of phage specificity that targets bacterial polysaccharide receptors, across 787,566 bacterial genomes from 5 virulent, antibiotic-resistant pathogens. Notably, 86.60% (143,200) of these proteins exhibited strong associations with specific bacterial polysaccharides. We discovered that phages with identical tailspike proteins can infect different bacterial species with similar polysaccharide receptors, underscoring the pivotal role of tailspike proteins in determining host range. The specificity is mainly attributed to the protein's C-terminal domain, which strictly correlates with host specificity during domain swapping in tailspike proteins. Importantly, our dataset-driven predictions of phage-host specificity closely match the phage-host pairs observed in real-world phage therapy cases we studied.

CONCLUSIONS

Our research provides a rich resource, including both the method and a database derived from a large-scale genomics survey. This substantially enhances understanding of phage specificity determinants at the strain level and offers a valuable framework for guiding phage selection in therapeutic applications.

Shigella flexneri remodeling and consumption of host lipids during infection

IMPORTANCE

Bacterial pathogens have vastly distinct sites that they inhabit during infection. This requires adaptation due to changes in nutrient availability and antimicrobial stress. The bacterial surface is a primary barrier, and here, we show that the bacterial pathogen Shigella flexneri increases its surface decorations when it transitions to an intracellular lifestyle. We also observed changes in bacterial and host cell fatty acid homeostasis. Specifically, intracellular S. flexneri increased the expression of their fatty acid degradation pathway, while the host cell lipid pool was significantly depleted. Importantly, bacterial proliferation could be inhibited by fatty acid supplementation of host cells, thereby providing novel insights into the possible link between human malnutrition and susceptibility to S. flexneri.

Repeat-Unit Elongations To Produce Bacterial Complex Long Polysaccharide Chains, an O-Antigen Perspective

The O-antigen, a long polysaccharide that constitutes the distal part of the outer membrane-anchored lipopolysaccharide, is one of the critical components in the protective outer membrane of Gram-negative bacteria. Most species produce one of the structurally diverse O-antigens, with nearly all the polysaccharide components having complex structures made by the Wzx/Wzy pathway. This pathway produces repeat-units of mostly 3-8 sugars on the cytosolic face of the cytoplasmic membrane that is translocated by Wzx flippase to the periplasmic face and polymerized by Wzy polymerase to give long-chain polysaccharides. The Wzy polymerase is a highly diverse integral membrane protein typically containing 10-14 transmembrane segments. Biochemical evidence confirmed that Wzy polymerase is the sole driver of polymerization, and recent progress also began to demystify its interacting partner, Wzz, shedding some light to speculate how the proteins may operate together during polysaccharide biogenesis. However, our knowledge of how the highly variable Wzy proteins work as part of the O-antigen processing machinery remains poor. Here, we discuss the progress to the current understanding of repeat-unit polymerization and propose an updated model to explain the formation of additional short chain O-antigen polymers found in the lipopolysaccharide of diverse Gram-negative species and their importance in the biosynthetic process.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

Bacterial genome-wide association study substantiates papGII of Escherichia coli as a major risk factor for urosepsis

BACKGROUND

Urinary tract infections (UTIs) are among the most common bacterial infections worldwide, often caused by uropathogenic Escherichia coli. Multiple bacterial virulence factors or patient characteristics have been linked separately to progressive, more invasive infections. In this study, we aim to identify pathogen- and patient-specific factors that drive the progression to urosepsis by jointly analysing bacterial and host characteristics.

METHODS

We analysed 1076 E. coli strains isolated from 825 clinical cases with UTI and/or bacteraemia by whole-genome sequencing (Illumina). Sequence types (STs) were determined via srst2 and capsule loci via fastKaptive. We compared the isolates from urine and blood to confirm clonality. Furthermore, we performed a bacterial genome-wide association study (bGWAS) (pyseer) using bacteraemia as the primary clinical outcome. Clinical data were collected by an electronic patient chart review. We concurrently analysed the association of the most significant bGWAS hit and important patient characteristics with the clinical endpoint bacteraemia using a generalised linear model (GLM). Finally, we designed qPCR primers and probes to detect papGII-positive E. coli strains and prospectively screened E. coli from urine samples (n = 1657) at two healthcare centres.

RESULTS

Our patient cohort had a median age of 75.3 years (range: 18.00-103.1) and was predominantly female (574/825, 69.6%). The bacterial phylogroups B2 (60.6%; 500/825) and D (16.6%; 137/825), which are associated with extraintestinal infections, represent the majority of the strains in our collection, many of which encode a polysaccharide capsule (63.4%; 525/825). The most frequently observed STs were ST131 (12.7%; 105/825), ST69 (11.0%; 91/825), and ST73 (10.2%; 84/825). Of interest, in 12.3% (13/106) of cases, the E. coli pairs in urine and blood were only distantly related. In line with previous bGWAS studies, we identified the gene papGII (p-value < 0.001), which encodes the adhesin subunit of the E. coli P-pilus, to be associated with 'bacteraemia' in our bGWAS. In our GLM, correcting for patient characteristics, papGII remained highly significant (odds ratio = 5.27, 95% confidence interval = [3.48, 7.97], p-value < 0.001). An independent cohort of cases which we screened for papGII-carrying E. coli at two healthcare centres further confirmed the increased relative frequency of papGII-positive strains causing invasive infection, compared to papGII-negative strains (p-value = 0.033, chi-squared test).

CONCLUSIONS

This study builds on previous work linking papGII with invasive infection by showing that it is a major risk factor for progression from UTI to bacteraemia that has diagnostic potential.

Extensive Diversity in Escherichia coli Group 3 Capsules Is Driven by Recombination and Plasmid Transfer from Multiple Species

Bacterial capsules provide protection against environmental challenges and host immunity. Historically, Escherichia coli K serotyping scheme, which relies on the hypervariable capsules, has identified around 80 K forms that fall into four distinct groups. Based on recent work by us and others, we predicted that E. coli capsular diversity is grossly underestimated. We exploited group 3 capsule gene clusters, the best genetically defined capsule group in E. coli, to analyze publicly available E. coli sequences for overlooked capsular diversity within the species. We report the discovery of seven novel group 3 clusters that fall into two distinct subgroups (3A and 3B). The majority of the 3B capsule clusters were found on plasmids, contrary to the defining feature of group 3 capsule genes localizing at the serA locus on the E. coli chromosome. Other new group 3 capsule clusters were derived from ancestral sequences through recombination events between shared genes found within the serotype variable central region 2. Intriguingly, flanking regions 1 and 3, known to be conserved areas among capsule clusters, showed considerable intra-subgroup variation in clusters from the 3B subgroup, containing genes of shared ancestry with other Enterobacteriaceae species. Variation of group 3 kps clusters within dominant E. coli lineages, including multidrug-resistant pathogenic lineages, further supports that E. coli capsules are undergoing rigorous change. Given the pivotal role of capsular polysaccharides in phage predation, our findings raise attention to the need of monitoring kps evolutionary dynamics in pathogenic E. coli in supporting phage therapy. IMPORTANCE Capsular polysaccharides protect pathogenic bacteria against environmental challenges, host immunity, and phage predations. The historical Escherichia coli K typing scheme, which relies on the hypervariable capsular polysaccharide, has identified around 80 different K forms that fall into four distinct groups. Taking advantage of the supposedly compact and genetically well-defined group 3 gene clusters, we analyzed published E. coli sequences to identify seven new gene clusters and revealed an unexpected capsular diversity. Genetic analysis revealed that group 3 gene clusters shared closely related serotype-specific region 2 and were diversified through recombination events and plasmid transfer between multiple Enterobacteriaceae species. Overall, capsular polysaccharides in E. coli are undergoing rigorous change. Given the pivotal role capsules play in phage interactions, this work highlighted the need to monitor the evolutionary dynamics of capsules in pathogenic E. coli for effective phage therapy.

Shigella flexneri Adapts to Niche-Specific Stresses through Modifications in Cell Envelope Composition and Decoration

Shigella flexneri is the primary causative agent of worldwide shigellosis. As the pathogen transverses the distinct niches of the gastrointestinal tract it necessitates dynamic adaptation strategies to mitigate host antimicrobials such as dietary fatty acids (FAs) and the bile salt, deoxycholate (DOC). This study investigates the dynamics of the S. flexneri cell envelope, by interrogating adaptations following FA or DOC exposure. We deciphered the effects of FAs and DOC on bacterial membrane fatty acid and lipopolysaccharide (LPS) compositions. We identified novel LPS-based strategies by the pathogen to support resistance to these host compounds. In particular, expression of S. flexneri very-long O antigen (VL-Oag) LPS was found to play a central role in stress mitigation, as VL-Oag protects against antimicrobial FAs, but its presence rendered S. flexneri susceptible to DOC stress. Collectively, this work underpins the importance for S. flexneri to maintain appropriate regulation of cell envelope constituents, in particular VL-Oag LPS, to adequately adapt to diverse stresses during infection.

Hypermucoviscosity Regulator RmpD Interacts with Wzc and Controls Capsular Polysaccharide Chain Length

Klebsiella pneumoniae is a leading cause of nosocomial infections, including pneumonia, bacteremia, and urinary tract infections. Treatment options are increasingly restricted by the high prevalence of resistance to frontline antibiotics, including carbapenems, and the recently identified plasmid-conferred colistin resistance. The classical pathotype (cKp) is responsible for most of the nosocomial infections observed globally, and these isolates are often multidrug resistant. The hypervirulent pathotype (hvKp) is a primary pathogen capable of causing community-acquired infections in immunocompetent hosts. The hypermucoviscosity (HMV) phenotype is strongly associated with the increased virulence of hvKp isolates. Recent studies demonstrated that HMV requires capsule (CPS) synthesis and the small protein RmpD but is not dependent on the increased amount of capsule associated with hvKp. Here, we identified the structure of the capsular and extracellular polysaccharide isolated from hvKp strain KPPR1S (serotype K2) with and without RmpD. We found that the polymer repeat unit structure is the same in both strains and that it is identical to the K2 capsule. However, the chain length of CPS produced by strains expressing rmpD demonstrates more uniform length. This property was reconstituted in CPS from Escherichia coli isolates that possess the same CPS biosynthesis pathway as K. pneumoniae but naturally lack rmpD. Furthermore, we demonstrate that RmpD binds Wzc, a conserved capsule biosynthesis protein required for CPS polymerization and export. Based on these observations, we present a model for how the interaction of RmpD with Wzc could impact CPS chain length and HMV. IMPORTANCE Infections caused by Klebsiella pneumoniae continue to be a global public health threat; the treatment of these infections is complicated by the high frequency of multidrug resistance. K. pneumoniae produces a polysaccharide capsule required for virulence. Hypervirulent isolates also have a hypermucoviscous (HMV) phenotype that increases virulence, and we recently demonstrated that a horizontally acquired gene, rmpD, is required for HMV and hypervirulence but that the identity of the polymeric product(s) in HMV isolates is uncertain. Here, we demonstrate that RmpD regulates capsule chain length and interacts with Wzc, a part of the capsule polymerization and export machinery shared by many pathogens. We further show that RmpD confers HMV and regulates capsule chain length in a heterologous host (E. coli). As Wzc is a conserved protein found in many pathogens, it is possible that RmpD-mediated HMV and increased virulence may not be restricted to K. pneumoniae.

Deciphering Phage-Host Specificity Based on the Association of Phage Depolymerases and Bacterial Surface Glycan with Deep Learning

Phage tailspike proteins are depolymerases that target diverse bacterial surface glycans with high specificity, determining the host-specificity of numerous phages. To address the challenge of identifying tailspike proteins due to their sequence diversity, we developed SpikeHunter, an approach based on the ESM-2 protein language model. Using SpikeHunter, we successfully identified 231,965 tailspike proteins from a dataset comprising 8,434,494 prophages found within 165,365 genomes of five common pathogens. Among these proteins, 143,035 tailspike proteins displayed strong associations with serotypes. Moreover, we observed highly similar tailspike proteins in species that share closely related serotypes. We found extensive domain swapping in all five species, with the C-terminal domain being significantly associated with host serotype highlighting its role in host range determination. Our study presents a comprehensive cross-species analysis of tailspike protein to serotype associations, providing insights applicable to phage therapy and biotechnology.

Evolutionary and functional history of the Escherichia coli K1 capsule

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

A point mutation in recC associated with subclonal replacement of carbapenem-resistant Klebsiella pneumoniae ST11 in China

Adaptation to selective pressures is crucial for clinically important pathogens to establish epidemics, but the underlying evolutionary drivers remain poorly understood. The current epidemic of carbapenem-resistant Klebsiella pneumoniae (CRKP) poses a significant threat to public health. In this study we analyzed the genome sequences of 794 CRKP bloodstream isolates collected in 40 hospitals in China between 2014 and 2019. We uncovered a subclonal replacement in the predominant clone ST11, where the previously prevalent subclone OL101:KL47 was replaced by O2v1:KL64 over time in a stepwise manner. O2v1:KL64 carried a higher load of mobile genetic elements, and a point mutation exclusively detected in the recC of O2v1:KL64 significantly promotes recombination proficiency. The epidemic success of O2v1:KL64 was further associated with a hypervirulent sublineage with enhanced resistance to phagocytosis, sulfamethoxazole-trimethoprim, and tetracycline. The phenotypic alterations were linked to the overrepresentation of hypervirulence determinants and antibiotic genes conferred by the acquisition of an rmpA-positive pLVPK-like virulence plasmid and an IncFII-type multidrug-resistant plasmid, respectively. The dissemination of the sublineage was further promoted by more frequent inter-hospital transmission. The results collectively demonstrate that the expansion of O2v1:KL64 is correlated to a repertoire of genomic alterations convergent in a subpopulation with evolutionary advantages.

Capsular Polysaccharide Is Essential for the Virulence of the Antimicrobial-Resistant Pathogen Enterobacter hormaechei

Nosocomial infections caused by multidrug-resistant (MDR) Enterobacter cloacae complex (ECC) pathogens are on the rise. However, the virulence strategies employed by these pathogens remain elusive. Here, we study the interaction of ECC clinical isolates with human serum to define how this pathogen evades the antimicrobial action of complement, one of the first lines of host-mediated immune defense. We identified a small number of serum-sensitive strains, including Enterobacter hormaechei strain NR3055, which we exploited for the in vitro selection of serum-resistant clones. Comparative genomics between the serum-sensitive NR3055 strain and the isolated serum-resistant clones revealed a premature stop codon in the wzy gene of the capsular polysaccharide biosynthesis locus of NR3055. The complementation of wzy conferred serum resistance to NR3055, prevented the deposition of complement proteins on the bacterial surface, inhibited phagocytosis by human neutrophils, and rendered the bacteria virulent in a mouse model of peritonitis. Mice exposed to a nonlethal dose of encapsulated NR3055 were protected from subsequent lethal infections by encapsulated NR3055, whereas mice that were previously exposed to unencapsulated NR3055 succumbed to infection. Thus, capsule is a key immune evasion determinant for E. hormaechei, and it is a potential target for prophylactics and therapeutics to combat these increasingly MDR human pathogens. IMPORTANCE Infections caused by antimicrobial resistant bacteria are of increasing concern, especially those due to carbapenem-resistant Enterobacteriaceae pathogens. Included in this group are species of the Enterobacter cloacae complex, regarding which there is a paucity of knowledge on the infection biology of the pathogens, despite their clinical relevance. In this study, we combine techniques in comparative genomics, bacterial genetics, and diverse models of infection to establish capsule as an important mechanism of Enterobacter pathogens to resist the antibacterial activity of serum, a first line of host defense against bacterial infections. We also show that immune memory targeting the Enterobacter capsule protects against lethal infection. The further characterization of Enterobacter infection biology and the immune response to infection are needed for the development of therapies and preventative interventions targeting these highly antibiotic resistant pathogens.

Sulfated polysaccharides as multi target molecules to fight COVID 19 and comorbidities

The majority of research to combat SARS-CoV-2 infection exploits the adaptive immune system, but innate immunity, the first line of defense against pathogenic microbes, is equally important in understanding and controlling infectious diseases. Various cellular mechanisms provide physiochemical barriers to microbe infection in mucosal membranes and epithelia, with extracellular polysaccharides, particularly sulfated polysaccharides, being among the most widespread and potent extracellular and secreted molecules blocking and deactivating bacteria, fungi, and viruses. New research reveals that a range of polysaccharides effectively inhibits COV-2 infection of mammalian cells in culture. This review provides an overview of sulfated polysaccharides nomenclature, its significance as immunomodulators, antioxidants, antitumors, anticoagulants, antibacterial, and as potent antivirals. It summarizes current research on various interactions of sulfated polysaccharide with a range of viruses, including SARS-CoV-2, and their application for potential treatments for COVID-19. These molecules interact with biochemical signaling in immune cell responses, by actions in oxidative reactions, cytokine signaling, receptor binding, and through antiviral and antibacterial toxicity. These properties provide the potential for the development of novel therapeutic treatments for SARS-CoV-2 and other infectious diseases from modified polysaccharides.

A dedicated C-6 β-hydroxyacyltransferase required for biosynthesis of the glycolipid anchor for Vi antigen capsule in typhoidal Salmonella

Vi antigen is an extracellular polysaccharide produced by Salmonella enterica Typhi, Citrobacter freundii, and some soil bacteria belonging to the Burkholderiales. In Salmonella Typhi, Vi-antigen capsule protects the bacterium against host defenses, and the glycan is used in a current glycoconjugate vaccine to protect against typhoid. Vi antigen is a glycolipid assembled in the cytoplasm and translocated to the cell surface by an export complex driven by an ABC transporter. In Salmonella Typhi, efficient export and cell-surface retention of the capsule layer depend on a reducing terminal acylated-HexNAc moiety. Although the precise structure and biosynthesis of the acylated terminus has not been resolved, it distinguishes Vi antigen from other known glycolipid substrates for bacterial ABC transporters. The genetic locus for Vi antigen-biosynthesis encodes a single acyltransferase candidate (VexE), which is implicated in the acylation process. Here, we determined the structure of the VexE in vitro reaction product by mass spectrometry and nuclear magnetic resonance spectroscopy, to reveal that VexE catalyzes β-hydroxyacyl-ACP dependent acylation of the activated sugar precursor, uridine-5'-diphospho-N-acetylglucosamine (UDP-GlcNAc), at C-6 to form UDP-6-O-[β-hydroxymyristoyl]-α-d-GlcNAc. VexE belongs to the lysophosphatidyl acyltransferase (LPLAT) family, and comparison of an Alphafold VexE model to solved LPLAT structures, together with modeling enzyme:substrate complexes, led us to predict an enzyme mechanism. This study provides new insight into Vi terminal structure, offers a new model substrate to investigate the mechanism of glycolipid ABC transporters, and adds biochemical understanding for a novel reaction used in synthesis of an important bacterial virulence factor.

Identification of the Shigella flexneri Wzy Domain Modulating WzzpHS-2 Interaction and Detection of the Wzy/Wzz/Oag Complex

Shigella flexneri implements the Wzy-dependent pathway to biosynthesize the O antigen (Oag) component of its surface lipopolysaccharide. The inner membrane polymerase Wzy_SF catalyzes the repeat addition of undecaprenol-diphosphate-linked Oag (Und-PP-RUs) to produce a polysaccharide, the length of which is tightly regulated by two competing copolymerase proteins, Wzz_SF (short-type Oag; 10 to 17 RUs) and Wzz_pHS-2 (very-long-type Oag; >90 RUs). The nature of the interaction between Wzy_SF and Wzz_SF/Wzz_pHS-2 in Oag polymerization remains poorly characterized, with the majority of the literature characterizing the individual protein constituents of the Wzy-dependent pathway. Here, we report instead a major investigation into the specific binding interactions of Wzy_SF with its copolymerase counterparts. For the first time, a region of Wzy_SF that forms a unique binding site for Wzz_pHS-2 has been identified. Specifically, this work has elucidated key Wzy_SF moieties at the N- and C-terminal domains (NTD and CTD) that form an intramolecular pocket modulating the Wzz_pHS-2 interaction. Novel copurification data highlight that disruption of residues within this NTD-CTD pocket impairs the interaction with Wzz_pHS-2 without affecting Wzz_SF binding, thereby specifically disrupting polymerization of longer polysaccharide chains. This study provides a novel understanding of the molecular interaction of Wzy_SF with Wzz_SF/Wzz_pHS-2 in the Wzy-dependent pathway and, furthermore, detects the Wzy/Wzz/Und-PP-Oag complex for the first time. Beyond S. flexneri, this work may be extended to provide insight into the interactions between protein homologues expressed by related species, especially members of Enterobacteriaceae, that produce dual Oag chain length determinants. IMPORTANCE Shigella flexneri is a pathogen causing significant morbidity and mortality, predominantly devastating the pediatric age group in developing countries. A major virulence factor contributing to S. flexneri pathogenesis is its surface lipopolysaccharide, which is comprised of three domains: lipid A, core oligosaccharide, and O antigen (Oag). The Wzy-dependent pathway is the most common biosynthetic mechanism implemented for Oag biosynthesis by Gram-negative bacteria, including S. flexneri. The nature of the interaction between the polymerase, Wzy_SF, and the polysaccharide copolymerases, Wzz_SF and Wzz_pHS-2, in Oag polymerization is poorly characterized. This study investigates the molecular interplay between Wzy_SF and its copolymerases, deciphering key interactions in the Wzy-dependent pathway that may be extended beyond S. flexneri, providing insight into Oag biosynthesis in Gram-negative bacteria.

Deep Population Genomics Reveals Systematic and Parallel Evolution at a Lipopolysaccharide Biosynthetic Locus in Xanthomonas Pathogens That Infect Rice and Sugarcane

The advent of high-throughput sequencing and population genomics has enabled researchers to investigate selection pressure at hypervariable genomic loci encoding pathogen-associated molecular pattern (PAMP) molecules like lipopolysaccharide (LPS). Xanthomonas is a model and a major group of phytopathogenic bacteria that infect hosts in tissue-specific manner. Our in-depth population-based genomic investigation revealed the emergence of major lineages in two Xanthomonas pathogens that infect xylem of rice and sugarcane is associated with the acquisition and later large-scale replacement by distinct type of LPS cassettes. In the population of the rice xylem pathogen, Xanthomonas oryzae pv. oryzae (Xoo) and sugarcane pathogens Xanthomonas sacchari (Xsac) and Xanthomonas vasicola (Xvv), the BXO8 type of LPS cassette is replaced by a BXO1 type of cassette in Xoo and by Xvv type LPS cassette in Xsac and Xvv. These findings suggest a wave of parallel evolution at an LPS locus mediated by horizontal gene transfer (HGT) events during its adaptation and emergence. Aside from xylem pathogens, two closely related lineages of Xoo that infect parenchyma of rice and Leersia hexandra grass have acquired an LPS cassette from Xanthomonas pathogens that infect parenchyma of citrus, walnut, and strawberries, indicating yet another instance of parallel evolution mediated by HGT at an LPS locus. Our targeted and megapopulation-based genome dynamic studies revealed the acquisition and dominance of specific types of LPS cassettes in adaptation and success of a major group of phytopathogenic bacteria. IMPORTANCE Lipopolysaccharide (LPS) is a major microbe associated molecular pattern and hence a major immunomodulator. As a major and outer member component, it is expected that LPS is a frontline defense mechanism to deal with different host responses. Limited studies have indicated that LPS loci are also highly variable at strain and species level in plant-pathogenic bacteria, suggesting strong selection pressure from plants and associated niches. The advent of high-throughput genomics has led to the availability of a large set of genomic resources at taxonomic and population levels. This provides an exciting and important opportunity to carryout megascale targeted and population-based comparative genomic/association studies at important loci like those encoding LPS biosynthesis to understand their role in the evolution of the host, tissue specificity, and also predominant lineages. Such studies will also fill major gap in understanding host and tissue specificity in pathogenic bacteria. Our pioneering study uses the Xanthomonas group of phytopathogens that are known for their characteristic host and tissue specificity. The present deep phylogenomics of diverse Xanthomonas species and its members revealed lineage association and dominance of distinct types of LPS in accordance with their origin, host, tissue specificity, and evolutionary success.

Adaptation to novel spatially-structured environments is driven by the capsule and alters virulence-associated traits

The extracellular capsule is a major virulence factor, but its ubiquity in free-living bacteria with large environmental breadths suggests that it shapes adaptation to novel niches. Yet, how it does so, remains unexplored. Here, we evolve three Klebsiella strains and their capsule mutants in parallel. Their comparison reveals different phenotypic and genotypic evolutionary changes that alter virulence-associated traits. Non-capsulated populations accumulate mutations that reduce exopolysaccharide production and increase biofilm formation and yield, whereas most capsulated populations become hypermucoviscous, a signature of hypervirulence. Hence, adaptation to novel environments primarily occurs by fine-tuning expression of the capsular locus. The same evolutionary conditions selecting for mutations in the capsular gene wzc leading to hypermucoviscosity also result in increased susceptibility to antibiotics by mutations in the ramA regulon. This implies that general adaptive processes outside the host can affect capsule evolution and its role in virulence and infection outcomes may be a by-product of such adaptation.

Phenotypic and genotypic evolution in worrisome Klebsiella spp. is influenced by the capsule. Here the authors show that adaptation outside the host can impact virulence-associated traits, including de novo emergence of hypermucoviscosity.

Resistance of Dickeya solani strain IPO 2222 to lytic bacteriophage ΦD5 results in fitness tradeoffs for the bacterium during infection

Resistance to bacteriophage infections protects bacteria in phage-replete environments, enabling them to survive and multiply in the presence of their viral predators. However, such resistance may confer costs for strains, reducing their ecological fitness as expressed as competitiveness for resources or virulence or both. There is limited knowledge about such costs paid by phage-resistant plant pathogenic bacteria in their natural habitats. This study analyzed the costs of phage resistance paid by the phytopathogenic pectinolytic bacterium Dickeya solani both in vitro and in potato (Solanum tuberosum L.) plants. Thirteen Tn5 mutants of D. solani IPO 2222 were identified that exhibited resistance to infection by lytic bacteriophage vB_Dsol_D5 (ΦD5). The genes disrupted in these mutants encoded proteins involved in the synthesis of bacterial envelope components (viz. LPS, EPS and capsule). Although phage resistance did not affect most of the phenotypes of ΦD5-resistant D. solani such as growth rate, production of effectors, swimming and swarming motility, use of various carbon and nitrogen sources and biofilm formation evaluated in vitro, all phage resistant mutants were significantly compromised in their ability to survive on leaf surfaces as well as to grow within and cause disease symptoms in potato plants.

Identification of distinct capsule types associated with Serratia marcescens infection isolates

Serratia marcescens is a versatile opportunistic pathogen that can cause a variety of infections, including bacteremia. Our previous work established that the capsule polysaccharide (CPS) biosynthesis and translocation locus contributes to the survival of S. marcescens in a murine model of bacteremia and in human serum. In this study, we determined the degree of capsule genetic diversity among S. marcescens isolates. Capsule loci (KL) were extracted from >300 S. marcescens genome sequences and compared. A phylogenetic comparison of KL sequences demonstrated a substantial level of KL diversity within S. marcescens as a species and a strong delineation between KL sequences originating from infection isolates versus environmental isolates. Strains from five of the identified KL types were selected for further study and electrophoretic analysis of purified CPS indicated the production of distinct glycans. Polysaccharide composition analysis confirmed this observation and identified the constituent monosaccharides for each strain. Two predominant infection-associated clades, designated KL1 and KL2, emerged from the capsule phylogeny. Bacteremia strains from KL1 and KL2 were determined to produce ketodeoxynonulonic acid and N-acetylneuraminic acid, two sialic acids that were not found in strains from other clades. Further investigation of KL1 and KL2 sequences identified two genes, designated neuA and neuB, that were hypothesized to encode sialic acid biosynthesis functions. Disruption of neuB in a KL1 isolate resulted in the loss of sialic acid and CPS production. The absence of sialic acid and CPS production also led to increased susceptibility to internalization by a human monocytic cell line, demonstrating that S. marcescens phagocytosis resistance requires CPS. Together, these results establish the capsule genetic repertoire of S. marcescens and identify infection-associated clades with sialic acid CPS components.

Kaptive 2.0: updated capsule and lipopolysaccharide locus typing for the Klebsiella pneumoniae species complex

The outer polysaccharide capsule and lipopolysaccharide (LPS) antigens are key targets for novel control strategies targeting Klebsiella pneumoniae and related taxa from the K. pneumoniae species complex (KpSC), including vaccines, phage and monoclonal antibody therapies. Given the importance and growing interest in these highly diverse surface antigens, we had previously developed Kaptive, a tool for rapidly identifying and typing capsule (K) and outer LPS (O) loci from whole genome sequence data. Here, we report two significant updates, now freely available in Kaptive 2.0 (https://github.com/katholt/kaptive): (i) the addition of 16 novel K locus sequences to the K locus reference database following an extensive search of >17 000 KpSC genomes; and (ii) enhanced O locus typing to enable prediction of the clinically relevant O2 antigen (sub)types, for which the genetic determinants have been recently described. We applied Kaptive 2.0 to a curated dataset of >12 000 public KpSC genomes to explore for the first time, to the best of our knowledge, the distribution of predicted O (sub)types across species, sampling niches and clones, which highlighted key differences in the distributions that warrant further investigation. As the uptake of genomic surveillance approaches continues to expand globally, the application of Kaptive 2.0 will generate novel insights essential for the design of effective KpSC control strategies.

A Polysaccharide Biosynthesis Locus in Vibrio parahaemolyticus Important for Biofilm Formation Has Homologs Widely Distributed in Aquatic Bacteria Mainly from Gammaproteobacteria

Vibrio parahaemolyticus is a seafood-borne pathogen that poses a great threat to public health worldwide. It is found in either a planktonic cell or a biofilm form in the natural environment. The cps locus has been the only extensively studied polysaccharide biosynthesis gene cluster involved in biofilm formation for this bacterium. In this study, we found that an additional polysaccharide biosynthesis locus, scv, is also necessary for biofilm maturation. The scv locus is composed of two operons, and a loss of their expression leads to a defective biofilm phenotype. The transcription of the scv locus is under the control of a sigma 54-dependent response regulator, ScvE. In contrast, the quorum-sensing regulator AphA stimulates the expression of the cps locus and the scvABCD operon found in the scv locus. Bioinformatic analyses demonstrated that scv loci are divergent and widely distributed among 28 genera, including 26 belonging to the Gammaproteobacteria and 2 within the Alphaproteobacteria. We also determined that all scv locus-positive species are water-dwelling. Some strains of Aeromonas, Aliivibrio salmonicida, Pseudomonas anguilliseptica, Vibrio breoganii, and Vibrio scophthalmi probably acquired scv loci through insertion sequences and/or integrase-mediated horizontal gene transfer. Gene duplication and fusion were also detected in some scv homologs. Together, our results suggest that the genome of V. parahaemolyticus harbors two distinct polysaccharide biosynthesis loci, which may play a role in fine-tuning biofilm development, and that scv loci likely evolved by horizontal gene transfer, gene loss, gene duplication, and fragment fusion. IMPORTANCE Polysaccharides are the major component of biofilms, which provide survival advantages for bacteria in aquatic environments. The seafood-borne pathogen V. parahaemolyticus possesses a functionally uncharacterized polysaccharide biosynthesis locus, scv. We demonstrated that the scv locus is important for biofilm maturation and that scv expression is positively regulated by ScvE. Strains from 148 aquatic bacterial species possess scv homolog loci. These bacterial species belong to 28 genera, most of which belong to the Gammaproteobacteria class. The evolution and diversification of scv loci are likely driven by horizontal gene transfer, gene loss, gene duplication, and fragment fusion. Our results provide new insights into the function and evolution of this widespread polysaccharide biosynthesis locus.

Capsules and Extracellular Polysaccharides in Escherichia coli and Salmonella

Escherichia coli and Salmonella isolates produce a range of different polysaccharide structures that play important roles in their biology. E. coli isolates often possess capsular polysaccharides (K antigens), which form a surface structural layer. These possess a wide range of repeat-unit structures. In contrast, only one capsular polymer (Vi antigen) is found in Salmonella, and it is confined to typhoidal serovars. In both genera, capsules are vital virulence determinants and are associated with the avoidance of host immune defenses. Some isolates of these species also produce a largely secreted exopolysaccharide called colanic acid as part of their complex Rcs-regulated phenotypes, but the precise function of this polysaccharide in microbial cell biology is not fully understood. E. coli isolates produce two additional secreted polysaccharides, bacterial cellulose and poly-N-acetylglucosamine, which play important roles in biofilm formation. Cellulose is also produced by Salmonella isolates, but the genes for poly-N-acetylglucosamine synthesis appear to have been lost during its evolution toward enhanced virulence. Here, we discuss the structures, functions, relationships, and sophisticated assembly mechanisms for these important biopolymers.

Enhanced Recovery of Microbial Genes and Genomes From a Marine Water Column Using Long-Read Metagenomics

Third-generation sequencing has penetrated little in metagenomics due to the high error rate and dependence for assembly on short-read designed bioinformatics. However, second-generation sequencing metagenomics (mostly Illumina) suffers from limitations, particularly in the assembly of microbes with high microdiversity and retrieval of the flexible (adaptive) fraction of prokaryotic genomes. Here, we have used a third-generation technique to study the metagenome of a well-known marine sample from the mixed epipelagic water column of the winter Mediterranean. We have compared PacBio Sequel II with the classical approach using Illumina Nextseq short reads followed by assembly to study the metagenome. Long reads allow for efficient direct retrieval of complete genes avoiding the bias of the assembly step. Besides, the application of long reads on metagenomic assembly allows for the reconstruction of much more complete metagenome-assembled genomes (MAGs), particularly from microbes with high microdiversity such as Pelagibacterales. The flexible genome of reconstructed MAGs was much more complete containing many adaptive genes (some with biotechnological potential). PacBio Sequel II CCS appears particularly suitable for cellular metagenomics due to its low error rate. For most applications of metagenomics, from community structure analysis to ecosystem functioning, long reads should be applied whenever possible. Specifically, for in silico screening of biotechnologically useful genes, or population genomics, long-read metagenomics appears presently as a very fruitful approach and can be analyzed from raw reads before a computationally demanding (and potentially artifactual) assembly step.

Ecology and Function of the Transmissible Locus of Stress Tolerance in Escherichia coli and Plant-Associated Enterobacteriaceae

The transmissible locus of stress tolerance (tLST) is a genomic island which confers resistance to heat and chlorine. In this study, we determined that the tLST is frequent in genomes of those Enterobacteriaceae that occur in association with plants as well as the intestines of humans and animals and are relevant as nosocomial pathogens, e.g., Klebsiella and Cronobacter species. The tLST is more frequent in environmental and clinical isolates of Klebsiella pneumoniae than in animal isolates, and heat and chlorine resistance of tLST-positive strains of K. pneumoniae matched the resistance of tLST-positive strains of Escherichia coli. The function of 13 tLST genes was determined by assessing the heat and chlorine resistance of E. coli MG1655 mutants. The deletion of sHsp20, clpK_GI, sHsp_GI, pscA, pscB, and hdeD_GI reduced both heat and chlorine resistance; deletion of kefB reduced only chlorine resistance. Genes coding for heat shock proteins sHsp20, clpK_GI, and sHsp_GI decreased the oxidation of cytoplasmic proteins, while kefB decreased the oxidation of membrane lipids. The fitness cost of the tLST for E. coli MG1655 was assessed by pairwise competition experiments with isogenic tLST-positive or tLST-negative strains. The tLST imposes a fitness cost that is compensated for by frequent and lethal challenges with chlorine. All core genes need to be present to maintain the ecological advantage relative to the fitness cost. Taken together, core tLST genes are necessary to provide protection for E. coli against heat and chlorine stress, and the selective pressure for the tLST maintains core genes.

IMPORTANCE The transmissible locus of stress tolerance (tLST) is a genomic island comprising 10 core genes that occurs in diverse Enterobacteriaceae and confers resistance to heat and chlorine. Experimentation described in the manuscript describes the physiological function of the core genes by characterization of the resistance of 13 single-knockout (KO) mutants and by characterization of protein and membrane oxidation in these strains after chlorine challenge. Results identify tLST resistance as a genomic island that is specific for those Enterobacteriaceae that occur in plant-associated habitats as well in the intestines of vertebrates. In addition, the ecological function of the genomic island was characterized by large-scale genomic analysis and competition experiments of wild-type and mutant strains. Results suggest that tLST-mediated resistance to chlorine may contribute to the persistence of nosocomial pathogens in hospitals.

Interplay between the cell envelope and mobile genetic elements shapes gene flow in populations of the nosocomial pathogen Klebsiella pneumoniae

Mobile genetic elements (MGEs) drive genetic transfers between bacteria using mechanisms that require a physical interaction with the cellular envelope. In the high-priority multidrug-resistant nosocomial pathogens (ESKAPE), the first point of contact between the cell and virions or conjugative pili is the capsule. While the capsule can be a barrier to MGEs, it also evolves rapidly by horizontal gene transfer (HGT). Here, we aim at understanding this apparent contradiction by studying the covariation between the repertoire of capsule genes and MGEs in approximately 4,000 genomes of Klebsiella pneumoniae (Kpn). We show that capsules drive phage-mediated gene flow between closely related serotypes. Such serotype-specific phage predation also explains the frequent inactivation of capsule genes, observed in more than 3% of the genomes. Inactivation is strongly epistatic, recapitulating the capsule biosynthetic pathway. We show that conjugative plasmids are acquired at higher rates in natural isolates lacking a functional capsular locus and confirmed experimentally this result in capsule mutants. This suggests that capsule inactivation by phage pressure facilitates its subsequent reacquisition by conjugation. Accordingly, capsule reacquisition leaves long recombination tracts around the capsular locus. The loss and regain process rewires gene flow toward other lineages whenever it leads to serotype swaps. Such changes happen preferentially between chemically related serotypes, hinting that the fitness of serotype-swapped strains depends on the host genetic background. These results enlighten the bases of trade-offs between the evolution of virulence and multidrug resistance and caution that some alternatives to antibiotics by selecting for capsule inactivation may facilitate the acquisition of antibiotic resistance genes (ARGs).

A study of how the complex interaction between capsules and mobile genetic elements shapes gene flow in populations of Klebsiella pneumoniae reveals that capsule inactivation by phage pressure facilitates its subsequent re-acquisition by conjugation, and this loss and re-gain process influences the gene flow towards other lineages whenever it leads to serotype changes.

Nutrient conditions are primary drivers of bacterial capsule maintenance in Klebsiella

The fitness cost associated with the production of bacterial capsules is considered to be offset by the protection provided by these extracellular structures against biotic aggressions or abiotic stress. However, it is unknown if the capsule contributes to fitness in the absence of these. Here, we explored conditions favouring the maintenance of the capsule in Klebsiella pneumoniae, where the capsule is known to be a major virulence factor. Using short-term experimental evolution on different Klebsiella strains, we showed that small environmental variations have a strong impact on the maintenance of the capsule. Capsule inactivation is frequent in nutrient-rich, but scarce in nutrient-poor media. Competitions between wild-type and capsule mutants in nine different strains confirmed that the capsule is costly in nutrient-rich media. Surprisingly, these results also showed that the presence of a capsule provides a clear fitness advantage in nutrient-poor conditions by increasing both growth rates and population yields. The comparative analyses of the wild-type and capsule mutants reveal complex interactions between the environment, genetic background and serotype even in relation to traits known to be relevant during pathogenesis. In conclusion, our data suggest there are novel roles for bacterial capsules yet to be discovered and further supports the notion that the capsule's role in virulence may be a by-product of its contribution to bacterial adaptation outside the host.

Genomic diversity and organization of complex polysaccharide biosynthesis clusters in the genus Dickeya

The pectinolytic genus Dickeya (formerly Erwinia chrysanthemi) comprises numerous pathogenic species which cause diseases in various crops and ornamental plants across the globe. Their pathogenicity is governed by complex multi-factorial processes of adaptive virulence gene regulation. Extracellular polysaccharides and lipopolysaccharides present on bacterial envelope surface play a significant role in the virulence of phytopathogenic bacteria. However, very little is known about the genomic location, diversity, and organization of the polysaccharide and lipopolysaccharide biosynthetic gene clusters in Dickeya. In the present study, we report the diversity and structural organization of the group 4 capsule (G4C)/O-antigen capsule, putative O-antigen lipopolysaccharide, enterobacterial common antigen, and core lipopolysaccharide biosynthesis clusters from 54 Dickeya strains. The presence of these clusters suggests that Dickeya has both capsule and lipopolysaccharide carrying O-antigen to their external surface. These gene clusters are key regulatory components in the composition and structure of the outer surface of Dickeya. The O-antigen capsule/group 4 capsule (G4C) coding region shows a variation in gene content and organization. Based on nucleotide sequence homology in these Dickeya strains, two distinct groups, G4C group I and G4C group II, exist. However, comparatively less variation is observed in the putative O-antigen lipopolysaccharide cluster in Dickeya spp. except for in Dickeya zeae. Also, enterobacterial common antigen and core lipopolysaccharide biosynthesis clusters are present mostly as conserved genomic regions. The variation in the O-antigen capsule and putative O-antigen lipopolysaccharide coding region in relation to their phylogeny suggests a role of multiple horizontal gene transfer (HGT) events. These multiple HGT processes might have been manifested into the current heterogeneity of O-antigen capsules and O-antigen lipopolysaccharides in Dickeya strains during its evolution.

VPsero: Rapid Serotyping of Vibrio parahaemolyticus Using Serogroup-Specific Genes Based on Whole-Genome Sequencing Data

Vibrio parahaemolyticus has emerged as a significant enteropathogen in human and marine habitats worldwide, notably in regions where aquaculture products constitute a major nutritional source. It is a growing cause of diseases including gastroenteritis, wound infections, and septicemia. Serotyping assays use commercially available antisera to identify V. parahaemolyticus strains, but this approach is limited by high costs, complicated procedures, cross-immunoreactivity, and often subjective interpretation. By leveraging high-throughput sequencing technologies, we developed an in silico method based on comparison of gene clusters for lipopolysaccharide (LPSgc) and capsular polysaccharide (CPSgc) by firstly using the unique-gene strategy. The algorithm, VPsero, which exploits serogroup-specific genes as markers, covers 43 K and all 12 O serogroups in serotyping assays. VPsero is capable of predicting serotypes from assembled draft genomes, outputting LPSgc/CPSgc sequences, and recognizing possible novel serogroups or populations. Our tool displays high specificity and sensitivity in prediction toward V. parahaemolyticus strains, with an average sensitivity in serogroup prediction of 0.910 for O and 0.961 for K serogroups and a corresponding average specificity of 0.990 for O and 0.998 for K serogroups.