Efficient inference, potential, and limitations of site-specific substitution models

Puller, Vadim and Sagulenko, Pavel and Neher, Richard A.. (2020) Efficient inference, potential, and limitations of site-specific substitution models. Virus Evolution, 6 (2). veaa066.

PDF - Published Version
Available under License CC BY (Attribution).


Official URL: https://edoc.unibas.ch/84945/

Downloads: Statistics Overview


Natural selection imposes a complex filter on which variants persist in a population resulting in evolutionary patterns that vary greatly along the genome. Some sites evolve close to neutrally, while others are highly conserved, allow only specific states, or only change in concert with other sites. On one hand, such constraints on sequence evolution can be to infer biological function, one the other hand they need to be accounted for in phylogenetic reconstruction. Phylogenetic models often account for this complexity by partitioning sites into a small number of discrete classes with different rates and/or state preferences. Appropriate model complexity is typically determined by model selection procedures. Here, we present an efficient algorithm to estimate more complex models that allow for different preferences at every site and explore the accuracy at which such models can be estimated from simulated data. Our iterative approximate maximum likelihood scheme uses information in the data efficiently and accurately estimates site-specific preferences from large data sets with moderately diverged sequences and known topology. However, the joint estimation of site-specific rates, and site-specific preferences, and phylogenetic branch length can suffer from identifiability problems, while ignoring variation in preferences across sites results in branch length underestimates. Site-specific preferences estimated from large HIV; pol; alignments show qualitative concordance with intra-host estimates of fitness costs. Analysis of these substitution models suggests near saturation of divergence after a few hundred years. Such saturation can explain the inability to infer deep divergence times of HIV and SIVs using molecular clock approaches and time-dependent rate estimates.
Faculties and Departments:05 Faculty of Science > Departement Biozentrum > Computational & Systems Biology > Computational Modeling of Biological Processes (Neher)
UniBasel Contributors:Neher, Richard A
Item Type:Article, refereed
Article Subtype:Research Article
Publisher:Oxford University Press
Note:Publication type according to Uni Basel Research Database: Journal article
Related URLs:
Identification Number:
edoc DOI:
Last Modified:24 Feb 2022 11:53
Deposited On:24 Feb 2022 11:53

Repository Staff Only: item control page