A new method for inferring timetrees from temporally sampled molecular sequences

Abstract

ABSTRACTPathogen timetrees are phylogenies scaled to time. They reveal the temporal history of a pathogen spread through the populations as captured in the evolutionary history of strains. These timetrees are inferred by using molecular sequences of pathogenic strains sampled at different times. That is, temporally sampled sequences enable the inference of sequence divergence times. Here, we present a new approach (RelTime with Dated Tips [RTDT]) to estimating pathogen timetrees based on the relative rate framework underlying the RelTime approach. RTDT does not require many of the priors demanded by Bayesian approaches, and it has light computing requirements. We found RTDT to be accurate on simulated datasets evolved under a variety of branch rates models. Interestingly, we found two non-Bayesian methods (RTDT and Least Squares Dating [LSD]) to perform similar to or better than the Bayesian approaches available in BEAST and MCMCTree programs. RTDT method was found to generally outperform all other methods for phylogenies in with autocorrelated evolutionary rates. In analyses of empirical datasets, RTDT produced dates that were similar to those from Bayesian analyses. Speed and accuracy of the new method, as compared to the alternatives, makes it appealing for analyzing growing datasets of pathogenic strains. Cross-platform MEGA X software, freely available from http://www.megasoftware.net, now contains the new method for use through a friendly graphical user interface and in high-throughput settings.AUTHOR SUMMARYPathogen timetrees trace the origins and evolutionary histories of strains in populations, hosts, and outbreaks. The tips of these molecular phylogenies often contain sampling time information because the sequences were generally obtained at different times during the disease outbreaks and propagation. We have developed a new method for inferring timetrees for phylogenies with tip dates, which improves on widely-used Bayesian methods (e.g., BEAST) in computational efficiency and does not require prior specification of population parameters, branch rate model, or clock model. We performed extensive computer simulation and found that RTDT performed better than the other methods for the estimation of divergence times at deep node in phylogenies where evolutionary rates were autocorrelated. The new method is available in the cross-platform MEGA software package that provides a graphical user interface, and allows use via a command line in scripting and high throughput analysis (www.megasoftware.net).