Distance-Based Genome Rearrangement Phylogeny

Evolution operates on whole genomes through direct rearrangements of genes, such as inversions, transpositions, and inverted transpositions, as well as through operations, such as duplications, losses, and transfers, that also affect the gene content of the genomes. Because these events are rare relative to nucleotide substitutions, gene order data offer the possibility of resolving ancient branches in the tree of life; the combination of gene order data with sequence data also has the potential to provide more robust phylogenetic reconstructions, since each can elucidate evolution at different time scales. Distance corrections greatly improve the accuracy of phylogeny reconstructions from DNA sequences, enabling distance-based methods to approach the accuracy of the more elaborate methods based on parsimony or likelihood at a fraction of the computational cost. This paper focuses on developing distance correction methods for phylogeny reconstruction from whole genomes. The main question we investigate is how to estimate evolutionary histories from whole genomes with equal gene content, and we present a technique, the empirically derived estimator (EDE), that we have developed for this purpose. We study the use of EDE on whole genomes with identical gene content, and we explore the accuracy of phylogenies inferred using EDE with the neighbor joining and minimum evolution methods under a wide range of model conditions. Our study shows that tree reconstruction under these two methods is much more accurate when based on EDE distances than when based on other distances previously suggested for whole genomes. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Martin Kreitman]     

Molecular evolution in Hawaiian drosophilids

The evolutionary relationships of the Hawaiian drosophilids have been studied by several molecular techniques. Immunological measurements have determined some of the distances between the major groups. More detailed studies have used protein polymorphism, DNA hybridization, DNA restriction enzyme analysis and DNA sequence studies of both nuclear and mitochondrial genomes. These studies indicate both the strengths and the weaknesses of the various methods and suggest new ways to examine phylogenetic relationships and question some of the accepted ph ylogenies.     

Molecular homology and DNA hybridization

Summary We reviewed the concept of homology, which can broadly be defined as a correspondence between characteristics that is caused by continuity of information (Van Valen 1982). The concept applies widely in molecular biology when correspondence is taken to mean a genetic relationship resulting from a unique heritable modification of a feature at some previous point in time. Such correspodence can be established for features within a single organism as well as between organisms, making paralogy a valid form of molecular homology under this definition. Molecular homology can be recognized at a variety of organizational levels, which are intedependent. For example, the recognition of homology at the site level involves a statement of homology at the sequence level, and vice versa. This hierarchy, the potential for nonhomologous identity at the site level, and such processes as sequence transposition combine to yield a molecular equivalent to complex structural homology at the anatomical level. As a result, statements of homology between heritable units can involve a valid sense of percent homology.We analyzed DNA hybridization with respect to the problems of recognizing homology and using it in phylogenetic inference. Under a model requiring continuous divergence among compared sequences, DNA hybridization distances embed evolutionary hierarchy, and groups inferred using pairwise methods of tree reconstruction are based on underlying patterns of apomorphic homology. Thus, symplesiomorphic homology will not confound DNA hybridization phylogenies. However, nonhomologous identities that act like apomorphic homologies can lead to inaccurate reconstructions. The main difference between methods of phylogenetic analysis of DNA sequences is that parsimony methods permit hypotheses of nonhomology, whereas distance methods do not.This article was presented at the C.S.E.O.L. Conference on DNA-DNA Hybridization and Evolution, Lake Arrowhead, California, May 11–14, 1989     

Sequencing of whole plastid genomes and nuclear ribosomal DNA of Diospyros species (Ebenaceae) endemic to New Caledonia: many species,little divergence

Background and Aims Some plant groups, especially on islands, have been shaped by strong ancestral bottlenecks and rapid, recent radiation of phenotypic characters. Single molecular markers are often not informative enough for phylogenetic reconstruction in such plant groups. Whole plastid genomes and nuclear ribosomal DNA (nrDNA) are viewed by many researchers as sources of information for phylogenetic reconstruction of groups in which expected levels of divergence in standard markers are low. Here we evaluate the usefulness of these data types to resolve phylogenetic relationships among closely related Diospyros species.Methods Twenty-two closely related Diospyros species from New Caledonia were investigated using whole plastid genomes and nrDNA data from low-coverage next-generation sequencing (NGS). Phylogenetic trees were inferred using maximum parsimony, maximum likelihood and Bayesian inference on separate plastid and nrDNA and combined matrices.Key Results The plastid and nrDNA sequences were, singly and together, unable to provide well supported phylogenetic relationships among the closely related New Caledonian Diospyros species. In the nrDNA, a 6-fold greater percentage of parsimony-informative characters compared with plastid DNA was found, but the total number of informative sites was greater for the much larger plastid DNA genomes. Combining the plastid and nuclear data improved resolution. Plastid results showed a trend towards geographical clustering of accessions rather than following taxonomic species.Conclusions In plant groups in which multiple plastid markers are not sufficiently informative, an investigation at the level of the entire plastid genome may also not be sufficient for detailed phylogenetic reconstruction. Sequencing of complete plastid genomes and nrDNA repeats seems to clarify some relationships among the New Caledonian Diospyros species, but the higher percentage of parsimony-informative characters in nrDNA compared with plastid DNA did not help to resolve the phylogenetic tree because the total number of variable sites was much lower than in the entire plastid genome. The geographical clustering of the individuals against a background of overall low sequence divergence could indicate transfer of plastid genomes due to hybridization and introgression following secondary contact.     

Phylogenetics beyond biology

Evolutionary processes have been described not only in biology but also for a wide range of human cultural activities including languages and law. In contrast to the evolution of DNA or protein sequences, the detailed mechanisms giving rise to the observed evolution-like processes are not or only partially known. The absence of a mechanistic model of evolution implies that it remains unknown how the distances between different taxa have to be quantified. Considering distortions of metric distances, we first show that poor choices of the distance measure can lead to incorrect phylogenetic trees. Based on the well-known fact that phylogenetic inference requires additive metrics, we then show that the correct phylogeny can be computed from a distance matrix \({\mathbf {D}}\) if there is a monotonic, subadditive function \(\zeta\) such that \(\zeta ^{-1}({\mathbf {D}})\) is additive. The required metric-preserving transformation \(\zeta\) can be computed as the solution of an optimization problem. This result shows that the problem of phylogeny reconstruction is well defined even if a detailed mechanistic model of the evolutionary process remains elusive.     

A fast algorithm for the multiple genome rearrangement problem with weighted reversals and transpositions

Background  

Due to recent progress in genome sequencing, more and more data for phylogenetic reconstruction based on rearrangement distances between genomes become available. However, this phylogenetic reconstruction is a very challenging task. For the most simple distance measures (the breakpoint distance and the reversal distance), the problem is NP-hard even if one considers only three genomes.     

Reconstruction of Certain Phylogenetic Networks from Their Tree-Average Distances

Trees are commonly utilized to describe the evolutionary history of a collection of biological species, in which case the trees are called phylogenetic trees. Often these are reconstructed from data by making use of distances between extant species corresponding to the leaves of the tree. Because of increased recognition of the possibility of hybridization events, more attention is being given to the use of phylogenetic networks that are not necessarily trees. This paper describes the reconstruction of certain such networks from the tree-average distances between the leaves. For a certain class of phylogenetic networks, a polynomial-time method is presented to reconstruct the network from the tree-average distances. The method is proved to work if there is a single reticulation cycle.     

DNA evolutionary rates in nine-primaried passerine birds

Differences in single-copy nuclear-DNA sequences among 13 species of passerine birds were measured using DNA-DNA hybridization. A matrix of pairwise dissimilarity values (delta mode distances) was constructed from analysis of fitted thermal dissociation curves. A least-squares method of phylogenetic estimation was used to construct two topologies from the distance matrix, one constraining branch lengths of sister taxa to be equal and the other permitting such lengths to vary. These topologies were identical in the pattern of branching of taxa, and the difference in their sums of squares was not statistically significant, suggesting that rates of DNA evolution in sister groups of nine- primaried oscines are equal. A nonparametric test for nonrandom variation in distances of sister groups to outgroup taxa revealed no statistically significant deviation from random variation that would be expected as a result of measurement error. However, the level of measurement error was such that rates of DNA evolution in sister taxa could vary by as much as 10% without being detected with the statistical methods used here.      

Relative patterns and rates of evolution in heron nuclear and mitochondrial DNA

Mitochondrial cytochrome b sequence data from 15 species of herons (Aves: Ardeidae), representing 13 genera, were compared with DNA hybridization data of single-copy nuclear DNA (scnDNA) from the same species in a taxonomic congruence assessment of heron phylogeny. The two data sets produced a partially resolved, completely congruent estimate of phylogeny with the following basic structure: (Tigrisoma, Cochlearius, (((Zebrilus, (Ixobrychus, Botaurus)), (((Ardea, Casmerodius), Bubulcus), ((Egretta thula, Egretta caerulea, Egretta tricolor), Syrigma), Butorides, Nycticorax, Nyctanassa)))). Because congruence indicated similar phylogenetic information in the two data sets, we used the relatively unsaturated DNA hybridization distances as surrogates of time to examine graphically the patterns and rates of change in cytochrome b distances. Cytochrome b distances were computed either from whole sequences or from partitioned sequences consisting of transitions, transversions, specific codon site positions, or specific protein-coding regions. These graphical comparisons indicated that unpartitioned cytochrome b has evolved at 5-10 times the rate of scnDNA. Third-position transversions appeared to offer the most useful sequence partition for phylogenetic analysis because of their relatively fast rate of substitution (two times that of scnDNA) and negligible saturation. We also examined lineage-based rates of evolution by comparing branch length patterns between the nuclear and cytochrome b trees. The degree of correlation in corresponding branch lengths between cytochrome b and DNA hybridization trees depended on DNA sequence partitioning. When cytochrome b sequences were not partitioned, branch lengths in the cytochrome b and DNA hybridization trees were not correlated. However, when cytochrome b sequences were reduced to third-position transversions (i.e., unsaturated, relatively fast changing data), branch lengths were correlated. This finding suggests that lineage-based rates of DNA evolution in nuclear and mitochondrial genomes are influenced by common causes.     

Distorted metrics on trees and phylogenetic forests

We study distorted metrics on binary trees in the context of phylogenetic reconstruction. Given a binary tree T on n leaves with a path metric d, consider the pairwise distances {d(u,v)} between leaves. It is well known that these determine the tree and the d length of all edges. Here, we consider distortions d of d such that, for all leaves u and v, it holds that |d(u,v)-dmacr(u,v)|1.....T₀ such that the true tree T may be obtained from that forest by adding alpha-1 edges and alpha-1les2^-Omega(M/g)n. Our distorted metric result implies a reconstruction algorithm of phylogenetic forests with a small number of trees from sequences of length logarithmic in the number of species. The reconstruction algorithm is applicable for the general Markov model. Both the distorted metric result and its applications to phylogeny are almost tight     

BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments

Background  

The quality of multiple sequence alignments plays an important role in the accuracy of phylogenetic inference. It has been shown that removing ambiguously aligned regions, but also other sources of bias such as highly variable (saturated) characters, can improve the overall performance of many phylogenetic reconstruction methods. A current scientific trend is to build phylogenetic trees from a large number of sequence datasets (semi-)automatically extracted from numerous complete genomes. Because these approaches do not allow a precise manual curation of each dataset, there exists a real need for efficient bioinformatic tools dedicated to this alignment character trimming step.     

DNA hybridization as a guide to phylogeny: Chemical and physical limits

Summary The technique of forming interspecific DNA heteroduplexes and estimating phylogenetic distances from the depression in their duplex melting temperature has several physical and chemical constraints. These constraints determine the maximum phylogenetic distance that may be estimated by this technique and the most appropriate method of analyzing that distance.Melting curves of self-renatured single copy primate DNAs reveal the presence of components absent from the renaturation products of exactly paired sequences. This observation, which confirms existing literature, challenges a fundamental assumption: that orthologous (i.e., corresponding) DNA sequences in the divergent species are being compared in DNA heteroduplex melting experiments.As a model system, the thermal stabilities of heteroduplexes formed between a human alpha-globin cDNA and four alpha-like globin genes isolated from chimpanzee are qualitatively compared. The results of this comparison show that the cross-hybrids of imperfectly matched gene duplicates from divergent species can contribute to the additional components that are present in renatured single copy DNAs. Single copy DNA, as usually defined, includes sequence duplicates that will obscure phylogenetic comparisons in a mass hybridization of genomes.     

DISCORDANCE BETWEEN NUCLEAR AND CHLOROPLAST PHYLOGENIES IN THE HEUCHERA GROUP (SAXIFRAGACEAE)

Various factors, including taxon density, sampling error, convergence, and heterogeneity of evolutionary rates, can potentially lead to incongruence between phylogenetic trees based on different genomes. Particularly at the generic level and below, chloroplast capture resulting from hybridization may distort organismal relationships in phylogenetic analyses based on the chloroplast genome, or genes included therein. However, the extent of such discord between chloroplast DNA (cpDNA) trees and those trees based on nuclear genes has rarely been assessed. We therefore used sequences of the internal transcribed spacer regions (ITS-1 and ITS-2) of nuclear ribosomal DNA (rDNA) to reconstruct phylogenetic relationships among members of the Heuchera group of genera (Saxifragaceae). The Heuchera group presents an important model for the analysis of chloroplast capture and its impact on phylogenetic reconstruction because hybridization is well documented within genera (e.g., Heuchera), and intergeneric hybrids involving six of the nine genera have been reported. An earlier study provided a well-resolved phylogenetic hypothesis for the Heuchera group based on cpDNA restriction-site variation. However, trees based on ITS sequences are discordant with the cpDNA-based tree. Evidence from both morphology and nuclear-encoded allozymes is consistent with the ITS trees, rather than the cpDNA tree, and several points of phylogenetic discord can clearly be attributed to chloroplast capture. Comparison of the organellar and ITS trees also raises the strong likelihood that ancient events of chloroplast capture occurred between lineages during the early diversification of the Heuchera group. Thus, despite the many advantages and widespread use of cpDNA data in phylogeny reconstruction, comparison of relationships based on cpDNA and ITS sequences for the Heuchera group underscores the need for caution in the use of organellar variation for retrieving phylogeny at lower taxonomic levels, particularly in groups noted for hybridization.     

Phylogenomics of Fargesia and Yushania reveals a history of reticulate evolution

Reticulate evolution is a common and important driving force in angiosperm evolution. In this study, we analyzed the phylogenetic signals of genomic regions with different inheritance patterns to understand the evolutionary process of organisms using species-rich Himalaya–Hengduan taxa of bamboos (Fargesia Franchet and Yushania Keng). We constructed phylogenetic trees using different sampling strategies and reconstruction methods based on genome skimming and double digest restriction-site-associated DNA sequencing data. We assessed the congruence of topologies generated from different datasets and employed several approaches to reveal the causes of phylogenetic incongruence, including the detection of hybridization and introgression using PhyloNetworks and the D-statistic test (ABBA-BABA test). We found that, in the plastome-based phylogeny, Fargesia bamboos can be clustered into three groups and Yushania was nested within one of them, which contradicts the nuclear–double digest restriction-site-associated DNA sequencing-based phylogeny. Moreover, the genetic variation of chloroplast DNA is significantly correlated with geographical distribution. The strong signal of incomplete lineage sorting, hybridization, introgression, and cytoplasmic gene flow found among genera and species suggests that reticulate evolution is the main cause for the phylogenetic incongruence between nuclear and chloroplast datasets. Our results add evidence that genomes with different inheritance patterns can reveal distinct evolutionary histories of species and suggest that reticulate evolution is prevalent in rapidly diversifying groups.     

Genome BLAST distance phylogenies inferred from whole plastid and whole mitochondrion genome sequences

Background  

Phylogenetic methods which do not rely on multiple sequence alignments are important tools in inferring trees directly from completely sequenced genomes. Here, we extend the recently described Genome BLAST Distance Phylogeny (GBDP) strategy to compute phylogenetic trees from all completely sequenced plastid genomes currently available and from a selection of mitochondrial genomes representing the major eukaryotic lineages. BLASTN, TBLASTX, or combinations of both are used to locate high-scoring segment pairs (HSPs) between two sequences from which pairwise similarities and distances are computed in different ways resulting in a total of 96 GBDP variants. The suitability of these distance formulae for phylogeny reconstruction is directly estimated by computing a recently described measure of "treelikeness", the so-called δ value, from the respective distance matrices. Additionally, we compare the trees inferred from these matrices using UPGMA, NJ, BIONJ, FastME, or STC, respectively, with the NCBI taxonomy tree of the taxa under study.     

Seven new dolphin mitochondrial genomes and a time-calibrated phylogeny of whales

Background  

The phylogeny of Cetacea (whales) is not fully resolved with substantial support. The ambiguous and conflicting results of multiple phylogenetic studies may be the result of the use of too little data, phylogenetic methods that do not adequately capture the complex nature of DNA evolution, or both. In addition, there is also evidence that the generic taxonomy of Delphinidae (dolphins) underestimates its diversity. To remedy these problems, we sequenced the complete mitochondrial genomes of seven dolphins and analyzed these data with partitioned Bayesian analyses. Moreover, we incorporate a newly-developed "relaxed" molecular clock to model heterogenous rates of evolution among cetacean lineages.     

Fast and accurate phylogenetic reconstruction from high-resolution whole-genome data and a novel robustness estimator

The rapid accumulation of whole-genome data has renewed interest in the study of genomic rearrangements. Comparative genomics, evolutionary biology, and cancer research all require models and algorithms to elucidate the mechanisms, history, and consequences of these rearrangements. However, even simple models lead to NP-hard problems, particularly in the area of phylogenetic analysis. Current approaches are limited to small collections of genomes and low-resolution data (typically a few hundred syntenic blocks). Moreover, whereas phylogenetic analyses from sequence data are deemed incomplete unless bootstrapping scores (a measure of confidence) are given for each tree edge, no equivalent to bootstrapping exists for rearrangement-based phylogenetic analysis. We describe a fast and accurate algorithm for rearrangement analysis that scales up, in both time and accuracy, to modern high-resolution genomic data. We also describe a novel approach to estimate the robustness of results-an equivalent to the bootstrapping analysis used in sequence-based phylogenetic reconstruction. We present the results of extensive testing on both simulated and real data showing that our algorithm returns very accurate results, while scaling linearly with the size of the genomes and cubically with their number. We also present extensive experimental results showing that our approach to robustness testing provides excellent estimates of confidence, which, moreover, can be tuned to trade off thresholds between false positives and false negatives. Together, these two novel approaches enable us to attack heretofore intractable problems, such as phylogenetic inference for high-resolution vertebrate genomes, as we demonstrate on a set of six vertebrate genomes with 8,380 syntenic blocks. A copy of the software is available on demand.     

DNA Hybridization as a Guide to Phylogenies: a Critical Analysis

Abstract— This article evaluates the use of DNA hybridization for estimating the extent of divergence among the single-copy fractions of vertebrate genomes. It focuses, in particular, on the nature and informational content of the melting profiles as a guide to phylogenetic relationships. While concluding that the DNA hybridization approach remains the best and most cost-effective guide to such relationships over its useful range, it demonstrates serious flaws in certain recent attempts to apply the method to specific cases among primates and birds. The major points are:

• 1 The T₅₀H statistic is flawed as a measure of mean sequence divergence, and also, therefore, as a measure of phylogenetic distance.
• 2 The T_mode statistic overcomes many of the problems inherent in interpreting thermal stabilities of DNA heteroduplexes for phylogenetic purposes.
• 3 The phylogenetic significance of ΔT_modes of > 15d? or so cannot be accurately assessed.
• 4 The putative slowdown in the rate of nuclear DNA sequence change among the lemurs is not justified by the data.
• 5 The claims of Sibley and Ahlquist to have resolved the human/chimpanzee/gorilla trichotomy are not supported by their data.
• 6 There are major problems in the published Sibley and Ahlquist avian phylogenies; in particular, with those containing evolutionary “staircases” of nodes separated by less than 1d? from one another.
• 7 There would appear to be a lineage misplacement involving a ΔT of at least 4d? in a recent publication on avian phylogeny.
• 8 Certain of the published ΔT₅₀ H values seem not to be representative of the actual data on which they are based.
• 9 Most important, it is recommended that no phytogenies based on DNA hybridization comparisons should be presented without being accompanied by the data relevant to each claim of a resolved lineage.

     

Recombination in Mitochondrial DNA of European Mussels <Emphasis Type="Italic">Mytilus</Emphasis>

Mitochondrial DNA was long believed to be purely clonal and free from recombination. Major phylogenetic studies still depend on such assumptions. The peculiar genetic system of marine mussels Mytilus in which two divergent mitochondrial genomes exist provides a unique opportunity to study mtDNA recombination. Previous reports showed the existence of a few haplotypes having very strong recombination signal in the control region of mtDNA. Those recombinant variants have been found in a Baltic Sea population of Mytilus trossulus as well as in Mytilus galloprovincialis from the Black Sea. In both cases the mosaic genomes switched their transmission route and have been inherited paternally. In the present study rearranged mtDNA genomes found in all three European Mytilus species are described. The structure of their control region is a result of intra- and intermolecular recombination between mitochondrial genomes. Together with the phylogenetic reconstruction and geographic distribution, this suggests that two interlineage recombination events have occurred in the control region of mtDNA of European mussels Mytilus. Contrary to earlier observations, some of the mosaic genomes do not show any gender bias, which has important implications regarding the transmission and evolution of blue mussel mitochondrial genomes.