Phylogenetics. Which mammalian supertree to bark up?

A supertree of Northern European macromoths

Ecological and life-history data on the Northern European macromoth (Lepidoptera: Macroheterocera) fauna is widely available and ideal for use in answering phylogeny-based research questions: for example, in comparative biology. However, phylogenetic information for such studies lags behind. Here, as a synthesis of all currently available phylogenetic information on the group, we produce a supertree of 114 Northern European macromoth genera (in four superfamilies, with Geometroidea considered separately), providing the most complete phylogenetic picture of this fauna available to date. In doing so, we assess those parts of the phylogeny that are well resolved and those that are uncertain. Furthermore, we identify those genera for which phylogenetic information is currently too poor to include in such a supertree, or entirely absent, as targets for future work. As an aid to studies involving these genera, we provide information on their likely positions within the macromoth tree. With phylogenies playing an ever more important role in the field, this supertree should be useful in informing future ecological and evolutionary studies.

A Total-Group Phylogenetic Metatree for Cetacea and the Importance of Fossil Data in Diversification Analyses

Phylogenetic trees provide a powerful framework for testing macroevolutionary hypotheses, but it is becoming increasingly apparent that inferences derived from extant species alone can be highly misleading. Trees incorporating living and extinct taxa are are needed to address fundamental questions about the origins of diversity and disparity but it has proved challenging to generate robust, species-rich phylogenies that include large numbers of fossil taxa. As a result, most studies of diversification dynamics continue to rely on molecular phylogenies. Here, we extend and apply a recently developed meta-analytic approach for synthesizing previously published phylogenetic studies to infer a well-resolved set of species level, time-scaled phylogenetic hypotheses for extinct and extant cetaceans (whales, dolphins and allies). Our trees extend sampling from the ∼ 90 extant species to over 500 living and extinct species, and therefore allow for more robust inference of macroevolutionary dynamics. While the diversification scenarios we recover are broadly concordant with those inferred from molecular phylogenies they differ in critical ways, notably in the relative contributions of extinction and speciation rate shifts in driving rapid radiations. The metatree approach provides the most immediate route for generating higher level phylogenies of extinct taxa, and opens the door to re-evaluation of macroevolutionary hypotheses derived only from extant taxa.

Evolutionary Models for the Diversification of Placental Mammals Across the KPg Boundary

Deciphering the timing of the placental mammal radiation is a longstanding problem in evolutionary biology, but consensus on the tempo and mode of placental diversification remains elusive. Nevertheless, an accurate timetree is essential for understanding the role of important events in Earth history (e.g., Cretaceous Terrestrial Revolution, KPg mass extinction) in promoting the taxonomic and ecomorphological diversification of Placentalia. Archibald and Deutschman described three competing models for the diversification of placental mammals, which are the Explosive, Long Fuse, and Short Fuse Models. More recently, the Soft Explosive Model and Trans-KPg Model have emerged as additional hypotheses for the placental radiation. Here, we review molecular and paleontological evidence for each of these five models including the identification of general problems that can negatively impact divergence time estimates. The Long Fuse Model has received more support from relaxed clock studies than any of the other models, but this model is not supported by morphological cladistic studies that position Cretaceous eutherians outside of crown Placentalia. At the same time, morphological cladistics has a poor track record of reconstructing higher-level relationships among the orders of placental mammals including the results of new pseudoextinction analyses that we performed on the largest available morphological data set for mammals (4,541 characters). We also examine the strengths and weaknesses of different timetree methods (node dating, tip dating, and fossilized birth-death dating) that may now be applied to estimate the timing of the placental radiation. While new methods such as tip dating are promising, they also have problems that must be addressed if these methods are to effectively discriminate among competing hypotheses for placental diversification. Finally, we discuss the complexities of timetree estimation when the signal of speciation times is impacted by incomplete lineage sorting (ILS) and hybridization. Not accounting for ILS results in dates that are older than speciation events. Hybridization, in turn, can result in dates than are younger or older than speciation dates. Disregarding this potential variation in "gene" history across the genome can distort phylogenetic branch lengths and divergence estimates when multiple unlinked genomic loci are combined together in a timetree analysis.

Fin modules: an evolutionary perspective on appendage disparity in basal vertebrates

BACKGROUND

Fishes are extremely speciose and also highly disparate in their fin configurations, more specifically in the number of fins present as well as their structure, shape, and size. How they achieved this remarkable disparity is difficult to explain in the absence of any comprehensive overview of the evolutionary history of fish appendages. Fin modularity could provide an explanation for both the observed disparity in fin configurations and the sequential appearance of new fins. Modularity is considered as an important prerequisite for the evolvability of living systems, enabling individual modules to be optimized without interfering with others. Similarities in developmental patterns between some of the fins already suggest that they form developmental modules during ontogeny. At a macroevolutionary scale, these developmental modules could act as evolutionary units of change and contribute to the disparity in fin configurations. This study addresses fin disparity in a phylogenetic perspective, while focusing on the presence/absence and number of each of the median and paired fins.

RESULTS

Patterns of fin morphological disparity were assessed by mapping fin characters on a new phylogenetic supertree of fish orders. Among agnathans, disparity in fin configurations results from the sequential appearance of novel fins forming various combinations. Both median and paired fins would have appeared first as elongated ribbon-like structures, which were the precursors for more constricted appendages. Among chondrichthyans, disparity in fin configurations relates mostly to median fin losses. Among actinopterygians, fin disparity involves fin losses, the addition of novel fins (e.g., the adipose fin), and coordinated duplications of the dorsal and anal fins. Furthermore, some pairs of fins, notably the dorsal/anal and pectoral/pelvic fins, show non-independence in their character distribution, supporting expectations based on developmental and morphological evidence that these fin pairs form evolutionary modules.

CONCLUSIONS

Our results suggest that the pectoral/pelvic fins and the dorsal/anal fins form two distinct evolutionary modules, and that the latter is nested within a more inclusive median fins module. Because the modularity hypotheses that we are testing are also supported by developmental and variational data, this constitutes a striking example linking developmental, variational, and evolutionary modules.

The placental mammal ancestor and the post-K-Pg radiation of placentals

To discover interordinal relationships of living and fossil placental mammals and the time of origin of placentals relative to the Cretaceous-Paleogene (K-Pg) boundary, we scored 4541 phenomic characters de novo for 86 fossil and living species. Combining these data with molecular sequences, we obtained a phylogenetic tree that, when calibrated with fossils, shows that crown clade Placentalia and placental orders originated after the K-Pg boundary. Many nodes discovered using molecular data are upheld, but phenomic signals overturn molecular signals to show Sundatheria (Dermoptera + Scandentia) as the sister taxon of Primates, a close link between Proboscidea (elephants) and Sirenia (sea cows), and the monophyly of echolocating Chiroptera (bats). Our tree suggests that Placentalia first split into Xenarthra and Epitheria; extinct New World species are the oldest members of Afrotheria.

Phylogenetic analysis of genome rearrangements among five mammalian orders

Evolutionary relationships among placental mammalian orders have been controversial. Whole genome sequencing and new computational methods offer opportunities to resolve the relationships among 10 genomes belonging to the mammalian orders Primates, Rodentia, Carnivora, Perissodactyla and Artiodactyla. By application of the double cut and join distance metric, where gene order is the phylogenetic character, we computed genomic distances among the sampled mammalian genomes. With a marsupial outgroup, the gene order tree supported a topology in which Rodentia fell outside the cluster of Primates, Carnivora, Perissodactyla, and Artiodactyla. Results of breakpoint reuse rate and synteny block length analyses were consistent with the prediction of random breakage model, which provided a diagnostic test to support use of gene order as an appropriate phylogenetic character in this study. We discussed the influence of rate differences among lineages and other factors that may contribute to different resolutions of mammalian ordinal relationships by different methods of phylogenetic reconstruction.

The historical biogeography of Mammalia

Palaeobiogeographic reconstructions are underpinned by phylogenies, divergence times and ancestral area reconstructions, which together yield ancestral area chronograms that provide a basis for proposing and testing hypotheses of dispersal and vicariance. Methods for area coding include multi-state coding with a single character, binary coding with multiple characters and string coding. Ancestral reconstruction methods are divided into parsimony versus Bayesian/likelihood approaches. We compared nine methods for reconstructing ancestral areas for placental mammals. Ambiguous reconstructions were a problem for all methods. Important differences resulted from coding areas based on the geographical ranges of extant species versus the geographical provenance of the oldest fossil for each lineage. Africa and South America were reconstructed as the ancestral areas for Afrotheria and Xenarthra, respectively. Most methods reconstructed Eurasia as the ancestral area for Boreoeutheria, Euarchontoglires and Laurasiatheria. The coincidence of molecular dates for the separation of Afrotheria and Xenarthra at approximately 100 Ma with the plate tectonic sundering of Africa and South America hints at the importance of vicariance in the early history of Placentalia. Dispersal has also been important including the origins of Madagascar's endemic mammal fauna. Further studies will benefit from increased taxon sampling and the application of new ancestral area reconstruction methods.

Extensive intron gain in the ancestor of placental mammals

Background

Genome-wide studies of intron dynamics in mammalian orthologous genes have found convincing evidence for loss of introns but very little for intron turnover. Similarly, large-scale analysis of intron dynamics in a few vertebrate genomes has identified only intron losses and no gains, indicating that intron gain is an extremely rare event in vertebrate evolution. These studies suggest that the intron-rich genomes of vertebrates do not allow intron gain. The aim of this study was to search for evidence of de novo intron gain in domesticated genes from an analysis of their exon/intron structures.

Results

A phylogenomic approach has been used to analyse all domesticated genes in mammals and chordates that originated from the coding parts of transposable elements. Gain of introns in domesticated genes has been reconstructed on well established mammalian, vertebrate and chordate phylogenies, and examined as to where and when the gain events occurred. The locations, sizes and amounts of de novo introns gained in the domesticated genes during the evolution of mammals and chordates has been analyzed. A significant amount of intron gain was found only in domesticated genes of placental mammals, where more than 70 cases were identified. De novo gained introns show clear positional bias, since they are distributed mainly in 5' UTR and coding regions, while 3' UTR introns are very rare. In the coding regions of some domesticated genes up to 8 de novo gained introns have been found. Intron densities in Eutheria-specific domesticated genes and in older domesticated genes that originated early in vertebrates are lower than those for normal mammalian and vertebrate genes. Surprisingly, the majority of intron gains have occurred in the ancestor of placentals.

Conclusions

This study provides the first evidence for numerous intron gains in the ancestor of placental mammals and demonstrates that adequate taxon sampling is crucial for reconstructing intron evolution. The findings of this comprehensive study slightly challenge the current view on the evolutionary stasis in intron dynamics during the last 100 - 200 My. Domesticated genes could constitute an excellent system on which to analyse the mechanisms of intron gain in placental mammals.

Reviewers: this article was reviewed by Dan Graur, Eugene V. Koonin and Jürgen Brosius.

Heterochrony and developmental modularity of cranial osteogenesis in lipotyphlan mammals

BACKGROUND

Here we provide the most comprehensive study to date on the cranial ossification sequence in Lipotyphla, the group which includes shrews, moles and hedgehogs. This unique group, which encapsulates diverse ecological modes, such as terrestrial, subterranean, and aquatic lifestyles, is used to examine the evolutionary lability of cranial osteogenesis and to investigate the modularity of development.

RESULTS

An acceleration of developmental timing of the vomeronasal complex has occurred in the common ancestor of moles. However, ossification of the nasal bone has shifted late in the more terrestrial shrew mole. Among the lipotyphlans, sequence heterochrony shows no significant association with modules derived from developmental origins (that is, neural crest cells vs. mesoderm derived parts) or with those derived from ossification modes (that is, dermal vs. endochondral ossification).

CONCLUSIONS

The drastic acceleration of vomeronasal development in moles is most likely coupled with the increased importance of the rostrum for digging and its use as a specialized tactile surface, both fossorial adaptations. The late development of the nasal in shrew moles, a condition also displayed by hedgehogs and shrews, is suggested to be the result of an ecological reversal to terrestrial lifestyle and reduced functional importance of the rostrum. As an overall pattern in lipotyphlans, our results reject the hypothesis that ossification sequence heterochrony occurs in modular fashion when considering the developmental patterns of the skull. We suggest that shifts in the cranial ossification sequence are not evolutionarily constrained by developmental origins or mode of ossification.

What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny

BACKGROUND

Efforts to solve higher-level evolutionary relationships within the class Insecta by using mitochondrial genomic data are hindered due to fast sequence evolution of several groups, most notably Hymenoptera, Strepsiptera, Phthiraptera, Hemiptera and Thysanoptera. Accelerated rates of substitution on their sequences have been shown to have negative consequences in phylogenetic inference. In this study, we tested several methodological approaches to recover phylogenetic signal from whole mitochondrial genomes. As a model, we used two classical problems in insect phylogenetics: The relationships within Paraneoptera and within Holometabola. Moreover, we assessed the mitochondrial phylogenetic signal limits in the deeper Eumetabola dataset, and we studied the contribution of individual genes.

RESULTS

Long-branch attraction (LBA) artefacts were detected in all the datasets. Methods using Bayesian inference outperformed maximum likelihood approaches, and LBA was avoided in Paraneoptera and Holometabola when using protein sequences and the site-heterogeneous mixture model CAT. The better performance of this method was evidenced by resulting topologies matching generally accepted hypotheses based on nuclear and/or morphological data, and was confirmed by cross-validation and simulation analyses. Using the CAT model, the order Strepsiptera was recovered as sister to Coleoptera for the first time using mitochondrial sequences, in agreement with recent results based on large nuclear and morphological datasets. Also the Hymenoptera-Mecopterida association was obtained, leaving Coleoptera and Strepsiptera as the basal groups of the holometabolan insects, which coincides with one of the two main competing hypotheses. For the Paraneroptera, the currently accepted non-monophyly of Homoptera was documented as a phylogenetic novelty for mitochondrial data. However, results were not satisfactory when exploring the entire Eumetabola, revealing the limits of the phylogenetic signal that can be extracted from Insecta mitogenomes. Based on the combined use of the five best topology-performing genes we obtained comparable results to whole mitogenomes, highlighting the important role of data quality.

CONCLUSION

We show for the first time that mitogenomic data agrees with nuclear and morphological data for several of the most controversial insect evolutionary relationships, adding a new independent source of evidence to study relationships among insect orders. We propose that deeper divergences cannot be inferred with the current available methods due to sequence saturation and compositional bias inconsistencies. Our exploratory analysis indicates that the CAT model is the best dealing with LBA and it could be useful for other groups and datasets with similar phylogenetic difficulties.

Chemical characterization of milk oligosaccharides of the island flying fox (Pteropus hypomelanus) (Chiroptera: Pteropodidae)

Although a considerable amount of information has accumulated about oligosaccharides in the milk and colostrum of representatives of various mammalian orders, nothing is so far known concerning these sugars in the milk of any bat species (order Chiroptera). In this study, we determined that the following oligosaccharides occur in milk of the island flying fox, Pteropus hypomelanus (Chiroptera: Pteropidae): Gal(α1-3)Gal(β1-4)Glc (isoglobotriose), Gal(β1-4)GlcNAc(β1-3)Gal(β1-4)Glc (lacto-N-neotetraose), Gal(β1-4)GlcNAc(β1-3)[Gal(β1-4)GlcNAc(β1-6)]Gal(β1-4)Glc (lacto-N-neohexaose) and Neu5Gc(α2-3)Gal(β1-4)Glc (3'-NGc-SL). However, lactose was found to be the dominant saccharide in this milk, as in most eutherian mammals. The biologic importance of oligosaccharides in Chiropteran milks warrants further study.

The performance of the Congruence Among Distance Matrices (CADM) test in phylogenetic analysis

BACKGROUND

CADM is a statistical test used to estimate the level of Congruence Among Distance Matrices. It has been shown in previous studies to have a correct rate of type I error and good power when applied to dissimilarity matrices and to ultrametric distance matrices. Contrary to most other tests of incongruence used in phylogenetic analysis, the null hypothesis of the CADM test assumes complete incongruence of the phylogenetic trees instead of congruence. In this study, we performed computer simulations to assess the type I error rate and power of the test. It was applied to additive distance matrices representing phylogenies and to genetic distance matrices obtained from nucleotide sequences of different lengths that were simulated on randomly generated trees of varying sizes, and under different evolutionary conditions.

RESULTS

Our results showed that the test has an accurate type I error rate and good power. As expected, power increased with the number of objects (i.e., taxa), the number of partially or completely congruent matrices and the level of congruence among distance matrices.

CONCLUSIONS

Based on our results, we suggest that CADM is an excellent candidate to test for congruence and, when present, to estimate its level in phylogenomic studies where numerous genes are analysed simultaneously.

Comparative performance of supertree algorithms in large data sets using the soapberry family (Sapindaceae) as a case study

For the last 2 decades, supertree reconstruction has been an active field of research and has seen the development of a large number of major algorithms. Because of the growing popularity of the supertree methods, it has become necessary to evaluate the performance of these algorithms to determine which are the best options (especially with regard to the supermatrix approach that is widely used). In this study, seven of the most commonly used supertree methods are investigated by using a large empirical data set (in terms of number of taxa and molecular markers) from the worldwide flowering plant family Sapindaceae. Supertree methods were evaluated using several criteria: similarity of the supertrees with the input trees, similarity between the supertrees and the total evidence tree, level of resolution of the supertree and computational time required by the algorithm. Additional analyses were also conducted on a reduced data set to test if the performance levels were affected by the heuristic searches rather than the algorithms themselves. Based on our results, two main groups of supertree methods were identified: on one hand, the matrix representation with parsimony (MRP), MinFlip, and MinCut methods performed well according to our criteria, whereas the average consensus, split fit, and most similar supertree methods showed a poorer performance or at least did not behave the same way as the total evidence tree. Results for the super distance matrix, that is, the most recent approach tested here, were promising with at least one derived method performing as well as MRP, MinFlip, and MinCut. The output of each method was only slightly improved when applied to the reduced data set, suggesting a correct behavior of the heuristic searches and a relatively low sensitivity of the algorithms to data set sizes and missing data. Results also showed that the MRP analyses could reach a high level of quality even when using a simple heuristic search strategy, with the exception of MRP with Purvis coding scheme and reversible parsimony. The future of supertrees lies in the implementation of a standardized heuristic search for all methods and the increase in computing power to handle large data sets. The latter would prove to be particularly useful for promising approaches such as the maximum quartet fit method that yet requires substantial computing power.

An application of supertree methods to Mammalian mitogenomic sequences

TWO DIFFERENT APPROACHES CAN BE USED IN PHYLOGENOMICS: combined or separate analysis. In the first approach, different datasets are combined in a concatenated supermatrix. In the second, datasets are analyzed separately and the phylogenetic trees are then combined in a supertree. The supertree method is an interesting alternative to avoid missing data, since datasets that are analyzed separately do not need to represent identical taxa. However, the supertree approach and the corresponding consensus methods have been highly criticized for not providing valid phylogenetic hypotheses. In this study, congruence of trees estimated by consensus and supertree approaches were compared to model trees obtained from a combined analysis of complete mitochondrial sequences of 102 species representing 93 mammal families. The consensus methods produced poorly resolved consensus trees and did not perform well, except for the majority rule consensus with compatible groupings. The weighted supertree and matrix representation with parsimony methods performed equally well and were highly congruent with the model trees. The most similar supertree method was the least congruent with the model trees. We conclude that some of the methods tested are worth considering in a phylogenomic context.

Nomenclature and placental mammal phylogeny

An issue arising from recent progress in establishing the placental mammal Tree of Life concerns the nomenclature of high-level clades. Fortunately, there are now several well-supported clades among extant mammals that require unambiguous, stable names. Although the International Code of Zoological Nomenclature does not apply above the Linnean rank of family, and while consensus on the adoption of competing systems of nomenclature does not yet exist, there is a clear, historical basis upon which to arbitrate among competing names for high-level mammalian clades. Here, we recommend application of the principles of priority and stability, as laid down by G.G. Simpson in 1945, to discriminate among proposed names for high-level taxa. We apply these principles to specific cases among placental mammals with broad relevance for taxonomy, and close with particular emphasis on the Afrotherian family Tenrecidae. We conclude that no matter how reconstructions of the Tree of Life change in years to come, systematists should apply new names reluctantly, deferring to those already published and maximizing consistency with existing nomenclature.

Phylogenetic evidence for early hemochorial placentation in eutheria

The eutherian placenta is remarkable for its structural and functional variability. In order to construct and test comparative hypotheses relating ecological, behavioral and physiological traits to placental characteristics it is first necessary to reconstruct the historical course of placental evolution. Previous attempts to do so have yielded inconsistent results, particularly with respect to the early evolution of structural relationships between fetal and maternal circulatory systems. Here, we bring a battery of phylogenetic methods - including parsimony, likelihood and Bayesian approaches - to bear on the question of placental evolution. All of these approaches are consistent in indicating that highly invasive hemochorial placentation, as found in human beings and numerous other taxa, was an early evolutionary innovation present in the most ancient ancestors of the living placental mammals.

[Combining phylogenetic information: concept, methodology, and challenges]

The DNA sequences, morphological and other homologous characters can be used to infer the origins and histories of biological taxa. Combining all the phylogenetic information available can produce more inclusive phylogenies, improve our understanding of living organisms, and enable biologists to prompt and test hypotheses on a larger scale and with stronger statistical power. In this article, the concept of combining phylogenetic information and its comparison with traditional analysis were reviewed. The most popular approaches of supertree and supermatrix were discussed in detail, and novel ways were presented. Although the combining analysis is facing rigid challenges from data and foundation, it is currently the only approach for realization of the Tree(Net) of Life, and its development will definitely expand our knowledge of evolution on the earth and contribute to the progress of evolutionary related disciplines.

Phylogenetically informed analysis of the allometry of Mammalian Basal metabolic rate supports neither geometric nor quarter-power scaling

The form of the relationship between the basal metabolic rate (BMR) and body mass (M) of mammals has been at issue for almost seven decades, with debate focusing on the value of the scaling exponent (b, where BMR is proportional to M(b)) and the relative merits of b= 0.67 (geometric scaling) and b= 0.75 (quarter-power scaling). However, most analyses are not phylogenetically informed (PI) and therefore fail to account for the shared evolutionary history of the species they consider. Here, we reanalyze the most rigorously selected and comprehensive mammalian BMR dataset presently available, and investigate the effects of data selection and phylogenetic method (phylogenetic generalized least squares and independent contrasts) on estimation of the scaling exponent relating mammalian BMR to M. Contrary to the results of a non-PI analysis of these data, which found an exponent of 0.67-0.69, we find that most of the PI scaling exponents are significantly different from both 0.67 and 0.75. Similarly, the scaling exponents differ between lineages, and these exponents are also often different from 0.67 or 0.75. Thus, we conclude that no single value of b adequately characterizes the allometric relationship between body mass and BMR.

Supertrees and the Tree of Life: generating a metaphylogeny for a diverse invertebrate family (Insecta:Diptera:Therevidae) using constraint trees and the parsimony ratchet to overcome low taxon overlap

The dipteran family Therevidae (stiletto flies) is cosmopolitan and has been the focus of many taxonomic and phylogenetic studies over the last 25 years. Despite this work, questions remain concerning the relationships between subfamilies, genera and generic groups and membership of those groups. We use the supertree method to produce an inclusive phylogeny for the family Therevidae from 24 phylogenetic studies using matrix representation with parsimony (MRP) analysis. The supertree method, one of the most common approaches to calculating globally inclusive phylogenies from smaller more exclusive analyses, produced the therevid metaphylogeny despite only 34% of the terminal taxa being found in more than one source tree. We describe a method for handling low taxon overlap in supertree analyses, in combination with the parsimony ratchet and constraint tree techniques. The supertree presented here is an overarching phylogenetic hypothesis of the Therevidae, incorporating extensive sampling of major lineages and summarising past phylogenetic work on the family. The inclusive metaphylogeny for 362 therevid taxa robustly retrieves the subfamilies Agapophytinae, Phycinae, Therevinae and Xestomyzinae, and the tribes Cyclotelini and Therevini. The Phycinae and Xestomyzinae form a clade, sister to the remaining Therevidae. The Australasian and South American Taenogera Kröber genus-group is monophyletic and sister to a clade of Therevinae and the Australian endemic Agapophytinae. The Therevinae consists of the Anabarhynchus Macquart genus-group of Australian, South American, New Caledonian and New Zealand taxa as sister to the non-Australasian ‘higher Therevinae’, which contains the tribes Cyclotelini and Therevini. The Therevini includes the Hoplosathe Lyneborg & Zaitzev, Litolinga Irwin & Lyneborg, Baryphora Loew, Pandivirilia Irwin & Lyneborg and Thereva Latreille generic-groups. MRP supertree methods can be used to produce inclusive metaphylogenies in situations where source trees have poor data overlap and low taxon overlap, and are therefore valuable in species-rich groups such as arthropods. These methods may be necessary for constructing the ‘Tree of Life’, representing phylogenetic relationships among the millions of known species. However, our analyses show that in situations of source tree conflict, MRP supertree analyses present only the majority signal. We also show that conflict between source trees can be hidden in MRP supertrees, thus our results emphasise the need to evaluate the resulting clades with reference to the source trees.

Confirming the phylogeny of mammals by use of large comparative sequence data sets

The ongoing generation of prodigious amounts of genomic sequence data from myriad vertebrates is providing unparalleled opportunities for establishing definitive phylogenetic relationships among species. The size and complexities of such comparative sequence data sets not only allow smaller and more difficult branches to be resolved but also present unique challenges, including large computational requirements and the negative consequences of systematic biases. To explore these issues and to clarify the phylogenetic relationships among mammals, we have analyzed a large data set of over 60 megabase pairs (Mb) of high-quality genomic sequence, which we generated from 41 mammals and 3 other vertebrates. All sequences are orthologous to a 1.9-Mb region of the human genome that encompasses the cystic fibrosis transmembrane conductance regulator gene (CFTR). To understand the characteristics and challenges associated with phylogenetic analyses of such a large data set, we partitioned the sequence data in several ways and utilized maximum likelihood, maximum parsimony, and Neighbor-Joining algorithms, implemented in parallel on Linux clusters. These studies yielded well-supported phylogenetic trees, largely confirming other recent molecular phylogenetic analyses. Our results provide support for rooting the placental mammal tree between Atlantogenata (Xenarthra and Afrotheria) and Boreoeutheria (Euarchontoglires and Laurasiatheria), illustrate the difficulty in resolving some branches even with large amounts of data (e.g., in the case of Laurasiatheria), and demonstrate the valuable role that very large comparative sequence data sets can play in refining our understanding of the evolutionary relationships of vertebrates.

A likelihood look at the supermatrix-supertree controversy

Supermatrix and supertree methods are two strategies advocated for phylogenetic analysis of sequence data from multiple gene loci, especially when some species are missing at some loci. The supermatrix method concatenates sequences from multiple genes into a data supermatrix for phylogenetic analysis, and ignores differences in evolutionary dynamics among the genes. The supertree method analyzes each gene separately and assembles the subtrees estimated from individual genes into a supertree for all species. Most algorithms suggested for supertree construction lack statistical justifications and ignore uncertainties in the subtrees. Instead of supermatrix or supertree, we advocate the use of likelihood function to combine data from multiple genes while accommodating their differences in the evolutionary process. This combines the strengths of the supermatrix and supertree methods while avoiding their drawbacks. We conduct computer simulation to evaluate the performance of the supermatrix, supertree, and maximum likelihood methods applied to two phylogenetic problems: molecular-clock dating of species divergences and reconstruction of species phylogenies. The results confirm the theoretical superiority of the likelihood method. Supertree or separate analyses of data of multiple genes may be useful in revealing the characteristics of the evolutionary process of multiple gene loci, and the information may be used to formulate realistic models for combined analysis of all genes by likelihood.

Mammalian evolution and biomedicine: new views from phylogeny

Recent progress resolving the phylogenetic relationships of the major lineages of mammals has had a broad impact in evolutionary biology, comparative genomics and the biomedical sciences. Novel insights into the timing and historical biogeography of early mammalian diversification have resulted from a new molecular tree for placental mammals coupled with dating approaches that relax the assumption of the molecular clock. We highlight the numerous applications to come from a well-resolved phylogeny and genomic prospecting in multiple lineages of mammals, from identifying regulatory elements in mammalian genomes to assessing the functional consequences of mutations in human disease loci and those driving adaptive evolution.

Phylogenomic data analyses provide evidence that Xenarthra and Afrotheria are sister groups

The phylogenetic positions of the 4 clades, Euarchontoglires, Laurasiatheria, Afrotheria, and Xenarthra, have been major issues in the recent discussion of basal relationships among placental mammals. However, despite considerable efforts these relationships, crucial to the understanding of eutherian evolution and biogeography, have remained essentially unresolved. Euarchontoglires and Laurasiatheria are generally joined into a common clade (Boreoeutheria), whereas the position of Afrotheria and Xenarthra relative to Boreoeutheria has been equivocal in spite of the use of comprehensive amounts of nuclear encoded sequences or the application of genome-level characters such as retroposons. The probable reason for this uncertainty is that the divergences took place long time ago and within a narrow temporal window, leaving only short common branches. With the aim of further examining basal eutherian relationships, we have collected conserved protein-coding sequences from 11 placental mammals, a marsupial and a bird, whose nuclear genomes have been largely sequenced. The length of the alignment of homologous sequences representing each individual species is 2,168,859 nt. This number of sites, representing 2840 protein-coding genes, exceeds by a considerable margin that of any previous study. The phylogenetic analysis joined Xenarthra and Afrotheria on a common branch, Atlantogenata. This topology was found to fit the data significantly better than the alternative trees.

Supertree bootstrapping methods for assessing phylogenetic variation among genes in genome-scale data sets

Nonparamtric bootstrapping methods may be useful for assessing confidence in a supertree inference. We examined the performance of two supertree bootstrapping methods on four published data sets that each include sequence data from more than 100 genes. In "input tree bootstrapping," input gene trees are sampled with replacement and then combined in replicate supertree analyses; in "stratified bootstrapping," trees from each gene's separate (conventional) bootstrap tree set are sampled randomly with replacement and then combined. Generally, support values from both supertree bootstrap methods were similar or slightly lower than corresponding bootstrap values from a total evidence, or supermatrix, analysis. Yet, supertree bootstrap support also exceeded supermatrix bootstrap support for a number of clades. There was little overall difference in support scores between the input tree and stratified bootstrapping methods. Results from supertree bootstrapping methods, when compared to results from corresponding supermatrix bootstrapping, may provide insights into patterns of variation among genes in genome-scale data sets.

Retroposed elements as archives for the evolutionary history of placental mammals

Reconstruction of the placental mammalian (eutherian) evolutionary tree has undergone diverse revisions, and numerous aspects remain hotly debated. Initial hierarchical divisions based on morphology contained many misgroupings due to features that evolved independently by similar selection processes. Molecular analyses corrected many of these misgroupings and the superordinal hierarchy of placental mammals was recently assembled into four clades. However, long or rapid evolutionary periods, as well as directional mutation pressure, can produce molecular homoplasies, similar characteristics lacking common ancestors. Retroposed elements, by contrast, integrate randomly into genomes with negligible probabilities of the same element integrating independently into orthologous positions in different species. Thus, presence/absence analyses of these elements are a superior strategy for molecular systematics. By computationally scanning more than 160,000 chromosomal loci and judiciously selecting from only phylogenetically informative retroposons for experimental high-throughput PCR applications, we recovered 28 clear, independent monophyly markers that conclusively verify the earliest divergences in placental mammalian evolution. Using tests that take into account ancestral polymorphisms, multiple long interspersed elements and long terminal repeat element insertions provide highly significant evidence for the monophyletic clades Boreotheria (synonymous with Boreoeutheria), Supraprimates (synonymous with Euarchontoglires), and Laurasiatheria. More importantly, two retropositions provide new support for a prior scenario of early mammalian evolution that places the basal placental divergence between Xenarthra and Epitheria, the latter comprising all remaining placentals. Due to its virtually homoplasy-free nature, the analysis of retroposon presence/absence patterns avoids the pitfalls of other molecular methodologies and provides a rapid, unequivocal means for revealing the evolutionary history of organisms.

The authors identified and sequenced retroposons in mammalian genomes. The presence and absence of these retroposons provided evolutionary markers from which the authors reconstructed the phylogenetic history of placental mammals.

A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla)

Despite the biological and economic importance of the Cetartiodactyla, the phylogeny of this clade remains controversial. Using the supertree approach of matrix representation with parsimony, we present the first phylogeny to include all 290 extant species of the Cetacea (whales and dolphins) and Artiodactyla (even-toed hoofed mammals). At the family-level, the supertree is fully resolved. For example, the relationships among the Ruminantia appear as (((Cervidae, Moschidae) Bovidae) (Giraffidae, Antilocapridae) Tragulidae). However, due to either lack of phylogenetic study or contradictory information, polytomies occur within the clades Sus, Muntiacus, Cervus, Delphinidae, Ziphiidae and Bovidae. Complete species-level phylogenies are necessary for both illustrating and analysing biological, geographical and ecological patterns in an evolutionary framework. The present species-level tree of the Cetartiodactyla provides the first opportunity to examine comparative hypotheses across entirely aquatic and terrestrial species within a single mammalian order.

A complete estimate of the phylogenetic relationships in Ruminantia: a dated species-level supertree of the extant ruminants

This paper presents the first complete estimate of the phylogenetic relationships among all 197 species of extant and recently extinct ruminants combining morphological, ethological and molecular information. The composite tree is derived by applying matrix representation using parsimony analysis to 164 previous partial estimates, and is remarkably well resolved, containing 159 nodes (> 80 % of the potential nodes in the completely resolved phylogeny). Bremer decay index has been used to indicate the degree of certainty associated with each clade. The ages of over 80% of the clades in the tree have been estimated from information in the literature. The supertree for Ruminantia illustrates which areas of ruminant phylogeny are still only roughly known because of taxa with controversial relationships (e.g. Odocoileini, Antilopinae) or not studied in great detail (e.g. Muntiacus). It supports the monophyly of the ruminant families and Pecora. According to this analysis Antilocapridae and Giraffidae constitute the superfamily Giraffoidea, which is the sister group of a clade clustering Bovoidea and Cervoidea. The position of several taxa whose systematic positions have remained controversial in the past (Saiga, Pelea, Aepycerus, Pantholops, Ammotragus, Pseudois) is unambiguously established. Nevertheless, the position of Neotragus and Oreotragus within the original radiation of the non-bovine bovids remains unresolved in the present analysis. It also shows that six successive rapid cladogenesis events occurred within the infraorder Pecora during the Oligocene to middle Pliocene, which coincided with periods of global climatic change. Finally, the presented supertree will be a useful framework for comparative and evolutionary biologists interested in studies involving the ruminants.

Evolution of somatosensory and motor cortex in primates

Inferences about how the complex somatosensory systems of anthropoid primates evolved are based on comparative studies of such systems in extant mammals. Experimental studies of members of the major clades of extant mammals suggest that somatosensory cortex of early mammals consisted of only a few areas, including a primary area, S1, bordered by strip-like rostral and caudal somatosensory fields, SR and SC. In addition, the second somatosensory area, S2, and the parietal ventral area, PV, were probably present. S1, S2, and PV were activated independently via parallel projections from the ventroposterior nucleus, VP. Little posterior parietal cortex existed, and it was unlikely that a separate primary motor area, M1, existed until placental mammals evolved. Early primates retained this basic organization and also had a larger posterior parietal region that mediated sensorimotor functions via connections with motor and premotor areas. The frontal cortex included M1, dorsal and ventral premotor areas, supplementary motor area, and cingulate motor fields. Ventroposterior superior and ventroposterior inferior nuclei were distinct from the ventroposterior nucleus in the thalamus. In early anthropoid primates, areas S1, SR, and SC had differentiated into the fields now recognized as areas 3b, 3a, and 1. Areas 3b and 1 contained parallel mirror-image representations of cutaneous receptors and a parallel representation in area 2 was probable. Serial processing became dominant, so that neurons in areas 1, S2, and PV became dependent on area 3b for activation. Posterior parietal cortex expanded into more areas that related to frontal cortex. Less is known about changes that might have occurred with the emergence of apes and humans, but their brains were larger and posed scaling problems most likely solved by increasing the number of cortical areas and reducing the proportion of long connections.

Phylogeny and life histories of the 'Insectivora': controversies and consequences

The evolutionary relationships of the eutherian order Insectivora (Lipotyphla sensu stricto) are the subject of considerable debate. The difficulties in establishing insectivore phylogeny stem from their lack of many shared derived characteristics. The grouping is therefore something of a 'wastebasket' taxon. Most of the older estimates of phylogeny, based on morphological evidence, assumed insectivore monophyly. More recently, molecular phylogenies argue strongly against monophyly, although they differ in the extent of polyphyly inferred for the order. I review the history of insectivore phylogenetics and systematics, focussing on the relationships between the six extant families (Erinaceidae--hedgehogs and moonrats, Talpidae - moles and desmans, Soricidae - shrews, Solenodontidae--solenodons, Tenrecidae--tenrecs and otter-shrews and Chrysochloridae--golden moles). I then examine how these various phylogenetic hypotheses influence the results of comparative analyses and our interpretation of insectivore life-history evolution. I assess which particular controversies have the greatest effect on results, and discuss the implications for comparative analyses where the phylogeny is controversial. I also explore and suggest explanations for certain insectivore life-history trends: increased gestation length and litter size in tenrecs, increased encephalization in moles, and the mixed fast and slow life-history strategies in solenodons. Finally, I consider the implications for comparative analyses of the recent strongly supported phylogenetic hypothesis of an endemic African clade of mammals that includes the insectivore families of tenrecs and golden moles.

Areal specialization of pyramidal cell structure in the visual cortex of the tree shrew: a new twist revealed in the evolution of cortical circuitry

Cortical pyramidal cells, while having a characteristic morphology, show marked phenotypic variation in primates. Differences have been reported in their size, branching structure and spine density between cortical areas. In particular, there is a systematic increase in the complexity of the structure of pyramidal cells with anterior progression through occipito-temporal cortical visual areas. These differences reflect area-specific specializations in cortical circuitry, which are believed to be important for visual processing. However, it remains unknown as to whether these regional specializations in pyramidal cell structure are restricted to primates. Here we investigated pyramidal cell structure in the visual cortex of the tree shrew, including the primary (V1), second (V2) and temporal dorsal (TD) areas. As in primates, there was a trend for more complex branching structure with anterior progression through visual areas in the tree shrew. However, contrary to the trend reported in primates, cells in the tree shrew tended to become smaller with anterior progression through V1, V2 and TD. In addition, pyramidal cells in V1 of the tree shrew are more than twice as spinous as those in primates. These data suggest that variables that shape the structure of adult cortical pyramidal cells differ among species.

A supertree approach to shorebird phylogeny

BACKGROUND

Order Charadriiformes (shorebirds) is an ideal model group in which to study a wide range of behavioural, ecological and macroevolutionary processes across species. However, comparative studies depend on phylogeny to control for the effects of shared evolutionary history. Although numerous hypotheses have been presented for subsets of the Charadriiformes none to date include all recognised species. Here we use the matrix representation with parsimony method to produce the first fully inclusive supertree of Charadriiformes. We also provide preliminary estimates of ages for all nodes in the tree.

RESULTS

Three main lineages are revealed: i) the plovers and allies; ii) the gulls and allies; and iii) the sandpipers and allies. The relative position of these clades is unresolved in the strict consensus tree but a 50% majority-rule consensus tree indicates that the sandpiper clade is sister group to the gulls and allies whilst the plover group is placed at the base of the tree. The overall topology is highly consistent with recent molecular hypotheses of shorebird phylogeny.

CONCLUSION

The supertree hypothesis presented herein is (to our knowledge) the only complete phylogenetic hypothesis of all extant shorebirds. Despite concerns over the robustness of supertrees (see Discussion), we believe that it provides a valuable framework for testing numerous evolutionary hypotheses relating to the diversity of behaviour, ecology and life-history of the Charadriiformes.

The evolution of supertrees

Supertrees result from combining many smaller, overlapping phylogenetic trees into a single, more comprehensive tree. As such, supertree construction is probably as old as the field of systematics itself, and remains our only way of visualizing the Tree of Life as a whole. Over the past decade, supertree construction has gained a more formal, objective footing, and has become an area of active theoretical and practical research. Here, I review the history of the supertree approach, focusing mainly on its current implementation. The supertrees of today represent some of the largest, complete phylogenies available for many groups, but are not without their critics. I conclude by arguing that the ever-growing molecular revolution will result in supertree construction taking on a new role and implementation in the future for analyzing large DNA sequence matrices as part of a divide-and-conquer phylogenetic approach.

A new phylogenetic marker, apolipoprotein B, provides compelling evidence for eutherian relationships

Higher-level relationships within, and the root of Placentalia, remain contentious issues. Resolution of the placental tree is important to the choice of mammalian genome projects and model organisms, as well as for understanding the biogeography of the eutherian radiation. We present phylogenetic analyses of 63 species representing all extant eutherian mammal orders for a new molecular phylogenetic marker, a 1.3kb portion of exon 26 of the apolipoprotein B (APOB) gene. In addition, we analyzed a multigene concatenation that included APOB sequences and a previously published data set (Murphy et al., 2001b) of three mitochondrial and 19 nuclear genes, resulting in an alignment of over 17kb for 42 placentals and two marsupials. Due to computational difficulties, previous maximum likelihood analyses of large, multigene concatenations for placental mammals have used quartet puzzling, less complex models of sequence evolution, or phylogenetic constraints to approximate a full maximum likelihood bootstrap. Here, we utilize a Unix load sharing facility to perform maximum likelihood bootstrap analyses for both the APOB and concatenated data sets with a GTR+Gamma+I model of sequence evolution, tree-bisection and reconnection branch-swapping, and no phylogenetic constraints. Maximum likelihood and Bayesian analyses of both data sets provide support for the superordinal clades Boreoeutheria, Euarchontoglires, Laurasiatheria, Xenarthra, Afrotheria, and Ostentoria (pangolins+carnivores), as well as for the monophyly of the orders Eulipotyphla, Primates, and Rodentia, all of which have recently been questioned. Both data sets recovered an association of Hippopotamidae and Cetacea within Cetartiodactyla, as well as hedgehog and shrew within Eulipotyphla. APOB showed strong support for an association of tarsier and Anthropoidea within Primates. Parsimony, maximum likelihood and Bayesian analyses with both data sets placed Afrotheria at the base of the placental radiation. Statistical tests that employed APOB to examine a priori hypotheses for the root of the placental tree rejected rooting on myomorphs and hedgehog, but did not discriminate between rooting at the base of Afrotheria, at the base of Xenarthra, or between Atlantogenata (Xenarthra+Afrotheria) and Boreoeutheria. An orthologous deletion of 363bp in the aligned APOB sequences proved phylogenetically informative for the grouping of the order Carnivora with the order Pholidota into the superordinal clade Ostentoria. A smaller deletion of 237-246bp was diagnostic of the superordinal clade Afrotheria.

Evaluating placental inter-ordinal phylogenies with novel sequences including RAG1, gamma-fibrinogen, ND6, and mt-tRNA, plus MCMC-driven nucleotide, amino acid, and codon models

It is essential to test a priori scientific hypotheses with independent data, not least to partly negate factors such as gene-specific base composition biases misleading our models. Seven new gene segments and sequences plus Bayesian likelihood phylogenetic methods were used to compare and test five recent placental phylogenies. These five phylogenies are similar to each other, yet quite different from Fthose of previously proposed trees, and span Waddell et al. [Syst. Biol. 48 (1999) 1] to Murphy et al. [Science 294 (2001b) 2348]. Trees for RAG1, gamma-fibrinogen, ND6, mt-tRNA, mt-RNA, c-MYC, epsilon -globin, and GHR are significantly congruent with the four main groups of mammals common to the five phylogenies, i.e., Afrotheria, Laurasiatheria, Euarchontoglires, Xenarthra plus Boreoeutheria (Laurasiatheria plus Euarchontoglires). Where these five a priori phylogenies differ, remain areas generally hard to resolve with the new sequences. The root remains ambiguous and does not reject a basal Afrotheria (the Exafroplacentalia hypothesis), Afrotheria plus Xenarthra together with basal (Atlantogenata), or Epitheria (Xenarthra basal) convincingly. Good evidence is found that Eulipotyphla is monophyletic and is located at the base of Laurasiatheria. The shrew mole, Uropsilus, is found to cluster consistently with other moles, while Solenodon may be the sister taxa to all other eulipotyphlans. Support is found for a probable sister pairing of just hedgehogs/gymnures and shrews. Relationships within Afrotheria, except the Paenungulata clade, remain hard to resolve, although there is congruent support for Afroinsectiphillia (aardvark, elephant shrews, golden moles, and tenrecs). A first-time use is made of MCMC enacted general time-reversible (GTR) amino acid and codon-based models for general tree selection. Even with ND6, a GTR amino acid model provided resolution of fine features, such as the sister group relationship of walrus to Otatriidae, and with BRCA a more reasonable rooting. An extensive analysis of GHR sequences reveals strong congruence with prior phylogenies, including strong support for Eulipotyphla, and good resolution within Rodentia. A codon model gives a worse likelihood than a nucleotide model and sometimes switches support, e.g., with RAG1+gamma-fibrinogen from a hyrax-sirenian association to support for Tethytheria. An analysis of the concatenated data is in accordance with well-resolved features of the gene trees. Taken all together, this work suggests that we are on the right path finding strong confirmation of prior phylogenies. However, with the use of robust criteria for assessing trees (i.e., not Bayesian posteriors), it is apparent parts of the tree remain hard to resolve. Since our current models are far from fitting the sequence data, we should continue with our exploratory analyses to arrive at a refined set of hypotheses for future testing using more model independent characters (e.g., rare indels, gene rearrangement, and SINE data).