Enterovirus A71 (EV-A71) causes hand, foot and mouth disease epidemics with neurological complications and fatalities. However, the neuropathogenesis of EV-A71 remains poorly understood. In mice, adaptation and virulence determinants have been mapped to mutations at VP2-149, VP1-145 and VP1-244. We investigate how these amino acids alter heparin-binding phenotype and shapes EV-A71 virulence in one-day old mice. We constructed six viruses with varying residues at VP1-98, VP1-145 (which are both heparin-binding determinants) and VP2-149 (based on the wild type 149K/98E/145Q, termed KEQ) to generate KKQ, KKE, KEE, IEE and IEQ variants. We demonstrated that the weak heparin-binder IEE was highly lethal in mice. The initially strong heparin-binding IEQ variant acquired an additional mutation VP1-K244E, which confers weak heparin-binding phenotype resulting in elevated viremia and increased virus antigens in mice brain, with subsequent high virulence. IEE and IEQ-244E variants inoculated into mice disseminated efficiently and displayed high viremia. Increasing polymerase fidelity and impairing recombination of IEQ attenuated the virulence, suggesting the importance of population diversity in EV-A71 pathogenesis in vivo. Combining in silico docking and deep sequencing approaches, we inferred that virus population diversity is shaped by electrostatic interactions at the five-fold axis of the virus surface. Electrostatic surface charges facilitate virus adaptation by generating poor heparin-binding variants for better in vivo dissemination in mice, likely due to reduced adsorption to heparin-rich peripheral tissues, which ultimately results in increased neurovirulence. The dynamic switching between heparin-binding and weak heparin-binding phenotype in vivo explained the neurovirulence of EV-A71.
Enterovirus A71 (EV-A71) has emerged as the most important cause of large outbreaks of severe and sometimes fatal hand, foot and mouth disease (HFMD) across the Asia-Pacific region. EV-A71 outbreaks have been associated with (sub)genogroup switches, sometimes accompanied by recombination events. Understanding EV-A71 population dynamics is therefore essential for understanding this emerging infection, and may provide pivotal information for vaccine development. Despite the public health burden of EV-A71, relatively few EV-A71 complete-genome sequences are available for analysis and from limited geographical localities. The availability of an efficient procedure for whole-genome sequencing would stimulate effort to generate more viral sequence data. Herein, we report for the first time the development of a next-generation sequencing based protocol for whole-genome sequencing of EV-A71 directly from clinical specimens. We were able to sequence viruses of subgenogroup C4 and B5, while RNA from culture materials of diverse EV-A71 subgenogroups belonging to both genogroup B and C was successfully amplified. The nature of intra-host genetic diversity was explored in 22 clinical samples, revealing 107 positions carrying minor variants (ranging from 0 to 15 variants per sample). Our analysis of EV-A71 strains sampled in 2013 showed that they all belonged to subgenogroup B5, representing the first report of this subgenogroup in Vietnam. In conclusion, we have successfully developed a high-throughput next-generation sequencing-based assay for whole-genome sequencing of EV-A71 from clinical samples.