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...
Main Authors: | , , , , , , |
---|---|
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 |