Estimate of within population incremental selection through branch imbalance in lineage trees

Incremental selection within a population, defined as limited fitness changes following mutation, is an important aspect of many evolutionary processes. Strongly advantageous or deleterious mutations are detected using the synonymous to non-synonymous mutations ratio. However, there are currently no precise methods to estimate incremental selection. We here provide for the first time such a detailed method and show its precision in multiple cases of micro-evolution. The proposed method is a novel mixed lineage tree/sequence based method to detect within population selection as defined by the effect of mutations on the average number of offspring. Specifically, we propose to measure the log of the ratio between the number of leaves in lineage trees branches following synonymous and non-synonymous mutations. The method requires a high enough number of sequences, and a large enough number of independent mutations. It assumes that all mutations are independent events. It does not require of a baseline model and is practically not affected by sampling biases. We show the method''s wide applicability by testing it on multiple cases of micro-evolution. We show that it can detect genes and inter-genic regions using the selection rate and detect selection pressures in viral proteins and in the immune response to pathogens.     

Bayesian inference of character evolution

Much recent progress in evolutionary biology is based on the inference of ancestral states and past transformations in important traits on phylogenetic trees. These exercises often assume that the tree is known without error and that ancestral states and character change can be mapped onto it exactly. In reality, there is often considerable uncertainty about both the tree and the character mapping. Recently introduced Bayesian statistical methods enable the study of character evolution while simultaneously accounting for both phylogenetic and mapping uncertainty, adding much needed credibility to the reconstruction of evolutionary history.     

Flammability and serotiny as strategies: correlated evolution in pines

Fire may act as a selective force on plants both through its direct effects by killing or wounding susceptible individuals and through its effect on the environment: the post-fire environment may select specific physiological traits or life histories. We used phylogenetic independent contrasts to test the hypothesis that fire has selected for correlated evolution among alternative suites of traits in pines: a survival/avoidance suite characterized by thick bark, height, and self-pruning of dead branches; and a fire-embracing strategy in which plants invest little into survival, exhibit traits which enhance flammability, and use fire as a means to cue seedling establishment to the post-fire environment through serotinous cones. We created a set of alternative 'supertree' phylogenies for the genus Pinus from published sources. Using these alternative phylogenies, published ecological data for 38 pine species, and newly collected morphological data, we demonstrate that much variation in trait evolution occurs along a fire-surviving/fire-embracing axis. Pines vary in their susceptibility to ignition since a tree that retains dead branches is more likely to carry a fire into the canopy than a tree that self-prunes. The evolution of increased flammability may have altered evolutionary trajectories prompting an evolutionary switch from a fire-surviving to a fire-embracing life history. Alternatively, the fire-embracing strategy may in fact select for increased flammability to ensure canopy ignition and the realization of serotinous seed-release.     

Irrational exuberance for resolved species trees

Phylogenomics has largely succeeded in its aim of accurately inferring species trees, even when there are high levels of discordance among individual gene trees. These resolved species trees can be used to ask many questions about trait evolution, including the direction of change and number of times traits have evolved. However, the mapping of traits onto trees generally uses only a single representation of the species tree, ignoring variation in the gene trees used to construct it. Recognizing that genes underlie traits, these results imply that many traits follow topologies that are discordant with the species topology. As a consequence, standard methods for character mapping will incorrectly infer the number of times a trait has evolved. This phenomenon, dubbed “hemiplasy,” poses many problems in analyses of character evolution. Here we outline these problems, explaining where and when they are likely to occur. We offer several ways in which the possible presence of hemiplasy can be diagnosed, and discuss multiple approaches to dealing with the problems presented by underlying gene tree discordance when carrying out character mapping. Finally, we discuss the implications of hemiplasy for general phylogenetic inference, including the possible drawbacks of the widespread push for “resolved” species trees.     

PHYLOGENETIC ANALYSIS OF THE EVOLUTIONARY CORRELATION USING LIKELIHOOD

Many evolutionary processes can lead to a change in the correlation between continuous characters over time or on different branches of a phylogenetic tree. Shifts in genetic or functional constraint, in the selective regime, or in some combination thereof can influence both the evolution of continuous traits and their relation to each other. These changes can often be mapped on a phylogenetic tree to examine their influence on multivariate phenotypic diversification. We propose a new likelihood method to fit multiple evolutionary rate matrices (also called evolutionary variance–covariance matrices) to species data for two or more continuous characters and a phylogeny. The evolutionary rate matrix is a matrix containing the evolutionary rates for individual characters on its diagonal, and the covariances between characters (of which the evolutionary correlations are a function) elsewhere. To illustrate our approach, we apply the method to an empirical dataset consisting of two features of feeding morphology sampled from 28 centrarchid fish species, as well as to data generated via phylogenetic numerical simulations. We find that the method has appropriate type I error, power, and parameter estimation. The approach presented herein is the first to allow for the explicit testing of how and when the evolutionary covariances between characters have changed in the history of a group.     

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.     

Positive selection of synonymous mutations in vesicular stomatitis virus

Prevailing evolutionary forces are typically deduced from the pattern of differences in synonymous and non-synonymous mutations, under the assumption of neutrality in the absence of amino acid change. We determined the complete sequence of ten vesicular stomatitis virus populations evolving under positive selection. A significant number of the mutations occurred independently in two or more strains, a process known as parallel evolution, and a substantial fraction of the parallel mutations were silent. Parallel evolution was also identified in non-coding regions. These results indicate that silent mutations can significantly contribute to adaptation in RNA viruses, and relative frequencies of synonymous and non-synonymous substitutions may not be useful to resolve their evolutionary history.     

Episodic sitewise positive selection on the signal recognition particle protein Ffh in Actinobacteria

In this work, we report a case of episodic sitewise positive selection acting on the highly conserved SRP protein Ffh in Actinobacteria. An elevated non-synonymous to synonymous mutation ratio (ω) was observed along the non-terminal branches of the species tree, which contained 11 Actinobacteria species, where positively selected residues were frequently observed within the signal sequence-binding domain. Together with the estimated lineage-specific ω ratios for several core components in the major protein translocation systems, our data suggest that this positive selection might be partially driven by the diversification of signal sequences.     

Using Plant Functional Traits and Phylogenies to Understand Patterns of Plant Community Assembly in a Seasonal Tropical Forest in Lao PDR

Plant functional traits reflect different evolutionary responses to environmental variation, and among extant species determine the outcomes of interactions between plants and their environment, including other plant species. Thus, combining phylogenetic and trait-based information can be a powerful approach for understanding community assembly processes across a range of spatial scales. We used this approach to investigate tree community composition at Phou Khao Khouay National Park (18°14’-18°32’N; 102°38’- 102°59’E), Laos, where several distinct forest types occur in close proximity. The aim of our study was to examine patterns of plant community assembly across the strong environmental gradients evident at our site. We hypothesized that differences in tree community composition were being driven by an underlying gradient in soil conditions. Thus, we predicted that environmental filtering would predominate at the site and that the filtering would be strongest on sandier soil with low pH, as these are the conditions least favorable to plant growth. We surveyed eleven 0.25 ha (50x50 m) plots for all trees above 10 cm dbh (1221 individual trees, including 47 families, 70 genera and 123 species) and sampled soils in each plot. For each species in the community, we measured 11 commonly studied plant functional traits covering both the leaf and wood economic spectrum traits and we reconstructed a phylogenetic tree for 115 of the species in the community using rbcL and matK sequences downloaded from Genebank (other species were not available). Finally we compared the distribution of trait values and species at two scales (among plots and 10x10m subplots) to examine trait and phylogenetic community structures. Although there was strong evidence that an underlying soil gradient was determining patterns of species composition at the site, our results did not support the hypothesis that the environmental filtering dominated community assembly processes. For the measured plant functional traits there was no consistent pattern of trait dispersion across the site, either when traits were considered individually or when combined in a multivariate analysis. However, there was a significant correlation between the degree of phylogenetic dispersion and the first principle component axis (PCA1) for the soil parameters. Moreover, the more phylogenetically clustered plots were on sandier soils with lower pH. Hence, we suggest that the community assembly processes across our site may reflect the influence of more conserved traits that we did not measure. Nevertheless, our results are equivocal and other interpretations are possible. Our study illustrates some difficulties in combining trait and phylogenetic approaches that may result from the complexities of integrating spatial and evolutionary processes that vary at different scales.     

Evolutionary change of restriction cleavage sites and phylogenetic inference for man and apes

A mathematical theory for the evolutionary change of restriction endonuclease cleavage sites is developed, and the probabilities of various types of restriction-site changes are evaluated. A computer simulation is also conducted to study properties of the evolutionary change of restriction sites. These studies indicate that parsimony methods of constructing phylogenetic trees often make erroneous inferences about evolutionary changes of restriction sites unless the number of nucleotide substitutions per site is less than 0.01 for all branches of the tree. This introduces a systematic error in estimating the number of mutational changes for each branch and, consequently, in constructing phylogenetic trees. Therefore, parsimony methods should be used only in cases where nucleotide sequences are closely related. Reexamination of Ferris et al.'s data on restriction-site differences of mitochondrial DNAs does not support Templeton's conclusions regarding the phylogenetic tree for man and apes and the molecular clock hypothesis. Templeton's claim that Nei and Li's method of estimating the number of nucleotide substitutions per site is seriously affected by parallel losses and loss-gains of restriction sites is also unsupported.      

Codon-based tests of positive selection,branch lengths,and the evolution of mammalian immune system genes

Using basic probability theory, we show that there is a substantial likelihood that even in the presence of strong purifying selection, there will be a number of codons in which the number of synonymous nucleotide substitutions per site (d (S)) exceeds the number of non-synonymous nucleotide substitutions per site (d (N)). In an empirical study, we examined the numbers of synonymous (b (S)) and non-synonymous substitutions (b (N)) along branches of the phylogenies of 69 single-copy orthologous genes from seven species of mammals. A pattern of b (N) > b (S) was most commonly seen in the shortest branches of the tree and was associated with a high coefficient of variation in both b (N) and b (S), suggesting that high stochastic error in b (N) and b (S) on short branches, rather than positive Darwinian selection, is the explanation of most cases where b (N) is greater than b (S) on a given branch. The branch-site method of Zhang et al. (Zhang, Nielsen, Yang, Mol Biol Evol, 22:2472-2479, 2005) identified 117 codons on 35 branches as "positively selected," but a majority of these codons lacked synonymous substitutions, while in the others, synonymous and non-synonymous differences per site occurred in approximately equal frequencies. Thus, it was impossible to rule out the hypothesis that chance variation in the pattern of mutation across sites, rather than positive selection, accounted for the observed pattern. Our results showed that b (N)/b (S) was consistently elevated in immune system genes, but neither the search for branches with b (N) > b (S) nor the branch-site method revealed this trend.     

Evaluation of methods for determination of a reconstructed history of gene sequence evolution

With whole-genome sequences being completed at an increasing rate, it is important to develop and assess tools to analyze them. Following annotation of the protein content of a genome, one can compare sequences with previously characterized homologous genes to detect novel functions within specific proteins in the evolution of the newly sequenced genome. One common statistical method to detect such changes is to compare the ratios of nonsynonymous (K(a)) to synonymous (K(s)) nucleotide substitution rates. Here, the effects of several parameters that can influence this calculation (sequence reconstruction method, phylogenetic tree branch length weighting, GC content, and codon bias) are examined. Also, two new alternative measures of adaptive evolution, the point accepted mutations (PAM)/neutral evolutionary distance (NED) ratio and the sequence space assessment (SSA) statistic are presented. All of these methods are compared using two sequence families: the recent divergence of leptin orthologs in primates, and the more ancient divergence of the deoxyribonucleoside kinase family. The examination of these and other measures to detect changes of gene function along branches of a phylogenetic tree will become increasingly important in the postgenomic era.     

Environmental variation is a major predictor of global trait turnover in mammals

Aim

To evaluate how environment and evolutionary history interact to influence global patterns of mammal trait diversity (a combination of 14 morphological and life‐history traits).

Location

The global terrestrial environment.

Taxon

Terrestrial mammals.

Methods

We calculated patterns of spatial turnover for mammalian traits and phylogenetic lineages using the mean nearest taxon distance. We then used a variance partitioning approach to establish the relative contribution of trait conservatism, ecological adaptation and clade specific ecological preferences on global trait turnover.

Results

We provide a global scale analysis of trait turnover across mammalian terrestrial assemblages, which demonstrates that phylogenetic turnover by itself does not predict trait turnover better than random expectations. Conversely, trait turnover is consistently more strongly associated with environmental variation than predicted by our null models. The influence of clade‐specific ecological preferences, reflected by the shared component of phylogenetic turnover and environmental variation, was considerably higher than expectations. Although global patterns of trait turnover are dependent on the trait under consideration, there is a consistent association between trait turnover and environmental predictive variables, regardless of the trait considered.

Main conclusions

Our results suggest that changes in phylogenetic composition are not always coupled with changes in trait composition on a global scale and that environmental conditions are strongly associated with patterns of trait composition across species assemblages, both within and across phylogenetic clades.     

Plant persistence traits in fire-prone ecosystems of the Mediterranean basin: a phylogenetic approach

The two main fire response traits found in the Mediterranean basin are the resprouting capacity (R) and the propagule-persistence capacity (P). Previous studies suggested that these two traits might be correlated. In this paper we first test whether R and P have evolved independently. Then, we ask if the correlation occurs because (a) one trait is not the target of selection but it is genetically linked to the other trait which is the one under selection pressure (indirect selection), or (b) because different evolutionary responses to the same selective pressure are acting in parallel on populations at different genetic starting points (parallel selection). Finally, we test to what extent resprouting is associated with some vegetative and reproductive traits.
To answer these questions we used a traits database for the eastern Iberian Peninsula and we assembled the phylogenetic tree on the basis of published information. The results indicate that the two traits are negatively associated and support the parallel selection scenario in which changes in R precedes changes in P. The phylogenetic–informed associations of resprouting with other traits (plant height, age at maturity) support the existence of allocation tradeoffs.
The results are consistent with the biogeographical history of the Mediterranean basin flora where most of lineages already resprouted to persist after a disturbance during the Tertiary, thus making it improbable that an additional costly persistence strategy would evolve under the Quaternary climatic conditions.     

Phylocom: software for the analysis of phylogenetic community structure and trait evolution

MOTIVATION: The increasing availability of phylogenetic and trait data for communities of co-occurring species has created a need for software that integrates ecological and evolutionary analyses. Capabilities: Phylocom calculates numerous metrics of phylogenetic community structure and trait similarity within communities. Hypothesis testing is implemented using several null models. Within the same framework, it measures phylogenetic signal and correlated evolution for species traits. A range of utility functions allow community and phylogenetic data manipulation, tree and trait generation, and integration into scientific workflows. Availability: Open source at: http://phylodiversity.net/phylocom/.     

Spatial patterns of wood traits in China are controlled by phylogeny and the environment

Aim Wood properties are related to tree physiology and mechanical stability and are influenced by both phylogeny and the environment. However, it remains unclear to what extent geographical gradients in wood traits are shaped by either phylogeny or the environment. Here we aimed to disentangle the influences of phylogeny and the environment on spatial trends in wood traits. Location China. Methods We compiled a data set of 11 wood properties for 618 tree species from 98 sampling sites in China to assess their phylogenetic and spatial patterns, and to determine how many of the spatial patterns in wood properties are attributable to the environment after correction for phylogenetic influences. Result All wood traits examined exhibited significant phylogenetic signal. The widest divergence in wood traits was observed between gymnosperms and angiosperms, Rosids and Asterids, Magnoiliids and Eudicots, and in Lamiales. For most wood traits, the majority of trait variation was observed at genus and species levels. The mechanical properties of wood showed correlated evolution with wood density. Most of the mechanical properties of wood exhibited significant latitudinal variation but limited or no altitudinal variation, and were positively correlated with mean annual precipitation based on both Pearson's correlation analysis and the phylogenetic comparative method. Correlations at family level between mean annual temperature and wood density, compression strength, cross‐section hardness, modulus of elasticity and volumetric shrinkage coefficient became significant after correction for phylogenetic influences. Main conclusions Phylogeny interacted with the environment in shaping the spatial patterns of wood traits of trees across China because most wood properties showed strong phylogenetic conservatism and thus affected environmental tolerances and distributions of tree species. Mean annual precipitation was a key environmental factor explaining the spatial patterns of wood traits. Our study provides valuable insights into the geographical patterns in productivity, distribution and ecological strategy of trees linking to wood traits.     

Host functional and phylogenetic composition rather than host diversity structure plant–herbivore networks

Declining plant diversity alters ecological networks, such as plant–herbivore interactions. However, our knowledge of the potential mechanisms underlying effects of plant species loss on plant–herbivore network structure is still limited. We used DNA barcoding to identify herbivore–host plant associations along declining levels of tree diversity in a large‐scale, subtropical biodiversity experiment. We tested for effects of tree species richness, host functional and phylogenetic diversity, and host functional (leaf trait) and phylogenetic composition on species, phylogenetic and network composition of herbivore communities. We found that phylogenetic host composition and related palatability/defence traits but not tree species richness significantly affected herbivore communities and interaction network complexity at both the species and community levels. Our study indicates that evolutionary dependencies and functional traits of host plants determine the composition of higher trophic levels and corresponding interaction networks in species‐rich ecosystems. Our findings highlight that characteristics of the species lost have effects on ecosystem structure and functioning across trophic levels that cannot be predicted from mere reductions in species richness.     

Fire-adapted traits of Pinus arose in the fiery Cretaceous

? The mapping of functional traits onto chronograms is an emerging approach for the identification of how agents of natural selection have shaped the evolution of organisms. Recent research has reported fire-dependent traits appearing among flowering plants from 60 million yr ago (Ma). Although there are many records of fossil charcoal in the Cretaceous (65-145 Ma), evidence of fire-dependent traits evolving in that period is lacking. ? We link the evolutionary trajectories for five fire-adapted traits in Pinaceae with paleoatmospheric conditions over the last 250 million yr to determine the time at which fire originated as a selective force in trait evolution among seed plants. ? Fire-protective thick bark originated in Pinus c. 126 Ma in association with low-intensity surface fires. More intense crown fires emerged c. 89 Ma coincident with thicker bark and branch shedding, or serotiny with branch retention as an alternative strategy. These innovations appeared at the same time as the Earth's paleoatmosphere experienced elevated oxygen levels that led to high burn probabilities during the mid-Cretaceous. ? The fiery environments of the Cretaceous strongly influenced trait evolution in Pinus. Our evidence for a strong correlation between the evolution of fire-response strategies and changes in fire regime 90-125 Ma greatly backdates the key role that fire has played in the evolution of seed plants.     

Joint effect of phylogenetic relatedness and trait selection on the elevational distribution of Rhododendron species

Congeneric species may coexist at fine spatial scales through niche differentiation, however, the magnitude to which the effects of functional traits and phylogenetic relatedness contribute to their distribution along elevational gradients remains understudied. To test the hypothesis that trait and elevational range overlap can affect local speciesʼ coexistence, we first compared phylogenetic relatedness and trait (including morphological traits and leaf elements) divergence among closely related species of Rhododendron L. on Yulong Mountain, China. We then assessed relationships between the overlap of multiple functional traits and the degree of elevational range overlap among species pairs in a phylogenetic context. We found that phylogeny was a good predictor for most functional traits, where closely related species showed higher trait similarity and occupied different elevational niches at our study site. Species pairs of R. subgen. Hymenanthes (Blume) K. Koch showed low elevational range overlap and some species pairs of R. subgen. Rhododendron showed obvious niche differentiation. Trait divergence is greater for species in R. subgen. Rhododendron, and it plays an important role between species pairs with low elevational range overlap. Trait convergent selection takes place between co-occurring closely related species that have high elevational range overlap, which share more functional trait space due to environmental filtering or ecological adaptation in more extreme habitats. Our results highlight the importance of evolutionary history and trait selection for species coexistence at fine ecological scales along environmental gradients.     

Independent evolution of leaf and root traits within and among temperate grassland plant communities

In this study, we used data from temperate grassland plant communities in Alberta, Canada to test two longstanding hypotheses in ecology: 1) that there has been correlated evolution of the leaves and roots of plants due to selection for an integrated whole-plant resource uptake strategy, and 2) that trait diversity in ecological communities is generated by adaptations to the conditions in different habitats. We tested the first hypothesis using phylogenetic comparative methods to test for evidence of correlated evolution of suites of leaf and root functional traits in these grasslands. There were consistent evolutionary correlations among traits related to plant resource uptake strategies within leaf tissues, and within root tissues. In contrast, there were inconsistent correlations between the traits of leaves and the traits of roots, suggesting different evolutionary pressures on the above and belowground components of plant morphology. To test the second hypothesis, we evaluated the relative importance of two components of trait diversity: within-community variation (species trait values relative to co-occurring species; α traits) and among-community variation (the average trait value in communities where species occur; β traits). Trait diversity was mostly explained by variation among co-occurring species, not among-communities. Additionally, there was a phylogenetic signal in the within-community trait values of species relative to co-occurring taxa, but not in their habitat associations or among-community trait variation. These results suggest that sorting of pre-existing trait variation into local communities can explain the leaf and root trait diversity in these grasslands.