pith. sign in

arxiv: 2308.03714 · v1 · submitted 2023-08-07 · 🧬 q-bio.PE · cond-mat.stat-mech· q-bio.GN· q-bio.QM

The role of APOBEC3-induced mutations in the differential evolution of monkeypox virus

Pith reviewed 2026-05-24 07:57 UTC · model grok-4.3

classification 🧬 q-bio.PE cond-mat.stat-mechq-bio.GNq-bio.QM
keywords monkeypox virusAPOBEC3SNPsGTR modelmutation biasevolutionary editinglineage inference
0
0 comments X

The pith

A generalized GTR model shows recent monkeypox SNPs arise from APOBEC3 editing rather than genetic drift.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests whether the mutation pattern in newly sampled monkeypox virus genomes matches the specific editing signature produced by the human APOBEC3 enzyme or instead reflects neutral drift in animal reservoirs. It introduces a statistical test that extends the standard GTR substitution model to include an APOBEC3-biased rate matrix and compares the fit to observed single-nucleotide changes. The method extracts consistent lineage signals and counts discrete evolutionary events once the editing component is isolated. A sympathetic reader would care because the distinction decides whether the 2022 outbreak strains accumulated changes under host immune pressure or through ordinary population processes. The approach supplies a practical way to separate one known mutational process from background evolution in poxvirus data.

Core claim

We develop a simple method based on a generalization of the General-Time-Reversible (GTR) model to show that the observed SNPs are likely the result of APOBEC3-induced editing. The statistical features allow us to extract lineage information and estimate evolutionary events.

What carries the argument

A generalization of the General-Time-Reversible (GTR) model that adds an APOBEC3-specific rate matrix to separate editing signatures from neutral substitution patterns.

If this is right

  • Recent MPXV single-nucleotide changes are produced by APOBEC3 editing rather than reservoir drift.
  • Lineage relationships among outbreak strains can be recovered once the editing component is removed.
  • The number of discrete evolutionary events since the last common ancestor can be counted from the residual substitution pattern.
  • Differential evolution between early and recent MPXV clades is explained by the action of APOBEC3.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same rate-matrix extension could be used on other poxviruses or DNA viruses known to encounter APOBEC3.
  • If editing is confirmed as the dominant recent process, it implies measurable human-host pressure on circulating MPXV.
  • Future reservoir sampling campaigns could test whether the same bias appears in animal sequences collected before human spillover.

Load-bearing premise

The added APOBEC3 rate parameters correctly isolate editing effects without being confounded by unmodeled reservoir biases or other mutational processes.

What would settle it

Sequencing of additional pre-2022 reservoir isolates that display the same TC-to-TT bias in the absence of human APOBEC3 exposure would falsify the attribution.

Figures

Figures reproduced from arXiv: 2308.03714 by Otto O. Yang, Sara Habibipour, Tom Chou, Xiangting Li.

Figure 1
Figure 1. Figure 1: Schematic of the stochastic evolution model. (a) There are different potential synonymous mutation [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Simulations of different scenarios. (a-b) Four representative trajectories of an iid substitution process [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Simulations of scenarios that include mutation and selection following Eqs. 11-13. We used [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Relative molecular clocks for different types of mutations with respect to an arbitrarily chosen reference [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) The number of observed APOBEC3-relevant mutations [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The number of synonymous APOBEC3-induced mutations plotted against the collection date of the [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 4
Figure 4. Figure 4: The tight linear proportionality observed is also worth reporting and further investigation. It may reflect [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 7
Figure 7. Figure 7: Different evolutionary environments shape the shared features of pre- and post-2016 groups of MPXV [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
read the original abstract

Recent studies show that newly sampled monkeypox virus (MPXV) genomes exhibit mutations consistent with Apolipoprotein B mRNA Editing Catalytic Polypeptide-like3 (APOBEC3)-mediated editing, compared to MPXV genomes collected earlier. It is unclear whether these single nucleotide polymorphisms (SNPs) result from APOBEC3-induced editing or are a consequence of genetic drift within one or more MPXV animal reservoirs. We develop a simple method based on a generalization of the General-Time-Reversible (GTR) model to show that the observed SNPs are likely the result of APOBEC3-induced editing. The statistical features allow us to extract lineage information and estimate evolutionary events.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

4 major / 2 minor

Summary. The manuscript develops a generalization of the General-Time-Reversible (GTR) model applied to SNPs in MPXV genomes. It claims this method demonstrates that mutations observed in recently sampled human-clade sequences are attributable to APOBEC3-mediated editing rather than neutral genetic drift in unsampled animal reservoirs, while also enabling extraction of lineage information and estimation of evolutionary events.

Significance. If validated, the approach could clarify the contribution of host editing enzymes to poxvirus evolution and improve phylogenetic resolution for outbreak tracking. The paper does not report machine-checked proofs, reproducible code, or falsifiable predictions that would strengthen this assessment.

major comments (4)
  1. [Abstract / Methods] Abstract and Methods: the generalized GTR is presented as distinguishing APOBEC3 editing from drift, yet no dinucleotide-context terms (e.g., TC/GG preference) are described; without them the model reduces to a re-parameterization of the six exchangeability rates and cannot isolate the APOBEC3 signature.
  2. [Methods / Results] Methods / Results: no forward simulations under neutral drift with reservoir-specific effective population sizes or host mutation-rate heterogeneity are reported; excess C-to-T transitions could therefore be mis-attributed without this control.
  3. [Results] Results: the manuscript supplies no quantitative fit statistics, bootstrap error estimates, data-exclusion criteria, or cross-validation against held-out sequences, leaving the central attribution claim without measurable support.
  4. [Methods] Methods: the fitting procedure appears to use the same tip sequences both to estimate the generalized rate matrix and to label the resulting excess transitions as APOBEC3 evidence, creating a circularity risk that is not addressed.
minor comments (2)
  1. [Methods] Notation for the generalized rate matrix is introduced without an explicit equation; adding Eq. (X) would clarify how the six GTR parameters are extended.
  2. [Figures] Figure legends do not state the number of sequences or alignment length used, hindering reproducibility.

Simulated Author's Rebuttal

4 responses · 0 unresolved

We thank the referee for their careful reading and valuable comments on our manuscript. We address each of the major comments below and indicate the revisions we will make to strengthen the paper.

read point-by-point responses
  1. Referee: [Abstract / Methods] Abstract and Methods: the generalized GTR is presented as distinguishing APOBEC3 editing from drift, yet no dinucleotide-context terms (e.g., TC/GG preference) are described; without them the model reduces to a re-parameterization of the six exchangeability rates and cannot isolate the APOBEC3 signature.

    Authors: Our generalization of the GTR model incorporates additional parameters to capture the elevated rate of C-to-T transitions characteristic of APOBEC3 editing, allowing us to distinguish it from neutral drift by comparing to standard GTR expectations. However, we acknowledge that explicit modeling of dinucleotide contexts would provide a more precise isolation of the APOBEC3 signature. We will revise the Methods section to include a discussion of dinucleotide preferences and how our model approximates the APOBEC3 effect through the observed transition biases. revision: partial

  2. Referee: [Methods / Results] Methods / Results: no forward simulations under neutral drift with reservoir-specific effective population sizes or host mutation-rate heterogeneity are reported; excess C-to-T transitions could therefore be mis-attributed without this control.

    Authors: We agree that forward simulations would be a useful addition to rule out alternative explanations such as drift in reservoirs. Our current approach uses the deviation from the standard GTR model fitted to the data as evidence for APOBEC3, supported by the known biology of the enzyme. In the revision, we will add a section discussing potential alternative explanations and include results from simple neutral simulations to address this concern. revision: yes

  3. Referee: [Results] Results: the manuscript supplies no quantitative fit statistics, bootstrap error estimates, data-exclusion criteria, or cross-validation against held-out sequences, leaving the central attribution claim without measurable support.

    Authors: The results in the manuscript are presented through the estimated rate parameters showing the excess C-to-T. To provide stronger quantitative support, we will include likelihood ratio tests comparing the generalized model to standard GTR, bootstrap confidence intervals for the rate estimates, and specify the data exclusion criteria used. Cross-validation will be discussed as a future direction if space permits. revision: yes

  4. Referee: [Methods] Methods: the fitting procedure appears to use the same tip sequences both to estimate the generalized rate matrix and to label the resulting excess transitions as APOBEC3 evidence, creating a circularity risk that is not addressed.

    Authors: The rate matrix is estimated from the sequences, but the attribution to APOBEC3 is based on the specific pattern of increased C-to-T transitions, which matches the known biochemical preference of APOBEC3 enzymes rather than being derived solely from the fit. We will clarify this distinction in the revised Methods section to address the potential circularity concern. revision: partial

Circularity Check

0 steps flagged

No circularity detected in provided text

full rationale

The abstract describes developing a generalization of the GTR model to attribute SNPs to APOBEC3 editing, but supplies no equations, parameter definitions, or derivation steps. Without explicit model equations, fitted quantities, or self-citations that reduce the central claim to its inputs by construction, no load-bearing circular step can be exhibited. The derivation chain is therefore treated as self-contained pending the actual manuscript equations.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no explicit list of free parameters, axioms, or invented entities; the generalized GTR model is presumed to rest on standard evolutionary assumptions whose details are unavailable here.

pith-pipeline@v0.9.0 · 5662 in / 997 out tokens · 25725 ms · 2026-05-24T07:57:30.540155+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

76 extracted references · 76 canonical work pages

  1. [1]

    An APOBEC3 molecular clock to estimate the date of emergence of hMPXV

    ´Aine O’Toole and Andrew Rambaut. An APOBEC3 molecular clock to estimate the date of emergence of hMPXV. https://virological.org/t/an-apobec3-molecular-clock-to-estimate-the-date-of- emergence-of-hmpxv/885, 2022. (Accessed on 10/16/2022)

  2. [2]

    Initial observations about putative APOBEC3 deaminase editing driv- ing short-term evolution of MPXV since 2017

    ´Aine O’Toole and Andrew Rambaut. Initial observations about putative APOBEC3 deaminase editing driv- ing short-term evolution of MPXV since 2017. https://virological.org/t/initial-observations- about-putative-apobec3-deaminase-editing-driving-short-term-evolution-of-mpxv-since- 2017/830, 2022. (Accessed on 10/16/2022)

  3. [3]

    Update to observations about putative APOBEC3 deaminase editing in the light of new genomes from USA

    ´Aine O’Toole and Andrew Rambaut. Update to observations about putative APOBEC3 deaminase editing in the light of new genomes from USA. https://virological.org/t/update-to-observations-about- putative-apobec3-deaminase-editing-in-the-light-of-new-genomes-from-usa/847 , 2022. (Ac- cessed on 10/16/2022)

  4. [4]

    Vallejo, and Germ´ an Bou

    Pablo Aja-Macaya, Soraya Rumbo-Feal, Margarita Poza, Angelina Ca˜ nizares, Juan A. Vallejo, and Germ´ an Bou. A new and efficient enrichment method for metagenomic sequencing of monkeypox virus. medRxiv, 2022

  5. [5]

    Henderson

    Isao Arita and Donald A. Henderson. Monkeypox and whitepox viruses in west and central africa. Bulletin of the World Health Organization , 53(4):347, 1976

  6. [6]

    Baliere, V

    C. Baliere, V. Hourdel, A. Kwasiborski, Q. Grassin, M. Feher, D. Hoinard, J. Vanhomwegen, F. Taieb, P.-H. Consigny, J.-C. Manuguerra, I. Leclercq, C. Batejat, and V. Caro. Complete genome sequence of Monkey- pox virus from a French clinical sample and the corresponding isolated strain using nanopore sequencing. Unpublished, 2023

  7. [7]

    Barbian, S

    H. Barbian, S. Green, K. Kunstman, F. Araujo Perez, S. Bobrovska, M. Hayden, and N. Moore. Direct sub- mission. Submitted (19-DEC-2022) Regional Innovative Public Health Laboratory (RIPHL), Rush University Meical Center, 1750 W Harrison St., Chicago, Illinois 60612, USA , 2022

  8. [8]

    Rupert C. L. Beale, Svend K. Petersen-Mahrt, Ian N. Watt, Reuben S. Harris, Cristina Rada, and Michael S. Neuberger. Comparison of the differential context-dependence of DNA deamination by APOBEC enzymes: Correlation with mutation spectra in vivo. Journal of Molecular Biology , 337(3):585–596, 2004

  9. [9]

    Betancur, C.-A

    I.I.B. Betancur, C.-A. C. V. Velarde-Hoyos, R. R. G. Gomez, and M. M. R. Mercado-Reyes. Direct submis- sion. Submitted (16-DEC-2022) Antioquia, Laboratorio Departamental de Salud Publica de Antioquia, Cra 72A 141 78b, Robledo, Medellin, Antioquia 050036, Colombia , 2023

  10. [10]

    Jeff Bezanson, Alan Edelman, Stefan Karpinski, and Viral B. Shah. Julia: A fresh approach to numerical computing. SIAM Review, 59(1):65–98, 2017

  11. [11]

    Jesse D. Bloom. An experimentally informed evolutionary model improves phylogenetic fit to divergent lactamase homologs. Molecular Biology and Evolution , 31(10):2753–2769, 2014

  12. [12]

    Breman, M

    Joel G. Breman, M. V. Steniowski, E. Zanotto, A. I. Gromyko, I. Arita, et al. Human monkeypox, 1970-79. Bulletin of the World Health Organization , 58(2):165, 1980

  13. [13]

    Possible adaption of the 2022 Monkeypox virus to the human host through gene duplication and loss

    Annika Brinkmann, Claudia Kohl, Katharina Pape, Daniel Bourquain, Andrea Thuermer, Janine Michel, Lars Schaade, and Andreas Nitsche. Possible adaption of the 2022 Monkeypox virus to the human host through gene duplication and loss. bioRxiv, 2022

  14. [14]

    The modern molecular clock

    Lindell Bromham and David Penny. The modern molecular clock. Nature Reviews Genetics, 4(3):216–224, 2003

  15. [15]

    The role of altered nucleotide excision repair and UVB-induced DNA damage in melanomagenesis

    Timothy Budden and Nikola Bowden. The role of altered nucleotide excision repair and UVB-induced DNA damage in melanomagenesis. International Journal of Molecular Sciences , 14(1):1132–1151, 2013. 18

  16. [16]

    Bunge, Bernard Hoet, Liddy Chen, Florian Lienert, Heinz Weidenthaler, Lorraine R

    Eveline M. Bunge, Bernard Hoet, Liddy Chen, Florian Lienert, Heinz Weidenthaler, Lorraine R. Baer, and Robert Steffen. The changing epidemiology of human monkeypox–a potential threat? a systematic review. PLoS Neglected Tropical Diseases, 16(2):e0010141, 2022

  17. [17]

    Burns, Lela Lackey, Michael A

    Michael B. Burns, Lela Lackey, Michael A. Carpenter, Anurag Rathore, Allison M. Land, Brandon Leonard, Eric W. Refsland, Delshanee Kotandeniya, Natalia Tretyakova, Jason B. Nikas, Douglas Yee, Nuri A. Temiz, Duncan E. Donohue, Rebecca M. McDougle, William L. Brown, Emily K. Law, and Reuben S. Harris. APOBEC3b is an enzymatic source of mutation in breast c...

  18. [18]

    Christensen, Ciar´ an O’Mara, Alexey Stukalov, Kenta Sato, Jakob Nybo Nissen, Sabrina Jaye Ward, Kevin Bonham, Daniel C

    Thomas A. Christensen, Ciar´ an O’Mara, Alexey Stukalov, Kenta Sato, Jakob Nybo Nissen, Sabrina Jaye Ward, Kevin Bonham, Daniel C. Jones, Tim Holy, Morten Piibeleht, and tanhevg. Biojulia/bioalignments.jl: v3.0.0, 2022

  19. [19]

    Identification and whole-genome sequencing of a Monkeypox virus strain isolated in Israel

    Inbar Cohen-Gihon, Ofir Israeli, Ohad Shifman, Noam Erez, Sharon Melamed, Nir Paran, Adi Beth-Din, and Anat Zvi. Identification and whole-genome sequencing of a Monkeypox virus strain isolated in Israel. Microbiology Resource Announcements, 9(10):e01524–19, 2020

  20. [20]

    UVB-induced cyclobutane pyrimidine dimer frequency correlates with skin cancer mutational hotspots in p53

    Regen Drouin and Jean-Philippe Therrien. UVB-induced cyclobutane pyrimidine dimer frequency correlates with skin cancer mutational hotspots in p53. Photochemistry and Photobiology, 66(5):719–726, 1997

  21. [21]

    Evaluation of methods to detect shifts in directional selection at the genome scale

    Louis Duchemin, Vincent Lanore, Philippe Veber, and Bastien Boussau. Evaluation of methods to detect shifts in directional selection at the genome scale. Molecular Biology and Evolution , 40(2), 2022

  22. [22]

    Evolutionary trees from DNA sequences: A maximum likelihood approach

    Joseph Felsenstein. Evolutionary trees from DNA sequences: A maximum likelihood approach. Journal of Molecular Evolution, 17(6):368–376, 1981

  23. [23]

    Taking variation of evolutionary rates between sites into account in inferring phylogenies

    Joseph Felsenstein. Taking variation of evolutionary rates between sites into account in inferring phylogenies. Journal of Molecular Evolution , 53(4-5):447–455, 2001

  24. [24]

    Fitch and Charles H

    Walter M. Fitch and Charles H. Langley. Evolutionary rates in proteins: Neutral mutations and the molecular clock. In Molecular Anthropology, pages 197–219. Springer US, 1976

  25. [25]

    About monkeypox

    Centers for Disease Control and Prevention. About monkeypox. https://www.cdc.gov/poxvirus/ monkeypox/about/index.html, 2022. (Accessed on 10/16/2022)

  26. [26]

    Risk assessment: Monkeypox multi-country outbreak

    European Centre for Disease Prevention and Control. Risk assessment: Monkeypox multi-country outbreak. https://www.ecdc.europa.eu/en/publications-data/risk-assessment-monkeypox-multi- country-outbreak, 2019. (Accessed on 10/16/2022)

  27. [27]

    Foster, Edward W

    Stanley O. Foster, Edward W. Brink, Deane L. Hutchins, John M. Pifer, Bernard Lourie, Claude R. Moser, Evelyn C. Cummings, O. E. K. Kuteyi, Reginald E. A. Eke, J. B Titus, et al. Human monkeypox. Bulletin of the World Health Organization , 46(5):569, 1972

  28. [28]

    K. T. Gagnon. Direct submission. Submitted (26-SEP-2022) Molecular Biology, Microbiology, and Biochem- istry, Southern Illinois University, 1245 Lincoln Drive, Neckers Building, Rm 208, Carbondale, IL 62901, USA, 2022

  29. [29]

    Gigante, Bette Korber, Matthew H

    Crystal M. Gigante, Bette Korber, Matthew H. Seabolt, Kimberly Wilkins, Whitni Davidson, Agam K. Rao, Hui Zhao, Todd G. Smith, Christine M. Hughes, Faisal Minhaj, et al. Multiple lineages of monkeypox virus detected in the United States, 2021–2022. Science, pages 560–565, 2022

  30. [30]

    Gillespie

    Daniel T. Gillespie. Exact stochastic simulation of coupled chemical reactions. The Journal of Physical Chemistry, 81(25):2340–2361, dec 1977

  31. [31]

    Goldman and Z

    N. Goldman and Z. Yang. A codon-based model of nucleotide substitution for protein-coding DNA sequences. Molecular Biology and Evolution , 1994

  32. [32]

    Balick, Elizabeth R

    Sidhartha Goyal, Daniel J. Balick, Elizabeth R. Jerison, Richard A. Neher, Boris I. Shraiman, and Michael M. Desai. Dynamic mutation-selection balance as an evolutionary attractor. Genetics, 191(4):1309– 1319, 2012. 19

  33. [33]

    Monkeypox: epidemiology, pathogenesis, treatment and prevention

    Yong Huang, Li Mu, and Wei Wang. Monkeypox: epidemiology, pathogenesis, treatment and prevention. Signal Transduction and Targeted Therapy, 7(1), 2022

  34. [34]

    Hultquist, Joy A

    Judd F. Hultquist, Joy A. Lengyel, Eric W. Refsland, Rebecca S. LaRue, Lela Lackey, William L. Brown, and Reuben S. Harris. Human and rhesus APOBEC3d, APOBEC3f, APOBEC3g, and APOBEC3h demonstrate a conserved capacity to restrict Vif-deficient HIV-1. Journal of Virology, 85(21):11220–11234, 2011

  35. [35]

    Mahmoud R. Hussein. Ultraviolet radiation and skin cancer: molecular mechanisms. Journal of Cutaneous Pathology, 32(3):191–205, 2005

  36. [36]

    Y. J. Hutin, R. Joel Williams, Philippe Malfait, Richard Pebody, Vladamir N. Loparev, Susan L. Ropp, Mariangelli Rodriguez, Janice C. Knight, Florimont K. Tshioko, Ali S. Khan, et al. Outbreak of human monkeypox, Democratic Republic of Congo, 1996 to 1997. Emerging infectious diseases, 7(3):434, 2001

  37. [37]

    Ingvarsson

    P¨ ar K. Ingvarsson. Molecular evolution of synonymous codon usage in populus.BMC Evolutionary Biology, 8(1):307, 2008

  38. [38]

    Phylogenomic characteriza- tion and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus

    Joana Isidro, V´ ıtor Borges, Miguel Pinto, Daniel Sobral, Jo˜ ao Dourado Santos, Alexandra Nunes, Ver´ onica Mix˜ ao, Rita Ferreira, Daniela Santos, Silvia Duarte, Lu´ ıs Vieira, Maria Jos´ e Borrego, Sofia N´ uncio, Is- abel Lopes de Carvalho, Ana Pelerito, Rita Cordeiro, and Jo˜ ao Paulo Gomes. Phylogenomic characteriza- tion and signs of microevolutio...

  39. [39]

    Jukes and Charles R

    Thomas H. Jukes and Charles R. Cantor. Evolution of protein molecules. InMammalian Protein Metabolism, pages 21–132. Elsevier, 1969

  40. [40]

    Raveendran, Doss Ganesh, Sampath Kumar N

    Sudeesh Karumathil, Nimal T. Raveendran, Doss Ganesh, Sampath Kumar N. S., Rahul R. Nair, and Vijaya R. Dirisala. Evolution of synonymous codon usage bias in West African and Central African strains of Monkeypox virus. Evolutionary Bioinformatics, 14:117693431876136, jan 2018

  41. [41]

    BLAST tree: Fast filtering for genomic sequence classification

    Stuart King, Yanni Sun, James Cole, and Sakti Pramanik. BLAST tree: Fast filtering for genomic sequence classification. In 2010 IEEE International Conference on BioInformatics and BioEngineering . IEEE, 2010

  42. [42]

    Komarova, and Dominik Wodarz

    Jesse Kreger, Natalia L. Komarova, and Dominik Wodarz. A hybrid stochastic-deterministic approach to explore multiple infection and evolution in HIV. PLOS Computational Biology , 17(12):1–26, 2021

  43. [43]

    Kugelman, Sara C

    Jeffrey R. Kugelman, Sara C. Johnston, Prime M. Mulembakani, Neville Kisalu, Michael S. Lee, Galina Koroleva, Sarah E. McCarthy, Marie C. Gestole, Nathan D. Wolfe, Joseph N. Fair, et al. Genomic variability of monkeypox virus among humans, Democratic Republic of the Congo. Emerging infectious diseases , 20(2):232, 2014

  44. [44]

    Lewontin, K¯ ostas V

    Richard C. Lewontin, K¯ ostas V. Krimpas, et al. Evolutionary Genetics: From Molecules to Morphology , volume 1. Cambridge University Press, 2000

  45. [45]

    Ellsworth, Julia Krushkal, Benny H.-J

    Wen-Hsiung Li, Darrell L. Ellsworth, Julia Krushkal, Benny H.-J. Chang, and David Hewett-Emmett. Rates of nucleotide substitution in primates and rodents and the generation-time effect hypothesis. Molecular Phylogenetics and Evolution , 5(1):182–187, 1996

  46. [46]

    Lord and Michael B

    Jennifer S. Lord and Michael B. Bonsall. The evolutionary dynamics of viruses: virion release strategies, time delays and fitness minima. Virus Evolution, 7(1), 2021. veab039

  47. [47]

    Raney, Bernard B

    Jian Ma, Aakrosh Ratan, Brian J. Raney, Bernard B. Suh, Webb Miller, and David Haussler. The infinite sites model of genome evolution. Proceedings of the National Academy of Sciences , 105(38):14254–14261, 2008

  48. [48]

    Mauldin, Andrea M

    Matthew R. Mauldin, Andrea M. McCollum, Yoshinori J. Nakazawa, Anna Mandra, Erin R. Whitehouse, Whitni Davidson, Hui Zhao, Jinxin Gao, Yu Li, Jeffrey Doty, Adesola Yinka-Ogunleye, Afolabi Akinpelu, Olusola Aruna, Dhamari Naidoo, Kuiama Lewandowski, Babak Afrough, Victoria Graham, Emma Aarons, Roger Hewson, Richard Vipond, Jake Dunning, Meera Chand, Colin ...

  49. [49]

    Milhaud, M

    C. Milhaud, M. Klein, and J. J. E. A. Virat. Analyse d’un cas de variole du singe (monkeypox) chez le chimpanz´ e (pan troglogdytes).Experimentation animale, 2:121–35, 1969

  50. [50]

    Amorim, Peter F

    Priya Moorjani, Carlos Eduardo G. Amorim, Peter F. Arndt, and Molly Przeworski. Variation in the molecular clock of primates. Proceedings of the National Academy of Sciences , 113(38):10607–10612, 2016

  51. [51]

    Kosakovsky Pond, and Konrad Scheffler and

    Ben Murrell, Tulio de Oliveira, Chris Seebregts, Sergei L. Kosakovsky Pond, and Konrad Scheffler and. Mod- eling HIV-1 drug resistance as episodic directional selection. PLoS Computational Biology , 8(5):e1002507, 2012

  52. [52]

    Lash, Jinxin Gao, Hui Zhao, Yu Li, Jean-Jacques Muyembe, Placide Kingebeni, Okito Wemakoy, Jean Malekani, Kevin Karem, Inger Damon, and Darin Carroll

    Yoshinori Nakazawa, Matthew Mauldin, Ginny Emerson, Mary Reynolds, R. Lash, Jinxin Gao, Hui Zhao, Yu Li, Jean-Jacques Muyembe, Placide Kingebeni, Okito Wemakoy, Jean Malekani, Kevin Karem, Inger Damon, and Darin Carroll. A phylogeographic investigation of African Monkeypox.Viruses, 7(4):2168–2184, 2015

  53. [53]

    Richard A. Neher. Contributions of adaptation and purifying selection to SARS-CoV-2 evolution, 2022

  54. [54]

    Schmidt, Arndt von Haeseler, and Bui Quang Minh

    Lam-Tung Nguyen, Heiko A. Schmidt, Arndt von Haeseler, and Bui Quang Minh. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molecular Biology and Evo- lution, 32(1):268–274, 2014

  55. [55]

    Nnaemeka, O

    N. Nnaemeka, O. Akanbi, C. Chukwu, O. Ayansola, A. Ahmad, A. Yinka-Ogunleye, A. Akinpelu, O. Ba- batunde, C. Ihekweazu, I. Adetifa, M. R. Mauldin, C. M. Gigante, D. Batra, M. Weigand, D. Howard, M. Sheth, M. H. Seabolt, and Y. Li. 2017-2018 Nigeria Monkeypox outbreak samples. Unpublished, 2022

  56. [56]

    Evolution of molecular phenotypes under stabilizing selection

    Armita Nourmohammad, Stephan Schiffels, and Michael L¨ assig. Evolution of molecular phenotypes under stabilizing selection. Journal of Statistical Mechanics: Theory and Experiment , 2013(01):P01012, 2013

  57. [57]

    Monkeypox fact sheet

    World Health Organization. Monkeypox fact sheet. https://www.who.int/news-room/fact-sheets/ detail/monkeypox, 2022. (Accessed on 10/16/2022)

  58. [58]

    N. C. Palmateer. Direct submission. Submitted (21-NOV-2022) Public Health and Environmental Labora- tories, New Jersey Department of Health, 3 Schwarzkopf Dr, Ewing, NJ 08628, USA , 2022

  59. [59]

    Codon usage bias

    Sujatha Thankeswaran Parvathy, Varatharajalu Udayasuriyan, and Vijaipal Bhadana. Codon usage bias. Molecular Biology Reports, 49(1):539–565, 2021

  60. [60]

    Gao, and Yi Shi

    Qi Peng, Yufeng Xie, Lu Kuai, Han Wang, Jianxun Qi, George F. Gao, and Yi Shi. Structure of monkeypox virus DNA polymerase holoenzyme. Science, 379(6627):100–105, 2023

  61. [61]

    Phylogenetic models of molecular evolution: Next-generation data, fit, and performance

    David Posada. Phylogenetic models of molecular evolution: Next-generation data, fit, and performance. Journal of Molecular Evolution , 76(6):351–352, 2013

  62. [62]

    Lam, Luiz Max Carvalho, and Oliver G

    Andrew Rambaut, Tommy T. Lam, Luiz Max Carvalho, and Oliver G. Pybus. Exploring the temporal structure of heterochronous sequences using TempEst (formerly path-o-gen). Virus Evolution, 2(1):vew007, jan 2016

  63. [63]

    Rouzine, ´Eric Brunet, and Claus O

    Igor M. Rouzine, ´Eric Brunet, and Claus O. Wilke. The traveling-wave approach to asexual evolution: M¨ uller’s ratchet and speed of adaptation.Theoretical Population Biology, 73(1):24–46, 2008

  64. [64]

    Human APOBEC3 variations and viral infection

    Shiva Sadeghpour, Saeideh Khodaee, Mostafa Rahnama, Hamzeh Rahimi, and Diako Ebrahimi. Human APOBEC3 variations and viral infection. Viruses, 13(7):1366, 2021

  65. [65]

    Jaydee Sereewit, Nicole A. P. Lieberman, Hong Xie, Shah A. K. Mohamed Bakhash, B. Ethan Nunley, Benjamin Chung, Margaret G. Mills, Pavitra Roychoudhury, and Alexander L. Greninger. ORF-Interrupting Mutations in Monkeypox Virus Genomes from Washington and Ohio, 2022. Viruses, 14(11):2393, 2022. 21

  66. [66]

    Spielman and Claus O

    Stephanie J. Spielman and Claus O. Wilke. Extensively parameterized mutation-selection models reliably capture site-specific selective constraint. Molecular Biology and Evolution , 33(11):2990–3002, 2016

  67. [67]

    Stenglein, Michael B

    Mark D. Stenglein, Michael B. Burns, Ming Li, Joy Lengyel, and Reuben S. Harris. APOBEC3 proteins mediate the clearance of foreign DNA from human cells. Nature Structural & Molecular Biology, 17(2):222– 229, 2010

  68. [68]

    Some probabilistic and statistical problems in the analysis of DNA sequences

    Simon Tavar´ e. Some probabilistic and statistical problems in the analysis of DNA sequences. Lect Math Life Sci (Am Math Soc) , 17:57–86, 1986

  69. [69]

    Edward W. J. Wallace, Edoardo M. Airoldi, and D. Allan Drummond. Estimating selection on synonymous codon usage from noisy experimental data. Molecular Biology and Evolution , 30(6):1438–1453, 2013

  70. [70]

    Beatty, Mang Shi, Edward C

    Xiuwan Wang, Maura Carrai, Kate Van Brussel, Shuo Feng, Julia A. Beatty, Mang Shi, Edward C. Holmes, Jun Li, and Vanessa R. Barrs. Low intrahost and interhost genetic diversity of carnivore protoparvovirus 1 in domestic cats during a feline panleukopenia outbreak. Viruses, 14(7), 2022

  71. [71]

    Estimating the pattern of nucleotide substitution

    Ziheng Yang. Estimating the pattern of nucleotide substitution. Journal of Molecular Evolution , 39(1), jul 1994

  72. [72]

    Molecular Evolution

    Ziheng Yang. Molecular Evolution. Oxford University PressOxford, 2014

  73. [73]

    Outbreak of human monkeypox in Nigeria in 2017–18: a clinical and epidemiological report

    Adesola Yinka-Ogunleye, Olusola Aruna, Mahmood Dalhat, Dimie Ogoina, Andrea McCollum, Yahyah Disu, Ibrahim Mamadu, Afolabi Akinpelu, Adama Ahmad, Joel Burga, et al. Outbreak of human monkeypox in Nigeria in 2017–18: a clinical and epidemiological report. The Lancet Infectious Diseases , 19(8):872–879, 2019

  74. [74]

    Y.-H. You. Cyclobutane pyrimidine dimers form preferentially at the major p53 mutational hotspot in UVB-induced mouse skin tumors. Carcinogenesis, 21(11):2113–2117, 2000

  75. [75]

    Coffin, and Nathaniel R

    Qin Yu, Renate K¨ onig, Satish Pillai, Kristopher Chiles, Mary Kearney, Sarah Palmer, Douglas Richman, John M. Coffin, and Nathaniel R. Landau. Single-strand specificity of APOBEC3G accounts for minus- strand deamination of the HIV genome. Nature Structural & Molecular Biology , 11(5):435–442, 2004

  76. [76]

    other mutations

    J. Zhang. Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level. Molecular Biology and Evolution , 22(12):2472–2479, 2005. Data availability The computer code produced in this study are available at: • Analysis computer scripts: Github ( github.com/hsianktin/monkeypox). This study includes no data ...