Incorporating heterogeneous sampling probabilities in continuous phylogeographic inference — Application to H5N1 spread in the Mekong region

Abstract

Motivation The potentially low precision associated with the geographic origin of sampled sequences represents an important limitation for spatially explicit (i.e. continuous) phylogeographic inference of fast-evolving pathogens such as RNA viruses. A substantial proportion of publicly available sequences is geo-referenced at broad spatial scale such as the administrative unit of origin, rather than more precise locations (e.g. geographic coordinates). Most frequently, such sequences are either discarded prior to continuous phylogeographic inference or arbitrarily assigned to the geographic coordinates of the centroid of their administrative area of origin for lack of a better alternative. Results We here implement and describe a new approach that allows to incorporate heterogeneous prior sampling probabilities over a geographic area. External data, such as outbreak locations, are used to specify these prior sampling probabilities over a collection of sub-polygons. We apply this new method to the analysis of highly pathogenic avian influenza H5N1 clade data in the Mekong region. Our method allows to properly include, in continuous phylogeographic analyses, H5N1 sequences that are only associated with large administrative areas of origin and assign them with more accurate locations. Finally, we use continuous phylogeographic reconstructions to analyse the dispersal dynamics of different H5N1 clades and investigate the impact of environmental factors on lineage dispersal velocities. Availability and implementation Our new method allowing heterogeneous sampling priors for continuous phylogeographic inference is implemented in the open-source multi-platform software package BEAST 1.10. Supplementary information Supplementary data are available at Bioinformatics online.

Accounting for spatial sampling patterns in Bayesian phylogeography

“…The signal conveyed by the data about this parameter (as well as the exponential growth parameter) is weak though as its posterior distribution is heavily influenced by the prior (SI Appendix, section 3). While dispersal parameter estimates are similar to that obtained in previous studies (39,40), values of that parameter inferred under the survey scheme are noticeably larger than under the detection scheme, thereby suggesting long dispersal events in short time frames (Fig. 1E).…”

Genomic Epidemiology, Evolution, and Transmission Dynamics of Porcine Deltacoronavirus

“…Because the sampling locations of several available sequences were only known at the province level, we used a sampling prior approach to define a potential area of origin for these sequences ( Nylinder et al. 2014 ; Dellicour et al. 2020 ).…”

Integrating animal movements with phylogeography to model the spread of PRRSV in the USA

“…Primarily, these approaches have enabled the inclusion of spatial and environmental data in phylodynamic analysis to understand the spread of different epidemics, e.g. Ebola, influenza and HIV in humans ( Dudas et al, 2017 ; Müller, Rasmussen, and Stadler 2018 ; Rasmussen et al, 2018 ), foot-and-mouth disease in livestock ( Duchatel, Bronsvoort, and Lycett 2019 ; Munsey et al, 2021 ), and infectious wildlife diseases ( Fountain-Jones et al, 2017 ; Yang et al, 2019 ; Dellicour et al, 2020 ). Given the role of long-distance animal movements in the spread of PRRS and other livestock diseases ( Mortensen et al, 2002 ; Carlsson et al, 2009 ; Amirpour et al, 2017 ; Neira et al, 2017 ; VanderWaal et al, 2020 ; Makau et al, 2021 ), the inclusion of empirical animal movement data in these models could increase the accuracy and robustness of phylodynamic inference and is essential for capturing how host population connectivity interacts with spatial drivers of transmission.…”