The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans has been monitored at an unprecedented level due to the public health crisis, yet the stochastic dynamics underlying such a process is dubious. Here, considering the number of acquired mutations as the displacement of the viral particle from the origin, we performed biostatistical analyses from numerous whole genome sequences on the basis of a time-dependent probabilistic mathematical model. We showed that a model with a constant variant-dependent evolution rate and nonlinear mutational variance with time (i.e., anomalous diffusion) explained the SARS-CoV-2 evolutionary motion in humans during the first 120 wk of the pandemic in the United Kingdom. In particular, we found subdiffusion patterns for the Primal, Alpha, and Omicron variants but a weak superdiffusion pattern for the Delta variant. Our findings indicate that non-Brownian evolutionary motions occur in nature, thereby providing insight for viral phylodynamics.
Keywords: anomalous diffusion; dynamic systems; stochastic process; virus evolution rate; whole genome sequencing.