Exploring data interaction and nucleotide alignment in a multiple gene analysis of Ips (Coleoptera: Scolytinae)

The possibility of gene tree incongruence in a species-level phylogenetic analysis of the genus Ips (Coleoptera: Scolytidae) was investigated based on mitochondrial 16S rRNA (16S) and nuclear elongation factor-1 alpha (EF-1 alpha) sequences, and existing cytochrome oxidase I (COI) and nonmolecular data sets. Separate cladistic analyses of the data partitions resulted in partially discordant most-parsimonious trees but revealed only low conflict of the phylogenetic signal. Interactions among data partitions, which differed in the extent of sequence divergence (COI > 16S > EF-1 alpha), base composition, and homoplasy, revealed that much of the branch support emerges only in the simultaneous analysis, particularly for deeper nodes in the tree, which are almost entirely supported through "hidden support" (sensu Gatesy et al., Cladistics 15:271-313, 1999). Apparent incongruence between data partitions is in part due to suboptimal alignments and bias of character transformations, but little evidence supports invoking incongruent phylogenetic histories of genetic loci. There is also no justification for eliminating or downweighting gene partitions on the basis of their apparent homoplasy or incongruence with other partitions, because the signal emerges only in the interaction of all data. In comparison with traditional taxonomy, the pini, plastographus, and perturbatus groups are polyphyletic, whereas the grandicollis group is monophyletic except for inclusion of the (monophyletic) calligraphus group. The latidens group and some European species are distantly related and closer to other genera within Ipini. Our robust cladogram was used to revise the classification of Ips. We provide new diagnoses for Ips and four subgeneric taxa.     

Trees within trees: genes and species, molecules and morphology

The construction and interpretation of gene trees is fundamental in molecular systematics. If the gene is defined in a historical (coalescent) sense, there can be multiple gene trees within the single contiguous set of nucleotides, and attempts to construct a single tree for such a sequence must deal with homoplasy created by conflict among divergent histories. On a larger scale, incongruence is expected among gene tree topologies at different loci of individuals within sexually reproducing species, and it has been suggested that this discordance can be used to delimit species. A practical concern for such topological methods is that polymorphisms may be maintained through numerous cladogenic events; this polymorphism problem is less of a concern for nontopological approaches to species delimitation using molecular data. Although a central theoretical concern in molecular systematics is discordance between a given gene tree and the true "species tree," the primary empirical problem faced in reconstructing taxic phylogeny is incongruence among the trees inferred from different sequences. Linkage relationships limit character independence and thus have important implications for handling multiple data sets in phylogenetic analysis, particularly at the species level, where incongruence among different historically associated loci is expected. Gene trees can also be reconstructed for loci that influence phenotypic characters, but there is at best a tenuous relationship between phenotypic homoplasy and homoplasy in such gene trees. Nevertheless, expression patterns and orthology relationships of genes involved in the expression of phenotypes can in theory provide criteria for homology assessment of morphological characters.     

Sources of character conflict in a clade of water striders (Heteroptera: Gerridae)

Incongruence among trees reconstructed with different data may stem from historical (gene tree‐species tree conflict) or process (character change biases) phenomena. Regardless of the source, incongruent data, as determined with “global” measures of homoplasy, have often been excluded from parsimony analysis of the combined data. Recent studies suggest that these homoplasy measures do not predict the contribution of each character to overall tree structure. Branch support measures identify, on a character to node basis, sources of support and conflict resulting from a simultaneous analysis of the data. We implement these branch support measures to identify sources of character conflict in a clade of water striders consisting of Gerris Fabricius, Aquarius Schellenberg, and Limnoporus Stål species. Separate analyses of morphology, mitochondrial cytochrome oxidase I (COI), large mitochondrial ribosomal subunit (16SrRNA), and elongation factor‐1α (EF‐1α) data resulted in cladograms that varied in resolution and topological concordance. Simultaneous analysis of the data resulted in two trees that were unresolved for one node in a strict consensus. The topology agreed with current classification except for the placements of Aquarius chilensis and the Aquarius remigis species group closer to Gerris than to congeneric species. Branch support measures indicated that support derived from each data set varied among nodes, but COI had an overall negative effect on branch support. However, Spearman rank correlation of partitioned branch support values indicated no negative associations of branch support between any data sets and a positive association between EF‐1α and 16SrRNA. Thus incongruence among data sets was not drastic and the gene‐tree versus species tree phenomenon was not implicated. Biases in character change were a more likely reason for incongruence, although saturation curves and incongruence length difference for COI indicated little potential for homoplasy. However, a posteriori inspection of COI nucleotide change with reference to the simultaneous analysis tree revealed AT and codon biases. These biases were not associated with branch support measures. Therefore, it is difficult to predict incongruence or identify its cause. Exclusion of data is ill advised because every character is potentially parsimony informative.     

Mitochondrial DNA,morphology, and the phylogenetic relationships of Antarctic icefishes (Notothenioidei: Channichthyidae)

The Channichthyidae is a lineage of 16 species in the Notothenioidei, a clade of fishes that dominate Antarctic near-shore marine ecosystems with respect to both diversity and biomass. Among four published studies investigating channichthyid phylogeny, no two have produced the same tree topology, and no published study has investigated the degree of phylogenetic incongruence between existing molecular and morphological datasets. In this investigation we present an analysis of channichthyid phylogeny using complete gene sequences from two mitochondrial genes (ND2 and 16S) sampled from all recognized species in the clade. In addition, we have scored all 58 unique morphological characters used in three previous analyses of channichthyid phylogenetic relationships. Data partitions were analyzed separately to assess the amount of phylogenetic resolution provided by each dataset, and phylogenetic incongruence among data partitions was investigated using incongruence length difference (ILD) tests. We utilized a parsimony-based version of the Shimodaira-Hasegawa test to determine if alternative tree topologies are significantly different from trees resulting from maximum parsimony analysis of the combined partition dataset. Our results demonstrate that the greatest phylogenetic resolution is achieved when all molecular and morphological data partitions are combined into a single maximum parsimony analysis. Also, marginal to insignificant incongruence was detected among data partitions using the ILD. Maximum parsimony analysis of all data partitions combined results in a single tree, and is a unique hypothesis of phylogenetic relationships in the Channichthyidae. In particular, this hypothesis resolves the phylogenetic relationships of at least two species (Channichthys rhinoceratus and Chaenocephalus aceratus), for which there was no consensus among the previous phylogenetic hypotheses. The combined data partition dataset provides substantial statistical power to discriminate among alternative hypotheses of channichthyid relationships. These findings suggest the optimal strategy for investigating the phylogenetic relationships of channichthyids is one that uses all available phylogenetic data in analyses of combined data partitions.     

Multiple sources of character information and the phylogeny of Hawaiian drosophilids

Relationships among representatives of the five major Hawaiian Drosophila species groups were examined using data from eight different gene regions. A simultaneous analysis of these data resulted in a single most-parsimonious tree that (1) places the adiastola picture-winged subgroup as sister taxon to the other picture-winged subgroups, (2) unites the modified-tarsus species group with flies from the Antopocerus species group, and (3) places the white-tip scutellum species group as the most basal taxon. Because of the different gene sources used in this study, numerous process partitions can be erected within this data set. We examined the incongruence among these various partitions and the ramifications of these data for the taxonomic consensus, prior agreement, and simultaneous analysis approaches to phylogenetic reconstruction. Separate analyses and taxonomic consensus appear to be inadequate methods for dealing with the partitions in this study. Although detection of incongruence is possible and helps elucidate particular areas of disagreement among data sets, separation of partitions on the basis of incongruence is problematic for many reasons. First, analyzing all genes separately and then either presenting them all as possible hypotheses or taking their consensus provides virtually no information concerning the relationships among these flies. Second, despite some evidence of incongruence, there are no clear delineations among the various gene partitions that separate only heterogeneous data. Third, to the extent that problematic genes can be identified, these genes have nearly the same information content, within a combined analysis framework, as the remaining nonproblematic genes. Our data suggest that significant incongruence among data partitions may be isolated to specific relationships and the "false" signal creating this incongruence is most likely to be overcome by a simultaneous analysis. We present a new method, partitioned Bremer support, for examining the contribution of a particular data partition to the topological support of the simultaneous analysis tree.     

Tempo and mode of evolution in an orthologous Can SINE

Tempo and mode of nucleotide change were examined in an orthologous carnivoran nuclear repetitive DNA element (Can SINE), and compared with those of the transthyretin intron I (TR-i-I) sequence in which it is embedded, by using a phylogenetic framework. The Can SINE is found in representatives of all living caniform carnivoran families, but no living feliform families. This suggests insertion 40–65 MYA, after the two lineages split, but before the caniform radiation. Despite linkage and a long shared evolutionary history, both parsimony and likelihood analyses showed the Can SINE to be significantly different from TR-i-I in rates of evolution and phylogenetic hypotheses supported. The substitution rate is significantly higher in Can SINE than in TR-i-I, and this is attributable to the tRNA-related region of the insertion. While the incongruence length difference test revealed significant conflict between the Can SINE and TR-i-I partitions, the test was shown to be sensitive to the distribution of homoplasy within partitions. The conflicting phylogenies are likely the result of differences in phylogenetic accuracy (homoplasy distribution) rather than in phylogenetic history (gene trees). The base composition of Can SINE contains a significantly higher GC percentage than TR-i-I. Our results indicate that differences between the two partitions may be the result of homoplasy introduced by an increased substitution rate in the tRNA-related region of Can SINE owing to CpG hypermutability. Received: 21 April 2000 / Accepted: 24 August 2000     

Failure of the ILD to determine data combinability for slow loris phylogeny

Tests for incongruence as an indicator of among-data partition conflict have played an important role in conditional data combination. When such tests reveal significant incongruence, this has been interpreted as a rationale for not combining data into a single phylogenetic analysis. In this study of lorisiform phylogeny, we use the incongruence length difference (ILD) test to assess conflict among three independent data sets. A large morphological data set and two unlinked molecular data sets--the mitochondrial cytochrome b gene and the nuclear interphotoreceptor retinoid binding protein (exon 1)--are analyzed with various optimality criteria and weighting mechanisms to determine the phylogenetic relationships among slow lorises (Primates, Loridae). When analyzed separately, the morphological data show impressive statistical support for a monophyletic Loridae. Both molecular data sets resolve the Loridae as paraphyletic, though with different branching orders depending on the optimality criterion or character weighting used. When the three data partitions are analyzed in various combinations, an inverse relationship between congruence and phylogenetic accuracy is observed. Nearly all combined analyses that recover monophyly indicate strong data partition incongruence (P = 0.00005 in the most extreme case), whereas all analyses that recover paraphyly indicate lack of significant incongruence. Numerous lines of evidence verify that monophyly is the accurate phylogenetic result. Therefore, this study contributes to a growing body of information affirming that measures of incongruence should not be used as indicators of data set combinability.     

Understanding phylogenetic incongruence: lessons from phyllostomid bats

All characters and trait systems in an organism share a common evolutionary history that can be estimated using phylogenetic methods. However, differential rates of change and the evolutionary mechanisms driving those rates result in pervasive phylogenetic conflict. These drivers need to be uncovered because mismatches between evolutionary processes and phylogenetic models can lead to high confidence in incorrect hypotheses. Incongruence between phylogenies derived from morphological versus molecular analyses, and between trees based on different subsets of molecular sequences has become pervasive as datasets have expanded rapidly in both characters and species. For more than a decade, evolutionary relationships among members of the New World bat family Phyllostomidae inferred from morphological and molecular data have been in conflict. Here, we develop and apply methods to minimize systematic biases, uncover the biological mechanisms underlying phylogenetic conflict, and outline data requirements for future phylogenomic and morphological data collection. We introduce new morphological data for phyllostomids and outgroups and expand previous molecular analyses to eliminate methodological sources of phylogenetic conflict such as taxonomic sampling, sparse character sampling, or use of different algorithms to estimate the phylogeny. We also evaluate the impact of biological sources of conflict: saturation in morphological changes and molecular substitutions, and other processes that result in incongruent trees, including convergent morphological and molecular evolution. Methodological sources of incongruence play some role in generating phylogenetic conflict, and are relatively easy to eliminate by matching taxa, collecting more characters, and applying the same algorithms to optimize phylogeny. The evolutionary patterns uncovered are consistent with multiple biological sources of conflict, including saturation in morphological and molecular changes, adaptive morphological convergence among nectar‐feeding lineages, and incongruent gene trees. Applying methods to account for nucleotide sequence saturation reduces, but does not completely eliminate, phylogenetic conflict. We ruled out paralogy, lateral gene transfer, and poor taxon sampling and outgroup choices among the processes leading to incongruent gene trees in phyllostomid bats. Uncovering and countering the possible effects of introgression and lineage sorting of ancestral polymorphism on gene trees will require great leaps in genomic and allelic sequencing in this species‐rich mammalian family. We also found evidence for adaptive molecular evolution leading to convergence in mitochondrial proteins among nectar‐feeding lineages. In conclusion, the biological processes that generate phylogenetic conflict are ubiquitous, and overcoming incongruence requires better models and more data than have been collected even in well‐studied organisms such as phyllostomid bats.     

Molecular phylogenetic analysis of tropical freshwater mussels (Mollusca: Bivalvia: Unionoida) resolves the position of Coelatura and supports a monophyletic Unionidae

In previous molecular phylogenetic analyses of the freshwater mussel family Unionidae (Bivalvia: Unionoida), the Afrotropical genus Coelatura had been recovered in various positions, generally indicating a paraphyletic Unionidae. However that result was typically poorly supported and in conflict with morphology-based analyses. We set out to test the phylogenetic position of Coelatura by sampling tropical lineages omitted from previous studies. Forty-one partial 28S nuclear rDNA and partial COI mtDNA sequences (1130 total aligned nucleotides) were analyzed separately and in combination under both maximum parsimony and likelihood, as well as Bayesian inference. There was significant phylogenetic incongruence between the character sets (partition homogeneity test, p < 0.01), but a novel heuristic for comparing bootstrap values among character sets analyzed separately and in combination illustrated that the observed conflict was due to homoplasy rather than separate gene histories. Phylogenetic analyses robustly supported a monophyletic Unionidae, with Coelatura recovered as part of a well-supported Africa–India clade (= Parreysiinae). The implications of this result are discussed in the context of Afrotropical freshwater mussel evolution and the classification of the family Unionidae.     

Exploring gene tree incongruence at the origin of ants and bees (Hymenoptera)

The fact that different phylogenomic data sets can lead to highly supported but inconsistent results suggest that conflict among gene trees in real data sets could be severe. We provide here a detailed exploration of gene tree space to investigate the relationships in Hymenoptera based on data obtained by Johnson et al. (Current Biology, 2013, 23, 2058), in which ants and Apoidea (bees and spheciform wasps) were recovered as sister groups, contradicting previous studies. We found high levels of topological variation among gene trees, several of them disagreeing with previously published hypotheses. To profile the dynamics of emerging support versus conflicting signal in combined analysis of data, we employed a novel method based on the incremental addition of randomized data to coalescence‐based phylogenetic inference. Although the monophyly of Aculeata and of Formicidae were consistently recovered using as little as 6.5% of the 308 available markers, signal for the Formicidae + Apoidea clade prevailed only after more than 50% of the loci were sampled. Still, non‐negligible support for alternative hypotheses remained until all genes were added to the analysis. Our results suggest that phylogenetic conflict is rather pervasive and not scattered as noise across individual gene trees because alternative topologies were recovered not from a specific subset, but from several random combinations of loci. Thus, even though phylogenetic signal recovered from full gene data sets was already dominant in much smaller ensembles, large amounts of data may be indeed necessary to overcome phylogenetic conflict.     

Phylogeny of the Liliales (Monocotyledons) with special emphasis on data partition congruence and RNA editing

A phylogenetic analysis of the monocot order Liliales was performed using sequence data from three mitochondrial (atp1, cob, nad5) and two plastid genes (rbcL, ndhF). The complete data matrix includes 46 terminals representing all 10 families currently included in Liliales. The two major partitions, mitochondrial and plastid data, were congruent, and parsimony analysis resulted in 50 equally parsimonious trees and a well resolved consensus tree confirming monophyly of all families. Mitochondrial genes are known to include RNA edited sites, and in some cases unprocessed genes are replaced by retro‐processed gene copies, that is processed paralogs. To test the effects on phylogeny reconstruction of predicted edited sites and potentially unintentionally sampled processed paralogs, a number of analyses were performed using subsets of the complete data matrix. In general, predicted edited sites were more homoplasious than the other characters and increased incongruence among most data partitions. The predicted edited sites have a non‐random phylogenetic signal in conflict with the signal of the non‐edited sites. The potentially misleading signal was caused partially by the apparent presence of processed paralogs in Galanthus (Amaryllidaceae), part of the outgroup, but also by a deviating evolutionary pattern of predicted edited sites in Liliaceae compared with the remainder of the Liliales. Despite the problems that processed paralogs may cause, we argue that they should not a priori be excluded from phylogenetic analysis.     

Phylogenetic analysis of Themira (Sepsidae: Diptera): sensitivity analysis, alignment, and indel treatment in a multigene study

In this study we use sensitivity analysis sensu Wheeler (1995 ) for a matrix entirely composed of DNA sequences. We propose that not only congruence but also phylogenetic structure, as measured by character resampling, should be used to choose among competing weighting regimes. An extensive analysis of a five‐gene data set for Themira (Sepsidae: Diptera) reveals that even with different ways of partitioning the data, measures of topological congruence, character incongruence, and phylogenetic structure favor similar weighting regimes involving the down‐weighting of transitions. We furthermore use sensitivity analysis for obtaining empirical evidence that allows us to select weights for third positions, deciding between treating indels as fifth character states or missing values, and choosing between manual and computational alignments. For our data, sensitivity analysis favors manual alignment over a Clustal‐generated numerical alignment, the treatment of indels as fifth character states over considering them missing values, and equal weights for all positions in protein‐encoding genes over the down‐weighting of third positions. Among the topological congruence measures compared, symmetric tree distance performed best. Partitioned Bremer Support analysis reveals that COI contributes the largest amount of support for our phylogenetic tree for Themira. © The Willi Hennig Society 2005.     

Phylogeny of bumble bees in the New World subgenus Fervidobombus (Hymenoptera: Apidae): congruence of molecular and morphological data

We present new DNA sequence data (12S, 16S, and opsin gene fragments) and morphological characters of the male genitalia for a phylogenetic analysis of the bumble bee subgenus Fervidobombus. There is no significant incongruence between the three molecular data sets, and little incongruence between the DNA and morphology. Simultaneous analysis of all the data partitions resulted in a tree that was entirely congruent with the All-DNA tree. Optimization of the geographic locations of the taxa onto the tree topology using dispersal/vicariance analysis suggests a complex picture of spread and diversification of Fervidobombus from the Old World into the southern New World. There is a phylogenetic component to their spread into tropical rain forest, as the two primary rain forest species (Bombus transversalis and Bombus pullatus) comprise a monophyletic clade, along with a third species, Bombus atratus, which is widely distributed in South America, including lowland subtropical habitats.     

Weighting and Congruence: A Case Study Based on Three Mitochondrial Genes in Pitvipers

Weighting strategies in a total evidence approach are often conducted in order to remove homoplasy, with the implicit hope to increase congruence between data partitions. Incongruence was measured using the test of Farriset al.(J. S. Farris, M. Källersjö, A. G. Kluge, and C. Bult, 1995,Cladistics10, 315–319) within and between three mitochondrial genes from pitvipers (Crotalinae) in partitioning each codon position for the coding genes. Incongruence between codon positions of a gene can be worse than incongruence between genes at analogous positions. Third positions of codons behave very differently in terms of incongruence from one gene to another while showing similar patterns in saturation tests. Instead of removing characters in order to discard homoplasy, which is hopeless and does not increase general congruence, we advocate for the removal of those substitutions that are incongruent with the rest. The genusCalloselasmaand its sister group the genusHypnaleare the most basal Crotalinae. Asiatic pitvipers are paraphyletic, while American pitvipers are monophyletic.     

The phylogenetic position of the Pterosauria within the Archosauromorpha re-examined

A previous analysis of the phylogenetic position of the Pterosauria argued that pterosaurs were not closely related to dinosaurs as is generally accepted, but rather were outside the crown group Archosauria. However, that study was dismissed for the use of inappropriate methods. Here, the data set from that analysis was divided into five partitions: one with characters associated with cursorial digitigrade bipedal locomotion and the other four with characters from the skull and mandible, postcranial axial skeleton, forelimb and hindlimb, respectively. The partitions were subjected to homogeneity testing, and the Cursorial partition was found to be incongruent with other partitions and all other characters at the α = 0.01 probability level. Deletion of the Pterosauria removed all significant incongruence, demonstrating that the incongruence results from the coding of pterosaurs for the cursorial characters. The cause of the incongruence was interpreted as homoplasy in hindlimb morphology, and after re-evaluating and reformulating the characters of the Cursorial partition, the revised data set was tested for homogeneity and no significant incongruence was found. Lastly, the data set was updated with additional characters and taxa from recent analyses, tested as before, and when analysed suggested that the Pterosauria were basal archosauriforms well outside the crown group Archosauria.     

Phylogenetic relationships of typical antbirds (Thamnophilidae) and test of incongruence based on Bayes factors

Background  

The typical antbirds (Thamnophilidae) form a monophyletic and diverse family of suboscine passerines that inhabit neotropical forests. However, the phylogenetic relationships within this assemblage are poorly understood. Herein, we present a hypothesis of the generic relationships of this group based on Bayesian inference analyses of two nuclear introns and the mitochondrial cytochrome b gene. The level of phylogenetic congruence between the individual genes has been investigated utilizing Bayes factors. We also explore how changes in the substitution models affected the observed incongruence between partitions of our data set.     

Attributes and congruence of three molecular data sets: inferring phylogenies among Septoria-related species from woody perennial plants

To improve our understanding of phylogenetic relationships within the anamorphic genus Septoria, three molecular data sets representing 2,417 bp of nuclear and mitochondrial genes were evaluated. Separate gene analyses and combined analyses were performed using first, the maximum parsimony criterion and second, a Bayesian framework. The homogeneity of data partitions was evaluated via a combination of homogeneity partition tests and tree topology incongruence tests before conducting combined analyses. A last incongruence re-evaluation using partitioned Bremer support was performed on the combined tree, which corroborated the previous estimates. After each separate data set attributes were examined, simple explanations were advocated as the causes of the significant incongruences detected. The analysis of multiple gene partitions showed unprecedented phylogenetic resolution within the genus Septoria that supported the results from previously published single gene phylogenies. Specifically, we have delimited distinct but closely related species representing monophyletic groups that frequently correlated with their respective host families. Conversely, the occurrence of well-supported groups including closely related but distinct molecular taxa sampled on unrelated host-plants allowed us to reject, in these particular cases, the co-evolutionary concept expected between a parasite and its host and to discuss alternative evolutionary models recently proposed for these pathogens.     

iGTP: A software package for large-scale gene tree parsimony analysis

Background  

The ever-increasing wealth of genomic sequence information provides an unprecedented opportunity for large-scale phylogenetic analysis. However, species phylogeny inference is obfuscated by incongruence among gene trees due to evolutionary events such as gene duplication and loss, incomplete lineage sorting (deep coalescence), and horizontal gene transfer. Gene tree parsimony (GTP) addresses this issue by seeking a species tree that requires the minimum number of evolutionary events to reconcile a given set of incongruent gene trees. Despite its promise, the use of gene tree parsimony has been limited by the fact that existing software is either not fast enough to tackle large data sets or is restricted in the range of evolutionary events it can handle.     

A phylogenomic perspective on gene tree conflict and character evolution in Caprifoliaceae using target enrichment data,with Zabelioideae recognized as a new subfamily

The use of diverse data sets in phylogenetic studies aiming for understanding evolutionary histories of species can yield conflicting inference. Phylogenetic conflicts observed in animal and plant systems have often been explained by hybridization, incomplete lineage sorting (ILS), or horizontal gene transfer. Here, we used target enrichment data, species tree, and species network approaches to infer the backbone phylogeny of the family Caprifoliaceae, while distinguishing among sources of incongruence. We used 713 nuclear loci and 46 complete plastome sequence data from 43 samples representing 38 species from all major clades to reconstruct the phylogeny of the family using concatenation and coalescence approaches. We found significant nuclear gene tree conflict as well as cytonuclear discordance. Additionally, coalescent simulations and phylogenetic species network analyses suggested putative ancient hybridization among subfamilies of Caprifoliaceae, which seems to be the main source of phylogenetic discordance. Ancestral state reconstruction of six morphological characters revealed some homoplasy for each character examined. By dating the branching events, we inferred the origin of Caprifoliaceae at approximately 66.65 Ma in the late Cretaceous. By integrating evidence from molecular phylogeny, divergence times, and morphology, we here recognize Zabelioideae as a new subfamily in Caprifoliaceae. This work shows the necessity of using a combination of multiple approaches to identify the sources of gene tree discordance. Our study also highlights the importance of using data from both nuclear and plastid genomes to reconstruct deep and shallow phylogenies of plants.