How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?

One of the central objectives in the field of phylodynamics is the quantification of population dynamic processes using genetic sequence data or in some cases phenotypic data. Phylodynamics has been successfully applied to many different processes, such as the spread of infectious diseases, within-h...

Full description

Bibliographic Details
Main Authors: Stadler, Tanja, Vaughan, Timothy G., Gavryushkin, Alex, Guindon, Stephane, Kühnert, Denise, Leventhal, Gabriel E., Drummond, Alexei J.
Format: Online
Language:English
Published: The Royal Society 2015
Online Access:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4426635/
id pubmed-4426635
recordtype oai_dc
spelling pubmed-44266352015-05-21 How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics? Stadler, Tanja Vaughan, Timothy G. Gavryushkin, Alex Guindon, Stephane Kühnert, Denise Leventhal, Gabriel E. Drummond, Alexei J. Research Articles One of the central objectives in the field of phylodynamics is the quantification of population dynamic processes using genetic sequence data or in some cases phenotypic data. Phylodynamics has been successfully applied to many different processes, such as the spread of infectious diseases, within-host evolution of a pathogen, macroevolution and even language evolution. Phylodynamic analysis requires a probability distribution on phylogenetic trees spanned by the genetic data. Because such a probability distribution is not available for many common stochastic population dynamic processes, coalescent-based approximations assuming deterministic population size changes are widely employed. Key to many population dynamic models, in particular epidemiological models, is a period of exponential population growth during the initial phase. Here, we show that the coalescent does not well approximate stochastic exponential population growth, which is typically modelled by a birth–death process. We demonstrate that introducing demographic stochasticity into the population size function of the coalescent improves the approximation for values of R0 close to 1, but substantial differences remain for large R0. In addition, the computational advantage of using an approximation over exact models vanishes when introducing such demographic stochasticity. These results highlight that we need to increase efforts to develop phylodynamic tools that correctly account for the stochasticity of population dynamic models for inference. The Royal Society 2015-05-07 /pmc/articles/PMC4426635/ /pubmed/25876846 http://dx.doi.org/10.1098/rspb.2015.0420 Text en http://creativecommons.org/licenses/by/4.0/ © 2015 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/by/4.0/, which permits unrestricted use, provided the original author and source are credited.
repository_type Open Access Journal
institution_category Foreign Institution
institution US National Center for Biotechnology Information
building NCBI PubMed
collection Online Access
language English
format Online
author Stadler, Tanja
Vaughan, Timothy G.
Gavryushkin, Alex
Guindon, Stephane
Kühnert, Denise
Leventhal, Gabriel E.
Drummond, Alexei J.
spellingShingle Stadler, Tanja
Vaughan, Timothy G.
Gavryushkin, Alex
Guindon, Stephane
Kühnert, Denise
Leventhal, Gabriel E.
Drummond, Alexei J.
How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
author_facet Stadler, Tanja
Vaughan, Timothy G.
Gavryushkin, Alex
Guindon, Stephane
Kühnert, Denise
Leventhal, Gabriel E.
Drummond, Alexei J.
author_sort Stadler, Tanja
title How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
title_short How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
title_full How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
title_fullStr How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
title_full_unstemmed How well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
title_sort how well can the exponential-growth coalescent approximate constant-rate birth–death population dynamics?
description One of the central objectives in the field of phylodynamics is the quantification of population dynamic processes using genetic sequence data or in some cases phenotypic data. Phylodynamics has been successfully applied to many different processes, such as the spread of infectious diseases, within-host evolution of a pathogen, macroevolution and even language evolution. Phylodynamic analysis requires a probability distribution on phylogenetic trees spanned by the genetic data. Because such a probability distribution is not available for many common stochastic population dynamic processes, coalescent-based approximations assuming deterministic population size changes are widely employed. Key to many population dynamic models, in particular epidemiological models, is a period of exponential population growth during the initial phase. Here, we show that the coalescent does not well approximate stochastic exponential population growth, which is typically modelled by a birth–death process. We demonstrate that introducing demographic stochasticity into the population size function of the coalescent improves the approximation for values of R0 close to 1, but substantial differences remain for large R0. In addition, the computational advantage of using an approximation over exact models vanishes when introducing such demographic stochasticity. These results highlight that we need to increase efforts to develop phylodynamic tools that correctly account for the stochasticity of population dynamic models for inference.
publisher The Royal Society
publishDate 2015
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4426635/
_version_ 1613221594157547520