MOTMOT: models of trait macroevolution on trees
|
|
- Erica Harrison
- 5 years ago
- Views:
Transcription
1 Methods in Ecology and Evolution 2012, 3, doi: /j X x APPLICATION MOTMOT: models of trait macroevolution on trees Gavin H. Thomas 1 * and Robert P. Freckleton 2 1 School of Biological Sciences, University of Bristol, Woodland Road, Bristol BS8 1UG, UK; and 2 Department of Animal & Plant Sciences, University of Sheffield, Sheffield S10 2TN, UK Summary 1. Models of trait macroevolution on trees (MOTMOT) is a new software package that tests for variation in the tempo and mode of continuous character evolution on phylogenetic trees. MOT- MOT provides tools to fit a range of models of trait evolution with emphasis on variation in the rate of evolution between clades and character states. 2. We introduce a new method, trait MEDUSA, to identify the location of major changes in the rate of evolution of continuous traits on phylogenetic trees. We demonstrate trait MEDUSA and the other main functions of MOTMOT, using body size of Anolis lizards. 3. MOTMOT is open source software written in the R language and is freely available from CRAN ( Key-words: comparative analysis, evolutionary biology, macroevolution Introduction Phylogenetic trees describe the evolutionary relationships among taxa but also contain information on variation in the tempo and mode of evolution among lineages and through time. Identifying variation in rates of phenotypic evolution on phylogenies provides an important step in revealing the mechanisms that have generated differences in the temporal, geographic and taxonomic distribution of biological diversity over large spatial and temporal scales. We introduce models of trait macroevolution on trees (MOTMOT), a new package written in the R programming language (R Development Core Team 2010), that provides methods for assessing patterns in the tempo of continuous trait evolution on phylogenetic trees. MOTMOT uses transforms of phylogenetic branch lengths to test for non-brownian trait evolution and complements several other R packages notably GEIGER (Harmon et al. 2008), ape (Paradis, Claude, & Strimmer 2004) and caper (Orme et al. 2011). We demonstrate some of the main uses of MOTMOT using body size evolution in Caribbean Anolis. Brownian motion and models of trait evolution MOTMOT uses statistical models of trait evolution based on Brownian motion. Brownian motion has been described in the context of phylogenetic trait evolution numerous times (e.g. *Correspondence author. gavin.thomas@bristol.ac.uk Correspondence site: Felsenstein 1973; Freckleton, Harvey, & Pagel 2002; O Meara et al. 2006), so we keep the description here brief. The central tenet is that the magnitude and direction of changes in trait values are independent of the current state of the character and have an expected mean change of zero. Trait variance accumulates linearly with time, and the expected covariances of related species are proportional to the amount of shared evolutionary history (phylogenetic branch length). Hansen & Martins (1996) showed that the Brownian model adequately describes the expected covariances among related species under some microevolutionary models including random genetic drift, random punctuated change, directional selection and stabilising selection with a shifting optimum. Under such models, the Brownian variance (r 2 ) can be taken as a measure of the average rate of evolution (Garland 1992; McPeek 1995; O Meara et al. 2006; Thomas, Freckleton, & Sze kely 2006). The unbiased estimate of r 2 can be obtained using generalised least squares (GLS; e.g., Pagel 1997; Garland & Ives 2000; Freckleton, Harvey, & Pagel 2002), r 2 ¼ 1 ðn 1Þ ðy ^axþt V 1 ðy ^axþ eqn 1 where n is the number of tips, y is an n 1 vector of trait values at the tips, ^a is an n 1 vector of the phylogenetic mean for the trait (under Brownian motion, this is the estimated state at the root), V is the variance covariance matrix representation of the phylogenetic tree and X is an n 1 vector of ones in the constant rate Brownian model but can be treated as a design matrix in heterogeneous Ó 2011 The Authors. Methods in Ecology and Evolution Ó 2011 British Ecological Society
2 146 G. H. Thomas & R. P. Freckleton rate models where rate varies according to the state of a discrete trait (see Thomas, Freckleton, & Sze kely 2006; Thomas, Meiri, & Phillimore 2009). The fit of the Brownian model to trait data and a phylogeny can be assessed with the log-likelihood equation lnðlþ¼ 1 2 nlogð2pr2 ÞþlogjVjþ ðy ^axþt V 1 ðy ^axþ r 2! : eqn 2 Evaluation of the log likelihood with eqn 2 can be prohibitively slow for large data sets. However,fortheBrownian model, GLS approaches are exactly equivalent to phylogenetically independent contrasts (Felsenstein 1973; Garland & Ives 2000; Freckleton & Jetz 2009). Following Freckleton & Jetz (2009), eqn 2 then becomes lnðlþ ¼ 1! 2 n logð2pr2 Þþ Xn log V i þ u2 i V i r 2 eqn 3 i¼0 where u i is the unstandardised contrast at node i and V i is the corresponding expected variance following Felsenstein (1985). Phylogenetically independent contrasts allow rapid evaluation of the log likelihood, even for trees with thousands of tips. It is possible to extend the contrasts method to multivariate data sets, " #! L ¼ 1 2 nk log ð 2pÞþn log j W jþxn i¼0 k log V i þ ut W 1 u i i V i eqn 4 where W is the k k covariance matrix for the standardised contrasts with elements calculated as follows: w jk ¼ 1 X n u j;i u k;i ; eqn 5 n 1 V i¼0 i in which u j,i is the contrast for trait j at node i. The likelihood equations above assess the fit of the Brownian motion model to the data and phylogeny; however, we are generally more interested in identifying non-brownian evolution. The expected covariances among species can be altered by adjusting the phylogenetic branch lengths, and the fit of alternative models can be used to test for the different modes of evolution (Grafen 1989; Garland, Harvey, & Ives 1992; Pagel 1992, 1994, 1997, 1999; Mooers & Heard 1997; Blomberg, Garland, & Ives 2003; Thomas, Freckleton, & Székely 2006; Harmon et al. 2010; Ingram 2011). MOTMOT functionality MOTMOT fits a range of models of trait evolution (see Table 1) and has two main functions for parameter estimation. The first, transformphylo.ml, optimises the likelihood of uni- or multivariate data using the independent contrasts method (eqns 3 and 4). This allows rapid fitting of the models of trait evolution listed above (also Table 1) that assume a single mean value across all taxa. We use the L-BFGS-B routine (Byrd et al. 1995) in the optim function (stats library) for optimisation. This allows upper and lower bounds on parameters. The free model (Mooers, Vamosi, & Schluter 1999) is over-parameterised, so finding the optimum branch scaling is difficult and slow. The free model in MOTMOT should only be used to give a rough visualisation of rate variation on the phylogeny. The second maximum likelihood function, ML.RatePhylo, implements methods described by Thomas, Freckleton, & Székely (2006), Thomas, Meiri, & Phillimore (2009) in which rates of trait diversification are expected to vary according to the state of a discrete trait. There are several ways in which the rate of evolution may vary among species within clades or groups of taxa (Fig. 1). First, the relative rate of evolution differs because of trait divergence among species (disparification) of one group compared to one or more other groups (Fig. 1a and b). However, changes in rate may also occur in the absence of any net divergence among focal species, for example, if there is a major change in the direction of trait evolution in the root branch but no subsequent divergence among taxa in the clade (Fig. 1c compared with Fig. 1a). Alternatively, there may be continued directional change along all branches; this change is Table 1. Models of trait evolution available in models of trait macroevolution on trees (MOTMOT) Model Tests for MOTMOT Function Multivariate k 1 Phylogenetic signal transformphylo.ml Yes a 2 Constraints or strength of selection transformphylo.ml Yes j 1 Speciational vs. gradual evolution transformphylo.ml Yes w 3 Speciational vs. gradual evolution transformphylo.ml Yes d 1 Temporal change in rate transformphylo.ml Yes Hypothesised clade or Rate variation between lineages transformphylo.ml Yes trait-dependent rate shift (single mean) 4 6 ML.RatePhylo No Hypothesised clade or Rate variation between lineages ML.RatePhylo No trait-dependent rate shift (multiple means) 4 6 Clade rate shifts 4 7 Location and magnitude of rate shifts transformphylo.ml Yes Free 8 Independent branch rates transformphylo.ml Yes 1 Pagel (1997); 2 Hansen (1997); 3 Ingram (2011); 4 Thomas, Freckleton, & Székely (2006); Thomas, Meiri, & Phillimore (2009, 2009); 6 O Meara et al. (2006); 7 Alfaro et al. (2009); 8 Mooers, Vamosi, & Schluter (1999).
3 MOTMOT 147 in the same direction in each branch but with no net divergence among taxa compared with the null (Fig. 1d compared with Fig. 1a). Fig. 1c and d are characterised by a change in the mean trait value but no net change in divergence among taxa. The function ML.RatePhylo allows the user to test for different forms of rate shift by explicitly modelling differences in mean between groups. This cannot be readily achieved using independent contrasts because it also requires specification of a design matrix (see Thomas, Meiri, & Phillimore 2009). We therefore use the GLS log-likelihood function(eqn 2)and optimise using the L-BFGS-B routine optim in ML.RatePhylo. (a) (b) (c) (d) Trait divergence Fig. 1. Forms of rate variation of traits. The rate of evolution among species is faster in (b) than in (a) but the mean is identical. The rate of evolution of among species evolution in (c) and (d) is similar to (a) but the means differ. In (c), the mean differs because of directional change in the stem lineage, whereas in (d), the mean differs because of directional change in all branches. In addition to hypothesis testing, MOTMOT s transformphylo.sim function also allows the user to simulate trait evolution under all of the models described above. IDENTIFYING RATE SHIFTS While rate variation models are typically used as hypothesis tests, we frequently have no aprioriexpectation on the location or number of rate shifts. We take inspiration from a recent method [Modelling Evolutionary Diversification Using Stepwise AIC (MEDUSA)] for identifying the position and magnitude of shifts in the rate of lineage diversification (Rabosky 2006; Alfaro et al. 2009) and apply an analogous approach to rate shifts in continuous character diversification. The trait medusa algorithm 1 (tm1 algorithm) in the transformphylo.ml function proceeds as follows: (1) evaluate the likelihood of a single-rate Brownian model; (2) fit rate heterogeneous models and evaluate the likelihood at each node in the phylogeny (one at a time) where the second rate is permitted in all branches descended from the focal node plus the stem branch leading to that node (new shifts are inferred only as monophyletic groups); (3) select the best-fitting two-rate model; (4) fit rate heterogeneous models with two-rate shifts where one of the shifts mustoccuratthenodeidentifiedinstep3(eachadditional shift defines a monophyletic set of taxa but rate shifts can be nested); (5) continue this procedure until a user-defined maximum number of rate shifts has been reached; (6) select the preferred model from the set of best-rate heterogeneous models using the Akaike Information Criterion (AIC). The number of parameters is equal to the total number of rates estimated plus the number of rate shifts inferred. For example, a model with two-rate shifts has three rates and hence five parameters in total. In the interests of computational efficiency, the tm1 algorithm uses independent contrasts. As a result, the tm1 algorithm searches for clade-wide rate shifts but cannot model differences in means between groups in the same way as ML.RatePhylo. However, we only infer clade (rather than polyphyletic group) shifts, so the possibility of directional trait change rather than net trait divergence can be tested by allowing rate shifts on the stem branch leading to a clade. The clade option in transformphylo.ml permits post hoc tests of single-branch vs. clade-wide rate shifts. Alternatively, the tm2 algorithm provides a more thorough model fitting procedure. The tm2 algorithm proceeds in the same way as tm1; however, at each node, it evaluates the fit of the standard clade-wide shift (as described for tm1 above) and a single-branch shift where the rate shift occurs only on the stem branch leading to the focal node. At each step, the best model is chosen by comparison of all clade-wide and single-branch models. In subsequent steps (three-rate, four-rate and so on), it is possible to revisit, for example, a node with a previously inferred clade-wide rate shift and test additionally for a single-branch shift. Consequently, the final model can contain any combination of single-branch and clade-wide shifts including both shift types at thesamenodes.
4 148 G. H. Thomas & R. P. Freckleton (a) delta AICc (b) delta AICc tm1 tm Minimum clade size tm1 tm2 For large trees, the number of possible rate shifts is high. MOTMOT allows the user to restrict search space by (1) setting an upper limit on the number of rate shifts to infer; (2) specifying a minimum clade species richness such that the rate shift search does not proceed within smaller clades; and (3) listing specific clades within which no further rate shifts should be inferred. The final preferred multiple rate model is determined by comparison of the Akaike Information Criterion (AIC, or its small sample alternative, AICc; Akaike 1973; Burnham & Anderson 2002). The choice of a suitable cut-off value for DAIC is not trivial. This is because the number of potential nodes at which rate shifts could occur will increase with increasing tree size and if the user decreases the minimum clade size. Consequently, the probability of falsely inferring at least one rate shift is likely to be high and to scale with the number of available nodes in the tree. We conducted simulations to assess a suitable DAIC and the accuracy of shift localities. Full details of the simulations are provided as Supporting Information. In brief, a suitable cut-off for DAICc (i.e. large enough to reject the null at a =0Æ05) declines as the number of possible shift locations declines (Fig. 2a). Even with a high DAICc cut-off, the power to detect real rate shifts remains reasonably high (60 70%; Fig. 2b) and the location of shifts is generally inferred accurately (Fig. 2c), although the latter is affected by minimum clade size (see Supporting Information). (c) Frequency Minimum clade size Proportion nested = 0 82 Proportion non nested = 0 18 Body size evolution in Greater Antillean Anolis The radiation of Anolis lizards in the Caribbean has become a classic model system for studying adaptive radiations and macroevolutionary trends with particular interest in the repeated evolution of ecomorphs on different islands (Losos 2009). We use a recently published data set of maximum adult snout-vent length (hereafter body size) measurements for male Anolis (Thomas, Meiri, & Phillimore 2009) and a molecular phylogeny (pruned from the larger tree in Nicholson et al. 2005) to demonstrate some of the main functionality of MOTMOT. The data, phylogeny and annotated R scripts are provided as online Supporting Information Node distance Fig. 2. Simulation results for trait Modelling Evolutionary Diversification Using Stepwise AIC. (a) Estimated DAIC values at a =0.05 for tm1 and tm2 algorithms with minimum clade sizes of 1, 2, 5, 10 and 20 based on data simulated under a constant rate Brownian model. (b) Violin plots showing the frequency distribution of estimated DAIC for data with a single-rate shift five times faster than the background rate. The area below the red line indicates the proportion of simulated data sets where the constant rate model could not be rejected. (c) Frequency distribution showing the distance (in number of nodes) between the inferred rate shift node and the actual rate shift node for the tm1 algorithm with a minimum clade size of five. The figure shows data for the 82% of simulated data sets where the inferred node nested (positive values) or was nested within (negative values) the true-rate shift node. Table 2. Parameter estimates and model fit for Anolis male body size Model ML parameter estimate (approx. conf. int.) Maximum likelihood bm 83Æ391 k 0Æ991 (0Æ953 1) 84Æ453 2Æ125 a 0Æ000 (0Æ000 0Æ005) 83Æ391 0Æ000 j 1Æ040 (0Æ775 1Æ256) 83Æ443 0Æ105 w 0Æ019 (0Æ000 0Æ103) 84Æ349 1Æ916 d 0Æ696 (0Æ145 1Æ426) 83Æ772 0Æ763 LRT vs. bm Model comparison with the Brownian (bm) null is by likelihood ratio test (LRT) with one degree of freedom. All LRTs were nonsignificant.
5 MOTMOT 149 MODELS OF TRAIT EVOLUTION We compared the fit of the Brownian model, transformphylo.ml (data, phy=phylogeny, model='bm'),tok (model='lambda'), a (model='ou'), j (model='kappa'), w (model='psi') and d (model='delta') for log body size. Although all models, with the exception of OU, fit the data better than the Brownian model, the improvement in likelihood is not sufficient to reject Brownian motion (Table 2). (a) (b) Relative rates Node tm1 tm Fig. 3. Body size rates shifts in Caribbean Anolis. (a) Rate shifts inferred from the tm1 algorithm. (b) Rate shifts inferred from the tm2 algorithm and tm1 with post hoc single-branch testing. Branch lengths are scaled according to the maximum likelihood relative rate estimates (inset table) and coloured according to low (blue), high (red) or background (black) rates.
6 150 G. H. Thomas & R. P. Freckleton RATE VARIATION We tested for rate variation in male body size first in relation to geography and second using the tm1 and tm2 algorithms. For the geographic analyses, we compared rates between species from the four main islands of the Greater Antilles (Cuba, Hispaniola, Jamaica, Puerto Rico) using the ML.RatePhylo function. We used the ancestral state reconstruction of island occupancy in Mahler et al. (2010) to assign all internal branches to one of the four islands where each internal node has a single fixed value (unlike Mahler et al. who allowed uncertainty in ancestral states in their models). We fitted a single-rate model while assuming a single mean across all taxa, ML.RatePhylo(rateData, fixed=c(1,false, FALSE, FALSE), pretty=true, common.mean=true), and found no evidence for rate heterogeneity (log likelihood = 84Æ650, LRT vs. Brownian = 2Æ518, P = 0Æ472, d.f. = 3). We then fitted a single-rate model while allowing means to differ between islands: ML.RatePhylo(rateData, fixed=c(1,false, FALSE, FALSE), pretty=true, common.mean= FALSE). We again found no evidence for rate heterogeneity across species from the four islands (log likelihood = 86Æ997 LRT vs. Brownian with multiple means = 2Æ209, P = 0Æ530, d.f. = 3). For the tm1 and tm2 algorithms, we set the maximum number of possible rate shifts to 10 but did not impose a minimum clade size: transformphylo.ml(data, phy=phylogeny, model='tm1', mincladesize=1, nsplits=10). Thebest-rateshiftmodel(basedonanAICccut-offofnineas determined by the simulations described above) identified two-rate shifts (Fig. 3a) for the tm1 algorithm and three-rate shifts for the tm2 algorithm (Fig. 3b). The first two-rate shifts using the tm2 algorithm occur at the same nodes as the two shifts identified using tm1. The tm2 algorithm additionally identified a substantial rate shift on a single branch. We note that if we use a more liberal AICc cut-off of four, the tm1 method also identifies a rate shift at this node but is limited to inference of a clade-wide rate shift (results not shown). We also note that the first rate shift using tm1 and tm2 (node 1 in Fig. 3) pertains to two taxa with identical trait values. It seems unlikely that the species genuinely have identical body sizes, and this may reflect crude measurement and or rounding error. Conclusions R is rapidly becoming the default platform for phylogenetic comparative analyses. MOTMOT complements the growing number of phylogenetic R packages. Here, we have demonstrated the main functionality of MOTMOT using Anolis of the Caribbean as an example. This functionality can readily be extended to incorporate any trait evolution model based on transformations to the branch lengths of the phylogeny. Acknowledgements We thank David Orme, Luke Harmon and Emanuel Paradis whose R packages (caper, GEIGER and ape) underpin much of the code in MOTMOT; Luke Harmon, Matt Pennell, David Orme and two anonymous reviewers for their insightful reviews; and Travis Ingram, Albert Phillimore, Annat Haber and Chris Vendetti for additional helpful discussion and comments on the methods, code and manuscript. This work was funded by a Natural Environment Research Council (NERC, UK) Postdoctoral Research Fellowship (grant number NE G ) and Royal Society Grant (grant number RG ) to GHT. References Akaike, H. (1973) Information theory as an extension of the maximum likelihood principle. Second International Symposium on Information Theory (eds B.N. Petrov & F. Csaki), pp Akademiai Kiado, Budapest. Alfaro, M.E., Santini, F., Brock, C., Alamillo, H., Dornburg, A., Rabosky, D.L., Carnevale, G. & Harmon, L.J. (2009) Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. Proceedings of the National Academy of Sciences of the United States of America, 106, Blomberg, S.P., Garland, T.J. & Ives, A.R. (2003) Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution, 57, Burnham, K.P. & Anderson, D.R. (2002) Model Selection and Mulimodel Inference: A Practical Information-theoretic Approach. Springer-Verlag, New York. Byrd, R.H., Lu, P., Nocedal, J. & Zhu, C. (1995) A limited memory algorithm for bound constrained optimization. SIAM Journal of Scientific Computing, 16, Felsenstein, J. (1973) Maximum-likelihood estimation of evolutionary trees from continuous characters. American Journal of Human Genetics, 25, Felsenstein, J. (1985) Phylogenies and the comparative method. The American Naturalist, 125, Freckleton, R.P., Harvey, P.H. & Pagel, M. (2002) Phylogenetic analysis and comparative data: a test and review of the evidence. The American Naturalist, 125,1 15. Freckleton, R.P. & Jetz, W. (2009) Space versus phylogeny: disentangling phylogenetic and spatial signals in comparative data. Proceedings of the Royal Society of London Series B, 276, Garland, T. (1992) Rate tests for phenotypic evolution using phylogenetically independent contrasts. The American Naturalist, 140, Garland, T.J., Harvey, P.H. & Ives, A.R. (1992) Procedures for the analysis of comparative data using phylogenetically independent contrasts. Systematic Biology, 41, Garland, T.J. & Ives, A.R. (2000) Using the past to predict the present: Confidence intervals for regression equations in phylogenetic comparative methods. The American Naturalist, 155, Grafen, A. (1989) The phylogenetic regression. Philosophical Transactions of the Royal Society B, 326, Hansen, T.F. (1997) Stabilizing selection and the comparative analysis of adaptation. Evolution, 51, Hansen, T.F. & Martins, E.P. (1996) Translating between microevolutionary process and macroevolutionary patterns: the correlation structure of interspecific data. Evolution, 50, Harmon, L.J., Weir, J., Brock, C., Glor, R.E. & Challenger, W. (2008) GEI- GER: investigating evolutionary radiations. Bioinformatics, 24, Harmon, L.J., Losos, J.B., Davies, T.J., Gillespie, R.G., Gittleman, J.L., Jennings, W.B., Kozak, K.H., McPeek, M.A., Moreno-Roark, F., Near, T.J., Purvis, A., Ricklefs, R.E., Schluter, D., Schulte II, J.A., Seehausen, O., Sidlauskas, B.L., Torres-Carvajal, O., Weir, J.T. & Mooers, A.Ø. (2010) Early bursts of body size and shape evolution are rare in comparative data. Evolution, 64, Ingram, T. (2011) Speciation along a depth gradient in a marine adaptive radiation. Proceedings of the Royal Society of London Series B, 278, Losos, J.B. (2009) Lizards in an Evolutionary Tree: Ecology and Adaptive Radiation of Anoles. University of California Press, Berkeley. Mahler, D.L., Revell, L.J., Glor, R.E. & Losos, J.B. (2010) Ecological opportunity and the rate of morphological evolution in the diversification of Greater antillean anoles. Evolution, 64,
7 MOTMOT 151 McPeek, M.A. (1995) Testing hypotheses about evolutionary change on single branches of phylogeny using evolutionary contrasts. The American Naturalist, 145, Mooers, A.O. & Heard, S.B. (1997) Inferring evolutionary process from phylogenetic tree shape. Quarterly Review of Biology, 72, Mooers, A.O., Vamosi, S.M. & Schluter, D. (1999) Using phylogenies to test macroevolutionary hypotheses of trait evolution in cranes (Gruidae). The American Naturalist, 154, Nicholson, K.E., Glor, R.E., Kolbe, J.J., Larson, A., Blair Hedges, S. & Losos, J. (2005) Mainland colonization by island lizards. Journal of Biogeography, 32, O Meara, B.C., Ané, C., Sanderson, M.J. & Wainwright, P.C. (2006) Testing for different rates of continuous trait evolution using likelihood. Evolution, 60, Orme, C.D.L., Freckleton, R.P., Thomas, G.H., Petzoldt, T. & Fritz, S.A. (2011) caper: Comparative Analyses of Phylogenetics and Evolution in R ( Pagel, M. (1992) A method for the analysis of comparative data. Journal of Theoretical Biology, 156, Pagel, M. (1994) Detecting correlated evolution on phylogenies: a general method for the comparative analysis of discrete characters. Proceedings of the Royal Society of London Series B, 255, Pagel, M. (1997) Inferring evolutionary processes from phylogenies. Zoologica Scripta, 26, Pagel, M. (1999) Inferring the historical patterns of biological evolution. Nature, 401, Paradis, E., Claude, J. & Strimmer, K. (2004) APE: analyses of Phylogenetics and Evolution in R language. Bioinformatics, 20, R Development Core Team (2010) R: A Language and Environment for Statistical Computing, Vienna, Austria, ISBN , Rabosky, D.L. (2006) Likelihood methods for detecting temporal rate shifts in diversification rates. Evolution, 60, Thomas, G.H., Freckleton, R.P. & Sze kely, T. (2006) Comparative analyses of the influence of developmental mode on phenotypic diversification rates in shorebirds. Proceedings of the Royal Society of London Series B, 273, Thomas, G.H., Meiri, S. & Phillimore, A.B. (2009) Body size diversification in Anolis: novel environments and island effects. Evolution, 63, Received 6 January 2011; accepted 13 May 2011 Handling Editor: David Orme Supporting Information Additional Supporting Information may be found in the online version of this article. Data S1. Full simulation results for trait MEDUSA and comparison of run times to fit models of trait evolution in MOTMOT and GEI- GER. As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials may be re-organised for online delivery, but are not copy-edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.
Phylogenetic comparative methods, CEB Spring 2010 [Revised April 22, 2010]
Phylogenetic comparative methods, CEB 35300 Spring 2010 [Revised April 22, 2010] Instructors: Andrew Hipp, Herbarium Curator, Morton Arboretum. ahipp@mortonarb.org; 630-725-2094 Rick Ree, Associate Curator,
More informationPOSITIVE CORRELATION BETWEEN DIVERSIFICATION RATES AND PHENOTYPIC EVOLVABILITY CAN MIMIC PUNCTUATED EQUILIBRIUM ON MOLECULAR PHYLOGENIES
doi:10.1111/j.1558-5646.2012.01631.x POSITIVE CORRELATION BETWEEN DIVERSIFICATION RATES AND PHENOTYPIC EVOLVABILITY CAN MIMIC PUNCTUATED EQUILIBRIUM ON MOLECULAR PHYLOGENIES Daniel L. Rabosky 1,2,3 1 Department
More informationAPPENDIX S1: DESCRIPTION OF THE ESTIMATION OF THE VARIANCES OF OUR MAXIMUM LIKELIHOOD ESTIMATORS
APPENDIX S1: DESCRIPTION OF THE ESTIMATION OF THE VARIANCES OF OUR MAXIMUM LIKELIHOOD ESTIMATORS For large n, the negative inverse of the second derivative of the likelihood function at the optimum can
More informationModels of Continuous Trait Evolution
Models of Continuous Trait Evolution By Ricardo Betancur (slides courtesy of Graham Slater, J. Arbour, L. Harmon & M. Alfaro) How fast do animals evolve? C. Darwin G. Simpson J. Haldane How do we model
More informationEvolutionary models with ecological interactions. By: Magnus Clarke
Evolutionary models with ecological interactions By: Magnus Clarke A thesis submitted in partial fulfilment of the requirements for the degree of Doctor of Philosophy The University of Sheffield Faculty
More informationBAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees
Methods in Ecology and Evolution 2, 5, 71 77 doi: 1.1111/241-21X.12199 APPLICATION BAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees Daniel L. Rabosky*, Michael Grundler,
More informationChapter 6: Beyond Brownian Motion Section 6.1: Introduction
Chapter 6: Beyond Brownian Motion Section 6.1: Introduction Detailed studies of contemporary evolution have revealed a rich variety of processes that influence how traits evolve through time. Consider
More informationBRIEF COMMUNICATIONS
Evolution, 59(12), 2005, pp. 2705 2710 BRIEF COMMUNICATIONS THE EFFECT OF INTRASPECIFIC SAMPLE SIZE ON TYPE I AND TYPE II ERROR RATES IN COMPARATIVE STUDIES LUKE J. HARMON 1 AND JONATHAN B. LOSOS Department
More informationpaleotree: an R package for paleontological and phylogenetic analyses of evolution
Methods in Ecology and Evolution 2012, 3, 803 807 doi: 10.1111/j.2041-210X.2012.00223.x APPLICATION paleotree: an R package for paleontological and phylogenetic analyses of evolution David W. Bapst* Department
More informationRapid evolution of the cerebellum in humans and other great apes
Rapid evolution of the cerebellum in humans and other great apes Article Accepted Version Barton, R. A. and Venditti, C. (2014) Rapid evolution of the cerebellum in humans and other great apes. Current
More informationBiol 206/306 Advanced Biostatistics Lab 11 Models of Trait Evolution Fall 2016
Biol 206/306 Advanced Biostatistics Lab 11 Models of Trait Evolution Fall 2016 By Philip J. Bergmann 0. Laboratory Objectives 1. Explore how evolutionary trait modeling can reveal different information
More informationLaboratoire d Écologie Alpine (LECA), University of Grenoble Alpes, F-38000, Grenoble, France 2
Reports Ecology, 97(2), 2016, pp. 286 293 2016 by the Ecological Society of America Improving phylogenetic regression under complex evolutionary models Florent Mazel, 1,2,6 T. Jonathan Davies, 3,4 Damien
More informationPhylogenetic Comparative Analysis: A Modeling Approach for Adaptive Evolution
vol. 164, no. 6 the american naturalist december 2004 Phylogenetic Comparative Analysis: A Modeling Approach for Adaptive Evolution Marguerite A. Butler * and Aaron A. King Department of Ecology and Evolutionary
More informationTitle: Area, climate heterogeneity, and the response of climate niches to ecological opportunity in island radiations of Anolis lizards.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 Title: Area, climate heterogeneity, and the response of climate niches to ecological opportunity in island radiations
More informationDNA-based species delimitation
DNA-based species delimitation Phylogenetic species concept based on tree topologies Ø How to set species boundaries? Ø Automatic species delimitation? druhů? DNA barcoding Species boundaries recognized
More information20 Grenoble, France ; CNRS, Laboratoire d Écologie Alpine (LECA), F Grenoble, France
1/23 1 Improving phylogenetic regression under complex evolutionary models 2 Florent Mazel, T. Jonathan Davies, Damien Georges, Sébastien Lavergne, Wilfried Thuiller and 3 Pedro R. Peres-Neto 4 Running
More informationASYMMETRIES IN PHYLOGENETIC DIVERSIFICATION AND CHARACTER CHANGE CAN BE UNTANGLED
doi:1.1111/j.1558-5646.27.252.x ASYMMETRIES IN PHYLOGENETIC DIVERSIFICATION AND CHARACTER CHANGE CAN BE UNTANGLED Emmanuel Paradis Institut de Recherche pour le Développement, UR175 CAVIAR, GAMET BP 595,
More informationThe Tempo of Macroevolution: Patterns of Diversification and Extinction
The Tempo of Macroevolution: Patterns of Diversification and Extinction During the semester we have been consider various aspects parameters associated with biodiversity. Current usage stems from 1980's
More informationDoes convergent evolution always result from different lineages experiencing similar
Different evolutionary dynamics led to the convergence of clinging performance in lizard toepads Sarin Tiatragula 1 *, Gopal Murali 2 *, James T. Stroud 3 1 Department of Biological Sciences, Auburn University,
More information"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200B Spring 2011 University of California, Berkeley
"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200B Spring 2011 University of California, Berkeley B.D. Mishler Feb. 1, 2011. Qualitative character evolution (cont.) - comparing
More informationIntegrative Biology 200 "PRINCIPLES OF PHYLOGENETICS" Spring 2018 University of California, Berkeley
Integrative Biology 200 "PRINCIPLES OF PHYLOGENETICS" Spring 2018 University of California, Berkeley B.D. Mishler Feb. 14, 2018. Phylogenetic trees VI: Dating in the 21st century: clocks, & calibrations;
More informationBiology 559R: Introduction to Phylogenetic Comparative Methods Topics for this week:
Biology 559R: Introduction to Phylogenetic Comparative Methods Topics for this week: Course general information About the course Course objectives Comparative methods: An overview R as language: uses and
More informationAppendix from L. J. Revell, On the Analysis of Evolutionary Change along Single Branches in a Phylogeny
008 by The University of Chicago. All rights reserved.doi: 10.1086/588078 Appendix from L. J. Revell, On the Analysis of Evolutionary Change along Single Branches in a Phylogeny (Am. Nat., vol. 17, no.
More informationGreater host breadth still not associated with increased diversification rate in the Nymphalidae A response to Janz et al.
doi:10.1111/evo.12914 Greater host breadth still not associated with increased diversification rate in the Nymphalidae A response to Janz et al. Christopher A. Hamm 1,2 and James A. Fordyce 3 1 Department
More informationTesting quantitative genetic hypotheses about the evolutionary rate matrix for continuous characters
Evolutionary Ecology Research, 2008, 10: 311 331 Testing quantitative genetic hypotheses about the evolutionary rate matrix for continuous characters Liam J. Revell 1 * and Luke J. Harmon 2 1 Department
More informationPackage OUwie. August 29, 2013
Package OUwie August 29, 2013 Version 1.34 Date 2013-5-21 Title Analysis of evolutionary rates in an OU framework Author Jeremy M. Beaulieu , Brian O Meara Maintainer
More informationEXTINCTION RATES SHOULD NOT BE ESTIMATED FROM MOLECULAR PHYLOGENIES
ORIGINAL ARTICLE doi:10.1111/j.1558-5646.2009.00926.x EXTINCTION RATES SHOULD NOT BE ESTIMATED FROM MOLECULAR PHYLOGENIES Daniel L. Rabosky 1,2,3,4 1 Department of Ecology and Evolutionary Biology, Cornell
More information"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200 Spring 2018 University of California, Berkeley
"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200 Spring 2018 University of California, Berkeley D.D. Ackerly Feb. 26, 2018 Maximum Likelihood Principles, and Applications to
More informationPhylogenetic Signal, Evolutionary Process, and Rate
Syst. Biol. 57(4):591 601, 2008 Copyright c Society of Systematic Biologists ISSN: 1063-5157 print / 1076-836X online DOI: 10.1080/10635150802302427 Phylogenetic Signal, Evolutionary Process, and Rate
More informationMeasurement errors should always be incorporated in phylogenetic comparative analysis
Methods in Ecology and Evolution 25, 6, 34 346 doi:./24-2x.2337 s should always be incorporated in phylogenetic comparative analysis Daniele Silvestro,2,3, Anna Kostikova 2,3, Glenn Litsios 2,3, Peter
More informationPRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION
"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200 Spring 2018 University of California, Berkeley D. Ackerly April 2, 2018 Lineage diversification Reading: Morlon, H. (2014).
More informationEcological Limits on Clade Diversification in Higher Taxa
vol. 173, no. 5 the american naturalist may 2009 Ecological Limits on Clade Diversification in Higher Taxa Daniel L. Rabosky * Department of Ecology and Evolutionary Biology, Cornell University, Ithaca,
More informationPhylogenetic approaches for studying diversification
Ecology Letters, (2014) doi: 10.1111/ele.12251 REVIEW AND SYNTHESIS Phylogenetic approaches for studying diversification Helene Morlon* Center for Applied Mathematics, Ecole Polytechnique, Palaiseau, Essonne,
More informationNonlinear relationships and phylogenetically independent contrasts
doi:10.1111/j.1420-9101.2004.00697.x SHORT COMMUNICATION Nonlinear relationships and phylogenetically independent contrasts S. QUADER,* K. ISVARAN,* R. E. HALE,* B. G. MINER* & N. E. SEAVY* *Department
More informationFrank T. Burbrink 1,2 *, Sara Ruane 1,2 and R. Alexander Pyron 3
Journal of Biogeography (J. Biogeogr.) (2012) 39, 465 475 ORIGINAL ARTICLE When are adaptive radiations replicated in areas? Ecological opportunity and unexceptional diversification in West Indian dipsadine
More informationPhenotypic Evolution. and phylogenetic comparative methods. G562 Geometric Morphometrics. Department of Geological Sciences Indiana University
Phenotypic Evolution and phylogenetic comparative methods Phenotypic Evolution Change in the mean phenotype from generation to generation... Evolution = Mean(genetic variation * selection) + Mean(genetic
More informationEvolutionary Physiology
Evolutionary Physiology Yves Desdevises Observatoire Océanologique de Banyuls 04 68 88 73 13 desdevises@obs-banyuls.fr http://desdevises.free.fr 1 2 Physiology Physiology: how organisms work Lot of diversity
More informationUser s Manual for. Continuous. (copyright M. Pagel) Mark Pagel School of Animal and Microbial Sciences University of Reading Reading RG6 6AJ UK
User s Manual for Continuous (copyright M. Pagel) Mark Pagel School of Animal and Microbial Sciences University of Reading Reading RG6 6AJ UK email: m.pagel@rdg.ac.uk (www.ams.reading.ac.uk/zoology/pagel/)
More informationQuantitative evolution of morphology
Quantitative evolution of morphology Properties of Brownian motion evolution of a single quantitative trait Most likely outcome = starting value Variance of the outcomes = number of step * (rate parameter)
More informationIntrinsic inference difficulties for trait evolution with Ornstein-Uhlenbeck models
Methods in Ecology and Evolution 204, 5, 33 46 doi: 0./204-20X.2285 Intrinsic inference difficulties for trait evolution with Ornstein-Uhlenbeck models Lam Si Tung Ho, 2* and Cecile Ane Department of Statistics,
More informationHow can molecular phylogenies illuminate morphological evolution?
How can molecular phylogenies illuminate morphological evolution? 27 October 2016. Joe Felsenstein UNAM How can molecular phylogenies illuminate morphological evolution? p.1/81 Where this lecture fits
More informationDiscrete & continuous characters: The threshold model
Discrete & continuous characters: The threshold model Discrete & continuous characters: the threshold model So far we have discussed continuous & discrete character models separately for estimating ancestral
More informationCONVERGENT EVOLUTION OF PHENOTYPIC INTEGRATION AND ITS ALIGNMENT WITH MORPHOLOGICAL DIVERSIFICATION IN CARIBBEAN ANOLIS ECOMORPHS
ORIGINAL ARTICLE doi:10.1111/j.1558-5646.2011.01416.x CONVERGENT EVOLUTION OF PHENOTYPIC INTEGRATION AND ITS ALIGNMENT WITH MORPHOLOGICAL DIVERSIFICATION IN CARIBBEAN ANOLIS ECOMORPHS Jason J. Kolbe, 1,2
More informationR eports. Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism
Ecology, 93(2), 2012, pp. 242 247 Ó 2012 by the Ecological Society of America Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism T. JONATHAN DAVIES, 1,6 NATHAN J. B.
More informationComparing the rates of speciation and extinction between phylogenetic trees
Received: 29 November 2017 Revised: 21 February 2018 Accepted: 8 March 2018 DOI: 10.1002/ece3.4030 ORIGINAL RESEARCH Comparing the rates of speciation and extinction between phylogenetic trees Liam J.
More informationSPECIATION. REPRODUCTIVE BARRIERS PREZYGOTIC: Barriers that prevent fertilization. Habitat isolation Populations can t get together
SPECIATION Origin of new species=speciation -Process by which one species splits into two or more species, accounts for both the unity and diversity of life SPECIES BIOLOGICAL CONCEPT Population or groups
More informationUoN, CAS, DBSC BIOL102 lecture notes by: Dr. Mustafa A. Mansi. The Phylogenetic Systematics (Phylogeny and Systematics)
- Phylogeny? - Systematics? The Phylogenetic Systematics (Phylogeny and Systematics) - Phylogenetic systematics? Connection between phylogeny and classification. - Phylogenetic systematics informs the
More informationFiSSE: A simple nonparametric test for the effects of a binary character on lineage diversification rates
ORIGINAL ARTICLE doi:0./evo.3227 FiSSE: A simple nonparametric test for the effects of a binary character on lineage diversification rates Daniel L. Rabosky,2, and Emma E. Goldberg 3, Department of Ecology
More informationAP Biology. Cladistics
Cladistics Kingdom Summary Review slide Review slide Classification Old 5 Kingdom system Eukaryote Monera, Protists, Plants, Fungi, Animals New 3 Domain system reflects a greater understanding of evolution
More informationThe phylogenetic effective sample size and jumps
arxiv:1809.06672v1 [q-bio.pe] 18 Sep 2018 The phylogenetic effective sample size and jumps Krzysztof Bartoszek September 19, 2018 Abstract The phylogenetic effective sample size is a parameter that has
More informationTesting the time-for-speciation effect in the assembly of regional biotas
Methods in Ecology and Evolution 2012, 3, 224 233 Testing the time-for-speciation effect in the assembly of regional biotas Daniel L. Rabosky doi: 10.1111/j.2041-210X.2011.00166.x Department of Integrative
More informationModeling Phylogenetic Comparative Methods with Hybridization
Modeling Phylogenetic Comparative Methods with Hybridization Tony Jhwueng NIMBioS Interdisciplinary Seminar Jan 25 2011 Outline: 1 Introduction: Phylognetic Comparative Methods (PCMs). 2 Modeling PCMs
More informationReconstructing the history of lineages
Reconstructing the history of lineages Class outline Systematics Phylogenetic systematics Phylogenetic trees and maps Class outline Definitions Systematics Phylogenetic systematics/cladistics Systematics
More informationChapter 7: Models of discrete character evolution
Chapter 7: Models of discrete character evolution pdf version R markdown to recreate analyses Biological motivation: Limblessness as a discrete trait Squamates, the clade that includes all living species
More informationSome of these slides have been borrowed from Dr. Paul Lewis, Dr. Joe Felsenstein. Thanks!
Some of these slides have been borrowed from Dr. Paul Lewis, Dr. Joe Felsenstein. Thanks! Paul has many great tools for teaching phylogenetics at his web site: http://hydrodictyon.eeb.uconn.edu/people/plewis
More informationAmira A. AL-Hosary PhD of infectious diseases Department of Animal Medicine (Infectious Diseases) Faculty of Veterinary Medicine Assiut
Amira A. AL-Hosary PhD of infectious diseases Department of Animal Medicine (Infectious Diseases) Faculty of Veterinary Medicine Assiut University-Egypt Phylogenetic analysis Phylogenetic Basics: Biological
More informationQuantitative Traits and Diversification
Syst. Biol. 59(6):619 633, 2010 c The Author(s) 2010. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oxfordjournals.org
More informationStatistical nonmolecular phylogenetics: can molecular phylogenies illuminate morphological evolution?
Statistical nonmolecular phylogenetics: can molecular phylogenies illuminate morphological evolution? 30 July 2011. Joe Felsenstein Workshop on Molecular Evolution, MBL, Woods Hole Statistical nonmolecular
More informationQuantifying and Comparing Phylogenetic Evolutionary Rates for Shape and Other High-Dimensional Phenotypic Data
Syst. Biol. 63(2):166 177, 2014 The Author(s) 2013. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oup.com
More informationAn integrative view of phylogenetic comparative methods: connections to population genetics, community ecology, and paleobiology
Ann. N.Y. Acad. Sci. ISSN 0077-8923 ANNALS OF THE NEW YORK ACADEMY OF SCIENCES Issue: The Year in Evolutionary Biology An integrative view of phylogenetic comparative methods: connections to population
More informationContrasts for a within-species comparative method
Contrasts for a within-species comparative method Joseph Felsenstein, Department of Genetics, University of Washington, Box 357360, Seattle, Washington 98195-7360, USA email address: joe@genetics.washington.edu
More informationECOLOGICAL CAUSES OF DECELERATING DIVERSIFICATION IN CARNIVORAN MAMMALS
ORIGINAL ARTICLE doi:10.1111/evo.12126 ECOLOGICAL CAUSES OF DECELERATING DIVERSIFICATION IN CARNIVORAN MAMMALS Antonin Machac, 1,2,3,4 David Storch, 2,3 and John J. Wiens 1,5 1 Department of Ecology and
More informationEXPLOSIVE EVOLUTIONARY RADIATIONS: DECREASING SPECIATION OR INCREASING EXTINCTION THROUGH TIME?
ORIGINAL ARTICLE doi:10.1111/j.1558-5646.2008.00409.x EXPLOSIVE EVOLUTIONARY RADIATIONS: DECREASING SPECIATION OR INCREASING EXTINCTION THROUGH TIME? Daniel L. Rabosky 1,2,3 and Irby J. Lovette 1,3 1 Department
More informationModel Inadequacy and Mistaken Inferences of Trait-Dependent Speciation
Syst. Biol. 64(2):340 355, 2015 The Author(s) 2015. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oup.com
More informationNon-independence in Statistical Tests for Discrete Cross-species Data
J. theor. Biol. (1997) 188, 507514 Non-independence in Statistical Tests for Discrete Cross-species Data ALAN GRAFEN* AND MARK RIDLEY * St. John s College, Oxford OX1 3JP, and the Department of Zoology,
More informationOptimum selection and OU processes
Optimum selection and OU processes 29 November 2016. Joe Felsenstein Biology 550D Optimum selection and OU processes p.1/15 With selection... life is harder There is the Breeder s Equation of Wright and
More informationPackage WGDgc. October 23, 2015
Package WGDgc October 23, 2015 Version 1.2 Date 2015-10-22 Title Whole genome duplication detection using gene counts Depends R (>= 3.0.1), phylobase, phyext, ape Detection of whole genome duplications
More informationGENETICS - CLUTCH CH.22 EVOLUTIONARY GENETICS.
!! www.clutchprep.com CONCEPT: OVERVIEW OF EVOLUTION Evolution is a process through which variation in individuals makes it more likely for them to survive and reproduce There are principles to the theory
More informationEvidence for adaptive radiation from a phylogenetic study of plant defenses
Evidence for adaptive radiation from a phylogenetic study of plant defenses Anurag A. Agrawal a,b,1, Mark Fishbein c,2, Rayko Halitschke a, Amy P. Hastings a, Daniel L. Rabosky a, and Sergio Rasmann a
More informationA NEW BAYESIAN METHOD FOR FITTING EVOLUTIONARY MODELS TO COMPARATIVE DATA WITH INTRASPECIFIC VARIATION
ORIGINAL ARTICLE doi:10.1111/j.1558-5646.01.01645.x A NEW BAYESIAN METHOD FOR FITTING EVOLUTIONARY MODELS TO COMPARATIVE DATA WITH INTRASPECIFIC VARIATION Liam J. Revell 1, and R. Graham Reynolds 1 1 Department
More informationHow to measure and test phylogenetic signal
Methods in Ecology and Evolution 2012, 3, 743 756 doi: 10.1111/j.2041-210X.2012.00196.x How to measure and test phylogenetic signal Tamara Mu nkemu ller 1 *, Se bastien Lavergne 1, Bruno Bzeznik 2, Ste
More informationCenter for Applied Mathematics, École Polytechnique, UMR 7641 CNRS, Route de Saclay,
1 Opinion Paper for Trends in Ecology and Evolution (TREE) 2 3 Title: Why does diversification slow down? 4 5 Authors: Daniel Moen 1 and Hélène Morlon 6 7 8 9 10 Center for Applied Mathematics, École Polytechnique,
More informationBio94 Discussion Activity week 3: Chapter 27 Phylogenies and the History of Life
Bio94 Discussion Activity week 3: Chapter 27 Phylogenies and the History of Life 1. Constructing a phylogenetic tree using a cladistic approach Construct a phylogenetic tree using the following table:
More informationA (short) introduction to phylogenetics
A (short) introduction to phylogenetics Thibaut Jombart, Marie-Pauline Beugin MRC Centre for Outbreak Analysis and Modelling Imperial College London Genetic data analysis with PR Statistics, Millport Field
More informationDr. Amira A. AL-Hosary
Phylogenetic analysis Amira A. AL-Hosary PhD of infectious diseases Department of Animal Medicine (Infectious Diseases) Faculty of Veterinary Medicine Assiut University-Egypt Phylogenetic Basics: Biological
More informationPhylogeny-Based Analyses of Evolution with a Paleo-Focus Part 2. Diversification of Lineages (Mar 02, 2012 David W. Bapst
Phylogeny-Based Analyses of Evolution with a Paleo-Focus Part 2. Diversification of Lineages (Mar 02, 2012 David W. Bapst Recommended Reading: Nee, S. 2006. Birth-Death Models in Macroevolution. Annual
More informationModeling continuous trait evolution in the absence of phylogenetic information
Modeling continuous trait evolution in the absence of phylogenetic information Mariia Koroliuk GU university, EM in complex systems under the supervision of Daniele Silvestro June 29, 2014 Contents 1 Introduction
More informationTrees, branches and (square) roots: why evolutionary relatedness is not linearly related to functional distance
Methods in Ecology and Evolution 2015, 6, 439 444 doi: 10.1111/2041-210X.12237 SPECIAL FEATURE NEW OPPORTUNITIES AT THE INTERFACE BETWEEN ECOLOGY AND STATISTICS Trees, branches and (square) roots: why
More informationPhylogeny and systematics. Why are these disciplines important in evolutionary biology and how are they related to each other?
Phylogeny and systematics Why are these disciplines important in evolutionary biology and how are they related to each other? Phylogeny and systematics Phylogeny: the evolutionary history of a species
More information3 Hours 18 / 06 / 2012 EXAMS OFFICE USE ONLY University of the Witwatersrand, Johannesburg Course or topic No(s) ANAT 4000
3 Hours 18 / 06 / 2012 EXAMS OFFICE USE ONLY University of the Witwatersrand, Johannesburg Course or topic No(s) ANAT 4000 Course or topic name(s) Paper Number & title HUMAN BIOLOGY HONOURS: PAPER 1 Examination
More informationPackage WGDgc. June 3, 2014
Package WGDgc June 3, 2014 Type Package Title Whole genome duplication detection using gene counts Version 1.1 Date 2014-06-03 Author Tram Ta, Charles-Elie Rabier, Cecile Ane Maintainer Tram Ta
More informationPolytomies and Phylogenetically Independent Contrasts: Examination of the Bounded Degrees of Freedom Approach
Syst. Biol. 48(3):547 558, 1999 Polytomies and Phylogenetically Independent Contrasts: Examination of the Bounded Degrees of Freedom Approach THEODORE GARLAND, JR. 1 AND RAMÓN DÍAZ-URIARTE Department of
More informationTrees of Unusual Size: Biased Inference of Early Bursts from Large Molecular Phylogenies
Trees of Unusual Size: Biased Inference of Early Bursts from Large Molecular Phylogenies Matthew W. Pennell 1,2,3 *, Brice A. J. Sarver 1,2,3, Luke J. Harmon 1,2,3 1 Institute for Bioinformatics and Evolutionary
More informationCAUGHT IN THE ACT APSQ 01
ANNEXE A.1 CAUGHT IN THE ACT APSQ 01 Document PowerPoint, disponible en téléchargement seulement à l adresse suivante : http://www.apsq.org/sautquantique/tresors.html CAUGHT IN TH E A CT : AG E NTS O F
More informationTESTING FOR DIFFERENT RATES OF CONTINUOUS TRAIT EVOLUTION USING LIKELIHOOD
Evolution, 60(5), 006, pp. 9 933 TESTIG FOR DIFFERET RATES OF COTIUOUS TRAIT EVOLUTIO USIG LIKELIHOOD BRIA C. O MEARA, 1 CÉCILE AÉ, MICHAEL J. SADERSO, 3,4 AD PETER C. WAIWRIGHT 3,5 1 Center for Population
More informationModel Inadequacy and Mistaken Inferences of Trait-Dependent Speciation
1 Model Inadequacy and Mistaken Inferences of Trait-Dependent Speciation Daniel L. Rabosky 1* and Emma E. Goldberg 2* 1 Museum of Zoology and Department of Ecology and Evolutionary Biology, University
More informationSupplementary material to Whitney, K. D., B. Boussau, E. J. Baack, and T. Garland Jr. in press. Drift and genome complexity revisited. PLoS Genetics.
Supplementary material to Whitney, K. D., B. Boussau, E. J. Baack, and T. Garland Jr. in press. Drift and genome complexity revisited. PLoS Genetics. Tree topologies Two topologies were examined, one favoring
More informationANALYSIS OF CHARACTER DIVERGENCE ALONG ENVIRONMENTAL GRADIENTS AND OTHER COVARIATES
ORIGINAL ARTICLE doi:10.1111/j.1558-5646.2007.00063.x ANALYSIS OF CHARACTER DIVERGENCE ALONG ENVIRONMENTAL GRADIENTS AND OTHER COVARIATES Dean C. Adams 1,2,3 and Michael L. Collyer 1,4 1 Department of
More information5/31/17. Week 10; Monday MEMORIAL DAY NO CLASS. Page 88
Week 10; Monday MEMORIAL DAY NO CLASS Page 88 Week 10; Wednesday Announcements: Family ID final in lab Today Final exam next Tuesday at 8:30 am here Lecture: Species concepts & Speciation. What are species?
More information8/23/2014. Phylogeny and the Tree of Life
Phylogeny and the Tree of Life Chapter 26 Objectives Explain the following characteristics of the Linnaean system of classification: a. binomial nomenclature b. hierarchical classification List the major
More informationBiology 211 (2) Week 1 KEY!
Biology 211 (2) Week 1 KEY Chapter 1 KEY FIGURES: 1.2, 1.3, 1.4, 1.5, 1.6, 1.7 VOCABULARY: Adaptation: a trait that increases the fitness Cells: a developed, system bound with a thin outer layer made of
More informationTESTING HYPOTHESES OF CORRELATED EVOLUTION USING PHYLOGENETICALLY INDEPENDENT CONTRASTS: SENSITIVITY TO DEVIATIONS. FROM BROWNIAN MOTION
Syst. Biol. 45(1):27-47, 1996 TESTING HYPOTHESES OF CORRELATED EVOLUTION USING PHYLOGENETICALLY INDEPENDENT CONTRASTS: SENSITIVITY TO DEVIATIONS. FROM BROWNIAN MOTION RAMON D~z-URIARTE~ AND THEODORE GARLAND,
More informationModes of Macroevolution
Modes of Macroevolution Macroevolution is used to refer to any evolutionary change at or above the level of species. Darwin illustrated the combined action of descent with modification, the principle of
More informationChapter 9: Beyond the Mk Model
Chapter 9: Beyond the Mk Model Section 9.1: The Evolution of Frog Life History Strategies Frog reproduction is one of the most bizarrely interesting topics in all of biology. Across the nearly 6,000 species
More information"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200B Spring 2011 University of California, Berkeley
"PRINCIPLES OF PHYLOGENETICS: ECOLOGY AND EVOLUTION" Integrative Biology 200B Spring 2011 University of California, Berkeley B.D. Mishler Feb. 3, 2011. Quantitative character evolution within a cladogram
More informationIntegrating Fossils into Phylogenies. Throughout the 20th century, the relationship between paleontology and evolutionary biology has been strained.
IB 200B Principals of Phylogenetic Systematics Spring 2011 Integrating Fossils into Phylogenies Throughout the 20th century, the relationship between paleontology and evolutionary biology has been strained.
More informationmvmorph:an R package for fitting multivariate evolutionary models to morphometric data
Methods in Ecology and Evolution 215, 6, 1311 1319 doi: 1.1111/241-21X.1242 APPLICATION mvmorph:an R package for fitting multivariate evolutionary models to morphometric data Julien Clavel 1,2 *, Gilles
More informationHominid Evolution What derived characteristics differentiate members of the Family Hominidae and how are they related?
Hominid Evolution What derived characteristics differentiate members of the Family Hominidae and how are they related? Introduction. The central idea of biological evolution is that all life on Earth shares
More informationPatterns of Evolution
Patterns of Evolution A tree that represents an estimate (hypothesis) of evolutionary relatedness is a phylogeny Classifications can be based on groupings within a phylogeny Groupings can be categorized
More informationBiology 559R: Introduction to Phylogenetic Comparative Methods Topics for this week (Jan 27 & 29):
Biology 559R: Introduction to Phylogenetic Comparative Methods Topics for this week (Jan 27 & 29): Statistical estimation of models of sequence evolution Phylogenetic inference using maximum likelihood:
More information