Detection of lateral gene transfer among microbial genomes

Systematic comparison of genome information processing and boundary recognition tools used for genomic island detection

Genomic islands are fragments of foreign DNA that are found in bacterial and archaeal genomes, and are typically associated with symbiosis or pathogenesis. While numerous genomic island detection methods have been proposed, there has been limited evaluation of the efficiency of the genome information processing and boundary recognition tools. In this study, we conducted a review of the statistical methods involved in genomic signatures, host signature extraction, informative signature selection, divergence measures, and boundary detection steps in genomic island prediction. We compared the performances of these methods on simulated experiments using alien fragments obtained from both artificial and real genomes. Our results indicate that among the nine genomic signatures evaluated, genomic signature frequency and full probability performed the best. However, their performance declined when normalized to their expectations and variances, such as Z-score and composition vector. Based on our experiments of the E. coli genome, we found that the confidence intervals of the window variances achieved the best performance in the signature extraction of the host, with the best confidence interval being 1.5-2 times the standard error. Ordered kurtosis was most effective in selecting informative signatures from a single genome, without requiring prior knowledge from other datasets. Among the three divergence measures evaluated, the two-sample t-test was the most successful, and a non-overlapping window with a small eye window (size 2) was best suited for identifying compositionally distinct regions. Finally, the maximum of the Markovian Jensen-Shannon divergence score, in terms of GC-content bias, was found to make boundary detection faster while maintaining a similar error rate.

Multi-model approaches to phylogenetics: Implications for idealization

Phylogenetic models traditionally represent the history of life as having a strictly-branching tree structure. However, it is becoming increasingly clear that the history of life is often not strictly-branching; lateral gene transfer, endosymbiosis, and hybridization, for example, can all produce lateral branching events. There is thus motivation to allow phylogenetic models to have a reticulate structure. One proposal involves the reconciliation of genealogical discordance. Briefly, this method uses patterns of disagreement - discordance - between trees of different genes to add lateral branching events to phylogenetic trees of taxa, and to estimate the most likely cause of these events. I use this practice to argue for: (1) a need for expanded accounts of multiple-models idealization, (2) a distinction between automatic and manual de-idealization, and (3) recognition that idealization may serve the meso-level aims of science in a different way than hitherto acknowledged.

Genomic Island Prediction via Chi-Square Test and Random Forest Algorithm

Genomic islands are related to microbial adaptation and carry different genomic characteristics from the host. Therefore, many methods have been proposed to detect genomic islands from the rest of the genome by evaluating its sequence composition. Many sequence features have been proposed, but many of them have not been applied to the identification of genomic islands. In this paper, we present a scheme to predict genomic islands using the chi-square test and random forest algorithm. We extract seven kinds of sequence features and select the important features with the chi-square test. All the selected features are then input into the random forest to predict the genome islands. Three experiments and comparison show that the proposed method achieves the best performance. This understanding can be useful to design more powerful method for the genomic island prediction.

2SigFinder: the combined use of small-scale and large-scale statistical testing for genomic island detection from a single genome

BACKGROUND

Genomic islands are associated with microbial adaptations, carrying genomic signatures different from the host. Some methods perform an overall test to identify genomic islands based on their local features. However, regions of different scales will display different genomic features.

RESULTS

We proposed here a novel method "2SigFinder ", the first combined use of small-scale and large-scale statistical testing for genomic island detection. The proposed method was tested by genomic island boundary detection and identification of genomic islands or functional features of real biological data. We also compared the proposed method with the comparative genomics and composition-based approaches. The results indicate that the proposed 2SigFinder is more efficient in identifying genomic islands.

CONCLUSIONS

From real biological data, 2SigFinder identified genomic islands from a single genome and reported robust results across different experiments, without annotated information of genomes or prior knowledge from other datasets. 2SigHunter identified 25 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats from 27 Pathogenicity, 1 tRNA, 2 Virulence and 2 Repeats, and detected 101 Phage and 28 HEG out of 130 Phage and 36 HEGs in S. enterica Typhi CT18, which shows that it is more efficient in detecting functional features associated with GIs.

Alignment-free inference of hierarchical and reticulate phylogenomic relationships

We are amidst an ongoing flood of sequence data arising from the application of high-throughput technologies, and a concomitant fundamental revision in our understanding of how genomes evolve individually and within the biosphere. Workflows for phylogenomic inference must accommodate data that are not only much larger than before, but often more error prone and perhaps misassembled, or not assembled in the first place. Moreover, genomes of microbes, viruses and plasmids evolve not only by tree-like descent with modification but also by incorporating stretches of exogenous DNA. Thus, next-generation phylogenomics must address computational scalability while rethinking the nature of orthogroups, the alignment of multiple sequences and the inference and comparison of trees. New phylogenomic workflows have begun to take shape based on so-called alignment-free (AF) approaches. Here, we review the conceptual foundations of AF phylogenetics for the hierarchical (vertical) and reticulate (lateral) components of genome evolution, focusing on methods based on k-mers. We reflect on what seems to be successful, and on where further development is needed.

MTGIpick allows robust identification of genomic islands from a single genome

Genomic islands (GIs) that are associated with microbial adaptations and carry sequence patterns different from that of the host are sporadically distributed among closely related species. This bias can dominate the signal of interest in GI detection. However, variations still exist among the segments of the host, although no uniform standard exists regarding the best methods of discriminating GIs from the rest of the genome in terms of compositional bias. In the present work, we proposed a robust software, MTGIpick, which used regions with pattern bias showing multiscale difference levels to identify GIs from the host. MTGIpick can identify GIs from a single genome without annotated information of genomes or prior knowledge from other data sets. When real biological data were used, MTGIpick demonstrated better performance than existing methods, as well as revealed potential GIs with accurate sizes missed by existing methods because of a uniform standard. Software and supplementary are freely available at http://bioinfo.zstu.edu.cn/MTGI or https://github.com/bioinfo0706/MTGIpick.

An efficient genomic signature ranking method for genomic island prediction from a single genome

Genomic islands that are associated with microbial adaptations and carry genomic signatures different from that of the host, and thus many methods have been proposed to select the informative genomic signatures from a range of organisms and discriminate genomic islands from the rest of the genome in terms of these signature biases. However, they are of limited use when closely related genomes are unavailable. In the present work, we proposed a kurtosis-based ranking method to select the informative genomic signatures from a single genome. In simulations with alien fragments from artificial and real genomes, the proposed kurtosis-based ranking method efficiently selected the informative genomic signatures from a single genome, without annotated information of genomes or prior knowledge from other datasets. This understanding can be useful to design more powerful method for genomic island detection.

Phylogenetic diversification and developmental implications of poly-(R)-3-hydroxyalkanoate gene cluster assembly in prokaryotes

Many polyhydroxyalkanoates (PHAs) system genes, such as phaC, phaA, phaB, phaR, phaP and phaZ, are often found to be organised in the form of operon-like clusters. In this study, a genome survey was performed to identify such clustered PHA systems among 256 prokaryotic organisms. These data were then used to generate a comprehensive 16S rRNA species tree depicting the phylogenetic distribution of the observed clusters with diverse gene arrangements. In addition, the gene occurrences and physical linkages between PHA system genes were quantitatively estimated. From this, we identified a centrally connected hub gene, i.e. the phaC gene of PHA. Furthermore, a comparative investigation was performed between the clusters of PHA and glycogen, which decoded the role of the hub gene in the cluster organisation of both systems. Together, these findings suggest that the highly connected hub gene might contribute substantively towards the organisation and maintenance of the gene network connectivity in the clusters, particularly in the storage reserve systems.

Robust Inference of Genetic Exchange Communities from Microbial Genomes Using TF-IDF

Bacteria and archaea can exchange genetic material across lineages through processes of lateral genetic transfer (LGT). Collectively, these exchange relationships can be modeled as a network and analyzed using concepts from graph theory. In particular, densely connected regions within an LGT network have been defined as genetic exchange communities (GECs). However, it has been problematic to construct networks in which edges solely represent LGT. Here we apply term frequency-inverse document frequency (TF-IDF), an alignment-free method originating from document analysis, to infer regions of lateral origin in bacterial genomes. We examine four empirical datasets of different size (number of genomes) and phyletic breadth, varying a key parameter (word length k) within bounds established in previous work. We map the inferred lateral regions to genes in recipient genomes, and construct networks in which the nodes are groups of genomes, and the edges natively represent LGT. We then extract maximum and maximal cliques (i.e., GECs) from these graphs, and identify nodes that belong to GECs across a wide range of k. Most surviving lateral transfer has happened within these GECs. Using Gene Ontology enrichment tests we demonstrate that biological processes associated with metabolism, regulation and transport are often over-represented among the genes affected by LGT within these communities. These enrichments are largely robust to change of k.

Impact of genomics on the understanding of microbial evolution and classification: the importance of Darwin's views on classification

Analyses of genome sequences, by some approaches, suggest that the widespread occurrence of horizontal gene transfers (HGTs) in prokaryotes disguises their evolutionary relationships and have led to questioning of the Darwinian model of evolution for prokaryotes. These inferences are critically examined in the light of comparative genome analysis, characteristic synapomorphies, phylogenetic trees and Darwin's views on examining evolutionary relationships. Genome sequences are enabling discovery of numerous molecular markers (synapomorphies) such as conserved signature indels (CSIs) and conserved signature proteins (CSPs), which are distinctive characteristics of different prokaryotic taxa. Based on these molecular markers, exhibiting high degree of specificity and predictive ability, numerous prokaryotic taxa of different ranks, currently identified based on the 16S rRNA gene trees, can now be reliably demarcated in molecular terms. Within all studied groups, multiple CSIs and CSPs have been identified for successive nested clades providing reliable information regarding their hierarchical relationships and these inferences are not affected by HGTs. These results strongly support Darwin's views on evolution and classification and supplement the current phylogenetic framework based on 16S rRNA in important respects. The identified molecular markers provide important means for developing novel diagnostics, therapeutics and for functional studies providing important insights regarding prokaryotic taxa.

Bioprospecting of Plant Growth Promoting Bacilli and Related Genera Prevalent in Soils of Pristine Sacred Groves: Biochemical and Molecular Approach

Bacillus spp. and related genera native to soils of the pristine sacred groves from Meghalaya, India were characterized using biochemical and 16S rRNA gene analysis which revealed dominance of Bacillus, Paenibacillus, Lysinibacillus and Viridibacillus in the groves. Biochemical estimation was carried out for in vitro testing of plant growth promoting traits present in these isolates. PCR screening were performed for plant growth-promoting related genes involved in the biosynthesis of acid phosphatase (AcPho), indolepyruvate decarboxylase (ipdC), 1-aminocyclopropane-1-carboxylate deaminase (accd) and siderophore biosynthesis protein (asbA). 76% of the sacred grove isolates gave an amplified fragment for AcPho. Three of the isolates gave an amplified fragment for IpdC gene. Apart from 2 isolates, all the other isolates including the reference strains were positive for the amplification of the accd gene indicating their potential to produce ACC deaminase enzyme. 42% of the isolates gave an amplified fragment for asbA gene indicating the potential ability of these isolates to produce the catechol type siderophore, petrobactin. Overall findings indicated multiple PGP genetic traits present in these isolates which suggested that these isolates are capable of expressing multiple PGP traits. Phylogenetic and sequence analysis of accd and asbA genes from the isolates revealed that asbA genes from Paenibacillus taichungiensis SG3 and Paenibacillus tylopili SG24 indicated the occurrence of intergeneric horizontal transfer between Paenibacillus and Bacillus.

Massive horizontal gene transfer, strictly vertical inheritance and ancient duplications differentially shape the evolution of Bacillus cereus enterotoxin operons hbl, cytK and nhe

Background

Bacillus cereus sensu lato comprises eight closely related species including the human pathogens Bacillus anthracis and Bacillus cereus. Within B. cereus sensu lato, chromosomally and plasmid-encoded toxins exist. While plasmid-mediated horizontal gene transfer of the emetic toxin, anthrax and insecticidal toxins is known, evolution of enterotoxin genes within the group has not been studied.

Results

We report draft genome assemblies of 25 strains, a phylogenetic network of 142 strains based on ANI derived from genome sequences and a phylogeny based on whole-genome SNP analysis. The data clearly support subdivision of B. cereus sensu lato into seven phylogenetic groups. While group I, V and VII represent B. pseudomycoides, B. toyonensis and B. cytotoxicus, which are distinguishable at species level (ANI border ≥ 96 %), strains ascribed to the other five species do not match phylogenic groups. The chromosomal enterotoxin operons nheABC and hblCDAB are abundant within B. cereus both isolated from infections and from the environment. While the duplicated hbl variant hbl_a is present in 22 % of all strains investigated, duplication of nheABC is extremely rare (0.02 %) and appears to be phylogenetically unstable. Distribution of toxin genes was matched to a master tree based on seven concatenated housekeeping genes, which depicts species relationships in B. cereus sensu lato as accurately as whole-genome comparisons. Comparison to the phylogeny of enterotoxin genes uncovered ample evidence for horizontal transfer of hbl, cytK and plcR, as well as frequent deletion of both toxins and duplication of hbl. No evidence for nhe deletion was found and stable horizontal transfer of nhe is rare. Therefore, evolution of B. cereus enterotoxin operons is shaped unexpectedly different for yet unknown reasons.

Conclusions

Frequent exchange of the pathogenicity factors hbl, cytK and plcR in B. cereus sensu lato appears to be an important mechanism of B. cereus virulence evolution, including so-called probiotic or non-pathogenic species, which might have consequences for risk assessment procedures. In contrast, exclusively vertical inheritance of nhe was observed, and since nhe-negative strains appear to be extremely rare, we suggest that fitness loss may be associated with deletion or horizontal transfer of the nhe operon.

Electronic supplementary material

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

Transfer of energy pathway genes in microbial enhanced biological phosphorus removal communities

BACKGROUND

Lateral gene transfer (LGT) is an important evolutionary process in microbial evolution. In sewage treatment plants, LGT of antibiotic resistance and xenobiotic degradation-related proteins has been suggested, but the role of LGT outside these processes is unknown. Microbial communities involved in Enhanced Biological Phosphorus Removal (EBPR) have been used to treat wastewater in the last 50 years and may provide insights into adaptation to an engineered environment. We introduce two different types of analysis to identify LGT in EBPR sewage communities, based on identifying assembled sequences with more than one strong taxonomic match, and on unusual phylogenetic patterns. We applied these methods to investigate the role of LGT in six energy-related metabolic pathways.

RESULTS

The analyses identified overlapping but non-identical sets of transferred enzymes. All of these were homologous with sequences from known mobile genetic elements, and many were also in close proximity to transposases and integrases in the EBPR data set. The taxonomic method had higher sensitivity than the phylogenetic method, identifying more potential LGTs. Both analyses identified the putative transfer of five enzymes within an Australian community, two in a Danish community, and none in a US-derived culture.

CONCLUSIONS

Our methods were able to identify sequences with unusual phylogenetic or compositional properties as candidate LGT events. The association of these candidates with known mobile elements supports the hypothesis of transfer. The results of our analysis strongly suggest that LGT has influenced the development of functionally important energy-related pathways in EBPR systems, but transfers may be unique to each community due to different operating conditions or taxonomic composition.

A new genome-wide method to track horizontally transferred sequences: application to Drosophila

Because of methodological breakthroughs and the availability of an increasing amount of whole-genome sequence data, horizontal transfers (HTs) in eukaryotes have received much attention recently. Contrary to similar analyses in prokaryotes, most studies in eukaryotes usually investigate particular sequences corresponding to transposable elements (TEs), neglecting the other components of the genome. We present a new methodological framework for the genome-wide detection of all putative horizontally transferred sequences between two species that requires no prior knowledge of the transferred sequences. This method provides a broader picture of HTs in eukaryotes by fully exploiting complete-genome sequence data. In contrast to previous genome-wide approaches, we used a well-defined statistical framework to control for the number of false positives in the results, and we propose two new validation procedures to control for confounding factors. The first validation procedure relies on a comparative analysis with other species of the phylogeny to validate HTs for the nonrepeated sequences detected, whereas the second one built upon the study of the dynamics of the detected TEs. We applied our method to two closely related Drosophila species, Drosophila melanogaster and D. simulans, in which we discovered 10 new HTs in addition to all the HTs previously detected in different studies, which underscores our method’s high sensitivity and specificity. Our results favor the hypothesis of multiple independent HTs of TEs while unraveling a small portion of the network of HTs in the Drosophila phylogeny.

The loose evolutionary relationships between transcription factors and other gene products across prokaryotes

Background

Tests for the evolutionary conservation of associations between genes coding for transcription factors (TFs) and other genes have been limited to a few model organisms due to the lack of experimental information of functional associations in other organisms. We aimed at surmounting this limitation by using the most co-occurring gene pairs as proxies for the most conserved functional interactions available for each gene in a genome. We then used genes predicted to code for TFs to compare their most conserved interactions against the most conserved interactions for the rest of the genes within each prokaryotic genome available.

Results

We plotted profiles of phylogenetic profiles, p-cubic, to compare the maximally scoring interactions of TFs against those of other genes. In most prokaryotes, genes coding for TFs showed lower co-occurrences when compared to other genes. We also show that genes coding for TFs tend to have lower Codon Adaptation Indexes compared to other genes.

Conclusions

The co-occurrence tests suggest that transcriptional regulation evolves quickly in most, if not all, prokaryotes. The Codon Adaptation Index analyses suggest quick gene exchange and rewiring of transcriptional regulation across prokaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-928) contains supplementary material, which is available to authorized users.

Comparative study of the marR genes within the family Enterobacteriaceae

marR genes are members of an ancient family originally identified in Escherichia coli. This family is widely distributed in archaea and bacteria. Homologues of this family have a conserved winged helix fold. MarR proteins are involved in non-specific resistance systems conferring resistance to multiple antibiotics. Extensive studies have shown the importance of MarR proteins in physiology and pathogenicity in Enterobacteria, but little is known about their origin or evolution. In this study, all the marR genes in 43 enterobacterial genomes representing 14 genera were identified, and the phylogenetic relationships and genetic parameters were analyzed. Several major findings were made. Three conserved marR genes originated earlier than Enterobacteriaceae and a geneloss event was found to have taken place in Yersinia pestis Antiqua. Three functional genes, rovA, hor, and slyA, were found to be clear orthologs among Enterobacteriaceae. The copy number of marR genes in Enterobacteriaceae was found to vary from 2 to 11. These marR genes exhibited a faster rate of nucleotide substitution than housekeeping genes did. Specifically, the regions of marR domain were found to be subject to strong purifying selection. The phylogenetic relationship and genetic parameter analyses were consistent with conservation and specificity of marR genes. These dual characters helped MarR to maintain a conserved binding motif and variable C-terminus, which are important to adaptive responses to a number of external stimuli in Enterobacteriaceae.

Horizontal gene transfer in eukaryotic plant pathogens

Gene transfer has been identified as a prevalent and pervasive phenomenon and an important source of genomic innovation in bacteria. The role of gene transfer in microbial eukaryotes seems to be of a reduced magnitude but in some cases can drive important evolutionary innovations, such as new functions that underpin the colonization of different niches. The aim of this review is to summarize published cases that support the hypothesis that horizontal gene transfer (HGT) has played a role in the evolution of phytopathogenic traits in fungi and oomycetes. Our survey of the literature identifies 46 proposed cases of transfer of genes that have a putative or experimentally demonstrable phytopathogenic function. When considering the life-cycle steps through which a pathogen must progress, the majority of the HGTs identified are associated with invading, degrading, and manipulating the host. Taken together, these data suggest HGT has played a role in shaping how fungi and oomycetes colonize plant hosts.

Quartet analysis of putative horizontal gene transfer in Crenarchaeota

Horizontal gene transfers (HGT) between four Crenarchaeota species (Metallosphaera cuprina Ar-4T, Acidianus hospitalis W1T, Vulcanisaeta moutnovskia 768-28T, and Pyrobaculum islandicum DSM 4184T) were investigated with quartet analysis. Strong support was found for individual genes that disagree with the phylogeny of the majority, implying genomic mosaicism. One such gene, a ferredoxin-related gene, was investigated further and incorporated into a larger phylogeny, which provided evidence for HGT of this gene from the Vulcanisaeta lineage to the Acidianus lineage. This is the first application of quartet analysis of HGT for the phylum Crenarchaeota. The results have shown that quartet analysis is a powerful technique to screen homologous sequences for putative HGTs and is useful in visually describing genomic mosaicism and HGT within four taxa.

Where does Neisseria acquire foreign DNA from: an examination of the source of genomic and pathogenic islands and the evolution of the Neisseria genus

Background

Pathogenicity islands (PAIs) or genomic islands (GEIs) are considered to be the result of a recent horizontal transfer. Detecting PAIs/GEIs as well as their putative source can provide insight into the organism’s pathogenicity within its host. Previously we introduced a tool called S-plot which provides a visual representation of the variation in compositional properties across and between genomic sequences. Utilizing S-plot and new functionality developed here, we examined 18 publicly available Neisseria genomes, including strains of both pathogenic and non-pathogenic species, in order to identify regions of unusual compositional properties (RUCPs) using both a sliding window as well as a gene-by-gene approach.

Results

Numerous GEIs and PAIs were identified including virulence genes previously found within the pathogenic Neisseria species. While some genes were conserved amongst all species, only pathogenic species, or an individual species, a number of genes were detected that are unique to an individual strain. While the majority of such genes have an origin unknown, a number of putative sources including pathogenic and capsule-containing bacteria were determined, indicative of gene exchange between Neisseria spp. and other bacteria within their microhabitat. Furthermore, we uncovered evidence that both N. meningitidis and N. gonorrhoeae have separately acquired DNA from their human host. Data suggests that all three Neisseria species have received horizontally transferred elements post-speciation.

Conclusions

Using this approach, we were able to not only find previously identified regions of virulence but also new regions which may be contributing to the virulence of the species. This comparative analysis provides a means for tracing the evolutionary history of the acquisition of foreign DNA within this genus. Looking specifically at the RUCPs present within the 18 genomes considered, a stronger similarity between N. meningitidis and N. lactamica is observed, suggesting that N. meningitidis arose before N. gonorrhoeae.

Occurrence of horizontal gene transfer of P(IB)-type ATPase genes among bacteria isolated from the uranium rich deposit of Domiasiat in North East India

Uranium (U) tolerant aerobic heterotrophs were isolated from the subsurface soils of one of the pre-mined U-rich deposits at Domiasiat located in the north-eastern part of India. On screening of genomic DNA from 62 isolates exhibiting superior U and heavy metal tolerance, 32 isolates were found to be positive for P(IB)-type ATPase genes. Phylogenetic incongruence and anomalous DNA base compositions revealed the acquisition of P(IB)-type ATPase genes by six isolates through horizontal gene transfer (HGT). Three of these instances of HGT appeared to have occurred at inter-phylum level and the other three instances indicated to have taken place at intra-phylum level. This study provides an insight into one of the possible survival strategies that bacteria might employ to adapt to environments rich in uranium and heavy metals.

The mercury resistance operon: from an origin in a geothermal environment to an efficient detoxification machine

Mercuric mercury (Hg[II]) is a highly toxic and mobile element that is likely to have had a pronounced and adverse effect on biology since Earth’s oxygenation ∼2.4 billion years ago due to its high affinity for protein sulfhydryl groups, which upon binding destabilize protein structure and decrease enzyme activity, resulting in a decreased organismal fitness. The central enzyme in the microbial mercury detoxification system is the mercuric reductase (MerA) protein, which catalyzes the reduction of Hg(II) to volatile Hg(0). In addition to MerA, mer operons encode for proteins involved in regulation, Hg binding, and organomercury degradation. Mer-mediated approaches have had broad applications in the bioremediation of mercury-contaminated environments and industrial waste streams. Here, we examine the composition of 272 individual mer operons and quantitatively map the distribution of mer-encoded functions on both taxonomic SSU rRNA gene and MerA phylogenies. The results indicate an origin and early evolution of MerA among thermophilic bacteria and an overall increase in the complexity of mer operons through evolutionary time, suggesting continual gene recruitment and evolution leading to an improved efficiency and functional potential of the Mer detoxification system. Consistent with a positive relationship between the evolutionary history and topology of MerA and SSU rRNA gene phylogenies (Mantel R = 0.81, p < 0.01), the distribution of the majority of mer functions, when mapped on these phylograms, indicates an overall tendency to inherit mer-encoded functions through vertical descent. However, individual mer functions display evidence of a variable degree of vertical inheritance, with several genes exhibiting strong evidence for acquisition via lateral gene transfer and/or gene loss. Collectively, these data suggest that (i) mer has evolved from a simple system in geothermal environments to a widely distributed and more complex and efficient detoxification system, and (ii) merA is a suitable biomarker for examining the functional diversity of Hg detoxification and for predicting the composition of mer operons in natural environments.

Evolutionary dynamics of small RNAs in 27 Escherichia coli and Shigella genomes

Small RNAs (sRNAs) are widespread in bacteria and play critical roles in regulating physiological processes. They are best characterized in Escherichia coli K-12 MG1655, where 83 sRNAs constitute nearly 2% of the gene complement. Most sRNAs act by base pairing with a target mRNA, modulating its translation and/or stability; many of these RNAs share only limited complementarity to their mRNA target, and require the chaperone Hfq to facilitate base pairing. Little is known about the evolutionary dynamics of bacterial sRNAs. Here, we apply phylogenetic and network analyses to investigate the evolutionary processes and principles that govern sRNA gene distribution in 27 E. coli and Shigella genomes. We identify core (encoded in all 27 genomes) and variable sRNAs; more than two-thirds of the E. coli K-12 MG1655 sRNAs are core, whereas the others show patterns of presence and absence that are principally due to genetic loss, not duplication or lateral genetic transfer. We present evidence that variable sRNAs are less tightly integrated into cellular genetic regulatory networks than are the core sRNAs, and that Hfq facilitates posttranscriptional cross talk between the E. coli–Shigella core and variable genomes. Finally, we present evidence that more than 80% of genes targeted by Hfq-associated core sRNAs have been transferred within the E. coli–Shigella clade, and that most of these genes have been transferred intact. These results suggest that Hfq and sRNAs help integrate laterally acquired genes into established regulatory networks.

Comparative genomics of Helicobacter pylori and the human-derived Helicobacter bizzozeronii CIII-1 strain reveal the molecular basis of the zoonotic nature of non-pylori gastric Helicobacter infections in humans

BACKGROUND

The canine Gram-negative Helicobacter bizzozeronii is one of seven species in Helicobacter heilmannii sensu lato that are detected in 0.17-2.3% of the gastric biopsies of human patients with gastric symptoms. At the present, H. bizzozeronii is the only non-pylori gastric Helicobacter sp. cultivated from human patients and is therefore a good alternative model of human gastric Helicobacter disease. We recently sequenced the genome of the H. bizzozeronii human strain CIII-1, isolated in 2008 from a 47-year old Finnish woman suffering from severe dyspeptic symptoms. In this study, we performed a detailed comparative genome analysis with H. pylori, providing new insights into non-pylori Helicobacter infections and the mechanisms of transmission between the primary animal host and humans.

RESULTS

H. bizzozeronii possesses all the genes necessary for its specialised life in the stomach. However, H. bizzozeronii differs from H. pylori by having a wider metabolic flexibility in terms of its energy sources and electron transport chain. Moreover, H. bizzozeronii harbours a higher number of methyl-accepting chemotaxis proteins, allowing it to respond to a wider spectrum of environmental signals. In this study, H. bizzozeronii has been shown to have high level of genome plasticity. We were able to identify a total of 43 contingency genes, 5 insertion sequences (ISs), 22 mini-IS elements, 1 genomic island and a putative prophage. Although H. bizzozeronii lacks homologues of some of the major H. pylori virulence genes, other candidate virulence factors are present. In particular, we identified a polysaccharide lyase (HBZC1_15820) as a potential new virulence factor of H. bizzozeronii.

CONCLUSIONS

The comparative genome analysis performed in this study increased the knowledge of the biology of gastric Helicobacter species. In particular, we propose the hypothesis that the high metabolic versatility and the ability to react to a range of environmental signals, factors which differentiate H. bizzozeronii as well as H. felis and H. suis from H. pylori, are the molecular basis of the of the zoonotic nature of H. heilmannii sensu lato infection in humans.

Evolution of patchily distributed proteins shared between eukaryotes and prokaryotes: Dictyostelium as a case study

Protein families are often patchily distributed in the tree of life; they are present in distantly related organisms, but absent in more closely related lineages. This could either be the result of lateral gene transfer between ancestors of organisms that encode them, or losses in the lineages that lack them. Here a novel approach is developed to study the evolution of patchily distributed proteins shared between prokaryotes and eukaryotes. Proteins encoded in the genome of cellular slime mold Dictyostelium discoideum and a restricted number of other lineages, including at least one prokaryote, were identified. Analyses of the phylogenetic distribution of 49 such patchily distributed protein families showed conflicts with organismal phylogenies; 25 are shared with the distantly related amoeboflagellate Naegleria (Excavata), whereas only two are present in the more closely related Entamoeba. Most protein families show unexpected topologies in phylogenetic analyses; eukaryotes are polyphyletic in 85% of the trees. These observations suggest that gene transfers have been an important mechanism for the distribution of patchily distributed proteins across all domains of life. Further studies of this exchangeable gene fraction are needed for a better understanding of the origin and evolution of eukaryotic genes and the diversification process of eukaryotes.

Towards more robust methods of alien gene detection

Because the properties of horizontally-transferred genes will reflect the mutational proclivities of their donor genomes, they often show atypical compositional properties relative to native genes. Parametric methods use these discrepancies to identify bacterial genes recently acquired by horizontal transfer. However, compositional patterns of native genes vary stochastically, leaving no clear boundary between typical and atypical genes. As a result, while strongly atypical genes are readily identified as alien, genes of ambiguous character are poorly classified when a single threshold separates typical and atypical genes. This limitation affects all parametric methods that examine genes independently, and escaping it requires the use of additional genomic information. We propose that the performance of all parametric methods can be improved by using a multiple-threshold approach. First, strongly atypical alien genes and strongly typical native genes would be identified using conservative thresholds. Genes with ambiguous compositional features would then be classified by examining gene context, including the class (native or alien) of flanking genes. By including additional genomic information in a multiple-threshold framework, we observed a remarkable improvement in the performance of several popular, but algorithmically distinct, methods for alien gene detection.

Evolution of northeastern and midwestern Borrelia burgdorferi, United States

Differences in animal ecology or human behavior may account for differences in human incidence in the 2 regions.

The per capita incidence of human Lyme disease in the northeastern United States is more than twice that in the Midwest. However, the prevalence of Borrelia burgdorferi, the bacterium that causes Lyme disease, in the tick vector is nearly identical in the 2 regions. The disparity in human Lyme disease incidence may result from a disparity in the human invasiveness of the bacteria in the Northeast and Midwest caused by fundamentally different evolutionary histories. B. burgdorferi populations in the Northeast and Midwest are geographically isolated, enabling evolutionary divergence in human invasiveness. However, we found that B. burgdorferi populations in the Northeast and Midwest shared a recent common ancestor, which suggests that substantial evolutionary divergence in human invasiveness has not occurred. We propose that differences in either animal ecology or human behavior are the root cause of the differences in human incidence between the 2 regions.

A benchmark of parametric methods for horizontal transfers detection

Horizontal gene transfer (HGT) has appeared to be of importance for prokaryotic species evolution. As a consequence numerous parametric methods, using only the information embedded in the genomes, have been designed to detect HGTs. Numerous reports of incongruencies in results of the different methods applied to the same genomes were published. The use of artificial genomes in which all HGT parameters are controlled allows testing different methods in the same conditions. The results of this benchmark concerning 16 representative parametric methods showed a great variety of efficiencies. Some methods work very poorly whatever the type of HGTs and some depend on the conditions or on the metrics used. The best methods in terms of total errors were those using tetranucleotides as criterion for the window methods or those using codon usage for gene based methods and the Kullback-Leibler divergence metric. Window methods are very sensitive but less specific and detect badly lone isolated gene. On the other hand gene based methods are often very specific but lack of sensitivity. We propose using two methods in combination to get the best of each category, a gene based one for specificity and a window based one for sensitivity.

Comparison of phylogenetic trees through alignment of embedded evolutionary distances

Background

The understanding of evolutionary relationships is a fundamental aspect of modern biology, with the phylogenetic tree being a primary tool for describing these associations. However, comparison of trees for the purpose of assessing similarity and the quantification of various biological processes remains a significant challenge.

Results

We describe a novel approach for the comparison of phylogenetic distance information based on the alignment of representative high-dimensional embeddings (xCEED: Comparison of Embedded Evolutionary Distances). The xCEED methodology, which utilizes multidimensional scaling and Procrustes-related superimposition approaches, provides the ability to measure the global similarity between trees as well as incongruities between them. We demonstrate the application of this approach to the prediction of coevolving protein interactions and demonstrate its improved performance over the mirrortree, tol-mirrortree, phylogenetic vector projection, and partial correlation approaches. Furthermore, we show its applicability to both the detection of horizontal gene transfer events as well as its potential use in the prediction of interaction specificity between a pair of multigene families.

Conclusions

These approaches provide additional tools for the study of phylogenetic trees and associated evolutionary processes. Source code is available at http://gomezlab.bme.unc.edu/tools.

Horizontal gene transfer in evolution: facts and challenges

The contribution of horizontal gene transfer to evolution has been controversial since it was suggested to be a force driving evolution in the microbial world. In this paper, I review the current standpoint on horizontal gene transfer in evolutionary thinking and discuss how important horizontal gene transfer is in evolution in the broad sense, and particularly in prokaryotic evolution. I review recent literature, asking, first, which processes are involved in the evolutionary success of transferred genes and, secondly, about the extent of horizontal gene transfer towards different evolutionary times. Moreover, I discuss the feasibility of reconstructing ancient phylogenetic relationships in the face of horizontal gene transfer. Finally, I discuss how horizontal gene transfer fits in the current neo-Darwinian evolutionary paradigm and conclude there is a need for a new evolutionary paradigm that includes horizontal gene transfer as well as other mechanisms in the explanation of evolution.

Evolutionary constraints permeate large metabolic networks

BACKGROUND

Metabolic networks show great evolutionary plasticity, because they can differ substantially even among closely related prokaryotes. Any one metabolic network can also effectively compensate for the blockage of individual reactions by rerouting metabolic flux through other pathways. These observations, together with the continual discovery of new microbial metabolic pathways and enzymes, raise the possibility that metabolic networks are only weakly constrained in changing their complement of enzymatic reactions.

RESULTS

To ask whether this is the case, I characterized pairwise and higher-order associations in the co-occurrence of genes encoding metabolic enzymes in more than 200 completely sequenced representatives of prokaryotic genera. The majority of reactions show constrained evolution. Specifically, genes encoding most reactions tend to co-occur with genes encoding other reaction(s). Constrained reaction pairs occur in small sets whose number is substantially greater than expected by chance alone. Most such sets are associated with single biochemical pathways. The respective genes are not always tightly linked, which renders horizontal co-transfer of constrained reaction sets an unlikely sole cause for these patterns of association.

CONCLUSION

Even a limited number of available genomes suffices to show that metabolic network evolution is highly constrained by reaction combinations that are favored by natural selection. With increasing numbers of completely sequenced genomes, an evolutionary constraint-based approach may enable a detailed characterization of co-evolving metabolic modules.

Phylogeny of prokaryotes: does it exist and why should we care?

Understanding microbial evolution is essential for gathering information on the most ancient events in the history of Life on our planet. Nevertheless, the idea that it is impossible to reconstruct the evolutionary history of prokaryotes because of horizontal gene transfer has become very popular. We review this important debate and how it can be solved.

Phylogenomic analysis demonstrates a pattern of rare and ancient horizontal gene transfer between plants and fungi

Horizontal gene transfer (HGT) describes the transmission of genetic material across species boundaries and is an important evolutionary phenomenon in the ancestry of many microbes. The role of HGT in plant evolutionary history is, however, largely unexplored. Here, we compare the genomes of six plant species with those of 159 prokaryotic and eukaryotic species and identify 1689 genes that show the highest similarity to corresponding genes from fungi. We constructed a phylogeny for all 1689 genes identified and all homolog groups available from the rice (Oryza sativa) genome (3177 gene families) and used these to define 14 candidate plant-fungi HGT events. Comprehensive phylogenetic analyses of these 14 data sets, using methods that account for site rate heterogeneity, demonstrated support for nine HGT events, demonstrating an infrequent pattern of HGT between plants and fungi. Five HGTs were fungi-to-plant transfers and four were plant-to-fungi HGTs. None of the fungal-to-plant HGTs involved angiosperm recipients. These results alter the current view of organismal barriers to HGT, suggesting that phagotrophy, the consumption of a whole cell by another, is not necessarily a prerequisite for HGT between eukaryotes. Putative functional annotation of the HGT candidate genes suggests that two fungi-to-plant transfers have added phenotypes important for life in a soil environment. Our study suggests that genetic exchange between plants and fungi is exceedingly rare, particularly among the angiosperms, but has occurred during their evolutionary history and added important metabolic traits to plant lineages.

Detection and quantitative assessment of horizontal gene transfer

This chapter discusses the pros and cons of the existing computational methods for the detection of horizontal (or lateral) gene transfer and highlights the genome-wide studies utilizing these methods. The impact of horizontal gene transfer (HGT) on prokaryote genome evolution is discussed.

Bioinformatic perspectives on NRPS/PKS megasynthases: advances and challenges

The increased understanding of both fundamental principles and mechanistic variations of NRPS/PKS megasynthases along with the unprecedented availability of microbial sequences has inspired a number of in silico studies of both enzyme families. The insights that can be extracted from these analyses go far beyond a rough classification of data and have turned bioinformatics into a frontier field of natural products research. As databases are flooded with NRPS/PKS gene sequence of microbial genomes and metagenomes, increasingly reliable structural prediction methods can help to uncover hidden treasures. Already, phylogenetic analyses have revealed that NRPS/PKS pathways should not simply be regarded as enzyme complexes, specifically evolved to product a selected natural product. Rather, they represent a collection of genetic opinions, allowing biosynthetic pathways to be shuffled in a process of perpetual chemical innovations and pathways diversification in nature can give impulses for specificities, protein interactions and genetic engineering of libraries of novel peptides and polyketides. The successful translation of the knowledge obtained from bioinformatic dissection of NRPS/PKS megasynthases into new techniques for drug discovery and design remain challenges for the future.

Horizontal gene transfer in cyanobacterial signature genes

Comparison of 15 phylogenetically diverse cyanobacterial genomes identified an updated list of 183 signature genes that are widely found in cyanobacteria but absent in non-cyanobacterial species. These signature genes comprise the unique portion of the core cyanobacterial phenotype, and their absence from other lineages implies that if they arose by horizontal gene transfer (HGT), it likely occurred before the last shared cyanobacterial ancestor. A remaining issue is whether or not these signature genes would be relatively immune to HGT within the cyanobacterial lineage. Phylogenetic trees for each signature gene were constructed and compared to cyanobacterial groupings based on 16S rRNA sequences, with clear incongruence considered indicative of HGT. Approximately 18% of the signature genes exhibited such anomalies, indicating that the incidence of inter-lineage HGT has been significant. A preliminary analysis of intra-lineage transfer was conducted using four Synechococcus/Prochlorococcus species. In this case, it was found that 13% of the signature genes had likely been involved in within group HGT. In order to compare this level of likely HGT to other gene types, the analysis was extended to 1380 genes shared by the four Synechococcus/Prochlorococcus species. Successful HGT events appear to be most frequent among genes involved in photosynthesis/respiration and genes of unknown function, many of which are signature genes. This is consistent with the hypothesis that genes that most directly effect competition and adaptation of similar species in neighboring niches would be most usefully transferred. Such genes may be more easily integrated into a new genomic environment due to close similarities in regulatory circuits. In summary, signature genes are not immune from HGT and in fact may be favored candidates for HGT among closely related cyanobacterial strains.

A unique combination of genetic systems for the synthesis of trehalose in Rubrobacter xylanophilus: properties of a rare actinobacterial TreT

Trehalose is the primary organic solute in Rubrobacter xylanophilus under all conditions tested, including those for optimal growth. We detected genes of four different pathways for trehalose synthesis in the genome of this organism, namely, the trehalose-6-phosphate synthase (Tps)/trehalose-6-phosphate phosphatase (Tpp), TreS, TreY/TreZ, and TreT pathways. Moreover, R. xylanophilus is the only known member of the phylum Actinobacteria to harbor TreT. The Tps sequence is typically bacterial, but the Tpp sequence is closely related to eukaryotic counterparts. Both the Tps/Tpp and the TreT pathways were active in vivo, while the TreS and the TreY/TreZ pathways were not active under the growth conditions tested and appear not to contribute to the levels of trehalose observed. The genes from the active pathways were functionally expressed in Escherichia coli, and Tps was found to be highly specific for GDP-glucose, a rare feature among these enzymes. The trehalose-6-phosphate formed was specifically dephosphorylated to trehalose by Tpp. The recombinant TreT synthesized trehalose from different nucleoside diphosphate-glucose donors and glucose, but the activity in R. xylanophilus cell extracts was specific for ADP-glucose. The TreT could also catalyze trehalose hydrolysis in the presence of ADP, but with a very high K(m). Here, we functionally characterize two systems for the synthesis of trehalose in R. xylanophilus, a representative of an ancient lineage of the actinobacteria, and discuss a possible scenario for the exceptional occurrence of treT in this extremophilic bacterium.