Characterization and inference of gene gain/loss along burkholderia evolutionary history

Analysis of multipartite bacterial genomes using alignment free and alignment-based pipelines

Since the discovery of second chromosome in Rhodobacter sphaeroides 2.4.1 in 1989, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This has shattered the unipartite (single chromosome) genome dogma in bacteria. Since then, many questions on various aspects of multipartite genomes in bacteria have been addressed. However, our understanding of how multipartite genomes emerge and evolve is still lacking. Importantly, the knowledge of genetic factors underlying the differences in multipartite and single-chromosome genomes is lacking. In this work, we have performed comparative evolutionary and functional genomics analyses to identify molecular factors that discriminate multipartite from unipartite bacteria, with the goal to decipher taxon-specific factors, and those that are prevalent across the taxa, underlying these traits. We assessed the roles of evolutionary mechanisms, specifically gene gain, in driving the divergence of bacteria with single and multiple chromosomes. In addition, we performed functional genomic analysis to garner support for our findings from comparative evolutionary analysis. We found genes such as those encoding conserved hypothetical proteins in Deinococcus radiodurans R1, and putative phage phi-C31 gp36 major capsid like and hypothetical proteins in Rhodobacter sphaeroides 2.4.1, which are located on accessory chromosomes in these bacteria but were not found in the inferred ancestral sequences, and on the primary chromosomes, as well as were not found in their closest relatives with single chromosome within the same clade. Our study shines a new light on the potential roles of the secondary chromosomes in helping bacteria with multipartite genomes to adapt to specialized environments or growth conditions.

Tal1NXtc01 in Xanthomonas translucens pv. cerealis Contributes to Virulence in Bacterial Leaf Streak of Wheat

Xanthomonas translucens pv. cerealis (Xtc) causes bacterial leaf streak (BLS) of important cereal crops, including wheat (Triticum aestivum) and barley (Hordeum vulgare). Transcription activator-like effectors (TALEs) play vital roles in many plant diseases caused by Xanthomonas spp., however, TALEs have not been previously characterized in Xtc. In this study, the whole genome of NXtc01, a virulent strain of Xtc from Xinjiang, China, was sequenced and compared with genomes of other Xanthomonas spp. Xtc NXtc01 consists of a single 4,622,298 bp chromosome that encodes 4,004 genes. Alignment of the NXtc01 sequence with the draft genome of Xtc strain CFBP 2541 (United States) revealed a single giant inversion and differences in the location of two tal genes, which were designated tal1 and tal2. In NXtc01, both tal genes are located on the chromosome, whereas tal2 is plasmid-encoded in CFBP 2541. The repeat variable diresidues (RVDs) at the 12th and 13th sites within Tal2 repeat units were identical in both strains, whereas Tal1 showed differences in the third RVD. Xtc NXtc01 and CFBP 2541 encoded 35 and 33 non-TALE type III effectors (T3Es), respectively. tal1, tal2, and tal-free deletion mutants of Xtc NXtc01 were constructed and evaluated for virulence. The tal1 and tal-free deletion mutants were impaired with respect to symptom development and growth in wheat, suggesting that tal1 is a virulence factor in NXtc01. This was confirmed in gain-of-function experiments that showed the introduction of tal1, but not tal2, restored virulence to the tal-free mutant. Furthermore, we generated a hrcC deletion mutant of NXtc01; the hrcC mutant was non-pathogenic on wheat and unable to elicit a hypersensitive response in the non-host Nicotiana benthamiana. Our data provide a platform for exploring the roles of both TALEs and non-TALEs in promoting BLS on wheat.

Complete Genome Sequence of Massilia oculi sp. nov. CCUG 43427T (=DSM 26321T), the Type Strain of M. oculi, and Comparison with Genome Sequences of Other Massilia Strains

Massilia oculi sp. nov. of type strain CCUG 43427^T is a Gram-negative, rod-shaped, nonspore-forming bacterium, which was recently isolated from the eye of a patient suffering from endophthalmitis and was described as novel species in Massilia genus. In this study, we present the complete genome sequence of this strain by using Pacbio SMRT cell platform and compare this sequence with the genomes of 30 Massilia representative strains. Also, a comprehensive search was conducted for genes and proteins involved in antibiotic resistance and pathogenicity. The genome of CCUG 43427^T is 5,844,653 bp with 65.55% GC content. This genome contains four prophages and four genomic islands (GIs). The cobalt/zinc/cadmium transporter locus CzcABCD is included in these GIs. This GI was predicted to play important role in bacterial heavy-metal tolerance. The in silico genome analysis also revealed that this strain contains a lot of antibiotic resistance and pathogenicity related genes. This result suggested that this strain may has evolved a wide arsenal of weapons for pathogenicity and survival. Genome comparison among CCUG 43427^T and other 30 Massilia strains revealed that more than 400 genes are unique in CCUG 43427^T. Among these, one gene cluster, which was annotated to be important for LOS biosynthesis, catalytic mechanism and the substrate specificity of the enzyme, was predicted to be horizontally transferred by using phylogenies and biased GC content.

Molecular signatures and phylogenomic analysis of the genus Burkholderia: proposal for division of this genus into the emended genus Burkholderia containing pathogenic organisms and a new genus Paraburkholderia gen. nov. harboring environmental species

The genus Burkholderia contains large number of diverse species which include many clinically important organisms, phytopathogens, as well as environmental species. However, currently, there is a paucity of biochemical or molecular characteristics which can reliably distinguish different groups of Burkholderia species. We report here the results of detailed phylogenetic and comparative genomic analyses of 45 sequenced species of the genus Burkholderia. In phylogenetic trees based upon concatenated sequences for 21 conserved proteins as well as 16S rRNA gene sequence based trees, members of the genus Burkholderia grouped into two major clades. Within these main clades a number of smaller clades including those corresponding to the clinically important Burkholderia cepacia complex (BCC) and the Burkholderia pseudomallei groups were also clearly distinguished. Our comparative analysis of protein sequences from Burkholderia spp. has identified 42 highly specific molecular markers in the form of conserved sequence indels (CSIs) that are uniquely found in a number of well-defined groups of Burkholderia spp. Six of these CSIs are specific for a group of Burkholderia spp. (referred to as Clade I in this work) which contains all clinically relevant members of the genus (viz. the BCC and the B. pseudomallei group) as well as the phytopathogenic Burkholderia spp. The second main clade (Clade II), which is composed of environmental Burkholderia species, is also distinguished by 2 identified CSIs that are specific for this group. Additionally, our work has also identified multiple CSIs that serve to clearly demarcate a number of smaller groups of Burkholderia spp. including 3 CSIs that are specific for the B. cepacia complex, 4 CSIs that are uniquely found in the B. pseudomallei group, 5 CSIs that are specific for the phytopathogenic Burkholderia spp. and 22 other CSI that distinguish two groups within Clade II. The described molecular markers provide highly specific means for the demarcation of different groups of Burkholderia spp. and they also offer novel and useful targets for the development of diagnostic assays for the clinically important members of the BCC or the pseudomallei groups. Based upon the results of phylogenetic analyses, the identified CSIs and the pathogenicity profile of Burkholderia species, we are proposing a division of the genus Burkholderia into two genera. In this new proposal, the emended genus Burkholderia will correspond to the Clade I and it will contain only the clinically relevant and phytopathogenic Burkholderia species. All other Burkholderia spp., which are primarily environmental, will be transferred to a new genus Paraburkholderia gen. nov.

The complete genome of Burkholderia phenoliruptrix strain BR3459a, a symbiont of Mimosa flocculosa: highlighting the coexistence of symbiotic and pathogenic genes

BACKGROUND

Burkholderia species play an important ecological role related to xenobiosis, the promotion of plant growth, the biocontrol of agricultural diseases, and symbiotic and non-symbiotic biological nitrogen fixation. Here, we highlight our study as providing the first complete genome of a symbiotic strain of B. phenoliruptrix, BR3459a (=CLA1), which was originally isolated in Brazil from nodules of Mimosa flocculosa and is effective in fixing nitrogen in association with this leguminous species.

RESULTS

Genomic comparisons with other pathogenic and non-pathogenic Burkholderia strains grouped B. phenoliruptrix BR3459a with plant-associated beneficial and environmental species, although it shares a high percentage of its gene repertoire with species of the B. cepacia complex (Bcc) and "pseudomallei" group. The genomic analyses showed that the bce genes involved in exopolysaccharide production are clustered together in the same genomic region, constituting part of the Group III cluster of non-pathogenic bacteria. Regarding environmental stresses, we highlight genes that might be relevant in responses to osmotic, heat, cold and general stresses. Furthermore, a number of particularly interesting genes involved in the machinery of the T1SS, T2SS, T3SS, T4ASS and T6SS secretion systems were identified. The xenobiotic properties of strain BR3459a were also investigated, and some enzymes involved in the degradation of styrene, nitrotoluene, dioxin, chlorocyclohexane, chlorobenzene and caprolactam were identified. The genomic analyses also revealed a large number of antibiotic-related genes, the most important of which were correlated with streptomycin and novobiocin. The symbiotic plasmid showed high sequence identity with the symbiotic plasmid of B. phymatum. Additionally, comparative analysis of 545 housekeeping genes among pathogenic and non-pathogenic Burkholderia species strongly supports the definition of a new genus for the second branch, which would include BR3459a.

CONCLUSIONS

The analyses of B. phenoliruptrix BR3459a showed key property of fixing nitrogen that together with genes for high tolerance to environmental stresses might explain a successful strategy of symbiosis in the tropics. The strain also harbours interesting sets of genes with biotechnological potential. However, the resemblance of certain genes to those of pathogenic Burkholderia raise concerns about large-scale applications in agriculture or for bioremediation.

Global patterns of protein domain gain and loss in superkingdoms

Domains are modules within proteins that can fold and function independently and are evolutionarily conserved. Here we compared the usage and distribution of protein domain families in the free-living proteomes of Archaea, Bacteria and Eukarya and reconstructed species phylogenies while tracing the history of domain emergence and loss in proteomes. We show that both gains and losses of domains occurred frequently during proteome evolution. The rate of domain discovery increased approximately linearly in evolutionary time. Remarkably, gains generally outnumbered losses and the gain-to-loss ratios were much higher in akaryotes compared to eukaryotes. Functional annotations of domain families revealed that both Archaea and Bacteria gained and lost metabolic capabilities during the course of evolution while Eukarya acquired a number of diverse molecular functions including those involved in extracellular processes, immunological mechanisms, and cell regulation. Results also highlighted significant contemporary sharing of informational enzymes between Archaea and Eukarya and metabolic enzymes between Bacteria and Eukarya. Finally, the analysis provided useful insights into the evolution of species. The archaeal superkingdom appeared first in evolution by gradual loss of ancestral domains, bacterial lineages were the first to gain superkingdom-specific domains, and eukaryotes (likely) originated when an expanding proto-eukaryotic stem lineage gained organelles through endosymbiosis of already diversified bacterial lineages. The evolutionary dynamics of domain families in proteomes and the increasing number of domain gains is predicted to redefine the persistence strategies of organisms in superkingdoms, influence the make up of molecular functions, and enhance organismal complexity by the generation of new domain architectures. This dynamics highlights ongoing secondary evolutionary adaptations in akaryotic microbes, especially Archaea.

Updated clusters of orthologous genes for Archaea: a complex ancestor of the Archaea and the byways of horizontal gene transfer

Background

Collections of Clusters of Orthologous Genes (COGs) provide indispensable tools for comparative genomic analysis, evolutionary reconstruction and functional annotation of new genomes. Initially, COGs were made for all complete genomes of cellular life forms that were available at the time. However, with the accumulation of thousands of complete genomes, construction of a comprehensive COG set has become extremely computationally demanding and prone to error propagation, necessitating the switch to taxon-specific COG collections. Previously, we reported the collection of COGs for 41 genomes of Archaea (arCOGs). Here we present a major update of the arCOGs and describe evolutionary reconstructions to reveal general trends in the evolution of Archaea.

Results

The updated version of the arCOG database incorporates 91% of the pangenome of 120 archaea (251,032 protein-coding genes altogether) into 10,335 arCOGs. Using this new set of arCOGs, we performed maximum likelihood reconstruction of the genome content of archaeal ancestral forms and gene gain and loss events in archaeal evolution. This reconstruction shows that the last Common Ancestor of the extant Archaea was an organism of greater complexity than most of the extant archaea, probably with over 2,500 protein-coding genes. The subsequent evolution of almost all archaeal lineages was apparently dominated by gene loss resulting in genome streamlining. Overall, in the evolution of Archaea as well as a representative set of bacteria that was similarly analyzed for comparison, gene losses are estimated to outnumber gene gains at least 4 to 1. Analysis of specific patterns of gene gain in Archaea shows that, although some groups, in particular Halobacteria, acquire substantially more genes than others, on the whole, gene exchange between major groups of Archaea appears to be largely random, with no major ‘highways’ of horizontal gene transfer.

Conclusions

The updated collection of arCOGs is expected to become a key resource for comparative genomics, evolutionary reconstruction and functional annotation of new archaeal genomes. Given that, in spite of the major increase in the number of genomes, the conserved core of archaeal genes appears to be stabilizing, the major evolutionary trends revealed here have a chance to stand the test of time.

Reviewers

This article was reviewed by (for complete reviews see the Reviewers’ Reports section): Dr. PLG, Prof. PF, Dr. PL (nominated by Prof. JPG).

An objective approach for Burkholderia pseudomallei strain selection as challenge material for medical countermeasures efficacy testing

Burkholderia pseudomallei is the causative agent of melioidosis, a rare disease of biodefense concern with high mortality and extreme difficulty in treatment. No human vaccines are available that protect against B. pseudomallei infection, and with the current limitations of antibiotic treatment, the development of new preventative and therapeutic interventions is crucial. Although clinical trials could be used to test the efficacy of new medical countermeasures (MCMs), the high mortality rates associated with melioidosis raises significant ethical issues concerning treating individuals with new compounds with unknown efficacies. The US Food and Drug Administration (FDA) has formulated a set of guidelines for the licensure of new MCMs to treat diseases in which it would be unethical to test the efficacy of these drugs in humans. The FDA “Animal Rule” 21 CFR 314 calls for consistent, well-characterized B. pseudomallei strains to be used as challenge material in animal models. In order to facilitate the efficacy testing of new MCMs for melioidosis using animal models, we intend to develop a well-characterized panel of strains for use. This panel will comprise of strains that were isolated from human cases, have a low passage history, are virulent in animal models, and are well-characterized phenotypically and genotypically. We have reviewed published and unpublished data on various B. pseudomallei strains to establish an objective method for selecting the strains to be included in the panel of B. pseudomallei strains with attention to five categories: animal infection models, genetic characterization, clinical and passage history, and availability of the strain to the research community. We identified 109 strains with data in at least one of the five categories, scored each strain based on the gathered data and identified six strains as candidate for a B. pseudomallei strain panel.