A new study published in Nature introduces Delphy, a scalable Bayesian phylogenetics platform designed from the ground up for near-real-time analysis of outbreak-scale datasets. The work addresses a bottleneck that has frustrated the genomics community since 2020: how to preserve the statistical rigor of Bayesian inference while making it fast enough to run continuously on datasets containing hundreds of thousands of genomes. The authors demonstrate that Delphy can produce dated, uncertainty-aware phylogenies from massive SARS-CoV-2 datasets on timescales of hours rather than days or weeks, opening the door to a style of genomic epidemiology in which phylogenetic estimates are refreshed as routinely as case counts.

The core challenge Delphy tackles is one of combinatorial explosion. Bayesian phylogenetic inference requires exploring an astronomically large space of possible trees, branch lengths, and substitution model parameters, typically via Markov chain Monte Carlo (MCMC) sampling. Each step of an MCMC chain proposes a small change to the current state of the tree and evaluates how well the resulting configuration explains the sequence data, a calculation whose cost scales with both the number of taxa and the length of the genome alignment. With a few hundred genomes, a desktop workstation can handle the workload. With hundreds of thousands, the likelihood evaluations become prohibitively expensive, and the mixing of the chain—the efficiency with which it explores the space of trees—deteriorates because the vast majority of proposed moves are rejected. Traditional implementations such as BEAST, which has served the field admirably for years, were architected for datasets orders of magnitude smaller than those generated during the pandemic.

Delphy’s designers approached the problem by rethinking the inference engine rather than simply optimizing existing code. The platform restructures the likelihood computation to take advantage of modern hardware accelerators, particularly graphics processing units (GPUs), which excel at the kind of massively parallel arithmetic that underlies phylogenetic likelihood calculations. In a phylogenetic likelihood evaluation, the probability of observing the sequence data at each site of the genome must be computed for every node of the tree, and these computations are largely independent of one another. By mapping them onto thousands of parallel processing cores, Delphy amortizes the cost of each likelihood evaluation across hardware that consumer laptops and cloud instances now routinely carry. The result is a dramatic speedup in the inner loop of Bayesian inference, the very component that dominated runtime in earlier tools.

Hardware alone, however, would not have been sufficient. The study also describes algorithmic innovations in how the chain explores tree space and how the analysis is partitioned. Rather than requiring the entire dataset to be re-analyzed from scratch each time new genomes arrive, Delphy is built to support incremental updates, in which an existing posterior sample of trees is extended and refined as additional sequences become available. This design mirrors the operational rhythm of outbreak response, where sequencing data arrive in a continuous stream and analysts need updated estimates on a daily or even hourly cadence. The ability to warm-start an analysis from a previous posterior, rather than beginning each run cold, converts what was historically a batch process into something closer to a living model of an epidemic.

Another distinguishing feature of the platform is its treatment of time. During the COVID-19 pandemic, the field converged on the value of time-calibrated phylogenies, in which branch lengths are expressed in units of calendar time using sampling dates as calibration points. These dated trees allow epidemiologists to estimate when lineages diverged, when a variant likely entered a region, and how quickly lineages were growing or declining in frequency. Delphy incorporates temporal structure directly into its Bayesian model, jointly estimating the evolutionary rate and the timing of divergence events alongside the tree topology. Crucially, it retains full posterior distributions over these quantities, so that an analyst can report not just a point estimate of when a lineage emerged but a credible interval that honestly reflects the limits of the data. In fast-moving outbreak situations, where decisions about border measures, vaccine updates, and public communications hinge on timing estimates, that quantified uncertainty is not a luxury; it is the difference between a defensible inference and an overconfident one.

The authors validate the platform on datasets drawn from the SARS-CoV-2 pandemic, the largest real-world stress test genomic epidemiology has ever faced. Global repositories accumulated many millions of viral genomes during that outbreak, and even curated subsets routinely exceeded the practical limits of legacy Bayesian tools. Delphy’s benchmarks show that analyses on datasets of this scale, which would have been effectively impossible with prior software, complete in hours on appropriate hardware, with posterior estimates that are consistent with those obtained by slower reference methods on smaller subsets. This kind of cross-validation matters: speed is worthless if the fast answer is wrong. By demonstrating agreement between its accelerated inference and established approaches on tractable data, the study builds the case that Delphy’s approximations and engineering choices preserve the statistical integrity that Bayesian methods are meant to guarantee.

The implications for public health practice extend well beyond coronaviruses. The same architecture applies to any rapidly evolving pathogen for which dense genomic sampling is available, including influenza, respiratory syncytial virus, dengue, mpox, and foodborne bacterial outbreaks. In each of these settings, the operational question is similar: can a public health laboratory turn a week’s worth of new genomes into an updated picture of transmission before the picture changes again? A tool that delivers dated phylogenies with quantified uncertainty on a near-real-time cadence shifts genomic epidemiology from retrospective analysis toward genuine situational awareness. It also lowers the barrier for smaller laboratories and public health agencies in low-resource settings, since the computational demands, while substantial, are within reach of cloud computing budgets that are far smaller than the supercomputing clusters previously imagined necessary for pandemic-scale Bayesian work.

The study also speaks to a broader lesson about scientific software. The pandemic exposed a gap between the methods developed in academic statistics and systematics communities and the operational needs of public health. Tools that were exemplary for studies of dozens or hundreds of taxa could not simply be scaled by throwing more processors at them; they required rethinking data structures, sampling strategies, and the software engineering practices needed for continuous, reliable operation. Delphy represents the kind of ground-up redesign that gap demands, combining modern statistical machinery with the performance engineering of contemporary high-performance computing. Its open availability means that the wider community can scrutinize, extend, and build upon the platform, which is essential if it is to become trusted infrastructure rather than a one-off demonstration.

Challenges remain, of course. Bayesian inference on enormous datasets still requires careful attention to model adequacy, because a fast and precise answer to a poorly specified model is of limited value. Recombination, which is substantial in coronaviruses, complicates the strictly tree-based assumptions underlying most phylogenetic software, and integrating recombination-aware inference at pandemic scale remains an open problem. Sampling bias—where genomes are sequenced unevenly across geography, time, and severity of illness—continues to shape what any phylogenetic method can honestly conclude. None of these caveats diminishes the significance of the advance, but they frame the work realistically: Delphy removes a computational barrier that has constrained the field for years, leaving the scientific community freer to focus on the biological and epidemiological subtleties that no amount of computing power can resolve on its own.

If the trajectory of genomic surveillance continues, the coming years will see sequencing embedded ever more deeply in routine public health, from wastewater monitoring to hospital diagnostics. In that world, the demand for phylogenetic analyses that are simultaneously rigorous, current, and scalable will only grow. Delphy’s demonstration that Bayesian phylogenetics can operate at near-real-time speed on outbreak-scale data marks a turning point in that trajectory, promising that the statistical rigor that made phylogenetics indispensable to evolutionary biology can finally be delivered at the tempo at which epidemics actually unfold.

Subject of Research: Scalable near-real-time Bayesian phylogenetic inference for infectious disease outbreak genomics

Article Title: Scalable near-real-time Bayesian phylogenetics for outbreaks with Delphy

Article References: Varilly, P., Schifferli, M., Yang, K., Cronan, P., Specht, I., Burcham, T., Glennon, O., Jacks, O., Laning, E., Marrs, L., Oba, K., Yeung, S., Zhao, K. W., Parker, E., Omah, I., Pekar, J. E., Luebbert, L., Andersen, K. G., Park, D. J., … Sabeti, P. C. (2026). Scalable near-real-time Bayesian phylogenetics for outbreaks with Delphy. Nature. https://doi.org/10.1038/s41586-026-11012-6

Image Credits: AI Generated

DOI: 10.1038/s41586-026-11012-6

Keywords: Bayesian phylogenetics, Delphy, outbreak genomics, viral evolution, SARS-CoV-2, genomic epidemiology, GPU computing, Markov chain Monte Carlo, phylodynamics, public health surveillance, time-calibrated phylogenies, real-time inference

Circulating recombinant forms, or CRFs, arise when two different HIV-1 subtypes recombine within a coinfected host, producing a mosaic genome that can then propagate through a population. CRF67_01B is one such mosaic, and its emergence illustrates a recurring problem in HIV control: once a recombinant strain establishes a foothold, its subsequent spread can escape the attention of surveillance systems that operate within provincial boundaries. Because each region sequences and reports viral genomes largely in isolation, the true geography of transmission can remain hidden until someone stitches the fragments back together. That is precisely what the Hangzhou team set out to do, drawing on viral genetic data that carry within them a record of who infected whom, when lineages diverged, and how the virus moved across the map.

The study assembled 183 CRF67_01B pol gene sequences sampled between 2011 and 2024 across six provinces. The researchers combined sequences generated locally in Hangzhou with publicly available sequences retrieved from the Los Alamos National Laboratory HIV Sequence Database, one of the largest repositories of HIV genetic data in the world. The pol gene, which encodes the virus’s essential enzymes and is routinely sequenced for drug-resistance testing, serves as a convenient and information-rich marker for molecular epidemiology. By comparing these sequences, the team could identify viruses that were so genetically similar that their carriers were likely connected directly or through a short chain of intermediate infections.

To construct those connections, the researchers built a distance-based molecular network, a technique in which each viral sequence becomes a node and pairs of sequences separated by fewer genetic differences than a chosen threshold are joined by edges. At a genetic distance threshold of 0.6 percent, the network resolved the 183 sequences into 24 putative transmission clusters, and 99 of the 183 sequences, an entry rate of 54.1 percent, fell into one of those clusters. That is a striking figure: it means that more than half of the people carrying this strain in the dataset are genetically linked to at least one other person in the network, indicating active and recent chains of transmission rather than scattered, independent introductions. The demographic pattern within those clusters was even more pronounced. Men who have sex with men accounted for 91.9 percent of the clustered individuals, confirming that this recombinant form has become firmly embedded in that community across the sampled provinces.

When the researchers examined the links that crossed provincial boundaries, one pairing stood out above all others: connections between Jiangsu Province and Hangzhou dominated the interprovincial portion of the network. This observation immediately suggested a direction of flow, with the virus moving from an established epidemic center in Jiangsu into the Zhejiang provincial capital and beyond. But a molecular network alone can only show that two viruses are related; it cannot say which lineage came first or reconstruct the historical routes the virus took. For that, the team turned to Bayesian phylogeography, a family of statistical methods that treats the geographic locations of viral ancestors as unknown variables to be inferred from the tree of genetic relationships.

The phylogeographic analysis, combined with a technique known as Bayesian stochastic search variable selection, or BSSVS, allowed the researchers to test which migration pathways between provinces were supported by the data rather than merely plausible. BSSVS works by allowing the model to switch migration rates on and off across the phylogenetic tree and retaining only the connections that measurably improve the explanation of the observed sampling locations. From this inference, the team estimated the time to the most recent common ancestor of the sampled CRF67_01B sequences, arriving at an origin in Jiangsu Province around 2007, roughly six years before the strain was first formally identified in China. The lag between true emergence and first detection is a familiar lesson of viral genomics, and it underscores how much transmission can occur before a novel recombinant attracts scientific attention.

The reconstruction also captured the epidemic’s tempo. After its origin around 2007, CRF67_01B entered a phase of exponential growth between 2009 and 2011, the classic signature of a lineage exploiting an under-saturated network of susceptible hosts. Growth then slowed, declined between 2017 and 2020, and subsequently leveled off into a plateau. That plateau does not mean the strain has been contained; rather, it suggests that the epidemic has reached a kind of equilibrium within its host population, with ongoing transmission balanced by the natural depletion of uninfected contacts, treatment-driven reductions in infectiousness, and behavioral changes. The inferred migration pathways tell a complementary story of movement: from Jiangsu into Hangzhou, from Hangzhou onward into Anhui, and from Jiangsu, Hangzhou, Anhui, and Ningbo converging into Guangdong, a province whose economic gravity draws large numbers of internal migrants.

Technically, the study demonstrates the value of combining two complementary molecular epidemiological tools. Distance-based networks are computationally inexpensive, transparent, and well suited to identifying recent transmission pairs and clusters, but they are blind to direction and time. Bayesian phylogeography is far more computationally demanding and depends on modeling assumptions, yet it recovers the temporal depth and directional flow that networks cannot. Used together, the two approaches allowed the authors to show not only that CRF67_01B sequences from Hangzhou cluster with those from Jiangsu, but that the common ancestor of the whole lineage most likely lived in Jiangsu around 2007 and that the branching order of migrations follows a coherent historical narrative matching the growth phases visible in the skyline plot of effective population size.

The public health implications are clear. An emerging recombinant strain circulating across provincial lines within a highly connected sexual network cannot be controlled by any single jurisdiction acting alone. The authors argue that their findings call for cross-regional joint prevention and control strategies, with interventions targeted at hub cities that serve as sources and conduits of viral spread, and at the high-risk MSM populations in which the strain is concentrated. In practical terms, this means coordinated partner notification across provincial boundaries, harmonized surveillance that shares sequence data in near real time, and prevention services, including testing, pre-exposure prophylaxis, and rapid treatment initiation, concentrated in the urban nodes the phylogeography identifies as transmission hubs. Hangzhou’s epidemic, on this reading, is best understood as epidemiologically linked to surrounding and earlier epicenters, and any strategy that treats it as a local problem will miss the upstream reservoirs that continually seed new infections.

The research also highlights the necessity of open genetic data. A substantial fraction of the sequences in this analysis came from a public database rather than from the authors’ own laboratory, and the cross-provincial picture would have been impossible without that shared resource. As HIV-1 continues to generate new recombinant forms worldwide, the combination of open sequence repositories, molecular network analysis, and Bayesian phylodynamic inference offers health authorities a way to see epidemics as they truly are: continuous, connected, and indifferent to administrative borders. For CRF67_01B, the genetic record now shows a lineage born in Jiangsu around 2007, exploding within the MSM community early in the following decade, and weaving an enduring web of transmission from the Yangtze delta to the south of China. The task that remains is to translate that molecular cartography into prevention that moves as fast and as far as the virus does.

Subject of Research: Spatiotemporal dynamics and cross-regional transmission of HIV-1 CRF67_01B in China

Article Title: Spatiotemporal dynamics and cross-regional transmission network of HIV-1 CRF67_01B in six provinces in China

Article References: Ye, L., Xu, K., Luo, W., Wu, S., Zhu, M., Zhang, X., & Sun, Z. (2026). Spatiotemporal dynamics and cross-regional transmission network of HIV-1 CRF67_01B in six provinces in China. Virology Journal. https://doi.org/10.1186/s12985-026-03302-2

Image Credits: AI Generated

DOI: 10.1186/s12985-026-03302-2

Keywords: HIV-1, CRF67_01B, circulating recombinant form, molecular transmission network, Bayesian phylogeography, cross-regional transmission, MSM, China, phylodynamics, viral epidemiology, Jiangsu, Hangzhou