Search

Model-based molecular phylogenetics plays an important role in comparisons of genomic data, and model selection is a key step in all such analyses. We present ModelFinder, a fast model-selection method that greatly improves the accuracy of phylogenetic estimates by incorporating a model of rate heterogeneity across sites not previously considered in this context and by allowing concurrent searches of model space and tree space.

Insects are the most speciose group of animals, but the phylogenetic relationships of many major lineages remain unresolved. We inferred the phylogeny of insects from 1478 protein-coding genes. Phylogenomic analyses of nucleotide and amino acid sequences, with site-specific nucleotide or domain-specific amino acid substitution models, produced statistically robust and congruent results resolving previously controversial phylogenetic relations hips. We dated the origin of insects to the Early Ordovician [~479 million years ago (Ma)], of insect flight to the Early Devonian (~406 Ma), of major extant lineages to the Mississippian (~345 Ma), and the major diversification of holometabolous insects to the Early Cretaceous. Our phylogenomic study provides a comprehensive reliable scaffold for future comparative analyses of evolutionary innovations among insects., (Copyright © 2014, American Association for the Advancement of Science.)

Background: Under a Markov model of evolution, recoding, or lumping, of the four nucleotides into fewer groups may permit analysis under simpler conditions but may unfortunately yield misleading results unless the evolutionary process of the recoded groups remains Markovian. If a Markov process is lumpable, then the evolutionary process of the recoded groups is Markovian., Results: We consider stationary, reversible, and homogeneous Markov processes on two taxa and compare three tests for lumpability: one using an ad hoc test statistic, which is based on an index that is evaluated using a bootstrap approximation of its distribution; one that is based on a test proposed specifically for Markov chains; and one using a likelihood-ratio test. We show that the likelihood-ratio test is more powerful than the index test, which is more powerful than that based on the Markov chain test statistic. We also show that for stationary processes on binary trees with more than two taxa, the tests can be applied to all pairs. Finally, we show that if the process is lumpable, then estimates obtained under the recoded model agree with estimates obtained under the original model, whereas, if the process is not lumpable, then these estimates can differ substantially. We apply the new likelihood-ratio test for lumpability to two primate data sets, one with a mitochondrial origin and one with a nuclear origin., Conclusions: Recoding may result in biased phylogenetic estimates because the original evolutionary process is not lumpable. Accordingly, testing for lumpability should be done prior to phylogenetic analysis of recoded data.

Although the status of Crinoidea (sea lilies and featherstars) as sister group to all other living echinoderms is well-established, relationships among crinoids, particularly extant forms, are debated. All living species are currently placed in Articulata, which is generally accepted as the only crinoid group to survive the Permian-Triassic extinction event. Recent classifications have recognized five major extant taxa: Isocrinida, Hyocrinida, Bourgueticrinina, Comatulidina and Cyrtocrinida, plus several smaller groups with uncertain taxonomic status, e.g., Guillecrinus, Proisocrinus and Caledonicrinus. Here we infer the phylogeny of extant Crinoidea using three mitochondrial genes and two nuclear genes from 59 crinoid terminals that span the majority of extant crinoid diversity. Although there is poor support for some of the more basal nodes, and some tree topologies varied with the data used and mode of analysis, we obtain several robust results. Cyrtocrinida, Hyocrinida, Isocrinida are all recovered as clades, but two stalked crinoid groups, Bourgueticrinina and Guillecrinina, nest among the featherstars, lending support to an argument that they are paedomorphic forms. Hence, they are reduced to families within Comatulida. Proisocrinus is clearly shown to be part of Isocrinida, and Caledonicrinus may not be a bourgueticrinid. Among comatulids, tree topologies show little congruence with current taxonomy, indicating that much systematic revision is required. Relaxed molecular clock analyses with eight fossil calibration points recover Articulata with a median date to the most recent common ancestor at 231-252mya in the Middle to Upper Triassic. These analyses tend to support the hypothesis that the group is a radiation from a small clade that passed through the Permian-Triassic extinction event rather than several lineages that survived. Our tree topologies show various scenarios for the evolution of stalks and cirri in Articulata, so it is clear that further data and taxon sampling are needed to recover a more robust phylogeny of the group., (Copyright © 2012 Elsevier Inc. All rights reserved.)

The selection of an optimal model for data analysis is an important component of model-based molecular phylogenetic studies. Owing to the large number of Markov models that can be used for data analysis, model selection is a combinatorial problem that cannot be solved by performing an exhaustive search of all possible models. Currently, model selection is based on a small subset of the available Markov models, namely those that assume the evolutionary process to be globally stationary, reversible, and homogeneous. This forces the optimal model to be time reversible even though the actual data may not satisfy these assumptions. This problem can be alleviated by including more complex models during the model selection. We present a novel heuristic that evaluates a small fraction of these complex models and identifies the optimal model.

Indoleamine 2,3-dioxygenase (IDO1) and tryptophan 2,3-dioxygenase (TDO) are tryptophan-degrading enzymes that catalyze the first step in tryptophan catabolism via the kynurenine pathway. TDO is widely distributed in both eukaryotes and bacteria. In contrast, IDO has been found only in mammals and yeast. In 2007, a third enzyme, indoleamine 2,3-dioxygenase-2 (IDO2), was discovered. IDO2 is found not only in mammals but also in lower vertebrates. Interestingly, the Km value of IDO2 for L-Trp was 500-1000 fold higher than that of IDO1. In this study, we isolated both IDO1 and IDO2 cDNA from a monotreme, the platypus (Ornithorhynchus anatinus), and a marsupial, the gray short-tailed opossum (Monodelphis domestica). We characterized the recombinant proteins and those of other known IDO1/IDO2 in intact cells and a cell-free system. It was found that methylene blue may not be suitable reductant for IDO2, hence resulting in an underestimation of recombinant IDO2 activity. In intact cells, the Km value of IDO2 for L-Trp was estimated to be much higher than that of IDO1 and this high Km value appears to have been conserved during the evolution of IDO2. The protein encoded by the ancestor gene of IDO1 and IDO2 is likely to have had properties more similar to present day IDO2 than to IDO1.

The delineation and phylogenetic relationships of the pulmonate family Helicarionidae are currently poorly understood. We have undertaken a phylogenetic analysis of Australian members of Helicarionidae. Three consecutive mitochondrial genes (COI, tRNA-val and 16S rRNA) were sequenced for 36 species from Helicarionidae and related groups. Helicarionidae grouped with Ariophantidae and Urocyclidae with good support in all trees, as did Microcystinae with Trochomorphidae, corresponding to the results of a previous morphological study, but the relationships of Euconulinae and Limacidae could not be resolved. Cystopeltidae never grouped with Helicarionidae and should be regarded as a separate family. The position of tRNA-val in the mitochondrial genome provided a new synapomorphy for Microcystinae. In addition, the COI and 16S rRNA sequences showed a high degree of compositional heterogeneity and incompatibility of phylogenetic signals, highlighting the importance of testing for the decay of the historical signal prior to the phylogenetic analysis.

Although recent studies indicate that estimating phylogenies from alignments of concatenated genes greatly reduces the stochastic error, the potential for systematic error still remains, heightening the need for reliable methods to analyze multigene data sets. Consensus methods provide an alternative, more inclusive, approach for analyzing collections of trees arising from multiple genes. We extend a previously described consensus network method for genome-scale phylogeny (Holland, B. R., K. T. Huber, V. Moulton, and P. J. Lockhart. 2004. Using consensus networks to visualize contradictory evidence for species phylogeny. Mol. Biol. Evol. 21:1459-1461) to incorporate additional information. This additional information could come from bootstrap analysis, Bayesian analysis, or various methods to find confidence sets of trees. The new methods can be extended to include edge weights representing genetic distance. We use three data sets to illustrate the approach: 61 genes from 14 angiosperm taxa and one gymnosperm, 106 genes from eight yeast taxa, and 46 members of a gene family from 15 vertebrate taxa.

Phylogenetic studies of dinoflagellates are often conducted using rDNA sequences. In analyses to date, the monophyly of some of the major lineages of dinoflagellates remain to be demonstrated. There are several reasons for this uncertainty, one of which may be the use of models of evolution that may not closely fit the data. We constructed and examined alignments of SSU and partial LSU rRNA along with a concatenated alignment of the two molecules. The alignments showed several characteristics that may confound phylogeny reconstruction: paired helix (stem) regions that contain non-independently evolving sites, high levels of compositional heterogeneity among some of the sequences, high levels of incompatibility (homoplasy), and rate heterogeneity among sites. Taking into account these confounding factors, we analysed the data and found that the Gonyaulacales, a well-supported clade, may be the most recently diverged order. Other supported orders were, in the analysis based on SSU, the Suessiales and the Dinophysiales; however, the Gymnodiniales and Prorocentrales appeared to be polyphyletic. The Peridiniales without Heterocapsa species appeared as a monophyletic group in the analysis based on LSU; however, the support was low. The concatenated alignment did not provide a better phylogenetic resolution than the single gene alignments.

Recent findings of new genes in fish related to the prion protein (PrP) gene PRNP, including our recent report of SPRN coding for Shadoo (Sho) protein found also in mammals, raise issues of their function and evolution. Here we report additional novel fish genes found in public databases, including a duplicated SPRN gene, SPRNB, in Fugu, Tetraodon, carp, and zebrafish encoding the Sho2 protein, and we use comparative genomic analysis to analyze the evolutionary relationships and to infer evolutionary trajectories of the complete data set. Phylogenetic footprinting performed on aligned human, mouse, and Fugu SPRN genes to define candidate regulatory promoter regions, detected 16 conserved motifs, three of which are known transcription factor-binding sites for a receptor and transcription factors specific to or associated with expression in brain. This result and other homology-based (VISTA global genomic alignment; protein sequence alignment and phylogenetics) and context-dependent (genomic context; relative gene order and orientation) criteria indicate fish and mammalian SPRN genes are orthologous and suggest a strongly conserved basic function in brain. Whereas tetrapod PRNPs share context with the analogous stPrP-2-coding gene in fish, their sequences are diverged, suggesting that the tetrapod and fish genes are likely to have significantly different functions. Phylogenetic analysis predicts the SPRN/SPRNB duplication occurred before divergence of fish from tetrapods, whereas that of stPrP-1 and stPrP-2 occurred in fish. Whereas Sho appears to have a conserved function in vertebrate brain, PrP seems to have an adaptive role fine-tuned in a lineage-specific fashion. An evolutionary model consistent with our findings and literature knowledge is proposed that has an ancestral prevertebrate SPRN-like gene leading to all vertebrate PrP-related and Sho-related genes. This provides a new framework for exploring the evolution of this unusual family of proteins and for searching for members in other fish branches and intermediate vertebrate groups.

We present a computer program to simulate the evolution of a nucleotide sequence on a phylogenetic tree with four tips. The program, Hetero, allows users to assign lineage-specific differences in the rate matrices used to describe the evolutionary process. It has a simple user interface and output, making it equally useful in the teaching and research of phylogenetics.

Intragenic duplication is an evolutionary process where segments of a gene become duplicated. While there has been much research into whole-gene or domain duplication, there have been very few studies of non-tandem intragenic duplication. The identification of intragenically replicated sequences may provide insight into the evolution of proteins, helping to link sequence data with structure and function. This paper describes a tool for autonomously modelling intragenic duplication. AMID provides: identification of modularly repetitive genes; an algorithm for identifying repeated modules; and a scoring system for evaluating the modules' similarity. An evaluation of the algorithms and use cases are presented.

Relationships among laurasiatherian clades represent one of the most highly disputed topics in mammalian phylogeny. In this study, we attempt to disentangle laurasiatherian interordinal relationships using two independent genome-level approaches: (1) quantifying retrotransposon presence/absence patterns, and (2) comparisons of exon datasets at the levels of nucleotides and amino acids. The two approaches revealed contradictory phylogenetic signals, possibly due to a high level of ancestral incomplete lineage sorting. The positions of Eulipotyphla and Chiroptera as the first and second earliest divergences were consistent across the approaches. However, the phylogenetic relationships of Perissodactyla, Cetartiodactyla, and Ferae, were contradictory. While retrotransposon insertion analyses suggest a clade with Cetartiodactyla and Ferae, the exon dataset favoured Cetartiodactyla and Perissodactyla. Future analyses of hitherto unsampled laurasiatherian lineages and synergistic analyses of retrotransposon insertions, exon and conserved intron/intergenic sequences might unravel the conflicting patterns of relationships in this major mammalian clade.

The general Markov model (GMM) of nucleotide substitution does not assume the evolutionary process to be stationary, reversible, or homogeneous. The GMM can be simplified by assuming the evolutionary process to be stationary. A stationary GMM is appropriate for analyses of phylogenetic data sets that are compositionally homogeneous; a data set is considered to be compositionally homogeneous if a statistical test does not detect significant differences in the marginal distributions of the sequences. Though the general time-reversible (GTR) model assumes stationarity, it also assumes reversibility and homogeneity. We propose two new stationary and nonhomogeneous models—one constrains the GMM to be reversible, whereas the other does not. The two models, coupled with the GTR model, comprise a set of nested models that can be used to test the assumptions of reversibility and homogeneity for stationary processes. The two models are extended to incorporate invariable sites and used to analyze a seven-taxon hominoid data set that displays compositional homogeneity. We show that within the class of stationary models, a nonhomogeneous model fits the hominoid data better than the GTR model. We note that if one considers a wider set of models that are not constrained to be stationary, then an even better fit can be obtained for the hominoid data. However, the methods for reducing model complexity from an extremely large set of nonstationary models are yet to be developed. [ABSTRACT FROM PUBLISHER]

Two hypotheses account for the evolution of the inner antenna light-harvesting proteins of oxygenic photosynthesis in cyanobacteria, algae, and plants: one in which the CP43 protein of photosytem II gave rise to the extrinsic CP43-like antennas of cyanobacteria (i.e. IsiA and Pcb proteins), as a late development, and the other in which CP43 and CP43-like proteins derive from an ancestral protein. In order to determine which of these hypotheses is most likely, we analyzed the family of antenna proteins by a variety of phylogenetic techniques, using alignments of the six common membrane-spanning helices, constructed using information on the antenna proteins’ three-dimensional structure, and surveyed for evidence of factors that might confound inference of a correct phylogeny. The first hypothesis was strongly supported. As a consequence, we conclude that the ancestral photosynthetic apparatus, with 11 membrane-spanning helices, split at an early stage during evolution to form, on the one hand, the reaction center of photosystem II and, on the other hand, the ancestor of inner antenna proteins, CP43 (PsbC) and CP47 (PsbB). Only much later in evolution did the CP43 lineage give rise to the CP43’ proteins (IsiA and Pcb) of cyanobacteria. [ABSTRACT FROM AUTHOR]

Alignments of nucleotide or amino acid sequences may contain a variety of different signals, one of which is the historical signal that we often try to recover by phylogenetic analysis. Other signals, such as those arising due to compositional heterogeneities, among-lineage and among-site rate heterogeneities, invariant sites, and covariotides, may interfere adversely with the recovery of the historical signal. The effect of the interaction of these signals on phylogenetic inference is not well understood and may, in many cases, even be underappreciated. In this study, we investigate this matter and present results based on Monte Carlo simulations. We explored the success of four phylogenetic methods in recovering the true tree from data that had evolved under conditions where the equilibrium base frequencies and substitution rates were allowed to vary among lineages. Seven scenarios with increasingly complex conditions were investigated. All of the methods tested, with the exception of neighbor-joining using LogDet distances, were sensitive to compositional convergence in nonsister lineages. Maximum parsimony was also susceptible to attraction between long edges. In many cases, however, phylogenetic inference methods can still recover the true tree when misleading signals are present, in some instances even when the historical signal is no longer dominant. These results highlight the growing need for simple methods to detect violation of the phylogenetic assumptions. [Compositional heterogeneity; edge lengths; Monte Carlo simulation; networks; phylogenetic signal; rate heterogeneity; substitutional saturation.] [ABSTRACT FROM AUTHOR]

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content. [ABSTRACT FROM AUTHOR]