Using supernetworks to distinguish hybridization from lineage-sorting
© Holland et al; licensee BioMed Central Ltd. 2008
Received: 16 December 2007
Accepted: 14 July 2008
Published: 14 July 2008
A simple and widely used approach for detecting hybridization in phylogenies is to reconstruct gene trees from independent gene loci, and to look for gene tree incongruence. However, this approach may be confounded by factors such as poor taxon-sampling and/or incomplete lineage-sorting.
Using coalescent simulations, we investigated the potential of supernetwork methods to differentiate between gene tree incongruence arising from taxon sampling and incomplete lineage-sorting as opposed to hybridization. For few hybridization events, a large number of independent loci, and well-sampled taxa across these loci, we found that it was possible to distinguish incomplete lineage-sorting from hybridization using the filtered Z-closure and Q-imputation supernetwork methods. Moreover, we found that the choice of supernetwork method was less important than the choice of filtering, and that count-based filtering was the most effective filtering technique.
Filtered supernetworks provide a tool for detecting and identifying hybridization events in phylogenies, a tool that should become increasingly useful in light of current genome sequencing initiatives and the ease with which large numbers of independent gene loci can be determined using new generation sequencing technologies.
In recent years there has been growing interest in the problem of building explicit models of reticulate evolution [1–12]. This work has to a large part been motivated by biological research highlighting the potential importance of hybridization in the origin of biotic diversity, biological invasion and rapid adaptation [13–27].
One simple and widely used approach for detecting hybridization has been to compare gene trees built from independent gene loci, and to consider gene tree incongruence as evidence for hybridization [1, 28, 29]. However, hybridization is not the only possible cause of gene tree incongruence. Other explanations include phylogenetic error , unrecognised gene duplication and loss , incomplete lineage-sorting  and lateral gene transfer .
In light of current genome sequencing efforts and the ease of sequencing large numbers of independent gene loci using new generation sequencing technologies, it is important to find ways to differentiate between various explanations of gene tree incongruence. Here we focus on distinguishing hybridization from incomplete lineage-sorting. In this regard, a helpful concept might be "principal trees", which are the trees displayed by a hybridization network (see the subsection Simulations for a more formal definition of principal trees). If a phylogeny contains no hybridization or lateral gene transfer, then the expectation is for one principal or "species" tree. However, if hybridization has occurred, then there will be multiple principal trees. In the absence of incomplete lineage-sorting, each principle tree will represent the evolutionary history for a large collection of loci, but where incomplete lineage-sorting occurs gene trees may differ from their underlying principal tree.
Here we investigate the potential of filtered Z-closure  and Q-imputation  supernetworks to distinguish phylogenetic signals arising from principal trees from phylogenetic signals caused by a combination of incomplete taxon-sampling and incomplete lineage-sorting in evolutionary histories involving hybridization. We consider these methods since they are not only designed to cope with conflicting phylogenetic signals, but also with data that current genome sequencing efforts can fail to gather (for example, for multi-gene data sets, such as those generated from expressed sequence tag (EST) databases, gene sequences are often missing for species of interest). In particular, using gene trees generated under a coalescent process, we test whether these methods can be used to filter out phylogenetic signals that do not correspond to edges in principal trees, signals that have the potential of confounding efforts to reconstruct complex phylogenies.
Analogous to the supertree framework  our input is a set of trees on overlapping but not necessarily identical taxa. We refer to the complete taxa set as the union of the taxa sets of the input trees; complete splits are bipartitions of the complete taxa set and trivial splits are splits where one part consists of precisely one element. Furthermore partial trees and partial splits are trees and splits on a subset of the complete taxa set. We denote a split of the taxa set X as A|B where A and B are both subsets of X, A ∪ B = X, and A ∩ B is the empty set.
Our overall approach is to first generate a collection of partial trees along a hybridization network in the presence of incomplete lineage-sorting (modelled by the coalescent process), to then apply a supernetwork method to this collection of partial trees to obtain a collection of complete splits, and to then apply a filter to reduce the complexity of this collection. The reduced collection of complete splits is then compared to the splits associated with the hybridization network to determine if they have been accurately recovered. We use this approach to study two supernetwork methods – Z-closure  and Q-imputation , and two types of filter – one counting-based and one homoplasy-based .
Note that the output of Z-closure is dependent on the order of elements in the list of partial splits, and so we repeat the procedure for 10 random orderings keeping a cumulative count of how many times each complete split appears. (Simulations indicate that this order dependence is not strong , so there would be little benefit in performing a larger number of random orderings.) Also note that the Z-closure implementation used for this paper differs slightly from that in SplitsTree4 in that it keeps track of multiple copies of partial splits and complete splits, as this information is required by the counting filter that we apply later.
The other supernetwork method we consider, Q-imputation , also uses partial trees as input but uses an alternative approach to generating complete splits that is based on the four-taxon subtrees (quartets) of the partial trees. For each partial tree with missing taxa – that is, taxa that are in the complete taxa set but not in the taxa set for that tree – the missing taxa are inserted in the tree. This is done by processing the missing taxa in a fixed order and placing each taxon within the partial tree to maximise the number of quartets that contain the missing taxon and are also quartets of the other partial input trees. Once all the trees have been completed the list of complete splits displayed by the completed trees is returned. (In the special case where all the trees are on identical taxa sets, the Q-imputation method reduces to the consensus network method ).
We apply two different kinds of filter to the lists of complete splits obtained from the two supernetwork methods, a homoplasy-based filter and a counting-based filter. The homoplasy filter  requires two user-defined parameters x and y. The level of homoplasy for each complete split and partial tree is determined by recoding the split as a binary character, reducing it to the same taxa set as the partial tree, and evaluating the number of character state changes required to explain the character on the partial tree (i.e. the parsimony score). Splits that have a parsimony score greater than a given number x in more than a given number y of the partial trees are filtered out. The counting filter has one user-defined parameter n and keeps the n splits that appear most frequently in the list of complete splits (ties are broken arbitrarily). Note that for Q-imputation this is equivalent to the filter described in  for some choice of threshold.
Accuracy of the supernetwork methods was determined by counting the number of false positive and false negative splits. A false positive is a split that is displayed by the supernetwork, but is not a true split of the hybridization network; a false negative is a split that is displayed by at least one of the principal trees of the hybridization network, but not displayed by the supernetwork. Note that these definitions differ from those used in  in that they measure accuracy with respect to an underlying hybridization network that has been used to generate the data.
Hybridization networks used in simulations.
Principle trees (given in Newick format)
We applied 4 different homoplasy filters to the lists of complete splits returned by Z-closure and Q-imputation:
• (HF1) keep only splits with no homoplasy (i.e. a parsimony score of 1) on all partial trees,
• (HF2) keep only splits with no homoplasy on 75% or more of the partial trees,
• (HF3) keep only splits with no homoplasy on 50% or more of the partial trees,
• (HF4) keep splits with a parsimony score of 2 or less on all partial trees.
We also applied the counting filter (CF) to both Z-closure and Q-imputation. For each hybridization network we selected n splits, where n was fixed to be the number of unique non-trivial splits associated with the principal trees of the network (Table 1).
Results and Discussion
False positives for hybridization network 7 using the counting filter to select the 8 highest weight splits.
m = 0
m = 1
m = 2
m = 3
m = 0
m = 1
m = 2
m = 3
b = ∞
g = 5
g = 5
g = 10
g = 10
g = 15
g = 15
g = 20
g = 20
m = 0
m = 1
m = 2
m = 3
m = 0
m = 1
m = 2
m = 3
b = 1
g = 5
g = 5
g = 10
g = 10
g = 15
g = 15
g = 20
g = 20
m = 0
m = 1
m = 2
m = 3
m = 0
m = 1
m = 2
m = 3
b = 0.5
g = 5
g = 5
g = 10
g = 10
g = 15
g = 15
g = 20
g = 20
Figure 4 shows the change in the average number of split false positives and false negatives with respect to the number of gene trees. The results are averaged over 100 repetitions and the 12 combinations of number of missing taxa m and coalescent branch lengths b.
As can be seen in Figure 4, the (HF1) filter is far too stringent in combination with either Z-closure or Q-imputation; it gives almost no false positives but false negatives increase with increasing g towards the maximum value of 8. The (HF4) filter is not stringent enough in combination with either Z-closure or Q-imputation; it gives almost no false negatives but false positives increase with increasing g. Moreover, (HF3) is ineffective in combination with Z-closure as the number of false positives increases with increasing number of partial trees, in combination with Q-imputation the average number of false positives stays close to 2 for all values of g. (HF2) is the most effective of the homoplasy filters, as for both Z-closure and Q-imputation both types of errors either decrease or stay reasonably constant with increasing g, a property that we would expect any filtered supernetwork method to satisfy. The counting filter also displays this property for both Z-closure and Q-imputation; both false positives and false negatives decrease with increasing number of input trees.
Figure 5 is similar to Figure 4 except that rather than averaging over all values of b and m we focus on the difficult case with the highest number of missing taxa and the most incongruence generated by incomplete lineage-sorting (m = 3, b = 0.5). While all the filtered supernetwork methods shown in Figure 5 control false negatives, false positives increase with increasing g for both Z-closure and Q-imputation using (HF3).
For the rest of this section, we restrict our attention to the best homoplasy filter (HF2) and the counting filter (CF).
Figure 6 shows the trends in the number of false positives and false negatives as the number of missing taxa per tree, m, increases from 0 to 3. Results are averaged over the 12 possible settings for number of partial trees g and coalescent branch lengths b. The (HF2) and (CF) filters exhibit very different behaviour. For (HF2) as m increases the number of false positives increases, in particular going from 2 to 3 missing taxa produces a dramatic increase. Conversely the number of false negatives decreases, presumably due to the fact that as the number of missing taxa gets large more splits meet the requirement of the filter. This effect was not observed for (CF) where the total number of splits is capped; here both false positives and negatives increase with growing m.
Recall that the parameter b affects the probability that the trees generated by COAL  will match the principal tree sampled from the hybridization network, b = 8 corresponds to trees that match exactly. Figure 7 shows the trends in the number of false positives and false negatives for different values of b. Results are averaged over the 16 possible settings for the number of partial trees g and the number of missing taxa m. As expected, both methods and filters perform better when b is large.
Table 2 shows the number of false positives for hybridization network 7 (Figure 2a), which has 2 hybridization events and 8 true splits, for m = 0, 1, 2 or 3, b = 0.5, 1 or 8, and g = 5, 10, 15 or 20 for both Z-closure and Q-imputation with (CF). If two out of three of the conditions (g, m, or b) are favourable (i.e. many input trees, few missing taxa, and the probability that the input trees are congruent with the principal trees is high) then both methods work well. However if two or three of the conditions are unfavourable then both methods start to break down.
We have evaluated the potential of Z-closure and Q-imputation filtered supernetworks to identify splits belonging to the sets of principal trees associated with hybridization networks. We have found that this approach can recover these splits when there are few hybridization events. However, our results imply that (1) if gene trees have many missing taxa then many gene trees are required; (2) if the gene trees are frequently incongruent with the principal trees of the hybridization network due to incomplete lineage-sorting then a large number of near complete gene trees is required; (3) and if there are few gene trees available they need to be both near complete and highly congruent with the principal trees.
In our simulations the counting filter picked the n best-supported splits, where n was chosen to be the known number of true non-trivial splits. Of course with real data n will not be known, although in practice n could be chosen by, for example, greedily introducing splits with highest support as long as the corresponding network does not become too complex to easily interpret. Approaches to do this are described in  for consensus networks. Note that by increasing n the risk of introducing false positive splits is increased, although the risk of failing to identify true positive splits is reduced.
Despite these limitations, with the potential now of obtaining large numbers of splits from independent gene loci using new generation sequencing technologies, our findings may nevertheless be applicable for tree-like phylogenies where some degree of hybridization is inferred . In such cases, filtered supernetworks can be used to identify the true splits of the underlying hybridization network. Once these are obtained, the method of  can be used to convert the split system into a hybridization scenario.
One of our most interesting findings is that the choice of whether to use Z-closure or Q-imputation seems to have much less impact on accuracy with regards to recovering the splits in the underlying hybridization network than the choice of filter. For both Q-imputation and Z-closure the counting filter (CF) has the desirable property that as the amount of data increases (more genes or more complete gene trees) the rate of both false positives and false negatives goes down. Several settings were tried for the homoplasy-based filter (HF1 – HF4). HF1 was too stringent, and HF3/HF4 tended to either suffer from increasing false positives or increasing false negatives as the number of gene trees increased. HF2 gave the best compromise between these extremes.
Using the HF2 filter, we found that Z-closure had a higher false positive rate than Q-imputation over a range of parameter combinations (Figure 8A). One explanation might be that Z-closure can potentially generate more splits than Q-imputation. For example, given g fully resolved gene trees on 8-m taxa (m = 1, 2, 3), Q-imputation can generate at most 5*g non-trivial complete splits, whereas Z-closure can produce at most 10*(5-m)*g non-trivial complete splits (where 10 is the number of random orderings). Hence the maximum number of splits that Z-closure could generate decreases as m grows, whereas the number of splits that Q-imputation could generate stays constant. What we observe for both methods and filters is that false positives increase with increasing m (Figure 6). Therefore, the maximum number of splits that Z-closure and Q-imputation could generate does not appear to explain the difference in false positive rates. We think a more likely explanation is that Q-imputation places missing taxa in such a way as to maximize agreement with the input trees, hence tending to produce multiple copies of the same splits. Conversely, Z-closure aims to find all possible complete splits that can be derived by extending partial splits using the Z-closure rule, a process that can yield many different splits. Hence we expect that Q-imputation would be likely to generate fewer false positives than Z-closure in general. This difference is not greatly reduced by HF2 as, in contrast to CF, it does not place a cap on the total number of splits.
We found that HF2 resulted in more false negatives than CF (Figure 8B). This may be due to the fact that this filter only selects splits that have no homoplasy when restricted to 75% of the input trees. Since the principal trees are obtained from a network, they can be different and in some cases may only agree on a small number of edges. Even a true split may have a high homoplasy score when restricted to a particular principal tree. In contrast, CF only selects splits that occur with high frequency, irrespective of whether they are in agreement with any of the input trees.
Although all the trees used in our simulations were fully resolved, both supernetwork methods considered here can be applied to partially resolved trees. Thus, when inferring gene trees to be used as input to a supernetwork method, it would probably be a reasonable approach to only retain those edges in the estimated gene trees that have high support (e.g. bootstrap support or posterior probability higher than some cut-off value).
In cases where there are many hybridization events, especially between individuals that are not closely related, there will be many principal trees and corresponding splits (as in hybridization network 10). Many of these splits will occur at low frequencies making them hard to distinguish from phylogenetic error. This means that phylogenetic inference will be limited, as gene-tree incongruence will be extensive. In such cases, rather than attempt to reconstruct a hybridization network, it may be more appropriate to formulate objective tests to better understand the complexity of the data and the extent to which hybridization contributes to this complexity. Joly, McLenachan and Lockhart (submitted manuscript) have recently proposed such a test.
An unexplored idea worthy of study is the investigation of model-based, rather than combinatorial, methods of filtering. One approach might be to consider posterior probability distributions on species trees . It will be interesting to investigate whether such posterior distributions can also be analysed for evidence of distinct principal trees in cases where evolutionary relationships are complex.
All authors thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, UK for hosting them in the context of their Phylogenetics Programme where part of the work presented in this paper was carried out. BRH and PJL thank Marsden Fund MAU0509 and MAU0510, VM thanks the Engineering and Physical Sciences Research Council (EPSRC), grant EP/D068800/1.
- Linder CR, Rieseberg L: Reconstructing patterns of reticulate evolution in plants. Amer J Botany. 2004, 91: 1700-1708.View ArticleGoogle Scholar
- Gusfield D, Eddhu S, Langley CH: Optimal, efficient reconstruction of phylogenetic networks with constrained recombination. J Bioinform Comput Biol. 2004, 2 (1): 173-213.View ArticlePubMedGoogle Scholar
- Nakhleh L, Warnow T, Linder CR, St. John K: Reconstructing reticulate evolution in species–theory and practice. J Comput Biol. 2005, 12 (6): 796-811.View ArticlePubMedGoogle Scholar
- Jansson J, Sung W: Inferring a level-1 phylogenetic network from a dense set of rooted triplets. Springer Lecture Notes in Computer Science. 2004, 3106: 462-471.View ArticleGoogle Scholar
- Huson DH, Kloepper T, Lockhart PJ, Steel M: Reconstruction of reticulate networks from gene trees. RECOMB 2005 LNBI 3500. Edited by: Miyano S. 2005, Springer-Verlag Heidelberg, 233-249.Google Scholar
- Song YS, Hein J: Constructing minimal ancestral recombination graphs. J Comput Biol. 2005, 12 (2): 147-169.View ArticlePubMedGoogle Scholar
- Baroni M, Semple C, Steel M: Hybrids in real time. Syst Biol. 2006, 55: 46-56.View ArticlePubMedGoogle Scholar
- Semple C, Steel M: Phylogenetics. 2003, Oxford: Oxford University PressGoogle Scholar
- Huber KT, Oxelman B, Lott M, Moulton V: Reconstructing the evolutionary history of polyploids from multi-labelled trees. Mol Biol Evo. 2006, 23: 1784-1791.View ArticleGoogle Scholar
- Willson S: Unique determination of some homoplasies at hybridization events. Bull Math Biol. 2007, 69: 1709-1725.View ArticlePubMedGoogle Scholar
- Cardona G, Rosselló F, Valiente G: Tripartitions do not always discriminate Phylogenetic Networks. Math Biosci. 211: 356-370.Google Scholar
- van Iersel L, Keijsper J, Kelk S, Stougie L: Constructing level-2 phylogenetic networks from triplets. 2007, arXiv:0707.2890v1 [q-bio.PE]Google Scholar
- Linnaeus C: Species plantarum, exhibentes plantas rite cognitas, ad genera relatas, cum differentiis specificis, nominibus trivialibus, synonymis selectis, locis natalibus, secundum systema sexuale digestas. 1753, Holmiae: Impensis Laurentii Salvii [L. Salvius, Stockholm]Google Scholar
- Anderson E: Hybridization of the habitat. Evolution. 1948, 2: 1-9.View ArticleGoogle Scholar
- Shaw KL, Stebbins GL: The role of hybridization in evolution. Proc Am Philos Soc. 1959, 103: 231-251.Google Scholar
- Ehrendorfer FR: Differentiation-hybridization cycles and polyploidy . Cold Spring Harb Symp Quant Biol. 1959, 24: 141-152.View ArticlePubMedGoogle Scholar
- Rattenbury JA: Cyclic hybridization as a survival mechanism in the New Zealand forest flora. Evolution. 1962, 16: 348-363.View ArticleGoogle Scholar
- Levin DA: Hybridization and evolution – a discussion. Taxon. 1970, 19: 167-171.View ArticleGoogle Scholar
- Grant V: Plant speciation. 1981, New York: Columbia University Press, 2Google Scholar
- Arnold ML: Natural hybidization and evolution. 1997, New York: Oxford University PressGoogle Scholar
- Barrier M, Baldwin BG, Robichaux RH, Purugganan MD: Interspecific hybrid ancestry of a plant adaptive radiation: allopolyploidy of the Hawaiian Silversword alliance (Asteraceae) inferred from floral homeotic gene duplications. Mol Biol Evol. 1999, 16 (8): 1105-1113.View ArticlePubMedGoogle Scholar
- Ellstrand NC, Schierenbeck KA: Hybridization as a stimulus for the evolution of invasivness in plants. Proc Natl Acad Sci USA. 2000, 97: 7043-7050.PubMed CentralView ArticlePubMedGoogle Scholar
- Barton NH: The role of hybridization in evolution. Molecular Ecology. 2001, 10: 551-568.View ArticlePubMedGoogle Scholar
- Rieseberg LH, Raymond O, Rosenthal DM, Lai Z, Livingstone K, Nakazato T, Durphy JL, Schwarzbach AE, Donovan LA, Lexer C: Major ecological transitions in annual sunflowers facilitated by hybridization. Science. 2003, 301: 1211-1216.View ArticlePubMedGoogle Scholar
- Seehausen O: Hybridization and adaptive radiation. Trends Ecol Evol. 2004, 19: 198-207.View ArticlePubMedGoogle Scholar
- Wissemann V: Plant evolution by means of hybridization. Systematics and Biodiversity. 2007, 5: 243-253.View ArticleGoogle Scholar
- Mallet J: Hybrid speciation. Nature. 2007, 446: 279-283.View ArticlePubMedGoogle Scholar
- Palmer JD, Shields CR, Cohen DB, Orton TJ: Chloroplast DNA evolution and the origin of amphidiploid Brassica species. Theor Appl Genet. 1983, 65: 181-189.View ArticlePubMedGoogle Scholar
- Sang T, Zhong Y: Testing hybridization hypotheses based on incongruent gene trees. Syst Biol. 2000, 49 (3): 422-434.View ArticlePubMedGoogle Scholar
- Cusimano N, Zhang L-B, Renner SS: Revaluation of the cox1 group I intron in Araceae and Angiosperms suggests a history dominated by loss rather than horizontal transfer. Mol Biol Evol.Google Scholar
- Hallet M, Lagergren J, Tofigh A: Simultaneous identification of duplications and lateral transfers. Proceedings of the eighth annual international conference on Resaerch in computational molecular biology San Diego; ACM. 2004, 347-356.Google Scholar
- Degnan JH, Salter LA: Gene tree distributions under the coalescent process. Evolution. 2005, 59 (1): 24-37.View ArticlePubMedGoogle Scholar
- Doolittle WF, Bapteste E: Pattern pluralism and the Tree of Life hypothesis. Proc Natl Acad Sci U S A. 2007, 104 (7): 2043-2049.PubMed CentralView ArticlePubMedGoogle Scholar
- Huson DH, Dezulian T, Kloepper TK, Steel MA: Phylogenetic supernetworks from partial trees. IEEE Trans Comput Biol Bioinform. 2004, 1: 151-158.View ArticleGoogle Scholar
- Holland BR, Conner G, Huber K, Moulton V: Imputing supertrees and supernetworks from quartets. Syst Biol. 2007, 56 (1): 57-67.View ArticlePubMedGoogle Scholar
- Bininda-Emonds ORP: The evolution of supertrees. Trends Ecol Evol. 2004, 19: 315-22.View ArticlePubMedGoogle Scholar
- Huson D, Steel MA, Whitfield J: Reducing distortion in phylogenetic networks. Proceedings of WABI (Workshop on algorithms in bioinformatics). Edited by: Buecher P, Moret BME. 2006, Springer-Verlage Berlin Heidelberg, Lecture Notes in Bioinformatics 475:150–161Google Scholar
- Huson DH, Bryant D: Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006, 23 (2): 254-267.View ArticlePubMedGoogle Scholar
- Holland BR, Huber KT, Moulton V, Lockhart PJ: Using consensus networks to visualize contradictory evidence for species phylogeny. Mol Biol Evol. 2004, 21 (7): 1459-1461.View ArticlePubMedGoogle Scholar
- Holland BR, Delsuc F, Moulton V: Visualizing conflicting evolutionary hypotheses in large collections of trees using consensus networks. Syst Biol. 2005, 54 (1): 66-76.View ArticlePubMedGoogle Scholar
- Lockhart PJ, McLenachan PA, Havell D, Glenny D, Huson D, Jensen U: Phylogeny, dispersal and radiation of New Zealand alpine buttercups: molecular evidence under split decomposition. Ann Missouri Bot Gard. 2001, 88: 458-477.View ArticleGoogle Scholar
- Edwards SV, Liu L, Pearl DK: High-resolution species trees without concatenation. Proc Natl Acad Sci U S A. 2007, 104 (14): 5936-5941.PubMed CentralView ArticlePubMedGoogle Scholar
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.