Inferring epidemiological parameters from phylogenies using regression-ABC: a comparative study

Archive ouverte

Saulnier, Emma | Gascuel, Olivier | Alizon, Samuel

Edité par CCSD ; PLOS -

International audience. Inferring epidemiological parameters such as the R0 from time-scaled phylogenies is a timely challenge. Most current approaches rely on likelihood functions, which raise specific issues that range from computing these functions to finding their maxima numerically. Here, we present a new regression-based Approximate Bayesian Computation (ABC) approach, which we base on a large variety of summary statistics intended to capture the information contained in the phylogeny and its corresponding lineage-through-time plot. The regression step involves the Least Absolute Shrinkage and Selection Operator (LASSO) method, which is a robust machine learning technique. It allows us to readily deal with the large number of summary statistics, while avoiding resorting to Markov Chain Monte Carlo (MCMC) techniques. To compare our approach to existing ones, we simulated target trees under a variety of epidemiological models and settings, and inferred parameters of interest using the same priors. We found that, for large phylogenies, the accuracy of our regression-ABC is comparable to that of likelihood-based approaches involving birth-death processes implemented in BEAST2. Our approach even outperformed these when inferring the host population size with a Susceptible-Infected-Removed epidemiological model. It also clearly outperformed a recent kernel-ABC approach when assuming a Susceptible-Infected epidemiological model with two host types. Lastly, by re-analyzing data from the early stages of the recent Ebola epidemic in Sierra Leone, we showed that regression-ABC provides more realistic estimates for the duration parameters (latency and infectiousness) than the likelihood-based method. Overall, ABC based on a large variety of summary statistics and a regression method able to perform variable selection and avoid overfitting is a promising approach to analyze large phylogenies.

Suggestions

Du même auteur

TiPS: rapidly simulating trajectories and phylogenies from compartmental models

Archive ouverte | Danesh, Gonché | CCSD

Abstract We introduce TiPS , an R package to generate trajectories and phylogenetic trees associated with a compartmental model. Trajectories are simulated using Gillespie’s exact or approximate stochastic simulation algorithm, or...

TiPS : Rapidly simulating trajectories and phylogenies from compartmental models

Archive ouverte | Danesh, Gonché | CCSD

International audience. Abstract Stochastic population dynamics simulations are essential for many ecological and epidemiological studies to generate time series and genealogies that capture the relatedness between ...

Phylodynamics of viral infections. Phylodynamique des infections virales

Archive ouverte | Alizon, Samuel | CCSD

International audience. Phylodynamics is a recent field that aims at using genetic sequence data to estimate epidemiological parameters such as the viral population growth rate, the number of infections in the popul...

Chargement des enrichissements...