Abstract
ABSTRACTThe coalescent is a powerful statistical framework that allows us to infer past population dynamics leveraging the ancestral relationships reconstructed from sampled molecular sequence data. In many biomedical applications, such as in the study of infectious diseases, cell development, and tumorgenesis, several distinct populations share evolutionary history and therefore become dependent. The inference of such dependence is a highly important, yet a challenging problem. With advances in sequencing technologies, we are well positioned to exploit the wealth of high-resolution biological data for tackling this problem. Here, we present a novel probabilistic model that relies on jointly distributed Markov random fields. We use this model to estimate past population dynamics of dependent populations and to quantify their degree of dependence. An essential feature of our approach is the ability to track the time-varying association between the populations while making minimal assumptions on their functional shapes via Markov random field priors. We provide nonparametric estimators, extensions of our base model that integrate multiple data sources, and fast scalable inference algorithms. We test our method using simulated data under various dependent population histories and demonstrate the utility of our model in shedding light on evolutionary histories of different variants of SARS-CoV-2.
Publisher
Cold Spring Harbor Laboratory
Reference42 articles.
1. Adams, R. P. , Murray, I. , and MacKay, D. J. Tractable nonparametric Bayesian inference in Poisson processes with Gaussian process intensities. In Proceedings of the 26th Annual International Conference on Machine Learning, pp. 9–16, 2009.
2. Bishop, C. M. Pattern recognition and machine learning. Springer, New York, NY, 2006.
3. BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis;PLoS Computational Biology,2019
4. Adaptive preferential sampling in phylodynamics with an application to SARS-CoV-2;Journal of Computational and Graphical Statistics,2021
5. Cappello, L. , Kim, J. , Liu, S. , and Palacios, J. A. Statistical challenges in tracking the evolution of SARS-CoV-2. arXiv, 2108.13362 [stat.AP], 2021.