Computational and phylogenetic validation of nematode horizontal gene transfer

Sequencing of expressed genes has shown that nematodes, particularly the plant-parasitic nematodes, have genes purportedly acquired from other kingdoms by horizontal gene transfer. The prevailing orthodoxy is that such transfer has been a driving force in the evolution of niche specificity, and a recent paper in BMC Evolutionary Biology that presents a detailed phylogenetic analysis of cellulase genes in the free-living nematode Pristionchus pacificus at the species, genus and family levels substantiates this hypothesis.     

Methods of horizontal gene transfer determination using phylogenetic data

A new approach for comparative analysis of multiple trees reconstructed for representative protein families is proposed. This approach is based on the hypothesis of gene duplication, gene loss and horizontal gene transfer and makes use of stochastic methods and optimization. We present a species tree of 40 prokaryotic organisms obtained by our algorithm on the basis of 132 clusters of orthologous groups of proteins (COGs) from the GenBank of the National Center for Biotechnology Information (USA). We also present a computer technology intended to determine horizontally transferred genes. Some application results of the technology, based on comparative analysis of protein and species trees, are given.     

SPRIT: Identifying horizontal gene transfer in rooted phylogenetic trees

Background  

Phylogenetic trees based on sequences from a set of taxa can be incongruent due to horizontal gene transfer (HGT). By identifying the HGT events, we can reconcile the gene trees and derive a taxon tree that adequately represents the species' evolutionary history. One HGT can be represented by a rooted Subtree Prune and Regraft (RSPR) operation and the number of RSPRs separating two trees corresponds to the minimum number of HGT events. Identifying the minimum number of RSPRs separating two trees is NP-hard, but the problem can be reduced to fixed parameter tractable. A number of heuristic and two exact approaches to identifying the minimum number of RSPRs have been proposed. This is the first implementation delivering an exact solution as well as the intermediate trees connecting the input trees.     

Identification of the horizontal gene transfer based on phylogenetic data

We suggest a new procedure to search for the genes with horizontal transfer events in their evolutionary history. The search is based on analysis of topology difference between the phylogenetic trees of gene (protein) groups and the corresponding phylogenetic species trees. Numeric values are introduced to measure the discrepancy between the trees. This approach was applied to analyze 40 prokaryotic genomes classified into 132 classes of orthologs. This resulted in a list of the candidate genes for which the hypothesis of horizontal transfer in evolution looks true.     

HGT-Gen: a tool for generating a phylogenetic tree with horizontal gene transfer

Horizontal gene transfer (HGT) is a common event in prokaryotic evolution. Therefore, it is very important to consider HGT in the study of molecular evolution of prokaryotes. This is true also for conducting computer simulations of their molecular phylogeny because HGT is known to be a serious disturbing factor for estimating their correct phylogeny. To the best of our knowledge, no existing computer program has generated a phylogenetic tree with HGT from an original phylogenetic tree. We developed a program called HGT-Gen that generates a phylogenetic tree with HGT on the basis of an original phylogenetic tree of a protein or gene. HGT-Gen converts an operational taxonomic unit or a clade from one place to another in a given phylogenetic tree. We have also devised an algorithm to compute the average length between any pair of branches in the tree. It defines and computes the relative evolutionary time to normalize evolutionary time for each lineage. The algorithm can generate an HGT between a pair of donor and acceptor lineages at the same evolutionary time. HGT-Gen is used with a sequence-generating program to evaluate the influence of HGT on the molecular phylogeny of prokaryotes in a computer simulation study.

Availability

The database is available for free at http://www.grl.shizuoka.ac.jp/˜thoriike/HGT-Gen.html     

Genetic and phylogenetic evidence for horizontal gene transfer among ecologically disparate groups of marine Vibrio

Vibrio represents a diverse bacterial genus found in different niches of the marine environment, including numerous genera of marine sponges (phylum Porifera), inhabiting different depths and regions of benthic seas, that are potentially important in driving adaptive change among Vibrio spp. Using 16S rRNA gene sequencing, a previous study showed that sponge‐derived (SD) vibrios clustered with their mainstream counterparts present in shallow, coastal ecosystems, suggesting a genetic relatedness between these populations. Sequences from the topA, ftsZ, mreB, rpoD, rctB and toxR genes were used to investigate the degree of relatedness existing between these two separate populations by examining their phylogenetic and genetic disparity. Phylogenies were constructed from the concatenated sequences of the six housekeeping genes using maximum‐parsimony, maximum‐likelihood and neighbour‐joining algorithms. Genetic recombination was evaluated using the incongruence length difference test, Split decomposition and measuring overall compatibility of sites. This combined technical approach provided evidence that SD Vibrio strains are largely genetically homologous to their shallow‐water counterparts. Moreover, the analyses conducted support the existence of extensive horizontal gene transfer between these two groups, supporting the idea of a single panmictic population structure among vibrios from two seemingly distinct, marine environments.     

Organization of butyrate synthetic genes in human colonic bacteria: phylogenetic conservation and horizontal gene transfer

Butyrate producers constitute an important bacterial group in the human large intestine. Butyryl-CoA is formed from two molecules of acetyl-CoA in a process resembling beta-oxidation in reverse. Three different arrangements of the six genes coding for this pathway have been found in low mol% G+C-content gram-positive human colonic bacteria using DNA sequencing and degenerate PCR. Gene arrangements were strongly conserved within phylogenetic groups defined by 16S rRNA gene sequence relationships. In the case of one of the genes, encoding beta-hydroxybutyryl-CoA dehydrogenase, however, sequence relationships were strongly suggestive of horizontal gene transfer between lineages. The newly identified gene for butyryl-CoA CoA-transferase, which performs the final step in butyrate formation in most known human colonic bacteria, was not closely linked to these central pathway genes.     

Automatic genome-wide reconstruction of phylogenetic gene trees

Gene duplication and divergence is a major evolutionary force. Despite the growing number of fully sequenced genomes, methods for investigating these events on a genome-wide scale are still in their infancy. Here, we present SYNERGY, a novel and scalable algorithm that uses sequence similarity and a given species phylogeny to reconstruct the underlying evolutionary history of all genes in a large group of species. In doing so, SYNERGY resolves homology relations and accurately distinguishes orthologs from paralogs. We applied our approach to a set of nine fully sequenced fungal genomes spanning 150 million years, generating a genome-wide catalog of orthologous groups and corresponding gene trees. Our results are highly accurate when compared to a manually curated gold standard, and are robust to the quality of input according to a novel jackknife confidence scoring. The reconstructed gene trees provide a comprehensive view of gene evolution on a genomic scale. Our approach can be applied to any set of sequenced eukaryotic species with a known phylogeny, and opens the way to systematic studies of the evolution of individual genes, molecular systems and whole genomes. Supplementary information: Supplementary data are available at Bioinformatics online.     

Detecting highways of horizontal gene transfer

In a horizontal gene transfer (HGT) event, a gene is transferred between two species that do not have an ancestor-descendant relationship. Typically, no more than a few genes are horizontally transferred between any two species. However, several studies identified pairs of species between which many different genes were horizontally transferred. Such a pair is said to be linked by a highway of gene sharing. We present a method for inferring such highways. Our method is based on the fact that the evolutionary histories of horizontally transferred genes disagree with the corresponding species phylogeny. Specifically, given a set of gene trees and a trusted rooted species tree, each gene tree is first decomposed into its constituent quartet trees and the quartets that are inconsistent with the species tree are identified. Our method finds a pair of species such that a highway between them explains the largest (normalized) fraction of inconsistent quartets. For a problem on n species and m input quartet trees, we give an efficient O(m + n(2))-time algorithm for detecting highways, which is optimal with respect to the quartets input size. An application of our method to a dataset of 1128 genes from 11 cyanobacterial species, as well as to simulated datasets, illustrates the efficacy of our method.     

Monitoring and modeling horizontal gene transfer

Monitoring efforts have failed to identify horizontal gene transfer (HGT) events occurring from transgenic plants into bacterial communities in soil or intestinal environments. The lack of such observations is frequently cited in biosafety literature and by regulatory risk assessment. Our analysis of the sensitivity of current monitoring efforts shows that studies to date have examined potential HGT events occurring in less than 2 g of sample material, when combined. Moreover, a population genetic model predicts that rare bacterial transformants acquiring transgenes require years of growth to out-compete wild-type bacteria. Time of sampling is there-fore crucial to the useful implementation of monitoring. A population genetic approach is advocated for elucidating the necessary sample sizes and times of sampling for monitoring HGT into large bacterial populations. Major changes in current monitoring approaches are needed, including explicit consideration of the population size of exposed bacteria, the bacterial generation time, the strength of selection acting on the transgene-carrying bacteria, and the sample size necessary to verify or falsify the HGT hypotheses tested.     

Eukaryotic signalling domain homologues in archaea and bacteria. Ancient ancestry and horizontal gene transfer.

Phyletic distributions of eukaryotic signalling domains were studied using recently developed sensitive methods for protein sequence analysis, with an emphasis on the detection and accurate enumeration of homologues in bacteria and archaea. A major difference was found between the distributions of enzyme families that are typically found in all three divisions of cellular life and non-enzymatic domain families that are usually eukaryote-specific. Previously undetected bacterial homologues were identified for# plant pathogenesis-related proteins, Pad1, von Willebrand factor type A, src homology 3 and YWTD repeat-containing domains. Comparisons of the domain distributions in eukaryotes and prokaryotes enabled distinctions to be made between the domains originating prior to the last common ancestor of all known life forms and those apparently originating as consequences of horizontal gene transfer events. A number of transfers of signalling domains from eukaryotes to bacteria were confidently identified, in contrast to only a single case of apparent transfer from eukaryotes to archaea.     

Deriving the genomic tree of life in the presence of horizontal gene transfer: conditioned reconstruction

The horizontal gene transfer (HGT) being inferred within prokaryotic genomes appears to be sufficiently massive that many scientists think it may have effectively obscured much of the history of life recorded in DNA. Here, we demonstrate that the tree of life can be reconstructed even in the presence of extensive HGT, provided the processes of genome evolution are properly modeled. We show that the dynamic deletions and insertions of genes that occur during genome evolution, including those introduced by HGT, may be modeled using techniques similar to those used to model nucleotide substitutions that occur during sequence evolution. In particular, we show that appropriately designed general Markov models are reasonable tools for reconstructing genome evolution. These studies indicate that, provided genomes contain sufficiently many genes and that the Markov assumptions are met, it is possible to reconstruct the tree of life. We also consider the fusion of genomes, a process not encountered in gene sequence evolution, and derive a method for the identification and reconstruction of genome fusion events. Genomic reconstructions of a well-defined classical four-genome problem, the root of the multicellular animals, show that the method, when used in conjunction with paralinear/logdet distances, performs remarkably well and is relatively unaffected by the recently discovered big genome artifact.     

Ancient phylogenetic relationships

Traditional views on deep evolutionary events have been seriously challenged over the last few years, following the identification of major pitfalls affecting molecular phylogeny reconstruction. Here we describe the principally encountered artifacts, notably long branch attraction, and their causes (i.e., difference in evolutionary rates, mutational saturation, compositional biases). Additional difficulties due to phenomena of biological nature (i.e., lateral gene transfer, recombination, hidden paralogy) are also discussed. Moreover, contrary to common beliefs, we show that the use of rare genomic events can also be misleading and should be treated with the same caution as standard molecular phylogeny. The universal tree of life, as described in most textbooks, is partly affected by tree reconstruction artifacts, e.g. (i) the bacterial rooting of the universal tree of life; (ii) the early emergence of amitochondriate lineages in eukaryotic phylogenies; and (iii) the position of hyperthermophilic taxa in bacterial phylogenies. We present an alternative view of this tree, based on recent evidence obtained from reanalyses of ancient data sets and from novel analyses of large combination of genes.     

Phylogenetic reconstruction and lateral gene transfer

Lateral gene transfer (LGT) is often seen as a form of noise, obscuring the phylogenetic signal with which we might hope to reconstruct the evolution of a group of organisms, or indeed the history of all life (the Tree of Life). Such reconstruction might still be possible if the subset of genes conserved among all genomes in a group (or common to all genomes) comprise a core that is relatively refractory to LGT. Several papers designed to test this notion have recently appeared, and here we re-analyze one, which claims that the core of single-copy orthologs shared by all sequenced genomes of the gammaproteobacteria is essentially free of LGT. This conclusion is unfortunately premature, and it is very hard to determine what fraction of this core has been affected by LGT. We discuss other difficulties with the core concept and suggest that, although the core idea must remain part of our understanding of phylogenetic relationships, it should not be the sole basis for defining such relationships, because these are not exclusively tree-like. We suggest instead a more complex but more natural framework for classification, which we call the Synthesis of Life.