Lumpy skin disease virus (LSDV) causes a disease of economic importance affecting cattle. Its global epidemiology is complex due to the combination of vector-borne and anthropogenic spread, the circulation of vaccine-like recombinants, and the use of vaccines. The slow molecular evolution of its DNA genome limits the utility of genetic variation for accurate tracing based on evolutionary analyses, but this limitation has not yet been formally assessed. Furthermore, until present, whole genome sequencing in affected areas has remained patchy. This study combines the first fine-grained sampling of LSDV whole genomes from a time-constrained (2015–2017) southeastern European (SEE) LSDV outbreak, which we analyze along with curated public genomes to investigate the global and regional viral dispersal dynamics. First, haplotype networks visualizing the limited genetic variability associated with the SEE LSDV outbreak show intense intermixing between countries. We also assess at which spatial scale a correlation between genetic and geographic distances can be detected for LSDV. On a global scale, we show the importance of accounting for recombination events that can impact phylogenetic and phylogeographic reconstructions. Following the assessment of the temporal signal in the recombination-free alignment, our time-scaled continuous phylogeographic analysis of Kenya-like and recent wild-type viruses confirms the origin and global dissemination history of LSDV. Our analyses highlight the importance of careful selection and application of phylodynamic approaches to DNA viruses, as well as the importance of whole genome sampling in endemic and outbreak areas to improve our understanding of the evolution, epidemiology, and transmission dynamics of DNA viruses.
IMPORTANCE
Lumpy skin disease virus (LSDV) has a complex epidemiology involving multiple strains, recombination, and vaccination. Its DNA genome provides limited genetic variation to trace outbreaks in space and time. Sequencing of LSDV whole genomes has also been patchy at global and regional scales. Here, we provide the first fine-grained whole genome sequence sampling of a constrained LSDV outbreak (southeastern Europe, 2015–2017), which we analyze along with global publicly available genomes. We formally evaluate the past occurrence of recombination events as well as the temporal signal that is required for calibrating molecular clock models and subsequently conduct a time-calibrated spatially explicit phylogeographic reconstruction. Our study further illustrates the importance of accounting for recombination events before reconstructing global and regional dynamics of DNA viruses. More LSDV whole genomes from endemic areas are needed to obtain a comprehensive understanding of global LSDV dispersal dynamics.