Posterior Summarization in Bayesian Phylogenetics Using Tracer 1.7 (original) (raw)

Abstract

Bayesian inference of phylogeny using Markov chain Monte Carlo (MCMC) plays a central role in understanding evolutionary history from molecular sequence data. Visualizing and analyzing the MCMC-generated samples from the posterior distribution is a key step in any non-trivial Bayesian inference. We present the software package Tracer (version 1.7) for visualizing and analyzing the MCMC trace files generated through Bayesian phylogenetic inference. Tracer provides kernel density estimation, multivariate visualization, demographic trajectory reconstruction, conditional posterior distribution summary, and more. Tracer is open-source and available at http://beast.community/tracer.

Keywords: Bayesian inference, Markov chain Monte Carlo, phylogenetics, visualization

Bayesian inference of phylogeny using Markov chain Monte Carlo (MCMC) (Rannala and Yang 1996; Mau et al. 1999; Drummond et al. 2002) flourishes as a popular approach to uncover the evolutionary relationships among taxa, such as genes, genomes, individuals, or species. MCMC approaches generate samples of model parameter values—including the phylogenetic tree—drawn from their posterior distribution given molecular sequence data and a selection of evolutionary models. Visualizing, tabulating, and marginalizing these samples are critical for approximating the posterior quantities of interest that one reports as the outcome of a Bayesian phylogenetic analysis. To facilitate this task, we have developed the Tracer (version 1.7) software package to process MCMC trace files containing parameter samples and to interactively explore the high-dimensional posterior distribution. Tracer works automatically with sample output from BEAST (Drummond et al. 2012), BEAST2 (Bouckaert et al. 2014), LAMARC (Kuhner 2006), Migrate (Beerli 2006), MrBayes (Ronquist et al. 2012), RevBayes (Höhna et al. 2016), and possibly other MCMC programs from other domains.

Design and Implementation

Tracer examines the posterior samples from all the available parameters—treating continuous, integer and categorical parameters appropriately—from a trace and presents statistical summaries and visualizations. Further, Tracer can analyze a single trace or combine samples from multiple files. Immediately apparent in the default Tracer view, the effective sample size (ESS) is one such statistic that allows users to assess the number of effectively independent draws from the posterior distribution the trace represents (Figure 1a). Color coding assists the user in determining potential MCMC mixing problems, with arbitrary cut-off values at 100 and 200.

Figure 1.

Figure 1.

Overview of Tracer functionality and individual parameter visualizations: a) Main Tracer panel upon loading a single trace file; b) boxplot representation of two continuous parameters; corresponding c) kernel density estimates; d) violin plots; e) the actual traces connecting the parameter values visited by the Markov chain.

Selecting multiple parameters from the “Traces” panel on the left generates a side-by-side comparison or an overlay of the selected parameters’ visualizations (Figure 1 b–e). Multiple trace files can be selected in a similar fashion to compare posterior samples between different replicates of an analysis. If multiple trace files contain the same collection of parameters, then a “Combined” trace appears automatically. Tracer generates four display panels for the selected parameters:

Tracer offers a solution of visualizing conditional posterior distributions as well. Selecting one continuous and one integer or categorical parameter generates side-by-side violin or boxplots under the Joint-Marginal panel. These plots present the continuous parameter distribution conditioned on the unique integer or categorical values. A typical use case involves Bayesian stochastic search variable selection (BSSVS), a form of model averaging, in which parameters influence the likelihood function only when a specific model is selected by a random indicator function. Under BSSVS, a posterior estimate of the parameter should only sample values from states where the indicator equals one. Discrete phylogeographic analyses frequently employ BSSVS due to the potentially large amount of transition rates that need to be estimated (Lemey et al. 2009), but this is also relevant when employing model averaging approaches, e.g., over relaxed molecular clocks (Li and Drummond 2012).

Finally, Tracer provides demographic reconstruction resulting in a graphical plot, often applied to reconstruct epidemic dynamics. Available models include, e.g., constant size, exponential and logistic growth (Drummond et al. 2002), and the non-parametric Bayesian skyline (Drummond et al. 2005; Heled and Drummond 2008), skyride (Minin et al. 2008), and skygrid (Gill et al. 2013).

Example

Cross-Species Dynamics of North American Bat Rabies

We use Tracer to infer the spatial dispersal and cross-species dynamics of rabies virus (RABV) in North American bats. The data set comprises 372 nucleoprotein gene sequences from 17 bat species, sampled between 1997 and 2006 across 14 states in the United States (Streicker et al. 2010; Faria et al. 2013). We estimate RABV ancestral locations and host-jumping history using a Bayesian discrete phylogeographic approach with BSSVS, while simultaneously estimating effective population sizes over time through a Bayesian skygrid coalescent model (Gill et al. 2013).

Phylogeographic BSSVS inference includes parameters of both integer (number of non-zero transition rates) and categorical (host or location-state) trace types. In Tracer, a bubble chart visualizes the joint probability distribution between two integer or categorical traces (see Figure 2a). Circle area is proportional to the joint probability, with a colored tile background if this probability reaches a nominal threshold to enhance visibility. Marginal density plots can also display multiple integer parameters, each with unique colour scales (see Figure 2b). With approximately equal numbers of transition rates, both figures suggest similar host and location trait model complexity. Tracer also provides popular visualizations for continuous parameters, including scatter plots for two parameters (see Figure 2c), and extensions for correlations between Inline graphic continuous parameters (Figure 2d; Murdoch and Chow 1996). Colour gradients indicate strength and direction of the correlation, from red (strong negative) to blue (strong positive). Ellipse shapes re-enforce the strength of correlation, with no correlation appearing as a circle and perfect (anti)correlation as a line.

Figure 2.

Figure 2.

Multi-parameter visualizations of: a) the joint probability distribution of two integer variables through a bubble chart; b) the marginal density of multiple integer or categorical variables through frequency plots; c) two continuous variables through a classic scatter or correlation plot; d) multiple (Inline graphic) continuous variables using large correlation matrices.

Tracer reconstructs the demographic history of RABV by drawing the effective population sizes over time (Figure 3). RABV has successfully established itself in North American bat species, with its effective population size rising steadily throughout recent centuries. Following a rapid decline at the end of last century, we observe a recent sharp increase in size.

Figure 3.

Figure 3.

Estimating the effective population sizes over time using a Bayesian skygrid demographic reconstruction for rabies virus in North America.

Other packages are available for the post-processing of MCMC samples. “coda” (Plummer et al. 2006) provides some of the functionality of Tracer within the R programming environment, while “AWTY” (Nylander et al. 2007) and “RWTY” (Warren et al. 2017) explore the convergence of the phylogenetic tree parameter itself across multiple MCMC runs. These alternative packages compute, e.g., Gelman–Rubin diagnostics (Gelman and Rubin 1992) that Tracer currently does not provide.

Availability

Tracer is open-source under the GNU lesser general public license and available in both source code (https://github.com/beast-dev/tracer) and executable (http://beast.community/tracer) forms. This latter page also serves up self-contained, step-by-step tutorials covering basic to advanced usage of Tracer to summarize posteriors under a variety of phylogenetic models using BEAST and diagnose MCMC chain convergence. Popular tutorials employ Tracer to generate marginal parameter summaries and to infer population dynamics trajectories over time. Tracer requires Java version 1.6 or greater.

Funding

This work was supported in part by the Wellcome Trust through project 206298/Z/17/Z (Artic Network), the European Union Seventh Framework Programme under [grant agreement no. 725422-RESERVOIRDOCS]; the National Science Foundation through grant DMS [1264153]; the National Institutes of Health under grants [R01 AI107034 and U19 AI135995]; Marsden grant contract [UOA1611 to A.J.D.]; Interne Fondsen KU Leuven / Internal Funds KU Leuven to G.B.

References

  1. Beerli P. (2006). Comparison of Bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics 22:341–345. [DOI] [PubMed] [Google Scholar]
  2. Bouckaert R.,, Heled J.,, Kühnert D.,, Vaughan T.,, Wu C.-H.,, Xie D.,, Suchard M.A.,, Rambaut A.,, Drummond A.J. (2014). BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comp. Biol. 10:e1003537. [DOI] [PMC free article] [PubMed] [Google Scholar]
  3. Drummond, A.J.,, Nicholls, G.K.,, Rodrigo, A.G.,, Solomon W. (2002). Estimating mutation parameters, population history and genealogy simultaneously from temporally spaced sequence data. Genetics 161:1307–1320. [DOI] [PMC free article] [PubMed] [Google Scholar]
  4. Drummond, A.J.,, Rambaut, A.,, Shapiro, B.,, Pybus, O.G. (2005). Bayesian coalescent inference of past population dynamics from molecular sequences. Mol. Biol. Evol. 22:1185–1192. [DOI] [PubMed] [Google Scholar]
  5. Drummond A.J.,, Suchard M.A., Xie D.,, Rambaut A. (2012). Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol. Evol. 29:1969–1973. [DOI] [PMC free article] [PubMed] [Google Scholar]
  6. Faria, N.,, Suchard M.,, Rambaut A.,, Streicker D.,, Lemey P. (2013). Simultaneously reconstructing viral cross-species transmission history and identifying the underlying constraint. Philos. Trans. R. Soc. Lon. B 368:20120196. [DOI] [PMC free article] [PubMed] [Google Scholar]
  7. Gelman A.,, Rubin D.B. (1992). Inference from iterative simulation using multiple sequences. Stat. Sci. 7:457–472. [Google Scholar]
  8. Gill, M.S.,, Lemey, P.,, Faria N.R.,, Rambaut A.,, Shapiro B.,, Suchard M.A. (2013). Improving Bayesian population dynamics inference: a coalescent-based model for multiple loci. Mol. Biol. Evol. 30:713–724. [DOI] [PMC free article] [PubMed] [Google Scholar]
  9. Heled J.,, Drummond A.J. (2008). Bayesian inference of population size history from multiple loci. BMC Evol. Biol. 8:289. [DOI] [PMC free article] [PubMed] [Google Scholar]
  10. Höhna S.,, Landis M.J.,, Heath T.A.,, Boussau B.,, Lartillot N.,, Moore B.R.,, Huelsenbeck J.P.,, Ronquist F. (2016). RevBayes: Bayesian phylogenetic inference using graphical models and an interactive model-specification language. Syst. Biol. 65:726–736. [DOI] [PMC free article] [PubMed] [Google Scholar]
  11. Kuhner M.K. (2006). LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics 22:768–770. [DOI] [PubMed] [Google Scholar]
  12. Lemey P.,, Rambaut A.,, Drummond A.,, Suchard M. (2009). Bayesian phylogeography finds its roots. PLoS Comp. Biol. 5:e1000520. [DOI] [PMC free article] [PubMed] [Google Scholar]
  13. Li W.L.S.,, Drummond A.J. (2012). Model averaging and Bayes factor calculation of relaxed molecular clocks in Bayesian phylogenetics. Mol. Biol. Evol. 29:751–761. [DOI] [PMC free article] [PubMed] [Google Scholar]
  14. Mau B.,, Newton M.A.,, Larget B. (1999). Bayesian phylogenetic inference via Markov chain Monte Carlo methods. Biometrics 55:1–12. [DOI] [PubMed] [Google Scholar]
  15. Minin V.N.,, Bloomquist E.W.,, Suchard M.A. (2008). Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics. Mol. Biol. Evol. 25:1459–1471. [DOI] [PMC free article] [PubMed] [Google Scholar]
  16. Murdoch D.,, Chow E. (1996). A graphical display of large correlation matrices. Am. Stat. 50:178–180. [Google Scholar]
  17. Nylander J.A.,, Wilgenbusch J.C.,, Warren D.L.,, Swofford D.L. (2007). AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics. Bioinformatics 24:581–583. [DOI] [PubMed] [Google Scholar]
  18. Plummer M.,, Best N.,, Cowles K.,, Vines K. (2006). CODA: convergence diagnosis and output analysis for MCMC. R News 6:7–11. [Google Scholar]
  19. Rannala B.,, Yang Z. (1996). Probability distribution of molecular evolutionary trees: a new method of phylogenetic inference. J. Mol. Evol. 43:304–311. [DOI] [PubMed] [Google Scholar]
  20. Ronquist F.,, Teslenko M.,, van der Mark P.,, Ayres D.L.,, Darling A.,, Höhna S.,, Larget B.,, Liu L.,, Suchard M.A.,, Huelsenbeck J.P. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61:539–542. [DOI] [PMC free article] [PubMed] [Google Scholar]
  21. Streicker D.,, Turmelle A.,, Vonhof M.,, Kuzmin I.,, McCracken G.F.,, Rupprecht C. (2010). Host phylogeny constrains cross-species emergence and establishment of rabies virus in bats. Science 329:676–679. [DOI] [PubMed] [Google Scholar]
  22. Warren D.L.,, Geneva A.J.,, Lanfear R. (2017). RWTY (R We There Yet): an R package for examining convergence of Bayesian phylogenetic analyses. Mol. Biol. Evol. 34:1016–1020. [DOI] [PubMed] [Google Scholar]