Phylogenomics with incomplete taxon coverage: the limits to inference
 Michael J Sanderson^{1}Email author,
 Michelle M McMahon^{2} and
 Mike Steel^{3}
DOI: 10.1186/1471214810155
© Sanderson et al; licensee BioMed Central Ltd. 2010
Received: 20 January 2010
Accepted: 25 May 2010
Published: 25 May 2010
Abstract
Background
Phylogenomic studies based on multilocus 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.
Results
We establish theoretical bounds on the impact of missing data on decisiveness. Results are derived for two contexts: a fixed taxon coverage pattern, such as that observed from an already assembled data set, and a randomly generated pattern derived from a process of sampling new data, such as might be observed in an ongoing comparative genomics sequencing project. Lower bounds on how many loci are needed for decisiveness are derived for the former case, and both lower and upper bounds for the latter. When data are not decisive for all trees, we estimate the probability of decisiveness and the chances that a given edge in the tree will be distinguishable. Theoretical results are illustrated using several empirical examples constructed by mining sequence databases, genomic libraries such as ESTs and BACs, and complete genome sequences.
Conclusion
Partial taxon coverage among loci can limit phylogenomic inference by making it impossible to distinguish among multiple alternative trees. However, even though lack of decisiveness is typical of many sparse phylogenomic data sets, it is often still possible to distinguish a large fraction of edges in the tree.
Background
Ready access to vast sequence databases and new data emerging from largescale sequencing efforts is transforming efforts to infer the phylogenetic history of species. Efforts to build ever larger subtrees of the "tree of life" are pushing the field to grapple with highly speciesrich data sets (plants: 2538 taxa [1]; legumes: 2228 taxa [2]; Asterales: 13,533 taxa [3]; eukaryotes: 73,060 taxa [4]). At the same time, some recent studies exploiting complete genome sequences or genomic libraries include 100s  1000s of loci [5–11]. The former tend to be data sets with many taxa and few loci; the latter few taxa and many loci. However, the distinction is blurring as larger data sets of both kinds have been constructed using semiautomated assembly pipelines [3, 10] and high performance phylogenetic inference packages (e.g. [12]).
Although such vast data sets have effectively solved certain phylogenetic problems, such as the broad outline of yeast phylogeny [5], in other taxa they have uncovered a surprisingly high level of incongruence and uncertainty (reviewed in [13]). Case studies in Drosophila [14], rice [1, 11, 15], house mouse [16], and primates [17] document processes generating incongruence among gene trees in different parts of the genome that contribute nearly as much to the signal as the true species history does. Phylogenetic methods are struggling to keep pace with inferential complications arising from these processes, such as the potential for statistical inconsistency when sequences are concatenated [18, 19]. Indeed, a minor paradigm shift appears to be occurring toward retaining the individuality of signals associated with small genomic regions ("gene trees") rather than simply concatenating all sequence data into an overall analysis [20–24]. This approach comprises socalled "supertree" methods [25] in the broadest sense of that term. In phylogenomic contexts this includes methods that explicitly try to build species trees while reconciling gene trees by modeling processes at the interface between phylogenetics and population genetics, such as lineage sorting [24, 26, 27], recombination [28] and hybridization [29, 30].
In addition to biological processes which evidently confound phylogenomic inference are complexities that stem from the origin of the data themselves. One of the hallmarks of typical phylogenomic data sets is partial taxon coverage among loci–by which we mean not merely missing nucleotides within an alignment but loci for which an entire sequence is missing for some taxa (Sanderson et al. [31]; also referred to as "taxon occupancy" by Hejnol et al. [10]). Missing data arise from causes including biases in sequence databases, the technicalities of sequencing strategies, and real evolutionary processes of sequence divergence or gene loss. Whenever sequences from genomic libraries such as ESTs or BACends provide only partial coverage of a genome, the mismatch between regions covered in different species will generate missing data. Whole genome shotgun sequencing using technologies that assemble at low coverage values will likewise entail missing data, especially at broader phylogenetic scales. Even finished whole genomes can fail to have homologs in a collection of species for various reasons ranging from accelerated sequence divergence to gene loss [32]. It is not unusual for finished genomes to have numerous genes with no detectable homologs in any other taxon (e.g., the maize genome has 2000 unique genes not found in close relatives: [33]).
The consequences of missing data for phylogenetic inference were first studied in the context of integrating fossil taxa, in which there are typically large subsets of missing morphological data, with more complete extant taxa [34–37]. More recently the impact of missing data on accuracy of inferred trees has been examined in the more general context of molecular and phylogenomic data sets [11, 38–41] where missing data arises for a variety of reasons. Collectively this body of work treats missing data at the scale of individual cells in a taxon by character data matrix or sequence alignment up to entirely missing sequences for some taxa in a supermatrix or supertree setting, which is the phylogenomic context we consider here. Conclusions about the negative impact of missing data on accuracy have ranged from optimistic [37, 42, 43] to somewhat more cautious [38, 41]; it is possibly fair to say that many empirical studies have argued that the impact of missing data tends to be unbiased and can be ameliorated by sufficient overall quantity of data (e.g., [10]), and most simulation work agrees (see exceptions: [35, 41]). However, there has been frequent speculation about the distribution of missing data in the matrix (e.g. its "evenness", [38]: page 546; [10]: SOM). General conclusions have remained elusive both because of lack of characterization of this distribution and the difficulty of separating the usual problems of phylogenetic inference–error and bias caused by the data–from the simple absence of data.
In the phylogenomic context of multiple loci, a natural question is when can phylogenetic information from different loci be combined to produce a good tree with all taxa present, even when taxa are missing for some loci? Ideally this question should be abstracted away from the details of tree reconstruction methods. We derived some formal results recently for a relatively simple notion of combinability  "decisiveness" ([44]; note, Goloboff [45] used this term with a very different definition). A pattern of taxon coverage is decisive for some tree on all the taxa in a multilocus data set if that tree is uniquely defined by combining the (correctly inferred and perfectly concordant) trees from the separate loci. Our previous paper [44] focused on a strong sense of decisiveness: decisiveness for all trees. In this paper we also consider weaker notions of decisiveness, derive further results on decisiveness in the important context of sampling from genomescale data, and examine the utility of the concept in real phylogenomic data sets. The theory will permit us to make relatively precise statements about matters such as the amount and distribution of missing entries, at least within bounds. The scale of phylogenomic data sets and the sampling properties of genomic sequencing protocols is likely to increase the relevance of our mathematical theory.
Methods
Overview
Questions about taxon coverage arise in at least two contexts. In the first, taxon coverage is given a priori. This might happen when a data matrix has already been assembled from GenBank, or in an examination of a data set from the literature. In this case, the true tree is unknown, and the key question is whether the data set is decisive for all trees, some trees, or no trees. In the second context, taxon coverage can be modeled as the outcome of a sampling process. This might result from a sequencing strategy, for example, that induces a certain quantity of missing data, but for which those data are more or less randomly distributed. Here the pertinent question is how many loci must be sampled to guarantee decisiveness for the unknown true tree.
Definitions
We begin by recalling some basic definitions from phylogenetic theory. Following Semple and Steel [46], given a set X of taxa, a phylogenetic Xtree is a (unrooted) tree T in which the leaves, or terminal taxa, of T consist of the label set X (taxon names) and all the remaining vertices (nodes) of T are unlabelled and have at least three incident edges. T is binary (completely bifurcating) if every internal vertex has exactly three incident edges. For a phylogenetic tree T and a subset Y of X, let TY denote the induced phylogenetic tree on leaf set Y (the tree obtained from the minimal subtree connecting Y by suppressing any vertices of degree two). We will let n = X throughout. A quartet tree is a binary phylogenetic tree on four leaves. For such a tree with leaves a, b, c, d, we write abcd if the internal edge of the tree separates the pair a, b from c, d.
A rooted phylogenetic tree is a tree obtained from a phylogenetic tree by deleting some leaf and its incident edge and regarding the other endpoint of that edge as a root vertex. Rather than defining the notions that follow for both rooted and unrooted phylogenetic trees separately, we will define them just for unrooted trees; the analogous definitions for rooted trees follow immediately by regarding any rooted tree on a set of n taxa as an unrooted phylogenetic tree on n+1 taxa, the additional taxon being some fixed reference outgroup (real or hypothetical).
We now define some key concepts. Let Y be the (label) set of taxa with nonmissing data for some locus (Y_{ i }if locus i is specified).

Whenever Y ⊆ X, we say that a binary phylogenetic Xtree T displays a binary phylogenetic Ytree, T' if TY = T'. More generally, for Y ⊆ X, we say that a phylogenetic Xtree T displays a phylogenetic Ytree, T' if TY either equals T' or is a resolution of that tree (i.e. all the splits of T' are contained in TY).

We say that a collection of phylogenetic trees T_{1},...,T_{ k }defines a phylogenetic Xtree T if X is the union of the leaf sets of the trees T_{1},...,T_{ k }and there exists one, and only one phylogenetic Xtree that displays these trees, and this tree is T. This implies T that must be binary.

An internal edge e of a binary tree T is said to be distinguished by some other tree if both trees display a quartet tree xywz but the tree T/e (the tree obtained from T by collapsing e) does not display this quartet (this means that x, y, w, z comes from each of the four subtrees of T incident with e). There is an analogous notion of 'distinguished' for rooted trees based on rooted triples rather than quartet trees. For example, in Figure 2A the internal edge below the clade C, D of the original tree (left) is not distinguished by either of the two subtrees (middle) shown; however in Figure 2B all three interior edges of the original tree (left) are distinguished by the subtrees (middle).

For a collection S = {Y_{1}, ..., Y_{ k }} of subsets of X, we say that S is the pattern of taxon coverage or just taxon coverage (Figure 1). Any taxon found in all k subsets is called a reference taxon. The taxon coverage density is .

We say that S is decisive for a phylogenetic tree T provided TY_{1},...,TY_{ k }defines (Figure 2).

We say that S is phylogenetically decisive for all trees if it satisfies the following property: If T and T' are binary phylogenetic Xtrees, with TY = T'Y for all Y ∈ S then T = T'. In other words, for any binary phylogenetic Xtree T, the collection of induced subtrees {TY:Y ∈ S} defines T.
These last two definitions are from [44], although we now adopt the convention of saying "decisive for all trees" instead of merely "decisive". Notice that a pattern of taxon coverage that is decisive for a given tree can involve far fewer loci than any pattern of taxon coverage that is decisive for all trees; for instance, for every binary phylogenetic tree T, there is a set S of just n  3 quadruples (set of four taxa) for which S is decisive for T [47] but any set of quadruples that is decisive for all trees must have size that is at least cubic in n (this follows from Theorem 1, below).
Notice that there is a subtle but fundamental distinction between questions about whether an arbitrary collection of input trees, T_{1},...,T_{ k }, defines a unique tree, T, on the one hand, and whether a pattern of taxon coverage, S = {Y_{1},...,Y_{ k }}, together with a tree, T, induces a set of subtrees, TY_{1},...,TY_{ k }, that uniquely defines T. Intuitively, in the first case the focus is on the subtrees, whereas in the second the focus is the pattern of taxon coverage. The term "define" refers to the first context; "decisive" to the second. What does it mean when taxon coverage is not decisive for all trees? Then there exists two or more distinct trees which cannot be differentiated from one another based only on the subtrees induced by the taxon coverage pattern, because these trees induce the same set of subtrees. Of course, these may or may not include the "true tree" for any given problem in phylogenetic inference, but surely it is risky to be unable to distinguish among possibilities in case the true tree is involved.
Context 1: Taxon coverage is fixed a priori
[50], where n_{ j }= Y_{ j }. As noted in this last paper, Inequality (2) is not tight in the following sense: certain nineleaf binary trees cannot be defined by any two sixleaf subtrees, even though Inequality (2) is satisfied since (63) + (63)≥ (93).
Similar bounds can be obtained for phylogenetic decisiveness for all trees.
Theorem 1 If a collection S = {Y_{1}, ..., Y_{ k }} of subsets of X is phylogenetically decisive for all trees then: .
(See Appendix for proof). This provides a necessary condition, or lower bound, for phylogenetic decisiveness for all trees. Roughly speaking, if some fraction of taxa, f = m/n, is present for each locus, then we need at least (1/f)^{3} loci for decisiveness for all trees, while at least (1/f) loci are required for decisiveness for a given tree.
Partial decisiveness
Determining when taxon coverage is decisive for only some trees involves an important open question. When exactly does a set of subtrees define a tree? The bound stated earlier does not specify how to tell if some set of n3 or more quartet trees defines a tree. The answer is known for rooted [47], but not for unrooted trees. A collection of rooted trees defines a tree if two conditions are met: each subtree is displayed by the tree, and each internal edge of the tree is distinguished by an internal edge of at least one of the input trees ([47]: see definition above). In contrast to the rooted setting, distinguishing all the edges of an unrooted tree is only a necessary (but not sufficient) condition for a collection of unrooted trees to define that tree. To proceed, we therefore restrict attention to the rooted case, which is equivalent to assuming there exists in the data set a reference taxon which can provisionally be used to root the input trees (in an artificial but useful sense of "rooting").
For rooted trees, the fraction of trees, ∏_{ D }, that are decisive for a given pattern of taxon coverage could be determined by checking these conditions for all possible trees. For small sets of taxa complete enumeration is feasible. For larger data sets, the true fraction of trees for which the taxon coverage is decisive can be estimated from the fraction in a large random sample of trees (under a uniform/equiprobable distribution of trees).
However, the condition of distinguishability of edges suggests a potentially more sensitive measure of "partial" decisiveness: the probability of an edge being distinguished on a tree for a given pattern of taxon coverage, ∏_{ d }. This is the fraction of edges on the tree that are distinguished by subtrees induced by the taxon coverage pattern, averaged across all trees, or a random sample of trees. We examine these measures of partial decisiveness in empirical examples below.
Context 2: Taxon coverage is the product of a random sampling process
In this context the taxon coverage is not fixed but instead is determined by a random process of sampling. This might be a reasonable model for data assembled from partial genomic coverage, for example. The question we want to address is how many loci must be sampled to guarantee with some high probability that the taxon coverage will be decisive for some unknown tree T. In this context we need not require that it be decisive for all trees, because the results we derive will hold for any tree T and the true tree will be one of these. This is a subtle difference compared to the first context of a fixed taxon coverage pattern, where it made sense to require decisiveness simultaneously for most if not all trees. One way to understand this is to realize that because the sampling scheme is generating random patterns of taxon coverage, there is a chance that they are decisive for any particular true tree, and the probability calculations will be derived on this assumption. [If there were reason to think that two or more trees are "correct", and we want to require the random taxon coverage will define all of them, we might insist on decisiveness for all trees in this context too–see Discussion for possible examples.]
Let T be an unknown binary phylogenetic Xtree, and suppose we independently sample k random subsets S_{1},...,S_{ k }of X according to the following process: For each i, taxon x is included in S_{ i }with probability p(x), independently of other taxa. We will first consider two cases:

Uniform coverage (UC): p(x) = p for all taxa x in X;

Uniform (+1 complete) coverage( UC^{+1}): p(y) = 1 for one taxon y in X, and p(x) = p for all other taxa x in X. This is the case in which a reference taxon is present. Later we will discuss possible extensions of this model to allow variable taxon coverage probabilities across a tree.
Recall that taxon coverage pattern leads to a set of induced subtrees for each locus, TS_{1},...,TS_{ k }. Previously we assumed there was no conflict among these induced trees–they could be reconstructed accurately, but here we can obtain some results even if these induced trees are not fully resolved. We use a simple random model for lack of resolution in the induced trees. It does not explicitly model the possible sources of this loss of resolution, such as multiple hits in sequence data, but it provides a tractable starting point. Assume for some fixed probability q that each edge in TS_{ i }(for each i) is retained with probability q or collapsed with probability 1q, and this is done independently across the edges of TS_{ i }, and across the different induced trees. Let denote the resulting set of partiallyresolved induced subtrees of T. Let E_{ k }(T) be the event that defines T.
We are interested in computing bounds on the probability P(E_{ k }(T)) of this event that depend just on p, q and the total number of taxa n =X. It is useful to distinguish two further cases:

Rooted: T (and thereby the induced subtrees) are rooted trees;

Unrooted: T (and thereby the induced subtrees) are unrooted trees.
Let R(n) denote the set of rooted binary trees on n leaf taxa, and U(n) the set of unrooted binary trees on n leaf taxa. Note that by rooting a tree on an arbitrary taxon we can view U(n+1) as equivalent to R(n). Moreover, in what follows the UC condition applied to R(n) is essentially equivalent to the UC^{+1} condition applied to U(n+1).
Theorem 2 (i). For any rooted tree T in R(n), under UC (or an unrooted tree T in U(n+1) under UC^{+1}) we have: P(E_{ k }(T))≥ 1(n2) · (1p^{3}q)^{ k }. Consequently, for ε > 0, if .
Moreover, for any unrooted tree T in U(n+1) under UC, P(E_{ k }(T)) ≥ 1  ε holds if the term p^{3} in the above inequality is replaced by p^{4}.
(ii). Conversely, for any ε > 0, if then there exist (exponentially many) choices of T ∈ R(n) under UC (and trees in U(n+1) under UC or UC^{+1}) for which P(E_{ k }(T)) < 1  ε (See Appendix for proof).
Together these conditions give (respectively) upper and lower bounds on the number of loci needed for decisiveness relative to a single (unknown) tree. The precise number depends on the shape of the tree.
Remarks
When q = 1 the event E_{ k }(T) is equivalent to the event that S = {S_{1},...,S_{ k }}is phylogenetically decisive for T, but for q < 1 the event is a stronger condition. Also, Theorem 2 readily extends to allow the probabilities of taxon coverage to vary across the tree, albeit at the expense of more complicated formulae. However, there is one simple case where the same formulae apply: Part (i) holds as stated if, in the Uniform Coverage conditions, the equality p(x) = p for all taxa x is weakened to p(x) ≥ p for all x (similarly, part (ii) holds if we impose an inequality in the reverse direction).
We can also consider models that explicitly reflect how taxon coverage is often much more thorough within a clade than across the entire tree, by considering the following cladespecific condition:

Variable coverage ( VC ): p(x) ≥ p for all taxa x in clade C, and p(x) ≥ p' for any taxon x not in C (where p' <p).
Suppose that we have a rooted tree T in R(n), and for each of its n2 nontrivial clades, we independently generate at least r taxon sets, where each taxon is included or not in each taxon set indepedently with a probability p(x) that satisfies the VC condition for that clade. Repeat this independently for all the n2 clades of T (thus the resulting total number k of taxon sets must be at least r(n2)). Then the proof of Theorem 2(i) easily modifies to show that P(E_{ k }(T)) ≥ 1 ε provided that .
Empirical examples
We explore these theoretical results in several data sets derived from a diversity of phylogenomic data. We assess whether necessary and sufficient conditions for decisiveness for all trees are satisfied using our theoretical results. We also estimate the partial decisiveness quantities ∏_{ D }and ∏_{ d }using simulations, because no analytical results are known. We estimated these quantities for each data set using 1000 replicate simulated random trees generated using PAUP*'s equiprobable option [51]. Since the two measures of partial decisiveness can only be estimated on rooted trees, we estimate these quantities in a "reference subset" of the original data set. A reference subset is found by first identifying the taxon with the most loci present in the taxon coverage pattern, pruning all loci not present in this taxon from the entire matrix, and then removing any taxa that have been rendered dataless by this operation. This generates a data set with a reference taxon in the sense defined earlier, which can then be used to artificially "root" the analysis. Finally, we assess the bounds on data set size implied under random sampling of loci. In those analyses we assume the unrooted case without a reference taxon and require an error probability no greater than 5%.
Data mining examples
Taxon coverage and decisiveness parameters for five data sets.
Taxon  Source  Taxa  Loci  Coverage density  Fraction of triples  ∏_{ D }  ∏_{ d } 

Green plants  Driskell et al. 2004  69 (69)  251^{a} (220)  0.15 (0.15)  0.84 (0.84)  0.019  0.95 
Pleurodira  PhyLoTA ^{b}  44 (37)  15 (12)  0.26 (0.33)  0.37 (0.63)  0.000  0.82 
Amorpheae  PhyLoTA ^{b}  64 (64)  3 (3)  0.61 (0.61)  0.84 (0.84)  0.017  0.95 
Cactaceae  PhyLoTA ^{b}  488 (486)  18 (11)  0.14 (0.17)  0.22 (0.23)  0.000  0.61 
Metazoa  Hejnol et al., 2009  94 (94)  1487 (1351)  0.18 (0.19)  0.96 (0.95)  0.340  0.99 
However, the structure of the largest matrix, for Cactaceae is sobering. If we were to continue sampling loci with the same coverage density (14%) in the hopes of constructing a data set with high probability of being decisive for the true tree (Theorem 2, with perfect resolution), 2932 loci would be minimally necessary, and 23,891 would suffice. Of course, these results presuppose something like a uniform distribution of sampling effort across taxa, which is more likely to be found in genomic data sets than database mining studies. Nonetheless, these bounds are suggestive.
Finally, in this category we turn to several studies with very large taxon sets, for which we can only apply some of the simple results in Theorem 1. Goloboff et al. [4] assembled a supermatrix of 13 loci and 70,060 taxa. The locus with largest taxon coverage included ~20,000 taxa. By Theorem 1 the number of loci would have had to exceed 43 for phylogenetic decisiveness for all trees. Since most of the loci were sampled for far fewer than this number of taxa, a more accurate lower bound is probably higher.
This analysis, with coverage density of 8%, may be characteristic of taxon rich phylogenetic analyses published recently in which overall density is low, but a few loci have high coverage densities. For example, in [2] (2228 taxa of Leguminosae) and [3] (4954 taxa of Asterales) the locus with best coverage, nuclear ribosomal ITS, actually covered more than 73% and 84% of the label sets respectively, such that by Theorem 1 the lower bound on the number of loci is small, less than 3. However, in both cases, the ITS sequences were aligned by relatively "heroic" methods, and a more conservative procedure might have kept many separate taxonomically overlapping ITS alignments instead of one large one. We reported results from such an analysis (socalled "sparse" analysis in [2]) that ultimately included 1794 taxa in 72 loci. By Theorem 1 the lower bound for decisiveness for all trees in this case, assuming n/m = 1794/286, where 286 is the largest cluster (trnL), would have been 247 loci, far more than we had.
Partial genome sequences from shallow libraries
Hejnol et al. [10] constructed a large, sparse phylogenomic matrix of 94 taxa and 1487 loci derived mainly from ESTs. This is an expanded analysis of the widely cited phylogenomic study of 70 metazoan taxa [8]. Mean taxon coverage density is 18% across all loci, with the maximum density of 83% in one locus and a long tail of very low coverage loci. The given taxon coverage is not decisive for all trees, because it only includes 96% of possible triples (Table 1). Homo sapiens is present in 1351 of the loci. If we restrict the data to this reference subset, we can "root" the analysis with Homo as a reference taxon and test for a sufficient condition for decisiveness for all trees (a reference taxon plus all triples present: see above). For these loci there are still 94 taxa and 96% of possible triples are present. This fraction might have increased to 1 if the number of taxa in the new matrix had been less. This data set stands out by having a relatively high probability of decisiveness for random trees (34%), and a 99% chance of distinguishing a random edge in a tree. Under random taxon coverage (perfect resolution), by Theorem 2 the lower and upper bounds on loci are 1095 and 7148, which spans the actual data set.
Partial genome sequences from deep genomic libraries
Some genomic libraries are "deep", providing 10s to 100s of thousands of sequences, but still only partially cover a genome. Cranston et al. [11] analyzed data from BACend sequences from 10 species of Oryza including the complete genome of O. sativa and nine deep libraries averaging ~90,000 sequences per species. Collectively, these covered 44% of the CDS (coding regions) in O. sativa, but because of the random nature of library construction, of course, the pattern of overlap varied from species to species. The data were further filtered to include only single copy exonic regions with homology to O. sativa and having taxon coverage of at least four taxa, a total of 9481 exonic regions.
A posteriori the pattern of taxon coverage is trivially decisive for all trees because there are 6 loci (out of 9481) that have complete coverage–which automatically satisfies the fourway condition described earlier. An additional 60 loci are only missing one taxon. However, on average the taxon coverage among these exons averaged 48.3%. By Theorem 2 this value under random sampling and full resolution implies lower and upper bounds of 35 and 91 loci for decisiveness. However, if we consider partial resolution and if only one internal branch is present on average out of the 8 internal branches, the bounds increase to 292 and 742. Among a subset of loci most relevant for 6 of the species (the AA genome group) fully one quarter produced completely unresolved trees (the mean length of aligned exons was only 350 nt), so this is clearly an issue.
Complete genome sequences
Phylogenomic data sets constructed from complete genome sequences should suffer least from taxon coverage problems. We used Hahn et al.'s [32] analysis of gene families constructed from the 12 Drosophila complete genome sequences then available. Data were downloaded from Hahn's FRB Database http://www.indiana.edu/~hahnlab/fly/DfamDB/drosophila_frb.html and parsed to quantify taxon coverage. Loci for which there was at most a single copy of a gene present were retained. Of these 13,677 loci, 6979 had some taxa missing. The average taxon coverage was 66.3%. Note that this number is somewhat inflated because all loci that had only one taxon sampled (singletons) were excluded (Hahn et al. suggested many of these represent annotation errors). Clearly with 6698 loci having complete taxon coverage, there is no problem with phylogenetic decisiveness for all trees. Moreover, at this rate of taxon coverage on average, by Theorem 2, the number of randomly sampled loci needed for decisiveness for the true tree is bounded by 12  25 loci for 12 taxa, and even if resolution among trees were dismal (say one branch only out of 9 internal branches) we would have bounds of 127 and 246 loci. Thus if one were setting out to build a phylogenomic matrix with these data, this simple analysis suggests there are plenty of loci available. Interestingly the fractional taxon coverage is so high that even were one able to do this across 1000 species of Drosophila, the number of loci is still bounded by 26  47.
Results and Discussion
Impact of missing data on multilocus phylogenetic inference
Assessing the effect of missing data on phylogenetic inference has received substantial attention from phylogeneticists over many years (reviewed in [53]), but unambiguous conclusions have been elusive. For some amounts and patterns of missing data, simulation and empirical studies support the conclusion that its impact can be reduced to tolerable levels [37, 38]; for other cases, significant problems remain [39–41]. Properties of the data themselves, such as rates of evolution, stationarity, etc., can have a strong impact on the accuracy of inference, irrespective of the amount or pattern of missing data, so attempts to understand the problem at a comprehensive level taking all impacts into account have had poor predictive success.
Our results relate the taxon coverage–the amount and pattern of missing data–to the ability of data to define an unambiguous phylogenetic tree for all the taxa in the input, by separating the impact of missing data from tree reconstruction errors per se. To simplify the problem, we have defined a notion of phylogenetic decisiveness [44], which is separate from any notion of correctness of the individual gene trees based on individual loci.
We infer limits (bounds) on decisiveness. In both the context of a fixed pattern of taxon coverage and a random pattern generated by sampling, there are lower bounds to how many loci must be included for a given fraction of missing data. For the sampling case it is also possible to identify upper bounds so that it is clear when sufficient data are present. For a fixed (i.e., observed) taxon coverage pattern, it is possible to characterize whether it is decisive for some or all possible trees, but computing the latter is nontrivial.
Our results may help clarify why some arguments in the literature appear at odds with each other. On the one hand, Wiens [37] has suggested that it is not the fraction of missing data that is most important, but the overall amount of character data present. The most direct support from this in our results is seen in the random sampling context. According to Theorem 2 (i), for any fraction of missing data (1  p), there will be a number of loci that is sufficient to guarantee decisiveness. Another suggestive bit of evidence is seen in Table 1: Cactaceae and Metazoa have about the same fraction of missing data overall, but because Metazoa have so many more loci, the fraction of triples seen in that coverage pattern is much closer to the minimum level necessary for decisiveness (96% vs. 22%). A caveat is that the number of loci required can be quite large if the fraction of missing data is large or the number of taxa large (when p << 1 and or q << 1, the lower bound on number of loci k scales as 1/(p^{3}q) (expanding the expression Theorem 2 part (ii) in a Taylor series approximation)).
On the other hand, Hartmann et al. [39] and Lemmon et al. [41] discuss cases in which missing data are more problematic. Both studies simulated missing data at the sequence level occurring in blocks, which in our context we can think of as "loci". Both attributed errors in phylogenetic inference in part to these patterns of taxon coverage, often in combination with other factors. Our results on phylogenetic decisiveness for a fixed pattern of taxon coverage make it clear that some patterns of taxon coverage can be pathologically bad. Once we cannot count on random taxon coverage to ultimately sample all 4way partitions, it is easy to imagine patterns in which certain important 4way partitions are never seen and thus ambiguity remains a real possibility even as the number of loci gets arbitrarily large. The obvious case is just to suppose that we sample many loci but all of them either have {A, B, C, D} or {A, B, C, F} taxon coverage. This will never be decisive according to the 4way partition condition. Some of Hartmann et al.'s [39] ESTlike patterns of missing data are reminiscent of this case because of blocks of missing data at the ends in many alignments.
Implications in the context of gene tree discordance
Evidence is mounting that high levels of gene tree discordance is driven in part by real biological processes of lineage sorting, hybridization and recombination [14, 16, 17, 28]. Our analyses have by design set aside the latter issue, and by necessity the former, but incongruence is obviously real. Intuitively, since the lower bounds on the number of loci needed to infer a single tree decisively under random taxon coverage were derived assuming just one tree (see Appendix), we expect it would require many more loci to simultaneously be able to recover more than one tree. Interestingly in the Cranston et al. [11] analysis of rice phylogenomics, 118 out of 236 possible rooted 5taxon trees were found among the collection of gene trees (some undoubtedly due to error rather than lineage sorting). If lineage sorting is this rampant in a data set, we would need to know bounds not on decisiveness for one tree under random sampling but decisiveness for all trees. In an a posteriori analysis of a given pattern of coverage (as with rice), we can check decisiveness for all trees after the fact, and gene tree discordance provides a strong motivation for requiring this property.
Implications for nextgen sequencing strategies
Next generation sequence technologies are currently characterized by assembly of a large number of small reads. In phylogenomic analyses in which de novo assembly (without a reference genome) with limited coverage may be unavoidable, partial taxon coverage in the final data set may be the rule rather than the exception [54]. For example, the 1KP project [55] aims to use a shortread highthroughput approach to sequence transcriptomes of 1000 green plant species across 1 billion years of evolutionary history. A combination of partial transcriptome coverage in the final assemblies and lack of homology across such deep divergences will certainly limit the overall taxon coverage densities. Given estimates of this fraction derived from pilot studies it should be possible to use Theorem 2 to predict how much of the transcriptome is needed to anticipate a decisive analysis.
Does conflict hurt or help?
A particular consequence of Inequality (1) is that a binary phylogenetic tree T for a large number of species cannot be defined by a small number of poorly resolved trees. For example, if T has 1000 leaves, then we require a total of 997 or more internal edges amongst any set of trees that defines T. Now suppose we use a supertree method like MRP (Matrix Representation with Parsimony, see review in [25]) on these input trees. An interesting question is whether Inequality (1) might still hold if we replace the condition that the trees "T_{1},...,T_{ k }define T" by "T_{1},...,T_{ k }have a unique MRP tree T".
This seems plausible since the two notions are equivalent when the input trees T_{1},...,T_{ k }are compatible, and it might be supposed that conflictfree data is in some sense 'optimal'. However, it turns out that such a modification to Inequality (1) fails. For example, there is an incompatible set of six trees for 12 taxa, each with just one internal edge, that has a unique MRP tree T on those 12 taxa (results not shown); by contrast if we wished to additionally require the input trees to be compatible, we would require at least nine trees (each having one internal edge), not six. In other words, allowing the trees to conflict can in some cases actually 'help' reduce the amount of data required to determine a unique supertree by a method such as MRP.
Relation to phylogenetic "groves"
Future work
Several problems are suggested by this work and its application to empirical studies. First, many real taxon coverage patterns are sparse enough that they are not decisive for all trees. It would be useful to know what data should be obtained to convert an indecisive to a decisive data set, or to convert one that is decisive with low probability to one that is decisive with high probability. In both cases there are natural optimality problems that could be stated that would try to minimize the cost in additional sequencing. Second, along the lines of phylogenetic groves, when a taxon coverage pattern is not decisive for the entire label set, it might well be decisive for a large subset of it. What is the size of the largest label set for which a given pattern of taxon coverage is decisive? In some sense this conveys the "effective tree size" which a taxon coverage pattern can address.
Conclusions
We have characterized the impact of partial taxon coverage in multilocus phylogenomic data sets in the context of both fixed and random patterns of coverage. For fixed coverage patterns we have determined lower bounds on the number of loci needed for a coverage pattern to be phylogenetically decisive for all trees–that is to allow reconstruction of a unique tree on all taxa irrespective of what the tree is. Most real data sets having reasonably large numbers of taxa are not decisive for all trees, but in these cases we developed methods for assessing the fraction of trees for which they are decisive and the fraction of edges in those trees that can be distinguished. In general, the first measure is still low, but the second can be quite high. For random taxon coverage patterns, such as might be seen in data sets derived from EST or BAC libraries, upper and lower bounds on the number of loci needed for decisive tree inference were obtained. Some real data sets required relatively few loci for decisiveness, particularly those with few taxa and modest levels of missing data, whereas other data sets required tens of thousands of loci when the number of taxa and fraction of missing data was high. Clearly, the pattern of taxon coverage is a factor in large scale phylogenetic inference, regardless of the details of the sequence data themselves or tree reconstruction method.
Appendix: Proofs
from which the results stated now follow.
since a tree in R(n) has precisely n  2 internal edges, as required. The resulting bound involving k now follows by routine algebra.
Turning to the case of an unrooted tree T ∈ U(n + 1), if condition UC^{+1} holds then we can directly apply the result for rooted trees under UC. However if the weaker condition UC holds for T ∈ U(n + 1)then we use the result (from [47]) that T is defined by any set of quartet trees that distinguishes every edge of T provided that each quartet contains a fixed leaf. We can then repeat the above argument but at the penalty of increasing the exponent on p from 3 to 4.
from which part (ii) follows upon applying the inequality: This completes the proof.
Declarations
Acknowledgements
This research was supported by a US NSF grant to MJS and MMM. MS thanks the Allan Wilson Centre for Molecular Ecology and Evolution, and the NZ Royal Society (James Cook Fellowship). We thank Tandy Warnow for discussion.
Authors’ Affiliations
References
 Kallersjo M, Albert VA, Farris JS: Homoplasy increases phylogenetic structure. Cladistics. 1999, 15 (1): 9193.Google Scholar
 McMahon MM, Sanderson MJ: Phylogenetic supermatrix analysis of GenBank sequences from 2228 papilionoid legumes. Syst Biol. 2006, 55: 818836. 10.1080/10635150600999150.View ArticlePubMedGoogle Scholar
 Smith SA, Beaulieu JM, Donoghue MJ: Megaphylogeny approach for comparative biology: an alternative to supertree and supermatrix approaches. BMC Evolutionary Biology. 2009, 9: 3710.1186/14712148937.PubMed CentralView ArticlePubMedGoogle Scholar
 Goloboff PA, Catalano SA, Mirande JM, Szumik CA, Arias JS, Kallersjo M, Farris JS: Phylogenetic analysis of 73 060 taxa corroborates major eukaryotic groups. Cladistics. 2009, 25 (3): 211230. 10.1111/j.10960031.2009.00255.x.View ArticleGoogle Scholar
 Rokas A, Williams B, King N, Carroll S: Genomescale approaches to resolving incongruence in molecular phylogenies. Nature. 2003, 425: 798804. 10.1038/nature02053.View ArticlePubMedGoogle Scholar
 Driskell AC, Ané C, Burleigh JG, McMahon MM, O'Meara B, Sanderson MJ: Prospects for building the tree of life from large sequence databases. Science (Washington D C). 2004, 306: 11721174. 10.1126/science.1102036.View ArticleGoogle Scholar
 Wildman D, Uddin M, Opazo J, Liu G, Lefort V, Guindon S, Gascuel O, Grossman L, Romero R, Goodman M: Genomics, biogeography, and the diversification of placental mammals. Proceedings of the National Academy of Sciences. 2007, 104 (36): 1439514400. 10.1073/pnas.0704342104.View ArticleGoogle Scholar
 Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, et al: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452 (7188): 745U745. 10.1038/nature06614.View ArticlePubMedGoogle Scholar
 Zou XH, Zhang FM, Zhang JG, Zang LL, Tang L, Wang J, Sang T, Ge S: Analysis of 142 genes resolves the rapid diversification of the rice genus. Genome Biology. 2008, 9 (3): 10.1186/gb200893r49.Google Scholar
 Hejnol A, Obst M, Stamatakis A, Ott M, Rouse GW, Edgecombe GD, Martinez P, Baguna J, Bailly X, Jondelius U, et al: Assessing the root of bilaterian animals with scalable phylogenomic methods. Proceedings of the Royal Society BBiological Sciences. 2009, 276 (1677): 42614270. 10.1098/rspb.2009.0896.PubMed CentralView ArticleGoogle Scholar
 Cranston KA, Hurwitz B, Sanderson MJ, Ware D, Wing RA, Stein L: Phylogenomic analysis from deep BACend sequence libraries in rice. Syst Bot. 2010Google Scholar
 Stamatakis A, Ott M: Efficient computation of the phylogenetic likelihood function on multigene alignments and multicore architectures. Philosophical Transactions of the Royal Society BBiological Sciences. 2008, 363 (1512): 39773984. 10.1098/rstb.2008.0163.PubMed CentralView ArticleGoogle Scholar
 Galtier N, Daubin V: Dealing with incongruence in phylogenomic analyses. Philosophical Transactions of the Royal Society BBiological Sciences. 2008, 363 (1512): 40234029. 10.1098/rstb.2008.0144.PubMed CentralView ArticleGoogle Scholar
 Pollard DA, Iyer VN, Moses AM, Eisen MB: Widespread discordance of gene trees with species tree in Drosophila: Evidence for incomplete lineage sorting. Plos Genetics. 2006, 2 (10): 16341647.Google Scholar
 Cranston KA, Hurwitz B, Ware D, Stein L, Wing RA: Species Trees from Highly Incongruent Gene Trees in Rice. Syst Biol. 2009, 58 (5): 489500. 10.1093/sysbio/syp054.View ArticlePubMedGoogle Scholar
 White MA, Ane C, Dewey CN, Larget BR, Payseur BA: FineScale Phylogenetic Discordance across the House Mouse Genome. Plos Genetics. 2009, 5 (11): 10.1371/journal.pgen.1000729.Google Scholar
 Ebersberger I, Galgoczy P, Taudien S, Taenzer S, Platzer M, von Haeseler A: Mapping human genetic ancestry. Mol Biol Evol. 2007, 24 (10): 22662276. 10.1093/molbev/msm156.View ArticlePubMedGoogle Scholar
 Kubatko LS, Degnan JH: Inconsistency of phylogenetic estimates from concatenated data under coalescence. Syst Biol. 2007, 56 (1): 1724. 10.1080/10635150601146041.View ArticlePubMedGoogle Scholar
 Liu L, Yu LL, Pearl DK: Maximum tree: a consistent estimator of the species tree. Journal of Mathematical Biology. 2010, 60 (1): 95106. 10.1007/s0028500902600.View ArticlePubMedGoogle Scholar
 Ren FR, Tanaka H, Yang ZH: A likelihood look at the supermatrixsupertree controversy. Gene. 2009, 441 (12): 119125PubMedGoogle Scholar
 Rannala B, Yang ZH: Phylogenetic inference using whole Genomes. Annual Review of Genomics and Human Genetics. 2008, 9: 217231. 10.1146/annurev.genom.9.081307.164407.View ArticlePubMedGoogle Scholar
 Edwards SV: Is a new and general theory of phylogenetic inference emerging?. Evolution. 2009, 63 (1): 119. 10.1111/j.15585646.2008.00549.x.View ArticlePubMedGoogle Scholar
 Degnan JH, Rosenberg NA: Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol Evol. 2009, 24 (6): 332340. 10.1016/j.tree.2009.01.009.View ArticlePubMedGoogle Scholar
 Knowles LL: Estimating Species Trees: Methods of Phylogenetic Analysis When There Is Incongruence across Genes. Syst Biol. 2009, 58 (5): 463467. 10.1093/sysbio/syp061.View ArticlePubMedGoogle Scholar
 BinindaEmonds ORP: The evolution of supertrees. Trends Ecol Evol. 2004, 19 (6): 315322. 10.1016/j.tree.2004.03.015.View ArticlePubMedGoogle Scholar
 Maddison WP, Knowles LL: Inferring phylogeny despite incomplete lineage sorting. Syst Biol. 2006, 55 (1): 2130. 10.1080/10635150500354928.View ArticlePubMedGoogle Scholar
 Edwards SV, Liu L, Pearl DK: Highresolution species trees without concatenation. Proc Natl Acad Sci USA. 2007, 104 (14): 59365941. 10.1073/pnas.0607004104.PubMed CentralView ArticlePubMedGoogle Scholar
 Siepel A: Phylogenomics of primates and their ancestral populations. Genome Res. 2009, 19 (11): 19291941. 10.1101/gr.084228.108.PubMed CentralView ArticlePubMedGoogle Scholar
 Kubatko LS: Identifying Hybridization Events in the Presence of Coalescence via Model Selection. Syst Biol. 2009, 58 (5): 478488. 10.1093/sysbio/syp055.View ArticlePubMedGoogle Scholar
 Meng C, Kubatko LS: Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: A model. Theoretical Population Biology. 2009, 75 (1): 3545. 10.1016/j.tpb.2008.10.004.View ArticlePubMedGoogle Scholar
 Sanderson MJ, Driskell AC, Ree RH, Eulenstein O, Langley S: Obtaining maximal concatenated phylogenetic data sets from large sequence databases. Mol Biol Evol. 2003, 20: 10361042. 10.1093/molbev/msg115.View ArticlePubMedGoogle Scholar
 Hahn MW, Han MV, Han SG: Gene family evolution across 12 Drosophila genomes. Plos Genetics. 2007, 3 (11): 21352146. 10.1371/journal.pgen.0030197.View ArticleGoogle Scholar
 Soderlund C, Descour A, Kudrna D, Bomhoff M, Boyd L, Currie J, Angelova A, Collura K, Wissotski M, Ashley E, et al: Sequencing, Mapping, and Analysis of 27,455 Maize FullLength cDNAs. Plos Genetics. 2009, 5 (11): 10.1371/journal.pgen.1000740.Google Scholar
 Donoghue MJ, Doyle JA, Gauthier J, Kluge AG, Rowe T: The Importance of Fossils in Phylogeny Reconstruction. Annual Review of Ecology and Systematics. 1989, 20: 431460. 10.1146/annurev.es.20.110189.002243.View ArticleGoogle Scholar
 Huelsenbeck JP: When are fossils better than extant taxa in phylogenetic analysis?. Syst Zool. 1991, 40 (4): 458469. 10.2307/2992240.View ArticleGoogle Scholar
 Kearney M: Fragmentary taxa, missing data, and ambiguity: Mistaken assumptions and conclusions. Syst Biol. 2002, 51 (2): 369381. 10.1080/10635150252899824.View ArticlePubMedGoogle Scholar
 Wiens JJ: Missing data, incomplete taxa, and phylogenetic accuracy. Syst Biol. 2003, 52 (4): 528538. 10.1080/10635150390218330.View ArticlePubMedGoogle Scholar
 Philippe H, Snell E, Bapteste E, Lopez P, Holland P, Casane D: Phylogenomics of eukaryotes: Impact of missing data on large alignments. Mol Biol Evol. 2004, 21: 17401752. 10.1093/molbev/msh182.View ArticlePubMedGoogle Scholar
 Hartmann S, Vision TJ: Using ESTs for phylogenomics: Can one accurately infer a phylogenetic tree from a gappy alignment?. BMC Evolutionary Biology. 2008, 8: 9510.1186/14712148895.PubMed CentralView ArticlePubMedGoogle Scholar
 Campbell V, Lapointe FJ: The Use and Validity of Composite Taxa in Phylogenetic Analysis. Syst Biol. 2009, 58 (6): 560572. 10.1093/sysbio/syp056.View ArticlePubMedGoogle Scholar
 Lemmon AR, Brown JM, StangerHall K, Lemmon EM: The Effect of Ambiguous Data on Phylogenetic Estimates Obtained by Maximum Likelihood and Bayesian Inference. Syst Biol. 2009, 58 (1): 130145. 10.1093/sysbio/syp017.View ArticlePubMedGoogle Scholar
 Philippe H, Delsuc F, Brinkmann H, Lartillot N: Phylogenomics. Annual Review of Ecology Evolution and Systematics. 2005, 36: 541562. 10.1146/annurev.ecolsys.35.112202.130205.View ArticleGoogle Scholar
 de Queiroz A, Gatesy J: The supermatrix approach to systematics. Trends Ecol Evol. 2007, 22 (1): 3441. 10.1016/j.tree.2006.10.002.View ArticlePubMedGoogle Scholar
 Steel M, Sanderson MJ: Characterizing phylogenetically decisive taxon coverage. Applied Mathematics Letters. 2010, 23 (1): 8286. 10.1016/j.aml.2009.08.009.View ArticleGoogle Scholar
 Goloboff PA: Homoplasy and the choice among cladograms. Cladisticsthe International Journal of the Willi Hennig Society. 1991, 7 (3): 215232. 10.1111/j.10960031.1991.tb00035.x.View ArticleGoogle Scholar
 Semple C, Steel M: Phylogenetics. 2003, New York: Oxford University PressGoogle Scholar
 Steel M: The complexity of reconstructing trees from qualitative characters and subtrees. J Classif. 1992, 9: 91116. 10.1007/BF02618470.View ArticleGoogle Scholar
 Humphries PJ: Combinatorial aspects of leaflabelled trees. 2008, Christchurch, New Zealand: University of CanterburyGoogle Scholar
 Böcker S, Bryant D, Dress AWM, Steel MA: Algorithmic aspects of tree amalgamation. J Algorithms. 2000, 37 (2): 522537. 10.1006/jagm.2000.1116.View ArticleGoogle Scholar
 Bryant D, Steel M: Extension operations on sets of leaflabelled trees. Advances in Applied Mathematics. 1995, 16 (4): 425453. 10.1006/aama.1995.1020.View ArticleGoogle Scholar
 Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). 2002, Sunderland, MA: Sinauer, 4Google Scholar
 Sanderson MJ, Boss D, Chen D, Cranston KA, Wehe A: The PhyLoTA Browser: Processing GenBank for Molecular Phylogenetics Research. Syst Biol. 2008, 57: 335346. 10.1080/10635150802158688.View ArticlePubMedGoogle Scholar
 Wiens JJ: Missing data and the design of phylogenetic analyses. Journal of Biomedical Informatics. 2006, 39 (1): 3442. 10.1016/j.jbi.2005.04.001.View ArticlePubMedGoogle Scholar
 Gibbons JG, Janson EM, Hittinger CT, Johnston M, Abbot P, Rokas A: Benchmarking NextGeneration Transcriptome Sequencing for Functional and Evolutionary Genomics. Mol Biol Evol. 2009, 26 (12): 27312744. 10.1093/molbev/msp188.View ArticlePubMedGoogle Scholar
 The 1kp Project. [http://www.onekp.com/essential.html]
 Ane C, Eulenstein O, PiaggioTalice R, Sanderson MJ: Groves of phylogenetic trees. Annals of Combinatorics. 2009, 13: 139167. 10.1007/s000260090017x.View ArticleGoogle Scholar
 Sanderson MJ, Ane C, Eulenstein O, FernandezBaca D, Kim J, McMahon MM, PiaggioTalice R: Fragmentation of large data sets in phylogenetic analysis. Reconstructing evolution: new mathematical and computational advances. Edited by: Gascuel O, Steel M. 2007, Oxford: Oxford University PressGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.