As with all data used to infer phylogeny, morphological characters should be independent, and character states homologous. For molecular sequence data, modeling approaches provide objective tests and qualification of these assumptions. Morphological data, however, are more intractable; anatomical character complexes are subject to comparatively greater levels of developmental and functional linkage, ecological convergence, and subjective interpretation ([@B52]; [@B74]; [@B18]; [@B42]; [@B21]; [@B48]; [@B78]; [@B66]; [@B37]). Nevertheless, morphology can be hugely informative under many circumstances. Not only can it increase resolution and reveal hidden support ([@B30]; [@B29]), but because it is usually the only class of data available from fossils, it is the best way to address historic patterns, break up long branches, and calibrate molecular clocks ([@B84]; [@B16]; [@B5]). It is, therefore, necessary to know whether some classes of morphology are more reliable than others. A better understanding of the distribution of homoplasy across different character types would enable more sophisticated treatments of morphological data and would be especially beneficial where the decisiveness of data is poor.

Here we focus on the dental and osteological data of mammals. Dental data have been used alone and alongside other morphological data to reconstruct a wide range of evolutionary transitions from the radiation of mammals ([@B6]), to the relationships of human ancestors ([@B80]). The focus and reliance on dental data are necessitated in large measure by the taphonomic filter particular to the group; the enhanced preservation potential of mammal teeth relative to bones means that many fossil species are known only from dental data. This makes testing for differences between readily fossilizable data (teeth) and less-fossilizable data (osteology) a priority because taphonomic biases can systematically undermine evolutionary inferences ([@B69], [@B70], [@B71], [@B73]; [@B69]). In addition to that, particular concerns about character linkage and overall levels of homoplasy in teeth have been raised on several grounds. Studies of tooth development in mice have revealed dramatic and correlated changes in the number and shape of cusps resulting from minor manipulation of developmental pathways ([@B42]). This suggests that dental characters are nonindependent and this is supported by subsequent developmental studies focusing on phenotypes and traits ([@B37]). Furthermore, measures of phenotypic complexity indicate that the morphology of mammal teeth is strongly tied to function such that form is highly homoplastic and contingent on diet ([@B21]). Aside from developmental and functional linkage, comparison of levels of dissimilarity in molecular and dental characters of phyllostomid bats found the latter to be oversaturated and potentially nonindependent ([@B17]), with the inference that this problem could be more widespread. Nondental anatomy is also expected to exhibit some degree of nonindependence and linkage; for example, serial homology is observed in tetrapod limbs (e.g., [@B65]; [@B88]). Nevertheless, it is mammal teeth for which developmental linkage of phylogenetic characters has been empirically demonstrated ([@B37]).

Investigations of morphological partitions of empirical data (dental, cranial, and postcranial) have found similar levels of homoplasy in some cases (e.g., [@B68]; [@B49]). As such, the characters partitions have been interpreted as equally informative. The same is true for isolated studies of primates ([@B85]) and Cetartiodactylia ([@B53]), although an earlier analysis of the latter's data found teeth to be "markedly different from the rest of the morphological data" ([@B50], p. 451) \[see [@B53]) and [@B50]) for further discussion\]. This raises a number of questions: Why is there an apparent discordance between developmental studies and empirical morphological data regarding the relative phylogenetic informativeness of tooth morphology? Are there differences between the phylogenetic signals contained within dental and osteological morphology? Can dental data alone be used reliably to reconstruct the evolutionary history of mammals? To answer these questions we compiled morphological data sets comprising dental and osteological characters for a broad range of mammal clades, both extant and extinct. Rather than use a single data set as a case study, we take a meta-analytical approach to maximize statistical power, maximize taxonomic coverage, and identify broad-scale patterns. While there are alternative partitions that might be applied to partitions of morphological data \[e.g., cranial vs. postcranial ([@B49]), axial vs. appendicular\], our focus on dental versus osteological characters allows us to address specific concerns raised regarding the developmental and functional linkage and oversaturation of dental characters, and the taphonomic biases particular to mammals (i.e., the enhanced preservation potential of enamel over bone).

It is impossible to know the evolutionary history of empirical taxa with certainty, but congruence between trees inferred from different sources of data offers a means for cross-validation. In this context, we used trees derived from molecular data to assess the congruence of different classes of morphological data. Sequence data provide a qualitatively different and vastly larger data source that serves as a suitable and well-validated benchmark. No single case study is compelling, but meta-analyses of combined data from different sources, clades, and authors seek broad-scale patterns and generalizations from a statistically meaningful sample of independent data (e.g., [@B59]). We, therefore, use a meta-analytical approach to address the following three questions: (1) Do morphological partitions of dental and osteological characters convey a homogenous phylogenetic signal? More specifically, we test the null hypothesis that the partitions do not exhibit significant partition heterogeneity according to the incongruence length difference test (ILD: [@B23], [@B24]) and incongruence relationships difference test (IRD: [@B49]).(2) Are dental and osteological partitions of morphological data sets equally consistent with trees derived from independent molecular sequence data? The specific null hypothesis is that dental and osteological partitions do not retain different levels of relative homoplasy as assessed by the ensemble retention index (RI: [@B22]).(3) Do the generalized taphonomic filters that occur during fossilization (i.e., loss of osteological data, retention of dental data) degrade the phylogentic signal any more or less than equivalent random filters? Specifically, are a similar number of nodes recovered by matrices missing osteological characters compared with matrices missing identical amounts of random characters (as assessed by the node recovery test [@B69]).

MATERIALS AND METHODS {#SEC1}
=====================

Morphological data matrices of mammals were compiled from published sources (*Google* and *Google Scholar* searches for "*clade* phylogeny $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$+$\end{document}$/$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$-$\end{document}$ morphology" and references and citations therein, and *MorphoBank* ([@B54]) from August 2013 to April 2015. Characters were categorized as either osteological or dental on the basis of tissue type (i.e., characters relating directly and explicitly to teeth vs. osteological tissues such as the mandible). Soft tissue characters were excluded. To ensure a balance in the distribution of missing data, taxa with greater than 50% missing entries in either partition were removed, and data sets with a difference of greater than 10% missing entries between partitions were edited further by removing characters with high proportions of missing data (after [@B69]). All uninformative characters were removed. Thresholds were set for the minimum data set dimensions (10 taxa, 30 characters) and minimum percentage of characters in the smallest partition (20%; [@B69]). To ensure independence of data, data sets with appreciable taxonomic overlap were eliminated by prioritizing the more recent source. Rejected data sets, and the reasons for their rejection, are listed in the Supplementary Materials available on Dryad at <http://dx.doi.org/10.5061/dryad.k23mq>.

We applied the ILD test ([@B23], [@B24]; Barker and Lutzoni 2002) to our matrix partitions to test for heterogeneity of signal. The ILD test has the null of partition homogeneity and assesses this by comparing the combined length of shortest trees inferred from each partition to combined lengths of shortest trees inferred from randomly allocated partitions of the same size as the original partitions ([Fig. 1a](#F1){ref-type="fig"}). We used scripts in TNT (Tree analysis with New Technology Goloboff et al. 2008), implementing 999 random replicates. The suitability of the ILD test has been questioned, particularly on the grounds of high type I error rate (false positives) ([@B19]; [@B39]; Planet and Sarkarm 2005; [@B60]; [@B63]; [@B49]), but it remains widely used nonetheless. In addition to the ILD test, we implemented the IRD test of [@B49]). This is also a randomization-based test, but rather than using differences in inferred tree length, partitions are compared via the distances between the optimal trees that result from them. There are many tree-to-tree distance metrics, but here we used the symmetric difference distance (RF; [@B64]; [@B57]) for reasons of familiarity and ease of computation. Sets of multiple optimal trees were compared by calculating the mean RF distance between each tree in one set and its nearest neighbor (the nearest neighbor distance) in the other set ([@B16]). As with the ILD test, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P$\end{document}$ values were approximated from 999 random partitions of the data to assess whether the value from the original partitions (either combined lengths of most parsimonious trees or distance between most parsimonious trees) falls significantly outside the range of values observed in random partitions of the same data.

![Analyses applied to morphological data. a) Partition heterogeneity tested by comparing combined tree length (ILD) or nearest neighbor tree-to-tree distances of trees (IRD) resulting from searches of each morphological character partition (osteological and dental) and random character partitions of the same size ([@B23], [@B24]; [@B49]). b) Molecular consistency tested by applying morphological data to a molecular tree and calculating retention indices of characters and partitions. c) Node recovery tested by comparing strict consensus trees resulting from searches using only dental characters and strict consensus trees resulting from searches using random subsets of characters in the same proportion (i.e., loss of signal with either systematic or random character removal). Signal recovery is assessed in terms of nodes shared with either the strict consensus tree using all characters or the molecular tree. Morphological data set example from [@B25]) with [@B26]) molecular tree.](syw116f1){#F1}

To assess the role of taphonomic biases, we used the node recovery test ([@B69]), also implemented in TNT. The strict consensus trees resulting from maximum parsimony searches using dental characters alone were compared with the strict consensus trees resulting from searches using the entire data set in terms of number shared nodes (i.e., percentage of original signal recovered). This percentage of nodes was then compared with the distribution of percentage of nodes recovered from similar analyses of 500 equivalent matrices, each comprising the same number of characters (i.e., the number of dental characters), but drawn randomly from the entire data set ([Fig. 1c](#F1){ref-type="fig"}).

Molecular trees were sourced in the same way as morphological data. In selecting molecular trees for each morphological data set, priority was given to the degree of taxonomic overlap, and subsequently, the underlying amount of sequence data, rather than the date of publication. Prioritizing taxonomic overlap maximized the data available for analyses because taxa not included in the molecular trees (principally extinct taxa) were removed, and consequently uninformative characters were deleted from all calculations. Both maximum parsimony and likelihood trees were used where available. In the event that multiple trees derived from the same method were presented in a single molecular study, priority was given to the tree derived from the greater amount of sequence data, and subsequently to that which showed greatest resolution. Details of all the molecular trees used are given in the Supplementary Material available on Dryad. Mesquite ([@B47]) was used to construct trees from published figures and to derive retention indices of morphological characters and partitions. Retention indices are influenced by the data set dimensions and, therefore, cannot be directly compared across matrices. Differences between the partitions were, therefore, assessed using $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$t$\end{document}$-tests of the ensemble retention index (RI i.e., for the partition as a whole) and mean character retention index (ri i.e., averaged for all characters in a partition), with values from partitions within the same data set paired.

Molecular data were also used for a modified node recovery test ([Fig. 1c](#F1){ref-type="fig"}). Instead of using the strict consensus tree derived from all morphological characters as the baseline for node recovery, a molecular tree was used. This enabled an assessment of whether trees derived using only dental data recovered significantly fewer molecularly compatible nodes than trees derived from a random subset of morphological characters of the same size drawn from across partitions. For this test, taxa for which molecular data were not available were removed, as were the characters subsequently rendered uninformative following the exclusion of those taxa.

RESULTS {#SEC2}
=======

Forty edited data sets (Supplementary Material available on Dryad) were compiled, constituting a combined total of 1234 taxa and 7403 characters and a near comprehensive sampling of the available morphological data for mammals ([Table 1](#T1){ref-type="table"}).

Table 1.Data sets analyzed and partition heterogeneity test valuesEdited matrixTest $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P$\end{document}$ valuesData sourceIngroupTaxaCharactersILDIRD[@B1])‡Procyonidae12660.2710.467[@B2])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Chrysochloridae30115**0.0040.004**[@B3])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Glires57221**0.004**0.062[@B6])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Mammalia37402**0.001**0.114[@B7])†Notoungulata411190.0720.152[@B8])†Odobenidae1785**0.001**0.090[@B9])†Hippopotamidae14360.1040.246[@B11])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Mustelidae2336**0.0240.042**[@B12])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Muroidea71420.443---[@B13])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Phyllostomidae42620.519**0.019**[@B14])†Rhinocerotidae4171**0.001**---[@B15])‡Otariidae24820.0640.433[@B20])†Sirenia31590.0730.559[@B25])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Arctoidea2575**0.0010.006**[@B27])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Chiroptera231300.114**0.012**[@B28])†Perissodactyls241100.3570.509[@B31])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Feliformia44261**0.003**0.293[@B32])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Bovidae271100.2220.555[@B33])†Proboscidae121130.4240.990[@B34])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Pteropodidae491650.5480.782[@B38])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Erinaceidae1273**0.0010.002**[@B43])†Allotheria1151**0.0020.018**[@B45])†Condylarthra1257**0.013**0.236[@B46])‡Muridae18370.204**0.014**[@B51])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Primatomorpha51893**0.0010.036**[@B53])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Theria462200**0.0010.004**[@B55])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Echimyidae21590.6240.255[@B56])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Didephinae51800.2610.225[@B62])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Macropodidae25810.1830.645[@B67])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Talpidae17152**0.028**0.253[@B75])†Apatemyidae21220**0.001**0.359[@B76])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Artiodactyla50502**0.001**0.371[@B77])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Marsupials1683**0.001**0.172[@B79])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Phyllotini49820.058**0.012**[@B80])‡Hominoidia141030.1630.756[@B81])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Pteromyinae2073**0.001**0.132[@B82])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Carnivora4786**0.001**0.287[@B83])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Oryzomyini52640.6560.333[@B87])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Thylacinidae2570**0.019**0.609[@B89])$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$^{\ast}$\end{document}$Canidae3277**0.007**0.082Totals12347403$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$21/40 <0.05$\end{document}$$\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$11/38<0.05$\end{document}$[^2][^3][^4]

Homogeneity of Signal {#SEC2.1}
---------------------

Application of the ILD test ([@B23], [@B24]) found significant ($\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P <0.05$\end{document}$) heterogeneity between the osteological and dental character partitions in 21 out of 40 data sets, and a highly significant difference overall (Fisher's combined probability $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 4 \times 10^{-29}$\end{document}$, [Table 1](#T1){ref-type="table"}). Applying more stringent thresholds for significance still found widespread heterogeneity between osteological and dental partitions (17 out of 40 data sets had $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P < 0.01$\end{document}$, and 12 out of 40 had $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P < 0.001$\end{document}$). Using the tree-to-tree distance based test (IRD; [@B49]) also identified widespread heterogeneity between dental and osteological partitions (11 out of 38 data sets have $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P <0.05$\end{document}$, Fisher's combined probability $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 1\times 10^{-8})$\end{document}$.

Relative Molecular Consistency {#SEC2.2}
------------------------------

Molecular data are unavailable for almost all extinct taxa \[an exception here being *Thylacinus* from [@B87])\]. Following the removal of fossil taxa, and the characters subsequently rendered uninformative, 14 morphological data sets were omitted from molecular comparisons (of which, 10 were composed largely of fossil taxa and 4 failed to meet the minimum dimensions for inclusion following taxon and character removal). This resulted in 26 data sets ([Table 1](#T1){ref-type="table"}) for which molecular data were available (seven of which molecular data were drawn directly from the same study as the morphological data). The subset of morphological matrices with molecular data available comprised a combined total of 698 taxa and 5589 morphological characters, and like the total data set, it also showed significant partition heterogeneity (Fisher's combined probabilities of $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P< 2 \times 10^{-14}$\end{document}$ for ILD tests and $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P <5 \times 10^{-6}$\end{document}$ for IRD tests). More importantly, osteological characters were found to have significantly higher retention indices than dental characters when optimized onto molecular trees ([Fig. 2](#F2){ref-type="fig"}). This difference was significant both for average retention indices of individual characters (paired $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$t$\end{document}$-test, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.003$\end{document}$) and for ensemble retention indices for whole partitions (paired $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$t$\end{document}$-test, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.036$\end{document}$). Furthermore, the difference was significant irrespective of the tree building methods that were used to infer molecular trees (i.e., maximum likelihood vs. maximum parsimony; see Supplementary Material available on Dryad). We note that it was not necessary to control for differences in numbers of characters in the dental and osteological partitions, using, for example, a jackknifing approach. While the Retention Index of most parsimonious trees is sensitive to differences in data set dimensions, our approach here was to optimize characters onto an existing tree (the molecular) in each case. The retention index (ri) for a given character is then the same regardless of whether this is calculated as part of the original partition, as part of a jackknifed subsample of this partition, or as a singly optimized character.

![Tests for molecular consistency of morphological characters found significant differences between osteological and dental partitions (paired $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$t$\end{document}$-test, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.003$\end{document}$ for average retention indices, ri; paired $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$t$\end{document}$-test, $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.036$\end{document}$ for ensemble retention indices, RI).](syw116f2){#F2}

Recovery of Nodes by Dental and Osteological Character Partitions {#SEC2.3}
-----------------------------------------------------------------

The node recovery test ([@B69]) found that strict consensus trees inferred from exclusively dental characters recovered only 36% of the nodes of strict consensus trees inferred using all morphological characters as aggregated across data sets. This is a significantly lower percentage than strict consensus trees inferred using the same number of characters as dental characters, but drawn at random from across both partitions (median of 43% of aggregated nodes across all studies) ($\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0002$\end{document}$, [Fig. 3](#F3){ref-type="fig"}, see Supplementary Material available on Dryad). Strict consensus trees recovered using only osteological characters also recovered significantly fewer of the original strict consensus nodes compared to the searches using the same number of characters drawn at random from across partitions (47% of aggregated nodes for osteological only searches vs. median of 56% from random partitions of the same size; $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.002$\end{document}$, see Supplementary Material available on Dryad). The modified node recovery test using molecular trees as a baseline for signal recovery yielded slightly different results. Strict consensus trees from searches using both osteological and dental characters of taxa for which molecular data were available recovered 29% of the nodes of molecular trees (as aggregated across data sets). Strict consensus trees from searches using only dental characters of these same taxa recovered just 19.8% of the molecular trees nodes; this is a significantly lower percentage of molecular tree nodes than strict consensus trees inferred using the same number of characters drawn at random from across both partitions (median of 22.2% of aggregated molecular nodes; $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.022$\end{document}$; [Fig. 3](#F3){ref-type="fig"}; Supplementary Material available on Dryad). However, strict consensus trees recovered using only osteological characters fell within the distribution of molecular node recovery seen in equivalent searches using the same number of characters drawn at random from across partitions (25.0% of aggregated molecular nodes for osteological searches vs. median of 25.2% from random partitions of the same size; $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P = 0.431$\end{document}$).

![Results of node recovery tests showing distribution of phylogenetic signal recovery using either total morphological signal as a benchmark (left, percentage aggregate nodes of strict consensus trees inferred using both osteological and dental characters recovered) or molecular data as a benchmark (right, percentage aggregate nodes of molecular trees recovered). Characters were drawn from across both partitions using either the proportion of dental characters (above) or the proportion of osteological characters (below). The node recovery of the original partitions (i.e., dental or osteological, black lines) is mapped onto the distributions of equivalent randomly sampled partitions with estimated $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$P$\end{document}$ values (using proportion of random replicates recovering less signal recovery).](syw116f3){#F3}

DISCUSSION {#SEC3}
==========

Within morphological data sets, osteological and dental character partitions exhibit significant heterogeneity as evidenced from both length-based tests (ILD) and tree-to-tree distance based tests (IRD). This indicates that a different phylogenetic signal is contained within each morphological partition; however, these tests alone do not indicate whether partitions differ in their ability to capture evolutionary history. Since it is impossible to know relationships with certainty, we used congruence with a qualitatively different and operationally distinct data source---namely molecules---as a benchmark of phylogenetic fidelity. Of the two morphological partitions, dental characters are a significantly poorer fit with independent molecular trees ([Fig. 2](#F2){ref-type="fig"}). We interpret this to indicate that teeth are relatively less phylogenetically informative, which bears out predictions from developmental studies ([@B42]; [@B37]) but contradicts findings from smaller scale empirical studies ([@B68]; [@B53];[@B85]). That dental characters are less congruent with molecular data is unfortunate from a palaeontological perspective given that teeth constitute so much of the mammalian fossil record. Here we demonstrate that the application of simplified paleontological filters (the removal of osteology in node recovery tests) yields a loss of signal far in excess of the amount predicted from random incompleteness ([Fig. 3](#F3){ref-type="fig"}). A similar pattern was found in an analysis of a single primate data set by [@B58]) whereby "artificial fossil templates" using largely dental characters were found to result in appreciably less signal recovery than "random templates". Indeed, they state that "for most fossils in our data set known for fewer than 180 sampled characters (i.e., $\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym} 
\usepackage{amsfonts} 
\usepackage{amssymb} 
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
}{}$>50%$\end{document}$ missing data), taphonomically induced bias decreases phylogenetic accuracy" ([@B58], p. 179). Their "random templates" are akin to the random missing data analysis used in our study (unlike their "random states" that we interpret as providing a baseline of node recovery for phylogenetically random data with the same general properties of state frequencies). [@B58]) found that sampling across morphological partitions performed better than sampling predominately from a single partition and interpreted this as supporting the principle that characters should be sampled from as many sources as possible ([@B44]; [@B53]; [@B49]). However, the empirical distribution of characters was such that artificial templates dominated by one partition were always dominated by dental characters. As such, insufficiency of a single partition equates to insufficiency of dental characters in this case. In our meta-analyses of 40 data sets (26 with molecular trees), we have explicitly tested both dental only partitions and osteology only partitions. Both dental only and osteology only partitions recover significantly less signal (percentage of original strict consensus nodes when using all morphology as benchmark) than partitions of equivalent size that sample characters randomly from across partitions. This in itself would imply that characters should be sampled widely from across both partitions. When molecular trees are used as a benchmark of signal recovery, however, dental only searches are found to perform significantly worse than random partitions of equivalent size, whereas osteological only searches do not ([Fig. 3](#F3){ref-type="fig"}). Application of molecular data, therefore, indicates that dental data are generally less-reliable indicators of phylogenetic history than osteological data. Not only do dental characters exhibit elevated levels of relative homoplasy, but they are also unlikely to be sufficient for accurate phylogeny reconstruction on their own.

When results from different studies are compared, a coherent picture emerges. Dental morphology is found to convey a phylogenetic signal that is different to that derived from osteology and is comparatively less consistent with molecular data. That dental characters are more homoplastic accords with the observation that they are oversaturated ([@B17]), as well as with developmental evidence for the relative ease and correlation of shape changes in teeth ([@B42]; [@B37]). Furthermore, loss of osteology (the common palaeontological case) causes significant and disproportionate degradation of phylogenetic signal, such that analyses of dental morphology alone are less likely to reconstruct relationships accurately. Notwithstanding all of the above considerations, we do not advocate the wholesale exclusion of dental morphology from phylogenetic and evolutionary studies. Not only would that exclude many taxa known only from teeth, but also neglect the reliable phylogenetic signal that must exist in at least some dental characters, as well as the hidden support that may exist between dental and other characters ([@B30]); the inclusion of dental characters alongside osteological recovers 29% of molecular tree nodes compared to 25% when using osteology only. Instead, it is necessary to identify subsets of morphology that are less subject to convergence and to directly address the issue of oversaturation and nonindependence ([@B17]) \[e.g., by identifying suites and cliques of convergent or correlated characters ([@B52]; [@B40])\]. The phenomenon appears to be widespread for all mammals; it is possible that it could extend to other vertebrates, although the condition of heterodonty in mammals may have compounded the issue, potentially through failure to properly account for serial homology. In all instances, we strongly advocate meta-analyses of cladistics data sets; it is only when data from a wide diversity of clades are comprehensively sampled and compared that large-scale patterns such as these become apparent ([@B69]; [@B49]).

We are grateful to all the authors who created the data that were included in the meta-analysis. Feedback was kindly provided by Manchester Evolution research group, Rob Asher (Cambridge), and three anonymous reviewers. The authors declare no conflict of interest.

SUPPLEMENTARY MATERIAL {#SEC4}
======================

Data available from the Dryad Digital Repository: <http://dx.doi.org/10.5061/dryad.k23mq>.

FUNDING {#SEC5}
=======

This work was supported by Biotechnology and Biological Sciences Research Council \[BBSRC Manchester DTP Summer research experience placement to T.W. and R.S.S., and grants BB/N015827/1 to R.S.S., BB/K015702/1 and BB/K006754/1 to M.A.W.\]; Natural Environment Research Council \[NERC Fellowship NE/I020253/1 to R. S.S.\]; Leverhulme Trust \[grant F/00351/Z to M. A. W.\]; and John Templeton Foundation \[grant 43915 to Mark Wilkinson and M.A.W.\].

[^1]: *Associate Editor: Mark Friedman*

[^2]: *Note*: Bold values, $\documentclass[12pt]{minimal}
    \usepackage{amsmath}
    \usepackage{wasysym} 
    \usepackage{amsfonts} 
    \usepackage{amssymb} 
    \usepackage{amsbsy}
    \usepackage{upgreek}
    \usepackage{mathrsfs}
    \setlength{\oddsidemargin}{-69pt}
    \begin{document}
    }{}$P < 0.05$\end{document}$.

[^3]: $\documentclass[12pt]{minimal}
    \usepackage{amsmath}
    \usepackage{wasysym} 
    \usepackage{amsfonts} 
    \usepackage{amssymb} 
    \usepackage{amsbsy}
    \usepackage{upgreek}
    \usepackage{mathrsfs}
    \setlength{\oddsidemargin}{-69pt}
    \begin{document}
    }{}$^{\ast}$\end{document}$Data sets for which molecular trees were available.

[^4]: †, ‡Data sets for which molecular data were unavailable: either because the ingroup were largely fossil taxa (b) or failed to meet minimum criteria following editing (c), for example, number of taxa, ratio of osteological to dental characters.
