
==== Front
Syst Biol
Syst Biol
sysbio
Systematic Biology
1063-5157
1076-836X
Oxford University Press US

38366939
10.1093/sysbio/syae003
syae003
Software for Systematics and Evolution
AcademicSubjects/SCI00960
AcademicSubjects/SCI01130
AcademicSubjects/SCI01130
Clockor2: Inferring Global and Local Strict Molecular Clocks Using Root-to-Tip Regression
https://orcid.org/0000-0002-8878-1758
Featherstone Leo A Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC 3010, Australia

https://orcid.org/0000-0003-4337-3707
Rambaut Andrew Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, UK

https://orcid.org/0000-0002-2863-0907
Duchene Sebastian Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC 3010, Australia
Department of Computational Biology, Institut Pasteur, Paris, France

https://orcid.org/0000-0001-7070-0078
Wirth Wytamma Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC 3010, Australia

Correspondence to be sent to: Leo A. Featherstone, Peter Doherty Institute for Infection and Immunity, University of Melbourne,Australia; E-mail: leo.featherstone@unimelb.edu.au
Sebastian Duchene and Wytamma Wirth supervised this work equally.

5 2024
17 2 2024
17 2 2024
73 3 623628
10 8 2023
12 1 2024
18 1 2024
02 3 2024
© The Author(s) 2024. Published by Oxford University Press, on behalf of the Society of Systematic Biologists.
2024
https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Molecular sequence data from rapidly evolving organisms are often sampled at different points in time. Sampling times can then be used for molecular clock calibration. The root-to-tip (RTT) regression is an essential tool to assess the degree to which the data behave in a clock-like fashion. Here, we introduce Clockor2, a client-side web application for conducting RTT regression. Clockor2 allows users to quickly fit local and global molecular clocks, thus handling the increasing complexity of genomic datasets that sample beyond the assumption of homogeneous host populations. Clockor2 is efficient, handling trees of up to the order of 104 tips, with significant speed increases compared with other RTT regression applications. Although clockor2 is written as a web application, all data processing happens on the client-side, meaning that data never leave the user’s computer. Clockor2 is freely available at https://clockor2.github.io/.

evolutionary rate heterogeneity
molecular clock
root-to-tip regression
Investissement d’Avenir ANR-16-CONV-0005 Australian National Health and Medical Research Council APP1157586 Australian Research Council 10.13039/501100000923 DE190100805
==== Body
pmcIntroduction

Phylodynamic analyses make use of genetic sequence data to understand the evolution, epidemiological, and ecological dynamics of a pathogen. Importantly, phyodynamics achieves its greatest value when generating insight about infectious disease dynamics beyond that offered by traditional epidemiological data. This frequently occurs at population interfaces, such as during transmission across host sub-populations, geographical boundaries, or host species. Despite the increased complexity of such datasets, the essential component to all phylodynamic modeling is the assumption of a molecular clock relating epidemiological and evolutionary timescales (Biek et al. 2015).

The simplest molecular clock model is the strict clock, which assumes a constant rate of substitution per unit time known as the “evolutionary rate” (Zuckerkandl and Pauling 1965). When the evolutionary rate is constant throughout a phylogenetic tree, the term global molecular clock is used. In contrast, a strict local clock refers to the situation where different substitution rates apply to different monophyletic groups within a tree (Ho and Duchêne 2014). The branches of local clocks are sometimes referred to as the “foreground” while the remaining branches are known as the “background,” such that there are foreground and background rates of evolution (Yoder and Yang 2000). The assumption of a local clock may, for example, correspond to sampling from different host populations, host species, or pathogen lineages (Worobey et al. 2014).

Several tools allow for the inference of strict molecular clocks via root-to-tip (RTT) regression, but none readily offers the ability to fit local clock models (Rambaut et al. 2016; Volz and Frost 2017; Hadfield et al. 2018; Sagulenko et al. 2018). Here, we introduce Clockor2, an RTT regression tool enabling rapid inference of global and local strict molecular clocks from phylogenetic trees where tips are annotated with sampling times and other relevant data. Clockor2 is bundled with an example from Dudas et al. (2018) where local clocks are fit to MERS-CoV samples from human and camel hosts, and another in the documentation from Porter et al. (2023) with SARS-CoV-2 samples from human and mink hosts.

Phylodynamic datasets are and will continue to grow in size and scope (Featherstone et al. 2022). For example, datasets of thousands to tens-of-thousands of samples have been used to understand the spread of SARS-CoV-2 at international scales, the emergence of variants of concern (VOC), and transmission between species (du Plessis et al. 2021; Hill et al. 2022; Nadeau et al. 2023; Porter et al. 2023). However, larger datasets are more likely to sample from distinct populations as a function of their size, making local clocks increasingly important. Currently, testing the fit of a local clock over alternative models, such as global or relaxed clocks, frequently requires intensive computational efforts using common Bayesian phylodynamic applications such as BEAST or RevBayes (Drummond and Rambaut 2007; Drummond et al. 2012; Höhna et al. 2016; Suchard et al. 2018; Bouckaert et al. 2019). Each generally requires hours to days of computation time. Clockor2 uniquely offers a scalable and accessible client-side web application for exploring the fit of local clocks, with results available in seconds to minutes to direct subsequent phylodynamic analysis.

Specifically, Clockor2 allows users to perform RTT regression for fitting global and local clocks (Fig. 1). The user begins by dropping or importing a rooted tree. Sampling dates and group identifiers can then be parsed from tip labels or separate files on input. Like other RTT regression applications, Clockor2 also allows users to infer the best fitting root based on the R2 value or residual mean square (RMS) of the RTT regression. Both are key indicators of clock-like evolution (Drummond et al. 2003; Rambaut et al. 2016). It also offers users a local clock-search function to explore assumptions about the number of local clocks in a dataset as well as the ability to add a local clock interactively.

Figure 1 Clockor2 presents the tree alongside RTT regression data. Users can toggle between local and global clocks and alter the appearance of the tree. In this case, a tree of MERS‑CoV samples from Dudas et al. (2018) is presented with light coloured (yellow in online version) points corresponding to samples from camel hosts and dark coloured points (dark blue in online version) corresponding to human hosts. The solid dark line in the regression panel refers to the global clock where all samples are pooled.

Methods

Documentation and Examples

The documentation for Clockor2 is available at https://clockor2.github.io/docs/. It includes demonstration of use with empirical datasets of MERS-CoV and SARS-CoV-2 from Dudas et al. (2018) and Porter et al. (2023). The MERS-CoV dataset can be loaded from the landing page and is also shown in Figure 1.

General Model for Global and Local Strict Clocks

RTT regression consists in modeling the evolutionary rate as the slope of a linear regression of the distance from the root to each tip (RTT distance), typically in units of substitutions per site (s⁢u⁢b⁢s/s⁢i⁢t⁢e), against the sampling date of each tip (Drummond et al. 2003). If we denote the evolutionary rate as r (usually in units of s⁢u⁢b⁢s/s⁢i⁢t⁢e/t⁢i⁢m⁢e), RTT distance as d (usually in units of s⁢u⁢b⁢s/s⁢i⁢t⁢e), o as the intercept (interpreted as origin), and sampling times as t, then the model for a global strict clock takes the form:

d=r⁢t+o+ϵ

where ϵ is an error term.

Clockor2 uses a generalization of this model to accommodate local clocks. For a given tree with a set of tips T, we define local clocks as pertaining to groups of tips gi and a rate parameter for each (ri). For a strict clock model with two local clocks, we then write

d={r1⁢t+o1+ϵ, if tip ∈g1r2⁢t+o2+ϵ, if tip ∈g2

We refer to groups instead of clades because while collections of tips belonging to one local clock necessarily share a common ancestor, they do not necessarily comprise a whole clade. This occurs when two or more local clocks are nested. The tips comprising the outer clock(s) then cannot comprise a whole clade if another local clock is nested within it. For example, local clock 1 shown in Figure 2b,d does not comprise a clade (i.e., is not monophyletic) because local clock 2 is nested within it.

Figure 2 Simulated examples of how local clocks appear in trees and RTT regression data. a) RTT regression data for two local clocks with similar rates separated by a long branch. b) A tree characteristic of two similar local clock rates separated by a long branch. c) RTT regression data where two local clocks have different evolutionary rates. d) A tree characteristic of two local clocks with different rates.

This general model then captures the two key scenarios where local clocks may be appropriate. The first is where rates are similar between local clocks, but separated by a long branch (Fig. 2a,b). For example, in the evolution of VOCs in SARS-CoV-2 or due to temporally sparse sampling in the case of ancient Yersisnia pestis samples (Hill et al. 2022; Tay et al. 2022; Eaton et al. 2023). The second scenario is where rates differ between local clocks (Fig. 2c,d). For example, this can occur when a pathogen spreads in different host species, such as has been observed for SARS-CoV-2 in mink and human hosts (Porter et al. 2023).

For each group of tips defining a local clock, we independently conduct an RTT regression to estimate the evolutionary rate (slope). R2 or RMS values are then an indication of clock-like behavior for each local clock. Clockor2 focuses on R2 and RMS as indicators of clock-like evolution because they offer the most straightforward interpretation of clock-like evolution. R2 values of 1 indicate perfect clock-like evolution, while values of 0 indicate a lack of a molecular clock. Likewise, lower RMS values indicate better fit of a strict clock.

Local clock and or global clock configurations can also be compared using an information criterion that combines the likelihood of each local clock’s RTT regression while penalizing the number of inferred parameters (three for each clock—slope, intercept, and variance). Clockor2 allows users to use either the Bayesian Information Cirterion (BIC), Aikake Information Criterion (AIC), or corrected Aikake Information Cirterion (AICc). We recommend using the BIC because it most heavily penalizes the addition of extra parameters, and local clocks in turn.

Derivations of the above information criteria for the local clock model are given in Supplementary material, Section 1. Briefly, these exploit the assumption of independent sampling to factor the likelihood across local clocks. Note, however, that this assumption is always flawed because samples necessarily share some ancestry by the assumption of a phylogenetic tree. In other words, ancestral branches are counted over many times in calculating the distance from root to tip for each sample (Duchêne et al. 2016). However, this is a limitation of the RTT regression approach generally, rather than of Clockor2 itself.

Local Clock Search: An Exploratory Feature

Where it is hypothesized that a dataset contains local clocks, Clockor2 provides functionality to corroborate this hypothesis by performing a search for local clocks in the tree. Briefly, the algorithm takes a maximum number of clocks and a minimum number of tips (group size) for each local clock as input parameters. It then iterates through all combinations of internal nodes from which local clocks could descend to induce corresponding local clock configurations. Importantly, the clock search algorithm tests for a number of clocks up to and including the maximum number so that it may find more parsimonious configurations with fewer clocks. Configurations are compared using the information criteria outlined above. Again, we recommend the BIC as it penalizes additional parameters (i.e., additional local clocks) most heavily. See here for an animation (https://github.com/LeoFeatherstone/clockor2Paper  /  blob / main / figures  /  clockSearchEg2Clocks.gif) of the clock search algorithm.

The clock-search algorithm operates in polynomial time (see Supplementary material, Section 2). Efficiency is improved by reducing the maximum number of local clocks in the search, increasing the minimum group size, and contingent on the topology of the underlying tree. However, the former two parameters exert a far greater effect on efficiency than topology.

Clock-Search User Guidelines

We stress that this algorithm is intended as an exploratory feature of Clockor2, rather than a formal test. It has a strong tendency to over-parameterise and select higher numbers of local clocks, even where these all have congruent evolutionary rates (see Supplementary material, Section 3 and Fig. S2). Based on this, we only suggest using the clock search if there is a biological hypothesis about why there may be a particular number of clocks. In this case, users should test up to and including the hypothesized number of clocks, but not more, to avoid the high-likelihood of over-fitting. The clock search is intended to help corroborate hypotheses about numbers of local clocks in a dataset, and we urge users to formally test these hypotheses using more rigorous methods, such as a marginal likelihood comparison if they wish to report results (Tay et al. 2023; Drummond and Suchard 2010). The Phylostems software also provides a web-based platform to explore local temporal signal in trees (Doizy et al. 2023). We point users to documentation on the current limitations of the clock search (https://clockor2.github.io/docs/examples/sars-cov-2/).

Finding the Best Fitting Root

Clockor2 selects the best fitting root based on the R2 or RMS of a global clock model for the input tree. It seeks a root minimizing among-lineage rate variation, which is equivalent to maximizing temporal signal. It follows the same algorithm as implemented in TempEst (Rambaut et al. 2016), but makes use of parallelization to improve speed for larger trees. Briefly, the tree is rooted along the branch leading to each internal node or tip, an RTT regression is performed, and the root position along the branch leading to the highest R2 or RMS is selected. When targeting R2, Clockor2 optimizes the root position using the golden-section search algorithm (Kiefer 1953). There is an analytical solution for the RMS (see Supplementary material, Section 4).

The best fitting root is inferred using a single, global clock because it presents the most parsimonious model of the evolutionary rate for a given tree. The fit of more elaborate local clock models can then be compared with this using information criteria and/or comparing the R2 or RMS values of each model. Clockor2 does not find the best fitting root for local clock models because the search space of best fitting roots and local clock configurations quickly becomes prohibitive and is possibly unidentifiable.

If a biologically informed root is available, such as with an outgroup, we suggest users retain it instead of the best fitting root. This is because the best fitting root essentially seeks to minimize among-lineage rate variation, which may contradict the biological reality.

Dependencies

Clockor2 has three key dependencies for handling, and plotting trees and RTT data. Trees are handled and manipulated using the PhyloJS (https://www.npmjs.com/package/phylojs) library. Phylocanvas is used to visualize trees and plotly.js is used to plot RTT data (Abudahab et al. 2021; Plotly-Technologies-Inc. 2015).

Results

Efficiency

Clockor2 can process trees of up to the order of 104 tips, and is thus fit for the expanding size and diversity of phylodynamic datasets. Finding the best fitting root makes use of parallelization to increase speed. Speedup is therefore proportional to the number of threads or cores available, in addition to the choice of browser and computer. For example, Clockor2 is faster than TempEst, v1.5.3, on a 2021 Macbook pro with 16 Gb of memory and 8 cores running Chrome, v113.0.5627.126 (Table 1). However, we found Clockor2 to be comparable or slower on other combinations of processor and browser, such as a Lenovo Thinkpad with an 11th Gen Intel i7 processor running Firefox v118.0.2.

Table 1 Time in seconds taken to find the best fitting root for test trees of 100, 500, 1000, 5000, and 10,000 tips using Clockor2 and TempEst v1.5.3 on a 2021 Macbook with 16 Gb of memory and 8 cores running Chrome v113.0.5627.126. Times vary with computer and browser. In general, the relative efficiency of Clockor2 will increase with the number of cores. Using the residual mean squared as an optimization target is faster because there is an analytical solution.

	R2	Residual mean squared	
Tips	Clockor2	TempEst	Clockor2	TempEst	
100	0.313	0.760	0.129	0.050	
500	1.370	2.400	0.502	1.500	
1000	3.476	10.280	1.514	5.430	
5000	78.013	272.290	24.370	122.000	
10,000	306.821	1310.340	94.992	951.000	
					

The user interface also remains responsive when working with large trees. This is in large part due to the use of WebGL in the tree and plotting components, which exploit GPU acceleration to render large and interactive trees and datasets through Phylocanvas and Plotly.js, respectively (Abudahab et al. 2021; Plotly-Technologies-Inc. 2015).

Discussion

Clockor2 provides a flexible and scalable front-end web platform for RTT regression. Its extension to fitting local clocks allows it to accommodate the growing complexity of phylodyanmic datasets as genomic epidemiology plays a growing role in infectious disease surveillance.

As a front end application, Clockor2 is also highly accessible with no installation steps required, although users have option of saving the site to run locally. Wherever there is a browser, it is possible to conduct an RTT regression using Clockor2 with the data never leaving the user’s computer. This is particularly valuable where data sharing restrictions apply, such as for patient confidentiality.

Future Directions

One future direction consists in finding the right information criterion to penalize the addition of local clocks in the clock search. This will help to transition the clock search from being an exploratory feature, to a formal test for finding local clocks. This may come in borrowing from the broader literature on continuously evolving traits (Khabbazian et al. 2016; Bastide et al. 2017). Non-parametric clustering approaches may also offer an alternative solution.

In the future, it will be possible to re-implement core functionality in increasingly popular and highly efficient programming languages that can compile to Web-Assembly format. For example, as the bioinformatics ecosystem in Rust continues to develop, it will be possible to further improve the efficiency of Clockor2 using packages such as Bio-Rust (Köster 2015).

Acknowledgments

We are grateful to Luiz Max Carvalho for helpful correspondence in ascertaining the analytical solution for the position of the best fitting root using the root mean squared. We are also grateful to Paul Bastide for helpful correspondence regarding the clock search.

Supplementary Material

Data available from the Dryad Digital Repository: http://dx.doi.org/10.5061/dryad.gxd2547sn.

Funding

This work was supported by the Inception program [Investissement d’Avenir grant ANR-16-CONV-0005], the Australian National Health and Medical Research Council [APP1157586], and the Australian Research Council [DE190100805].

Conflict of Interest

The authors declare no conflict of interest.

Data Availability

Clockor2 is available at https://clockor2.github.io/. All code required to replicate the simulation study and figures in the paper is available at https://github.com/LeoFeatherstone/clockor2Paper and https://datadryad.org / stash / dataset / doi:10.5061 / dryad.gxd2547sn. Thecode for Clockor2 is open source at https://github.com/clockor2/clockor2 and https://zenodo.org/records/10208495. No new empirical data are associated with this article.
==== Refs
References

Abudahab K. , UnderwoodA., TaylorB., YeatsC., AanensenD.M. 2021. Phylocanvas.gl: a WebGL-powered JavaScript library for large tree visualisation. Available from: URL https://osf.io/nfv6m/.
Bastide P. , MariadassouM., RobinS. 2017. Detection of adaptive shifts on phylogenies by using shifted stochastic processes on a tree. J. R. Stat. Soc. Ser. B, Methodol. 79 (4 ):1067–1093. doi:10.1111/rssb.12206
Biek R. , PybusO.G., Lloyd-SmithJ.O., DidelotX. 2015. Measurably evolving pathogens in the genomic era. Trends Ecol. Evol. 30 (6 ):306–313. doi:10.1016/j.tree.2015.03.009 25887947
Bouckaert R. , VaughanT.G., Barido-SottaniJ., DuchêneS., FourmentM., GavryushkinaA., HeledJ., JonesG., KühnertD., MaioN.D., MatschinerM., MendesF.K., MüllerN.F., OgilvieH.A., PlessisL., PopingaA., RambautA., RasmussenD., SiveroniI., SuchardM.A., WuC.-H., XieD., ZhangC., StadlerT., DrummondA.J. 2019. BEAST 2.5: An advanced software platform for bayesian evolutionary analysis. PLoS Comput. Biol. 15 (4 ):e1006650. doi:10.1371/journal.pcbi.1006650 30958812
Doizy A. , PrinA., CornuG., ChiroleuF., RieuxA. 2023. Phylostems: a new graphical tool to investigate temporal signal of heterochronous sequences datasets. Bioinform. Adv. 3 (1 ):vbad026. doi:10.1093/bioadv/vbad026 36936370
Drummond A. , PybusO.G., RambautA. 2003. Inference of viral evolutionary rates from molecular sequences. Adv. Parasitol. 54 :331–358. doi:10.1016/S0065-308X(03)54008-8 14711090
Drummond A.J. , RambautA. 2007. Beast: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7 (1 ):1–8.17214884
Drummond A.J. , SuchardM.A. 2010. Bayesian random local clocks, or one rate to rule them all. BMC Biol. 8 (1 ):114. doi:10.1186/1741-7007-8-114 20807414
Drummond A.J. , SuchardM.A., XieD., RambautA. 2012. Bayesian phylogenetics with beauti and the beast 1.7. Mol. Biol. Evol. 29 (8 ):1969–1973.22367748
du Plessis L. , McCroneJ.T., ZarebskiA.E., HillV., RuisC., GutierrezB., RaghwaniJ., AshworthJ., ColquhounR., ConnorT.R., FariaN.R., JacksonB., LomanN.J., O’TooleA., NichollsS.M., ParagK.V., ScherE., VasylyevaT.I., VolzE.M., WattsA., BogochI.I., KhanK., AanensenD.M., KraemerM.U.G., RambautA., PybusO.G.; COVID-19 Genomics UK (COG-UK) Consortium. 2021. Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK. Science. 371 (6530 ):708–712. doi:10.1126/science.abf2946 33419936
Duchêne S. , GeogheganJ.L., HolmesE.C., HoS.Y. 2016. Estimating evolutionary rates using time-structured data: a general comparison of phylogenetic methods. Bioinformatics 32 (22 ):3375–3379. doi:10.1093/bioinformatics/btw421 27412094
Dudas G. , CarvalhoL.M., RambautA., BedfordT. 2018. MERS-CoV spillover at the camel-human interface. eLife 7 :e31257. doi:10.7554/eLife.31257 29336306
Eaton K. , FeatherstoneL., DucheneS., CarmichaelA.G., VarlıkN., GoldingG.B., HolmesE.C., PoinarH.N. 2023. Plagued by a cryptic clock: insight and issues from the global phylogeny of Yersinia pestis. Commun. Biol. 6 (1 ):23.36658311
Featherstone L.A. , ZhangJ.M., VaughanT.G., DucheneS. 2022. Epidemiological inference from pathogen genomes: a review of phylodynamic models and applications. Virus Evol. 8 (1 ):veac045. doi:10.1093/ve/veac045 35775026
Hadfield J. , MegillC., BellS.M., HuddlestonJ., PotterB., CallenderC., SagulenkoP., BedfordT., NeherR.A. 2018. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. 34 (23 ):4121–4123. doi:10.1093/bioinformatics/bty407 29790939
Hill V. , Du PlessisL., PeacockT.P., AggarwalD., ColquhounR., CarabelliA.M., EllabyN., GallagherE., GrovesN., JacksonB., McCroneJ.T., O’TooleA., PriceA., SandersonT., ScherE., SouthgateJ., VolzE., BarclayW.S., BarrettJ.C., ChandM., ConnorT., GoodfellowI., GuptaR.K., HarrisonE.M., LomanN., MyersR., RobertsonD.L., PybusO.G., RambautA.; The COVID-19 Genomics UK (COG-UK) Consortium. 2022. The origins and molecular evolution of SARS-CoV-2 lineage b.1.1.7 in the UK. Virus Evol. 8 (2 ):veac080. doi:10.1093/ve/veac080 36533153
Ho S.Y. , DuchêneS. 2014. Molecular-clock methods for estimating evolutionary rates and timescales. Mol. Ecol. 23 (24 ):5947–5965.25290107
Höhna S. , LandisM.J., HeathT.A., BoussauB., LartillotN., MooreB.R., HuelsenbeckJ.P., RonquistF. 2016. RevBayes: Bayesian phylogenetic inference using graphical models and an interactive model-specification language. Syst. Biol. 65 (4 ):726–736. doi:10.1093/sysbio/syw021 27235697
Khabbazian M. , KriebelR., RoheK., AnéC. 2016. Fast and accurate detection of evolutionary shifts in Ornstein–-Uhlenbeck models. Methods Ecol. Evol. 7 (7 ):811–824. doi:10.1111/2041-210X.12534
Kiefer J. 1953. Sequntial minmax search for a maximum. Proc. Am. Math. Soc. 4 (3 ):502–506.
Köster J. 2015. Rust-Bio: a fast and safe bioinformatics library. Bioinformatics 32 (3 ):444–446. doi:10.1093/bioinformatics/btv573 26446134
Nadeau S.A. , VaughanT.G., BeckmannC., TopolskyI., ChenC., HodcroftE., SchärT., NissenI., SantacroceN., BurcklenE., FerreiraP., JablonskiK.P., Posada-CéspedesS., CapeceV., SeidelS., Santamaria de SouzaN., Martinez-GomezJ.M., ChengP., BosshardP.P., LevesqueM.P., KufnerV., SchmutzS., ZaheriM., HuberM., TrkolaA., CordeyS., LaubscherF., GonçalvesA.R., AebyS., PillonelT., JacotD., BertelliC., GreubG., LeuzingerK., StangeM., MariA., RoloffT., Seth-SmithH., HirschH.H., EgliA., RedondoM., KobelO., NoppenC., du PlessisL., BeerenwinkelN., NeherR.A., BeiselC., StadlerT. 2023. Swiss public health measures associated with reduced SARS-CoV-2 transmission using genome data. Sci. Transl. Med. 15 (680 ):eabn7979. doi:10.1126/scitranslmed.abn7979 36346321
Plotly-Technologies-Inc. 2015. Collaborative data science. Available from: URL https://plot.ly.
Porter A.F. , PurcellD.F., HowdenB.P., DucheneS. 2023. 2023. Evolutionary rate of SARS-CoV-2 increases during zoonotic infection of farmed mink. Virus Evol. 9 (1 ):vead002.36751428
Rambaut A. , LamT.T., Max CarvalhoL., PybusO.G. 2016. Exploring the temporal structure of heterochronous sequences using TempEst (formerly path-o-gen). Virus Evol. 2 (1 ):vew007. doi:10.1093/ve/vew007 27774300
Sagulenko P. , PullerV., NeherR.A. 2018. TreeTime: maximum-likelihood phylodynamic analysis. Virus Evol. 4 (1 ):vex042. doi:10.1093/ve/vex042 29340210
Suchard M.A. , LemeyP., BaeleG., AyresD.L., DrummondA.J., RambautA. 2018. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol. 4 (1 ):vey016. doi:10.1093/ve/vey016 29942656
Tay J.H. , PorterA.F., WirthW., DucheneS. 2022. The emergence of SARS-CoV-2 variants of concern is driven by acceleration of the substitution rate. Mol. Biol. Evol. 39 (2 ):msac013. doi:10.1093/molbev/msac013 35038741
Tay J.H. , BaeleG., DucheneS. 2023. Detecting episodic evolution through Bayesian inference of molecular clock models. Mol. Biol. Evol. 40 (10 ):msad212. doi:10.1093/molbev/msad212 37738550
Volz E.M. , FrostS.D.W. 2017. Scalable relaxed clock phylogenetic dating. Virus Evol. 3 (2 ):vex025. doi:10.1093/ve/vex025
Worobey M. , HanG.-Z., RambautA. 2014. A synchronized global sweep of the internal genes of modern avian influenza virus. Nature 508 (7495 ):254–257. doi:10.1038/nature13016 24531761
Yoder A.D. , YangZ. 2000. Estimation of primate speciation dates using local molecular clocks. Mol. Biol. Evol. 17 (7 ):1081–1090.10889221
Zuckerkandl E. , PaulingL. 1965. Evolutionary divergence and convergence in proteins. In: BrysonV., VogelH.J., editors. Evolving genes and proteins. New York: Academic Press. p. 97–166.
