Complete Genome Sequence of Prepandemic Vibrio parahaemolyticus BB22OP

Virulence phenotypes differ between toxigenic Vibrio parahaemolyticus isolated from western coasts of Europe

Vibrio parahaemolyticus is the leading bacterial cause of gastroenteritis associated with seafood consumption worldwide. Not all members of the species are thought to be pathogenic, thus identification of virulent organisms is essential to protect public health and the seafood industry. Correlations of human disease and known genetic markers (e.g. thermostable direct hemolysin (TDH), TDH-related hemolysin (TRH)) appear complex. Some isolates recovered from patients lack these factors, while their presence has become increasingly noted in isolates recovered from the environment. Here, we used whole-genome sequencing in combination with mammalian and insect models of infection to assess the pathogenic potential of V. parahaemolyticus isolated from European Atlantic shellfish production areas. We found environmental V. parahaemolyticus isolates harboured multiple virulence-associated genes, including TDH and/or TRH. However, carriage of these factors did not necessarily reflect virulence in the mammalian intestine, as an isolate containing TDH and the genes coding for a type 3 secretion system (T3SS) 2α virulence determinant, appeared avirulent. Moreover, environmental V. parahaemolyticus lacking TDH or TRH could be assigned to groups causing low and high levels of mortality in insect larvae, with experiments using defined bacterial mutants showing that a functional T3SS1 contributed to larval death. When taken together, our findings highlight the genetic diversity of V. parahaemolyticus isolates found in the environment, their potential to cause disease and the need for a more systematic evaluation of virulence in diverse V. parahaemolyticus to allow better genetic markers.

The effect of environmental calcium on gene expression, biofilm formation and virulence of Vibrio parahaemolyticus

Calcium (Ca2+) can regulate the swarming motility and virulence of Vibrio parahaemolyticus BB22. However, the effects of Ca2+ on the physiology of V. parahaemolyticus RIMD2210633, whose genomic composition is quite different with that of BB22, have not been investigated. In this study, the results of phenotypic assays showed that the biofilm formation, c-di-GMP production, swimming motility, zebrafish survival rate, cytoxicity against HeLa cells, and adherence activity to HeLa cells of V. parahaemolyticus RIMD2210633 were significantly enhanced by Ca2+. However, Ca2+ had no effect on the growth, swarming motility, capsular polysaccharide (CPS) phase variation and hemolytic activity. The RNA sequencing (RNA-seq) assay disclosed 459 significantly differentially expressed genes (DEGs) in response to Ca2+, including biofilm formation-associated genes and those encode virulence factors and putative regulators. DEGs involved in polar flagellum and T3SS1 were upregulated, whereas majority of those involved in regulatory functions and c-di-GMP metabolism were downregulated. The work helps us understand how Ca2+ affects the behavior and gene expression of V. parahaemolyticus RIMD2210633.

Biological Function of Prophage-Related Gene Cluster ΔVpaChn25_RS25055~ΔVpaChn25_0714 of Vibrio parahaemolyticus CHN25

Vibrio parahaemolyticus is the primary foodborne pathogen known to cause gastrointestinal infections in humans. Nevertheless, the molecular mechanisms of V. parahaemolyticus pathogenicity are not fully understood. Prophages carry virulence and antibiotic resistance genes commonly found in Vibrio populations, and they facilitate the spread of virulence and the emergence of pathogenic Vibrio strains. In this study, we characterized three such genes, VpaChn25_0713, VpaChn25_0714, and VpaChn25_RS25055, within the largest prophage gene cluster in V. parahaemolyticus CHN25. The deletion mutants ΔVpaChn25_RS25055, ΔVpaChn25_0713, ΔVpaChn25_0714, and ΔVpaChn25_RS25055-0713-0714 were derived with homologous recombination, and the complementary mutants ΔVpaChn25_0713-com, ΔVpaChn25_0714-com, ΔVpaChn25_RS25055-com, ΔVpaChn25_RS25055-0713-0714-com were also constructed. In the absence of the VpaChn25_RS25055, VpaChn25_0713, VpaChn25_0714, and VpaChn25_RS25055-0713-0714 genes, the mutants showed significant reductions in low-temperature survivability and biofilm formation (p < 0.001). The ΔVpaChn25_0713, ΔVpaChn25_RS25055, and ΔVpaChn25_RS25055-0713-0714 mutants were also significantly defective in swimming motility (p < 0.001). In the Caco-2 model, the above four mutants attenuated the cytotoxic effects of V. parahaemolyticus CHN25 on human intestinal epithelial cells (p < 0.01), especially the ΔVpaChn25_RS25055 and ΔVpaChn25_RS25055-0713-0714 mutants. Transcriptomic analysis showed that 15, 14, 8, and 11 metabolic pathways were changed in the ΔVpaChn25_RS25055, ΔVpaChn25_0713, ΔVpaChn25_0714, and ΔVpaChn25_RS25055-0713-0714 mutants, respectively. We labeled the VpaChn25_RS25055 gene with superfolder green fluorescent protein (sfGFP) and found it localized at both poles of the bacteria cell. In addition, we analyzed the evolutionary origins of the above genes. In summary, the prophage genes VpaChn25_0713, VpaChn25_0714, and VpaChn25_RS25055 enhance V. parahaemolyticus CHN25's survival in the environment and host. Our work improves the comprehension of the synergy between prophage-associated genes and the evolutionary process of V. parahaemolyticus.

Investigation of exopolysaccharide formation and its impact on anaerobic succinate production with Vibrio natriegens

Vibrio natriegens is an emerging host for biotechnology due to its high growth and substrate consumption rates. In industrial processes typically fed-batch processes are applied to obtain high space-time yields. In this study, we established an aerobic glucose-limited fed-batch fermentation with the wild type (wt) of V. natriegens which yielded biomass concentrations of up to 28.4 gX  L-1 . However, we observed that the viscosity of the culture broth increased by a factor of 800 at the end of the cultivation due to the formation of 157 ± 20 mg exopolysaccharides (EPS) L-1 . Analysis of the genomic repertoire revealed several genes and gene clusters associated with EPS formation. Deletion of the transcriptional regulator cpsR in V. natriegens wt did not reduce EPS formation, however, it resulted in a constantly low viscosity of the culture broth and altered the carbohydrate content of the EPS. A mutant lacking the cps cluster secreted two-fold less EPS compared to the wt accompanied by an overall low viscosity and a changed EPS composition. When we cultivated the succinate producer V. natriegens Δlldh Δdldh Δpfl Δald Δdns::pycCg (Succ1) under anaerobic conditions on glucose, we also observed an increased viscosity at the end of the cultivation. Deletion of cpsR and the cps cluster in V. natriegens Succ1 reduced the viscosity five- to six-fold which remained at the same level observed at the start of the cultivation. V. natriegens Succ1 ΔcpsR and V. natriegens Succ1 Δcps achieved final succinate concentrations of 51 and 46 g L-1 with a volumetric productivity of 8.5 and 7.7 gSuc L-1  h-1 , respectively. Both strains showed a product yield of about 1.4 molSuc molGlc -1 , which is 27% higher compared with that of V. natriegens Succ1 and corresponds to 81% of the theoretical maximum.

Insights into virulence: structure classification of the Vibrio parahaemolyticus RIMD mobilome

IMPORTANCE

The pandemic Vpar strain RIMD causes seafood-borne illness worldwide. Previous comparative genomic studies have revealed pathogenicity islands in RIMD that contribute to the success of the strain in infection. However, not all virulence determinants have been identified, and many of the proteins encoded in known pathogenicity islands are of unknown function. Based on the EOCD database, we used evolution-based classification of structure models for the RIMD proteome to improve our functional understanding of virulence determinants acquired by the pandemic strain. We further identify and classify previously unknown mobile protein domains as well as fast evolving residue positions in structure models that contribute to virulence and adaptation with respect to a pre-pandemic strain. Our work highlights key contributions of phage in mediating seafood born illness, suggesting this strain balances its avoidance of phage predators with its successful colonization of human hosts.

Vibrio parahaemolyticus Epidemiology and Pathogenesis: Novel Insights on an Emerging Foodborne Pathogen

The epidemiological dynamics of V. parahaemolyticus´ infections have been characterized by the abrupt appearance of outbreaks in remote areas where these diseases had not been previously detected, without knowing the routes of entry of the pathogens in the new area. However, there are recent studies that show the link between the appearance of epidemic outbreaks of Vibrio and environmental factors such as oceanic transport of warm waters, which has provided a possible mechanism for the dispersion of Vibrio diseases globally. Despite this evidence, there is little information on the possible routes of entry and transport of infectious agents from endemic countries to the entire world. In this sense, the recent advances in genomic sequencing tools are making it possible to infer possible biogeographical patterns of diverse pathogens with relevance in public health like V. parahaemolyticus. In this chapter, we will address several general aspects about V. parahaemolyticus, including their microbiological and genetic detection, main virulence factors, and the epidemiology of genotypes involved in foodborne outbreaks globally.

QsvR and OpaR coordinately repress biofilm formation by Vibrio parahaemolyticus

Mature biofilm formation by Vibrio parahaemolyticus requires exopolysaccharide (EPS), type IV pili, and capsular polysaccharide (CPS). Production of each is strictly regulated by various control pathways including quorum sensing (QS) and bis-(3'-5')-cyclic di-GMP (c-di-GMP). QsvR, an AraC-type regulator, integrates into the QS regulatory cascade via direct control of the transcription of the master QS regulators, AphA and OpaR. Deletion of qsvR in wild-type or opaR mutant backgrounds altered the biofilm formation by V. parahaemolyticus, suggesting that QsvR may coordinate with OpaR to control biofilm formation. Herein, we demonstrated both QsvR and OpaR repressed biofilm-associated phenotypes, c-di-GMP metabolism, and the formation of V. parahaemolyticus translucent (TR) colonies. QsvR restored the biofilm-associated phenotypic changes caused by opaR mutation, and vice versa. In addition, QsvR and OpaR worked coordinately to regulate the transcription of EPS-associated genes, type IV pili genes, CPS genes and c-di-GMP metabolism-related genes. These results demonstrated how QsvR works with the QS system to regulate biofilm formation by precisely controlling the transcription of multiple biofilm formation-associated genes in V. parahaemolyticus.

Vibrio parahaemolyticus T6SS2 effector repertoires

All strains of the marine bacterium Vibrio parahaemolyticus harbor a type VI secretion system (T6SS) named T6SS2, suggesting that this system plays an important role in the life cycle of this emerging pathogen. Although T6SS2 was recently shown to play a role in interbacterial competition, its effector repertoire remains unknown. Here, we employed proteomics to investigate the T6SS2 secretome of two V. parahaemolyticus strains, and we identified several antibacterial effectors encoded outside of the main T6SS2 gene cluster. We revealed two T6SS2-secreted proteins that are conserved in this species, indicating that they belong to the core secretome of T6SS2; other identified effectors are found only in subsets of strains, suggesting that they comprise an accessory effector arsenal of T6SS2. Remarkably, a conserved Rhs repeat-containing effector serves as a quality control checkpoint and is required for T6SS2 activity. Our results reveal effector repertoires of a conserved T6SS, including effectors that have no known activity and that have not been previously associated with T6SSs.

Genomic attributes differ between Vibrio parahaemolyticus environmental and clinical isolates including pathotypes

Vibrio parahaemolyticus is a marine bacterium and causes opportunistic gastroenteritis in humans. Clinical strains of V. parahaemolyticus contain haemolysin and type III secretion systems (T3SS) that define their pathotype. A growing number of strains isolated recently from the environment have acquired these virulence genes constituting a pool of potential pathogens. This study used comparative genomics to identify genetic factors that delineate environmental and clinical V. parahaemolyticus population and understand the similarities and differences between the T3SS2 phylotypes. The comparative analysis revealed the presence of a cluster of genes belonging to bacterial cellulose synthesis (bcs) in isolates of environmental origin. This cluster, previously unreported in V. parahaemolyticus, exhibit significant similarity to that of Aliivibrio fischeri, and might dictate a potentially new mechanism of its environmental adaptation and persistence. The study also identified many genes predicted in silico to be T3SS effectors that are unique to T3SS2β of tdh^- trh⁺ and tdh⁺ trh⁺ pathotype and having no identifiable homologue in tdh⁺ trh^- T3SS2α. Overall, these findings highlight the importance of understanding the genes and strategies V. parahaemolyticus utilize for the myriad interactions with its hosts, either marine invertebrates or humans.

A Rapid Fluorescence-Based Screen to Identify Regulators and Components of Interbacterial Competition Mechanisms in Bacteria

Contact-dependent antibacterial mechanisms enhance bacterial fitness as they enable bacteria to outcompete their rivals and thrive in diverse environments. Such systems also allow pathogenic bacteria to establish a niche inside a host, where they must compete with commensal microflora. In many cases, antibacterial systems are tightly regulated by complex sensor and signal transduction networks. Deciphering these regulatory networks, as well as identifying functional components of antibacterial mechanisms, are valuable objectives since essential regulators and components present possible targets for developing antivirulence therapies. Here we describe Bacterial Competition Fluorescence (BaCoF), a methodology that relies on a fluorescence signal to determine the outcome of bacterial competitions. This methodology enables screening of mutant libraries to identify genes that are essential for activating a contact-dependent antibacterial system of interest. Thus, this methodology can be applied to reveal essential regulators and components of antibacterial systems in bacterial pathogens.

OpaR Controls the Metabolism of c-di-GMP in Vibrio parahaemolyticus

Vibrio parahaemolyticus, the leading cause of seafood-associated gastroenteritis worldwide, has a strong ability to form biofilms on surfaces. Quorum sensing (QS) is a process widely used by bacteria to communicate with each other and control gene expression via the secretion and detection of autoinducers. OpaR is the master QS regulator of V. parahaemolyticus operating under high cell density (HCD). OpaR regulation of V. parahaemolyticus biofilm formation has been reported, but the regulatory mechanisms are still not fully understood. bis-(3'-5')-cyclic di-GMP (c-di-GMP) is an omnipresent intracellular second messenger that regulates diverse behaviors of bacteria including activation of biofilm formation. In this work, we showed that OpaR repressed biofilm formation and decreased the intracellular concentration of c-di-GMP in V. parahaemolyticus RIMD2210633. The OpaR box-like sequences were detected within the regulatory DNA regions of scrA, scrG, VP0117, VPA0198, VPA1176, VP0699, and VP2979, encoding a group of GGDEF and/or EAL-type proteins. The results of qPCR, LacZ fusion, EMSA, and DNase I footprinting assays demonstrated that OpaR bound to the upstream DNA regions of scrA, VP0117, VPA0198, VPA1176, and VP0699 to repress their transcription, whereas it positively and directly regulated the transcription of scrG and VP2979. Thus, transcriptional regulation of these genes by OpaR led directly to changes in the intracellular concentration of c-di-GMP. The direct association between QS and c-di-GMP metabolism in V. parahaemolyticus RIMD2210633 would be conducive to precise control of gene transcription and bacterial behaviors such as biofilm formation.

Genotypic Diversity and Pathogenic Potential of Clinical and Environmental Vibrio parahaemolyticus Isolates From Brazil

Vibrio parahaemolyticus strains recovered from human diarrheal stools (one in 1975 and two in 2001) and environmental sources (four, between 2008 and 2010) were investigated for the presence of virulence genes (trh, tdh, and vpadF), pandemic markers (orf8, toxRS_new), and with respect to their pathogenic potential in two systemic infection models. Based only on the presence or absence of these genetic markers, they were classified as follows: the environmental strains were non-pathogenic, whereas among the clinical strains, the one isolated in 1975 was pathogenic (non-pandemic), and the other two were pathogenic (pandemic). The pathogenic potential of the strains was evaluated in mice and Galleria mellonella larvae infection models, and except for the clinical (pathogenic, non-pandemic) isolate, the others produced lethal infection in both organisms, regardless of their source, serotype, and genotype (tdh, orf8, toxRS_new, and vpadF). Based on mice and larval mortality rates, the strains were then grouped according to virulence (high, intermediate, and avirulent), and remarkably similar results were obtained by using these models: The clinical strain (pathogenic and non-pandemic) was classified as avirulent, and other strains (four non-pathogenic and two pandemic) were considered of high or intermediate virulence. In summary, these findings demonstrate that G. mellonella larvae can indeed be used as an alternative model to study the pathogenicity of V. parahaemolyticus. Moreover, they raise doubts about the use of traditional virulence markers to predict pathogenesis of the species and show that reliable models are indispensable to determine the pathogenic potential of environmental isolates considered non-pathogenic, based on the absence of the long-standing virulence indicators.

Prophage-Related Gene VpaChn25_0724 Contributes to Cell Membrane Integrity and Growth of Vibrio parahaemolyticus CHN25

Vibrio parahaemolyticus is a leading seafood-borne pathogen that can cause acute gastroenteritis and even death in humans. In aquatic ecosystems, phages constantly transform bacterial communities by horizontal gene transfer. Nevertheless, biological functions of prophage-related genes in V. parahaemolyticus remain to be fully unveiled. Herein, for the first time, we studied one such gene VpaChn25_0724 encoding an unknown hypothetical protein in V. parahaemolyticus CHN25. This gene deletion mutant ΔVpaChn25_0724 was constructed by homologous recombination, and its complementary mutant ΔVpaChn25_0724-com was also obtained. The ΔVpaChn25_0724 mutant exhibited a sever defect in growth and swimming motility particularly at lower temperatures. Biofilm formation and cytotoxicity capacity of V. parahaemolyticus CHN25 was significantly lowered in the absence of VpaChn25_0724. Comparative secretomic analysis revealed an increase in extracellular proteins of ΔVpaChn25_0724, which likely resulted from its damaged cell membrane. Comparison of transcriptome data showed twelve significantly altered metabolic pathways in ΔVpaChn25_0724, suggesting inactive transport and utilization of carbon sources, repressed energy production and membrane biogenesis in ΔVpaChn25_0724. Comparative transcriptomic analysis also revealed several remarkably down-regulated key regulators in bacterial gene regulatory networks linked to the observed phenotypic variations. Overall, the results here facilitate better understanding of biological significance of prophage-related genes remaining in V. parahaemolyticus.

A comparative genomics methodology reveals a widespread family of membrane-disrupting T6SS effectors

Gram-negative bacteria deliver effectors via the type VI secretion system (T6SS) to outcompete their rivals. Each bacterial strain carries a different arsenal of effectors; the identities of many remain unknown. Here, we present an approach to identify T6SS effectors encoded in bacterial genomes of interest, without prior knowledge of the effectors’ domain content or genetic neighborhood. Our pipeline comprises a comparative genomics analysis followed by screening using a surrogate T6SS⁺ strain. Using this approach, we identify an antibacterial effector belonging to the T6SS1 of Vibrio parahaemolyticus, representing a widespread family of T6SS effectors sharing a C-terminal domain that we name Tme (Type VI membrane-disrupting effector). Tme effectors function in the periplasm where they intoxicate bacteria by disrupting membrane integrity. We believe our approach can be scaled up to identify additional T6SS effectors in various bacterial genera.

Gram-negative bacteria deliver effectors via the type VI secretion system (T6SS) to outcompete their rivals. Here, Fridman et al. present an approach to identify T6SS effectors encoded in bacterial genomes without prior knowledge of their domain content or genetic neighbourhood, and identify a new family of membrane-disrupting effectors.

Regulatory effects of Shewanella putrefaciens isolated from shrimp Penaeus orientalis on the virulence factors of Vibrio parahaemolyticus and evaluation of the role of quorum sensing in virulence factors regulation

As an aquatic pathogen widely present in aquatic food, Vibrio parahaemolyticus causes outbreaks of gastroenteritis across the globe. Virulence factors of V. parahaemolyticus increases with the amount of spoilage in aquatic organisms including shrimp, but mechanisms regulating its virulence factors are not well understood. In this study, five spoilage bacteria isolated from shrimp were investigated for their regulatory effects on the virulence factors including haemolysin and biofilm of V. parahaemolyticus. Among these isolates, Shewanella putrefaciens induced haemolytic activity in V. parahaemolyticus in a time-dose-temperature-dependent manner and we found the main component responsible for this effect to be the supernatant or cell-free extract of S. putrefaciens. Total haemolytic activity, expression of the thermostable direct haemolysin gene tdh and biofilm production of V. parahaemolyticus were significantly up-regulated by S. putrefaciens, but also by deletion of quorum-sensing luxM or luxS gene of V. parahaemolyticus. However, this regulation by S. putrefaciens was significantly impaired by deletion of the luxM gene, but not by deletion of the luxS gene. Further study showed that S. putrefaciens exhibited a strong degradation ability on the signalling molecule acylated homoserine lactone (AHL) synthesised by the LuxM enzyme. This study revealed a novel virulence regulatory mechanism that S. putrefaciens can significantly increase the virulence factors of V. parahaemolyticus via interfering with the luxM- type quorum-sensing signalling pathway through its AHL-degradation ability.

Cell-cell communication, chemotaxis and recruitment in Vibrio parahaemolyticus

Motile bacteria are proficient at finding optimal environments for colonization. Often, they use chemotaxis to sense nutrient availability and dangerous concentrations of toxic chemicals. For many bacteria, the repertoire of chemoreceptors is large, suggesting they possess a broad palate with respect to sensing. However, knowledge of the molecules detected by chemotaxis signal transduction systems is limited. Some bacteria, like Vibrio parahaemolyticus, are social and swarm in groups on surfaces. This marine bacterium and human pathogen secretes the S signal autoinducer, which cues degradation of intracellular c-di-GMP leading to transcription of the swarming program. Here, we report that the S signal also directs motility at a behavioral level by serving as a chemoattractant. The data demonstrate that V. parahaemolyticus senses the S signal using SscL and SscS, homologous methyl-accepting chemotaxis proteins. SscL is required by planktonic bacteria for S signal chemotaxis. SscS plays a role during swarming, and mutants lacking this chemoreceptor swarm faster and produce colonies with more deeply branched swarming fronts than the wild type or the sscL mutant. Other Vibrio species can swim toward the S signal, suggesting a recruitment role for this cell-cell communication molecule in the context of polymicrobial marine communities.

Comparative Genomic Analysis Reveals the Potential Risk of Vibrio parahaemolyticus Isolated From Ready-To-Eat Foods in China

Vibrio parahaemolyticus is a major foodborne pathogen associated with the consumption of aquatic products. The presence of this bacterium in ready-to-eat (RTE) foods has recently been reported. However, the genomic features and potential risks of V. parahaemolyticus isolated from RTE foods are poorly understood. To help understand the genome-wide characteristics of RTE food isolates, the complete genomes of 27 RTE food isolates were sequenced and compared to those of 20 clinical and 19 other environmental (e.g., water and aquatic product source) isolates using a comparative genomics approach. Analysis revealed that V. parahaemolyticus RTE food isolates had higher numbers of genes on average and possessed more accessory genes than isolates from other sources. Most RTE food isolates were positive for some known virulence-associated genes and pathogenicity islands (PAIs), and some of these isolates were genetically homologous to clinical isolates. Genome-wide association analysis revealed 79 accessory genes and 78 missense single-nucleotide polymorphisms that affected 11 protein-coding genes were significantly associated with RTE food sources. These genes were mostly involved in defense mechanisms and energy production and conversion according to functional annotation in the COG database. KEGG Pathway analysis showed that these genes mainly affected the biofilm formation of V. parahaemolyticus, and subsequent experiments confirmed that nearly all RTE food isolates possessed the ability to form biofilm. The biofilm formation can facilitate the persistence of V. parahaemolyticus in RTE foods, and the presence of virulence-associated genes poses a pathogenic potential to humans. Our findings highlight the potential risk of V. parahaemolyticus in Chinese RTE foods and illustrate the genomic basis for the persistence of these isolates. This study will aid in re-evaluating the food safety threats conferred by this bacterium.

Regulation of Thermostable Direct Hemolysin and Biofilm Formation of Vibrio parahaemolyticus by Quorum-Sensing Genes luxM and luxS

Vibrio parahaemolyticus is a seafood opportunistic pathogen. There are evidences suggesting that virulence skills, including hemolytic activity and biofilm formation, are regulated by the luxM/luxS-dependent quorum-sensing system in V. parahaemolyticus, and their regulatory mechanism is not well understood. To better understand the virulence regulatory mechanism of V. parahaemolyticus, the luxM deletion (△luxM) and luxS deletion (△luxS) mutants were constructed and their impacts on growth, hemolysin activity, and biofilm were investigated. Results show that both luxM and luxS are involved in the adaptation to environmental conditions in early adaptive-log phase growth of V. parahaemolyticus. Thermostable direct hemolysin gene (tdh) was negatively regulated by luxM and positively regulated by luxS. The biofilm formation was negatively regulated by both luxS and luxM. This study provides an insight into some aspects of V. parahaemolyticus virulence regulation by luxM/luxS-dependent quorum-sensing system.

Plasmid-Mediated Quinolone Resistance (PMQR) Genes and Class 1 Integrons in Quinolone-Resistant Marine Bacteria and Clinical Isolates of Escherichia coli from an Aquacultural Area

Antimicrobial usage in aquaculture selects for antimicrobial-resistant microorganisms in the marine environment. The relevance of this selection to terrestrial animal and human health is unclear. Quinolone-resistance genes qnrA, qnrB, and qnrS were chromosomally located in four randomly chosen quinolone-resistant marine bacteria isolated from an aquacultural area with heavy quinolone usage. In quinolone-resistant uropathogenic clinical isolates of Escherichia coli from a coastal area bordering the same aquacultural region, qnrA was chromosomally located in two E. coli isolates, while qnrB and qnrS were located in small molecular weight plasmids in two other E. coli isolates. Three quinolone-resistant marine bacteria and three quinolone-resistant E. coli contained class 1 integrons but without physical association with PMQR genes. In both marine bacteria and uropathogenic E. coli, class 1 integrons had similar co-linear structures, identical gene cassettes, and similarities in their flanking regions. In a Marinobacter sp. marine isolate and in one E. coli clinical isolate, sequences immediately upstream of the qnrS gene were homologous to comparable sequences of numerous plasmid-located qnrS genes while downstream sequences were different. The observed commonality of quinolone resistance genes and integrons suggests that aquacultural use of antimicrobials might facilitate horizontal gene transfer between bacteria in diverse ecological locations.

Environmental vibrios represent a source of antagonistic compounds that inhibit pathogenic Vibrio cholerae and Vibrio parahaemolyticus strains

With the overuse of antibiotics, many pathogens including Vibrio cholerae and Vibrio parahaemolyticus have evolved multidrug resistance making treatment more difficult. While understanding the mechanisms that underlie pathogenesis is crucial, knowledge of bacterial interactions of V. cholerae and V. parahaemolyticus could provide insight to their susceptibility outside of the human host. Based on previous work showing competition among environmental strains, we predict that marine-derived bacteria should inhibit Vibrio pathogens and may be a source of unique antibiotic compounds. We tested a collection of 3,456 environmental Vibrio isolates from diverse habitats against a panel of V. cholerae and V. parahaemolyticus, and identified 102 strains that inhibited the growth of these pathogens. Phylogenetic analysis revealed that 40 pathogen-inhibiting strains were unique at the hsp60 gene sequence while 62 of the isolates were identical suggesting clonal groups. Genomic comparisons of ten strains revealed diversity even between clonal isolates and were identified as being closely related to known Vibrio crassostreae, Vibrio splendidus, and Vibrio tasmaniensis strains. Further analysis revealed multiple biosynthetic gene clusters within all sequenced genomes that encoded secondary metabolites with potential antagonistic activity. Thus, environmental vibrios represent a source of compounds that inhibit Vibrio pathogens.

Genomic and transcriptomic analyses reveal distinct biological functions for cold shock proteins (VpaCspA and VpaCspD) in Vibrio parahaemolyticus CHN25 during low-temperature survival

Background

Vibrio parahaemolyticus causes serious seafood-borne gastroenteritis and death in humans. Raw seafood is often subjected to post-harvest processing and low-temperature storage. To date, very little information is available regarding the biological functions of cold shock proteins (CSPs) in the low-temperature survival of the bacterium. In this study, we determined the complete genome sequence of V. parahaemolyticus CHN25 (serotype: O5:KUT). The two main CSP-encoding genes (VpacspA and VpacspD) were deleted from the bacterial genome, and comparative transcriptomic analysis between the mutant and wild-type strains was performed to dissect the possible molecular mechanisms that underlie low-temperature adaptation by V. parahaemolyticus.

Results

The 5,443,401-bp V. parahaemolyticus CHN25 genome (45.2% G + C) consisted of two circular chromosomes and three plasmids with 4,724 predicted protein-encoding genes. One dual-gene and two single-gene deletion mutants were generated for VpacspA and VpacspD by homologous recombination. The growth of the ΔVpacspA mutant was strongly inhibited at 10 °C, whereas the VpacspD gene deletion strongly stimulated bacterial growth at this low temperature compared with the wild-type strain. The complementary phenotypes were observed in the reverse mutants (ΔVpacspA-com, and ΔVpacspD-com). The transcriptome data revealed that 12.4% of the expressed genes in V. parahaemolyticus CHN25 were significantly altered in the ΔVpacspA mutant when it was grown at 10 °C. These included genes that were involved in amino acid degradation, secretion systems, sulphur metabolism and glycerophospholipid metabolism along with ATP-binding cassette transporters. However, a low temperature elicited significant expression changes for 10.0% of the genes in the ΔVpacspD mutant, including those involved in the phosphotransferase system and in the metabolism of nitrogen and amino acids. The major metabolic pathways that were altered by the dual-gene deletion mutant (ΔVpacspAD) radically differed from those that were altered by single-gene mutants. Comparison of the transcriptome profiles further revealed numerous differentially expressed genes that were shared among the three mutants and regulators that were specifically, coordinately or antagonistically modulated by VpaCspA and VpaCspD. Our data also revealed several possible molecular coping strategies for low-temperature adaptation by the bacterium.

Conclusions

This study is the first to describe the complete genome sequence of V. parahaemolyticus (serotype: O5:KUT). The gene deletions, complementary insertions, and comparative transcriptomics demonstrate that VpaCspA is a primary CSP in the bacterium, while VpaCspD functions as a growth inhibitor at 10 °C. These results have improved our understanding of the genetic basis for low-temperature survival by the most common seafood-borne pathogen worldwide.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3784-5) contains supplementary material, which is available to authorized users.

Defining a Core Genome Multilocus Sequence Typing Scheme for the Global Epidemiology of Vibrio parahaemolyticus

Vibrio parahaemolyticus is an important human foodborne pathogen whose transmission is associated with the consumption of contaminated seafood, with a growing number of infections reported over recent years worldwide. A multilocus sequence typing (MLST) database for V. parahaemolyticus was created in 2008, and a large number of clones have been identified, causing severe outbreaks worldwide (sequence type 3 [ST3]), recurrent outbreaks in certain regions (e.g., ST36), or spreading to other regions where they are nonendemic (e.g., ST88 or ST189). The current MLST scheme uses sequences of 7 genes to generate an ST, which results in a powerful tool for inferring the population structure of this pathogen, although with limited resolution, especially compared to pulsed-field gel electrophoresis (PFGE). The application of whole-genome sequencing (WGS) has become routine for trace back investigations, with core genome MLST (cgMLST) analysis as one of the most straightforward ways to explore complex genomic data in an epidemiological context. Therefore, there is a need to generate a new, portable, standardized, and more advanced system that provides higher resolution and discriminatory power among V. parahaemolyticus strains using WGS data. We sequenced 92 V. parahaemolyticus genomes and used the genome of strain RIMD 2210633 as a reference (with a total of 4,832 genes) to determine which genes were suitable for establishing a V. parahaemolyticus cgMLST scheme. This analysis resulted in the identification of 2,254 suitable core genes for use in the cgMLST scheme. To evaluate the performance of this scheme, we performed a cgMLST analysis of 92 newly sequenced genomes, plus an additional 142 strains with genomes available at NCBI. cgMLST analysis was able to distinguish related and unrelated strains, including those with the same ST, clearly showing its enhanced resolution over conventional MLST analysis. It also distinguished outbreak-related from non-outbreak-related strains within the same ST. The sequences obtained from this work were deposited and are available in the public database (http://pubmlst.org/vparahaemolyticus). The application of this cgMLST scheme to the characterization of V. parahaemolyticus strains provided by different laboratories from around the world will reveal the global picture of the epidemiology, spread, and evolution of this pathogen and will become a powerful tool for outbreak investigations, allowing for the unambiguous comparison of strains with global coverage.

Molecular and Genomic Characterization of Vibrio mimicus Isolated from a Frozen Shrimp Processing Facility in Mexico

Vibrio mimicus is a gram-negative bacterium responsible for diseases in humans. Three strains of V. mimicus identified as V. mimicus 87, V. mimicus 92 and V. mimicus 93 were isolated from a shrimp processing facility in Guaymas, Sonora, Mexico. The strains were analyzed using several molecular techniques and according to the cluster analysis they were different, their similarities ranged between 51.3% and 71.6%. ERIC-PCR and RAPD (vmh390R) were the most discriminatory molecular techniques for the differentiation of these strains. The complete genomes of two strains (V. mimicus 87, renamed as CAIM 1882, and V. mimicus 92, renamed as CAIM 1883) were sequenced. The sizes of the genomes were 3.9 Mb in both strains, with 2.8 Mb in ChI and 1.1 Mb in ChII. A 12.7% difference was found in the proteome content (BLAST matrix). Several virulence genes were detected (e.g. capsular polysaccharide, an accessory colonization factor and genes involved in quorum-sensing) which were classified in 16 categories. Variations in the gene content between these genomes were observed, mainly in proteins and virulence genes (e.g., hemagglutinin, mobile elements and membrane proteins). According to these results, both strains were different, even when they came from the same source, giving an insight of the diversity of V. mimicus. The identification of various virulence genes, including a not previously reported V. mimicus gene (acfD) in ChI in all sequenced strains, supports the pathogenic potential of this species. Further analysis will help to fully understand their potential virulence, environmental impact and evolution.

Multicolor Melting Curve Analysis-Based Multilocus Melt Typing of Vibrio parahaemolyticus

Vibrio parahaemolyticus is the leading cause of seafood-borne gastroenteritis outbreaks. To track the source of these diseases in a timely manner, a high throughput typing method is critical. We hereby describe a novel genotyping method for V. parahaemolyticus, termed multilocus melt typing (MLMT), based on multilocus sequence typing (MLST). MLMT utilizes melting curve analysis to interrogate the allelic types of a set of informative single nucleotide polymorphisms (SNPs) derived from the housekeeping genes used in MLST. For each SNP, one allelic type generates distinct T_m values, which are converted into a binary code. Multiple SNPs thus generate a series of binary codes, forming a melt type (MT) corresponding with a sequence type (ST) of MLST. Using a set of 12 SNPs, the MLMT scheme could resolve 218 V.parahaemolyticus isolates into 50 MTs corresponding with 56 STs. The discriminatory power of MLMT and MLST was similar with Simpson’s index of diversity of 0.638 and 0.646, respectively. The global (adjusted Rand index = 0.982) and directional congruence (adjusted Wallace coefficient, MT→ST = 0.965; ST→MT = 1.000) between the two typing approaches was high. The entire procedure of MLMT could be finished within 3 h with negligible hands on time in a real-time PCR machine. We conclude that MLMT provides a reliable and efficient approach for V. parahaemolyticus genotyping and might also find use in other pathogens.

OpaR controls a network of downstream transcription factors in Vibrio parahaemolyticus BB22OP

Vibrio parahaemolyticus is an emerging world-wide human pathogen that is associated with food-borne gastroenteritis when raw or undercooked seafood is consumed. Expression of virulence factors in this organism is modulated by the phenomenon known as quorum sensing, which permits differential gene regulation at low versus high cell density. The master regulator of quorum sensing in V. parahaemolyticus is OpaR. OpaR not only controls virulence factor gene expression, but also the colony and cellular morphology associated with growth on a surface and biofilm formation. Whole transcriptome Next Generation sequencing (RNA-Seq) was utilized to determine the OpaR regulon by comparing strains BB22OP (opaR+, LM5312) and BB22TR (∆opaR1, LM5674). This work, using the published V. parahaemolyticus BB22OP genome sequence, confirms and expands upon a previous microarray analysis for these two strains that used an Affymetrix GeneChip designed from the closely related V. parahaemolyticus RIMD2210633 genome sequence. Overall there was excellent correlation between the microarray and RNA-Seq data. Eleven transcription factors under OpaR control were identified by both methods and further confirmed by quantitative reverse transcription PCR (qRT-PCR) analysis. Nine of these transcription factors were demonstrated to be direct OpaR targets via in vitro electrophoretic mobility shift assays with purified hexahistidine-tagged OpaR. Identification of the direct and indirect targets of OpaR, including small RNAs, will enable the construction of a network map of regulatory interactions important for the switch between the nonpathogenic and pathogenic states.

Use of Whole-Genome Phylogeny and Comparisons for Development of a Multiplex PCR Assay To Identify Sequence Type 36 Vibrio parahaemolyticus

Vibrio parahaemolyticus sequence type 36 (ST36) strains that are native to the Pacific Ocean have recently caused multistate outbreaks of gastroenteritis linked to shellfish harvested from the Atlantic Ocean. Whole-genome comparisons of 295 genomes of V. parahaemolyticus, including several traced to northeastern U.S. sources, were used to identify diagnostic loci, one putatively encoding an endonuclease (prp), and two others potentially conferring O-antigenic properties (cps and flp). The combination of all three loci was present in only one clade of closely related strains of ST36, ST59, and one additional unknown sequence type. However, each locus was also identified outside this clade, with prp and flp occurring in only two nonclade isolates and cps in four. Based on the distribution of these loci in sequenced genomes, prp identified clade strains with >99% accuracy, but the addition of one more locus increased accuracy to 100%. Oligonucleotide primers targeting prp and cps were combined in a multiplex PCR method that defines species using the tlh locus and determines the presence of both the tdh and trh hemolysin-encoding genes, which are also present in ST36. Application of the method in vitro to a collection of 94 clinical isolates collected over a 4-year period in three northeastern U.S. states and 87 environmental isolates revealed that the prp and cps amplicons were detected only in clinical isolates identified as belonging to the ST36 clade and in no environmental isolates from the region. The assay should improve detection and surveillance, thereby reducing infections.

Insights into the environmental reservoir of pathogenic Vibrio parahaemolyticus using comparative genomics

Vibrio parahaemolyticus is an aquatic halophilic bacterium that occupies estuarine and coastal marine environments, and is a leading cause of seafood-borne food poisoning cases. To investigate the environmental reservoir and potential gene flow that occurs among V. parahaemolyticus isolates, the virulence-associated gene content and genome diversity of a collection of 133 V. parahaemolyticus isolates were analyzed. Phylogenetic analysis of housekeeping genes, and pulsed-field gel electrophoresis, demonstrated that there is genetic similarity among V. parahaemolyticus clinical and environmental isolates. Whole-genome sequencing and comparative analysis of six representative V. parahaemolyticus isolates was used to identify genes that are unique to the clinical and environmental isolates examined. Comparative genomics demonstrated an O3:K6 environmental isolate, AF91, which was cultured from sediment collected in Florida in 2006, has significant genomic similarity to the post-1995 O3:K6 isolates. However, AF91 lacks the majority of the virulence-associated genes and genomic islands associated with these highly virulent post-1995 O3:K6 genomes. These findings demonstrate that although they do not contain most of the known virulence-associated regions, some V. parahaemolyticus environmental isolates exhibit significant genetic similarity to clinical isolates. This highlights the dynamic nature of the V. parahaemolyticus genome allowing them to transition between aquatic and host-pathogen states.

Characterization of Vibrio parahaemolyticus clinical strains from Maryland (2012-2013) and comparisons to a locally and globally diverse V. parahaemolyticus strains by whole-genome sequence analysis

Vibrio parahaemolyticus is the leading cause of foodborne illnesses in the US associated with the consumption of raw shellfish. Previous population studies of V. parahaemolyticus have used Multi-Locus Sequence Typing (MLST) or Pulsed Field Gel Electrophoresis (PFGE). Whole genome sequencing (WGS) provides a much higher level of resolution, but has been used to characterize only a few United States (US) clinical isolates. Here we report the WGS characterization of 34 genomes of V. parahaemolyticus strains that were isolated from clinical cases in the state of Maryland (MD) during 2 years (2012-2013). These 2 years saw an increase of V. parahaemolyticus cases compared to previous years. Among these MD isolates, 28% were negative for tdh and trh, 8% were tdh positive only, 11% were trh positive only, and 53% contained both genes. We compared this set of V. parahaemolyticus genomes to those of a collection of 17 archival strains from the US (10 previously sequenced strains and 7 from NCBI, collected between 1988 and 2004) and 15 international strains, isolated from geographically-diverse environmental and clinical sources (collected between 1980 and 2010). A WGS phylogenetic analysis of these strains revealed the regional outbreak strains from MD are highly diverse and yet genetically distinct from the international strains. Some MD strains caused outbreaks 2 years in a row, indicating a local source of contamination (e.g., ST631). Advances in WGS will enable this type of analysis to become routine, providing an excellent tool for improved surveillance. Databases built with phylogenetic data will help pinpoint sources of contamination in future outbreaks and contribute to faster outbreak control.

Complete Genome Sequence of Vibrio parahaemolyticus Environmental Strain UCM-V493

Vibrio parahaemolyticus is the leading bacterial cause of seafood-related gastroenteritis in the world. Here, we report the complete genome sequence and annotation of an environmental strain of V. parahaemolyticus, UCM-V493, with the aim of understanding the differences between the clinical and environmental isolates of the bacteria. We also make some preliminary sequence comparisons with the clinical strain RIMD2210633.

Distribution and dynamics of epidemic and pandemic Vibrio parahaemolyticus virulence factors

Vibrio parahaemolyticus, autochthonous to estuarine, marine, and coastal environments throughout the world, is the causative agent of food-borne gastroenteritis. More than 80 serotypes have been described worldwide, based on antigenic properties of the somatic (O) and capsular (K) antigens. Serovar O3:K6 emerged in India in 1996 and subsequently was isolated worldwide, leading to the conclusion that the first V. parahaemolyticus pandemic had taken place. Most strains of V. parahaemolyticus isolated from the environment or seafood, in contrast to clinical strains, do not produce a thermostable direct hemolysin (TDH) and/or a TDH-related hemolysin (TRH). Type 3 secretion systems (T3SSs), needle-like apparatuses able to deliver bacterial effectors into host cytoplasm, were identified as triggering cytotoxicity and enterotoxicity. Type 6 secretion systems (T6SS) predicted to be involved in intracellular trafficking and vesicular transport appear to play a role in V. parahaemolyticus virulence. Recent advances in V. parahaemolyticus genomics identified several pathogenicity islands (VpaIs) located on either chromosome in both epidemic and pandemic strains and comprising additional colonization factors, such as restriction-modification complexes, chemotaxis proteins, classical bacterial surface virulence factors, and putative colicins. Furthermore, studies indicate strains lacking toxins and genomic regions associated with pathogenicity may also be pathogenic, suggesting other important virulence factors remain to be identified. The unique repertoire of virulence factors identified to date, their occurrence and distribution in both epidemic and pandemic strains worldwide are described, with the aim of highlighting the complexity of V. parahaemolyticus pathogenicity as well as its dynamic genome.

Draft Genome Sequence of Vibrio parahaemolyticus V110, Isolated from Shrimp in Hong Kong

We report the whole-genome sequence of a tdh- and trh-negative Vibrio parahaemolyticus strain, V110, from shrimp. The major difference of V110 from clinical strains was its lack of the type III secretion system T3SS2, a key component of virulence. Further sequence comparison can shed light on the pathogenesis of V. parahaemolyticus.