Identifying heterogeneity in rates of morphological evolution: discrete character change in the evolution of lungfish (Sarcopterygii; Dipnoi)

Quantifying rates of morphological evolution is important in many macroevolutionary studies, and critical when assessing possible adaptive radiations and episodes of punctuated equilibrium in the fossil record. However, studies of morphological rates of change have lagged behind those on taxonomic diversification, and most authors have focused on continuous characters and quantifying patterns of morphological rates over time. Here, we provide a phylogenetic approach, using discrete characters and three statistical tests to determine points on a cladogram (branches or entire clades) that are characterized by significantly high or low rates of change. These methods include a randomization approach that identifies branches with significantly high rates and likelihood ratio tests that pinpoint either branches or clades that have significantly higher or lower rates than the pooled rate of the remainder of the tree. As a test case for these methods, we analyze a discrete character dataset of lungfish, which have long been regarded as "living fossils" due to an apparent slowdown in rates since the Devonian. We find that morphological rates are highly heterogeneous across the phylogeny and recover a general pattern of decreasing rates along the phylogenetic backbone toward living taxa, from the Devonian until the present. Compared with previous work, we are able to report a more nuanced picture of lungfish evolution using these new methods.     

Phylogenetic relationships and divergence times of Charadriiformes genera: multigene evidence for the Cretaceous origin of at least 14 clades of shorebirds

Comparative study of character evolution in the shorebirds is presently limited because the phylogenetic placement of some enigmatic genera remains unclear. We therefore used Bayesian methods to obtain a well-supported phylogeny of 90 recognized genera using 5 kb of mitochondrial and nuclear sequences. The tree comprised three major clades: Lari (gulls, auks and allies plus buttonquails) as sister to Scolopaci (sandpipers, jacanas and allies), and in turn sister to Charadrii (plovers, oystercatchers and allies), as in previous molecular studies. Plovers and noddies were not recovered as monophyletic assemblages, and the Egyptian plover Pluvianus is apparently not a plover. Molecular dating using multiple fossil constraints suggests that the three suborders originated in the late Cretaceous between 79 and 102 Mya, and at least 14 lineages of modern shorebirds survived the mass extinction at the K/T boundary. Previous difficulties in determining the phylogenetic relationships of enigmatic taxa reflect the fact that they are well-differentiated relicts of old, genus-poor lineages. We refrain from suggesting systematic revisions for shorebirds at this time because gene trees may fail to recover the species tree when long branches are connected to deep, shorter branches, as is the case for some of the enigmatic taxa.     

Biases in Maximum Likelihood and Parsimony: A Simulation Approach to a 10-Taxon Case

Biases present in maximum likelihood and parsimony are investigated through a simulation study in a 10-taxon case in which several long branches coexist with short branches in the modeled topology. The performance of these methods is explored while increasing the length of the long branches with different amounts of data. Also, simulations with different taxonomic sampling schemes are examined through this study. The presence of a strong bias in parsimony is corroborated: the well-known long-branch attraction. Likelihood performance is found to be sensitive to the mere presence extreme of branch length disparity, retrieving topologies compatible with long-branch attraction and long-branch repulsion, irrespective of the correctness of the model used.     

Size and number of the hypoglossal nerve foramina in the avian skull and their potential neuroanatomical significance

The cranial openings of nervus hypoglossus, the 12th cranial nerve, are for the first time studied across a broad range of higher avian taxa. This nerve plays an important role in the innervation of the syrinx and exits the skull through a variable number of foramina. Most previous authors described 2–3 foramina nervi hypoglossi (FNH) for neornithine birds, but the number, size, and arrangement of FNH is actually more variable than what is apparent from the literature. In the case of three foramina, there is usually a pair of caudal foramina and a rostral one, but even in closely related taxa, a great variability of the FNH pattern may exist. Many taxa of Neognathae exhibit a quadruple of symmetrically arranged FNH, in others four foramina are strung together in a line. A few taxa show more than four FNH, although in these cases the additional foramina are very small. Of particular interest is the occurrence of a very large caudal FNH in Trochilidae and many species of Passeriformes. This large foramen is suggestive of a correlation with the highly developed vocal capabilities of these birds and may transmit fibers from the tracheosyringeal portion of nucleus nervi hypoglossi, in which case it would be an osteological correlate of vocal capabilities. However, targeted neuroanatomical studies are required to determine how individual hypoglossal foramina relate to hypoglossal roots and their branches, and which of them receive fibers supplying axial, lingual, and syringeal muscles.     

A review of long-branch attraction

The history of long‐branch attraction, and in particular methods suggested to detect and avoid the artifact to date, is reviewed. Methods suggested to avoid LBA‐artifacts include excluding long‐branch taxa, excluding faster evolving third codon positions, using inference methods less sensitive to LBA such as likelihood, the Aguinaldo et al. approach, sampling more taxa to break up long branches and sampling more characters especially of another kind, and the pros and cons of these are discussed. Methods suggested to detect LBA are numerous and include methodological disconcordance, RASA, separate partition analyses, parametric simulation, random outgroup sequences, long‐branch extraction, split decomposition and spectral analysis. Less than 10 years ago it was doubted if LBA occurred in real datasets. Today, examples are numerous in the literature and it is argued that the development of methods to deal with the problem is warranted. A 16 kbp dataset of placental mammals and a morphological and molecular combined dataset of gall waSPS are used to illustrate the particularly common problem of LBA of problematic ingroup taxa to outgroups. The preferred methods of separate partition analysis, methodological disconcordance, and long branch extraction are used to demonstrate detection methods. It is argued that since outgroup taxa almost always represent long branches and are as such a hazard towards misplacing long branched ingroup taxa, phylogenetic analyses should always be run with and without the outgroups included. This will detect whether only the outgroup roots the ingroup or if it simultaneously alters the ingroup topology, in which case previous studies have shown that the latter is most often the worse. Apart from that LBA to outgroups is the major and most common problem; scanning the literature also detected the ill advised comfort of high support values from thousands of characters, but very few taxa, in the age of genomics. Taxon sampling is crucial for an accurate phylogenetic estimate and trust cannot be put on whole mitochondrial or chloroplast genome studies with only a few taxa, despite their high support values. The placental mammal example demonstrates that parsimony analysis will be prone to LBA by the attraction of the tenrec to the distant marsupial outgroups. In addition, the murid rodents, creating the classic “the guinea‐pig is not a rodent” hypothesis in 1996, are also shown to be attracted to the outgroup by nuclear genes, although including the morphological evidence for rodents and Glires overcomes the artifact. The gall wasp example illustrates that Bayesian analyses with a partition‐specific GTR + Γ + I model give a conflicting resolution of clades, with a posterior probability of 1.0 when comparing ingroup alone versus outgroup rooted topologies, and this is due to long‐branch attraction to the outgroup. © The Willi Hennig Society 2005.     

The cytochrome b gene as a phylogenetic marker: the limits of resolution for analyzing relationships among cichlid fishes

The mitochondrial cytochrome b (cyt-b) gene is widely used in systematic studies to resolve divergences at many taxonomic levels. The present study focuses mainly on the utility of cyt-b as a molecular marker for inferring phylogenetic relationship at various levels within the fish family Cichlidae. A total of 78 taxa were used in the present analysis, representing all the major groups in the family Cichlidae (72 taxa) and other families from the suborders Labroidei and Percoidei. Gene trees obtained from cyt-b are compared to a published total evidence tree derived from previous studies. Minimum evolution trees based on cyt-b data resulted in topologies congruent with all previous analyses. Parsimony analyses downweighting transitions relative to transversions (ts1:tv4) or excluding transitions at third codon positions resulted in more robust bootstrap support for recognized clades than unweighted parsimony. Relative rate tests detected significantly long branches for some taxa (LB taxa) which were composed mainly by dwarf Neotropical cichlids. An improvement of the phylogenetic signal, as shown by the four-cluster likelihood mapping analysis, and higher bootstrap values were obtained by excluding LB taxa. Despite some limitations of cyt-b as a phylogenetic marker, this gene either alone or in combination with other data sets yields a tree that is in agreement with the well-established phylogeny of cichlid fish. Received: 11 October 2000 / Accepted: 26 February 2001     

Escaping from the Felsenstein Zone by Detecting Long Branches in Phylogenetic Data

Long branches in a true phylogeny tend to disrupt hierarchical character covariation (phylogenetic signal) in the distribution of traits among organisms. The distortion of hierarchical structure in character-state matrices can lead to errors in the estimation of phylogenetic relationships and inconsistency of methods of phylogenetic inference. Examination of trees distorted by long-branch attraction will not reveal the identities of problematic taxa, in part because the distortion can mask long branches by reducing inferred branch lengths and through errors in branching order. Here we present a simple method for the detection of taxa whose placement in evolutionary trees is made difficult by the effects of long-branch attraction. The method is an extension of a tree-independent conceptual framework of phylogenetic data exploration (RASA). Taxa that are likely to attract are revealed because long branches leave distinct footprints in the distribution of character states among taxa, and these traces can be directly observed in the error structure of the RASA regression. Problematic taxa are identified using a new diagnostic plot called the taxon variance plot, in which the apparent cladistic and phenetic variances contributed by individual taxa are compared. The procedure for identifying long edges employs algorithms solved in polynomial time and can be applied to morphological, molecular, and mixed characters. The efficacy of the method is demonstrated using simulated evolution and empirical evidence of long branches in a set of recently published sequences. We show that the accuracy of evolutionary trees can be improved by detecting and combating the potentially misleading influences of long-branch taxa.     

On the analysis of evolutionary change along single branches in a phylogeny

Comparative biologists are sometimes interested in estimating the evolutionary rate along single branches in a phylogeny. I evaluate two methods by which the evolutionary rate along single branches can be compared with the evolutionary rate throughout the rest of the tree. The first is McPeek's contrasts method, and the second is a likelihood method presented independently in two recently published studies. Although the latter method was developed primarily for the comparison of rates among clades, the approach is equally suited for the analysis of evolutionary rate along single or isolated branches. I find that Type I error is acceptable in both methods but that power and parameter estimation are relatively poor in McPeek's method as it is typically applied.     

Can incomplete taxa rescue phylogenetic analyses from long-branch attraction?

Taxon sampling may be critically important for phylogenetic accuracy because adding taxa can help to subdivide misleading long branches. Although the idea that added taxa can break up long branches was exemplified by a study of "incomplete" fossil taxa, the issue of taxon completeness (i.e., proportion of missing data) has been largely ignored in most subsequent discussions of taxon sampling and long-branch attraction. In this article, I use simulations to test the ability of incomplete taxa to subdivide long branches and improve phylogenetic accuracy in situations of potential long-branch attraction. The results show that for most methods and conditions examined, adding taxa that are only 50% complete may provide similar benefits to adding the same number of complete taxa (suggesting that the advantages of increased taxon sampling may be obtained with less data than previously considered). For parsimony, taxa that are less complete (5% to 25% complete) may often have limited ability to rescue analyses from long-branch attraction. In contrast, highly incomplete taxa can be surprisingly beneficial when using model-based methods. The results also suggest the importance of model-based methods in phylogenetic analyses that combine molecular and fossil data.     

Long-Branch Abstractions

Recent attention has been focused on the sensitivities of various tree reconstructing algorithms to sequence rate heterogeneity (long-branch attraction). Phylogenetic conclusions from two recent empirical studies have been indicted as artifacts attributable to long-branch attraction. Siddall et al. (1995) concluded that Myxozoa are cnidarians and sister group to Polypodium based on 18S rDNA and morphology. Hanelt et al. (1996) argued that this result is due to long-branch attraction. Whiting et al. (1997) concluded that the Strepsiptera are sister group to Diptera based on parsimony analysis of 18S rDNA, 28S rDNA, and morphology. Huelsenbeck (1997) argued that this result also is attributable to long-branch attraction. We demonstrate that the analyses and arguments dismissing these results as the effects of long-branch attraction are fundamentally flawed. The criteria employed by these authors were applied arbitrarily by them to the groups that they did not want, and yet using those same criteria, there is more reason to exclude other taxa besides Polypodium and there is more reason to disbelieve monophyly of Diptera than monophyly of Strepsiptera with Diptera. Moreover, it is asserted, long-branch attraction cannot explain the presence of nematocysts in Myxozoa and halteres in Strepsiptera. For these reasons, and in light of the demonstration that long branches cannot attract each other in their mutual absence, we conclude that the monophyly of Myxozoa + Polypodium and Strepsiptera + Diptera is not due to long-branch attraction. We suggest that maximum likelihood methods are extremely sensitive to taxon and character sampling and that these data sets are demonstrative of the long-branch repulsion problem.     

Phylogenetic inference: how much evolutionary history is knowable?

In order to reconstruct phylogenetic trees from extremely dissimilar sequences it is necessary to estimate accurately the extent of sequence divergence. In this paper a new method of sequence analysis, Markov triple analysis, is developed for determining the relative frequencies of nucleotide substitutions within the three branches of a three-taxon dendrogram. Assuming that nucleotide sites are independently and identically distributed and assuming a Markov model for nucleotide (or protein) evolution, it is shown that the unique Markov matrices can be reconstructed given only the joint probability distribution relating three taxa. (In the much simpler case involving only two taxa and two character states, Markov matrices can also be reconstructed, provided symmetry assumptions are placed on the elements of the matrices.) The method is illustrated using sequence data from the combined first and second codon positions derived from complete human, mouse, and cow mitochondrial sequences.      

Discordance of species trees with their most likely gene trees: the case of five taxa

Under a coalescent model for within-species evolution, gene trees may differ from species trees to such an extent that the gene tree topology most likely to evolve along the branches of a species tree can disagree with the species tree topology. Gene tree topologies that are more likely to be produced than the topology that matches that of the species tree are termed anomalous, and the region of branch-length space that gives rise to anomalous gene trees (AGTs) is the anomaly zone. We examine the occurrence of anomalous gene trees for the case of five taxa, the smallest number of taxa for which every species tree topology has a nonempty anomaly zone. Considering all sets of branch lengths that give rise to anomalous gene trees, the largest value possible for the smallest branch length in the species tree is greater in the five-taxon case (0.1934 coalescent time units) than in the previously studied case of four taxa (0.1568). The five-taxon case demonstrates the existence of three phenomena that do not occur in the four-taxon case. First, anomalous gene trees can have the same unlabeled topology as the species tree. Second, the anomaly zone does not necessarily enclose a ball centered at the origin in branch-length space, in which all branches are short. Third, as a branch length increases, it is possible for the number of AGTs to increase rather than decrease or remain constant. These results, which help to describe how the properties of anomalous gene trees increase in complexity as the number of taxa increases, will be useful in formulating strategies for evading the problem of anomalous gene trees during species tree inference from multilocus data.     

Phylogenetic relationships and rates of allozyme evolution among the lineages of sceloporine sand lizards

Parsimony analysis of characters derived from an electrophoretic survey of allozyme variation in the sceloporine sand lizards indicates that Uma is outside of a clade formed by the rest of the sand lizards and that Cophosaurus and Holbrookia share a more recent common ancestor with one another than either does with Callisaurus. Previous electrophoretic studies used phenetic clustering based on genetic distance data to assess relationships among these taxa. The resulting dendrograms were used to argue that Holbrookia is the sister group of all other sand lizards and that Callisaurus and Cophosaurus are sister taxa. When reanalysed using parsimony methods, the data from these previous studies are found to support the conclusions of the present study, namely, that Uma rather than Holbrookia is the sister group of all other sand lizards and that Cophosaurus is the sister taxon of Holbrookia rather than of Callisaurus. Relative rate tests indicate that the incongruencies between branching diagrams derived from phenetic clustering of genetic distances versus those derived from parsimony analysis of electrophoretic characters are attributable to increased rates of protein evolution in the Holbrookia lineage.     

缺失数据和贝叶斯系统发育分析精确度之探索

The effect of missing data on phylogenetic methods is a potentially important issue in our attempts to reconstruct the Tree of Life. If missing data are truly problematic, then it may be unwise to include species in an analysis that lack data for some characters (incomplete taxa) or to include characters that lack data for some species. Given the difficulty of obtaining data from all characters for all taxa (e.g., fossils), missing data might seriously impede efforts to reconstruct a comprehensive phylogeny that includes all species. Fortunately, recent simulations and empirical analyses suggest that missing data cells are not themselves problematic, and that incomplete taxa can be accurately placed as long as the overall number of characters in the analysis is large. However, these studies have so far only been conducted on parsimony, likelihood, and neighbor-joining methods. Although Bayesian phylogenetic methods have become widely used in recent years, the effects of missing data on Bayesian analysis have not been adequately studied. Here, we conduct simulations to test whether Bayesian analyses can accurately place incomplete taxa despite extensive missing data. In agreement with previous studies of other methods, we find that Bayesian analyses can accurately reconstruct the position of highly incomplete taxa (i.e., 95% missing data), as long as the overall number of characters in the analysis is large. These results suggest that highly incomplete taxa can be safely included in many Bayesian phylogenetic analyses.     

Highly incomplete taxa can rescue phylogenetic analyses from the negative impacts of limited taxon sampling

Background

Phylogenies are essential to many areas of biology, but phylogenetic methods may give incorrect estimates under some conditions. A potentially common scenario of this type is when few taxa are sampled and terminal branches for the sampled taxa are relatively long. However, the best solution in such cases (i.e., sampling more taxa versus more characters) has been highly controversial. A widespread assumption in this debate is that added taxa must be complete (no missing data) in order to save analyses from the negative impacts of limited taxon sampling. Here, we evaluate whether incomplete taxa can also rescue analyses under these conditions (empirically testing predictions from an earlier simulation study).

Methodology/Principal Findings

We utilize DNA sequence data from 16 vertebrate species with well-established phylogenetic relationships. In each replicate, we randomly sample 4 species, estimate their phylogeny (using Bayesian, likelihood, and parsimony methods), and then evaluate whether adding in the remaining 12 species (which have 50, 75, or 90% of their data replaced with missing data cells) can improve phylogenetic accuracy relative to analyzing the 4 complete taxa alone. We find that in those cases where sampling few taxa yields an incorrect estimate, adding taxa with 50% or 75% missing data can frequently (>75% of relevant replicates) rescue Bayesian and likelihood analyses, recovering accurate phylogenies for the original 4 taxa. Even taxa with 90% missing data can sometimes be beneficial.

Conclusions

We show that adding taxa that are highly incomplete can improve phylogenetic accuracy in cases where analyses are misled by limited taxon sampling. These surprising empirical results confirm those from simulations, and show that the benefits of adding taxa may be obtained with unexpectedly small amounts of data. These findings have important implications for the debate on sampling taxa versus characters, and for studies attempting to resolve difficult phylogenetic problems.     

Phylogenetic analyses of small subunit ribosomal RNA coding regions reveal a monophyletic lineage of euglyphid testate amoebae (Order Euglyphida)

The Testaceafilosia includes amoebae with filopodia and with a proteinaceous, agglutinated or siliceous test. To explore the deeper phylogeny of this group, we sequenced the small subunit ribosomal RNA coding region of 13 species, including the first sequence of an amoeba with an agglutinated test, Pseudodifflugia sp. Phylogenetic analyses using maximum parsimony and maximum likelihood methods as well as neighbor joining method yielded the following results: the order Euglyphida forms a monophyletic lineage with the sarcomonads as sister group. The next related taxa are the Chlorarachnea and the unidentified filose strain N-Por. In agreement with the previous studies the Phytomyxea branch off at the base of this lineage. The Monadofilosa (Testaceafilosia and Sarcomonadea) appear monophyletic. The Testaceafilosia are polyphyletic, because Pseudodifflugia sp. is positioned as the sister taxon to the sarcomonads. Within the order Euglyphida Paulinella branches off first, together with Cyphoderia followed by Tracheleuglypha. In maximum likelihood and neighbor joining analyses, the genus Euglypha is monophyletic. The branching pattern within the order Euglyphida reflects the evolution of shell morphology from simple to complex built test.     

Multilocus analyses of an Antarctic fish species flock (Teleostei, Notothenioidei, Trematominae): phylogenetic approach and test of the early-radiation event

Clades that have undergone episodes of rapid cladogenesis are challenging from a phylogenetic point of view. They are generally characterised by short or missing internal branches in phylogenetic trees and by conflicting topologies among individual gene trees. This may be the case of the subfamily Trematominae, a group of marine teleosts of coastal Antarctic waters, which is considered to have passed through a period of rapid diversification. Despite much phylogenetic attention, the relationships among Trematominae species remain unclear. In contrast to previous studies that were mostly based on concatenated datasets of mitochondrial and/or single nuclear loci, we applied various single-locus and multilocus phylogenetic approaches to sequences from 11 loci (eight nuclear) and we also used several methods to assess the hypothesis of a radiation event in Trematominae evolution. Diversification rate analyses support the hypothesis of a period of rapid diversification during Trematominae history and only a few nodes in the hypothetical species tree were consistently resolved with various phylogenetic methods. We detected significant discrepancies among trees from individual genes of these species, most probably resulting from incomplete lineage sorting, suggesting that concatenation of loci is not the most appropriate way to investigate Trematominae species interrelationships. These data also provide information about the possible effects of historic climate changes on the diversification rate of this group of fish.     

Dramatically elevated rate of mitochondrial substitution in lice (Insecta: Phthiraptera)

Few estimates of relative substitution rates, and the underlying mutation rates, exist between mitochondrial and nuclear genes in insects. Previous estimates for insects indicate a 2-9 times faster substitution rate in mitochondrial genes relative to nuclear genes. Here we use novel methods for estimating relative rates of substitution, which incorporate multiple substitutions, and apply these methods to a group of insects (lice, Order: Phthiraptera). First, we use a modification of copath analysis (branch length regression) to construct independent comparisons of rates, consisting of each branch in a phylogenetic tree. The branch length comparisons use maximum likelihood models to correct for multiple substitution. In addition, we estimate codon-specific rates under maximum likelihood for the different genes and compare these values. Estimates of the relative synonymous substitution rates between a mitochondrial (COI) and nuclear (EF-1alpha) gene in lice indicate a relative rate of several 100 to 1. This rapid relative mitochondrial rate (>100 times) is at least an order of magnitude faster than previous estimates for any group of organisms. Comparisons using the same methods for another group of insects (aphids) reveals that this extreme relative rate estimate is not simply attributable to the methods we used, because estimates from aphids are substantially lower. Taxon sampling affects the relative rate estimate, with comparisons involving more closely related taxa resulting in a higher estimate. Relative rate estimates also increase with model complexity, indicating that methods accounting for more multiple substitution estimate higher relative rates.     

Attractive and repulsive functions of Slit are mediated by different receptors in the Drosophila trachea

Oxygen delivery in many animals is enabled by the formation of unicellular capillary tubes that penetrate target tissues to facilitate gas exchange. We show that the tortuous outgrowth of tracheal unicellular branches towards their target tissues is controlled by complex local interactions with target cells. Slit, a phylogenetically conserved axonal guidance signal, is expressed in several tracheal targets and is required both for attraction and repulsion of tracheal branches. Robo and Robo2 are expressed in different branches, and are both necessary for the correct orientation of branch outgrowth. At the CNS midline, Slit functions as a repellent for tracheal branches and this function is mediated primarily by Robo. Robo2 is necessary for the tracheal response to the attractive Slit signal and its function is antagonized by Robo. We propose that the attractive and repulsive tracheal responses to Slit are mediated by different combinations of Robo and Robo2 receptors on the cell surface.     

Tree trunks and branches as optimum mechanical supports of the crown: II. The branches

Using assumptions made in a previous paper, a theory of the shapes of primary branches is developed. Two cases are studied: a primary branch which has a portion denuded of secondary branches, and a primary branch with a continuous load of secondary branches. The first case leads to hyperelliptic integrals, the second—to polynomials of second degree.