Background

During outbreaks of novel infectious diseases or variant viruses, genomic sequence data serves as a critical resource for identifying transmission pathways and mutation rates. Health authorities reconstruct the evolutionary history of the virus using phylogenetic trees to estimate when and where the virus was introduced and how it spread. The technique most widely recognized as the most accurate statistical framework in this process is Bayesian phylogenetics. This is because it can simultaneously estimate key epidemiological metrics, such as evolutionary rate and reproduction number, while systematically accounting for uncertainty.

The problem is computational load. Existing Bayesian phylogenetic tools use Markov Chain Monte Carlo (MCMC) algorithms. As sample sizes exceed hundreds, computation time increases exponentially, often taking days or weeks to complete analysis. During a pandemic, when tens or hundreds of thousands of whole-genome sequencing datasets are rapidly generated, existing methods struggle to contribute in time to real-time epidemic prevention decision-making. Developing nations or regional health agencies lacking high-performance computing infrastructure are unable to even attempt to apply the latest phylogenetic analysis results directly to disease control efforts.

Key Findings

The research team has unveiled Delphy, a new computational framework capable of processing accumulating large-scale viral genomic data in near real-time. The core of Delphy is an online Bayesian update structure that organically integrates new sequences into the previously learned posterior probability distribution of the phylogenetic tree, rather than recalculating the entire tree from scratch every time new data arrives.

Thanks to algorithmic optimization, Delphy drastically reduces computational complexity even as samples accumulate. While existing tools took dozens of hours on high-performance clusters to process thousands of sequences, Delphy completed Bayesian analysis within dozens of minutes on standard workstation-class computers. It maintained the highest level of statistical accuracy while increasing speed. Verification through simulations and actual infectious disease genomic datasets demonstrated that the posterior probability distributions of divergence time estimates and phylogenetic tree topologies achieved accuracy consistent with traditional full MCMC methods. This means high-level phylodynamic analysis can be performed with minimal computational cost even in resource-limited environments.

Significance and Outlook

The emergence of Delphy marks a turning point in shifting the paradigm of genomic surveillance from centralized post-hoc analysis to decentralized, real-time local surveillance. This is because local laboratories can perform phylogenetic analysis immediately as data is generated and feed the results back into quarantine operations. It provides a foundation for public health agencies worldwide to independently operate standardized Bayesian precision analysis while maintaining their own data sovereignty.

Challenges remain. For pathogens with high recombination rates or complex insertion/deletion mutations, additional phylogenetic modeling beyond simple point mutation models is required. The risk of bias during long-term updates cannot be ruled out if quality variations or sequencing errors in new data accumulate. Follow-up research will require expanding pipelines to flexibly accommodate various mutation mechanisms and automate the data cleaning process.