Robert Dyrdak, Monika Mastafa, Emma Hodcroft, Richard Neher and Jan Albert
bioRxiv, vol. , 420836, 2018
Worldwide outbreaks of enterovirus D68 (EV-D68) in 2014 and 2016 have caused widespread serious respiratory and neurological disease. To investigate diversity, spread, and evolution of EV-D68 we performed near full-length deep sequencing in 54 samples obtained in Sweden during the 2014 and 2016 outbreaks. In most samples, intrapatient variability was low and dominated by rare synonymous variants, but three patients showed evidence of dual infections with distinct EV-D68 variants from the same sub-clade. Interpatient evolution showed a very strong temporal signal, with an evolutionary rate of 0.0038±0.0001 substitutions per site and year. Phylogenetic trees reconstructed from the sequences suggest that EV-D68 was introduced into Stockholm several times during the 2016 outbreak and then spread locally. Putative neutralization targets in the BC- and DE-loops of the VP1 protein are slightly more diverse within-host and tend to undergo more frequent substitution than other genomic regions. However, evolution in these loops does not appear to have driven the emergence of the 2016 B3-subclade directly from the 2014 B1-subclade. Instead, the most recent ancestor of both clades was dated to 2009. The study provides a comprehensive description of the intra- and interpatient evolution of EV-D68, including the first report of intrapatient diversity and dual infections. The new data along with publicly available EV-D68 sequences are included in an interactive phylodynamic analysis on nextstrain.org/enterovirus/d68 to facilitate timely EV-D68 tracking in the future.