Genomics of epidemic pathogens

Immunometabolism in Arthropod Vectors: Redefining Interspecies Relationships

Metabolism influences biochemical networks, and arthropod vectors are endowed with an immune system that affects microbial acquisition, persistence, and transmission to humans and other animals. Here, we aim to persuade the scientific community to expand their interests in immunometabolism beyond mammalian hosts and towards arthropod vectors. Immunometabolism investigates the interplay of metabolism and immunology. We provide a conceptual framework for investigators from diverse disciplines and indicate that relationships between microbes, mammalian hosts and their hematophagous arthropods may result in cost-effective (mutualism) or energetically expensive (parasitism) interactions. We argue that disparate resource allocations between species may partially explain why some microbes act as pathogens when infecting humans and behave as mutualistic or commensal organisms when colonizing arthropod vectors.

In Silico Identification of Three Types of Integrative and Conjugative Elements in Elizabethkingia anophelis Strains Isolated from around the World

Elizabethkingia anophelis is an opportunistic human pathogen, and the genetic diversity between strains from around the world becomes apparent as more genomes are sequenced. Genome comparison identified three types of putative ICEs in 31 of 36 strains. The diversity of ICEs suggests that they had different origins. One of the ICEs was discovered previously from a large E. anophelis outbreak in Wisconsin in the United States; this ICE has integrated into the mutY gene of the outbreak strain, creating a mutator phenotype. Similar to ICEs found in many bacterial species, ICEs in E. anophelis carry various cargo genes that enable recipients to resist antibiotics and adapt to various ecological niches. The adaptive immune CRISPR-Cas system is present in nine of 36 strains. An ICE-derived spacer was found in the CRISPR locus in a strain that has no ICE, suggesting a past encounter and effective defense against ICE.

Elizabethkingia anophelis is an emerging global multidrug-resistant opportunistic pathogen. We assessed the diversity among 13 complete genomes and 23 draft genomes of E. anophelis strains derived from various environmental settings and human infections from different geographic regions around the world from 1950s to the present. Putative integrative and conjugative elements (ICEs) were identified in 31/36 (86.1%) strains in the study. A total of 52 putative ICEs (including eight degenerated elements lacking integrases) were identified and categorized into three types based on the architecture of the conjugation module and the phylogeny of the relaxase, coupling protein, TraG, and TraJ protein sequences. The type II and III ICEs were found to integrate adjacent to tRNA genes, while type I ICEs integrate into intergenic regions or into a gene. The ICEs carry various cargo genes, including transcription regulator genes and genes conferring antibiotic resistance. The adaptive immune CRISPR-Cas system was found in nine strains, including five strains in which CRISPR-Cas machinery and ICEs coexist at different locations on the same chromosome. One ICE-derived spacer was present in the CRISPR locus in one strain. ICE distribution in the strains showed no geographic or temporal patterns. The ICEs in E. anophelis differ in architecture and sequence from CTnDOT, a well-studied ICE prevalent in Bacteroides spp. The categorization of ICEs will facilitate further investigations of the impact of ICE on virulence, genome epidemiology, and adaptive genomics of E. anophelis.

IMPORTANCEElizabethkingia anophelis is an opportunistic human pathogen, and the genetic diversity between strains from around the world becomes apparent as more genomes are sequenced. Genome comparison identified three types of putative ICEs in 31 of 36 strains. The diversity of ICEs suggests that they had different origins. One of the ICEs was discovered previously from a large E. anophelis outbreak in Wisconsin in the United States; this ICE has integrated into the mutY gene of the outbreak strain, creating a mutator phenotype. Similar to ICEs found in many bacterial species, ICEs in E. anophelis carry various cargo genes that enable recipients to resist antibiotics and adapt to various ecological niches. The adaptive immune CRISPR-Cas system is present in nine of 36 strains. An ICE-derived spacer was found in the CRISPR locus in a strain that has no ICE, suggesting a past encounter and effective defense against ICE.

Zinc-dependent substrate-level phosphorylation powers Salmonella growth under nitrosative stress of the innate host response

The metabolic processes that enable the replication of intracellular Salmonella under nitrosative stress conditions engendered in the innate response of macrophages are poorly understood. A screen of Salmonella transposon mutants identified the ABC-type high-affinity zinc uptake system ZnuABC as a critical determinant of the adaptation of Salmonella to the nitrosative stress generated by the enzymatic activity of inducible nitric oxide (NO) synthase of mononuclear phagocytic cells. NO limits the virulence of a znuB mutant in an acute murine model of salmonellosis. The ZnuABC transporter is crucial for the glycolytic function of fructose bisphosphate aldolase, thereby fueling growth of Salmonella during nitrosative stress produced in the innate response of macrophages. Our investigations demonstrate that glycolysis mediates resistance of Salmonella to the antimicrobial activity of NO produced in an acute model of infection. The ATP synthesized by substrate-level phosphorylation at the payoff phase of glycolysis and acetate fermentation powers the replication of Salmonella experiencing high levels of nitrosative stress. In contrast, despite its high potential for ATP synthesis, oxidative phosphorylation is a major target of inhibition by NO and contributes little to the antinitrosative defenses of intracellular Salmonella. Our investigations have uncovered a previously unsuspected conjunction between zinc homeostasis, glucose metabolism and cellular energetics in the adaptation of intracellular Salmonella to the reactive nitrogen species synthesized in the innate host response.

Microbial pathogens are exposed to multiple antimicrobial defenses during their associations with host cells. Nitric oxide generated in the innate response exerts widespread antimicrobial activity against a variety of pathogenic microorganisms. Nitric oxide has high affinity for metal groups of terminal cytochromes of the respiratory chain, and thus nitrosative stress exerts extreme deleterious actions against the cellular energetics that rely on oxidative phosphorylation. Intracellular Salmonella have resolved this dilemma by satisfying a significant portion of their energetic demands via substrate level phosphorylation in the payoff phase of glycolysis and acetate fermentation. A high affinity zinc uptake system promotes antinitrosative defense of intracellular Salmonella by in great part supporting the enzymatic activity of an essential enzyme in the preparatory phase of glycolysis. Our research provides novel insights into the metabolic and energetic adaptations that allow a bacterial pathogen to thrive in the midst of the innate host response of vertebrate cells.

Identification of constraints influencing the bacterial genomes evolution in the PVC super-phylum

BACKGROUND

Horizontal transfer plays an important role in the evolution of bacterial genomes, yet it obeys several constraints, including the ecological opportunity to meet other organisms, the presence of transfer systems, and the fitness of the transferred genes. Bacteria from the Planctomyctetes, Verrumicrobia, Chlamydiae (PVC) super-phylum have a compartmentalized cell plan delimited by an intracytoplasmic membrane that might constitute an additional constraint with particular impact on bacterial evolution. In this investigation, we studied the evolution of 33 genomes from PVC species and focused on the rate and the nature of horizontally transferred sequences in relation to their habitat and their cell plan.

RESULTS

Using a comparative phylogenomic approach, we showed that habitat influences the evolution of the bacterial genome's content and the flux of horizontal transfer of DNA (HT). Thus bacteria from soil, from insects and ubiquitous bacteria presented the highest average of horizontal transfer compared to bacteria living in water, extracellular bacteria in vertebrates, bacteria from amoeba and intracellular bacteria in vertebrates (with a mean of 379 versus 110 events per species, respectively and 7.6% of each genomes due to HT against 4.8%). The partners of these transfers were mainly bacterial organisms (94.9%); they allowed us to differentiate environmental bacteria, which exchanged more with Proteobacteria, and bacteria from vertebrates, which exchanged more with Firmicutes. The functional analysis of the horizontal transfers revealed a convergent evolution, with an over-representation of genes encoding for membrane biogenesis and lipid metabolism, among compartmentalized bacteria in the different habitats.

CONCLUSIONS

The presence of an intracytoplasmic membrane in PVC species seems to affect the genome's evolution through the selection of transferred DNA, according to their encoded functions.

Compartmentalization in PVC super-phylum: evolution and impact

BACKGROUND

The PVC super-phylum gathers bacteria from seven phyla (Planctomycetes, Verrucomicrobiae, Chlamydiae, Lentisphaera, Poribacteria, OP3, WWE2) presenting different lifestyles, cell plans and environments. Planctomyces and several Verrucomicrobiae exhibit a complex cell plan, with an intracytoplasmic membrane inducing the compartmentalization of the cytoplasm into two regions (pirellulosome and paryphoplasm). The evolution and function of this cell plan is still subject to debate. In this work, we hypothesized that it could play a role in protection of the bacterial DNA, especially against Horizontal Genes Transfers (HGT). Therefore, 64 bacterial genomes belonging to seven different phyla (whose four PVC phyla) were studied. We reconstructed the evolution of the cell plan as precisely as possible, thanks to information obtained by bibliographic study and electronic microscopy. We used a strategy based on comparative phylogenomic in order to determine the part occupied by the horizontal transfers for each studied genomes.

RESULTS

Our results show that the bacteria Simkania negevensis (Chlamydiae) and Coraliomargarita akajimensis (Verrucomicrobiae), whose cell plan were unknown before, are compartmentalized, as we can see on the micrographies. This is one of the first indication of the presence of an intracytoplasmic membrane in a Chlamydiae. The proportion of HGT does not seems to be related to the cell plan of bacteria, suggesting that compartmentalization does not induce a protection of bacterial DNA against HGT. Conversely, lifestyle of bacteria seems to impact the ability of bacteria to exchange genes.

CONCLUSIONS

Our study allows a best reconstruction of the evolution of intracytoplasmic membrane, but this structure seems to have no impact on HGT occurrences.

REVIEWERS

This article was reviewed by Mircea Podar and Olivier Tenaillon.

Emerging Roles of Toxin-Antitoxin Modules in Bacterial Pathogenesis

Toxin-antitoxin (TA) cassettes are encoded widely by bacteria. The modules typically comprise a protein toxin and protein or RNA antitoxin that sequesters the toxin factor. Toxin activation in response to environmental cues or other stresses promotes a dampening of metabolism, most notably protein translation, which permits survival until conditions improve. Emerging evidence also implicates TAs in bacterial pathogenicity. Bacterial persistence involves entry into a transient semi-dormant state in which cells survive unfavorable conditions including killing by antibiotics, which is a significant clinical problem. TA complexes play a fundamental role in inducing persistence by downregulating cellular metabolism. Bacterial biofilms are important in numerous chronic inflammatory and infectious diseases and cause serious therapeutic problems due to their multidrug tolerance and resistance to host immune system actions. Multiple TAs influence biofilm formation through a network of interactions with other factors that mediate biofilm production and maintenance. Moreover, in view of their emerging contributions to bacterial virulence, TAs are potential targets for novel prophylactic and therapeutic approaches that are required urgently in an era of expanding antibiotic resistance. This review summarizes the emerging evidence that implicates TAs in the virulence profiles of a diverse range of key bacterial pathogens that trigger serious human disease.

The contribution of genomics to the study of Q fever

Coxiella burnetii is the etiological agent of Q fever, a worldwide zoonosis that can result in large outbreaks. The birth of genomics and sequencing of C. burnetii strains has revolutionized many fields of study of this infection. Accurate genotyping methods and comparative genomic analysis have enabled description of the diversity of strains around the world and their link with pathogenicity. Genomics has also permitted the development of qPCR tools and axenic culture medium, facilitating the diagnosis of Q fever. Moreover, several pathophysiological mechanisms can now be predicted and therapeutic strategies can be determined thanks to in silico genome analysis. An extensive pan-genomic analysis will allow for a comprehensive view of the clonal diversity of C. burnetii and its link with virulence.

One cannot rule them all: Are bacterial toxins-antitoxins druggable?

Type II (proteic) toxin–antitoxin (TA) operons are widely spread in bacteria and archaea. They are organized as operons in which, usually, the antitoxin gene precedes the cognate toxin gene. The antitoxin generally acts as a transcriptional self-repressor, whereas the toxin acts as a co-repressor, both proteins constituting a harmless complex. When bacteria encounter a stressful environment, TAs are triggered. The antitoxin protein is unstable and will be degraded by host proteases, releasing the free toxin to halt essential processes. The result is a cessation of cell growth or even death. Because of their ubiquity and the essential processes targeted, TAs have been proposed as good candidates for development of novel antimicrobials. We discuss here the possible druggability of TAs as antivirals and antibacterials, with focus on the potentials and the challenges that their use may find in the ‘real’ world. We present strategies to develop TAs as antibacterials in view of novel technologies, such as the use of very small molecules (fragments) as inhibitors of protein–protein interactions. Appropriate fragments could disrupt the T:A interfaces leading to the release of the targeted TA pair. Possible ways of delivery and formulation of Tas are also discussed.

We consider various approaches to develop the toxins of the type II family as possible candidates to drug discovery; druggability of toxins-antitoxins could be possible as antivirals. As antibacterials, they might be considered as druggable but delivery and formulation may not be simple so far.

Recent molecular insights into rickettsial pathogenesis and immunity

Human infections with arthropod-borne Rickettsia species remain a major global health issue, causing significant morbidity and mortality. Epidemic typhus due to Rickettsia prowazekii has an established reputation as the 'scourge of armies', and as a major determinant of significant 'historical turning points'. No suitable vaccines for human use are currently available to prevent rickettsial diseases. The unique lifestyle features of rickettsiae include obligate intracellular parasitism, intracytoplasmic niche within the host cell, predilection for infection of microvascular endothelium in mammalian hosts, association with arthropods and the tendency for genomic reduction. The fundamental research in the field of Rickettsiology has witnessed significant recent progress in the areas of pathogen adhesion/invasion and host immune responses, as well as the genomics, proteomics, metabolomics, phylogenetics, motility and molecular manipulation of important rickettsial pathogens. The focus of this review article is to capture a snapshot of the latest developments pertaining to the mechanisms of rickettsial pathogenesis and immunity.

Whole genome sequencing of bacteria in cystic fibrosis as a model for bacterial genome adaptation and evolution

Cystic fibrosis (CF) airways harbor a wide variety of new and/or emerging multidrug resistant bacteria which impose a heavy burden on patients. These bacteria live in close proximity with one another, which increases the frequency of lateral gene transfer. The exchange and movement of mobile genetic elements and genomic islands facilitate the spread of genes between genetically diverse bacteria, which seem to be advantageous to the bacterium as it allows adaptation to the new niches of the CF lungs. Niche adaptation is one of the major evolutionary forces shaping bacterial genome composition and in CF the chronic strains adapt and become less virulent. The purpose of this review is to shed light on CF bacterial genome alterations. Next-generation sequencing technology is an exciting tool that may help us to decipher the genome architecture and the evolution of bacteria colonizing CF lungs.

LTQ-XL mass spectrometry proteome analysis expands the Pseudomonas aeruginosa AmpR regulon to include cyclic di-GMP phosphodiesterases and phosphoproteins, and identifies novel open reading frames

Pseudomonas aeruginosa is well known for its antibiotic resistance and intricate regulatory network, contributing to its success as an opportunistic pathogen. This study is an extension of our transcriptomic analyses (microarray and RNA-Seq) to understand the global changes in PAO1 upon deleting a gene encoding a transcriptional regulator AmpR, in the presence and absence of β-lactam antibiotic. This study was performed under identical conditions to explore the proteome profile of the ampR deletion mutant (PAOΔampR) using LTQ-XL mass spectrometry. The proteomic data identified ~53% of total PAO1 proteins and expanded the master regulatory role of AmpR in determining antibiotic resistance and multiple virulence phenotypes in P. aeruginosa. AmpR proteome analysis identified 853 AmpR-dependent proteins, which include 102 transcriptional regulators and 21 two-component system proteins. AmpR also regulates cyclic di-GMP phosphodiesterases (PA4367, PA4969, PA4781) possibly affecting major virulence systems. Phosphoproteome analysis also suggests a significant role for AmpR in Ser, Thr and Tyr phosphorylation. These novel mechanisms of gene regulation were previously not associated with AmpR. The proteome analysis also identified many unannotated and misannotated ORFs in the P. aeruginosa genome. Thus, our data sheds light on important virulence regulatory pathways that can potentially be exploited to deal with P. aeruginosa infections.

BIOLOGICAL SIGNIFICANCE

The AmpR proteome data not only confirmed the role of AmpR in virulence and resistance to multiple antibiotics, but also expanded the perimeter of AmpR regulon. The data presented here points to the role of AmpR in regulating cyclic di-GMP levels and phosphorylation of Ser, Thr and Tyr, adding another dimension to the regulatory functions of AmpR. We also identify some previously unannotated/misannotated ORFs in the P. aeruginosa genome, indicating the limitations of existing ORF analyses software. This study will contribute towards understanding complex genetic organization of P. aeruginosa. Whole genome proteomic picture of regulators at higher nodal positions in the regulatory network will not only help us link various virulence phenotypes but also design novel therapeutic strategies.

Characterization and comparison of 2 distinct epidemic community-associated methicillin-resistant Staphylococcus aureus clones of ST59 lineage

Sequence type (ST) 59 is an epidemic lineage of community-associated (CA) methicillin-resistant Staphylococcus aureus (MRSA) isolates. Taiwanese CA-MRSA isolates belong to ST59 and can be grouped into 2 distinct clones, a virulent Taiwan clone and a commensal Asian-Pacific clone. The Taiwan clone carries the Panton–Valentine leukocidin (PVL) genes and the staphylococcal chromosomal cassette mec (SCCmec) V_T, and is frequently isolated from patients with severe disease. The Asian-Pacific clone is PVL-negative, carries SCCmec IV, and a frequent colonizer of healthy children. Isolates of both clones were characterized by their ability to adhere to respiratory A549 cells, cytotoxicity to human neutrophils, and nasal colonization of a murine and murine sepsis models. Genome variation was determined by polymerase chain reaction of selected virulence factors and by multi-strain whole genome microarray. Additionally, the expression of selected factors was compared between the 2 clones. The Taiwan clone showed a much higher cytotoxicity to the human neutrophils and caused more severe septic infections with a high mortality rate in the murine model. The clones were indistinguishable in their adhesion to A549 cells and persistence of murine nasal colonization. The microarray data revealed that the Taiwan clone had lost the ø3-prophage that integrates into the β-hemolysin gene and includes staphylokinase- and enterotoxin P-encoding genes, but had retained the genes for human immune evasion, scn and chps. Production of the virulence factors did not differ significantly in the 2 clonal groups, although more α-toxin was expressed in Taiwan clone isolates from pneumonia patients. In conclusion, the Taiwan CA-MRSA clone was distinguished by enhanced virulence in both humans and an animal infection model. The evolutionary acquisition of PVL, the higher expression of α-toxin, and possibly the loss of a large portion of the β-hemolysin-converting prophage likely contribute to its higher pathogenic potential than the Asian-Pacific clone.

Investigation of antibiotic resistance in the genomic era of multidrug-resistant Gram-negative bacilli, especially Enterobacteriaceae, Pseudomonas and Acinetobacter

The increase and spread of multidrug-resistant (MDR) Gram-negative bacteria, including Enterobacteriaceae, Pseudomonas and Acinetobacter species, have become major concerns worldwide. Although the frequent misuse of antibiotic drugs has greatly contributed to worldwide antibiotic resistance by causing a large dispersal of resistance determinants, recent studies demonstrate that these resistance determinants could have emerged from ancient or environmental sources. Moreover, during the last 10 years, we have been witnessing the emergence and development of technologies for high-throughput sequencing, coinciding with an exponential increase in the number of bacterial genomes sequenced. These sequencing technologies allow a complete study of MDR bacterial genomes and are the best way to investigate the genetic determinants of antimicrobial resistance. Accordingly, studies using genome sequencing to decipher resistance determinants in Enterobacteriaceae, Pseudomonas and Acinetobacter species have demonstrated several advantages including, among others: an exhaustive identification of resistance determinants from any clinical, epidemiological or environmental MDR bacterium; identification of the acquisition mechanisms for resistance determinants exchanged between bacterial species through mobile genetic elements and elucidation and understanding, in record time (less than 1 week), of some critical or epidemic events caused by particular pathogenic bacteria. Therefore, it is clear today that the bacterial genome sequencing approach has revolutionized all fields of scientific research and represents a powerful tool to explore the world of microorganisms.

Integrated whole-genome sequencing and temporospatial analysis of a continuing Group A Streptococcus epidemic

Analysis of microbial epidemics has been revolutionized by whole-genome sequencing. We recently sequenced the genomes of 601 type emm59 Group A Streptococcus (GAS) organisms responsible for an ongoing epidemic of invasive infections in Canada and some of the United States. The epidemic has been caused by the emergence of a genetically distinct, hypervirulent clone that has genetically diversified. The ease of obtaining genomic data contrasts with the relatively difficult task of translating them into insightful epidemiological information. Here, we sequenced the genomes of 90 additional invasive Canadian emm59 GAS organisms, including 80 isolated recently in 2010–2011. We used an improved bioinformatics pipeline designed to rapidly process and analyze whole-genome data and integrate strain metadata. We discovered that emm59 GAS organisms are undergoing continued multiclonal evolutionary expansion. Previously identified geographic patterns of strain dissemination are being diluted as mixing of subclones over time and space occurs. Our integrated data analysis strategy permits prompt and accurate mapping of the dissemination of bacterial organisms in an epidemic wave, permitting rapid generation of hypotheses that inform public health and virulence studies.

A comparison of the endotoxin biosynthesis and protein oxidation pathways in the biogenesis of the outer membrane of Escherichia coli and Neisseria meningitidis

The Gram-negative bacterial cell envelope consists of an inner membrane (IM) that surrounds the cytoplasm and an asymmetrical outer-membrane (OM) that forms a protective barrier to the external environment. The OM consists of lipopolysaccahride (LPS), phospholipids, outer membrane proteins (OMPs), and lipoproteins. Oxidative protein folding mediated by periplasmic oxidoreductases is required for the biogenesis of the protein components, mainly constituents of virulence determinants such as pili, flagella, and toxins, of the Gram-negative OM. Recently, periplasmic oxidoreductases have been implicated in LPS biogenesis of Escherichia coli and Neisseria meningitidis. Differences in OM biogenesis, in particular the transport pathways for endotoxin to the OM, the composition and role of the protein oxidation, and isomerization pathways and the regulatory networks that control them have been found in these two Gram-negative species suggesting that although form and function of the OM is conserved, the pathways required for the biosynthesis of the OM and the regulatory circuits that control them have evolved to suit the lifestyle of each organism.

The rhizome of the multidrug-resistant Enterobacter aerogenes genome reveals how new "killer bugs" are created because of a sympatric lifestyle

Here, we sequenced the 5,419,609 bp circular genome of an Enterobacter aerogenes clinical isolate that killed a patient and was resistant to almost all current antibiotics (except gentamicin) commonly used to treat Enterobacterial infections, including colistin. Genomic and phylogenetic analyses explain the discrepancies of this bacterium and show that its core genome originates from another genus, Klebsiella. Atypical characteristics of this bacterium (i.e., motility, presence of ornithine decarboxylase, and lack of urease activity) are attributed to genomic mosaicism, by acquisition of additional genes, such as the complete 60,582 bp flagellar assembly operon acquired "en bloc" from the genus Serratia. The genealogic tree of the 162,202 bp multidrug-resistant conjugative plasmid shows that it is a chimera of transposons and integrative conjugative elements from various bacterial origins, resembling a rhizome. Moreover, we demonstrate biologically that a G53S mutation in the pmrA gene results in colistin resistance. E. aerogenes has a large RNA population comprising 8 rRNA operons and 87 cognate tRNAs that have the ability to translate transferred genes that use different codons, as exemplified by the significantly different codon usage between genes from the core genome and the "mobilome." On the basis of our findings, the evolution of this bacterium to become a "killer bug" with new genomic repertoires was from three criteria that are "opportunity, power, and usage" to indicate a sympatric lifestyle: "opportunity" to meet other bacteria and exchange foreign sequences since this bacteria was similar to sympatric bacteria; "power" to integrate these foreign sequences such as the acquisition of several mobile genetic elements (plasmids, integrative conjugative element, prophages, transposons, flagellar assembly system, etc.) found in his genome; and "usage" to have the ability to translate these sequences including those from rare codons to serve as a translator of foreign languages.

Unique core genomes of the bacterial family vibrionaceae: insights into niche adaptation and speciation

Background

The criteria for defining bacterial species and even the concept of bacterial species itself are under debate, and the discussion is apparently intensifying as more genome sequence data is becoming available. However, it is still unclear how the new advances in genomics should be used most efficiently to address this question. In this study we identify genes that are common to any group of genomes in our dataset, to determine whether genes specific to a particular taxon exist and to investigate their potential role in adaptation of bacteria to their specific niche. These genes were named unique core genes. Additionally, we investigate the existence and importance of unique core genes that are found in isolates of phylogenetically non-coherent groups. These groups of isolates, that share a genetic feature without sharing a closest common ancestor, are termed genophyletic groups.

Results

The bacterial family Vibrionaceae was used as the model, and we compiled and compared genome sequences of 64 different isolates. Using the software orthoMCL we determined clusters of homologous genes among the investigated genome sequences. We used multilocus sequence analysis to build a host phylogeny and mapped the numbers of unique core genes of all distinct groups of isolates onto the tree. The results show that unique core genes are more likely to be found in monophyletic groups of isolates. Genophyletic groups of isolates, in contrast, are less common especially for large groups of isolate. The subsequent annotation of unique core genes that are present in genophyletic groups indicate a high degree of horizontally transferred genes. Finally, the annotation of the unique core genes of Vibrio cholerae revealed genes involved in aerotaxis and biosynthesis of the iron-chelator vibriobactin.

Conclusion

The presented work indicates that genes specific for any taxon inside the bacterial family Vibrionaceae exist. These unique core genes encode conserved metabolic functions that can shed light on the adaptation of a species to its ecological niche. Additionally, our study suggests that unique core genes can be used to aid classification of bacteria and contribute to a bacterial species definition on a genomic level. Furthermore, these genes may be of importance in clinical diagnostics and drug development.