Coalescent simulations reveal hybridization and incomplete lineage sorting in Mediterranean Linaria

We examined the phylogenetic history of Linaria with special emphasis on the Mediterranean sect. Supinae (44 species). We revealed extensive highly supported incongruence among two nuclear (ITS, AGT1) and two plastid regions (rpl32-trnL(UAG), trnS-trnG). Coalescent simulations, a hybrid detection test and species tree inference in *BEAST revealed that incomplete lineage sorting and hybridization may both be responsible for the incongruent pattern observed. Additionally, we present a multilabelled *BEAST species tree as an alternative approach that allows the possibility of observing multiple placements in the species tree for the same taxa. That permitted the incorporation of processes such as hybridization within the tree while not violating the assumptions of the *BEAST model. This methodology is presented as a functional tool to disclose the evolutionary history of species complexes that have experienced both hybridization and incomplete lineage sorting. The drastic climatic events that have occurred in the Mediterranean since the late Miocene, including the Quaternary-type climatic oscillations, may have made both processes highly recurrent in the Mediterranean flora.     

Patterns and causes of incongruence between plastid and nuclear Senecioneae (Asteraceae) phylogenies

One of the longstanding questions in phylogenetic systematics is how to address incongruence among phylogenies obtained from multiple markers and how to determine the causes. This study presents a detailed analysis of incongruent patterns between plastid and ITS/ETS phylogenies of Tribe Senecioneae (Asteraceae). This approach revealed widespread and strongly supported incongruence, which complicates conclusions about evolutionary relationships at all taxonomic levels. The patterns of incongruence that were resolved suggest that incomplete lineage sorting (ILS) and/or ancient hybridization are the most likely explanations. These phenomena are, however, extremely difficult to distinguish because they may result in similar phylogenetic patterns. We present a novel approach to evaluate whether ILS can be excluded as an explanation for incongruent patterns. This coalescence-based method uses molecular dating estimates of the duration of the putative ILS events to determine if invoking ILS as an explanation for incongruence would require unrealistically high effective population sizes. For four of the incongruent patterns identified within the Senecioneae, this approach indicates that ILS cannot be invoked to explain the observed incongruence. Alternatively, these patterns are more realistically explained by ancient hybridization events.     

Coalescent histories on phylogenetic networks and detection of hybridization despite incomplete lineage sorting

Analyses of the increasingly available genomic data continue to reveal the extent of hybridization and its role in the evolutionary diversification of various groups of species. We show, through extensive coalescent-based simulations of multilocus data sets on phylogenetic networks, how divergence times before and after hybridization events can result in incomplete lineage sorting with gene tree incongruence signatures identical to those exhibited by hybridization. Evolutionary analysis of such data under the assumption of a species tree model can miss all hybridization events, whereas analysis under the assumption of a species network model would grossly overestimate hybridization events. These issues necessitate a paradigm shift in evolutionary analysis under these scenarios, from a model that assumes a priori a single source of gene tree incongruence to one that integrates multiple sources in a unifying framework. We propose a framework of coalescence within the branches of a phylogenetic network and show how this framework can be used to detect hybridization despite incomplete lineage sorting. We apply the model to simulated data and show that the signature of hybridization can be revealed as long as the interval between the divergence times of the species involved in hybridization is not too small. We reanalyze a data set of 106 loci from 7 in-group Saccharomyces species for which a species tree with no hybridization has been reported in the literature. Our analysis supports the hypothesis that hybridization occurred during the evolution of this group, explaining a large amount of the incongruence in the data. Our findings show that an integrative approach to gene tree incongruence and its reconciliation is needed. Our framework will help in systematically analyzing genomic data for the occurrence of hybridization and elucidating its evolutionary role.     

Testing hybridization hypotheses based on incongruent gene trees

Hybridization is an important evolutionary mechanism in plants and has been increasingly documented in animals. Difficulty in reconstruction of reticulate evolution, however, has been a long-standing problem in phylogenetics. Consequently, hybrid speciation may play a major role in causing topological incongruence between gene trees. The incongruence, in turn, offers an opportunity to detect hybrid speciation. Here we characterized certain distinctions between hybridization and other biological processes, including lineage sorting, paralogy, and lateral gene transfer, that are responsible for topological incongruence between gene trees. Consider two incongruent gene trees with three taxa, A, B, and C, where B is a sister group of A on gene tree 1 but a sister group of C on gene tree 2. With a theoretical model based on the molecular clock, we demonstrate that time of divergence of each gene between taxa A and C is nearly equal in the case of hybridization (B is a hybrid) or lateral gene transfer, but differs significantly in the case of lineage sorting or paralogy. After developing a bootstrap test to test these alternative hypotheses, we extended the model and test to account for incongruent gene trees with numerous taxa. Computer simulation studies supported the validity of the theoretical model and bootstrap test when each gene evolved at a constant rate. The computer simulation also suggested that the model remained valid as long as the rate heterogeneity was occurring proportionally in the same taxa for both genes. Although the model could not test hypotheses of hybridization versus lateral gene transfer as the cause of incongruence, these two processes may be distinguished by comparing phylogenies of multiple unlinked genes.     

Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: A model

The application of phylogenetic inference methods, to data for a set of independent genes sampled randomly throughout the genome, often results in substantial incongruence in the single-gene phylogenetic estimates. Among the processes known to produce discord between single-gene phylogenies, two of the best studied in a phylogenetic context are hybridization and incomplete lineage sorting. Much recent attention has focused on the development of methods for estimating species phylogenies in the presence of incomplete lineage sorting, but phylogenetic models that allow for hybridization have been more limited. Here we propose a model that allows incongruence in single-gene phylogenies to be due to both hybridization and incomplete lineage sorting, with the goal of determining the contribution of hybridization to observed gene tree incongruence in the presence of incomplete lineage sorting. Using our model, we propose methods for estimating the extent of the role of hybridization in both a likelihood and a Bayesian framework. The performance of our methods is examined using both simulated and empirical data.     

Multi‐locus phylogenetics,lineage sorting,and reticulation in Pinus subsection Australes

Premise of the Study

Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long‐lived, wind‐pollinated plants with large population sizes and weak reproductive barriers.

Methods

We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low‐copy‐number nuclear genes and high‐copy‐number plastomes (Hyb‐Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation.

Key Results

Concatenation‐ and coalescent‐based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed‐cone pines) and Oocarpae (the egg‐cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non‐monophyly of the Attenuatae in the plastid tree was equivocal.

Conclusions

Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.     

Estimating the species tree for Hawaiian Schiedea (Caryophyllaceae) from multiple loci in the presence of reticulate evolution

Schiedea (Caryophyllaceae) is a monophyletic genus of 34 species, all endemic to the Hawaiian Islands, that arose from a single colonization, providing one of the best examples of adaptive radiation in Hawai'i. Species utilize a range of habitats and exhibit a variety of growth forms and transitions in breeding systems from hermaphroditism toward dimorphism or autogamy. Our study included the most thorough sampling to date: 2-5 individuals per species and 4 independent genetic partitions: eight plastid and three low-copy nuclear loci (9217bps), allowing a three-locus BEST species tree. Despite incomplete resolution at the tips, our results support monophyly for each extant species. Gene trees revealed several clear cases of cytonuclear incongruence, likely created by interspecific introgression. Conflict occurs at the divergence of section Alphaschiedea as well as at the tips. Ages inferred from a BEAST analysis allow an original colonization onto either Nihoa or Kauaì and inform some aspects of inter-island migrations. We suggest that several hard polytomies on the species tree are biologically realistic, signifying either nearly simultaneous speciation or historical introgressive hybridization. Based on inferred node ages that exceed expected coalescent times, we propose that undetected nuclear introgression may play a larger role than incomplete lineage sorting in sections Schiedea and Mononeura.     

Complex evolution in Arundinarieae (Poaceae: Bambusoideae): incongruence between plastid and nuclear GBSSI gene phylogenies

The monophyly of tribe Arundinarieae (the temperate woody bamboos) has been unequivocally recovered in previous molecular phylogenetic studies. In a recent phylogenetic study, 10 major lineages in Arundinarieae were resolved based on eight non-coding plastid regions, which conflicted significantly with morphological classifications both at the subtribal and generic levels. Nevertheless, relationships among and within the 10 lineages remain unclear. In order to further unravel the evolutionary history of Arundinarieae, we used the nuclear GBSSI gene sequences along with those of eight plastid regions for phylogenetic reconstruction, with an emphasis on Chinese species. The results of the plastid analyses agreed with previous studies, whereas 13 primary clades revealed in the GBSSI phylogeny were better resolved at the generic level than the plastid phylogeny. Our analyses also revealed many inconsistencies between the plastid DNA and the nuclear GBSSI trees. These results implied that the nuclear genome and the plastid genome had different evolutionary trajectories. The patterns of incongruence suggested that lack of informative characters, incomplete lineage sorting, and/or hybridization (introgression) could be the causes. Seven putative hybrid species were hypothesized, four of which are discussed in detail on the basis of topological incongruence, chromosome numbers, morphology, and distribution patterns, and those taxa probably resulted from homoploid hybrid speciation. Overall, our study indicates that the tribe Arundinarieae has undergone a complex evolution.     

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.     

Phylogeny of Nolana (Solanaceae) of the Atacama and Peruvian deserts inferred from sequences of four plastid markers and the nuclear LEAFY second intron

The phylogeny of Nolana (Solanaceae), a genus primarily distributed in the coastal Atacama and Peruvian deserts with a few species in the Andes and one species endemic to the Galápagos Islands, was reconstructed using sequences of four plastid regions (ndhF, psbA-trnH, rps16-trnK and trnC-psbM) and the nuclear LEAFY second intron. The monophyly of Nolana was strongly supported by all molecular data. The LEAFY data suggested that the Chilean species, including Nolana sessiliflora, the N. acuminata group and at least some members of the Alona group, are basally diverged, supporting the Chilean origin of the genus. Three well-supported clades in the LEAFY tree were corroborated by the SINE (short interspersed elements) or SINE-like insertions. Taxa from Peru are grouped roughly into two clades. Nolana galapagensis from the Galápagos Island is most likely to have derived from a Peruvian ancestor. The monophyly of the morphologically well-diagnosed Nolana acuminata group (N. acuminata, N. baccata, N. paradoxa, N. parviflora, N. pterocarpa, N. rupicola and N. elegans) was supported by both plastid and LEAFY data. Incongruence between the plastid and the LEAFY data was detected concerning primarily the positions of N. sessiliflora, N. galapagensis, taxa of the Alona group and the two Peruvian clades. Such incongruence may be due to reticulate evolution or in some cases lineage sorting of plastid DNA. Incongruence between our previous GBSSI trees and the plastid-LEAFY trees was also detected concerning two well-supported major clades in the GBSSI tree. Duplication of the GBSSI gene may have contributed to this incongruence.     

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.     

The reticulate history of Medicago (Fabaceae)

The phylogenetic history of Medicago was examined for 60 accessions from 56 species using two nuclear genes (CNGC5 and beta-cop) and one mitochondrial region (rpS14-cob). The results of several analyses revealed that extensive robustly supported incongruence exists among the nuclear genes, the cause of which we seek to explain. After rejecting several processes, hybridization and lineage sorting of ancestral polymorphisms remained as the most likely factors promoting incongruence. Using coalescence simulations, we rejected lineage sorting alone as an explanation of the differences among gene trees. The results indicate that hybridization has been common and ongoing among lineages since the origin of Medicago. Coalescence provides a good framework to test the causes of incongruence commonly seen among gene trees but requires knowledge of effective population sizes and generation times. We estimated the effective population size at 240,000 individuals and assumed a generation time of 1 year in Medicago (many are annual plants). A sensitivity analysis showed that our conclusions remain unchanged using a larger effective population size and/or longer generation time.     

Widespread discordance of gene trees with species tree in Drosophila: evidence for incomplete lineage sorting

The phylogenetic relationship of the now fully sequenced species Drosophila erecta and D. yakuba with respect to the D. melanogaster species complex has been a subject of controversy. All three possible groupings of the species have been reported in the past, though recent multi-gene studies suggest that D. erecta and D. yakuba are sister species. Using the whole genomes of each of these species as well as the four other fully sequenced species in the subgenus Sophophora, we set out to investigate the placement of D. erecta and D. yakuba in the D. melanogaster species group and to understand the cause of the past incongruence. Though we find that the phylogeny grouping D. erecta and D. yakuba together is the best supported, we also find widespread incongruence in nucleotide and amino acid substitutions, insertions and deletions, and gene trees. The time inferred to span the two key speciation events is short enough that under the coalescent model, the incongruence could be the result of incomplete lineage sorting. Consistent with the lineage-sorting hypothesis, substitutions supporting the same tree were spatially clustered. Support for the different trees was found to be linked to recombination such that adjacent genes support the same tree most often in regions of low recombination and substitutions supporting the same tree are most enriched roughly on the same scale as linkage disequilibrium, also consistent with lineage sorting. The incongruence was found to be statistically significant and robust to model and species choice. No systematic biases were found. We conclude that phylogenetic incongruence in the D. melanogaster species complex is the result, at least in part, of incomplete lineage sorting. Incomplete lineage sorting will likely cause phylogenetic incongruence in many comparative genomics datasets. Methods to infer the correct species tree, the history of every base in the genome, and comparative methods that control for and/or utilize this information will be valuable advancements for the field of comparative genomics.     

Multilocus phylogenetics of a rapid radiation in the genus Thomomys (Rodentia: Geomyidae)

Species complexes undergoing rapid radiation present a challenge in molecular systematics because of the possibility that ancestral polymorphism is retained in component gene trees. Coalescent theory has demonstrated that gene trees often fail to match lineage trees when taxon divergence times are less than the ancestral effective population sizes. Suggestions to increase the number of loci and the number of individuals per taxon have been proposed; however, phylogenetic methods to adequately analyze these data in a coalescent framework are scarce. We compare two approaches to estimating lineage (species) trees using multiple individuals and multiple loci: the commonly used partitioned Bayesian analysis of concatenated sequences and a modification of a newly developed hierarchical Bayesian method (BEST) that simultaneously estimates gene trees and species trees from multilocus data. We test these approaches on a phylogeny of rapidly radiating species wherein divergence times are likely to be smaller than effective population sizes, and incomplete lineage sorting is known, in the rodent genus, Thomomys. We use seven independent noncoding nuclear sequence loci (total approximately 4300 bp) and between 1 and 12 individuals per taxon to construct a phylogenetic hypothesis for eight Thomomys species. The majority-rule consensus tree from the partitioned concatenated analysis included 14 strongly supported bipartitions, corroborating monophyletic species status of five of the eight named species. The BEST tree strongly supported only the split between the two subgenera and showed very low support for any other clade. Comparison of both lineage trees to individual gene trees revealed that the concatenation method appears to ignore conflicting signals among gene trees, whereas the BEST tree considers conflicting signals and downweights support for those nodes. Bayes factor analysis of posterior tree distributions from both analyses strongly favor the model underlying the BEST analysis. This comparison underscores the risks of overreliance on results from concatenation, and ignoring the properties of coalescence, especially in cases of recent, rapid radiations.     

From gene trees to species trees II: species tree inference by minimizing deep coalescence events

When gene copies are sampled from various species, the resulting gene tree might disagree with the containing species tree. The primary causes of gene tree and species tree discord include incomplete lineage sorting, horizontal gene transfer, and gene duplication and loss. Each of these events yields a different parsimony criterion for inferring the (containing) species tree from gene trees. With incomplete lineage sorting, species tree inference is to find the tree minimizing extra gene lineages that had to coexist along species lineages; with gene duplication, it becomes to find the tree minimizing gene duplications and/or losses. In this paper, we present the following results: 1) The deep coalescence cost is equal to the number of gene losses minus two times the gene duplication cost in the reconciliation of a uniquely leaf labeled gene tree and a species tree. The deep coalescence cost can be computed in linear time for any arbitrary gene tree and species tree. 2) The deep coalescence cost is always not less than the gene duplication cost in the reconciliation of an arbitrary gene tree and a species tree. 3) Species tree inference by minimizing deep coalescence events is NP-hard.     

Palaeophylogenomics of the vertebrate ancestor--impact of hidden paralogy on hagfish and lamprey gene phylogeny

In dissecting the transition from invertebrates to vertebrates at the molecular level, whole-genome duplications are recognized as a key event. This gave rise to more copies of genes in jawed vertebrates (gnathostomes), such as the four Hox clusters in the human, compared to the single ancestral cluster in invertebrates. To date, as the most early-branching lineages in vertebrates, cyclostomes (hagfishes and lampreys) have been used for comparative analyses of gene regulations and functions. However, assignment of orthology/paralogy for cyclostomes' genes is not unambiguously demonstrated. Thus, there is a high degree of incongruence in tree topologies between gene families, although whole genome duplications postulate uniform patterns in gene phylogeny. In this review, we demonstrate how expansion of an ancient genome before the cyclostome-gnathostome split, followed by reciprocal gene loss, can cause this incongruence. This is sometimes referred to as 'hidden paralogy'.     

The probability of a gene tree topology within a phylogenetic network with applications to hybridization detection

Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa.     

Introgressive Hybridization between Anciently Diverged Lineages of Silene (Caryophyllaceae)

Hybridization has played a major role during the evolution of angiosperms, mediating both gene flow between already distinct species and the formation of new species. Newly formed hybrids between distantly related taxa are often sterile. For this reason, interspecific crosses resulting in fertile hybrids have rarely been described to take place after more than a few million years after divergence. We describe here the traces of a reproductively successful hybrid between two ancestral species of Silene, diverged for about six million years prior to hybridization. No extant hybrids between the two parental lineages are currently known, but introgression of the RNA polymerase gene NRPA2 provides clear evidence of a temporary and fertile hybrid. Parsimony reconciliation between gene trees and the species tree, as well as consideration of clade ages, help exclude gene paralogy and lineage sorting as alternative hypotheses. This may represent one of the most extreme cases of divergence between species prior to introgressive hybridization discovered yet, notably at a homoploid level. Although species boundaries are generally believed to be stable after millions of years of divergence, we believe that this finding may indicate that gene flow between distantly related species is merely largely undetected at present.     

Random roots and lineage sorting

Lineage sorting has been suggested as a major force in generating incongruent phylogenetic signal when multiple gene partitions are examined. The degree of lineage sorting can be estimated using the coalescent process and simulation studies have also pointed to a major role for incomplete lineage sorting as a factor in phylogenetic inference. Some recent empirical studies point to an extreme role for this phenomenon with up to 50-60% of all informative genes showing incongruence as a result of lineage sorting. Here, we examine seven large multi-partition genome level data sets over a large range of taxonomic representation. We took the approach of examining outgroup choice and its impact on tree topology, by swapping outgroups into analyses with successively larger genetics distances to the ingroup. Our results indicate a linear relationship of outgroup distance with incongruence in the data sets we examined suggesting a strong random rooting effect. In addition, we attempted to estimate the degree of lineage sorting in several large genome level data sets by examining triads of very closely related taxa. This exercise resulted in much lower estimates of incongruent genes that could be the result of lineage sorting, with an overall estimate of around 10% of the total number of genes in a genome showing incongruence as a result of true lineage sorting. Finally we examined the behavior of likelihood and parsimony approaches on the random rooting phenomenon. Likelihood tends to stabilize incongruence as outgroups get further and further away from the ingroup. In one extreme case, likelihood overcompensates for sequence divergence but increases random rooting causing long branch repulsion.     

Perils of paralogy: using HSP70 genes for inferring organismal phylogenies

Conserved genes have found their way into the mainstream of molecular systematics. Many of these genes are members of multigene families. A difficulty with using single genes of multigene families for phylogenetic inference is that genes from one species may be paralogous to those from another taxon. We focus attention on this problem using heat shock 70 (HSP70) genes. Using polymerase chain reaction techniques with genomic DNA, we isolated and sequenced 123 distinct sequences from 12 species of sharks. Phylogenetic analysis indicated that the sequences cluster with constituitively expressed cytoplasmic heat shock-like genes. Three highly divergent gene clades were sampled. A number of similar sequences were sampled from each species within each distinct gene clade. Comparison of published species trees with an HSP70 gene tree inferred using Bayesian phylogenetic analysis revealed several cases of gene duplication and differential sorting of gene lineages within this group of sharks. Gene tree parsimony based on the objective criteria of duplication and losses showed that previously published hypotheses of species relationships and two novel hypothesis based on Bayesian phylogenetics were concordant with the history of HSP70 gene duplication and loss. By contrast, two published hypotheses based on morphological data were not significantly different from the null hypothesis of a random association between species relatedness and the HSP70 gene tree. These results suggest that gene tree parsimony using data from multigene families can be used for inferring species relationships or testing published alternative hypotheses. More importantly, the results suggest that systematic studies relying on phylogenetic inferences from HSP70 genes may by plagued by unrecognized paralogy of sampled genes. Our results underscore the distinction between gene and species trees and highlight an underappreciated source of discordance between gene trees and organismal phylogeny, i.e., unrecognized paralogy of sampled genes.