Genomic island identification in Vibrio vulnificus reveals significant genome plasticity in this human pathogen

A novel Diguanylate cyclase VdcR has multifaceted regulatory functions in the pathogenicity of Vibrio vulnificus

BACKGROUND

Vibrio vulnificus is a Gram-negative pathogen that infects humans through foodborne or wound infections. Victims of V. vulnificus infections face significant health risks, including cellulitis and septicemia, which have rapid disease progression and high mortality rates. Diguanylate cyclase is responsible for producing the secondary messenger cyclic di-GMP. It plays a crucial role in regulating various bacterial physiological processes, such as motility, toxicity, and pathogenicity, through transcriptional regulation and affecting cyclic di-GMP levels. However, the DGC-mediated pathogenicity regulation in V. vulnificus is still unclear.

METHODS

The vdcR gene in V. vulnificus was studied using a deletion strain (ΔVdcR) and an overexpression strain (oeVdcR) to understand its role in regulating the bacterium's pathogenicity. The electrophoretic mobility shift assay and RT-qPCR confirmed VdcR's impact on phosphodiesterase gene expression. To investigate how VdcR affects pathogenicity, V. vulnificus variant strains were assays for hemolysis, metalloprotease activity, cytotoxicity, resistance to phagocytosis, and lethality assays of the nematode Caenorhabditis elegans after infection.

RESULTS

This study discovered a virulence-associated diguanylate cyclase, VdcR, which serves as a transcriptional regulator to induce phosphodiesterases and reduce the accumulation of cyclic di-GMP. VdcR expression resulted in low hemolysis, metalloprotease, and cytotoxicity activity. It also improved the cell adhesion ability and anti-phagocytosis activity to infect the host cell and escape the macrophage phagocytosis. The constitutively expressed VdcR in V. vulnificus caused low mortality rates in Caenorhabditis elegans survival assays.

CONCLUSION

The above evidence demonstrated that VdcR suppresses the pathogenicity in V. vulnificus YJ016.

A novel genomic island encodes vibrioferrin synthesis in the marine pathogen Photobacterium damselae subsp. damselae

In this study, we identified and analyzed a novel genomic island (GI), named pddGI-1, located on chromosome II of certain strains of the marine pathogen Photobacterium damselae subsp. damselae (Pdd). This GI shares structural similarities with other GIs found in Vibrio species, such as the Vibrio seventh pandemic island-II (VSP-II) of V. cholerae. The pddGI-1 island is a mosaic of gene blocks that encode functions related to ROS defense, anaerobic energy metabolism, and restriction-modification (RM) systems. Notably, pddGI-1 also includes a complete vibrioferrin siderophore system, enabling the bacteria to thrive in low-iron environments. Vibrioferrin was chemically identified from cell-free supernatants of Pdd RG91. Additionally, a pvsD mutant deficient in vibrioferrin biosynthesis was generated and analyzed. The results suggest that Pdd strains harbouring pddGI-1 gain a distinct growth advantage under iron-limited conditions. These findings, along with previous research, highlight the significant heterogeneity in iron assimilation systems among Pdd strains.

Genome analysis of clinical genotype Vibrio vulnificus isolated from seafood in Mangaluru Coast, India provides insights into its pathogenicity

Vibrio vulnificus an opportunistic human pathogen native to marine/estuarine environment, is one of the leading causes of death due to seafood consumption and exposure of wounds to seawater worldwide. The present study involves the whole genome sequence analysis of an environmental strain of V. vulnificus (clinical genotype) isolated from seafood along the Mangaluru coast of India. The sequenced genome data was subjected to in-silico analysis of phylogeny, virulence genes, antimicrobial resistance determinants, and secretary proteins using suitable bioinformatics tools. The sequenced isolate had an overall genome length of 4.8 Mb and GC content of 46% with 4400 coding DNA sequences. The sequenced strain belongs to a new sequence type (Multilocus sequence typing) and was also found to branch with a phylogenetic lineage that groups the most infectious strains of V. vulnificus. The seafood isolate had complete genes involved in conferring serum resistance yet showed limited serum resistance. The study identified several genes against the antibiotics that are commonly used in their treatment, highlighting the need for alternative treatments. Also, the secretory protein analysis revealed genes associated with major pathways like ABC transporters, two-component systems, quorum sensing, biofilm formation, cationic antimicrobial peptide (CAMP) resistance, and others that play a critical role in the pathogenesis of the V. vulnificus. To the best of our knowledge, this is the first report of a detailed analysis of the genomic information of a V. vulnificus isolated from the Indian subcontinent and provides evidence that raises public health concerns about the safety of seafood.

Genomic Signature in Evolutionary Biology: A Review

Organisms are unique physical entities in which information is stored and continuously processed. The digital nature of DNA sequences enables the construction of a dynamic information reservoir. However, the distinction between the hardware and software components in the information flow is crucial to identify the mechanisms generating specific genomic signatures. In this work, we perform a bibliometric analysis to identify the different purposes of looking for particular patterns in DNA sequences associated with a given phenotype. This study has enabled us to make a conceptual breakdown of the genomic signature and differentiate the leading applications. On the one hand, it refers to gene expression profiling associated with a biological function, which may be shared across taxa. This signature is the focus of study in precision medicine. On the other hand, it also refers to characteristic patterns in species-specific DNA sequences. This interpretation plays a key role in comparative genomics, identifying evolutionary relationships. Looking at the relevant studies in our bibliographic database, we highlight the main factors causing heterogeneities in genome composition and how they can be quantified. All these findings lead us to reformulate some questions relevant to evolutionary biology.

Phylogeography of the marine pathogen, Vibrio vulnificus, revealed the ancestral scenarios of its evolution

Vibrio vulnificus is the leading cause of seafood‐associated deaths worldwide. Despite the growing knowledge about the population structure of V. vulnificus, the evolutionary history and the ancestral relationships of strains isolated from various regions around the world have not been determined. Using the largest collection of sequence and isolate data of V. vulnificus to date, we applied ancestral character reconstruction to study the phylogeography of V. vulnificus. Multilocus sequence typing data from 10 housekeeping genes were used for the inference of ancestral states and reconstruction of the evolutionary history. The findings showed that the common ancestor of all V. vulnificus populations originated from East Asia, and later evolved into two main clusters that spread with time and eventually evolved into distinct populations in different parts of the world. While we found no meaningful insights concerning the evolution of V. vulnificus populations in the Middle East; however, we were able to reconstruct the ancestral scenarios of its evolution in East Asia, North America, and Western Europe.

Occurrence and significance of pathogenicity and fitness islands in environmental vibrios

Pathogenicity islands (PAIs) are large genomic regions that contain virulence genes, which aid pathogens in establishing infections. While PAIs in clinical strains (strains isolated from a human infection) are well-studied, less is known about the occurrence of PAIs in strains isolated from the environment. In this study we describe three PAIs found in environmental Vibrio vulnificus and Vibrio parahaemolyticus strains, as well as a genomic fitness island found in a Vibrio diabolicus strain. All four islands had markedly different GC profiles than the rest of the genome, indicating that all of these islands were acquired via lateral gene transfer. Genes on the PAIs and fitness island were characterized. The PAI found in V. parahaemolyticus contained the tdh gene, a collagenase gene, and genes involved in the type 3 secretion system II (T3SS2). A V. vulnificus environmental strain contained two PAIs, a small 25 kbp PAI and a larger 143 kbp PAI. Both PAIs contained virulence genes. Toxin-antitoxin (TA) genes were found in all three species: on the V. diabolicus fitness island, and on the V. parahaemolyticus and V. vulnificus PAIs.

Patterns of triclosan resistance in Vibrionaceae

The antimicrobial additive triclosan has been used in personal care products widely across the globe for decades. Triclosan resistance has been noted among Vibrio spp., but reports have been anecdotal and the extent of phenotypic triclosan resistance across the Vibrionaceae family has not been established. Here, triclosan resistance was determined for Vibrionaceae strains across nine distinct clades. Minimum inhibitory concentrations (MIC) were determined for 70 isolates from clinical (n = 6) and environmental sources (n = 64); only two were susceptible to triclosan. The mean MIC for all resistant Vibrionaceae was 53 µg mL-1 (range 3.1-550 µg mL-1), but was significantly different between clades (p < 0.001). The highest mean triclosan MIC was observed in the Splendidus clade (200 µg mL-1; n = 3). Triclosan mean MICs were 68.8 µg mL-1 in the Damselae clade and 45.3 µg mL-1 in the Harveyi clade. The lowest mean MIC was observed in the Cholerae clade with 14.4 µg mL-1, which was primarily represented by clinical strains. There were no significant differences in triclosan MIC among individual species or among environmental strains isolated from different locations. Overall, phenotypic triclosan resistance appears to be widespread across multiple clades of Vibrionaceae.

Eco-evolutionary Dynamics Linked to Horizontal Gene Transfer in Vibrios

Vibrio is a genus of ubiquitous heterotrophic bacteria found in aquatic environments. Although they are a small percentage of the bacteria in these environments, vibrios can predominate during blooms. Vibrios also play important roles in the degradation of polymeric substances, such as chitin, and in other biogeochemical processes. Vibrios can be found as free-living bacteria, attached to particles, or associated with other organisms in a mutualistic, commensal, or pathogenic relationship. This review focuses on vibrio ecology and genome plasticity, which confers an ability to adapt to new niches and is driven, at least in part, by horizontal gene transfer (HGT). The extent of HGT and its role in pathogen emergence are discussed based on genomic studies of environmental and pathogenic vibrios, mobile genetically encoded virulence factors, and mechanistic studies on the different modes of HGT.

Post-Genomic Analysis of Members of the Family Vibrionaceae

Similar to other genera and species of bacteria, whole genomic sequencing has revolutionized how we think about and address questions of basic Vibrio biology. In this review we examined 36 completely sequenced and annotated members of the Vibrionaceae family, encompassing 12 different species of the genera Vibrio, Aliivibrio, and Photobacterium. We reconstructed the phylogenetic relationships among representatives of this group of bacteria by using three housekeeping genes and 16S rRNA sequences. With an evolutionary framework in place, we describe the occurrence and distribution of primary and alternative sigma factors, global regulators present in all bacteria. Among Vibrio we show that the number and function of many of these sigma factors differs from species to species. We also describe the role of the Vibrio-specific regulator ToxRS in fitness and survival. Examination of the biochemical capabilities was and still is the foundation of classifying and identifying new Vibrio species. Using comparative genomics, we examine the distribution of carbon utilization patterns among Vibrio species as a possible marker for understanding bacteria-host interactions. Finally, we discuss the significant role that horizontal gene transfer, specifically, the distribution and structure of integrons, has played in Vibrio evolution.

Pathogenicity Island Cross Talk Mediated by Recombination Directionality Factors Facilitates Excision from the Chromosome

Pathogenicity islands (PAIs) are mobile integrated genetic elements (MIGEs) that contain a diverse range of virulence factors and are essential in the evolution of pathogenic bacteria. PAIs are widespread among bacteria and integrate into the host genome, commonly at a tRNA locus, via integrase-mediated site-specific recombination. The excision of PAIs is the first step in the horizontal transfer of these elements and is not well understood. In this study, we examined the role of recombination directionality factors (RDFs) and their relationship with integrases in the excision of two PAIs essential for Vibrio cholerae host colonization: Vibrio pathogenicity island 1 (VPI-1) and VPI-2. VPI-1 does not contain an RDF, which allowed us to answer the question of whether RDFs are an absolute requirement for excision. We found that an RDF was required for efficient excision of VPI-2 but not VPI-1 and that RDFs can induce excision of both islands. Expression data revealed that the RDFs act as transcriptional repressors to both VPI-1- and VPI-2-encoded integrases. We demonstrated that the RDFs Vibrio excision factor A (VefA) and VefB bind at the attachment sites (overlapping the int promoter region) of VPI-1 and VPI-2, thus supporting this mode of integrase repression. In addition, V. cholerae RDFs are promiscuous due to their dual functions of promoting excision of both VPI-1 and VPI-2 and acting as negative transcriptional regulators of the integrases. This is the first demonstration of cross talk between PAIs mediated via RDFs which reveals the complex interactions that occur between separately acquired MIGEs.

IMPORTANCE

Deciphering the mechanisms of pathogenicity island excision is necessary for understanding the evolution and spread of these elements to their nonpathogenic counterparts. Such mechanistic insight would assist in predicting the mobility of uncharacterized genetic elements. This study identified extensive RDF-mediated cross talk between two nonhomologous VPIs and demonstrated the dual functionality of RDF proteins: (i) inducing PAI excision and (ii) acting as transcriptional regulators. Findings from this study may be implicated in determining the mobilome contribution of other bacteria with multiple MIGEs.

Draft Genome Sequence of the Pathogenic Bacterium Vibrio vulnificus V252 Biotype 1, Isolated in Israel

We report the genome sequence of the pathogenic Vibrio vulnificus biotype 1 clade B, which is suggested to have a common ancestor with biotype 3. This draft genome of the clinical strain V252, isolated in Israel, represents the clonal clade B group that contains both clinical and environmental strains.

Draft Genome Sequence of Environmental Bacterium Vibrio vulnificus CladeA-yb158

We report the genome sequence of the environmental Vibrio vulnificus biotype 1_cladeA. This draft genome of the CladeA-yb158 strain, isolated in Israel, represents this newly emerged clonal group that contains both clinical and environmental strains.

Genome-wide SNP-genotyping array to study the evolution of the human pathogen Vibrio vulnificus biotype 3

Vibrio vulnificus is an aquatic bacterium and an important human pathogen. Strains of V. vulnificus are classified into three different biotypes. The newly emerged biotype 3 has been found to be clonal and restricted to Israel. In the family Vibrionaceae, horizontal gene transfer is the main mechanism responsible for the emergence of new pathogen groups. To better understand the evolution of the bacterium, and in particular to trace the evolution of biotype 3, we performed genome-wide SNP genotyping of 254 clinical and environmental V. vulnificus isolates with worldwide distribution recovered over a 30-year period, representing all phylogeny groups. A custom single-nucleotide polymorphism (SNP) array implemented on the Illumina GoldenGate platform was developed based on 570 SNPs randomly distributed throughout the genome. In general, the genotyping results divided the V. vulnificus species into three main phylogenetic lineages and an additional subgroup, clade B, consisting of environmental and clinical isolates from Israel. Data analysis suggested that 69% of biotype 3 SNPs are similar to SNPs from clade B, indicating that biotype 3 and clade B have a common ancestor. The rest of the biotype 3 SNPs were scattered along the biotype 3 genome, probably representing multiple chromosomal segments that may have been horizontally inserted into the clade B recipient core genome from other phylogroups or bacterial species sharing the same ecological niche. Results emphasize the continuous evolution of V. vulnificus and support the emergence of new pathogenic groups within this species as a recurrent phenomenon. Our findings contribute to a broader understanding of the evolution of this human pathogen.

Sequence and expression divergence of an ancient duplication of the chaperonin groESEL operon in Vibrio species

Heat-shock proteins are molecular chaperones essential for protein folding, degradation and trafficking. The human pathogen Vibrio vulnificus encodes a copy of the groESEL operon in both chromosomes and these genes share <80 % similarity with each other. Comparative genomic analysis was used to determine whether this duplication is prevalent among Vibrionaceae specifically or Gammaproteobacteria in general. Among the Vibrionaceae complete genome sequences in the database (31 species), seven Vibrio species contained a copy of groESEL in each chromosome, including the human pathogens Vibrio cholerae, Vibrio parahaemolyticus and V. vulnificus. Phylogenetic analysis of GroEL among the Gammaproteobacteria indicated that GroESEL-1 encoded in chromosome I was the ancestral copy and GroESEL-2 in chromosome II arose by an ancient gene duplication event. Interestingly, outside of the Vibrionaceae within the Gammaproteobacteria, groESEL chromosomal duplications were rare among the 296 genomes examined; only five additional species contained two or more copies. Examination of the expression pattern of groEL from V. vulnificus cells grown under different conditions revealed differential expression between the copies. The data demonstrate that groEL-1 was more highly expressed during growth in exponential phase than groEL-2 and a similar pattern was also found in both V. cholerae and V. parahaemolyticus. Overall these data suggest that retention of both copies of groESEL in Vibrio species may confer an evolutionary advantage.

Insight into the evolution of Vibrio vulnificus biotype 3's genome

Vibrio vulnificus is an aquatic bacterium and an important human pathogen. Strains of V. vulnificus are biochemically classified into three biotypes. The newly emerged biotype 3 appears to be rather clonal and geographically restricted to Israel, where it caused an outbreak of wound infections and bacteremia. To understand the evolution of the bacterium's genome, we sequenced and analyzed the genome of biotype 3 strain VVyb1(BT3), and then conducted a microbial environmental survey of the hypothesized niche from which it probably evolved. The genome of this environmental isolate revealed higher similarity to the published biotype 1 genomes of clinical strains (90%) than to the environmental strains (87%), supporting the virulence of the biotype 3 group. Moreover, 214 of the total 5361 genes were found to be unique to strain VVyb1(BT3), having no sequence similarity to any of the known genomes of V. vulnificus; 35 of them function in DNA mobility and rearrangement, supporting the role of horizontal gene transfer in genome evolution. Interestingly, 29 of the “unique” genes had homologies among Shewanella species. In a survey conducted in aquaculture ponds in Israel, we successfully co-isolated Shewanella and V. vulnificus from the same niche, further supporting the probable contribution of Shewanella to the genome evolution of biotype 3. Indeed, one gene was found in a S. algae isolate. Surprisingly, molecular analysis revealed that some of the considered unique genes are harbored by non-sequenced biotype 1 strains isolated from the same environment. Finally, analyses of the biotype 3 genome together with the environmental survey suggested that its genome originated from a biotype 1 Israeli strain that acquired a rather small number of genes from other bacterial species in the niche, such as Shewanella. Therefore, aquaculture is likely to play a major role as a man-made ecological niche in bacterial evolution, leading the emergence of new pathogenic groups in V. vulnificus.

Fitness factors in vibrios: a mini-review

Vibrios are Gram-negative curved bacilli that occur naturally in marine, estuarine, and freshwater systems. Some species include human and animal pathogens, and some vibrios are necessary for natural systems, including the carbon cycle and osmoregulation. Countless in vivo and in vitro studies have examined the interactions between vibrios and their environment, including molecules, cells, whole animals, and abiotic substrates. Many studies have characterized virulence factors, attachment factors, regulatory factors, and antimicrobial resistance factors, and most of these factors impact the organism's fitness regardless of its external environment. This review aims to identify common attributes among factors that increase fitness in various environments, regardless of whether the environment is an oyster, a rabbit, a flask of immortalized mammalian cells, or a planktonic chitin particle. This review aims to summarize findings published thus far to encapsulate some of the basic similarities among the many vibrio fitness factors and how they frame our understanding of vibrio ecology. Factors representing these similarities include hemolysins, capsular polysaccharides, flagella, proteases, attachment factors, type III secretion systems, chitin binding proteins, iron acquisition systems, and colonization factors.

Genome Sequence of the Pathogenic Bacterium Vibrio vulnificus Biotype 3

We report the first genome sequence of the pathogenic Vibrio vulnificus biotype 3. This draft genome sequence of the environmental strain VVyb1(BT3), isolated in Israel, provides a representation of this newly emerged clonal group, which reveals higher similarity to the clinical strains of biotype 1 than to the environmental ones.

Genome sequence of the human-pathogenetic bacterium Vibrio vulnificus B2

Vibrio vulnificus, which can lead to rapidly expanding cellulitis or septicemia, is present in the marine environment. Here, we present the draft genome sequence of strain B2, which was isolated from a septicemia patient in 2010.

An extracellular serine protease produced by Vibrio vulnificus NCIMB 2137, a metalloprotease-gene negative strain isolated from a diseased eel

Vibrio vulnificus is a ubiquitous estuarine microorganism but causes fatal systemic infections in immunocompromised humans, cultured eels or shrimps. An extracellular metalloprotease VVP/VvpE has been reported to be a potential virulence factor of the bacterium; however, a few strains isolated from a diseased eel or shrimp were recently found to produce a serine protease termed VvsA, but not VVP/VvpE. In the present study, we found that these strains had lost the 80 kb genomic region including the gene encoding VVP/VvpE. We also purified VvsA from the culture supernatant through ammonium sulfate fractionation, gel filtration and ion-exchange column chromatography, and the enzyme was demonstrated to be a chymotrypsin-like protease, as well as those from some vibrios. The gene vvsA was shown to constitute an operon with a downstream gene vvsB, and several Vibrio species were found to have orthologues of vvsAB. These findings indicate that the genes vvp/vvpE and vvsAB might be mobile genetic elements.

Integrative genome-scale metabolic analysis of Vibrio vulnificus for drug targeting and discovery

Chromosome 1 of Vibrio vulnificus tends to contain larger portion of essential or housekeeping genes on the basis of the genomic analysis and gene knockout experiments performed in this study, while its chromosome 2 seems to have originated and evolved from a plasmid.

The genome-scale metabolic network model of V. vulnificus was reconstructed based on databases and literature, and was used to identify 193 essential metabolites.

Five essential metabolites finally selected after the filtering process are 2-amino-4-hydroxy-6-hydroxymethyl-7,8-dihydropteridine (AHHMP), D-glutamate (DGLU), 2,3-dihydrodipicolinate (DHDP), 1-deoxy-D-xylulose 5-phosphate (DX5P), and 4-aminobenzoate (PABA), which were predicted to be essential in V. vulnificus, absent in human, and are consumed by multiple reactions.

Chemical analogs of the five essential metabolites were screened and a hit compound showing the minimal inhibitory concentration (MIC) of 2 μg/ml and the minimal bactericidal concentration (MBC) of 4 μg/ml against V. vulnificus was identified.

Discovering new antimicrobial targets and consequently new antimicrobials is important as drug resistance of pathogenic microorganisms is becoming an increasingly serious problem in human healthcare management (Fischbach and Walsh, 2009). There clearly exists a gap between genomic studies and drug discovery as the accumulation of knowledge on pathogens at genome level has not successfully transformed into the development of effective drugs (Mills, 2006; Payne et al, 2007). In this study, we dissected the genome of a microbial pathogen in detail, and subsequently developed a systems biological strategy of employing genome-scale metabolic modeling and simulation together with metabolite essentiality analysis for effective drug targeting and discovery. This strategy was used for identifying new drug targets in an opportunistic pathogen Vibrio vulnificus CMCP6 as a model.

V. vulnificus is a Gram-negative halophilic bacterium that is found in estuarine waters, brackish ponds, or coastal areas, and its Biotype 1 is an opportunistic human pathogen that can attack immune-compromised patients, and causes primary septicemia, necrotized wound infections, and gastroenteritis. We previously found that many metabolic genes were specifically induced in vivo, suggesting that specific metabolic pathways are essential for in vivo survival and virulence of this pathogen (Kim et al, 2003; Lee et al, 2007). These results motivated us to carry out systems biological analysis of the genome and the metabolic network for new drug target discovery.

V. vulnificus CMCP6 has two chromosomes. We first re-sequenced genomic regions assembled in low quality and low depth, and subsequently re-annotated the whole genome of V. vulnificus. Horizontal gene transfer was suspected to be responsible for the diversification of each chromosome of V. vulnificus, and the presence of metabolic genes was more biased to chromosome 1 than chromosome 2. Further studies on V. vulnificus genome revealed that chromosome 2 is more prone to diversification for better adaptation to the environment than its chromosome 1, while chromosome 1 tends to expand their genetic repertoire while maintaining the core genes at a constant level.

Next, a genome-scale metabolic network VvuMBEL943 was reconstructed based on literature, databases and experiments for systematic studies on the metabolism of this pathogen and prediction of drug targets. The VvuMBEL943 model is composed of 943 reactions and 765 metabolites, and covers 673 genes. The model was validated by comparing its simulated cell growth phenotype obtained by constraints-based flux analysis with the V. vulnificus-specific experimental data previously reported in the literature. In this study, constraints-based flux analysis is an optimization-based simulation method that calculates intracellular fluxes under the specific genetic and environmental condition (Kim et al, 2008). As a result, 17 growth phenotypes were correctly predicted out of 18 cases, which demonstrate the validity of VvuMBEL943.

The main objective of constructing VvuMBEL943 in this study is to predict potential drug targets by system-wide analysis of the metabolic network for the effective treatment of V. vulnificus. To achieve this goal, a set of drug target candidates was predicted by taking a metabolite-centric approach. Metabolite essentiality analysis is a concept recently introduced for the study of cellular robustness to complement conventional reaction or gene-centric approach (Kim et al, 2007b). Metabolite essentiality analysis observes changes in flux distribution by removing each metabolite from the in silico metabolic network. Hence, metabolite essentiality predicts essential metabolites whose absence causes cell death. By selecting essential metabolites, it is possible to directly screen only their structural analogs, which substantially reduces the number of chemical compounds to screen from the chemical compound library. As a result of implementing this approach, 193 metabolites were initially identified to be essential to the cell. These essential metabolites were then further filtered based on the predetermined criteria, mainly organism specificity and multiple connectivity associated with each metabolite, in order to reduce the number of initial target candidates towards identifying the most effective ones.

Five essential metabolites finally selected are 2-amino-4-hydroxy-6-hydroxymethyl-7,8-dihydropteridine (AHHMP), D-glutamate (DGLU), 2,3-dihydrodipicolinate (DHDP), 1-deoxy-D-xylulose 5-phosphate (DX5P), and 4-aminobenzoate (PABA). Enzymes that consume these essential metabolites were experimentally verified to be essential, which indeed demonstrates the essentiality of these five metabolites. On the basis of the structural information of these five essential metabolites, whole-cell screening assay was performed using their analogs for possible antibacterial discovery. We screened 352 chemical analogs of the essential metabolites selected from the chemical compound library, and found a hit compound 24837, which shows the minimal inhibitory concentration (MIC) of 2 μg/ml and minimal bactericidal concentration (MBC) of 4 μg/ml, showing good antibacterial activity without further structural modification. Although this study demonstrates a proof-of-concept, the approaches and their rationale taken here should serve as a general strategy for discovering novel antibiotics and drugs based on systems-level analysis of metabolic networks.

Although the genomes of many microbial pathogens have been studied to help identify effective drug targets and novel drugs, such efforts have not yet reached full fruition. In this study, we report a systems biological approach that efficiently utilizes genomic information for drug targeting and discovery, and apply this approach to the opportunistic pathogen Vibrio vulnificus CMCP6. First, we partially re-sequenced and fully re-annotated the V. vulnificus CMCP6 genome, and accordingly reconstructed its genome-scale metabolic network, VvuMBEL943. The validated network model was employed to systematically predict drug targets using the concept of metabolite essentiality, along with additional filtering criteria. Target genes encoding enzymes that interact with the five essential metabolites finally selected were experimentally validated. These five essential metabolites are critical to the survival of the cell, and hence were used to guide the cost-effective selection of chemical analogs, which were then screened for antimicrobial activity in a whole-cell assay. This approach is expected to help fill the existing gap between genomics and drug discovery.

Integrating conjugative elements of the SXT/R391 family trigger the excision and drive the mobilization of a new class of Vibrio genomic islands

In vibrios and enterobacteria lateral gene transfer is often facilitated by integrating conjugative elements (ICEs) of the SXT/R391 family. SXT/R391 ICEs integrate by site-specific recombination into prfC and transfer by conjugation, a process that is initiated at a specific locus called the origin of transfer (oriT(SXT) ). We identified genomic islands (GIs) harbouring a sequence that shares >63% identity with oriT(SXT) in three species of Vibrio. Unlike SXT/R391 ICEs, these GIs are integrated into a gene coding for a putative stress-induced protein and do not appear to carry any gene coding for a conjugative machinery or for mobilization proteins. Our results show that SXT/R391 ICEs trigger the excision and mediate the conjugative transfer in trans of the three Vibrio GIs at high frequency. GIs' excision is independent of the ICE-encoded recombinase and is controlled by the ICE-encoded transcriptional activator SetCD, which is expressed during the host SOS response. Both mobI and traI, two ICE-borne genes involved in oriT recognition, are essential for GIs' transfer. We also found that SXT/R391 ICEs mobilize in trans over 1 Mb of chromosomal DNA located 5' of the GIs' integration site. Together these results support a novel mechanism of mobilization of GIs by ICEs of the SXT/R391 family.

Comparative genomics of the family Vibrionaceae reveals the wide distribution of genes encoding virulence-associated proteins

BACKGROUND

Species of the family Vibrionaceae are ubiquitous in marine environments. Several of these species are important pathogens of humans and marine species. Evidence indicates that genetic exchange plays an important role in the emergence of new pathogenic strains within this family. Data from the sequenced genomes of strains in this family could show how the genes encoded by all these strains, known as the pangenome, are distributed. Information about the core, accessory and panproteome of this family can show how, for example, genes encoding virulence-associated proteins are distributed and help us understand how virulence emerges.

RESULTS

We deduced the complete set of orthologs for eleven strains from this family. The core proteome consists of 1,882 orthologous groups, which is 28% of the 6,629 orthologous groups in this family. There were 4,411 accessory orthologous groups (i.e., proteins that occurred in from 2 to 10 proteomes) and 5,584 unique proteins (encoded once on only one of the eleven genomes). Proteins that have been associated with virulence in V. cholerae were widely distributed across the eleven genomes, but the majority was found only on the genomes of the two V. cholerae strains examined.

CONCLUSIONS

The proteomes are reflective of the differing evolutionary trajectories followed by different strains to similar phenotypes. The composition of the proteomes supports the notion that genetic exchange among species of the Vibrionaceae is widespread and that this exchange aids these species in adapting to their environments.

Vibrio vulnificus produces quorum sensing signals of the AHL-class

Vibrio vulnificus is an aquatic pathogenic bacterium that can cause vibriosis in humans and fish. The species is subdivided into three biotypes with the fish-virulent strains belonging to biotype 2. The quorum sensing (QS) phenomenon mediated by furanosyl borate diester or autoinducer 2 (AI-2) has been described in human strains of biotype 1, and here we show that the luxS gene which encodes AI-2 is present in all strains of V. vulnificus regardless of origin, biotype or serovar. In this study, we also demonstrate that V. vulnificus produces QS signals of the acylated homoserine lactone (AHL) class (AI-1). AHLs were detected in strains of biotype 1 and 2 from water, fish and human wound infections but not in strains isolated from human septicaemic cases. The AHL compound was identified as N-butanoyl-homoserine-lactone (C(4)-HL) by both reporter strains and by HPLC-high-resolution MS. C(4)-HL was detected when AHL-positive strains were grown in low-nutrient medium [modified sea water yeast extract (MSWYE)] but not in rich media (tryptic soy broth or brain-heart infusion) and its production was enhanced when blood factors were added to MSWYE. C(4)-HL was detected in vivo, in eels infected with AHL-positive biotype 2 strains. No known AHL-related gene was detected by PCR or Southern blot suggesting that AHL-related genes in V. vulnificus are different from those found in other Gram-negative bacteria.

Genomic islands are dynamic, ancient integrative elements in bacterial evolution

Acquisition of genomic islands plays a central part in bacterial evolution as a mechanism of diversification and adaptation. Genomic islands are non-self-mobilizing integrative and excisive elements that encode diverse functional characteristics but all contain a recombination module comprised of an integrase, associated attachment sites and, in some cases, a recombination directionality factor. Here, we discuss how a group of related genomic islands are evolutionarily ancient elements unrelated to plasmids, phages, integrons and integrative conjugative elements. In addition, we explore the diversity of genomic islands and their insertion sites among Gram-negative bacteria and discuss why they integrate at a limited number of tRNA genes.

Using Mahalanobis distance to compare genomic signatures between bacterial plasmids and chromosomes

Plasmids are ubiquitous mobile elements that serve as a pool of many host beneficial traits such as antibiotic resistance in bacterial communities. To understand the importance of plasmids in horizontal gene transfer, we need to gain insight into the ‘evolutionary history’ of these plasmids, i.e. the range of hosts in which they have evolved. Since extensive data support the proposal that foreign DNA acquires the host's nucleotide composition during long-term residence, comparison of nucleotide composition of plasmids and chromosomes could shed light on a plasmid's evolutionary history. The average absolute dinucleotide relative abundance difference, termed δ-distance, has been commonly used to measure differences in dinucleotide composition, or ‘genomic signature’, between bacterial chromosomes and plasmids. Here, we introduce the Mahalanobis distance, which takes into account the variance–covariance structure of the chromosome signatures. We demonstrate that the Mahalanobis distance is better than the δ-distance at measuring genomic signature differences between plasmids and chromosomes of potential hosts. We illustrate the usefulness of this metric for proposing candidate long-term hosts for plasmids, focusing on the virulence plasmids pXO1 from Bacillus anthracis, and pO157 from Escherichia coli O157:H7, as well as the broad host range multi-drug resistance plasmid pB10 from an unknown host.

On the origin and function of an insertion element VPaI-1 specific to post-1995 pandemic Vibrio parahaemolyticus strains

Since 1995, new virulent strains of Vibrio parahaemolyticus have emerged and spread throughout the world. These "pandemic" strains have four strain specific genomic islands (GIs), which are considered to be potential factors of the pandemicity. We investigated the origin and function of 24 genes in the so-called VPaI-1, one of the four GIs, by searching homologs in various species in Bacteria and Archaea. Of these 24 genes, two are found only in Vibrio vulnificus CMCP6 and Shewanella sp. MR-7. The genomic segment (- 8 kb) encompassing the two genes shows the synteny among the three species. Since many of the Shewanella species can grow at 4 degrees C, these two genes may be candidates of adaptation to temperature stress. Further, we found a candidate for a swarming gene, which is reported as the V. cholerae virulence gene. Based on these findings, we hypothesized the emergence of pandemicity and discuss the mechanism for how these strains spread throughout the world.

Three pathogenicity islands of Vibrio cholerae can excise from the chromosome and form circular intermediates

Vibrio pathogenicity island-2 (VPI-2) is a 57-kb region integrated at a transfer RNA (tRNA)-serine locus that encompasses VC1758 to VC1809 on the V. cholerae N16961 genome and is present in pandemic isolates. VPI-2 encodes a P4-like integrase, a restriction modification system, a Mu phage-like region, and a sialic acid metabolism region, as well as neuraminidase (VC1784), which is a glycosylhydrolase known to release sialic acid from sialoglycoconjugates to unmask GM1 gangliosides, the receptor for cholera toxin. We examined the tRNA-serine locus among the sequenced V. cholerae genomes and identified five variant VPI-2 regions, four of which retained the sialometabolism region. Three variant VPI-2 regions contained a type three secretion system. By using an inverse nested PCR approach, we found that the VPI-2 region can form an extrachromosomal circular intermediate (CI) molecule after precise excision from its tRNA-serine attachment site. We constructed a knockout mutant of VC1758 (int) with V. cholerae strain N16961 and found that no excision PCR product was produced, indicating that a functional cognate, VPI-2 integrase, is required for excision. The Vibrio seventh pandemic island-I (VSP-I) and VSP-II regions are present in V. cholerae O1 El Tor and O139 serogroup isolates. Novel regions are present at the VSP-I insertion site in strain MZO-3 and at the VSP-II insertion site in strain 623-39. VSP-II is a 27-kb region that integrates at a tRNA-methionine locus, is flanked by direct repeats, and encodes a P4-like integrase. We show that VSP-II can excise and form a CI and that the cognate VSP-II integrase is required for excision. Interestingly, VSP-I is not inserted at a tRNA locus and does encode a XerDC-like recombinase, but similar to VPI-2 and VSP-II, VSP-I does excise from the genome to form a CI. These results show that all three pathogenicity islands can excise from the chromosome, which is likely a first step in their horizontal transfer.

Emergence of a virulent clade of Vibrio vulnificus and correlation with the presence of a 33-kilobase genomic island

Vibrio vulnificus is a ubiquitous inhabitant of the marine coastal environment, and an important pathogen of humans. We characterized a globally distributed sample of environmental isolates from a range of habitats and hosts and compared these with isolates recovered from cases of human infection. Multilocus sequence typing data using six housekeeping genes divided 63 of the 67 isolates into the two main lineages previously noted for this species, and this division was also confirmed using the 16S rRNA and open reading frame VV0401 markers. Lineage I was comprised exclusively of biotype 1 isolates, whereas lineage II contained biotype 1 and all biotype 2 isolates. Four isolates did not cluster within either lineage: two biotype 3 and two biotype 1 isolates. The proportion of isolates recovered from a clinical setting was noted to be higher in lineage I than in lineage II. Lineage I isolates were also associated with a 33-kb genomic island (region XII), one of three regions identified by genome comparisons as unique to the species. Region XII contained an arylsulfatase gene cluster, a sulfate reduction system, two chondroitinase genes, and an oligopeptide ABC transport system, all of which are absent from the majority of lineage II isolates. Arylsulfatases and the sulfate reduction system, along with performing a scavenging role, have been hypothesized to play a role in pathogenic processes in other bacteria. Our data suggest that lineage I may have a higher pathogenic potential and that region XII, along with other regions, may give isolates a selective advantage either in the human host or in the aquatic environment or both.

Detection and transformation of genome segments that differ within a coastal population of Vibrio cholerae strains

Vibrio cholerae is an autochthonous member of diverse aquatic ecosystems around the globe. Collectively, the genomes of environmental V. cholerae strains comprise a large repository of encoded functions which can be acquired by individual V. cholerae lineages through uptake and recombination. To characterize the genomic diversity of environmental V. cholerae, we used comparative genome hybridization to study 41 environmental strains isolated from diverse habitats along the central California coast, a region free of endemic cholera. These data were used to classify genes of the epidemic V. cholerae O1 sequenced strain N16961 as conserved, variably present, or absent from the isolates. For the most part, absent genes were restricted to large mobile elements and have known functions in pathogenesis. Conversely, genes present in some, but not all, California isolates were in smaller contiguous clusters and were less likely to be near genes with functions in DNA mobility. Two such clusters of variable genes encoding different selectable metabolic phenotypes (mannose and diglucosamine utilization) were transformed into the genomes of environmental isolates by chitin-dependent competence, indicating that this mechanism of general genetic exchange is conserved among V. cholerae. The transformed DNA had an average size of 22.7 kbp, demonstrating that natural competence can mediate the movement of large chromosome fragments. Thus, whether variable genes arise through the acquisition of new sequences by horizontal gene transfer or by the loss of preexisting DNA though deletion, natural transformation provides a mechanism by which V. cholerae clones can gain access to the V. cholerae pan-genome.

A real-time PCR assay for the rapid determination of 16S rRNA genotype in Vibrio vulnificus

In a terminal restriction fragment polymorphism (T-RFLP) study, we recently reported a significant association between the type B 16S rRNA gene and clinical strains of Vibrio vulnificus associated with the consumption of raw oysters. In the present study we describe a real-time PCR assay for the rapid determination of the 16S rRNA type of V. vulnificus isolates. This assay was used to reexamine the 16S rRNA gene type in the strains studied previously by T-RFLP and additional isolates from selected sources. Analyses revealed that 15 of the strains (10 environmental and 5 clinical) previously found to be 16S rRNA type A actually appear to possess both the type A and B genes. The presence of both alleles was confirmed by cloning and sequencing both gene types from one strain. To our knowledge, this is the first report of 16S rRNA sequence heterogeneity within individual strains of V. vulnificus. The findings confirm the T-RFLP data that 16S rRNA type may be a useful marker for determining the clinical significance of V. vulnificus in disease in humans and cultured eels. The real-time PCR assay is much more rapid and less resource-intensive than T-RFLP, and should facilitate further study of the occurrence and distribution of the 16S rRNA genotypes of V. vulnificus. These studies should provide more definitive estimates of the risks associated with this organism and may lead to a better understanding of its virulence mechanism(s).

Patterns and mechanisms of genetic and phenotypic differentiation in marine microbes

Microbes in the ocean dominate biogeochemical processes and are far more diverse than anticipated. Thus, in order to understand the ocean system, we need to delineate microbial populations with predictable ecological functions. Recent observations suggest that ocean communities comprise diverse groups of bacteria organized into genotypic (and phenotypic) clusters of closely related organisms. Although such patterns are similar to metazoan communities, the underlying mechanisms for microbial communities may differ substantially. Indeed, the potential among ocean microbes for vast population sizes, extensive migration and both homologous and illegitimate genetic recombinations, which are uncoupled from reproduction, challenges classical population models primarily developed for sexually reproducing animals. We examine possible mechanisms leading to the formation of genotypic clusters and consider alternative population genetic models for differentiation at individual loci as well as gene content at the level of whole genomes. We further suggest that ocean bacteria follow at least two different adaptive strategies, which constrain rates and bounds of evolutionary processes: the 'opportunitroph', exploiting spatially and temporally variable resources; and the passive oligotroph, efficiently using low nutrient concentrations. These ecological lifestyle differences may represent a fundamental divide with major consequences for growth and predation rates, genome evolution and population diversity, as emergent properties driving the division of labour within microbial communities.

The genomic code: inferring Vibrionaceae niche specialization

The Vibrionaceae show a wide range of niche specialization, from free-living forms to those attached to biotic and abiotic surfaces, from symbionts to pathogens and from estuarine inhabitants to deep-sea piezophiles. The existence of complete genome sequences for closely related species from varied aquatic niches makes this group an excellent case study for genome comparison.

Four genomic islands that mark post-1995 pandemic Vibrio parahaemolyticus isolates

Background

Vibrio parahaemolyticus is an aquatic, halophilic, Gram-negative bacterium, first discovered in 1950 in Japan during a food-poisoning outbreak. Infections resulting from consumption of V. parahaemolyticus have increased globally in the last 10 years leading to the bacterium's classification as a newly emerging pathogen. In 1996 the first appearance of a pandemic V. parahaemolyticus clone occurred, a new O3:K6 serotype strain that has now been identified worldwide as a major cause of seafood-borne gastroenteritis.

Results

We examined the sequenced genome of V. parahaemolyticus RIMD2210633, an O3:K6 serotype strain isolated in Japan in 1996, by bioinformatic analyses to uncover genomic islands (GIs) that may play a role in the emergence and pathogenesis of pandemic strains. We identified 7 regions ranging in size from 10 kb to 81 kb that had the characteristics of GIs such as aberrant base composition compared to the core genome, presence of phage-like integrases, flanked by direct repeats and the absence of these regions from closely related species. Molecular analysis of worldwide clinical isolates of V. parahaemolyticus recovered over the last 33 years demonstrated that a 24 kb region named V. parahaemolyticus island-1 (VPaI-1) encompassing ORFs VP0380 to VP0403 is only present in new O3:K6 and related strains recovered after 1995. We investigated the presence of 3 additional regions, VPaI-4 (VP2131 to VP2144), VPaI-5 (VP2900 to VP2910) and VPaI-6 (VPA1254 to VPA1270) by PCR assays and Southern blot analyses among the same set of V. parahaemolyticus isolates. These 3 VPaI regions also gave similar distribution patterns amongst the 41 strains examined.

Conclusion

The 4 VPaI regions examined may represent DNA acquired by the pandemic group of V. parahaemolyticus isolates that increased their fitness either in the aquatic environment or in their ability to infect humans.