Heterotachy and tree building: a case study with plastids and eubacteria

The nature of heterotachy at the center of recent controversy over the relative performance of tree-building methods is different from the form of heterotachy that has been inferred in empirical studies. The latter have suggested that proportions of variable sites (p(var)) vary among orthologues and among paralogues. However, the strength of this inference, describing what may be one of the most important evolutionary properties of sequence data, has remained weak. Consequently, other models of sequence evolution have been proposed to explain some long-branch attraction (LBA) problems that could be attributed to differences in p(var). For an empirical case with plastid and eubacterial RNA polymerase sequences, we confirm using capture-recapture estimates and simulations that p(var) can differ among orthologues in anciently diverged evolutionary lineages. We find that parsimony and a least squares distance method that implements an overly simple model of sequence evolution are susceptible to LBA induced by this form of heterotachy. Although homogeneous maximum likelihood inference was found to be robust to model misspecification in our specific example, we caution against assuming that it will always be so.     

Phylogenomic analyses support the monophyly of Taphrinomycotina, including Schizosaccharomyces fission yeasts

Several morphologically dissimilar ascomycete fungi includingSchizosaccharomyces, Taphrina, Saitoella, Pneumocystis, andNeolecta have been grouped into the taxon Taphrinomycotina (Archiascomycotaor Archiascomycotina), originally based on rRNA phylogeny. Theseanalyses lack statistically significant support for the monophylyof this grouping, and although confirmed by more recent multigeneanalyses, this topology is contradicted by mitochondrial phylogenies.To resolve this inconsistency, we have assembled phylogenomicmitochondrial and nuclear data sets from four distantly relatedtaphrinomycotina taxa: Schizosaccharomyces pombe, Pneumocystiscarinii, Saitoella complicata, and Taphrina deformans. Our phylogenomicanalyses based on nuclear data (113 proteins) conclusively supportthe monophyly of Taphrinomycotina, diverging as a sister groupto Saccharomycotina + Pezizomycotina. However, despite the improvedtaxon sampling, Taphrinomycotina continue to be paraphyleticwith the mitochondrial data set (13 proteins): Schizosaccharomycesspecies associate with budding yeasts (Saccharomycotina) andthe other Taphrinomycotina group as a sister group to Saccharomycotina+ Pezizomycotina. Yet, as Schizosaccharomyces and Saccharomycotinaspecies are fast evolving, the mitochondrial phylogeny may beinfluenced by a long-branch attraction (LBA) artifact. Afterremoval of fast-evolving sequence positions from the mitochondrialdata set, we recover the monophyly of Taphrinomycotina. Ourcombined results suggest that Taphrinomycotina is a legitimatetaxon, that this group of species diverges as a sister groupto Saccharomycotina + Pezizomycotina, and that phylogeneticpositioning of yeasts and fission yeasts with mitochondrialdata is plagued by a strong LBA artifact.     

Evaluating the performance of a successive-approximations approach to parameter optimization in maximum-likelihood phylogeny estimation

Almost all studies that estimate phylogenies from DNA sequencedata under the maximum-likelihood (ML) criterion employ an approximateapproach. Most commonly, model parameters are estimated on someinitial phylogenetic estimate derived using a rapid method (neighbor-joiningor parsimony). Parameters are then held constant during a treesearch, and ideally, the procedure is repeated until convergenceis achieved. However, the effectiveness of this approximationhas not been formally assessed, in part because doing so requirescomputationally intensive, full-optimization analyses. Here,we report both indirect and direct evaluations of the effectivenessof successive approximations. We obtained an indirect evaluationby comparing the results of replicate runs on real data thatuse random trees to provide initial parameter estimates. Forsix real data sets taken from the literature, all replicateiterative searches converged to the same joint estimates oftopology and model parameters, suggesting that the approximationis not starting-point dependent, as long as the heuristic searchesof tree space are rigorous. We conducted a more direct assessmentusing simulations in which we compared the accuracy of phylogeniesestimated using full optimization of all model parameters oneach tree evaluated to the accuracy of trees estimated via successiveapproximations. There is no significant difference between theaccuracy of the approximation searches relative to full-optimizationsearches. Our results demonstrate that successive approximationis reliable and provide reassurance that this much faster approachis safe to use for ML estimation of topology.     

Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia

Almost a decade ago, a new phylogeny of bilaterian animals was inferred from small-subunit ribosomal RNA (rRNA) that claimed the monophyly of two major groups of protostome animals: Ecdysozoa (e.g., arthropods, nematodes, onychophorans, and tardigrades) and Lophotrochozoa (e.g., annelids, molluscs, platyhelminths, brachiopods, and rotifers). However, it received little additional support. In fact, several multigene analyses strongly argued against this new phylogeny. These latter studies were based on a large amount of sequence data and therefore showed an apparently strong statistical support. Yet, they covered only a few taxa (those for which complete genomes were available), making systematic artifacts of tree reconstruction more probable. Here we expand this sparse taxonomic sampling and analyze a large data set (146 genes, 35,371 positions) from a diverse sample of animals (35 species). Our study demonstrates that the incongruences observed between rRNA and multigene analyses were indeed due to long-branch attraction artifacts, illustrating the enormous impact of systematic biases on phylogenomic studies. A refined analysis of our data set excluding the most biased genes provides strong support in favor of the new animal phylogeny and in addition suggests that urochordates are more closely related to vertebrates than are cephalochordates. These findings have important implications for the interpretation of morphological and genomic data.     

Phylogenetic analyses of mode of larval development

Phylogenies based on morphological or molecular characters have been used to provide an evolutionary context for analysis of larval evolution. Studies of gastropods, bivalves, tunicates, sea stars, sea urchins, and polychaetes have revealed massive parallel evolution of similar larval forms. Some of these studies were designed to test, and have rejected, the species selection hypothesis for evolutionary trends in the frequency of derived larvae or life history traits. However, the lack of well supported models of larval character evolution leave some doubt about the quality of inferences of larval evolution from phylogenies of living taxa. Better models based on maximum likelihood methods and known prior probabilities of larval character state changes will improve our understanding of the history of larval evolution.     

An empirical codon model for protein sequence evolution

In the past, 2 kinds of Markov models have been considered to describe protein sequence evolution. Codon-level models have been mechanistic with a small number of parameters designed to take into account features, such as transition-transversion bias, codon frequency bias, and synonymous-nonsynonymous amino acid substitution bias. Amino acid models have been empirical, attempting to summarize the replacement patterns observed in large quantities of data and not explicitly considering the distinct factors that shape protein evolution. We have estimated the first empirical codon model (ECM). Previous codon models assume that protein evolution proceeds only by successive single nucleotide substitutions, but our results indicate that model accuracy is significantly improved by incorporating instantaneous doublet and triplet changes. We also find that the affiliations between codons, the amino acid each encodes and the physicochemical properties of the amino acids are main factors driving the process of codon evolution. Neither multiple nucleotide changes nor the strong influence of the genetic code nor amino acids' physicochemical properties form a part of standard mechanistic models and their views of how codon evolution proceeds. We have implemented the ECM for likelihood-based phylogenetic analysis, and an assessment of its ability to describe protein evolution shows that it consistently outperforms comparable mechanistic codon models. We point out the biological interpretation of our ECM and possible consequences for studies of selection.     

Phylogenetic inference of calyptrates,with the first mitogenomes for Gasterophilinae (Diptera: Oestridae) and Paramacronychiinae (Diptera: Sarcophagidae)

The complete mitogenome of the horse stomach bot fly Gasterophilus pecorum (Fabricius) and a near-complete mitogenome of Wohlfahrt''s wound myiasis fly Wohlfahrtia magnifica (Schiner) were sequenced. The mitogenomes contain the typical 37 mitogenes found in metazoans, organized in the same order and orientation as in other cyclorrhaphan Diptera. Phylogenetic analyses of mitogenomes from 38 calyptrate taxa with and without two non-calyptrate outgroups were performed using Bayesian Inference and Maximum Likelihood. Three sub-analyses were performed on the concatenated data: (1) not partitioned; (2) partitioned by gene; (3) 3rd codon positions of protein-coding genes omitted. We estimated the contribution of each of the mitochondrial genes for phylogenetic analysis, as well as the effect of some popular methodologies on calyptrate phylogeny reconstruction. In the favoured trees, the Oestroidea are nested within the muscoid grade. Relationships at the family level within Oestroidea are (remaining Calliphoridae (Sarcophagidae (Oestridae, Pollenia + Tachinidae))). Our mito-phylogenetic reconstruction of the Calyptratae presents the most extensive taxon coverage so far, and the risk of long-branch attraction is reduced by an appropriate selection of outgroups. We find that in the Calyptratae the ND2, ND5, ND1, COIII, and COI genes are more phylogenetically informative compared with other mitochondrial protein-coding genes. Our study provides evidence that data partitioning and the inclusion of conserved tRNA genes have little influence on calyptrate phylogeny reconstruction, and that the 3rd codon positions of protein-coding genes are not saturated and therefore should be included.     

Towards a phylogeny of chitons (Mollusca, Polyplacophora) based on combined analysis of five molecular loci

This study represents the first phylogenetic analysis of the molluscan class Polyplacophora using DNA sequence data. We employed DNA from a nuclear protein-coding gene (histone H3), two nuclear ribosomal genes (18S rRNA and the D3 expansion fragment of 28S rRNA), one mitochondrial protein-coding gene (cytochrome c oxidase subunit I), and one mitochondrial ribosomal gene (16S rRNA). A series of analyses were performed on independent and combined data sets. All these analyses were executed using direct optimization with parsimony as the optimality criterion, and analyses were repeated for nine combinations of parameters affecting indel and transversion/transition cost ratios. Maximum likelihood was also explored for the combined molecular data set, also using the direct optimization method, with a model equivalent to GTR + I + Γ that accommodates gaps. The results of all nine parameter sets for the combined parsimony analysis of all molecular data (as well as ribosomal data) and the maximum-likelihood analysis of all molecular data support monophyly of Polyplacophora. The resulting topologies mostly agree with a division of Polyplacophora into two major lineages: Lepidopleuridae and Chitonida (sensu Sirenko 1993). In our analyses the genus Callochiton is positioned as the sister group to Lepidopleuridae, and not as sister group to the remaining Chitonida (sensu Buckland-Nicks & Hodgson 2000), nor as the sister group to the remaining Chitonina (sensu Buckland-Nicks 1995). Chitonida (excluding Callochiton) is monophyletic, but conventional subgroupings of Chitonida are not supported. Acanthochitonina (sensu Sirenko 1993) is paraphyletic, or alternatively monophyletic, and is split into two clades, both with abanal gills only and cupules in the egg hull, but one has simple cupules whereas the other has more strongly hexagonal cupules. Sister to the Acanthochitonina clades is Chitonina, including taxa with adanal gills and a spiny egg hull. Schizochiton, the only genus with adanal gills that has an egg hull with cupules, is the sister-taxon to one of the Acanthochitonina clades plus Chitonina, or alternatively basal to Chitonina. Support values for either position are low, leaving this relationship unsettled. Our results refute several aspects of conventional classifications of chitons that are based primarily on shell characters, reinforcing the idea that chiton classification should be revised using additional characters.     

Root of the Eukaryota tree as inferred from combined maximum likelihood analyses of multiple molecular sequence data

Extensive studies aiming to establish the structure and root of the Eukaryota tree by phylogenetic analyses of molecular sequences have thus far not resulted in a generally accepted tree. To re-examine the eukaryotic phylogeny using alternative genes, and to obtain a more robust inference for the root of the tree as well as the relationship among major eukaryotic groups, we sequenced the genes encoding isoleucyl-tRNA and valyl-tRNA synthetases, cytosolic-type heat shock protein 90, and the largest subunit of RNA polymerase II from several protists. Combined maximum likelihood analyses of 22 protein-coding genes including the above four genes clearly demonstrated that Diplomonadida and Parabasala shared a common ancestor in the rooted tree of Eukaryota, but only when the fast-evolving sites were excluded from the original data sets. The combined analyses, together with recent findings on the distribution of a fused dihydrofolate reductase-thymidylate synthetase gene, narrowed the possible position of the root of the Eukaryota tree on the branch leading to Opisthokonta or to the common ancestor of Diplomonadida/Parabasala. However, the analyses did not agree with the position of the root located on the common ancestor of Opisthokonta and Amoebozoa, which was argued by Stechmann and Cavalier-Smith [Curr. Biol. 13:R665-666, 2003] based on the presence or absence of a three-gene fusion of the pyrimidine biosynthetic pathway: carbamoyl-phosphate synthetase II, dihydroorotase, and aspartate carbamoyltransferase. The presence of the three-gene fusion recently found in the Cyanidioschyzon merolae (Rhodophyta) genome sequence data supported our analyses against the Stechmann and Cavalier-Smith-rooting in 2003.     

Sweet or salty? The origin of freshwater gastrotrichs (Gastrotricha,Chaetonotida) revealed by molecular phylogenetic analysis

Chaetonotidae is the most diverse and widely distributed family of the order Chaetonotida (Gastrotricha) and includes both marine and freshwater species. Although the family is regarded as a sister taxon to the exclusively marine Xenotrichulidae, the type of environment, marine or freshwater, where Chaetonotidae originated is still not known. Here, we reconstructed the phylogeny of the family based on molecular sequence data and mapped both morphological and ecological characters to determine the ancestral environment of the first members of the family. Our results revealed that the freshwater genus Bifidochaetus is the earliest branching lineage in the paraphyletic Chaetonotidae (encompassing Dasydytidae and Neogosseidae). Moreover, we reconstructed Lepidochaetus-Cephalionotus clade as a monophyletic sister group to the remaining chaetonotids, which supports Kisielewski's morphological based hypothesis concerning undifferentiated type of body scales as a most primary character in Chaetonotidae. We also found that reversals to marine habitats occurred independently in different Chaetonotidae lineages, thus marine species in the genera Heterolepidoderma, Halichaetonotus, Aspidiophorus and subgenera Chaetonotus (Schizochaetonotus) or Chaetonotus (Marinochaetus) should be assumed as having secondarily invaded the marine environment. Character mapping revealed a series of synapomorphies that define the clade that includes Chaetonotidae (with Dasydytidae and Neogosseidae), the most important of which may be those linked to reproduction.     

Statistical methods for detecting molecular adaptation

The past few years have seen the development of powerful statistical methods for detecting adaptive molecular evolution. These methods compare synonymous and nonsynonymous substitution rates in protein-coding genes, and regard a nonsynonymous rate elevated above the synonymous rate as evidence for darwinian selection. Numerous cases of molecular adaptation are being identified in various systems from viruses to humans. Although previous analyses averaging rates over sites and time have little power, recent methods designed to detect positive selection at individual sites and lineages have been successful. Here, we summarize recent statistical methods for detecting molecular adaptation, and discuss their limitations and possible improvements.