Quartet decomposition server: a platform for analyzing phylogenetic trees

Evaluating Compression-Based Phylogeny Estimation in the Presence of Incomplete Lineage Sorting

This study assesses characteristics of the normalized compression distance (NCD) technique for building phylogenetic trees from molecular data. We examined results from a mammalian biological data set as well as a collection of simulated data with varying levels of incomplete lineage sorting. The implementation of NCD we analyze is a concatenation-based, distance-based, alignment-free, and model-free phylogeny estimation method, which takes concatenated unaligned sequence data as input and outputs a matrix of distances. We compare the NCD phylogeny estimation method with various other methods, including coalescent- and concatenation-based methods.

A New Quartet-Based Statistical Method for Comparing Sets of Gene Trees Is Developed Using a Generalized Hoeffding Inequality

Extracting the strength of the tree signal that is encompassed by a collection of gene trees is an exceptionally challenging problem in phylogenomics. Often, this problem not only involves the construction of individual phylogenies based on different genes, which may be a difficult endeavor on its own, but is also exacerbated by many factors that create conflicts between the evolutionary histories of different gene families, such as duplications or losses of genes; hybridization events; incomplete lineage sorting; and horizontal gene transfer, the latter two play central roles in the evolution of eukaryotes and prokaryotes, respectively. In this work, we tackle the aforementioned problem by focusing on quartet trees, which are the most basic unit of information in the context of unrooted phylogenies. In the first part, we show how a theorem of Janson that generalizes the classical Hoeffding inequality can be used to develop a statistical test involving quartets. In the second part, we study real and simulated data using this theoretical advancement, thus demonstrating how the significance of the differences between sets of quartets can be assessed. Our results are particularly intriguing since they nonstandardly require the analysis of dependent random variables.

Efficient Quartet Representations of Trees and Applications to Supertree and Summary Methods

Quartet trees displayed by larger phylogenetic trees have long been used as inputs for species tree and supertree reconstruction. Computational constraints prevent the use of all displayed quartets in many practical problems with large numbers of taxa. We introduce the notion of an Efficient Quartet System (EQS) to represent a phylogenetic tree with a subset of the quartets displayed by the tree. We show mathematically that the set of quartets obtained from a tree via an EQS contains all of the combinatorial information of the tree itself. Using performance tests on simulated datasets, we also demonstrate that using an EQS to reduce the number of quartets in both summary method pipelines for species tree inference as well as methods for supertree inference results in only small reductions in accuracy.

Genomic and phenotypic differentiation among Methanosarcina mazei populations from Columbia River sediment

Methanogenic archaea are genotypically and phenotypically diverse organisms that are integral to carbon cycling in anaerobic environments. Owing to their genetic tractability and ability to be readily cultivated, Methanosarcina spp. have become a powerful model system for understanding methanogen biology at the cellular systems level. However, relatively little is known of how genotypic and phenotypic variation is partitioned in Methanosarcina populations inhabiting natural environments and the possible ecological and evolutionary implications of such variation. Here, we have identified how genomic and phenotypic diversity is partitioned within and between Methanosarcina mazei populations obtained from two different sediment environments in the Columbia River Estuary (Oregon, USA). Population genomic analysis of 56 M. mazei isolates averaging <1% nucleotide divergence revealed two distinct clades, which we refer to as 'mazei-T' and 'mazei-WC'. Genomic analyses showed that these clades differed in gene content and fixation of allelic variants, which point to potential differences in primary metabolism and also interactions with foreign genetic elements. This hypothesis of niche partitioning was supported by laboratory growth experiments that revealed significant differences in trimethylamine utilization. These findings improve our understanding of the ecologically relevant scales of genomic variation in natural systems and demonstrate interactions between genetic and ecological diversity in these easily cultivable and genetically tractable model methanogens.

Evidence for extensive gene flow and Thermotoga subpopulations in subsurface and marine environments

Oil reservoirs represent a nutrient-rich ecological niche of the deep biosphere. Although most oil reservoirs are occupied by microbial populations, when and how the microbes colonized these environments remains unanswered. To address this question, we compared 11 genomes of Thermotoga maritima-like hyperthermophilic bacteria from two environment types: subsurface oil reservoirs in the North Sea and Japan, and marine sites located in the Kuril Islands, Italy and the Azores. We complemented our genomes with Thermotoga DNA from publicly available subsurface metagenomes from North America and Australia. Our analysis revealed complex non-bifurcating evolutionary history of the isolates' genomes, suggesting high amounts of gene flow across all sampled locations, a conjecture supported by numerous recombination events. Genomes from the same type of environment tend to be more similar, and have exchanged more genes with each other than with geographically close isolates from different types of environments. Hence, Thermotoga populations of oil reservoirs do not appear isolated, a requirement of the 'burial and isolation' hypothesis, under which reservoir bacteria are descendants of the isolated communities buried with sediments that over time became oil reservoirs. Instead, our analysis supports a more complex view, where bacteria from subsurface and marine populations have been continuously migrating into the oil reservoirs and influencing their genetic composition. The Thermotoga spp. in the oil reservoirs in the North Sea and Japan probably entered the reservoirs shortly after they were formed. An Australian oil reservoir, on the other hand, was likely colonized very recently, perhaps during human reservoir development.

Quantifying homologous replacement of loci between haloarchaeal species

In vitro studies of the haloarchaeal genus Haloferax have demonstrated their ability to frequently exchange DNA between species, whereas rates of homologous recombination estimated from natural populations in the genus Halorubrum are high enough to maintain random association of alleles between five loci. To quantify the effects of gene transfer and recombination of commonly held (relaxed core) genes during the evolution of the class Halobacteria (haloarchaea), we reconstructed the history of 21 genomes representing all major groups. Using a novel algorithm and a concatenated ribosomal protein phylogeny as a reference, we created a directed horizontal genetic transfer (HGT) network of contemporary and ancestral genomes. Gene order analysis revealed that 90% of testable HGTs were by direct homologous replacement, rather than nonhomologous integration followed by a loss. Network analysis revealed an inverse log-linear relationship between HGT frequency and ribosomal protein evolutionary distance that is maintained across the deepest divergences in Halobacteria. We use this mathematical relationship to estimate the total transfers and amino acid substitutions delivered by HGTs in each genome, providing a measure of chimerism. For the relaxed core genes of each genome, we conservatively estimate that 11–20% of their evolution occurred in other haloarchaea. Our findings are unexpected, because the transfer and homologous recombination of relaxed core genes between members of the class Halobacteria disrupts the coevolution of genes; however, the generation of new combinations of divergent but functionally related genes may lead to adaptive phenotypes not available through cumulative mutations and recombination within a single population.