Examining the molecular clock hypothesis for the contemporary evolution of the rabies virus
Authors: Durrant Rowan , Cobbold Christina A. , Brunker Kirstyn , Campbell Kathryn , Dushoff Jonathan , Ferguson Elaine A. , Jaswant Gurdeep , Lugelo Ahmed , Lushasi Kennedy , Sikana Lwitiko , Hampson Katie 
Published: September 2023

The molecular clock hypothesis assumes that mutations accumulate on an organism’s genome at a constant rate over time, but this assumption does not always hold true. While modelling approaches exist to accommodate deviations from a strict molecular clock, assumptions about rate variation may not fully represent the underlying evolutionary processes. There is considerable variability in rabies virus (RABV) incubation periods, ranging from days to over a year, during which viral replication may be reduced. This prompts the question of whether modelling RABV on a per infection generation basis might be more appropriate. We investigate how variable incubation periods affect root-to-tip divergence under per-unit time and per-generation models of mutation. Additionally, we assess how well these models represent root-to-tip divergence in time-stamped RABV sequences. We find that at low substitution rates (<1 substitution per genome per generation) divergence patterns between these models are difficult to distinguish, while above this threshold differences become apparent across a range of sampling rates. Using a Tanzanian RABV dataset, we calculate the mean substitution rate to be 0.17 substitutions per genome per generation. At RABV’s substitution rate, the per-generation substitution model is unlikely to represent rabies evolution substantially differently than the molecular clock model when examining contemporary outbreaks; over enough generations for any divergence to accumulate, extreme incubation periods average out. However, measuring substitution rates per-generation holds potential in applications such as inferring transmission trees and predicting lineage emergence.

Author Summary

Rabies is a neglected disease that kills around 60,000 people each year. After entering the body, the incubation period of the virus is usually less than one month, but can sometimes span months to years. While we normally assume a virus accumulates mutations at a constant rate, it is possible that rabies’ occasional long incubation periods mean that mutations accumulate at varying rates if the virus replicates (and thus mutates) more slowly during the incubation period. We compared how the rabies virus evolves over time using two simulation models where mutations either occur per unit time or per infection generation. We also calculated the mean substitution rate per infection generation, which can be useful for inferring linkage between related rabies cases. We found that at realistic substitution rates for the rabies virus, we could not distinguish between the two models. Our calculations show that in most generations no mutations are expected to occur. Thus, over a time period long enough to observe genetic divergence, occasional long incubation periods would be “cancelled out” by shorter than average incubation periods, meaning that the two models are almost equivalent. However our work suggests that modelling substitution rates per generation may be useful for epidemiological inference.