Molecular systematics of the old world monkey tribe papionini: analysis of the total available genetic sequences

The phylogenetic relationships among the genera of the tribe Papionini are inferred using a taxonomic congruence approach in which gene trees derived for eight unlinked genetic sequence datasets are compared. Population genetics theory predicts that species relationships will be revealed with greater probability when the topology of gene trees from many unlinked loci are found to be congruent. The theory underlying this approach is described. Monophyly of the mangabeys is not supported by any of the gene trees; instead, they are polyphyletic with Cercocebus found to be the sister taxon to Mandrillus in five gene trees (with no conflicting trees), and Lophocebus found to be closely related to Papio and/or Theropithecus in all trees. Theropithecus and Papio are not strongly supported as sister taxa (present in one or two trees only);Lophocebus and Papio are supported as sister taxa in the majority of trees. A close relationship between Mandrillus and Papio is not supported in any of the trees.The relationships among Papio, Lophocebus, and Theropithecus cannot be resolved by congruence, probably due to the short time interval estimated between their divergences. The mtDNA COII sequences are used to estimate divergence dates within the papionins. The internode between the divergences of these species is estimated to be between 290 ka and 370 ka. Lastly, the evolution of morphological features such as long faces, suborbital facial fossae, and terrestrial skeletal adaptations is discussed.     

Phylogenomics with incomplete taxon coverage: the limits to inference

Background  

Phylogenomic studies based on multi-locus sequence data sets are usually characterized by partial taxon coverage, in which sequences for some loci are missing for some taxa. The impact of missing data has been widely studied in phylogenetics, but it has proven difficult to distinguish effects due to error in tree reconstruction from effects due to missing data per se. We approach this problem using a explicitly phylogenomic criterion of success, decisiveness, which refers to whether the pattern of taxon coverage allows for uniquely defining a single tree for all taxa.     

Missing data,incomplete taxa,and phylogenetic accuracy

The problem of missing data is often considered to be the most important obstacle in reconstructing the phylogeny of fossil taxa and in combining data from diverse characters and taxa for phylogenetic analysis. Empirical and theoretical studies show that including highly incomplete taxa can lead to multiple equally parsimonious trees, poorly resolved consensus trees, and decreased phylogenetic accuracy. However, the mechanisms that cause incomplete taxa to be problematic have remained unclear. It has been widely assumed that incomplete taxa are problematic because of the proportion or amount of missing data that they bear. In this study, I use simulations to show that the reduced accuracy associated with including incomplete taxa is caused by these taxa bearing too few complete characters rather than too many missing data cells. This seemingly subtle distinction has a number of important implications. First, the so-called missing data problem for incomplete taxa is, paradoxically, not directly related to their amount or proportion of missing data. Thus, the level of completeness alone should not guide the exclusion of taxa (contrary to common practice), and these results may explain why empirical studies have sometimes found little relationship between the completeness of a taxon and its impact on an analysis. These results also (1) suggest a more effective strategy for dealing with incomplete taxa, (2) call into question a justification of the controversial phylogenetic supertree approach, and (3) show the potential for the accurate phylogenetic placement of highly incomplete taxa, both when combining diverse data sets and when analyzing relationships of fossil taxa.     

Identifying optimal incomplete phylogenetic data sets from sequence databases

We introduce a new method for identifying optimal incomplete data sets from large sequence databases based on the graph theoretic concept of alpha-quasi-bicliques. The quasi-biclique method searches large sequence databases to identify useful phylogenetic data sets with a specified amount of missing data while maintaining the necessary amount of overlap among genes and taxa. The utility of the quasi-biclique method is demonstrated on large simulated sequence databases and on a data set of green plant sequences from GenBank. The quasi-biclique method greatly increases the taxon and gene sampling in the data sets while adding only a limited amount of missing data. Furthermore, under the conditions of the simulation, data sets with a limited amount of missing data often produce topologies nearly as accurate as those built from complete data sets. The quasi-biclique method will be an effective tool for exploiting sequence databases for phylogenetic information and also may help identify critical sequences needed to build large phylogenetic data sets.     

A novel supermatrix approach improves resolution of phylogenetic relationships in a comprehensive sample of danthonioid grasses

Phylogeny reconstruction is challenging when branch lengths vary and when different genetic loci show conflicting signals. The number of DNA sequence characters required to obtain robust support for all the nodes in a phylogeny becomes greater with denser taxon sampling. We test the usefulness of an approach mixing densely sampled, variable non-coding sequences (trnL-F; rpl16; atpB-rbcL; ITS) with sparsely sampled, more conservative protein coding and ribosomal sequences (matK; ndhF; rbcL; 26S), for the grass subfamily Danthonioideae. Previous phylogenetic studies of Danthonioideae revealed extensive generic paraphyly, but were often impeded by insufficient character and taxon sampling and apparent inter-gene conflict. Our variably-sampled supermatrix approach allowed us to represent 79% of the species with up to c. 9900 base pairs for taxa representing the major clades. A 'taxon duplication' approach for taxa with conflicting phylogenetic signals allowed us to combine the data whilst representing the differences between chloroplast and nuclear encoded gene trees. This approach efficiently improves resolution and support whilst maximising representation of taxa and their sometimes composite evolutionary histories, resulting in a phylogeny of the Danthonioideae that will be useful both for a wide range of evolutionary studies and to inform forthcoming realignment of generic delimitations in the subfamily.     

Phylogenetic inference: how much evolutionary history is knowable?

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

Semi-strict supertrees

A method to calculate semi‐strict supertrees is proposed. The semi‐strict supertrees are calculated by creating the matrix that represents all the groups in the source trees (as done in already existing techniques), and then finding the trees determined by the ultra‐clique. The ultra‐clique is defined as the set of characters where each possible subset is compatible with each possible subset from the entire matrix. Finding the ultra‐clique is computationally complex (since in most cases many of the characters have missing entries), but a heuristic method yields reliable results. When the trees have no conflict, or when there are only two trees, the method produces the exact result for any ordering of the input trees and any ordering of the groups within them; when there are more than two trees and they have conflict, a single ordering or sequence can create some spurious groups, but doing multiple sequences eliminates the spurious groups. The method uses only state set operations, and is thus easily implemented in computer programs. Unlike any existing type of supertree, semi‐strict supertrees display all the groups, and only those groups, that are implied by at least some combination of the input trees and contradicted by none. The idea that supertrees should take into account the number of occurences of a given group, so as to retain some groups even in the case of conflict, is discussed; it is argued that a conceptual equivalent of the majority rule consensus is not possible when the sets of taxa differ among trees. Also, when pruning taxa from a set of trees, the supertree can display groups that contradict the consensus for the entire trees, suggesting that supertrees for matrices with very dissimilar sets of taxa should be interpreted with caution. If (for any valid reason) the data cannot be combined in a single matrix, it is advisable that the taxon sets in the matrices be as similar as possible.     

Phylogeny of extant and fossil Juglandaceae inferred from the integration of molecular and morphological data sets

It is widely acknowledged that integrating fossils into data sets of extant taxa is imperative for proper placement of fossils, resolution of relationships, and a better understanding of character evolution. The importance of this process has been further magnified because of the crucial role of fossils in dating divergence times. Outstanding issues remain, including appropriate methods to place fossils in phylogenetic trees, the importance of molecules versus morphology in these analyses, as well as the impact of potentially large amounts of missing data for fossil taxa. In this study we used the angiosperm clade Juglandaceae as a model for investigating methods of integrating fossils into a phylogenetic framework of extant taxa. The clade has a rich fossil record relative to low extant diversity, as well as a robust molecular phylogeny and morphological database for extant taxa. After combining fossil organ genera into composite and terminal taxa, our objectives were to (1) compare multiple methods for the integration of the fossils and extant taxa (including total evidence, molecular scaffolds, and molecular matrix representation with parsimony [MRP]); (2) explore the impact of missing data (incomplete taxa and characters) and the evidence for placing fossils on the topology; (3) simulate the phylogenetic effect of missing data by creating "artificial fossils"; and (4) place fossils and compare the impact of single and multiple fossil constraints in estimating the age of clades. Despite large and variable amounts of missing data, each of the methods provided reasonable placement of both fossils and simulated "artificial fossils" in the phylogeny previously inferred only from extant taxa. Our results clearly show that the amount of missing data in any given taxon is not by itself an operational guideline for excluding fossils from analysis. Three fossil taxa (Cruciptera simsonii, Paleoplatycarya wingii, and Platycarya americana) were placed within crown clades containing living taxa for which relationships previously had been suggested based on morphology, whereas Polyptera manningii, a mosaic taxon with equivocal affinities, was placed firmly as sister to two modern crown clades. The position of Paleooreomunnea stoneana was ambiguous with total evidence but conclusive with DNA scaffolds and MRP. There was less disturbance of relationships among extant taxa using a total evidence approach, and the DNA scaffold approach did not provide improved resolution or internal support for clades compared to total evidence, whereas weighted MRP retained comparable levels of support but lost crown clade resolution. Multiple internal minimum age constraints generally provided reasonable age estimates, but the use of single constraints provided by extinct genera tended to underestimate clade ages.     

Optimal taxonomic groups for biodiversity assessment: a meta‐analytic approach

A fundamental decision in biodiversity assessment is the selection of one or more study taxa, a choice that is often made using qualitative criteria such as historical precedent, ease of detection, or available technical or taxonomic expertise. A more robust approach would involve selecting taxa based on the a priori expectation that they will provide the best possible information on unmeasured groups, but data to inform such hypotheses are often lacking. Using a global meta‐analysis, we quantified the proportion of variability that each of 12 taxonomic groups (at the Order level or above) explained in the richness or composition of other taxa. We then applied optimization to matrices of pairwise congruency to identify the best set of complementary surrogate groups. We found that no single taxon was an optimal surrogate for both the richness and composition of unmeasured taxa if we used simple methods to aggregate congruence data between studies. In contrast, statistical methods that accounted for well‐known drivers of cross‐taxon congruence (spatial extent, grain size, and latitude) lead to the prioritization of similar surrogates for both species richness and composition. Advanced statistical methods were also more effective at describing known ecological relationships between taxa than simple methods, and show that congruence is typically highest between taxonomically and functionally dissimilar taxa. Birds and vascular plants were most frequently selected by our algorithm as surrogates for other taxonomic groups, but the extent to which any one taxon was the ‘optimal’ choice of surrogate for other biodiversity was highly context‐dependent. In the absence of other information – such as in data‐poor areas of the globe, and under limited budgets for monitoring or assessment – ecologists can use our results to assess which taxa are most likely to reflect the distribution of the richness or composition of ‘total’ biodiversity.     

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

Background

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

Methodology/Principal Findings

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

Conclusions

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

The complete sequence and gene organization of the mitochondrial genome of the gadilid scaphopod Siphonondentalium lobatum (Mollusca)

Comparisons of mitochondrial gene sequences and gene arrangements can be informative for reconstructing high-level phylogenetic relationships. We determined the complete sequence of the mitochondrial genome of Siphonodentalium lobatum, (Mollusca, Scaphopoda). With only 13,932 bases, it is the shortest molluscan mitochondrial genome reported so far. The genome contains the usual 13 protein-coding genes, two rRNA and 22 tRNA genes. The ATPase subunit 8 gene is exceptionally short. Several transfer RNAs show truncated TpsiC arms or DHU arms. The gene arrangement of S. lobatum is markedly different from all other known molluscan mitochondrial genomes and shows low similarity even to an unpublished gene order of a dentaliid scaphopod. Phylogenetic analyses of all available complete molluscan mitochondrial genomes based on amino acid sequences of 11 protein-coding genes yield trees with low support for the basal branches. None of the traditionally accepted molluscan taxa and phylogenies are recovered in all analyses, except for the euthyneuran Gastropoda. S. lobatum appears as the sister taxon to two of the three bivalve species. We conclude that the deep molluscan phylogeny is probably beyond the resolution of mitochondrial protein sequences. Moreover, assessing the phylogenetic signal in gene order data requires a much larger taxon sample than is currently available, given the exceptional diversity of this character set in the Mollusca.     

?-tubulin Paralogs Provide a Qualitative Test for a Phylogeny of Cyst Nematodes

Evolutionary relationships among cyst nematodes based on predicted ß-tubulin amino acid and DNA sequence data were compared with phylogenies inferred from ribosomal DNA (ITS1, 5.8S gene, ITS2). The ß-tubulin amino acid data were highly conserved and not useful for phylogenetic inference at the taxonomic level of genus and species. Phylogenetic trees based on ß-tubulin DNA sequence data were better resolved, but the relationships at lower taxonomic levels could not be inferred with confidence. Sequences from single species often appeared in more than one monophyletic clade, indicating the presence of ß-tubulin paralogs (confirmed by Southern blot analysis). For a subset of taxa, good congruence between the two data sets was revealed by the presence of the same putative ß-tubulin gene paralogs in monophyletic groups on the rDNA tree, corroborating the taxon relationships inferred from ribosomal DNA data.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Phylogenetic relationships of spatangoid sea urchins (Echinoidea): taxon sampling density and congruence between morphological and molecular estimates

A phylogeny for 21 species of spatangoid sea urchins is constructed using data from three genes and results compared with morphology-based phylogenies derived for the same taxa and for a much larger sample of 88 Recent and fossil taxa. Different data sets and methods of analysis generate different phylogenetic hypotheses, although congruence tests show that all molecular approaches produce trees that are congruent with each other. By contrast, the trees generated from morphological data differ significantly according to taxon sampling density and only those with dense sampling (after a posteriori weighting) are congruent with molecular estimates. With limited taxon sampling, secondary reversals in deep-water taxa are interpreted as plesiomorphies, pulling them to a basal position. The addition of fossil taxa with their unique character combinations reveals hidden homoplasy and generates a phylogeny that is compatible with molecular estimates. As homoplasy levels were found to be broadly similar across different anatomical structures in the echinoid test, no one suite of morphological characters can be considered to provide more reliable phylogenetic information. Some traditional groupings are supported, including the grouping of Loveniidae, Brissidae and Spatangidae within the Micrasterina, but the Asterostomatidae is shown to be polyphyletic with members scattered amongst at least five different clades. As these are mostly deep-sea taxa, this finding implies multiple independent invasions into the deep sea.     

POLYMORPHIC TAXA, MISSING VALUES AND CLADISTIC ANALYSIS

Abstract Missing values have been used in cladistic analyses when data are unavailable, inapplicable or sometimes when character states are variable within terminal taxa. The practice of scoring taxa as having "missing values" for polymorphic characters introduces errors into the calculation of cladogram lengths and consistency indices because some character change is hidden within terminals. Because these hidden character steps are not counted, the set of most parsimonious cladograms may differ from those that would be found if polymorphic taxa had been broken into monomorphic subunits. In some cases, the trees found when polymorphisms are scored as missing values may not include any of the most parsimonious trees found when the data are scored properly. Additionally, in some cases, polymorphic taxa may be found to be polyphyletic when broken into monomorphic subunits; this is undetected when polymorphisms are treated as missing. Because of these problems, terminal units in cladistic analysis should be based on unique, fixed combinations of characters. Polymorphic taxa should be subdivided into subunits that are monomorphic for each character used in the analysis. Disregarding errors in topology, the additional hidden steps in a cladogram in which polymorphisms are scored as missing can be calculated by a simple formula, based on the observation that if it is assumed that polymorphic terminals include all combinations of character states, 2 ^p − 1 additional steps are required for each taxon in which p polymorphic binary characters are scored as missing values. Thus, when several polymorphisms are scored as missing in the same taxon, very large errors can be introduced into the calculation of tree length.     

How wide to cast the net? Cross‐taxon congruence of species richness,community similarity and indicator taxa in ponds

1. Broad‐scale assessments of biodiversity often rely on the use of surrogate taxa, whose reliability has rarely been tested, particularly in freshwater systems. Here we use data from 46 ponds in two regions of the U.K. to explore the performance of macroinvertebrate taxa as surrogates for the rapid assessment of pond biodiversity. For the four dominant taxonomic groups in these ponds (Chironomidae, Coleoptera, Gastropoda and Trichoptera) we explore cross‐taxon species richness relationships in each of the two regions, and also determine the degree of concordance between the different taxa in accurately representing the similarity relationships between pond assemblages. 2. Patterns of cross‐taxon congruence in species richness were highly variable among taxa and study sites, making the use of a single taxon as a predictor of overall macroinvertebrate species richness problematic. In contrast, all four taxa show >70% congruence with the pattern of community similarity between sites resulting from the entire macroinvertebrate dataset, this result being consistent within and between regions. Canonical correspondence analysis demonstrated that all taxa were related in a similar manner to measured environmental parameters, meaning that limited additional ecological information is gained by including a wider range of pond taxa in rapid site assessment. 3. Single taxonomic groups can, therefore, perform consistently as indicators of community similarity between ponds, and no one taxon dramatically outperforms any other in this respect. The relative merits of the four focal taxa as surrogates for pond invertebrate assemblage composition are discussed with reference to ease of survey, ease of identification and ecological range occupied. 4. It is suggested that Coleoptera have a number of advantages as a surrogate taxon, being diverse, easily sampled, readily identified, taxonomically stable, ecologically well understood and occurring across a wide spectrum of pond types. They are therefore recommended for use as a focal group in rapid pond biodiversity assessments, employing an approach such as ours, which examines patterns of assemblage similarity, rather than species richness alone.     

No speed dating please! Patterns of social preference in male and female house mice

Glass sponges (Class Hexactinellida) are important components of deep-sea ecosystems and are of interest from geological and materials science perspectives. The reconstruction of their phylogeny with molecular data has only recently begun and shows a better agreement with morphology-based systematics than is typical for other sponge groups, likely because of a greater number of informative morphological characters. However, inconsistencies remain that have far-reaching implications for hypotheses about the evolution of their major skeletal construction types (body plans). Furthermore, less than half of all described extant genera have been sampled for molecular systematics, and several taxa important for understanding skeletal evolution are still missing. Increased taxon sampling for molecular phylogenetics of this group is therefore urgently needed. However, due to their remote habitat and often poorly preserved museum material, sequencing all 126 currently recognized extant genera will be difficult to achieve. Utilizing morphological data to incorporate unsequenced taxa into an integrative systematics framework therefore holds great promise, but it is unclear which methodological approach best suits this task. Here, we increase the taxon sampling of four previously established molecular markers (18S, 28S, and 16S ribosomal DNA, as well as cytochrome oxidase subunit I) by 12 genera, for the first time including representatives of the order Aulocalycoida and the type genus of Dactylocalycidae, taxa that are key to understanding hexactinellid body plan evolution. Phylogenetic analyses suggest that Aulocalycoida is diphyletic and provide further support for the paraphyly of order Hexactinosida; hence these orders are abolished from the Linnean classification. We further assembled morphological character matrices to integrate so far unsequenced genera into phylogenetic analyses in maximum parsimony (MP), maximum likelihood (ML), Bayesian, and morphology-based binning frameworks. We find that of these four approaches, total-evidence analysis using MP gave the most plausible results concerning congruence with existing phylogenetic and taxonomic hypotheses, whereas the other methods, especially ML and binning, performed more poorly. We use our total-evidence phylogeny of all extant glass sponge genera for ancestral state reconstruction of morphological characters in MP and ML frameworks, gaining new insights into the evolution of major hexactinellid body plans and other characters such as different spicule types. Our study demonstrates how a comprehensive, albeit in some parts provisional, phylogeny of a larger taxon can be achieved with an integrative approach utilizing molecular and morphological data, and how this can be used as a basis for understanding phenotypic evolution. The datasets and associated trees presented here are intended as a resource and starting point for future work on glass sponge evolution.     

Phylogeny of the Polytrichales (Bryophyta) based on simultaneous analysis of molecular and morphological data

Phylogenetic analyses of Polytrichales were conducted using morphology and sequence data from the chloroplast genes rbcL and rps4 plus the trnL-F gene region, part of the mitochondrial nad5 and the nuclear-encoded 18S rDNA. Our analyses included 46 species representing all genera of Polytrichales. Phylogenetic trees were constructed with simultaneous parsimony analyses of all sequences plus morphology and separate combinations of sequence data only. Results lend support for recognition of Polytrichales as a monophyletic entity. Oedipodium griffithianum appears as a sister taxon to Polytrichales or as a sister taxon of all mosses excluding Sphagnales and Andreaeles. Within Polytrichales, Alophosia and Atrichopsis, species without the adaxial lamellae (in Atrichopsis present but poorly developed on male gametophyte) otherwise typical of the group are sister to the remaining species followed by a clade including Bartramiopsis and Lyellia, species with adaxial lamellae covering only the central portion of the leaves. Six taxa with an exclusively Southern Hemisphere distribution form a grade between the basal lineages and a clade including genera that are mostly confined to the Northern Hemisphere.     

Does adding characters with missing data increase or decrease phylogenetic accuracy?

Missing data are a widely recognized nuisance factor in phylogenetic analyses, and the fear of missing data may deter systematists from including characters that are highly incomplete. In this paper, I used simulations to explore the consequences of including sets of characters that contain missing data. More specifically, I tested whether the benefits of increasing the number of characters outweigh the costs of adding missing data cells to a matrix. The results show that the addition of a set of characters with missing data is generally more likely to increase phylogenetic accuracy than decrease it, but the potential benefits of adding these characters quickly disappear as the proportion of missing data increases. Furthermore, despite the overall trend, adding characters with missing data does decrease accuracy in some cases. In these situations, the missing data entries are not themselves misleading, but their presence may mimic the effects of limited taxon sampling, which can positively mislead. Criteria are discussed for predicting whether adding characters with missing data may increase or decrease accuracy. The results of this study also suggest that accuracy can be increased to a surprising degree by (1) "filling the holes" in a data matrix as much as possible (even when relatively few taxa are missing data), and (2) adding fewer characters scored for all taxa rather than adding a larger number of characters known for fewer taxa. Missing data can also be eliminated from an analysis through the exclusion of incomplete taxa rather than incomplete characters, but this approach may reduce the usefulness of the analysis and (in some cases) the accuracy of the estimated trees.     

The phylogeny of Cetartiodactyla: the importance of dense taxon sampling, missing data, and the remarkable promise of cytochrome b to provide reliable species-level phylogenies

We perform Bayesian phylogenetic analyses on cytochrome b sequences from 264 of the 290 extant cetartiodactyl mammals (whales plus even-toed ungulates) and two recently extinct species, the 'Mouse Goat' and the 'Irish Elk'. Previous primary analyses have included only a small portion of the species diversity within Cetartiodactyla, while a complete supertree analysis lacks resolution and branch lengths limiting its utility for comparative studies. The benefits of using a single-gene approach include rapid phylogenetic estimates for a large number of species. However, single-gene phylogenies often differ dramatically from studies involving multiple datasets suggesting that they often are unreliable. However, based on recovery of benchmark clades-clades supported in prior studies based on multiple independent datasets-and recovery of undisputed traditional taxonomic groups, Cytb performs extraordinarily well in resolving cetartiodactyl phylogeny when taxon sampling is dense. Missing data, however, (taxa with partial sequences) can compromise phylogenetic accuracy, suggesting a tradeoff between the benefits of adding taxa and introducing question marks. In the full data, a few species with a short sequences appear misplaced, however, sequence length alone seems a poor predictor of this phenomenon as other taxa with equally short sequences were not conspicuously misplaced. Although we recommend awaiting a better supported phylogeny based on more character data to reconsider classification and taxonomy within Cetartiodactyla, the new phylogenetic hypotheses provided here represent the currently best available tool for comparative species-level studies within this group. Cytb has been sequenced for a large percentage of mammals and appears to be a reliable phylogenetic marker as long as taxon sampling is dense. Therefore, an opportunity exists now to reconstruct detailed phylogenies of most of the major mammalian clades to rapidly provide much needed tools for species-level comparative studies.