A multispecies coalescent model for quantitative traits
Figures
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig1-v2.tif/full/617,/0/default.jpg)
Three distinct discrete traits with their states mapped to the gene trees they evolved on, and to the species tree.
Hollow and filled shapes represent the ancestral and derived states, respectively, with character-state transitions being indicated by filled shapes along internal branches. Traits 1 and 2 undergo a single character-state transition in their evolutionary history, but when the states are resolved on the species tree, a homoplasy-like (yet not truly homoplastic) pattern emerges (i.e. hemiplasy). Trait three has evolved along a gene tree that matches the species tree in topology, and so no hemiplasy occurs.
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig2-v2.tif/full/617,/0/default.jpg)
Three- and five-taxon trees with branch lengths indicated.
(a) Three-species phylogeny (and a concordant gene tree within it) and its corresponding T matrix entries. (b) Five-species phylogeny used in coalescent simulations for PCM analyses. Branch lengths are indicated in units of 2N generations.
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig3-v2.tif/full/617,/0/default.jpg)
A single continuous traits controlled by five loci, four of which have discordant gene trees.
(a) Genealogies of the loci controlling the trait. Asterisks represent mutations at a given site in each of the five loci. Ancestral alleles (0) have no effect on the trait value. Derived alleles (1) have their random effects on the trait value drawn from a mutational effect distribution (see (b)). (b) Mutational effect distribution of derived alleles. The distribution has a mean of zero and unit variance. (c) The outcome of a simulation consists of one trait value per species, which correspond to the sum of all derived allele mutational effects coming from all loci controlling the trait.
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig4-v2.tif/full/617,/0/default.jpg)
Quantitative trait summaries under varying levels of genealogical discordance, for a three-taxon species tree.
(a) Expected and observed variances in trait values of species A in each of the five ILS conditions. Expected values come from Equation 2. (b) Expected and observed covariances between species A and B in each of the five ILS conditions. Expected values come from Equation 3.
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig5-v2.tif/full/617,/0/default.jpg)
Quantitative trait summaries under varying levels of genealogical discordance, indicating the effect of hemiplasy for a three-taxon species tree.
(a) Expected and observed covariances between species A and C in each of the five ILS conditions. Expected values come from Equation 4. (b) Observed covariances between a pair of species normalized by the variance in species A, for all five ILS conditions.
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig6-v2.tif/full/617,/0/default.jpg)
Mean parameter estimates under varying levels of genealogical discordance and number of loci controlling quantitative trait.
(a) Mean evolutionary rate for different number of loci controlling the simulated continuous trait and different levels of discordance. (b) Mean value of Pagel’s λ for different number of loci controlling the simulated continuous trait and different levels of discordance. (c) Mean evolutionary rate when 100 loci control the trait (‘Without λ’ is the same as shown in (a); in ‘With λ’, the rate was estimated with Pagel’s λ).
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig7-v2.tif/full/617,/0/default.jpg)
Mean number of independent trait-value shifts (i.e. significant phylogenetic ANOVA tests) among all possible groupings of two versus three species.
https://doi.org/10.7554/eLife.36482.008![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-fig8-v2.tif/full/617,/0/default.jpg)
Frequency of incongruent trait patterns (out of all informative trait patterns) for threshold traits.
Each combination of level of discordance and number of loci was simulated 1000 times.
![](https://iiif.elifesciences.org/lax/36482%2Felife-36482-app1-fig1-v2.tif/full/617,/0/default.jpg)
Divergence times, (, and for the species tree depicted in gray), and coalescence times, (, and for genealogy depicted in black).
https://doi.org/10.7554/eLife.36482.013Additional files
-
Source code file 1
R script for plotting figures.
- https://doi.org/10.7554/eLife.36482.010
-
Transparent reporting form
- https://doi.org/10.7554/eLife.36482.011