An independent heterotachy model and its implications for phylogeny and divergence time estimation

Resolving the evolution of the mammalian middle ear using Bayesian inference

BACKGROUND

The minute, finely-tuned ear ossicles of mammals arose through a spectacular evolutionary transformation from their origins as a load-bearing jaw joint. This involved detachment from the postdentary trough of the mandible, and final separation from the dentary through resorption of Meckel's cartilage. Recent parsimony analyses of modern and fossil mammals imply up to seven independent postdentary trough losses or even reversals, which is unexpected given the complexity of these transformations. Here we employ the first model-based, probabilistic analysis of the evolution of the definitive mammalian middle ear, supported by virtual 3D erosion simulations to assess for potential fossil preservation artifacts.

RESULTS

Our results support a simple, biologically plausible scenario without reversals. The middle ear bones detach from the postdentary trough only twice among mammals, once each in the ancestors of therians and monotremes. Disappearance of Meckel's cartilage occurred independently in numerous lineages from the Late Jurassic to the Late Cretaceous. This final separation is recapitulated during early development of extant mammals, while the earlier-occurring disappearance of a postdentary trough is not.

CONCLUSIONS

Our results therefore suggest a developmentally congruent and directional two-step scenario, in which the parallel uncoupling of the auditory and feeding systems in northern and southern hemisphere mammals underpinned further specialization in both lineages. Until ~168 Ma, all known mammals retained attached middle ear bones, yet all groups that diversified from ~163 Ma onwards had lost the postdentary trough, emphasizing the adaptive significance of this transformation.

A test for heterotachy using multiple pairs of sequences

Heterotachy is a general term to describe positions that evolve at different rates in different lineages. Heterotachy also can generally be viewed as multivariate rates-across-sites variation, which can be described as randomly drawing rates (or branch lengths) from a multivariate distribution for each branch at each site (Wu J, Susko E. 2009. General heterotachy and distance method adjustments. Mol Biol Evol. 26:2689-2697). Motivated by this result, we propose three new distance-based tests: a heterogeneity test, a heterotachy test, and a within-gene heterotachy test and demonstrate with simulations that they perform well under a wide range of conditions. We also applied the first two tests to two real data sets and found that although all these data sets showed significant evidence of heterotachy, there were subtrees for which the data were consistent with an equal rates or rates-across-sites model.heterogeneity, heterotachy, within-gene heterotachy, covarion model, distance method, hypothesis test.

Mitochondrial genome evolution in fire ants (Hymenoptera: Formicidae)

Background

Complete mitochondrial genome sequences have become important tools for the study of genome architecture, phylogeny, and molecular evolution. Despite the rapid increase in available mitogenomes, the taxonomic sampling often poorly reflects phylogenetic diversity and is often also biased to represent deeper (family-level) evolutionary relationships.

Results

We present the first fully sequenced ant (Hymenoptera: Formicidae) mitochondrial genomes. We sampled four mitogenomes from three species of fire ants, genus Solenopsis, which represent various evolutionary depths. Overall, ant mitogenomes appear to be typical of hymenopteran mitogenomes, displaying a general A+T-bias. The Solenopsis mitogenomes are slightly more compact than other hymentoperan mitogenomes (~15.5 kb), retaining all protein coding genes, ribosomal, and transfer RNAs. We also present evidence of recombination between the mitogenomes of the two conspecific Solenopsis mitogenomes. Finally, we discuss potential ways to improve the estimation of phylogenies using complete mitochondrial genome sequences.

Conclusions

The ant mitogenome presents an important addition to the continued efforts in studying hymenopteran mitogenome architecture, evolution, and phylogenetics. We provide further evidence that the sampling across many taxonomic levels (including conspecifics and congeners) is useful and important to gain detailed insights into mitogenome evolution. We also discuss ways that may help improve the use of mitogenomes in phylogenetic analyses by accounting for non-stationary and non-homogeneous evolution among branches.

Limited effects of among-lineage rate variation on the phylogenetic performance of molecular markers

Variation in substitution rates among evolutionary lineages (among-lineage rate variation or ALRV) has been reported to negatively affect the estimation of phylogenies. When the substitution processes underlying ALRV are modeled inadequately, non-sister taxa with similar substitution rates are estimated incorrectly as sister species due to long-branch attraction. Recent advances in modeling site-specific rate variation (heterotachy) have reduced the impacts of ALRV on phylogeny estimation in several empirical and simulated datasets. However, the addition of parameters to the substitution model reduces power to estimate each parameter correctly, which can also lead to incorrect phylogeny estimation. A potential solution to this problem is to identify the levels of ALRV that negatively impact phylogeny estimation such that molecular markers with non-deleterious levels of ALRV can be identified. To this end, we used analyses of empirical and simulated gene datasets to evaluate whether levels of ALRV identified in a mitochondrial genomic dataset for salamanders negatively impacted phylogeny estimation. We simulated data with and without ALRV, holding all other evolutionary parameters constant, and compared the phylogenetic performance of both simulated and empirical datasets. Overall, we found limited, positive effects of ALRV on phylogeny estimation in this dataset, the majority of which resulted from an increase in substitution rate on short branches. We conclude that ALRV does not always negatively impact phylogeny estimation. Therefore, ALRV can likely be disregarded as a criterion for marker selection in comparable phylogenetic studies.

General heterotachy and distance method adjustments

Heterotachy is a general term to describe positions in a sequence that evolve at different rates in different lineages. Kolaczkowski and Thornton (2004. Performance of maximum parsimony and likelihood phylogenetics when evolution is heterogeneous. Nature 431:980-984.) recently described an intriguing heterotachy model that leads to topological bias for likelihood-based methods and parsimony methods. In this article, we show that heterotachy can generally be viewed as multivariate rates-across-sites variation, which can be described as randomly drawing rates (or branch lengths) from a multivariate distribution for each branch at each site. Motivated by this idea, we propose a pairwise alpha heterotachy adjustment model, which gives us much improved topological estimation in the settings by Kolaczkowski and Thornton (2004).