Comparative genome analysis of rice-pathogenic Burkholderia provides insight into capacity to adapt to different environments and hosts

Background

In addition to human and animal diseases, bacteria of the genus Burkholderia can cause plant diseases. The representative species of rice-pathogenic Burkholderia are Burkholderia glumae, B. gladioli, and B. plantarii, which primarily cause grain rot, sheath rot, and seedling blight, respectively, resulting in severe reductions in rice production. Though Burkholderia rice pathogens cause problems in rice-growing countries, comprehensive studies of these rice-pathogenic species aiming to control Burkholderia-mediated diseases are only in the early stages.

Results

We first sequenced the complete genome of B. plantarii ATCC 43733^T. Second, we conducted comparative analysis of the newly sequenced B. plantarii ATCC 43733^T genome with eleven complete or draft genomes of B. glumae and B. gladioli strains. Furthermore, we compared the genome of three rice Burkholderia pathogens with those of other Burkholderia species such as those found in environmental habitats and those known as animal/human pathogens. These B. glumae, B. gladioli, and B. plantarii strains have unique genes involved in toxoflavin or tropolone toxin production and the clustered regularly interspaced short palindromic repeats (CRISPR)-mediated bacterial immune system. Although the genome of B. plantarii ATCC 43733^T has many common features with those of B. glumae and B. gladioli, this B. plantarii strain has several unique features, including quorum sensing and CRISPR/CRISPR-associated protein (Cas) systems.

Conclusions

The complete genome sequence of B. plantarii ATCC 43733^T and publicly available genomes of B. glumae BGR1 and B. gladioli BSR3 enabled comprehensive comparative genome analyses among three rice-pathogenic Burkholderia species responsible for tissue rotting and seedling blight. Our results suggest that B. glumae has evolved rapidly, or has undergone rapid genome rearrangements or deletions, in response to the hosts. It also, clarifies the unique features of rice pathogenic Burkholderia species relative to other animal and human Burkholderia species.

Electronic supplementary material

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

Diversities in Virulence,Antifungal Activity,Pigmentation and DNA Fingerprint among Strains of Burkholderia glumae

Burkholderia glumae is the primary causal agent of bacterial panicle blight of rice. In this study, 11 naturally avirulent and nine virulent strains of B. glumae native to the southern United States were characterized in terms of virulence in rice and onion, toxofalvin production, antifungal activity, pigmentation and genomic structure. Virulence of B. glumae strains on rice panicles was highly correlated to virulence on onion bulb scales, suggesting that onion bulb can be a convenient alternative host system to efficiently determine the virulence of B. glumae strains. Production of toxoflavin, the phytotoxin that functions as a major virulence factor, was closely associated with the virulence phenotypes of B. glumae strains in rice. Some strains of B. glumae showed various levels of antifungal activity against Rhizoctonia solani, the causal agent of sheath blight, and pigmentation phenotypes on casamino acid-peptone-glucose (CPG) agar plates regardless of their virulence traits. Purple and yellow-green pigments were partially purified from a pigmenting strain of B. glumae, 411gr-6, and the purple pigment fraction showed a strong antifungal activity against Collectotrichum orbiculare. Genetic variations were detected among the B. glumae strains from DNA fingerprinting analyses by repetitive element sequence-based PCR (rep-PCR) for BOX-A1R-based repetitive extragenic palindromic (BOX) or enterobacterial repetitive intergenic consensus (ERIC) sequences of bacteria; and close genetic relatedness among virulent but pigment-deficient strains were revealed by clustering analyses of DNA fingerprints from BOX-and ERIC-PCR.     

Type VI secretion systems of plant-pathogenic Burkholderia glumae BGR1 play a functionally distinct role in interspecies interactions and virulence

In the environment, bacteria show close association, such as interspecies interaction, with other bacteria as well as host organisms. The type VI secretion system (T6SS) in gram-negative bacteria is involved in bacterial competition or virulence. The plant pathogen Burkholderia glumae BGR1, causing bacterial panicle blight in rice, has four T6SS gene clusters. The presence of at least one T6SS gene cluster in an organism indicates its distinct role, like in the bacterial and eukaryotic cell targeting system. In this study, deletion mutants targeting four tssD genes, which encode the main component of T6SS needle formation, were constructed to functionally dissect the four T6SSs in B. glumae BGR1. We found that both T6SS group_4 and group_5, belonging to the eukaryotic targeting system, act independently as bacterial virulence factors toward host plants. In contrast, T6SS group_1 is involved in bacterial competition by exerting antibacterial effects. The ΔtssD1 mutant lost the antibacterial effect of T6SS group_1. The ΔtssD1 mutant showed similar virulence as the wild-type BGR1 in rice because the ΔtssD1 mutant, like the wild-type BGR1, still has key virulence factors such as toxin production towards rice. However, metagenomic analysis showed different bacterial communities in rice infected with the ΔtssD1 mutant compared to wild-type BGR1. In particular, the T6SS group_1 controls endophytic plant-associated bacteria such as Luteibacter and Dyella in rice plants and may have an advantage in competing with endophytic plant-associated bacteria for settlement inside rice plants in the environment. Thus, B. glumae BGR1 causes disease using T6SSs with functionally distinct roles.     

Comparative and bioinformatics analyses of pathogenic bacterial secretomes identified by mass spectrometry in <Emphasis Type="Italic">Burkholderia</Emphasis> species

Secreted proteins (secretomes) play crucial roles during bacterial pathogenesis in both plant and human hosts. The identification and characterization of secretomes in the two plant pathogens Burkholderia glumae BGR1 and B. gladioli BSR3, which cause diseases in rice such as seedling blight, panicle blight, and grain rot, are important steps to not only understand the disease-causing mechanisms but also find remedies for the diseases. Here, we identified two datasets of secretomes in B. glumae BGR1 and B. gladioli BSR3, which consist of 118 and 111 proteins, respectively, using mass spectrometry approach and literature curation. Next, we characterized the functional properties, potential secretion pathways and sequence information properties of secretomes of two plant pathogens in a comparative analysis by various computational approaches. The ratio of potential non-classically secreted proteins (NCSPs) to classically secreted proteins (CSPs) in B. glumae BGR1 was greater than that in B. gladioli BSR3. For CSPs, the putative hydrophobic regions (PHRs) which are essential for secretion process of CSPs were screened in detail at their N-terminal sequences using hidden Markov model (HMM)–based method. Total 31 pairs of homologous proteins in two bacterial secretomes were indicated based on the global alignment (identity ≥ 70%). Our results may facilitate the understanding of the species-specific features of secretomes in two plant pathogenic Burkholderia species.     

Characterization of the core and accessory genomes of Pseudomonas aeruginosa using bioinformatic tools Spine and AGEnt

Background

Pseudomonas aeruginosa is an important opportunistic pathogen responsible for many infections in hospitalized and immunocompromised patients. Previous reports estimated that approximately 10% of its 6.6 Mbp genome varies from strain to strain and is therefore referred to as “accessory genome”. Elements within the accessory genome of P. aeruginosa have been associated with differences in virulence and antibiotic resistance. As whole genome sequencing of bacterial strains becomes more widespread and cost-effective, methods to quickly and reliably identify accessory genomic elements in newly sequenced P. aeruginosa genomes will be needed.

Results

We developed a bioinformatic method for identifying the accessory genome of P. aeruginosa. First, the core genome was determined based on sequence conserved among the completed genomes of twelve reference strains using Spine, a software program developed for this purpose. The core genome was 5.84 Mbp in size and contained 5,316 coding sequences. We then developed an in silico genome subtraction program named AGEnt to filter out core genomic sequences from P. aeruginosa whole genomes to identify accessory genomic sequences of these reference strains. This analysis determined that the accessory genome of P. aeruginosa ranged from 6.9-18.0% of the total genome, was enriched for genes associated with mobile elements, and was comprised of a majority of genes with unknown or unclear function. Using these genomes, we showed that AGEnt performed well compared to other publically available programs designed to detect accessory genomic elements. We then demonstrated the utility of the AGEnt program by applying it to the draft genomes of two previously unsequenced P. aeruginosa strains, PA99 and PA103.

Conclusions

The P. aeruginosa genome is rich in accessory genetic material. The AGEnt program accurately identified the accessory genomes of newly sequenced P. aeruginosa strains, even when draft genomes were used. As P. aeruginosa genomes become available at an increasingly rapid pace, this program will be useful in cataloging the expanding accessory genome of this bacterium and in discerning correlations between phenotype and accessory genome makeup. The combination of Spine and AGEnt should be useful in defining the accessory genomes of other bacterial species as well.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-737) contains supplementary material, which is available to authorized users.     

Comparison of 61 Sequenced Escherichia coli Genomes

Escherichia coli is an important component of the biosphere and is an ideal model for studies of processes involved in bacterial genome evolution. Sixty-one publically available E. coli and Shigella spp. sequenced genomes are compared, using basic methods to produce phylogenetic and proteomics trees, and to identify the pan- and core genomes of this set of sequenced strains. A hierarchical clustering of variable genes allowed clear separation of the strains into clusters, including known pathotypes; clinically relevant serotypes can also be resolved in this way. In contrast, when in silico MLST was performed, many of the various strains appear jumbled and less well resolved. The predicted pan-genome comprises 15,741 gene families, and only 993 (6%) of the families are represented in every genome, comprising the core genome. The variable or ‘accessory’ genes thus make up more than 90% of the pan-genome and about 80% of a typical genome; some of these variable genes tend to be co-localized on genomic islands. The diversity within the species E. coli, and the overlap in gene content between this and related species, suggests a continuum rather than sharp species borders in this group of Enterobacteriaceae.     

Evidence of a large novel gene pool associated with prokaryotic genomic islands

Microbial genes that are “novel” (no detectable homologs in other species) have become of increasing interest as environmental sampling suggests that there are many more such novel genes in yet-to-be-cultured microorganisms. By analyzing known microbial genomic islands and prophages, we developed criteria for systematic identification of putative genomic islands (clusters of genes of probable horizontal origin in a prokaryotic genome) in 63 prokaryotic genomes, and then characterized the distribution of novel genes and other features. All but a few of the genomes examined contained significantly higher proportions of novel genes in their predicted genomic islands compared with the rest of their genome (Paired t test = 4.43E-14 to 1.27E-18, depending on method). Moreover, the reverse observation (i.e., higher proportions of novel genes outside of islands) never reached statistical significance in any organism examined. We show that this higher proportion of novel genes in predicted genomic islands is not due to less accurate gene prediction in genomic island regions, but likely reflects a genuine increase in novel genes in these regions for both bacteria and archaea. This represents the first comprehensive analysis of novel genes in prokaryotic genomic islands and provides clues regarding the origin of novel genes. Our collective results imply that there are different gene pools associated with recently horizontally transmitted genomic regions versus regions that are primarily vertically inherited. Moreover, there are more novel genes within the gene pool associated with genomic islands. Since genomic islands are frequently associated with a particular microbial adaptation, such as antibiotic resistance, pathogen virulence, or metal resistance, this suggests that microbes may have access to a larger “arsenal” of novel genes for adaptation than previously thought.     

Involvement of a Quorum-Sensing-Regulated Lipase Secreted by a Clinical Isolate of Burkholderia glumae in Severe Disease Symptoms in Rice

Burkholderia glumae is an emerging rice pathogen in several areas around the world. Closely related Burkholderia species are important opportunistic human pathogens for specific groups of patients, such as patients with cystic fibrosis and patients with chronic granulomatous disease. Here we report that the first clinical isolate of B. glumae, strain AU6208, has retained its capability to be very pathogenic to rice. As previously reported for rice isolate B. glumae BGR1 (and also for the clinical isolate AU6208), TofI or TofR acyl homoserine lactone (AHL) quorum sensing played a pivotal role in rice virulence. We report that AHL quorum sensing in B. glumae AU6208 regulates secreted LipA lipase and toxoflavin, the phytotoxin produced by B. glumae. B. glumae AU6208 lipA mutants were no longer pathogenic to rice, indicating that the lipase is an important virulence factor. It was also established that type strain B. glumae ATCC 33617 did not produce toxoflavin and lipase and was nonpathogenic to rice. It was determined that in strain ATCC 33617 the LuxR family quorum-sensing sensor/regulator TofR was inactive. Introducing the tofR gene of B. glumae AU6208 in strain ATCC 33617 restored its ability to produce toxoflavin and the LipA lipase. This study extends the role of AHL quorum sensing in rice pathogenicity through the regulation of a lipase which was demonstrated to be a virulence factor. It is the first report of a clinical B. glumae isolate retaining strong rice pathogenicity and finally determined that B. glumae can undergo phenotypic conversion through a spontaneous mutation in the tofR regulator.     

A novel strategy for the identification of genomic islands by comparative analysis of the contents and contexts of tRNA sites in closely related bacteria

We devised software tools to systematically investigate the contents and contexts of bacterial tRNA and tmRNA genes, which are known insertion hotspots for genomic islands (GIs). The strategy, based on MAUVE-facilitated multigenome comparisons, was used to examine 87 Escherichia coli MG1655 tRNA and tmRNA genes and their orthologues in E.coli EDL933, E.coli CFT073 and Shigella flexneri Sf301. Our approach identified 49 GIs occupying ~1.7 Mb that mapped to 18 tRNA genes, missing 2 but identifying a further 30 GIs as compared with Islander [Y. Mantri and K. P. Williams (2004), Nucleic Acids Res., 32, D55–D58]. All these GIs had many strain-specific CDS, anomalous GC contents and/or significant dinucleotide biases, consistent with foreign origins. Our analysis demonstrated marked conservation of sequences flanking both empty tRNA sites and tRNA-associated GIs across all four genomes. Remarkably, there were only 2 upstream and 5 downstream deletions adjacent to the 328 loci investigated. In silico PCR analysis based on conserved flanking regions was also used to interrogate hotspots in another eight completely or partially sequenced E.coli and Shigella genomes. The tools developed are ideal for the analysis of other bacterial species and will lead to in silico and experimental discovery of new genomic islands.     

Sequence and functional analyses of Haemophilus spp. genomic islands

Background

A major part of horizontal gene transfer that contributes to the diversification and adaptation of bacteria is facilitated by genomic islands. The evolution of these islands is poorly understood. Some progress was made with the identification of a set of phylogenetically related genomic islands among the Proteobacteria, recognized from the investigation of the evolutionary origins of a Haemophilus influenzae antibiotic resistance island, namely ICEHin1056. More clarity comes from this comparative analysis of seven complete sequences of the ICEHin1056 genomic island subfamily.

Results

These genomic islands have core and accessory genes in approximately equal proportion, with none demonstrating recent acquisition from other islands. The number of variable sites within core genes is similar to that found in the host bacteria. Furthermore, the GC content of the core genes is similar to that of the host bacteria (38% to 40%). Most of the core gene content is formed by the syntenic type IV secretion system dependent conjugative module and replicative module. GC content and lack of variable sites indicate that the antibiotic resistance genes were acquired relatively recently. An analysis of conjugation efficiency and antibiotic susceptibility demonstrates that phenotypic expression of genomic island-borne genes differs between different hosts.

Conclusion

Genomic islands of the ICEHin1056 subfamily have a longstanding relationship with H. influenzae and H. parainfluenzae and are co-evolving as semi-autonomous genomes within the 'supragenomes' of their host species. They have promoted bacterial diversity and adaptation through becoming efficient vectors of antibiotic resistance by the recent acquisition of antibiotic resistance transposons.     

Comparative genomic analysis of Acinetobacter baumannii clinical isolates reveals extensive genomic variation and diverse antibiotic resistance determinants

Background

Acinetobacter baumannii is an important nosocomial pathogen that poses a serious health threat to immune-compromised patients. Due to its rapid ability to develop multidrug resistance (MDR), A. baumannii has increasingly become a focus of attention worldwide. To better understand the genetic variation and antibiotic resistance mechanisms of this bacterium at the genomic level, we reported high-quality draft genome sequences of 8 clinical isolates with various sequence types and drug susceptibility profiles.

Results

We sequenced 7 MDR and 1 drug-sensitive clinical A. baumannii isolates and performed comparative genomic analysis of these draft genomes with 16 A. baumannii complete genomes from GenBank. We found a high degree of variation in A. baumannii, including single nucleotide polymorphisms (SNPs) and large DNA fragment variations in the AbaR-like resistance island (RI) regions, the prophage and the type VI secretion system (T6SS). In addition, we found several new AbaR-like RI regions with highly variable structures in our MDR strains. Interestingly, we found a novel genomic island (designated as GI_BJ4) in the drug-sensitive strain BJ4 carrying metal resistance genes instead of antibiotic resistance genes inserted into the position where AbaR-like RIs commonly reside in other A. baumannii strains. Furthermore, we showed that diverse antibiotic resistance determinants are present outside the RIs in A. baumannii, including antibiotic resistance-gene bearing integrons, the bla_OXA-23-containing transposon Tn2009, and chromosomal intrinsic antibiotic resistance genes.

Conclusions

Our comparative genomic analysis revealed that extensive genomic variation exists in the A. baumannii genome. Transposons, genomic islands and point mutations are the main contributors to the plasticity of the A. baumannii genome and play critical roles in facilitating the development of antibiotic resistance in the clinical isolates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1163) contains supplementary material, which is available to authorized users.     

Identification and comparative analysis of a genomic island in Mycobacterium avium subsp. hominissuis

Mycobacterium avium subsp. hominissuis (MAH) is an environmental bacterium causing opportunistic infections. The objective of this study was to identify flexible genome regions in MAH isolated from different sources. By comparing five complete and draft MAH genomes we identified a genomic island conferring additional flexibility to the MAH genomes. The island was absent in one of the five strains and had sizes between 16.37 and 84.85 kb in the four other strains. The genes present in the islands differed among strains and included phage- and plasmid-derived genes, integrase genes, hypothetical genes, and virulence-associated genes like mmpL or mce genes.     

Bacterial classification using genomic fingerprints obtained by virtual hybridization

In silico genomic fingerprints were produced by virtual hybridization of 191 fully sequenced bacterial genomes using a set of 15,264 13-mer probes specially designed to produce universal whole genome fingerprints. A novel approach for constructing phylogenetic trees, based on comparative analysis of genomic fingerprints, was developed. The resultant bacterial phylogenetic tree had strong similarities to those produced from the alignment of conserved sequences. Notably, the trees derived from the alignment of other conserved COG genes divided the Bacillus and Corynebacterium genera into the same subgroups produced by the novel bacterial tree. A number of discrepancies between both techniques were observed for the grouping of some Lactobacillus species. However, a detailed analysis of the alignment of these genomes using other bioinformatics tools revealed that the grouping of these organisms in the novel tree was more satisfactory than the groupings from previous classifications, which used only a few conserved genes. All these data suggest that the bacterial taxonomy produced by genomic fingerprints is satisfactory, but sometimes different from classical taxonomies. Discrepancies probably arise because the fingerprinting technique analyzes genomic sequences and reveals more information than previously used approaches.     

Bacterial classification using genomic fingerprints obtained by virtual hybridization

In silico genomic fingerprints were produced by virtual hybridization of 191 fully sequenced bacterial genomes using a set of 15,264 13-mer probes specially designed to produce universal whole genome fingerprints. A novel approach for constructing phylogenetic trees, based on comparative analysis of genomic fingerprints, was developed. The resultant bacterial phylogenetic tree had strong similarities to those produced from the alignment of conserved sequences. Notably, the trees derived from the alignment of other conserved COG genes divided the Bacillus and Corynebacterium genera into the same subgroups produced by the novel bacterial tree. A number of discrepancies between both techniques were observed for the grouping of some Lactobacillus species. However, a detailed analysis of the alignment of these genomes using other bioinformatics tools revealed that the grouping of these organisms in the novel tree was more satisfactory than the groupings from previous classifications, which used only a few conserved genes. All these data suggest that the bacterial taxonomy produced by genomic fingerprints is satisfactory, but sometimes different from classical taxonomies. Discrepancies probably arise because the fingerprinting technique analyzes genomic sequences and reveals more information than previously used approaches.     

Strain-specific genes of Helicobacter pylori: genome evolution driven by a novel type IV secretion system and genomic island transfer

The availability of multiple bacterial genome sequences has revealed a surprising extent of variability among strains of the same species. The human gastric pathogen Helicobacter pylori is known as one of the most genetically diverse species. We have compared the genome sequence of the duodenal ulcer strain P12 and six other H. pylori genomes to elucidate the genetic repertoire and genome evolution mechanisms of this species. In agreement with previous findings, we estimate that the core genome comprises about 1200 genes and that H. pylori possesses an open pan-genome. Strain-specific genes are preferentially located at potential genome rearrangement sites or in distinct plasticity zones, suggesting two different mechanisms of genome evolution. The P12 genome contains three plasticity zones, two of which encode type IV secretion systems and have typical features of genomic islands. We demonstrate for the first time that one of these islands is capable of self-excision and horizontal transfer by a conjugative process. We also show that excision is mediated by a protein of the XerD family of tyrosine recombinases. Thus, in addition to its natural transformation competence, conjugative transfer of genomic islands has to be considered as an important source of genetic diversity in H. pylori.     

Differential regulation of toxoflavin production and its role in the enhanced virulence of Burkholderia gladioli

Burkholderia gladioli is a causal agent of bacterial panicle blight and sheath/grain browning in rice in many countries. Many strains produce the yellow pigment toxoflavin, which is highly toxic to plants, fungi, animals and microorganisms. Although there have been several studies on the toxoflavin biosynthesis system of B. glumae, it is still unclear how B. gladioli activates toxoflavin biosynthesis. In this study, we explored the genomic organization of the toxoflavin system of B. gladioli and its biological functions using comparative genomic analysis between toxoflavin‐producing strains (B. glumae BGR1 and B. gladioli BSR3) and a strain not producing toxoflavin (B. gladioli KACC11889). The latter exhibits normal physiological characteristics similar to other B. gladioli strains. Burkholderia gladioli KACC11889 possesses all the genes involved in toxoflavin biosynthesis, but lacks the quorum‐sensing (QS) system that functions as an on/off switch for toxoflavin biosynthesis. These data suggest that B. gladioli has evolved to use the QS signalling cascade of toxoflavin production (TofI/TofR of QS → ToxJ or ToxR → tox operons) similar to that in B. glumae. However, some strains may have evolved to eliminate toxoflavin production through deletion of the QS genes. In addition, we demonstrate that the toxoflavin biosynthetic system enhances the virulence of B. gladioli. These findings provide another line of evidence supporting the differential regulation of the toxoflavin system in Burkholderia strains.     

Dissection of Quorum-Sensing Genes in Burkholderia glumae Reveals Non-Canonical Regulation and the New Regulatory Gene tofM for Toxoflavin Production

Burkholderia glumae causes bacterial panicle blight of rice and produces major virulence factors, including toxoflavin, under the control of the quorum-sensing (QS) system mediated by the luxI homolog, tofI, and the luxR homolog, tofR. In this study, a series of markerless deletion mutants of B. glumae for tofI and tofR were generated using the suicide vector system, pKKSacB, for comprehensive characterization of the QS system of this pathogen. Consistent with the previous studies by other research groups, ΔtofI and ΔtofR strains of B. glumae did not produce toxoflavin in Luria-Bertani (LB) broth. However, these mutants produced high levels of toxoflavin when grown in a highly dense bacterial inoculum (∼ 10¹¹ CFU/ml) on solid media, including LB agar and King’s B (KB) agar media. The ΔtofI/ΔtofR strain of B. glumae, LSUPB201, also produced toxoflavin on LB agar medium. These results indicate the presence of previously unknown regulatory pathways for the production of toxoflavin that are independent of tofI and/or tofR. Notably, the conserved open reading frame (locus tag: bglu_2g14480) located in the intergenic region between tofI and tofR was found to be essential for the production of toxoflavin by tofI and tofR mutants on solid media. This novel regulatory factor of B. glumae was named tofM after its homolog, rsaM, which was recently identified as a novel negative regulatory gene for the QS system of another rice pathogenic bacterium, Pseudomonas fuscovaginae. The ΔtofM strain of B. glumae, LSUPB286, produced a less amount of toxoflavin and showed attenuated virulence when compared with its wild type parental strain, 336gr-1, suggesting that tofM plays a positive role in toxoflavin production and virulence. In addition, the observed growth defect of the ΔtofI strain, LSUPB145, was restored by 1 µM N-octanoyl homoserine lactone (C8-HSL).     

Insights on genomic diversity of Vibrio spp. through Pan-genome analysis

Purpose

The aquaculture sector is a major contributor to the economic and nutritional security for a number of countries. India’s total seafood exports for the year 2017–2018 accounted for US$ Million 7082. One of the major setbacks in this sector is the frequent outbreaks of diseases often due to bacterial pathogens. Vibriosis is one of the major diseases caused by bacteria of Vibrio spp., causing significant economic loss to the aquaculture sector. The objective of this study was to understand the genetic composition of Vibrio spp.

Methods

Thirty-five complete genomes were downloaded from GenBank comprising seven vibrio species, namely, Vibrio alginolyticus, V. anguillarum, V. campbellii, V. harveyi, V. furnissii, V. parahaemolyticus, and V. vulnificus. Pan-genome analysis was carried out with coding sequences (CDS) generated from all the Vibrio genomes. In addition, genomes were mined for genes coding for toxin-antitoxin systems, antibiotic resistance, genomic islands, and virulence factors.

Results

Results revealed an open pan-genome comprising of 2004 core, 8249 accessory, and 6780 unique genes. Downstream analysis of genomes and the identified unique genes resulted in 312 antibiotic resistance genes, 430 genes coding for toxin and antitoxin systems along with 4802, and 4825 putative virulent genes from genomic island regions and unique gene sets, respectively.

Conclusion

Pan-genome and other downstream analytical procedures followed in this study have the potential to predict strain-specific genes and their association with habitat and pathogenicity.