Hybridization and genome duplication for early evolutionary success in the Asian Palmate group of Araliaceae

The phenomenal advances in sequencing techniques and analytical development during the last decade have provided a unique opportunity to unravel the evolutionary history of lineages under complex patterns of evolution. This is the case of the largest clade of the ginseng family (Araliaceae), the Asian Palmate group (AsPG), where the large internal polytomies and genome incongruences detected in previous studies pointed to a scenario of radiation with hybridization events between genera for the early evolution of the group. In this study, we aim to obtain well-resolved nuclear and plastid phylogenies of the AsPG using Hyb-Seq to evaluate the radiation hypothesis and assess the role of hybridization in the early evolution of the group. We performed concatenated- and coalescent-based phylogenetic analyses from the 936 targeted nuclear loci and 261 plastid loci obtained for 72 species representing 20 genera of the AsPG and the main clades of Araliaceae. The impact of hybridization and incomplete lineage sorting (ILS) was assessed with SNaQ, and genome duplications were evaluated with ChromEvol. Our nuclear and plastid phylogenies are compatible with a scenario of early radiation in the AsPG. Also, the identification of extensive signals of hybridization and ILS behind the genome incongruences supports hybridization as a major driving force during the early radiation. We hypothesize a whole-genome duplication event at the origin of the AsPG, followed by a radiation that led to extensive ILS, which, alongside the early inter-genera hybridization, is obscuring the phylogenetic signal in the early evolution of this major clade.     

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.     

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.     

Can three incongruence tests predict when data should be combined?

Advocates of conditional combination have argued that testing for incongruence between data partitions is an important step in data exploration. Unless the partitions have had distinct histories, as in horizontal gene transfer, incongruence means that one or more data support the wrong phylogeny. This study examines the relationship between incongruence and phylogenetic accuracy using three tests of incongruence. These tests were applied to pairs of mitochondrial DNA data partitions from two well-corroborated vertebrate phylogenies. Of the three tests, the most useful was the incongruence length difference test (ILD, also called the partition homogeneity test). This test distinguished between cases in which combining the data generally improved phylogenetic accuracy (P > 0.01) and cases in which accuracy of the combined data suffered relative to the individual partitions (P < 0.001). In contrast, in several cases, the Templeton and Rodrigo tests detected highly significant incongruence (P < 0.001) even though combining the incongruent partitions actually increased phylogenetic accuracy. All three tests identified cases in which improving the reconstruction model would improve the phylogenetic accuracy of the individual partitions.      

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.     

A biogeographic–ecological approach to disentangle reticulate evolution in the Triatoma phyllosoma species group (Heteroptera: Triatominae), vectors of Chagas disease

Introgression and incomplete lineage sorting (ILS) are two of the main sources of gene‐tree incongruence; both can confound the assessment of phylogenetic relationships among closely related species. The Triatoma phyllosoma species group is a clade of partially co‐distributed and cross‐fertile Chagas disease vectors. Despite previous efforts, the phylogeny of this group remains unresolved, largely because of substantial gene‐tree incongruence. Here, we sequentially address introgression and ILS to provide a robust phylogenetic hypothesis for the T. phyllosoma species group. To identify likely instances of introgression prior to molecular scrutiny, we assessed biogeographic data and information on fertility of inter‐specific crosses. We first derived a few explicit hybridization hypotheses by considering the degree of spatial overlap within each species pair. Then, we assessed the plausibility of these hypotheses in the light of each species pair's cross‐fertility. Using this contextual information, we evaluated mito‐nuclear (cyt b, ITS‐2) gene‐tree incongruence and found evidence suggesting introgression within two species pairs. Finally, we modeled ILS using a Bayesian multispecies coalescent approach and either (a) a “complete” dataset with all the specimens in our sample, or (b) a “filtered” dataset without putatively introgressed specimens. The “filtered tree” had higher posterior‐probability support, as well as more plausible topology and divergence times, than the “complete tree.” Detecting and filtering out introgression and modeling ILS allowed us to derive an improved phylogenetic hypothesis for the T. phyllosoma species group. Our results illustrate how biogeographic and ecological‐reproductive contextual information can help clarify the systematics and evolution of recently diverged taxa prone to introgression and ILS.     

A data based parsimony method of cophylogenetic analysis

Phylogenies of closely interacting groups, such as hosts and parasites, are seldom completely congruent. Incongruence can arise from biologically meaningful differences in the histories of the two groups, or can be generated by artifactual differences that are merely the result of incorrect phylogenies with weakly supported nodes. We present a method that distinguishes between these sources of incongruence and identifies lineages that are responsible for significant differences between phylogenies. We use the logic of conditional combination in that we first test for statistically significant incongruence using the partition homogeneity test. Then we remove all possible combinations of taxa until a non-significant result of this test is achieved. Finally, we construct a 'combined evidence' phylogeny and then reposition the incongruent taxa. This method produces trees for final comparison using reconciliation methods, but it includes only as many incongruence events as can be statistically justified from the data sets. We apply this method to a host–parasite (gopher–louse) data set and identify many fewer incongruence events than do topology based analyses alone. Our method is broadly applicable to comparisons of phylogenies of interacting taxa, such as hosts and parasites, or mutualists. The method should also be useful for other problems involving comparisons of phylogenies, such as multiple gene trees or cladistic biogeography.     

Dealing with incongruence in the quest for the species tree: a case study from the orchid genus Satyrium

A species tree was reconstructed for the mainly African terrestrial orchid genus Satyrium. Separate phylogenetic analysis of both plastid and ribosomal nuclear DNA sequences for 63 species, revealed extensive topological conflict. Here we describe a detailed protocol to deal with incongruence involving three steps: identifying incongruence and testing its significance, assessing the cause of incongruence, and reconstructing the species tree. The Incongruence Length Difference test revealed that many cases of incongruence were non-significant. For the remaining significant cases, results from taxon jack-knifing experiments and parametric bootstrap suggested that non-biological artefacts such as sparse taxon sampling and long-branch attraction could be excluded as causes for the observed incongruence. In order to evaluate biological causes, such as orthology/paralogy conflation, lineage sorting, and hybridization, the number of events was counted that needs to be invoked a-posteriori to explain the observed pattern. In most cases where incongruence was significant, this resulted in a similar number of events for each of these different causes. Only for the three species from south east Asia, that form a monophyletic clade, hybridization was favoured over the alternative causes. This conclusion is based on the large number of events that needs to be invoked, in order for either orthology/paralogy conflation or lineage sorting to have been the cause of the incongruence+morphological evidence. The final species tree presented here is the product of the combined analysis of plastid and ITS sequences for all non-incongruent species and a-posteriori grafting of the incongruent clades or accessions onto the tree.     

Congruence of Morphological and Molecular Phylogenies

When phylogenetic trees constructed from morphological and molecular evidence disagree (i.e. are incongruent) it has been suggested that the differences are spurious or that the molecular results should be preferred a priori. Comparing trees can increase confidence (congruence), or demonstrate that at least one tree is incorrect (incongruence). Statistical analyses of 181 molecular and 49 morphological trees shows that incongruence is greater between than within the morphological and molecular partitions, and this difference is significant for the molecular partition. Because the level of incongruence between a pair of trees gives a minimum bound on how much error is present in the two trees, our results indicate that the level of error may be underestimated by congruence within partitions. Thus comparisons between morphological and molecular trees are particularly useful for detecting this incongruence (spurious or otherwise). Molecular trees have higher average congruence than morphological trees, but the difference is not significant, and both within- and between-partition incongruence is much lower than expected by chance alone. Our results suggest that both molecular and morphological trees are, in general, useful approximations of a common underlying phylogeny and thus, when molecules and morphology clash, molecular phylogenies should not be considered more reliable a priori.     

Detecting phylogenetic incongruence using BIONJ: an improvement of the ILD test

The problem of testing for congruence between phylogenetic data has long been debated among phylogeneticists, but reaches a critical point with the availability of large amount of biological sequences. Notably in prokaryotes, where the amount of lateral transfers is believed to be important, the inference of phylogenies using multiple genes requires testing for incongruence before concatenating the genes. On another scale, incongruence tests can be used to detect recombination points within single gene alignments. The incongruence length difference test (ILD), based on parsimony, has been proved to be useful for finding incongruent data sets, but its application remains limited to small data sets for computational time reasons. Here, we have adapted the principle of ILD to the BIONJ algorithm. This algorithm is based on a tree length minimisation criterion and is suitable to replace parsimony in this test when used with uncorrected distance (model-free approach). We show that this new test, ILD-BIONJ, while being much faster, is often more accurate than the ILD test, especially when the alignments compared are simulated under different evolutionary models.     

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.     

SINEs of speciation: tracking lineages with retroposons

The value of short interspersed elements (SINEs) for diagnosing common ancestry is being expanded to examine the differential sorting of lineages through the course of speciation events. Because most SINEs are neutral markers of identical descent, are not precisely excised from the genome and have a known ancestral condition, they are advantageous for reconciling gene trees and species trees with minimal phylogenetic error. A population perspective on SINE evolution combined with coalescence theory provides a context for investigating the phenomenon of ancestral polymorphism and its role in producing incongruent SINE insertion patterns among multiple loci. Studies of human Alu repeats demonstrate the value of young polymorphic SINEs for assessing human genomic diversity and tracking ancient demographics of human populations, whereas incongruent insertion patterns revealed by older fixed SINE loci, such as those in African cichlid fishes, contain information that might help identify ancient radiations that are otherwise obscured by accumulated mutations in sequence data. Here, we review the utility of retroposons for inferring common ancestry, discuss limits to the method, and clarify confusion by providing examples from the literature that illustrate how discordant multi-locus insertion patterns of retroelements can indicate lineage-sorting events that should not be misinterpreted as phylogenetic noise.     

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.     

Phylogeny of the complex Vanguerieae (Rubiaceae) genera Fadogia, Rytigynia, and Vangueria with close relatives and a new circumscription of Vangueria

The phylogeny of the Vanguerieae genera Fadogia, Rytigynia, Vangueria, and closely related genera is investigated using nuclear rDNA ITS sequences and plastid trnT-F and rpl16 sequences. Individual and combined analyses reveal several strongly supported clades. There are indications that Fadogia, Rytigynia, Tapiphyllum, and Vangueria are para- or polyphyletic and only Multidentia is strongly supported as monophyletic. Several taxa are found to have incongruent positions in the ITS and chloroplast phylogenies, and possible reasons behind these incongruencies are discussed. For Ancylanthos rubiginosus a chloroplast capture event can explain the incongruent position. In the Fadogia-Rytigynia group the incongruence is more widespread and cannot be attributed to a single taxon or a few taxa, but hybridization and introgression is the most likely explanation for the incongruence. It is concluded that the genera Ancylanthos, Lagynias, Pachystigma, and Tapiphyllum and the three species Fadogia agrestis, Rytigynia fuscosetulosa, and Rytigynia induta should be transferred to Vangueria.     

A priori assumptions about characters as a cause of incongruence between molecular and morphological hypotheses of primate interrelationships

When molecules and morphology produce incongruent hypotheses of primate interrelationships, the data are typically viewed as incompatible, and molecular hypotheses are often considered to be better indicators of phylogenetic history. However, it has been demonstrated that the choice of which taxa to include in cladistic analysis as well as assumptions about character weighting, character state transformation order, and outgroup choice all influence hypotheses of relationships and may positively influence tree topology, so that relationships between extant taxa are consistent with those found using molecular data. Thus, the source of incongruence between morphological and molecular trees may lie not in the morphological data themselves but in assumptions surrounding the ways characters evolve and their impact on cladistic analysis. In this study, we investigate the role that assumptions about character polarity and transformation order play in creating incongruence between primate phylogenies based on morphological data and those supported by multiple lines of molecular data. By releasing constraints imposed on published morphological analyses of primates from disparate clades and subjecting those data to parsimony analysis, we test the hypothesis that incongruence between morphology and molecules results from inherent flaws in morphological data. To quantify the difference between incongruent trees, we introduce a new method called branch slide distance (BSD). BSD mitigates many of the limitations attributed to other tree comparison methods, thus allowing for a more accurate measure of topological similarity. We find that releasing a priori constraints on character behavior often produces trees that are consistent with molecular trees. Case studies are presented that illustrate how congruence between molecules and unconstrained morphological data may provide insight into issues of polarity, transformation order, homology, and homoplasy.     

Rampant horizontal transfer and duplication of rubisco genes in eubacteria and plastids

Previous work has shown that molecular phylogenies of plastids, cyanobacteria, and proteobacteria based on the rubisco (ribulose-1,5- bisphosphate carboxylase/oxygenase) genes rbcL and rbcS are incongruent with molecular phylogenies based on other genes and are also incompatible with structural and biochemical information. Although it has been much speculated that this is the consequence of a single horizontal gene transfer (of a proteobacterial or mitochondrial rubisco operon into plastids of rhodophytic and chromophytic algae), neither this hypothesis nor the alternative hypothesis of ancient gene duplication have been examined in detail. We have conducted phylogenetic analyses of all available bacterial rbcL sequences, and representative plastid sequences, in order to explore these alternative hypothesis and fully examine the complexity of rubisco gene evolution. The rbcL phylogeny reveals a surprising number of gene relationships that are fundamentally incongruent with organismal relationships as inferred from multiple lines of other molecular evidence. On the order of six horizontal gene transfers are implied by the form I (L8S8) rbcL phylogeny, two between cyanobacteria and proteobacteria, one between proteobacteria and plastids, and three within proteobacteria. Alternatively, a single ancient duplication of the form I rubisco operon, followed by repeated and pervasive differential loss of one operon or the other, would account for much of this incongruity. In all probability, the rubisco operon has undergone multiple events of both horizontal gene transfer and gene duplication in different lineages.      

A model of horizontal gene transfer and the bacterial phylogeny problem

How much horizontal gene transfer (HGT) between species influences bacterial phylogenomics is a controversial issue. This debate, however, lacks any quantitative assessment of the impact of HGT on phylogenies and of the ability of tree-building methods to cope with such events. I introduce a Markov model of genome evolution with HGT, accounting for the constraints on time -- an HGT event can only occur between concomitantly living species. This model is used to simulate multigene sequence data sets with or without HGT. The consequences of HGT on phylogenomic inference are analyzed and compared to other well-known phylogenetic artefacts. It is found that supertree methods are quite robust to HGT, keeping high levels of performance even when gene trees are largely incongruent with each other. Gene tree incongruence per se is not indicative of HGT. HGT, however, removes the (otherwise observed) positive relationship between sequence length and gene tree congruence to the estimated species tree. Surprisingly, when applied to a bacterial and a eukaryotic multigene data set, this criterion rejects the HGT hypothesis for the former, but not the latter data set.     

Hybridization and polyploidy of an aquatic plant, Ruppia (Ruppiaceae), inferred from plastid and nuclear DNA phylogenies

? Premise of the study: The monogeneric family Ruppiaceae is found primarily in brackish water and is widely distributed on all continents, many islands, and from subartic to tropical zones. Ruppia taxonomy has been confusing because of its simplified morphology yet high phenotypic plasticity and the existence of polyploidy and putative hybrids. This study addresses the current classification of species in the genus, the origin of putative hybrids and polyploids, and the distribution of Ruppia species. ? Methods: Separate molecular phylogenetic analyses using plastid DNA and nuclear-encoded PHYB data sets were performed after chromosome observations. ? Key results: The resultant trees were largely congruent between genomes, but were incongruent in two respects: the first incongruence may be caused by long outgroup branches and their effect on ingroup rooting, and the second is caused by the existence of heterogeneous PHYB sequences for several accessions that may reflect several independent hybridization events. Several morphological species recognized in previous taxonomic revisions appear paraphyletic in plastid DNA and PHYB trees. ? Conclusions: Given the molecular phylogenies, and considering chromosome number and morphology, three species and one species complex comprising six lineages were discerned. A putative allotriploid, an allotetraploid, and a lineage of hybrid origin were identified within the species complex, and a hybrid was found outside the species complex, and their respective putative parental taxa were inferred. With respect to biogeography, a remarkably discontinuous distribution was identified in two cases, for which bird-mediated seed dispersal may be a reasonable explanation.     

Testing the utility of internally transcribed spacer sequences in coral phylogenetics

Reef-building corals often possess high levels of intraindividual and intraspecific ribosomal DNA (rDNA) variation that is largely polyphyletic between closely related species. Polyphyletic rDNA phylogenies coupled with high intraindividual rDNA variation have been taken as evidence of introgressive hybridization in corals. Interpreting the data is problematic because the rDNA cluster evolves in a complex fashion and polyphyletic lineages can be generated by a variety of processes--such as incomplete lineage sorting and slow concerted evolution--in addition to hybridization. Using the genetically characterized Caribbean Acropora hybridization system, we evaluate how well rDNA data perform in revealing patterns of recent introgressive hybridization in contrast to genetic data from four single-copy loci. While the rDNA data are broadly consistent with the unidirectional introgression seen in other loci, we show that the phylogenetic signature of recent introgressive hybridization is obscured in the Caribbean Acropora by ancient shared rDNA lineages that predate the divergence of the species.     

Phylogenetic inconsistency caused by ancient sex-biased gene migration

Inferences of ancient sex-biased migration patterns using sex-linked genetic markers are usually difficult because of a stochastic process of allele fixation. Nevertheless, incongruent phylogenetic trees between different sex-linked markers and between sex-linked and autosomal markers are frequently interpreted as a signature of sex-biased migration without further statistical evaluation. I investigated the types of incongruent phylogenetic trees from which past sex-biased migration events can be statistically supported under the coalescent model. In the case of mammals, detecting a sex-biased migration pattern is not guaranteed by comparing the phylogenetic pattern of mitochondrial and Y-chromosomal loci. Likewise, evidence of introgression at a mitochondrial locus, but not at autosomal loci, does not support the hypothesis of an ancient female-biased migration pattern with statistical significance. In contrast, evidence of introgression at ≥ 5 unlinked autosomal loci, but not at a Y-chromosomal locus, would reject the null hypothesis of a sexually equal migration rate with statistical significance. A similar argument can be made to infer a male-biased migration pattern. Furthermore, the investigation of many recombining sex-biased markers such as X-chromosomal loci in mammals has the potential to efficiently detect ancient sex-biased demographic patterns.