Ancient lateral gene transfer in the evolution of Bdellovibrio bacteriovorus

Leveraging comparative genomics to uncover alien genes in bacterial genomes

A significant challenge in bacterial genomics is to catalogue genes acquired through the evolutionary process of horizontal gene transfer (HGT). Both comparative genomics and sequence composition-based methods have often been invoked to quantify horizontally acquired genes in bacterial genomes. Comparative genomics methods rely on completely sequenced genomes and therefore the confidence in their predictions increases as the databases become more enriched in completely sequenced genomes. Recent developments including in microbial genome sequencing call for reassessment of alien genes based on information-rich resources currently available. We revisited the comparative genomics approach and developed a new algorithm for alien gene detection. Our algorithm compared favourably with the existing comparative genomics-based methods and is capable of detecting both recent and ancient transfers. It can be used as a standalone tool or in concert with other complementary algorithms for comprehensively cataloguing alien genes in bacterial genomes.

Bdellovibrio reynosensis sp. nov., from a Mexico soil sample

A novel predatory bacterium, strain LBG001T, has been isolated from Reynosa, Mexico. The 16S rRNA shares approximately 97 % sequence identity with many reported strains in the genus Bdellovibrio including the type strain Bdellovibrio bacteriovorus HD100T. Phylogenetic trees based on the 16S rRNA gene and on 30 concatenated housekeeping genes or core genes showed that LBG001T is on a separate branch from the B. bacteriovorus group. LBG0001T has a genome size of 3 582 323 bp with a G+C content of 43.1 mol %. The average nucleotide identity, average amino acid identity and digital DNA-DNA hybridization values with other members of the genus Bdellovibrio (<79, <72 and <17 %, respectively) qualifies the strain to represent a new species in the genus. Strain LBG001T formed visible plaques on all 10 tested Gram-negative bacterial species. The phenotypic characteristics, phylogenetic analysis and genomic taxonomic studies support the classification of the strain as representing a new species for which the name Bdellovibrio reynosensis sp. nov. is proposed. The type strain is LBG001T(=ATCC TSD-288T =CM-CNRG 0932T).

Myxobacterial Genomics and Post-Genomics: A Review of Genome Biology, Genome Sequences and Related 'Omics Studies

Myxobacteria are fascinating and complex microbes. They prey upon other members of the soil microbiome by secreting antimicrobial proteins and metabolites, and will undergo multicellular development if starved. The genome sequence of the model myxobacterium Myxococcus xanthus DK1622 was published in 2006 and 15 years later, 163 myxobacterial genome sequences have now been made public. This explosion in genomic data has enabled comparative genomics analyses to be performed across the taxon, providing important insights into myxobacterial gene conservation and evolution. The availability of myxobacterial genome sequences has allowed system-wide functional genomic investigations into entire classes of genes. It has also enabled post-genomic technologies to be applied to myxobacteria, including transcriptome analyses (microarrays and RNA-seq), proteome studies (gel-based and gel-free), investigations into protein–DNA interactions (ChIP-seq) and metabolism. Here, we review myxobacterial genome sequencing, and summarise the insights into myxobacterial biology that have emerged as a result. We also outline the application of functional genomics and post-genomic approaches in myxobacterial research, highlighting important findings to emerge from seminal studies. The review also provides a comprehensive guide to the genomic datasets available in mid-2021 for myxobacteria (including 24 genomes that we have sequenced and which are described here for the first time).

Bdellovibrio and like organisms: current understanding and knowledge gaps of the smallest cellular hunters of the microbial world

Almost sixty years ago, Bdellovibrio and like organisms (BALOs) were discovered as the first obligate bacterial predators of other bacteria known to science. Since then, they were shown to be diverse and ubiquitous in the environment, and to bear astonishing ecological, physiological, and metabolic capabilities. The last decade has seen important strides made in understanding the mechanistic basis of their life cycle, the dynamics of their interactions with prey, along with significant developments towards their use in medicine, agriculture, and industry. This review details these achievements, identify current understanding and knowledge gaps to encourage and guide future BALO research.

Using modern low-oxygen marine ecosystems to understand the nitrogen cycle of the Paleo- and Mesoproterozoic oceans

During the productive Paleoproterozoic (2.4-1.8 Ga) and less productive Mesoproterozoic (1.8-1.0 Ga), the ocean was suboxic to anoxic and multicellular organisms had not yet evolved. Here, we link geologic information about the Proterozoic ocean to microbial processes in modern low-oxygen systems. High iron concentrations and rates of Fe cycling in the Proterozoic are the largest differences from modern oxygen-deficient zones. In anoxic waters, which composed most of the Paleoproterozoic and ~40% of the Mesoproterozoic ocean, nitrogen cycling dominated. Rates of N₂ production by denitrification and anammox were likely linked to sinking organic matter fluxes and in situ primary productivity under anoxic conditions. Additionally autotrophic denitrifiers could have used reduced iron or methane. 50% of the Mesoproterozoic ocean may have been suboxic, promoting nitrification and metal oxidation in the suboxic water and N₂ O and N₂ production by partial and complete denitrification in anoxic zones in organic aggregates. Sulfidic conditions may have composed ~10% of the Mesoproterozoic ocean focused along continental margins. Due to low nitrate concentrations in offshore regions, anammox bacteria likely dominated N₂ production immediately above sulfidic zones, but in coastal regions, higher nitrate concentrations probably promoted complete S-oxidizing autotrophic denitrification at the sulfide interface.

Variation in genome content and predatory phenotypes between Bdellovibrio sp. NC01 isolated from soil and B. bacteriovorus type strain HD100

Defining phenotypic and associated genotypic variation among Bdellovibrio may further our understanding of how this genus attacks and kills different Gram-negative bacteria. We isolated Bdellovibrio sp. NC01 from soil. Analysis of 16S rRNA gene sequences and average amino acid identity showed that NC01 belongs to a different species than the type species bacteriovorus. By clustering amino acid sequences from completely sequenced Bdellovibrio and comparing the resulting orthologue groups to a previously published analysis, we defined a 'core genome' of 778 protein-coding genes and identified four protein-coding genes that appeared to be missing only in NC01. To determine how horizontal gene transfer (HGT) may have impacted NC01 genome evolution, we performed genome-wide comparisons of Bdellovibrio nucleotide sequences, which indicated that eight NC01 genomic regions were likely acquired by HGT. To investigate how genome variation may impact predation, we compared protein-coding gene content between NC01 and the B. bacteriovorus type strain HD100, focusing on genes implicated as important in successful killing of prey. Of these, NC01 is missing ten genes that may play roles in lytic activity during predation. Compared to HD100, NC01 kills fewer tested prey strains and kills Escherichia coli ML35 less efficiently. NC01 causes a smaller log reduction in ML35, after which the prey population recovers and the NC01 population decreases. In addition, NC01 forms turbid plaques on lawns of E. coli ML35, in contrast to clear plaques formed by HD100. Linking phenotypic variation in interactions between Bdellovibrio and Gram-negative bacteria with underlying Bdellovibrio genome variation is valuable for understanding the ecological significance of predatory bacteria and evaluating their effectiveness in clinical applications.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

Whole-Genome Sequencing and Comparative Genome Analysis Provided Insight into the Predatory Features and Genetic Diversity of Two Bdellovibrio Species Isolated from Soil

Bdellovibrio spp. are predatory bacteria with great potential as antimicrobial agents. Studies have shown that members of the genus Bdellovibrio exhibit peculiar characteristics that influence their ecological adaptations. In this study, whole genomes of two different Bdellovibrio spp. designated SKB1291214 and SSB218315 isolated from soil were sequenced. The core genes shared by all the Bdellovibrio spp. considered for the pangenome analysis including the epibiotic B. exovorus were 795. The number of unique genes identified in Bdellovibrio spp. SKB1291214, SSB218315, W, and B. exovorus JJS was 1343, 113, 857, and 1572, respectively. These unique genes encode hydrolytic, chemotaxis, and transporter proteins which might be useful for predation in the Bdellovibrio strains. Furthermore, the two Bdellovibrio strains exhibited differences based on the % GC content, amino acid identity, and 16S rRNA gene sequence. The 16S rRNA gene sequence of Bdellovibrio sp. SKB1291214 shared 99% identity with that of an uncultured Bdellovibrio sp. clone 12L 106 (a pairwise distance of 0.008) and 95-97% identity (a pairwise distance of 0.043) with that of other culturable terrestrial Bdellovibrio spp., including strain SSB218315. In Bdellovibrio sp. SKB1291214, 174 bp sequence was inserted at the host interaction (hit) locus region usually attributed to prey attachment, invasion, and development of host independent Bdellovibrio phenotypes. Also, a gene equivalent to Bd0108 in B. bacteriovorus HD100 was not conserved in Bdellovibrio sp. SKB1291214. The results of this study provided information on the genetic characteristics and diversity of the genus Bdellovibrio that can contribute to their successful applications as a biocontrol agent.

Keeping it together: absence of genetic variation and DNA incorporation by the predatory bacteria Micavibrio aeruginosavorus and Bdellovibrio bacteriovorus during predation

The use of predatory bacteria as a potential live therapeutic to control human infection is gaining increased attention. Earlier work with Micavibrio spp. and Bdellovibrio spp. has demonstrated the ability of these predators to control drug-resistant Gram-negative pathogens, Tier-1 select agents and biofilms. Additional studies also confirmed that introducing high doses of the predators into animals does not negatively impact animal well-being and might assist in reducing bacterial burden in vivo. The survival of predators requires extreme proximity to the prey cell, which might bring about horizontal transfer of genetic material, such as genes encoding for pathogenic genetic islands that would indirectly facilitate the spread of genetic material to other organisms. In this study, we examined the genetic makeup of several lab isolates of the predators Bdellovibriobacteriovorus and Micavibrioaeruginosavorus that were cultured repeatedly and stored over a course of 13 years. We also conducted controlled experiments in which the predators were sequentially co-cultured on Klebsiella pneumoniae followed by genetic analysis of the predator. In both cases, we saw little genetic variation and no evidence of horizontally transferred chromosomal DNA from the prey during predator-prey interaction. Culturing the predators repeatedly did not cause any change in predation efficacy.

Geosmithia-Ophiostoma: a New Fungus-Fungus Association

In Europe as in North America, elms are devastated by Dutch elm disease (DED), caused by the alien ascomycete Ophiostoma novo-ulmi. Pathogen dispersal and transmission are ensured by local species of bark beetles, which established a novel association with the fungus. Elm bark beetles also transport the Geosmithia fungi genus that is found in scolytids' galleries colonized by O. novo-ulmi. Widespread horizontal gene transfer between O. novo-ulmi and Geosmithia was recently observed. In order to define the relation between these two fungi in the DED pathosystem, O. novo-ulmi and Geosmithia species from elm, including a GFP-tagged strain, were grown in dual culture and mycelial interactions were observed by light and fluorescence microscopy. Growth and sporulation of O. novo-ulmi in the absence or presence of Geosmithia were compared. The impact of Geosmithia on DED severity was tested in vivo by co-inoculating Geosmithia and O. novo-ulmi in elms. A close and stable relation was observed between the two fungi, which may be classified as mycoparasitism by Geosmithia on O. novo-ulmi. These results prove the existence of a new component in the complex of organisms involved in DED, which might be capable of reducing the disease impact.

Prey Range and Genome Evolution of Halobacteriovorax marinus Predatory Bacteria from an Estuary

Predatory bacteria attack and digest other bacteria and therefore may play a role in shaping microbial communities. To investigate phenotypic and genotypic variation in saltwater-adapted predatory bacteria, we isolated Halobacteriovorax marinus BE01 from an estuary in Rhode Island, assayed whether it could attack different prey bacteria, and sequenced and analyzed its genome. We found that BE01 is a prey generalist, attacking bacteria from different phylogenetic groups and environments. Gene order and amino acid sequences are highly conserved between BE01 and the H. marinus type strain, SJ. By comparative genomics, we detected two regions of gene content difference that likely occurred via horizontal gene transfer events. Acquired genes encode functions such as modification of DNA, membrane synthesis and regulation of gene expression. Understanding genome evolution and variation in predation phenotypes among predatory bacteria will inform their development as biocontrol agents and clarify how they impact microbial communities.

Halobacteriovorax strains are saltwater-adapted predatory bacteria that attack Gram-negative bacteria and may play an important role in shaping microbial communities. To understand how Halobacteriovorax strains impact ecosystems and develop them as biocontrol agents, it is important to characterize variation in predation phenotypes and investigate Halobacteriovorax genome evolution. We isolated Halobacteriovorax marinus BE01 from an estuary in Rhode Island using Vibrio from the same site as prey. Small, fast-moving, attack-phase BE01 cells attach to and invade prey cells, consistent with the intraperiplasmic predation strategy of the H. marinus type strain, SJ. BE01 is a prey generalist, forming plaques on Vibrio strains from the estuary, Pseudomonas from soil, and Escherichia coli. Genome analysis revealed extremely high conservation of gene order and amino acid sequences between BE01 and SJ, suggesting strong selective pressure to maintain the genome in this H. marinus lineage. Despite this, we identified two regions of gene content difference that likely resulted from horizontal gene transfer. Analysis of modal codon usage frequencies supports the hypothesis that these regions were acquired from bacteria with different codon usage biases than H. marinus. In one of these regions, BE01 and SJ carry different genes associated with mobile genetic elements. Acquired functions in BE01 include the dnd operon, which encodes a pathway for DNA modification, and a suite of genes involved in membrane synthesis and regulation of gene expression that was likely acquired from another Halobacteriovorax lineage. This analysis provides further evidence that horizontal gene transfer plays an important role in genome evolution in predatory bacteria.

IMPORTANCE Predatory bacteria attack and digest other bacteria and therefore may play a role in shaping microbial communities. To investigate phenotypic and genotypic variation in saltwater-adapted predatory bacteria, we isolated Halobacteriovorax marinus BE01 from an estuary in Rhode Island, assayed whether it could attack different prey bacteria, and sequenced and analyzed its genome. We found that BE01 is a prey generalist, attacking bacteria from different phylogenetic groups and environments. Gene order and amino acid sequences are highly conserved between BE01 and the H. marinus type strain, SJ. By comparative genomics, we detected two regions of gene content difference that likely occurred via horizontal gene transfer events. Acquired genes encode functions such as modification of DNA, membrane synthesis and regulation of gene expression. Understanding genome evolution and variation in predation phenotypes among predatory bacteria will inform their development as biocontrol agents and clarify how they impact microbial communities.

Determination of the Predatory Capability of Bdellovibrio bacteriovorus HD100

Bdellovibrio bacteriovorus HD100 is an obligate predator that preys upon a wide variety of Gram negative bacteria. The biphasic growth cycle of Bdellovibrio includes a free-swimming attack phase and an intraperiplasmic growth phase, where the predator replicates its DNA and grows using the prey as a source of nutrients, finally dividing into individual cells (Sockett, 2009). Due to its obligatory predatory lifestyle, manipulation of Bdellovibrio requires two-member culturing techniques using selected prey microorganisms ( Lambert et al., 2003 ). In this protocol, we describe a detailed workflow to grow and quantify B. bacteriovorus HD100 and its predatory ability, to easily carry out these laborious and time-consuming techniques.

Analysis of gene gain and loss in the evolution of predatory bacteria

Predatory bacteria are ubiquitously distributed in nature in including in aquatic environments, sewage, intestinal tracts of animals and humans, rhizophere and, soils. However, our understanding of their evolutionary history is limited. Results of recent studies have shown that acquiring novel genes is a major force driving bacterial evolution. Therefore, to gain a better understanding of the impact of gene gain and loss in the evolution of bacterial predators, this study employed comparative genomic approaches to identify core-set gene families and species-specific gene families, and model gene gain and loss events among 11 genomes that represented diverse lineages. In total, 1977 gene families were classified. Of these 509 (pattern 11111111111) were present all of the 11 species. Among the non-core set gene families, 52 were present only in saltwater bacteria predators and had no ortholog in the other genomes. Similarly 109 and 44 were present only in the genomes of Micavibrio spp. and Bdellovibrio spp., respectively. In this study, the gain loss mapping engine GLOOME was selected to analyze and estimate the expectations and probabilities of both gain and loss events in the predatory bacteria. In total, 354 gene families were involved in significant gene gain events, and 407 gene families were classified into gene loss events with high supported value. Moreover, 18 families from the core set gene family were identified as putative genes under positive selection. The results of this study suggest that acquisition of particular genes that encode functional proteins in metabolism and cellular processes and signaling, especially ABC systems, may help bacterial predators adapt to surrounding environmental changes and present different predation strategies for survival in their habitats.

Bacterial predation: 75 years and counting!

The first documented study on bacterial predation was carried out using myxobacteria three quarters of a century ago. Since then, many predatory strains, diverse hunting strategies, environmental consequences and potential applications have been reported by groups all over the world. Now we know that predatory bacteria are distributed in a wide variety of environments and that interactions between predatory and non-predatory populations seem to be the most important factor in bacterial selection and mortality in some ecosystems. Bacterial predation has now been proposed as an evolutionary driving force. The structure and diversity of the predatory bacterial community is beginning to be recognized as an important factor in biodiversity due to its potential role in controlling and modelling bacterial populations in the environment. In this paper, we review the current understanding of bacterial predation, going over the strategies used by the main predatory bacteria to kill their prey. We have also reviewed and integrated the accumulated advances of the last 75 years with the interesting new insights that are provided by the analyses of genomes, predatomes, predatosomes and other comparative genomics studies, focusing on potential applications that derive from all of these areas of study.

Ankyrin-mediated self-protection during cell invasion by the bacterial predator Bdellovibrio bacteriovorus

Predatory Bdellovibrio bacteriovorus are natural antimicrobial organisms, killing other bacteria by whole-cell invasion. Self-protection against prey-metabolizing enzymes is important for the evolution of predation. Initial prey entry involves the predator's peptidoglycan DD-endopeptidases, which decrosslink cell walls and prevent wasteful entry by a second predator. Here we identify and characterize a self-protection protein from B. bacteriovorus, Bd3460, which displays an ankyrin-based fold common to intracellular pathogens of eukaryotes. Co-crystal structures reveal Bd3460 complexation of dual targets, binding a conserved epitope of each of the Bd3459 and Bd0816 endopeptidases. Complexation inhibits endopeptidase activity and cell wall decrosslinking in vitro. Self-protection is vital - ΔBd3460 Bdellovibrio deleteriously decrosslink self-peptidoglycan upon invasion, adopt a round morphology, and lose predatory capacity and cellular integrity. Our analysis provides the first mechanistic examination of self-protection in Bdellovibrio, documents protection-multiplicity for products of two different genomic loci, and reveals an important evolutionary adaptation to an invasive predatory bacterial lifestyle.

Genome-wide analysis of myxobacterial two-component systems: genome relatedness and evolutionary changes

BACKGROUND

Two-component systems (TCSs) are abundant prokaryotic signaling pathways, whose evolution is of particular importance because of their role in bacterial pathogenicity. Comparative genomics can provide important insights into the evolution of these genes, but inferences are dependent on the relatedness of the compared genomes. This study investigated the relationship between evolutionary distance and TCS evolution in myxobacterial genomes, of which there are several sequenced examples, of varying relatedness, and which encode large numbers of TCSs.

METHODS

Myxobacterial TCS gene sets were compared, orthologues defined, and changes in TCS properties such as gene organisation, domain architecture and size identified.

RESULTS

Genome relatedness/evolutionary distance was found to have a large effect on the apparent frequency of evolutionary events affecting TCS genes, but not on the relative dominance of different types of mutations. Large (≥1 gene) indels were the most common changes, often giving rise to gene organisation changes. Smaller indels were also common, sometimes changing domain architecture, and/or leading to pseudogene formation. Individuality of myxobacterial TCS gene sets seems primarily due to lineage specific gene loss. However, there is also evidence of extensive acquisition of genes by lateral transfer, with gene duplication also creating new TCS genes.

CONCLUSIONS

This study provides catalogues of myxobacterial TCS gene sets and their orthology relationships, benchmarked against genome relatedness. It also provides insights into the relationship between evolutionary distance and the inference of TCS estudies of TCS evolution beyond the myxobacteriavolution, which may be important for studies of TCS evolutiThe online version of this articleon beyond the myxobacteria.

Inferring horizontal gene transfer

Horizontal or Lateral Gene Transfer (HGT or LGT) is the transmission of portions of genomic DNA between organisms through a process decoupled from vertical inheritance. In the presence of HGT events, different fragments of the genome are the result of different evolutionary histories. This can therefore complicate the investigations of evolutionary relatedness of lineages and species. Also, as HGT can bring into genomes radically different genotypes from distant lineages, or even new genes bearing new functions, it is a major source of phenotypic innovation and a mechanism of niche adaptation. For example, of particular relevance to human health is the lateral transfer of antibiotic resistance and pathogenicity determinants, leading to the emergence of pathogenic lineages. Computational identification of HGT events relies upon the investigation of sequence composition or evolutionary history of genes. Sequence composition-based ("parametric") methods search for deviations from the genomic average, whereas evolutionary history-based ("phylogenetic") approaches identify genes whose evolutionary history significantly differs from that of the host species. The evaluation and benchmarking of HGT inference methods typically rely upon simulated genomes, for which the true history is known. On real data, different methods tend to infer different HGT events, and as a result it can be difficult to ascertain all but simple and clear-cut HGT events.

The futalosine pathway played an important role in menaquinone biosynthesis during early prokaryote evolution

Menaquinone (MK) is an important component of the electron-transfer system in prokaryotes. One of its precursors, 1,4-dihydroxy-2-naphthoate, can be synthesized from chorismate by the classical MK pathway. Interestingly, in some bacteria, chorismate can also be converted to 1,4-dihydroxy-6-naphthoate by four enzymes encoded by mqnABCD in an alternative futalosine pathway. In this study, six crucial enzymes belonging to these two independent nonhomologous pathways were identified in the predicted proteomes of prokaryotes representing a broad phylogenetic distribution. Although the classical MK pathway was found in 32.1% of the proteomes, more than twice the proportion containing the futalosine pathway, the latter was found in a broader taxonomic range of organisms (18 of 31 phyla). The prokaryotes equipped with the classical MK pathway were almost all aerobic or facultatively anaerobic, but those with the futalosine pathway were not only aerobic or facultatively anaerobic but also anaerobic. Phylogenies of enzymes of the classical MK pathway indicated that its genes in archaea were probably acquired by an ancient horizontal gene transfer from bacterial donors. Therefore, the organization of the futalosine pathway likely predated that of the classical MK pathway in the evolutionary history of prokaryotes.

A global transcriptional switch between the attack and growth forms of Bdellovibrio bacteriovorus

Bdellovibrio bacteriovorus is an obligate predator of bacteria ubiquitously found in the environment. Its life cycle is composed of two essential phases: a free-living, non-replicative, fast swimming attack phase (AP) wherein the predator searches for prey; and a non-motile, actively dividing growth phase (GP) in which it consumes the prey. The molecular regulatory mechanisms governing the switch between AP and GP are largely unknown. We used RNA-seq to generate a single-base-resolution map of the Bdellovibrio transcriptome in AP and GP, revealing a specific "AP" transcriptional program, which is largely mutually exclusive of the GP program. Based on the expression map, most genes in the Bdellovibrio genome are classified as "AP only" or "GP only". We experimentally generated a genome-wide map of 140 AP promoters, controlling the majority of AP-specific genes. This revealed a common sigma-like DNA binding site highly similar to the E. coli flagellar genes regulator sigma28 (FliA). Further analyses suggest that FliA has evolved to become a global AP regulator in Bdellovibrio. Our results also reveal a non-coding RNA that is massively expressed in AP. This ncRNA contains a c-di-GMP riboswitch. We suggest it functions as an intracellular reservoir for c-di-GMP, playing a role in the rapid switch from AP to GP.

Being Aquifex aeolicus: Untangling a hyperthermophile's checkered past

Lateral gene transfer (LGT) is an important factor contributing to the evolution of prokaryotic genomes. The Aquificae are a hyperthermophilic bacterial group whose genes show affiliations to many other lineages, including the hyperthermophilic Thermotogae, the Proteobacteria, and the Archaea. Previous phylogenomic analyses focused on Aquifex aeolicus identified Thermotogae and Aquificae either as successive early branches or sisters in a rooted bacterial phylogeny, but many phylogenies and cellular traits have suggested a stronger affiliation with the Epsilonproteobacteria. Different scenarios for the evolution of the Aquificae yield different phylogenetic predictions. Here, we outline these scenarios and consider the fit of the available data, including three sequenced Aquificae genomes, to different sets of predictions. Evidence from phylogenetic profiles and trees suggests that the Epsilonproteobacteria have the strongest affinities with the three Aquificae analyzed. However, this pattern is shown by only a minority of encoded proteins, and the Archaea, many lineages of thermophilic bacteria, and members of genus Clostridium and class Deltaproteobacteria also show strong connections to the Aquificae. The phylogenetic affiliations of different functional subsystems showed strong biases: Most but not all genes implicated in the core translational apparatus tended to group Aquificae with Thermotogae, whereas a wide range of metabolic and cellular processes strongly supported the link between Aquificae and Epsilonproteobacteria. Depending on which sets of genes are privileged, either Thermotogae or Epsilonproteobacteria is the most plausible adjacent lineage to the Aquificae. Both scenarios require massive sharing of genes to explain the history of this enigmatic group, whose history is further complicated by specific affinities of different members of Aquificae to different partner lineages.

By their genes ye shall know them: genomic signatures of predatory bacteria

Predatory bacteria are taxonomically disparate, exhibit diverse predatory strategies and are widely distributed in varied environments. To date, their predatory phenotypes cannot be discerned in genome sequence data thereby limiting our understanding of bacterial predation, and of its impact in nature. Here, we define the 'predatome,' that is, sets of protein families that reflect the phenotypes of predatory bacteria. The proteomes of all sequenced 11 predatory bacteria, including two de novo sequenced genomes, and 19 non-predatory bacteria from across the phylogenetic and ecological landscapes were compared. Protein families discriminating between the two groups were identified and quantified, demonstrating that differences in the proteomes of predatory and non-predatory bacteria are large and significant. This analysis allows predictions to be made, as we show by confirming from genome data an over-looked bacterial predator. The predatome exhibits deficiencies in riboflavin and amino acids biosynthesis, suggesting that predators obtain them from their prey. In contrast, these genomes are highly enriched in adhesins, proteases and particular metabolic proteins, used for binding to, processing and consuming prey, respectively. Strikingly, predators and non-predators differ in isoprenoid biosynthesis: predators use the mevalonate pathway, whereas non-predators, like almost all bacteria, use the DOXP pathway. By defining predatory signatures in bacterial genomes, the predatory potential they encode can be uncovered, filling an essential gap for measuring bacterial predation in nature. Moreover, we suggest that full-genome proteomic comparisons are applicable to other ecological interactions between microbes, and provide a convenient and rational tool for the functional classification of bacteria.

Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria

Background

Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as “obligate predators” because only by mutations, often in gene bd0108, are 1 in ~1x10⁷ of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome.

Results

To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired.

Conclusions

The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

Identification and biochemical evidence of a medium-chain-length polyhydroxyalkanoate depolymerase in the Bdellovibrio bacteriovorus predatory hydrolytic arsenal

The obligate predator Bdellovibrio bacteriovorus HD100 shows a large set of proteases and other hydrolases as part of its hydrolytic arsenal needed for its predatory life cycle. We present genetic and biochemical evidence that open reading frame (ORF) Bd3709 of B. bacteriovorus HD100 encodes a novel medium-chain-length polyhydroxyalkanoate (mcl-PHA) depolymerase (PhaZ(Bd)). The primary structure of PhaZ(Bd) suggests that this enzyme belongs to the α/β-hydrolase fold family and has a typical serine hydrolase catalytic triad (serine-histidine-aspartic acid) in agreement with other PHA depolymerases and lipases. PhaZ(Bd) has been extracellularly produced using different hypersecretor Tol-pal mutants of Escherichia coli and Pseudomonas putida as recombinant hosts. The recombinant PhaZ(Bd) has been characterized, and its biochemical properties have been compared to those of other PHA depolymerases. The enzyme behaves as a serine hydrolase that is inhibited by phenylmethylsulfonyl fluoride. It is also affected by the reducing agent dithiothreitol and nonionic detergents like Tween 80. PhaZ(Bd) is an endoexohydrolase that cleaves both large and small PHA molecules, producing mainly dimers but also monomers and trimers. The enzyme specifically degrades mcl-PHA and is inactive toward short-chain-length polyhydroxyalkanoates (scl-PHA) like polyhydroxybutyrate (PHB). These studies shed light on the potentiality of these predators as sources of new biocatalysts, such as an mcl-PHA depolymerase, for the production of enantiopure hydroxyalkanoic acids and oligomers as building blocks for the synthesis of biobased polymers.

Characterization of the gastrointestinal microbiota in health and inflammatory bowel disease

The enteric bacterial flora play a key role in maintaining health. Inflammatory bowel disease is associated with quantitative and qualitative alterations in the microbiota. Early characterization of the microbiota involved culture-dependent techniques. The advent of metagenomic techniques, however, allows for structural and functional characterization using culture-independent methods. Changes in diversity, together with quantitative alterations in specific bacterial species, have been identified. The functional significance of these changes, and their pathogenic role, remain to be elucidated.

Analysis of the genome and proteome composition of Bdellovibrio bacteriovorus: indication for recent prey-derived horizontal gene transfer

The genome/proteome composition of Bdellovibrio bacteriovorus, the predatory microorganism that preys on other Gram-negative bacteria, has been analyzed. The study elucidates that translational selection plays a major role in genome compositional variation with higher intensity compared to other deltaproteobacteria. Other sources of variations having relatively minor contributions are local GC-bias, horizontal gene transfer and strand-specific mutational bias. The study identifies a group of AT-rich genes with distinct codon composition that is presumably acquired by Bdellovibrio recently from Gram-negative prey-bacteria other than deltaproteobacteria. The proteome composition of this species is influenced by various physico-chemical factors, viz, alcoholicity, residue-charge, aromaticity and hydropathy. Cell-wall-surface-anchor-family (CSAPs) and transporter proteins with distinct amino acid composition and specific secondary-structure also contribute notably to proteome compositional variation. CSAPs, which are low molecular-weight, outer-membrane proteins with highly disordered secondary-structure, have preference toward polar-uncharged residues and cysteine that presumably help in prey-predator interaction by providing particular bonds of attachment.

Assembling networks of microbial genomes using linear programming

Background

Microbial genomes exhibit complex sets of genetic affinities due to lateral genetic transfer. Assessing the relative contributions of parent-to-offspring inheritance and gene sharing is a vital step in understanding the evolutionary origins and modern-day function of an organism, but recovering and showing these relationships is a challenging problem.

Results

We have developed a new approach that uses linear programming to find between-genome relationships, by treating tables of genetic affinities (here, represented by transformed BLAST e-values) as an optimization problem. Validation trials on simulated data demonstrate the effectiveness of the approach in recovering and representing vertical and lateral relationships among genomes. Application of the technique to a set comprising Aquifex aeolicus and 75 other thermophiles showed an important role for large genomes as 'hubs' in the gene sharing network, and suggested that genes are preferentially shared between organisms with similar optimal growth temperatures. We were also able to discover distinct and common genetic contributors to each sequenced representative of genus Pseudomonas.

Conclusions

The linear programming approach we have developed can serve as an effective inference tool in its own right, and can be an efficient first step in a more-intensive phylogenomic analysis.

A database of phylogenetically atypical genes in archaeal and bacterial genomes, identified using the DarkHorse algorithm

Background

The process of horizontal gene transfer (HGT) is believed to be widespread in Bacteria and Archaea, but little comparative data is available addressing its occurrence in complete microbial genomes. Collection of high-quality, automated HGT prediction data based on phylogenetic evidence has previously been impractical for large numbers of genomes at once, due to prohibitive computational demands. DarkHorse, a recently described statistical method for discovering phylogenetically atypical genes on a genome-wide basis, provides a means to solve this problem through lineage probability index (LPI) ranking scores. LPI scores inversely reflect phylogenetic distance between a test amino acid sequence and its closest available database matches. Proteins with low LPI scores are good horizontal gene transfer candidates; those with high scores are not.

Description

The DarkHorse algorithm has been applied to 955 microbial genome sequences, and the results organized into a web-searchable relational database, called the DarkHorse HGT Candidate Resource . Users can select individual genomes or groups of genomes to screen by LPI score, search for protein functions by descriptive annotation or amino acid sequence similarity, or select proteins with unusual G+C composition in their underlying coding sequences. The search engine reports LPI scores for match partners as well as query sequences, providing the opportunity to explore whether potential HGT donor sequences are phylogenetically typical or atypical within their own genomes. This information can be used to predict whether or not sufficient information is available to build a well-supported phylogenetic tree using the potential donor sequence.

Conclusion

The DarkHorse HGT Candidate database provides a powerful, flexible set of tools for identifying phylogenetically atypical proteins, allowing researchers to explore both individual HGT events in single genomes, and large-scale HGT patterns among protein families and genome groups. Although the DarkHorse algorithm cannot, by itself, provide definitive proof of horizontal gene transfer, it is a flexible, powerful tool that can be combined with slower, more rigorous methods in situations where these other methods could not otherwise be applied.

Genome evolution and the emergence of fruiting body development in Myxococcus xanthus

Background

Lateral gene transfer (LGT) is thought to promote speciation in bacteria, though well-defined examples have not been put forward.

Methodology/Principle Findings

We examined the evolutionary history of the genes essential for a trait that defines a phylogenetic order, namely fruiting body development of the Myxococcales. Seventy-eight genes that are essential for Myxococcus xanthus development were examined for LGT. About 73% of the genes exhibit a phylogeny similar to that of the 16S rDNA gene and a codon bias consistent with other M. xanthus genes suggesting vertical transmission. About 22% have an altered codon bias and/or phylogeny suggestive of LGT. The remaining 5% are unique. Genes encoding signal production and sensory transduction were more likely to be transmitted vertically with clear examples of duplication and divergence into multigene families. Genes encoding metabolic enzymes were frequently acquired by LGT. Myxobacteria exhibit aerobic respiration unlike most of the δ Proteobacteria. M. xanthus contains a unique electron transport pathway shaped by LGT of genes for succinate dehydrogenase and three cytochrome oxidase complexes.

Conclusions/Significance

Fruiting body development depends on genes acquired by LGT, particularly those involved in polysaccharide production. We suggest that aerobic growth fostered innovation necessary for development by allowing myxobacteria access to a different gene pool from anaerobic members of the δ Proteobacteria. Habitat destruction and loss of species diversity could restrict the evolution of new bacterial groups by limiting the size of the prospective gene pool.

Predatory mechanisms of Bdellovibrio and like organisms

Bdellovibrio and like organisms (BALOs) are predatory, Gram-negative delta-proteobacteria with a complex developmental lifecycle. In the free-living attack phase they are highly motile and seek out prey bacteria that they invade. The ensuing intracellular growth and replication is characterized by the development of a long filament that septates into individual cells that differentiate further into the flagellated attack-phase bacterium. The prey bacterium is lysed and motile predators are released. BALOs have recently been considered to have potential as living antibiotics. The idea of using predatory bacteria as therapeutic agents to combat pathogenic Gram-negative bacteria is intriguing, as they can prey upon human pathogenic bacteria including Salmonella, Pseudomonas and Escherichia coli. However, our current knowledge of the amazing biology of these prokaryotes that cause the systematic destruction of Gram-negative bacteria is still rather limited. More has to be learned about their predatory lifestyle before their application as therapeutic agents will become feasible.

DarkHorse: a method for genome-wide prediction of horizontal gene transfer

DarkHorse is a new approach to rapid, genome-wide identification and ranking of horizontal transfer candidate proteins.

A new approach to rapid, genome-wide identification and ranking of horizontal transfer candidate proteins is presented. The method is quantitative, reproducible, and computationally undemanding. It can be combined with genomic signature and/or phylogenetic tree-building procedures to improve accuracy and efficiency. The method is also useful for retrospective assessments of horizontal transfer prediction reliability, recognizing orthologous sequences that may have been previously overlooked or unavailable. These features are demonstrated in bacterial, archaeal, and eukaryotic examples.

Evolution of sensory complexity recorded in a myxobacterial genome

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.