Phylogenomic assessment of the role of hybridization and introgression in trait evolution

Trait evolution among a set of species—a central theme in evolutionary biology—has long been understood and analyzed with respect to a species tree. However, the field of phylogenomics, which has been propelled by advances in sequencing technologies, has ushered in the era of species/gene tree incongruence and, consequently, a more nuanced understanding of trait evolution. For a trait whose states are incongruent with the branching patterns in the species tree, the same state could have arisen independently in different species (homoplasy) or followed the branching patterns of gene trees, incongruent with the species tree (hemiplasy). Another evolutionary process whose extent and significance are better revealed by phylogenomic studies is gene flow between different species. In this work, we present a phylogenomic method for assessing the role of hybridization and introgression in the evolution of polymorphic or monomorphic binary traits. We apply the method to simulated evolutionary scenarios to demonstrate the interplay between the parameters of the evolutionary history and the role of introgression in a binary trait’s evolution (which we call xenoplasy). Very importantly, we demonstrate, including on a biological data set, that inferring a species tree and using it for trait evolution analysis in the presence of gene flow could lead to misleading hypotheses about trait evolution.     

Evolution of myrmecophytism in western Malesian Macaranga (Euphorbiaceae)

Plants inhabited by ants (myrmecophytes) have evolved in a diversity of tropical plant lineages. Macaranga includes approximately 300 paleotropical tree species; in western Malesia there are 26 myrmecophytic species that vary in morphological specializations for ant association. The origin and diversification of myrmecophytism in Macaranga was investigated using phylogenetic analyses of morphological and nuclear ITS DNA characters and studies of character evolution. Despite low ITS variation, the combined analysis resulted in a well-supported hypothesis of relationships. Mapping myrmecophytism on all most parsimonious trees resulting from the combined analysis indicated that the trait evolved independently between two and four times and was lost between one and three times (five changes). This hypothesis was robust when tested against trees constrained to have three or fewer evolutionary transformations, although increased taxon sampling for the ITS analysis is required to confirm this. Mapping morphological traits on the phylogeny indicated that myrmecophytism was not homologous among lineages; each independent origin involved a suite of different specializations for ant-plant association. There was no evidence that myrmecophytic traits underwent sequential change through evolution; self-hollowing domatia evolved independently from ant-excavated domatia, and different food-body production types evolved in different lineages. The multiple origins of myrmecophytism in Macaranga were restricted to one small, exclusively western Malesian lineage of an otherwise large and nonmyrmecophytic genus. Although the evolution of aggregated food-body production and the formation of domatia coincided with the evolution of myrmecophytism in all cases, several morphological, ecological, and biogeographic factors appear to have facilitated and constrained this radiation of ant-plants.     

LIKELIHOOD OF ANCESTOR STATES IN ADAPTIVE RADIATION

Theories of ecological diversification make predictions about the timing and ordering of character state changes through history. These theories are testable by “reconstructing” ancestor states using phylogenetic trees and measurements of contemporary species. Here we use maximum likelihood to estimate and evaluate the accuracy of ancestor reconstructions. We present likelihoods of discrete ancestor states and derive probability distributions for continuous ancestral traits. The methods are applied to several examples: diets of ancestral Darwin's finches; origin of inquilinism in gall wasps; microhabitat partitioning and body size evolution in scrubwrens; digestive enzyme evolution in artiodactyl mammals; origin of a sexually selected male trait, the sword, in platies and swordtails; and evolution of specialization in Anolis lizards. When changes between discrete character states are rare, the maximum-likelihood results are similar to parsimony estimates. In this case the accuracy of estimates is often high, with the exception of some nodes deep in the tree. If change is frequent then reconstructions are highly uncertain, especially of distant ancestors. Ancestor states for continuous traits are typically highly uncertain. We conclude that measures of uncertainty are useful and should always be provided, despite simplistic assumptions about the probabilistic models that underlie them. If uncertainty is too high, reconstruction should be abandoned in favor of approaches that fit different models of trait evolution to species data and phylogenetic trees, taking into account the range of ancestor states permitted by the data.     

Evolution and losses of spines in slug caterpillars (Lepidoptera: Limacodidae)

Larvae of the cosmopolitan family Limacodidae, commonly known as “slug” caterpillars, are well known because of the widespread occurrence of spines with urticating properties, a morpho‐chemical adaptive trait that has been demonstrated to protect the larvae from natural enemies. However, while most species are armed with rows of spines (“nettle” caterpillars), slug caterpillars are morphologically diverse with some species lacking spines and thus are nonstinging. It has been demonstrated that the evolution of spines in slug caterpillars may have a single origin and that this trait is possibly derived from nonstinging slug caterpillars, but these conclusions were based on limited sampling of mainly New World taxa; thus, the evolution of spines and other traits within the family remains unresolved. Here, we analyze morphological variation in slug caterpillars within an evolutionary framework to determine character evolution of spines with samples from Asia, Australia, North America, and South America. The phylogeny of the Limacodidae was reconstructed based on a multigene dataset comprising five molecular markers (5.6 Kbp: COI, 28S, 18S, EF‐1α, and wingless) representing 45 species from 40 genera and eight outgroups. Based on this phylogeny, we infer that limacodids evolved from a common ancestor in which the larval type possessed spines, and then slug caterpillars without spines evolved independently multiple times in different continents. While larvae with spines are well adapted to avoiding generalist predators, our results imply that larvae without spines may be suited to different ecological niches. Systematic relationships of our dataset indicate six major lineages, several of which have not previously been identified.     

Evaluating evolutionary history in the face of high gene tree discordance in Australian Gehyra (Reptilia: Gekkonidae)

Species tree methods have provided improvements for estimating species relationships and the timing of diversification in recent radiations by allowing for gene tree discordance. Although gene tree discordance is often observed, most discordance is attributed to incomplete lineage sorting rather than other biological phenomena, and the causes of discordance are rarely investigated. We use species trees from multi-locus data to estimate the species relationships, evolutionary history and timing of diversification among Australian Gehyra—a group renowned for taxonomic uncertainty and showing a large degree of gene tree discordance. We find support for a recent Asian origin and two major clades: a tropically adapted clade and an arid adapted clade, with some exceptions, but no support for allopatric speciation driven by chromosomal rearrangement in the group. Bayesian concordance analysis revealed high gene tree discordance and comparisons of Robinson–Foulds distances showed that discordance between gene trees was significantly higher than that generated by topological uncertainty within each gene. Analysis of gene tree discordance and incomplete taxon sampling revealed that gene tree discordance was high whether terminal taxon or gene sampling was maximized, indicating discordance is due to biological processes, which may be important in contributing to gene tree discordance in many recently diversified organisms.     

When should we expect early bursts of trait evolution in comparative data? Predictions from an evolutionary food web model

Conceptual models of adaptive radiation predict that competitive interactions among species will result in an early burst of speciation and trait evolution followed by a slowdown in diversification rates. Empirical studies often show early accumulation of lineages in phylogenetic trees, but usually fail to detect early bursts of phenotypic evolution. We use an evolutionary simulation model to assemble food webs through adaptive radiation, and examine patterns in the resulting phylogenetic trees and species' traits (body size and trophic position). We find that when foraging trade-offs result in food webs where all species occupy integer trophic levels, lineage diversity and trait disparity are concentrated early in the tree, consistent with the early burst model. In contrast, in food webs in which many omnivorous species feed at multiple trophic levels, high levels of turnover of species' identities and traits tend to eliminate the early burst signal. These results suggest testable predictions about how the niche structure of ecological communities may be reflected by macroevolutionary patterns.     

Power of the concentrated changes test for correlated evolution

The concentrated changes test (CCT) calculates the probability that changes in a binary character are distributed randomly on the branches of a cladogram. This test is used to examine hypotheses of correlated evolution, especially cases where changes in the state of one character influence changes in the state of another character. The test may be sensitive to the addition of branches that lack either trait of interest (white branches). To examine the effects of the proportion of white branches and of tree topology on the CCT probability, we conducted a simulation analysis using a series of randomly generated 100-taxon trees, in addition to a nearly perfectly balanced (symmetrical) and a completely imbalanced (asymmetrical) 100-taxon tree. Using two models of evolution (gains only, or gains and losses), we evolved character pairs randomly onto these trees to simulate cases where (1) characters evolve independently (i.e., no correlation among the traits) or (2) all changes in the dependent character occur on branches containing the independent trait (i.e., a strong correlation among the traits). This allowed us to evaluate the sensitivity of the CCT to type I and type II errors, respectively. In the simulations, the CCT did not appear to be overly sensitive to the inclusion of white branches (low likelihood of type I error with both CCT probabilities < 0.05 and < 0.01). However, the CCT was susceptible to type II error when the proportion of white branches was < 20%. The test was also sensitive to tree shape and was positively correlated to Colless's tree imbalance statistic I. Finally, the CCT responded differently for simulations where only gains were allowed and those where both gains and losses were permitted. These results indicate that the CCT is unlikely to detect a correlation between characters when no such correlation exists. However, when a trait can be gained but not lost, the CCT is conservative and may fail to detect true correlations among traits (increased type II error). Determination of the sampling universe (the taxa included in the comparative analysis) can strongly influence the probability of making such type II errors. We suggest guidelines to circumvent these limitations.     

PHYLOGENIES AND THE ANALYSIS OF EVOLUTIONARY SEQUENCES,WITH EXAMPLES FROM SEED PLANTS

Studies of character evolution have frequently relied on ahistorical correlations rather than on phylogenies. However, correlations do not estimate the number of times that a trait evolved, and they are insensitive to the direction or the temporal sequence of character transformation. In contrast, cladograms can provide this information. A cladistic test of the hypothesis that the evolution of dioecy is favored in animal-dispersed plants indicates that dioecy may have originated somewhat more often in such lineages. Nevertheless, differences in rates of speciation or extinction must largely account for the observed species-level correlation between dispersal and breeding system. In considering the evolution of individual traits, cladograms help identify the context in which a feature evolved and specify which organisms should be compared in evaluating the causes of character change. Determining whether a feature and a performance advantage were strictly historically correlated or followed one another in sequence helps to distinguish whether the trait is an adaptation or an exaptation for the function. For example, cladograms of seed plants suggest that double fertilization arose incidentally prior to the origin of angiosperms and that the resulting product was later co-opted and elaborated as a nutritive tissue for the developing embryo. The order of character assembly in a lineage also bears on the evolution of functional and developmental interdependencies. In particular, it may be possible to trace the evolution of a character's “burden” from an initial period, during which change is more likely, through later stages, wherein successful modification is less likely owing to the evolution of dependent characters. The evolution of vessels and of floral phyllotaxis in angiosperms may exemplify this pattern. Recognition that the likelihood of character transformation may change during the evolution of a group warns against character weighting in phylogenetic analysis.     

Behavioural evolution in penguins does not reflect phylogeny

Over the past two decades, behavioural biologists and ecologists have made effective use of the comparative method, but have often stopped short of adopting an explicitly phylogenetic approach. We examined 68 behaviour and life history (BLH) traits of 15 penguin species to: (i) infer penguin phylogeny, (ii) assess homology of behavioural characters, and (iii) evaluate hypotheses about character evolution and ancestral states. Parsimony analysis of the BLH dataset found either two shortest trees (characters coded as unordered) or a single shortest tree (characters coded as a combination of unordered and Dollo). The BLH data had significant structure. Kishino–Hasegawa tests indicated that BLH trees were significantly different from most previous estimates of penguin phylogeny. The BLH phylogeny generated from Dollo characters appeared to be less accurate than the tree derived from the completely unordered dataset. Dividing BLH data into display and non‐display traits resulted in no significant differences in level of homoplasy and no difference in the accuracy of phylogeny. Tests for homology of BLH traits were performed by mapping the characters onto a molecular tree. Assuming that independent gains are less likely than losses of character states, 65 of the 68 characters were likely to be homologous across taxa, and at least several characters appeared to have been stable since the origin of modern penguins around 30 Myr. Finally, the likely BLH traits of the most recent common ancestor of extant penguins were reconstructed from character states along the internal branch leading to the penguins. This analysis suggested that the “proto‐penguin” probably had a similar life history to current temperate penguins but few ritualized behaviours. A southern, cool‐temperate origin of penguins is suggested.     

Comparing two Bayesian methods for gene tree/species tree reconstruction: simulations with incomplete lineage sorting and horizontal gene transfer

With the increasing interest in recognizing the discordance between gene genealogies, various gene tree/species tree reconciliation methods have been developed. We present here the first attempt to assess and compare two such Bayesian methods, Bayesian estimation of species trees (BEST) and BUCKy (Bayesian untangling of concordance knots), in the presence of several known processes of gene tree discordance. DNA alignments were simulated under the influence of incomplete lineage sorting (ILS) and of horizontal gene transfer (HGT). BEST and BUCKy both account for uncertainty in gene tree estimation but differ substantially in their assumptions of what caused gene tree discordance. BEST estimates a species tree using the coalescent model, assuming that all gene tree discordance is due to ILS. BUCKy does not assume any specific biological process of gene tree discordance through the use of a nonparametric clustering of concordant genes. BUCKy estimates the concordance factor (CF) of a clade, which is defined as the proportion of genes that truly have the clade in their trees. The estimated concordance tree is then built from clades with the highest estimated CFs. Because of their different assumptions, it was expected that BEST would perform better in the presence of ILS and that BUCKy would perform better in the presence of HGT. As expected, the species tree was more accurately reconstructed by BUCKy in the presence of HGT, when the HGT events were unevenly placed across the species tree. BUCKy and BEST performed similarly in most other cases, including in the presence of strong ILS and of HGT events that were evenly placed across the tree. However, BUCKy was shown to underestimate the uncertainty in CF estimation, with short credibility intervals. Despite this, the discordance pattern estimated by BUCKy could be compared with the signature of ILS. The resulting test for the adequacy of the coalescent model proved to have low Type I error. It was powerful when HGT was the major source of discordance and when HGT events were unevenly placed across the species tree.     

The structure of cynipid oak galls: patterns in the evolution of an extended phenotype

Galls are highly specialized plant tissues whose development is induced by another organism. The most complex and diverse galls are those induced on oak trees by gallwasps (Hymenoptera: Cynipidae: Cynipini), each species inducing a characteristic gall structure. Debate continues over the possible adaptive significance of gall structural traits; some protect the gall inducer from attack by natural enemies, although the adaptive significance of others remains undemonstrated. Several gall traits are shared by groups of oak gallwasp species. It remains unknown whether shared traits represent (i) limited divergence from a shared ancestral gall form, or (ii) multiple cases of independent evolution. Here we map gall character states onto a molecular phylogeny of the oak cynipid genus Andricus, and demonstrate three features of the evolution of gall structure: (i) closely related species generally induce galls of similar structure; (ii) despite this general pattern, closely related species can induce markedly different galls; and (iii) several gall traits (the presence of many larval chambers in a single gall structure, surface resins, surface spines and internal air spaces) of demonstrated or suggested adaptive value to the gallwasp have evolved repeatedly. We discuss these results in the light of existing hypotheses on the adaptive significance of gall structure.     

Maximum tree: a consistent estimator of the species tree

We propose a model based approach to use multiple gene trees to estimate the species tree. The coalescent process requires that gene divergences occur earlier than species divergences when there is any polymorphism in the ancestral species. Under this scenario, speciation times are restricted to be smaller than the corresponding gene split times. The maximum tree (MT) is the tree with the largest possible speciation times in the space of species trees restricted by available gene trees. If all populations have the same population size, the MT is the maximum likelihood estimate of the species tree. It can be shown the MT is a consistent estimator of the species tree even when the MT is built upon the estimates of the true gene trees if the gene tree estimates are statistically consistent. The MT converges in probability to the true species tree at an exponential rate.     

壳斗科雌花不同发育时期比较胚胎学特征初探

壳斗科是北半球暖温带森林的优势树种之一,具有重要的生态功能和经济开发价值。该研究根据文献资料和野外考察及标本观察,对壳斗科雌花不同发育时期的重要性状进行整理和对比分析,并在系统发育框架下对雌花重要性状的演化趋势进行分析,探讨了上述性状在系统发育树上的演化趋势。选取壳斗科6个属13种雌花发育的5个性状,利用Mesquite在依据OhManos基于核基因CRC和ITS构建的系统发育树上,将性状进行编码和赋值,其中合子形成滞后时间按照数量性状进行分析,其余性状按照质量性状进行分析,并用最大简约法重建雌花形态及受精特征的祖先性状。结果表明:壳斗科雌花的4个重要性状子房室无毛被、胚珠具明显珠柄、子房与花柱比值大和花粉在子房室内休眠是祖征;合子形成滞后时间的性状是平行演化特征。针对壳斗科胚囊成熟时珠心组织残余结构的解剖术语混乱的现象,建议以后统一采用"珠心柱"指示该结构。该研究结果揭示了壳斗科花及胚胎学性状的演化趋势,也为探讨这些性状演化提供基础资料。     

iGLASS: an improvement to the GLASS method for estimating species trees from gene trees

Several methods have been designed to infer species trees from gene trees while taking into account gene tree/species tree discordance. Although some of these methods provide consistent species tree topology estimates under a standard model, most either do not estimate branch lengths or are computationally slow. An exception, the GLASS method of Mossel and Roch, is consistent for the species tree topology, estimates branch lengths, and is computationally fast. However, GLASS systematically overestimates divergence times, leading to biased estimates of species tree branch lengths. By assuming a multispecies coalescent model in which multiple lineages are sampled from each of two taxa at L independent loci, we derive the distribution of the waiting time until the first interspecific coalescence occurs between the two taxa, considering all loci and measuring from the divergence time. We then use the mean of this distribution to derive a correction to the GLASS estimator of pairwise divergence times. We show that our improved estimator, which we call iGLASS, consistently estimates the divergence time between a pair of taxa as the number of loci approaches infinity, and that it is an unbiased estimator of divergence times when one lineage is sampled per taxon. We also show that many commonly used clustering methods can be combined with the iGLASS estimator of pairwise divergence times to produce a consistent estimator of the species tree topology. Through simulations, we show that iGLASS can greatly reduce the bias and mean squared error in obtaining estimates of divergence times in a species tree.     

Gene trees, species trees, and morphology converge on a similar phylogeny of living gars (Actinopterygii: Holostei: Lepisosteidae), an ancient clade of ray-finned fishes

Extant gars represent the remaining members of a formerly diverse assemblage of ancient ray-finned fishes and have been the subject of multiple phylogenetic analyses using morphological data. Here, we present the first hypothesis of phylogenetic relationships among living gar species based on molecular data, through the examination of gene tree heterogeneity and coalescent species tree analyses of a portion of one mitochondrial (COI) and seven nuclear (ENC1, myh6, plagl2, S7 ribosomal protein intron 1, sreb2, tbr1, and zic1) genes. Individual gene trees displayed varying degrees of resolution with regards to species-level relationships, and the gene trees inferred from COI and the S7 intron were the only two that were completely resolved. Coalescent species tree analyses of nuclear genes resulted in a well-resolved and strongly supported phylogenetic tree of living gar species, for which Bayesian posterior node support was further improved by the inclusion of the mitochondrial gene. Species-level relationships among gars inferred from our molecular data set were highly congruent with previously published morphological phylogenies, with the exception of the placement of two species, Lepisosteus osseus and L. platostomus. Re-examination of the character coding used by previous authors provided partial resolution of this topological discordance, resulting in broad concordance in the phylogenies inferred from individual genes, the coalescent species tree analysis, and morphology. The completely resolved phylogeny inferred from the molecular data set with strong Bayesian posterior support at all nodes provided insights into the potential for introgressive hybridization and patterns of allopatric speciation in the evolutionary history of living gars, as well as a solid foundation for future examinations of functional diversification and evolutionary stasis in a "living fossil" lineage.     

Discordance of species trees with their most likely gene trees

Because of the stochastic way in which lineages sort during speciation, gene trees may differ in topology from each other and from species trees. Surprisingly, assuming that genetic lineages follow a coalescent model of within-species evolution, we find that for any species tree topology with five or more species, there exist branch lengths for which gene tree discordance is so common that the most likely gene tree topology to evolve along the branches of a species tree differs from the species phylogeny. This counterintuitive result implies that in combining data on multiple loci, the straightforward procedure of using the most frequently observed gene tree topology as an estimate of the species tree topology can be asymptotically guaranteed to produce an incorrect estimate. We conclude with suggestions that can aid in overcoming this new obstacle to accurate genomic inference of species phylogenies.     

Efficient FPT Algorithms for (Strict) Compatibility of Unrooted Phylogenetic Trees

In phylogenetics, a central problem is to infer the evolutionary relationships between a set of species X; these relationships are often depicted via a phylogenetic tree—a tree having its leaves labeled bijectively by elements of X and without degree-2 nodes—called the “species tree.” One common approach for reconstructing a species tree consists in first constructing several phylogenetic trees from primary data (e.g., DNA sequences originating from some species in X), and then constructing a single phylogenetic tree maximizing the “concordance” with the input trees. The obtained tree is our estimation of the species tree and, when the input trees are defined on overlapping—but not identical—sets of labels, is called “supertree.” In this paper, we focus on two problems that are central when combining phylogenetic trees into a supertree: the compatibility and the strict compatibility problems for unrooted phylogenetic trees. These problems are strongly related, respectively, to the notions of “containing as a minor” and “containing as a topological minor” in the graph community. Both problems are known to be fixed parameter tractable in the number of input trees k, by using their expressibility in monadic second-order logic and a reduction to graphs of bounded treewidth. Motivated by the fact that the dependency on k of these algorithms is prohibitively large, we give the first explicit dynamic programming algorithms for solving these problems, both running in time \(2^{O(k^2)} \cdot n\), where n is the total size of the input.