Electrostatic interactions at the five-fold axis alter heparin-binding phenotype and drive enterovirus A71 virulence in mice

Abstract

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.

About the Authors:

Han Kang Tee

Roles Conceptualization, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Chee Wah Tan

Roles Conceptualization, Formal analysis, Investigation, Methodology, Supervision, Validation, Writing – review & editing

Affiliation: Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Thinesshwary Yogarajah

Roles Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Michelle Hui Pheng Lee

Roles Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

ORCID logo http://orcid.org/0000-0001-8897-9732

Hann Juang Chai

Roles Formal analysis, Investigation, Writing – review & editing

Affiliation: Department of Pharmacology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

ORCID logo http://orcid.org/0000-0002-4818-5564

Nur Aziah Hanapi

Roles Formal analysis, Investigation, Methodology, Writing – review & editing

Affiliation: Centre for Drug Research, Universiti Sains Malaysia, Penang, Malaysia

Siti R. Yusof

Roles Formal analysis, Investigation, Methodology, Validation, Writing – review & editing

Affiliation: Centre for Drug Research, Universiti Sains Malaysia, Penang, Malaysia

ORCID logo http://orcid.org/0000-0003-2257-9933

Kien Chai Ong

Roles Formal analysis, Investigation, Methodology, Validation, Writing – review & editing

Affiliation: Department of Biomedical Science, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

Vannajan Sanghiran Lee

Roles Formal analysis, Investigation, Methodology, Validation, Writing – review & editing

Affiliation: Department of Chemistry, Faculty of Science, University of Malaya, Kuala Lumpur, Malaysia

I-Ching Sam

Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

ORCID logo http://orcid.org/0000-0002-1296-7898

Yoke Fun Chan

Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected], [email protected]

Affiliation: Department of Medical Microbiology, Faculty of Medicine, University of Malaya, Kuala Lumpur, Malaysia

ORCID logo http://orcid.org/0000-0001-7089-0510

Introduction

Enterovirus A71 (EV-A71) causes cyclical outbreaks of hand, foot and mouth disease (HFMD) in the Asia-Pacific region [1]. HFMD primarily affects children younger than 5 years old. The clinical manifestations are usually mild and characterized by fever, oral ulcers and skin rashes on hands and feet [2, 3]. In some cases, infection also results in severe neurological complications, including encephalitis, aseptic meningitis, acute flaccid paralysis and death [4]. There are no licensed antivirals, and licensed vaccines are only available in China. The EV-A71 genome encodes a polyprotein with a single open reading frame (ORF) flanked by 5’ and 3’ untranslated regions [5]. The polyprotein is cleaved into four capsid proteins (VP1 to VP4) and seven nonstructural proteins (2A, 2B, 2C and 3A to 3D). The capsid proteins form a protomer, and five protomers form a pentamer, and then twelve pentamers assemble around a genome forming a provirion. The five-fold axis symmetry is formed by VP1 and surrounded by a canyon. Virus-receptor binding at the canyon or other physical alterations, such as heat, will displace the lipid pocket factor and trigger viral uncoating [6].

Multiple receptors including human scavenger receptor class B2 (SCARB2) [7], P-selectin glycoprotein ligand-1 (PSGL-1) [8], heparan sulfate [9], sialylated glycan [10, 11], annexin II [12], vimentin [13] and nucleolin [14] have roles in EV-A71 attachment or entry. Following infection, EV-A71 can disseminate to different organs and invade the human central nervous system (CNS) either through retrograde axonal transport or hematogenous spread [15–18]. However, the neuropathogenesis of EV-A71 remains unclear and virulence determinants are not well elucidated. Recent reports have shown that EV-A71 strains with VP1-145G/Q were more frequently isolated from severe HFMD with neurological complications [19–23]. However, the results contradicted in vivo studies which showed that VP1-145E strains but not 145G/Q exhibit higher virulence and lethality in murine models [24–27]. These findings were also consistent with studies in cynomolgus monkeys, in which strong selection of VP1-145E over VP1-145G was observed [28]. VP1-145G/Q residues but not VP1-145E are responsible for binding to the receptor PSGL-1 [29]. Both VP1-145E and VP2-149I mutations conferred Chinese Hamster Ovary (CHO) cell adaptation and increased mouse virulence [24, 26, 30].

Neurotropism and neurovirulence have been previously related to the virus binding to heparan sulfate (heparin), a negatively-charged glycosaminoglycan (GAG) found abundantly on most cell surfaces. Strong affinity for heparin was reported to cause attenuation in viruses such as Theiler's murine encephalomyelitis virus [31], Japanese encephalitis virus [32], Murray Valley encephalitis virus [32], West Nile virus [33], yellow fever virus [34] and tick-borne encephalitis virus [35]. In contrast, virus variants with weak heparin-binding resulted in higher mortality in mice, as shown for Sindbis virus [36] and eastern equine encephalitis virus (EEEV) [37]. Interestingly, variants with strong heparin-binding can contribute to higher neurovirulence if inoculated directly into the CNS, as demonstrated with EEEV [38].

EV-A71 also utilizes heparin as an attachment receptor [9]. Amino acids clustered around the five-fold symmetry axis, specifically VP1-98, VP1-145, VP1-242 and VP1-244 modulate positive charges required for heparin-binding [39]. The VP1-145G/Q residues are associated with stronger heparin-binding, but not negatively-charged VP1-145E. The propensity of EV-A71 to acquire positively-charged residues at the five-fold axis is associated with increased heparin-binding upon in vitro culture adaptation [39]. In contrast, experimental studies of VP1-145E in mice and monkey studies showed higher virulence and lethality.

In this current study, we sought to delineate the discrepancy exhibited between in vitro cytopathogenicity due to heparin-binding and association of weak heparin-binding with in vivo virulence. We demonstrated that acquisition of negative charges at the five-fold axis reduced heparin-binding. This resulted in high viremia and enhanced lethality in mice.

Results

Construction and rescue of clone-derived virus variants

We have previously described the role of VP1-98 and VP1-145 as modulators of heparin-binding in cell culture [39]. In the present study, we elucidated the role of heparin-binding in shaping EV-A71 virulence in mice. Five EV-A71 variants were engineered from a laboratory-adapted EV-A71 strain (5865/SIN/00009, subgenogroup B4), which posesses VP2-149K, VP1-98E and VP1-145Q residues (designated as KEQ, following amino acid residue order). Site-directed mutagenesis was performed to generate the five variants denoted as KKQ, KKE, KEE, IEE and IEQ (Fig 1A). Some variants had the mouse-adaptive determinant, K149I [24, 26, 40] along with substitutions in the heparin-binding determinants, VP1-98 and VP1-145. After a single passage, all EV-A71 variants were viable with comparable plaque morphology (Fig 1B) and were confirmed by sequencing (S1 Fig).

[Figure omitted. See PDF.]

Fig 1. In vitro characterization of constructed EV-A71 variants.

(A) Schematic illustration of the EV-A71 genome and the infectious clone constructs. Different amino acids were substituted at VP2-149, VP1-98 and VP1-145 (labeled in bold), with reference to the wild type (WT) strain KEQ. (B) The clone-derived EV-A71 variants were propagated in RD cells and showed comparable plaque morphologies (C) The binding affinity of EV-A71 variants to heparin sepharose beads was analyzed. (D) Inhibitory effect of heparin on EV-A71 variants was evaluated by pretreating the viruses with soluble heparin before infection of RD cells. Results are presented as mean ± SD (n = 3). Error bars indicate standard deviations from triplicates. Statistical significances are denoted with *P < 0.05, **P < 0.01 as compared to the WT.

https://doi.org/10.1371/journal.ppat.1007863.g001

Heparin-binding phenotypes of EV-A71 in in vitro cell culture

We further investigated the heparin-binding of the EV-A71 variants. EV-A71 variants with VP1-145E (KKE, KEE and IEE) displayed significant reduction of 29.7%, 18.8% and 37.1% in heparin-binding, respectively (Fig 1C). Similar findings were also observed in a heparin inhibition assay. Pre-treated heparin inhibition of these three VP1-145E variants significantly reduced RD cell viability to 51.5%, 45.5% and 34.3%, respectively (Fig 1D). Based on these heparin-binding results, the EV-A71 variants were categorized into two groups: strong heparin binder (KEQ, KKQ and IEQ) and weak heparin binder (KKE, KEE and IEE).

Potential association of weak heparin-binding variants with in vivo virulence

To determine the role of heparin-binding in EV-A71 in vivo virulence, we performed intraperitoneal (i.p.) infection of one-day old suckling mice with 1 × 105 PFU of each EV-A71 variant. The clinical score and survival analysis of infected mice are shown in Fig 2A. IEE-infected mice exhibited the highest virulence in vivo, with all the infected mice dying by day 4 post-infection (n = 10). Notably, 66.7% (n = 8) of the IEQ-infected mice died by day 12 post-infection, followed by 20% (n = 2) mortality in KEE-infected mice. Consistent with previous findings [24–26, 41], EV-A71 variants with VP1-145E were associated with an increased virulence phenotype in animal models, except for KKE. Mice infected with other EV-A71 variants showed no apparent clinical signs and survived beyond day 13 post-infection. These data suggest that in vivo pathogenicity in mice was increased by an additional VP2-K149I mutation, as demonstrated in the IEE-infected mice.

[Figure omitted. See PDF.]

Fig 2. Clinical scores and survival rates of suckling mice infected with EV-A71 variants.

One-day old suckling mice (n = 9–12) were inoculated with 1 × 105 PFU of different EV-A71 variants by i.p. injection. (A) A litter of mock-infected mice was used as a control group receiving PBS injection. The clinical scores and percentage of survival of the infected mice groups were monitored daily for 13 days. The severity of clinical symptoms was scored as follows: 0, healthy; 1, weak or less active; 2 hunched posture and lethargy; 3, one-limb paralysis; 4, two-limb paralysis; 5, moribund or dead. Significant differences compared to KEQ are labelled as * (P < 0.05) and ***P (< 0.001). KEQ, KKQ and KKE curves are identical to that of the mock-infected group. The representative VP1 sequence chromatograms of IEQ (B) and IEE (C) populations isolated from infected animal organs (n = 3) are shown, highlighting VP1-98, VP1-145 and VP1-244. The emergence of IEQ-244E virus isolated from brains and limbs of IEQ-infected moribund mice (n = 3) are shown.

https://doi.org/10.1371/journal.ppat.1007863.g002

We further asked why the IEQ variant, a strong heparin binder, exhibited a relatively high virulence in mice. Viral genomic RNA from the brains and hind limbs of dead IEQ-infected mice (n = 3) were harvested for genome sequencing, and revealed that the IEQ variant had acquired a VP1-K244E mutation (Fig 2B). This mutation however, was not present in the IEE-infected mice (Fig 2C).

Emergence of IEQ-K244E variant resulted in abolished heparin-binding and regained virulence

Since VP1-244K is a key determinant of heparin-binding [39], we speculated that the emergence of VP1-244E had abolished the heparin-binding of the IEQ variant. To investigate the role of VP1-K244E mutation in heparin-binding and virulence in mice, we introduced this mutation into the IEQ variant through site-directed mutagenesis. The IEQ variant with K244E mutation (termed IEQ-244E) however failed to achieve high virus yield in cell culture for subsequent in vivo experiments. We thus collected IEQ-244E and IEE from the brain homogenates (indicated with +) for subsequent experiments, after confirming the sequences using Sanger sequencing. IEQ-244E+ displayed significant reduction of heparin-binding compared to clone-derived IEQ and IEE, and IEE+ (Fig 3A). The IEE input was generated after one passage in RD cells, and next-generation sequencing (NGS) revealed a mixed population of 98E, 98K and 145E. In contrast, the IEE+ from brains and muscles showed dominance of 98E and 145E, which explained the weak heparin binding observed (S1 Table).

[Figure omitted. See PDF.]

Fig 3. Characterization of virulent phenotype of IEE and IEQ-244E+ variants.

(A) EV-A71 variants collected from the brain homogenates of dead infected mice were compared to clone-derived variants for heparin-binding properties. (B) One-day old suckling mice (n = 8–12) were infected with IEE, IEQ-244E+ and IEQ+ through i.p. injection. A control group receiving PBS injection was also included. The clinical scores and percentage of survival are shown over 13 days post-infection. Significant differences compared to WT are labelled as ** (P < 0.01) and ***P (< 0.001). Tissue samples of mice which succumbed to IEE (n = 2) and IEQ-244E+ (n = 2) infection were subjected to IHC and H&E staining. Virus antigen was seen in muscle cells (C), along with increased inflammatory infiltrates. In brains (D), antigen-positive neurons were seen in the pons (representative image from an IEE-infected brain) and dentate nucleus (representative image from an IEQ -244E+-infected brain), with mild inflammation in the cortex shown by H&E staining. Magnification for IHC staining: X40; H&E staining: X20.

https://doi.org/10.1371/journal.ppat.1007863.g003

To determine the association of heparin-binding and in vivo virulence, we then infected one-day old suckling mice with clone-derived IEE and IEQ-244E+ by i.p. administration. Brain homogenate from IEQ-infected surviving mice was also harvested and used as a negative control (IEQ+; viral RNA not detected in RT-PCR). At day 4 post-infection, 100% mortality was observed in IEE and IEQ-244E+-infected mice but none of the mice succumbed to IEQ+ infection (Fig 3B).

The hind limb and brain samples from IEE and IEQ-244E+-infected mice were then processed for histopathological analysis, and results supported earlier findings. Immunohistochemical (IHC) examination revealed massive localization of viral antigens in both IEE and IEQ-244E+-infected muscles, indicating that skeletal muscle is an important replication site (Fig 3C). Inflammation and extensive muscle damage were also observed in the haematoxylin and eosin (H&E)-stained sections of muscle. Mild inflammation was detected in the cortex. Viral antigens were detected in neuron bodies and their axon dendrites mainly distributed in the pons and dentate nucleus (Fig 3D). Mononuclear cell infiltrates were also evident in the cortices. In contrast, no distinctive histopathological change was observed in the mock-infected organ samples.

High lethality of weak heparin-binding EV-A71 variants correlates with high viremia

Strong heparin-binding confers the advantage of promoting virus attachment on the cell surface, thus increasing the probability of virus-functional receptor interaction in vitro [9]. However, we have demonstrated that a strong heparin-binding phenotype is deleterious to virus pathogenesis in vivo. To unravel the discrepancy of cytopathogenicity in vitro and in vivo virulence, we investigated virus dissemination in mice. Following i.p. infection, five mice were sacrificed for viral load quantitation in brain and hind limbs. IEE and IEQ-244E+ variants replicated to higher titers than IEQ in both hind limbs and brain, at day 2 and 4 post-infection (Fig 4A).

[Figure omitted. See PDF.]

Fig 4. Viral load quantitation from harvested organs and viremia level induced by EV-A71 variants.

(A) At selected time points, IEQ, IEE and IEQ-244E+-infected mice (n = 5) were sacrificed and viral loads were determined from harvested hind limbs and brains using qRT-PCR. Significant differences between viral variants are labelled as ** (P < 0.01) and *** (P < 0.001). (B) Virus clearance from blood was quantitated using qRT-PCR following intravenous inoculation of KEQ, KKE, KEE, IEE, IEQ or IEQ-244E+ into 3–4 weeks old mice. The viral RNA copies were quantitated at selected timepoints up to 30 minutes. The data are presented as Vt/V0, indicating the fraction of viral RNA copies at each timepoint (Vt) over the initial viral RNA copies (V0). Significant differences between KEQ and IEQ-244E+ are labelled as * (P < 0.05) and ** (P < 0.01).

https://doi.org/10.1371/journal.ppat.1007863.g004

To investigate the association of viremia with in vivo pathogenesis, three- to four-week-old mice were infected intravenously with KEQ, KKE, KEE, IEE, IEQ and IEQ-244E+ variants. Blood samples were then collected at 5, 15 and 30 min post-inoculation for viral load quantitation. Both strong heparin-binding variants KEQ and IEQ showed rapid viral clearance, with approximately 70% cleared from the bloodstream at 30 minutes (Fig 4B). Weak heparin-binding variants such as KKE, KEE and IEE displayed mild reduction of viremia throughout 30 minutes. Only about 23% of IEE had been cleared by 30 minutes post-infection. A sustained viremia level with minimal clearance was observed in IEQ-244E+-infected mice. Our data indicated stronger heparin-binding resulted in low viremia levels in the host, which may be due to adsorption and sequestration of viruses in surrounding tissues. The low viremia level may render the virus less efficient in disseminating to other organs, as seen in KEQ and IEQ variants.

High fidelity and impaired recombination attenuate in vivo virulence of IEQ

The emergence of a weak heparin-binding variant with VP1-244E mutation is the key determinant of in vivo adaptation and pathogenesis of IEQ. We hypothesized that IEQ is avirulent without the acquisition of the VP1-244E mutation in vivo. To reduce mutation rates and restrict generation of viral quasispecies, we engineered the viral RNA-dependent RNA polymerase (RdRp) of IEQ to harbor previously identified high-fidelity mutations G64R and L123F (abbreviated as HF in Fig 5A) [42–44] and recombination-deficient mutation Y276H (labelled as Rec-) [45]. We employed a luciferase-based replicon system to assess the impact of these mutations on genome replication. As demonstrated in Fig 5B, no significant differences in luciferase activities were observed between wild type EV-A71 Nluc Rep, EV-A71 Nluc Rep-HF and EV-A71 Nluc Rep-Rec-, suggesting that these mutated RdRp variants replicate as efficiently as the wild type. Next, we generated and rescued the IEQ-HF and IEQ-Rec- virus variants. These IEQ-HF and IEQ-Rec- variants were genetically stable with no reversion of mutations nor emergence of VP1-244E observed after a few passages, in addition to indistinguishable plaque morphology to IEQ (S2 Fig).

[Figure omitted. See PDF.]

Fig 5. High fidelity and impairing recombination attenuate in vivo virulence of IEQ.

The RdRp of IEQ was modified to harbor the high fidelity G64R and L123F mutations (HF) or the recombination-deficient Y276H mutation (Rec-) and virulence was analyzed in mice. (A) EV-A71 subgenomic replicon (EV-A71 Nluc Rep), EV-A71 Nluc Rep-HF (HF mutations are indicated by yellow triangles), EV-A71 Nluc-Rep-Rec- (Rec- mutation is denoted by orange triangle) and truncated replicon (EV-A71 Nluc Rep-del3D) were transfected separately into RD cells. (B) Luciferase activities were determined up to 8 hours post-transfection. Significant differences between replicon variants and WT are labelled as *** (P < 0. 001). (C) Virulence of IEQ-HF, IEQ-Rec- and IEQ following i.p. infection of mice was measured by clinical scores and percentage of survival. IEQ-HF and IEQ-Rec- curves are identical to that of the mock-infected PBS group. Significant difference between viral variants and IEQ is labelled as ** (P < 0.01).

https://doi.org/10.1371/journal.ppat.1007863.g005

To characterize the impact of increased fidelity and recombination deficiency on in vivo virulence, one-day old suckling mice were infected with IEQ, IEQ-HF and IEQ-Rec-. Half of the IEQ-infected mice died by day 12 post-infection (Fig 5C), while none of the mice died following IEQ-HF and IEQ-Rec- infection. Attempts to recover IEQ-HF and IEQ-Rec- viruses from all the surviving mice were unsuccessful. IEQ-HF and IEQ-Rec- did not result in any infections, and attenuated in vivo virulence of IEQ.

Emergence of VP1-244E is important for neuroinvasion

To determine if emergence of the VP1-244E mutation is critical for systemic dissemination, we next examined neurovirulence (the ability to directly infect the CNS) of all the EV-A71 variants following direct intracerebral inoculation. A dose of 1 × 105 PFU of each of the EV-A71 variants was intracerebrally injected into separate litters of one-day old mice (Fig 6). Similar mortality rates (100% mortality at day 4 post-infection) were shown by the IEE and IEQ-244E+ variants, the former having been earlier shown to be highly lethal following i.p. infection (Fig 2A). Notably, IEQ infection, which caused 66.7% mortality when inoculated intraperitoneally, now showed a drop of virulence to 8.3% (n = 1) mortality following i.c. infection. Detection of the VP1-K244E mutation by sequencing in organ samples from dead IEQ-infected mice further confirmed its importance as a neuroinvasion determinant. No virus was detected in the remaining mice with no disease symptoms. Taken together, our data indicates that the critical mutation conferring neuroinvasive phenotype to the IEQ variant is VP1-244E, which mainly arises during systemic dissemination when IEQ is inoculated intraperitoneally and not directly into the brain.

[Figure omitted. See PDF.]

Fig 6. IEQ virulence was attenuated when inoculated intracerebrally.

A dose of 1 × 105 PFU of each EV-A71 variant was administered intracerebrally into one-day old suckling mice. Clinical scores and percentage of survival over 13 days post-infection are shown. Significant differences between viral variants and WT are labelled as *** (P < 0.001).

https://doi.org/10.1371/journal.ppat.1007863.g006

Uncharged VP1-244 intermediate variants emerged during transition to IEQ-244E

To investigate how the VP1-K244E mutation emerges during in vivo infection, 14 suckling mice were intraperitoneally infected with IEQ. The mice were sacrificed at days 3, 7 and 11 post-infection to harvest hind limbs and brains for NGS of the virus population diversity. At day 9 post-infection, two moribund mice were collected. We first screened all the harvested samples using RT-PCR. At day 3 post-infection, none of the collected organs were positive for EV-A71 (Fig 7A). All five muscle samples collected were positive for EV-A71 at day 7 post-infection, suggesting that viruses were replicating in skeletal muscles. Brain samples were negative suggesting that limited virus had disseminated to the brain at this time point. As expected, both muscle and brain samples collected from the moribund mice at day 9 post-infection were positive for EV-A71. None of the surviving mice collected at day 11 post-infection were positive for EV-A71 RNA in both hind limbs and brain.

[Figure omitted. See PDF.]

Fig 7. Sequential emergence of K244E mutation involves intermediate transition variants.

One-day old suckling mice were intraperitoneally injected with 1 × 105 PFU of IEQ. (A) At days 3 and 7 post-infection, mice were euthanized to harvest their hind limb muscles and brains. At day 9 post-infection, two moribund mice were sacrificed for processing (highlighted in red). At the end of the experiment (day 11 post-infection), the remaining three healthy mice were sacrificed and categorized as ‘survivors’. The boxes indicate RT-PCR results for EV-A71 for each muscle and brain sample. A total of nine samples were positive, comprising five muscle samples from day 7 post-infection, and two muscle and brain samples from day 9 post-infection. (B) NGS results showing the frequency of different variants at the VP1-145 and VP1-244 sites. Note that 7 dpi-M5 sample had poor sequencing coverage. M indicates muscle; B indicates brain.

https://doi.org/10.1371/journal.ppat.1007863.g007

Deep sequencing of the EV-A71-positive organ samples revealed changes mainly at position VP1-145 and VP1-244 (Fig 7B). Interestingly, Q145E mutation was detected in hind leg muscle samples 7 dpi-M1 and M2 with frequency of over 80%. We observed a sequential transition of 244K to 244E from day 7 to day 9 post-infection. At day 7 post-infection, VP1-244K was predominant in two out of five hind limb samples (M1 and M2), while only a single sample (M4) showed VP1-244E as the dominant viral population. Notably, we identified a substitution, VP1-K244T from the hind limb of sample 7dpi-M3 with high frequency of over 90%. Sample 7 dpi-M5 also showed 244T and 244E, but due to poor sequencing coverage it was eliminated from the analysis. The emergence of an uncharged 244T intermediate appears to be important during the transition from positively-charged 244K to negatively-charged 244E, which is a radical change in physicochemical properties. As infection progressed, VP1-244E was solely detected in the organs harvested at day 9 post-infection, further affirming the contribution of VP1-244E to in vivo virulence. Apart from VP1-145 and 244, we also detected other non-synonymous mutations previously linked to heparin-binding from the harvested samples, including VP1-L97R, N104S and E167G (Table 1).

[Figure omitted. See PDF.]

Table 1. Non-synonymous mutations related to heparin binding detected from different organ samples of IEQ-infected mice.

https://doi.org/10.1371/journal.ppat.1007863.t001

To test if VP1-244E variant will develop specific mutation or reversion during cell culture adaptation, the IEQ-244E+ homogenates harvested from muscle and brain were propagated in RD cells. No significant CPE was observed in passage 1, and therefore the cell supernatants were collected for further propagation. At passage 2, mild CPE was observed with three out of four samples showing VP1-244E had reverted very quickly to 244K (S3 Fig). Clearly, reversion of E244K frequently emerged during in vitro culture adaptation.

Weak heparin-binding is due to loss of electrostatic interactions at the five-fold axis

IEE was experimentally proven to be highly lethal in mice. However, we observed that IEQ selectively acquires VP1-K244E over the VP1-Q145E mutation to gain neuroinvasiveness and neurovirulence (Figs 3B & 6). We reasoned that VP1-244E exhibits weaker heparin-binding ability compared to VP1-145E, and therefore, could be favorably selected in vivo. We employed in silico analysis to characterize the heparin-binding affinity of VP1-145E and VP1-244E. VP1-98, 145 and 244 are located around the five-fold axis of the EV-A71 pentamer (Fig 8A). Based on the electrostatic maps, the five-fold axis of the IEQ variant is highly positive-charged (Fig 8B), implying strong affinity to heparin. With the VP1-Q145E substitution, IEE has lower positive charges at its five-fold axis. Since the VP1-244K residue is protruding from the surface of the five-fold axis, substitution of a positively-charged lysine residue to a negatively-charged glutamic acid at this position greatly reduces the electrostatic potential on the five-fold axis and could change the capsid conformation.

[Figure omitted. See PDF.]

Fig 8. Comparison of electrostatic surface properties of the EV-A71 structure.

The analyzed amino acid residues are labelled in different colors in the EV-A71 capsid pentamer structure (PDB ID: 4AED). (A) The front views of the EV-A71 pentamer are displayed along with a magnified view of the five-fold axis using Chimera 1.10.1. (B) Electrostatic surface properties of IEQ, IEE and IEQ-244E variants were examined. Electrostatic charges are contoured from red (-5 kcal/mol·e, negatively-charged) to blue (+5 kcal/mol·e, positively-charged). Magnified views of the five-fold axis are shown for each variant. Total interaction energy (C) and individual interaction energy of residues (D) within a 4Å radius of EV-A71 variants and VP1-244 variants docked to heparin were evaluated. Total interaction energy (E) and individual interaction energy of residues (F) within a 4Å radius of other mutations detected from NGS docked to heparin were examined.

https://doi.org/10.1371/journal.ppat.1007863.g008

Simulated docking of 12-mer heparin to EV-A71 residues within a 4Å radius revealed no notable change in interaction energies between KEQ (-730.76 kcal/mol) and IEQ (-703.71 kcal/mol) (Fig 8C). The VP1-Q145E mutation of IEE resulted in a drop of 32.34 kcal/mol interaction energy when compared to KEQ. The IEQ-244E variant showed the weakest heparin-binding ability as the interaction energy drastically dropped to -476.32 kcal/mol. Compared to the IEQ-244E variant, the uncharged IEQ-244T variant exhibited higher interaction energies of -552.21 kcal/mol and -582.89 kcal/mol, respectively. The energy change was mainly contributed by the VP1-244 residue, with a descending energy order of 244K, 244T and 244E, which correlates with the scale of strong to weak heparin-binding (Fig 8D).

When examining the heparin docking of non-synonymous variants detected by NGS, the additional VP1-E167G mutation had no effect on interaction energy compared to IEQ (Fig 8E). Interestingly, variants with VP1-L97R and VP1-N104S mutations showed slight increase in the interaction energy of IEQ-244E variants, suggesting an enhanced effect of heparin-binding. Although IEQ-244E-97R showed very weak binding strength similar to IEQ-244E, a compensatory effect was seen in the 97R site which singly contributes to a more negative interaction energy, i.e. strong heparin-binding (Fig 8F). Meanwhile, the additional VP1-N104S mutation has also increased the heparin-binding of IEQ-244E variant. Taken together, different compensatory mutations could emerge to overcome the capsid instability and alter the virus fitness.

Discussion

As heparin-binding phenotype has been implicated in virulence of some neurotropic viruses, we studied the relationship between heparin-binding phenotype and mouse neurovirulence in EV-A71. Our data highlighted the key role of electrostatic interactions in shaping heparin-binding to confer virulence in mice. Among the weak heparin-binding variants used in this study (KKE, KEE and IEE), only IEE was associated with increased virulence and virus fitness in vivo. Strikingly, IEQ, which should be a strong heparin binder, showed high virulence, and we showed that this was due to the mutation VP1-244E, which conferred weak heparin-binding. This VP1-244E mutation has been previously identified as a mouse virulence determinant [46, 47]. Increasing polymerase fidelity or impairing recombination of IEQ abolishes in vivo neurovirulence possibly by restricting the emergence of VP1-244E, suggesting the importance of viral population diversity in EV-A71 pathogenesis, as reported for poliovirus [48, 49]. We showed that selection of adaptive mutations around the five-fold axis with roles in heparin-binding impacts viral pathogenicity.

Six EV-A71 variants were engineered with different amino acids at VP2-149, VP1-98 and VP1-145. VP1-98 and VP1-145 also have important roles in binding to PSGL-1 found in blood cells [8]. Both VP1-145 and VP2-149 have been implicated in mouse adaptation and virulence [24, 25, 27, 30, 46]. The VP1-98 and VP1-145 residues also act as modulators of heparin-binding in cell culture [39]. Heparin-binding may impact virulence outcome through neuroinvasion or neurovirulence [50]. Unlike SCARB2, heparin has not been showed to have roles in EV-A71 viral uncoating and internalization [7, 51–53], suggesting that heparin-dependent virulence may not to be due to direct virus-functional receptor interaction. In mice, mSCARB2 and mPSGL-1 are known to poorly support EV-A71 infection [51, 54, 55]. Therefore, an unidentified mouse receptor could be utilized to achieve high viremia and dissemination. Many reports have shown that the mutations in our engineered variants could result in better binding to mouse cells after adaptation [25, 56, 57]. As different cells may have different expression of EV-A71 receptors, tissue tropism and receptor availability will warrant further investigations. Also, KKE a poor heparin-binder, is avirulent in vivo. Heparin-binding is not solely controlled by a single mutation, and compensatory or complementary mutations may change weak heparin binders to strong binders. Other unknown receptors and additional mechanisms may also be involved in in vivo virulence.

Establishment of viremia is crucial for further dissemination to other target tissues such as skin and invasion into CNS [58]. During in vivo dissemination, the strong binding affinities of heparin binders KEQ and IEQ increase the likelihood of virus being sequestered by tissue GAG, resulting in rapid virus clearance from blood circulation [34, 36, 59]. This gives rise to a low viremia level with a substantial reduction of virulence in mice. Kobayashi and colleagues reported that KEG virus is less virulent compared to KEE due to the former harbouring the VP1-145G residue, enabling the virus to adsorb more strongly to heparin, resulting in attenuated virulence in SCARB2-expressing transgenic mice [60]. Similarly, IEE, KEE and KKE were cleared slowly in the bloodstream, and the poor heparin-binder IEQ-244E+ remained in the blood. Our findings that a strong heparin-binding phenotype attenuates virulence is also observed in other viruses, including Sindbis virus [36], Venezuelan equine encephalitis virus [59], West Nile virus [33], yellow fever virus [34] and Japanese encephalitis virus [61]. Using a monkey model, Zhang et al. showed that the establishment of viremia was strongly correlated with EV-A71 neuroinvasion into CNS [16]. A clinical study correlating prolonged viremia in EV-A71 patients with severe CNS involvement further suggests the importance of viremia in determining severity outcome [62].

We performed i.c. infection to bypass peripheral dissemination and neuroinvasion, therefore directly measuring the neurovirulence of each variant. Unlike the weak heparin binders IEE and IEQ-244E+, intracerebrally inoculated IEQ failed to exhibit the same neurovirulent phenotype as it did following i.p. administration, for which we propose two possible explanations. First, virus replication in brain cells may be restricted. This is supported by histopathological studies showing very few neurons in brain are infected with EV-A71 [63, 64]. IEQ viruses failed to acquire the VP1-K244E mutation, presumably due to lower infectivity of neurons and suboptimal replication following direct inoculation into brain, resulting in low virulence in mice. Secondly, viral multiplication in extraneural tissues such as hind limb skeletal muscle plays a key role in neuropathogenesis [65]. Mice intracerebrally inoculated with IEE showed higher viral load in hind limb muscles (S4 Fig), implying that the virus spread to peripheral tissues and underwent further extraneural replication (especially in the skeletal muscles) before re-entering the brain at high titers [16, 27, 65]. Retrograde axonal transport is the main transmission route for neuroinvasion [17, 64]. High replication will result in muscle damage that increases retrograde axonal transport and virus trafficking to the CNS [17, 66, 67]. The weak heparin binders IEE and IEQ-244E+ remain lethal since they disseminate effectively and establish high viremia prior to neuroinvasion. Using an in vitro porcine blood brain barrier (BBB) model (S1 Text), we found no correlation between heparin-binding and neuroinvasion across the BBB through tight junction leakages, as demonstrated by IEE and IEQ which induce poorer permeability compared to wild type KEQ (S5 Fig). This suggests that hematogenous spread is not the main route of spread to brain [68]; rather, establishment of high viremia appears crucial for virus dissemination to other target organs which support high levels of replication.

The mechanism of pathogenesis associated with heparin-binding is driven by electrostatic interactions at the five-fold axis of the virus surface. Alteration of the five-fold axis may affect capsid instability resulting in conformational changes which trigger genome uncoating, bypassing the need for receptor-virus binding. Mouse adaptation in poliovirus is controlled by a balance between capsid plasticity during uncoating and thermostability of the virion [69]. Similarly, electrostatic repulsion around the five-fold axis which results in capsid instability is observed in naturally thermo-labile foot-and-mouth-disease virus [70]. Dynamic switching between weak and strong heparin-binding phenotypes has been observed in vitro [39, 71]. In IEQ, VP1-145E was detected in samples harvested at early timepoints (7 dpi) but was 145Q at later timepoints (9 dpi). Multiple VP1-244 variants (244K, 244T and 244E) were generated with different effects on electrostatic interactions with heparin (K>T>E). As K to E represents a non-conservative substitution, the virus has evolved to transition through a non-charged residue (T) which has lower fitness cost during viral dissemination [72–74]. As GAGs are ubiquitous in tissues, a natural selective pressure thus exists to revert weak heparin-binding variants to heparin-binding variants which can then attach to and infect a wide range of tissues. This observed dynamic switching between weak and strong heparin-binding phenotypes highlights the importance for the virus to maintain an optimal electrostatic interaction for stable capsid conformation.

Three complementary mutations, VP1-L97R, VP1-N104S and VP1-E167G, were detected along with VP1-K244E from IEQ-infected moribund mice. Given that these mutations are near the VP1-244 residue site and located at loops (S6A Fig), they could be selectively utilized by the virus to stabilize the conformational structure of the VP1-K244E mutation at the five-fold axis. The VP1-N104S mutation was frequently associated with the VP1-244 variants. Root mean square fluctuation (RMSF) was measured to examine the dynamic movement of residues (S6B Fig); the RMSF value of the BC loop (where VP1-104 is located) was relatively higher in IEQ-244E compared to other virus variants. This indicates that residues within the VP1 BC loop are flexible and may bind poorly to heparin. The emergence of the VP1-N104S mutation could contribute to the conformational stability of VP1-K244E, resulting in enhanced viral infectivity [75]. The VP1-L97R mutation, within the VP1 BC loop, has been reported to enhance heparin-binding [76]. This mutation was first detected from EV-A71 in the blood, CSF and stool from an immunocompromised patient [77]. However, we found this mutation in both the skeletal muscle and brain of an IEQ-infected moribund mouse. Similarly, the presence of the VP1-L97R could stabilize the conformational structure of VP1-244E. Interestingly, both VP1-97R and 104S did not co-exist in the same sample. The VP1-D31G [78] and VP1-E167G [77] mutations previously implicated in neurotropism in humans were also observed in the brain samples but at very low frequencies. Furthermore, restricting polymerase fidelity and impairing recombination rendered the virus avirulent, suggesting that restriction of population diversity alters fitness of the virus in vivo. The transition of different virus subpopulations is pivotal for the virus to gain fitness and adaptation in vivo.

We speculate that our findings are not only limited to the murine model but may reflect neuropathogenesis in humans. Strong heparin-binding variants (KEQ and EGK) have been more frequently detected from sequencing of virus cultures than from direct sequencing of clinical specimens, suggesting that heparin-binding phenotypes are a consequence of adaptation to cell culture (S2 Table, S1 Appendix). Detection of VP1-145E in the sequences of a fatal encephalitis autopsy specimen [79] further reinforces our view that weak heparin-binding is associated with virulence in humans, as we have shown in mice in the present study. EV-A71 infection is usually mild and limited to HFMD, with neurological complications seen in 0.1–1.1% and deaths in 0.01% to 0.03% [4, 80–83]. The frequent reversion between heparin-binding and weak heparin-binding variants fulfils the trade-off hypothesis in which the virus juggles between the virulence and transmission. IEE and IEQ-244E are completely lethal whereas the IEQ variant demonstrated intermediate virulence in mice.

We proposed a hypothetical EV-A71 pathogenesis model to show the importance of heparin-binding in human infection (Fig 9). Three determinant factors of EV-A71 virulence are virus entry, dissemination and neuroinvasion. Both strong and weak heparin binders could infect humans, probably at different susceptibility. Primary viremia is established upon virus entry in primary replication sites such as tonsils and oropharynx [15]. Strong heparin binders are absorbed to GAG-rich tissues and readily removed from blood circulation due to their high affinity to heparin, thereby reducing the viremia level. Weak heparin binders are not adsorbed into tissues, but remain in the bloodstream, undergoing further extraneural replication in skeletal muscles, giving rise to high viremia. Upon overcoming the immune system, the high viremia results in better replication and dissemination to skeletal muscles, and peripheral motor nerves, through which the virus invades the CNS by retrograde axonal transport [15, 68]. High viremia alone does not result in direct CNS invasion as the virus cannot traverse the BBB. Many research questions however remain. What additional host factors could alter the heparin-binding phenotypes? Could the host immune responses influence the fate of the viruses if activated early enough before the high viremia stage?

[Figure omitted. See PDF.]

Fig 9. Hypothesized model of EV-A71 heparin-dependent pathogenesis in human.

Three major factors are responsible for EV-A71 virulence determination, namely virus entry, peripheral dissemination and neuroinvasion. (A) Both strong and weak heparin binders infect humans at the same rate, using the same inoculation route and receptor. (B) Viremia is established upon virus entry. Strong heparin binders are more readily removed from the blood circulation by binding to peripheral tissues due to their high affinity to heparin. Meanwhile, weak heparin binders give rise to higher viremia with better dissemination to other organs. (C) Neuroinvasion occurs when virus travels from peripheral motor nerves to the CNS via retrograde axonal transport.

https://doi.org/10.1371/journal.ppat.1007863.g009

Our current results have several limitations. We used i.p. and i.c. routes, and both are not the natural route of infection in humans. It would be interesting to use the oral route and examine the virulence in mice in future. The animal models used here were 1-day old suckling mice lacking of mature immune system, and both hSCARB2 and hPSGL-1 receptors. The use of hybrid mouse models such as hSCARB2/mSCARB2, hSCARB2 with stat-1 (interferon signalling) knock-out can provide efficient virus replication and mimics human pathogenesis [84, 85]. The use of original clinical isolates without cell culture propagation for in vivo infection may help to mimic infection in human. Also in the study, only VP1 was sequenced, and there may be other genomic changes within the virus genome.

In summary, we showed that weak heparin-binding EV-A71 is highly virulent in mice, in contrast with strong heparin binders which show higher replication in vitro due to culture adaptation. This study shows that weak heparin-binding EV-A71 is preferentially selected to disseminate via the bloodstream; in contrast, strong heparin-binding EV-A71 is adsorbed to peripheral tissues and rapidly cleared. The electrostatic surface charges at the VP1 capsid shape heparin-binding and hence EV-A71 virulence. Our findings provide the mechanistic action of heparin-dependent virulence, and have potential therapeutic implications for viruses which utilize heparin as an attachment receptor and are dependent on high viremia levels to cause infection.

Materials and methods

Ethics statement

The animal experiments were carried out in accordance with the rules and guidelines of the Animal Experimental Unit (AEU) in University of Malaya. The protocols were reviewed and approved by the Institutional Animal Care and Use Committee of the Faculty of Medicine, University of Malaya (reference number: 2016-190908/R/TCW).

Cell lines and viruses

Human rhabdomyosarcoma (RD, ATCC no.: CCL-136) was propagated in Dulbecco’s Modified Eagle’s Medium (DMEM) (Life Technologies). All cells were supplemented with 10% fetal bovine serum (FBS). Infected cells were maintained in media containing 2% FBS. All cells were maintained at 37°C in 5% CO2.

EV-A71 strain 41 (5865/SIN/000009, GenBank accession no. AF316321; subgenogroup B4) was used for construction of infectious clones using a DNA-launched strategy as reported previously [86]. Unlike its original sequence, this lab strain had been propagated in tissue culture previously and had acquired a VP1-E145Q mutation. We have denoted this strain as KEQ.

Different mutations were incorporated into the EV-A71 infectious clone plasmid using Q5 high-fidelity DNA polymerase (NEB) PCR site-directed mutagenesis with primers listed in S3 Table. The purified PCR products were treated with T4 polynucleotide kinase, T4 ligase and DpnI (NEB) for 1 hour at room temperature. The ligation mixture was then transformed into E. coli XL-10 GOLD ultracompetent cells (Agilent Technologies). The plasmids were transfected into RD cells using Lipofectamine LTX (life Technologies) as reported previously [39]. The transfected virus stock (P0) was further propagated in RD to generate P1 stock for all experiments.

Binding of EV-A71 particles to immobilized heparin sepharose beads

A binding assay was performed using columns with immobilized heparin sepharose beads as previously reported [39]. In brief, 200 μl of heparin sepharose (Abcam, UK) was aliquoted into a Pierce Spin Cup with cellulose acetate filter (Thermo Scientific, USA). The heparin sepharose beads were washed twice with binding buffer (0.02 M Tris-HCl, 0.14 M NaCl, pH 7.4), before addition of each virus variant (1 × 105 PFU in 600 μl). The columns containing viruses were incubated for 30 minutes at 4°C, and this was followed by centrifugation and 5 washing steps. The heparin-bound viruses were collected after eluted with elution buffer (0.02 M Tris-HCl, 2M NaCl, pH 7.4). Both virus input and output fractions were quantitated using real-time PCR, and the heparin binding efficiency was normalized by dividing the output viral RNA copy number over the input viral RNA copy number.

Evaluation of inhibitory effect of soluble heparin on EV-A71 variants

To determine the inhibitory effect of soluble heparin on EV-A71 variants, a virus inactivation assay was performed as previously described [9]. In brief, viruses were incubated with 2.5 mg/ml of soluble heparin (Sigma, USA) for an hour at 37°C. The treated viruses were inoculated onto pre-seeded RD cells and incubated at 37°C. Two days later, the cell viability of each infected virus variants was measured using CellTiter 96 Aqueous One solution Cell Proliferation Assay (Promega, USA). The relative cell viability was calculated with the following formula:

Mice infection experiments

Groups of one-day old ICR suckling mice (n = 9 to 12) were obtained from AEU. Each group of suckling mice were either intraperitoneally or intracerebrally inoculated with 1 × 105 PFU of each EV-A71 variant or PBS alone. All infected mice were monitored daily for weight change and health status up to 13 days post-infection. A clinical score was recorded using the following grades: 0, healthy; 1, weak or less active; 2, hunched posture and lethargy; 3, one-limb paralysis; 4, two-limb paralysis; 5, moribund or dead. Moribund mice were sacrificed and removed along with any mice found dead. Harvested mice organs were homogenized using hard tissue homogenizing mix (Omni International, USA). RNA was extracted from the homogenates with QIAamp viral RNA mini kit (Qiagen, Germany). The viral loads in organs were determined using TaqMan fast virus 1-step master mix (ABI, USA). One step RT-PCR was also performed to amplify viral RNA from the organs using MyTaq One-Step RT-PCR kit (Bioline, UK) for Sanger sequencing or deep sequencing. Illumina Miseq (Illumina, USA) was performed with 150 nucleotides and 250 nucleotides paired end reads. All listed mutations had an average coverage of at least 20,000 reads unless otherwise stated. The NGS reads were analyzed using CLC Bio Genomic Workbench (Qiagen) and Geneious Prime (Biomatters Ltd, New Zealand). Only variants with frequency >1% were reported.

For the virus clearance assay, 3- to 4 weeks old ICR mice (weighing within 25-35g) were i.p. injected with ketamine/xylazine cocktail prior to infection. Anesthetized animals were then intravenously inoculated with 5 × 105 PFU of each EV-A71 variant via the tail vein. At certain timepoints, blood was collected from anaesthetized mice through the retro-orbital plexus with the use of a sodium heparinized hematocrit capillary (Hirsschmann, Germany). The collected whole blood was then used for viral RNA quantitation.

Immunohistochemistry

Immunohistochemistry (IHC) were performed by the standard ENVISION technique as described previously [87]. Briefly, deparaffinised and rehydrated tissue sections were blocked using standard immunoperoxidase procedure before antigen retrieval (30 minutes, 99°C, Tris EDTA buffer with 0.05% Tween-20). Tissues were then incubated with rabbit polyclonal EV-A71 VP1 (GeneTex, USA) at 4°C overnight. After washing, tissues were then incubated with goat-anti rabbit HRP-conjugate (Dako, Denmark) for 30 minutes at room temperature. Tissues were stained using DAB (Dako) and counterstained with hematoxylin (Dako). The tissues were mounted using DPX (Dako) prior to examination under a light microscope. The negative control tissues for IHC included mock-infected ICR mice brains (n = 2) and hind-limb muscle tissues (n = 2). Isotype control antibodies or normal rabbit immunoglobulin fractions (Dako) were also used to exclude non-specific staining.

Electrostatic surface charge analysis of EV-A71 structure

The EV-A71 structure was visualized using Chimera software (UCSF Chimera version 1.10.1, USA). Electrostatic surface potentials of virus capsid were analyzed using the ‘Coulombic surface coloring’ function in which the capsid residues were labelled with different colors based on their electrostatic charges. Positively-charged residues were colored blue while negatively-charged residues were colored red.

Molecular docking simulation of EV-A71 VP1 and 12-mer heparin

Molecular docking simulation of EV-A71 crystal structure (PDB ID: 4AED) was performed using CDOCKER (CHARMm-based DOCKER) [88], as previously described [39].

Statistics

All experiments were performed with at least two biological duplicates. Data are shown with error bars indicating standard deviations. Student’s t-test was performed for all in vitro experiments as well as viral load quantitation from mice organs. Survival of mice was evaluated using Kaplan-Meier analysis. GraphPad Prism version 5.03 (GraphPad Software, USA) was used for statistical analyses with a P value of < 0.05 indicating significance.

Supporting information

S1 Fig. EV-A71 variants were confirmed by Sanger sequencing.

https://doi.org/10.1371/journal.ppat.1007863.s001

(TIF)

S2 Fig. Comparison of plaque morphologies for IEQ, IEQ-HF and IEQ Rec-.

https://doi.org/10.1371/journal.ppat.1007863.s002

(TIF)

S3 Fig. Reversion of IEQ-244E+ variant sequences after propagation in RD cells.

A sequential reversion of VP1-244E to 244K was observed in IEQ-244E+ after cell culture propagation. Note that the 244T variant was also present in the mixed population. The mixed populations of K244 (red), T244 (pink) and E244 (maroon) are shown in different proportions. Only variants with frequency >1% are reported. Note that samples B1 at P0 and B2 at P0 had poor sequencing coverages (between 110–5100). M indicates muscle, B indicates brain.

https://doi.org/10.1371/journal.ppat.1007863.s003

(TIF)

S4 Fig. Comparison of viral loads from muscles and brain in mice infected through the i.c. route.

One-day old suckling mice (n = 3) were infected with IEE through i.c. route of administration. At day 4 post-infection, muscles and brains were harvested and viral loads were quantitated using qRT-PCR. Significant comparisons are labelled *** (P < 0.001).

https://doi.org/10.1371/journal.ppat.1007863.s004

(TIF)

S5 Fig. Effects of EV-A71 exposure on a porcine in vitro BBB model.

Illustration of the porcine in vitro BBB model which simulates the movement of virus particles through an in vivo BBB, in which the luminal side represents the blood capillary while the abluminal compartment represents the brain (A). The in vitro model was exposed to different EV-A71 variants with titer of 1 × 105 PFU. The BBB permeability induced by EV-A71 variants were assessed in terms of transendothelial electrical resistance (TEER), with a greater reduction of TEER indicating greater permeability of BBB through tight junction leakages. The TEER was recorded at 2 hours (B) and 6 hours post-exposure (C) along with non-infected cell controls (white bars) and normalized with TEER values measured before virus exposure. Results are presented as mean ± SD (n = 6). Significant differences between viral variants and WT (black bars) are labelled as * (P < 0.05) and **(P < 0.01), using the Student’s t test.

https://doi.org/10.1371/journal.ppat.1007863.s005

(TIF)

S6 Fig. Structural modelling of IEQ-244E variant and root mean square fluctuation analysis of different EV-A71 variants.

(A) Structural modelling of VP1 amino acid residues of IEQ-244E (left panel). Each important amino acid is labelled with different colors: VP1-244E in red, VP1-244K in green, VP1-97L in orange and VP1-104N in magenta. Note that VP1-145E is not visible from this angle. The surface of IEQ-244E (right panel) is displayed corresponding to the structural model. (B) Root mean square fluctuation (RMSF) value of VP1 and VP2 amino acids are displayed for different variants. VP1 comprises residues 1–297 whereas VP2 consists of residues 298–542. BC, EF and GH loops of VP1 are labelled accordingly.

https://doi.org/10.1371/journal.ppat.1007863.s006

(TIF)

S1 Table. Non-synonymous mutations related to heparin binding detected from different organ samples of IEE-infected mice.

https://doi.org/10.1371/journal.ppat.1007863.s007

(DOCX)

S2 Table. Comparison of EV-A71 isolate sequences of primary specimens and passaged isolates.

Strong heparin binders (denoted with asterisks) were more frequently identified from sequencing of passaged EV-A71 (at least one passage) than from direct sequencing of primary specimens, suggesting that the virus isolates have undergone heparin-binding adaptation in cell culture (P = 0.00012, chi-square test).

https://doi.org/10.1371/journal.ppat.1007863.s008

(DOCX)

S3 Table. Primer sequences used for RT-PCR and qRT-PCR.

Primer sets used for EV-A71 VP1 sequencing and qRT-PCR are shown.

https://doi.org/10.1371/journal.ppat.1007863.s009

(DOCX)

S1 Appendix. Details of published sequences obtained from direct sequencing of clinical specimens and sequencing from cell cultures.

https://doi.org/10.1371/journal.ppat.1007863.s010

(PDF)

S1 Text. Establishment of in vitro blood-brain barrier model.

https://doi.org/10.1371/journal.ppat.1007863.s011

(DOCX)

Citation: Tee HK, Tan CW, Yogarajah T, Lee MHP, Chai HJ, Hanapi NA, et al. (2019) Electrostatic interactions at the five-fold axis alter heparin-binding phenotype and drive enterovirus A71 virulence in mice. PLoS Pathog 15(11): e1007863. https://doi.org/10.1371/journal.ppat.1007863

References

1. NikNadia N, Sam IC, Rampal S, WanNorAmalina W, NurAtifah G, Verasahib K, et al. Cyclical patterns of hand, foot and mouth disease caused by Enterovirus A71 in Malaysia. PLoS Negl Trop Dis. 2016;10(3):e0004562. Epub 2016/03/25. pmid:27010319.

2. Schmidt NJ, Lennette EH, Ho HH. An apparently new enterovirus isolated from patients with disease of the central nervous system. J Infect Dis. 1974;129(3):304–9. Epub 1974/03/01. pmid:4361245.

3. Ooi MH, Wong SC, Lewthwaite P, Cardosa MJ, Solomon T. Clinical features, diagnosis, and management of enterovirus 71. The Lancet Neurology. 2010;9(11):1097–105. Epub 2010/10/23. pmid:20965438.

4. McMinn PC. An overview of the evolution of enterovirus 71 and its clinical and public health significance. FEMS microbiology reviews. 2002;26(1):91–107. Epub 2002/05/15. pmid:12007645.

5. Brown BA, Pallansch MA. Complete nucleotide sequence of enterovirus 71 is distinct from poliovirus. Virus Res. 1995;39(2–3):195–205. Epub 1995/12/01. pmid:8837884.

6. Shingler KL, Yoder JL, Carnegie MS, Ashley RE, Makhov AM, Conway JF, et al. The enterovirus 71 A-particle forms a gateway to allow genome release: a CryoEM study of picornavirus uncoating. PLoS Pathog. 2013;9(3):e1003240. pmid:23555253

7. Yamayoshi S, Yamashita Y, Li J, Hanagata N, Minowa T, Takemura T, et al. Scavenger receptor B2 is a cellular receptor for enterovirus 71. Nat Med. 2009;15(7):798–801. Epub 2009/06/23. pmid:19543282.

8. Nishimura Y, Shimojima M, Tano Y, Miyamura T, Wakita T, Shimizu H. Human P-selectin glycoprotein ligand-1 is a functional receptor for enterovirus 71. Nat Med. 2009;15(7):794–7. Epub 2009/06/23. pmid:19543284.

9. Tan CW, Poh CL, Sam IC, Chan YF. Enterovirus 71 uses cell surface heparan sulfate glycosaminoglycan as an attachment receptor. J Virol. 2013;87(1):611–20. Epub 2012/10/26. pmid:23097443.

10. Yang F, Ren L, Xiong Z, Li J, Xiao Y, Zhao R, et al. Enterovirus 71 outbreak in the People's Republic of China in 2008. J Clin Microbiol. 2009;47(7):2351–2. Epub 2009/05/15. pmid:19439545.

11. Su PY, Liu YT, Chang HY, Huang SW, Wang YF, Yu CK, et al. Cell surface sialylation affects binding of enterovirus 71 to rhabdomyosarcoma and neuroblastoma cells. BMC Microbiol. 2012;12:162. Epub 2012/08/03. pmid:22853823.

12. Yang B, Solakyildirim K, Chang Y, Linhardt RJ. Hyphenated techniques for the analysis of heparin and heparan sulfate. Anal Bioanal Chem. 2011;399(2):541–57. Epub 2010/09/21. pmid:20853165.

13. Du N, Cong H, Tian H, Zhang H, Zhang W, Song L, et al. Cell surface vimentin is an attachment receptor for enterovirus 71. J Virol. 2014;88(10):5816–33. Epub 2014/03/14. pmid:24623428.

14. Su PY, Wang YF, Huang SW, Lo YC, Wang YH, Wu SR, et al. Cell surface nucleolin facilitates enterovirus 71 binding and infection. J Virol. 2015;89(8):4527–38. Epub 2015/02/13. pmid:25673703.

15. Ong KC, Wong KT. Understanding enterovirus 71 neuropathogenesis and its impact on other neurotropic enteroviruses. Brain Pathol. 2015;25(5):614–24. Epub 2015/08/16. pmid:26276025.

16. Zhang Y, Cui W, Liu L, Wang J, Zhao H, Liao Y, et al. Pathogenesis study of enterovirus 71 infection in rhesus monkeys. Lab Invest. 2011;91(9):1337–50. Epub 2011/05/11. pmid:21555996.

17. Chen CS, Yao YC, Lin SC, Lee YP, Wang YF, Wang JR, et al. Retrograde Axonal Transport: a Major Transmission Route of Enterovirus 71 in Mice. J Virol 2007;81(17):8996–9003. pmid:17567704

18. Feng M, Guo S, Fan S, Zeng X, Zhang Y, Liao Y, et al. The Preferential Infection of Astrocytes by Enterovirus 71 Plays a Key Role in the Viral Neurogenic Pathogenesis. Front Cell Infect Mi. 2016;6:192-. pmid:28066727.

19. Li Y, Zhu R, Qian Y, Deng J, Sun Y, Liu L, et al. Comparing Enterovirus 71 with Coxsackievirus A16 by analyzing nucleotide sequences and antigenicity of recombinant proteins of VP1s and VP4s. BMC Microbiol. 2011;11:246. Epub 2011/11/05. pmid:22050722.

20. Chang SC, Li WC, Chen GW, Tsao KC, Huang CG, Huang YC, et al. Genetic characterization of enterovirus 71 isolated from patients with severe disease by comparative analysis of complete genomes. J Med Virol. 2012;84(6):931–9. Epub 2012/04/14. pmid:22499017.

21. Zhang B, Wu X, Huang K, Li L, Zheng L, Wan C, et al. The variations of VP1 protein might be associated with nervous system symptoms caused by enterovirus 71 infection. BMC Infect Dis. 2014;14:243. Epub 2014/06/03. pmid:24886383.

22. Liu Y, Fu C, Wu S, Chen X, Shi Y, Zhou B, et al. A novel finding for enterovirus virulence from the capsid protein VP1 of EV71 circulating in mainland China. Virus Genes. 2014;48(2):260–72. Epub 2014/01/21. pmid:24442718.

23. Huang YP, Lin TL, Lin TH, Wu HS. Antigenic and genetic diversity of human enterovirus 71 from 2009 to 2012, Taiwan. PLoS One. 2013;8(11):e80942. Epub 2013/12/19. pmid:24348916.

24. Chua BH, Phuektes P, Sanders SA, Nicholls PK, McMinn PC. The molecular basis of mouse adaptation by human enterovirus 71. J Gen Virol. 2008;89(Pt 7):1622–32. Epub 2008/06/19. pmid:18559932.

25. Huang SW, Wang YF, Yu CK, Su IJ, Wang JR. Mutations in VP2 and VP1 capsid proteins increase infectivity and mouse lethality of enterovirus 71 by virus binding and RNA accumulation enhancement. Virology. 2012;422(1):132–43. Epub 2011/11/15. pmid:22078110.

26. Zaini Z, McMinn P. A single mutation in capsid protein VP1 (Q145E) of a genogroup C4 strain of human enterovirus 71 generates a mouse-virulent phenotype. J Gen Virol. 2012;93(Pt 9):1935–40. Epub 2012/06/01. pmid:22647370.

27. Zhu J, Chen N, Zhou S, Zheng K, Sun L, Zhang Y, et al. Severity of enterovirus A71 infection in a human SCARB2 knock-in mouse model is dependent on infectious strain and route. Emerg Microbes Infect. 2018;7(1):205. Epub 2018/12/07. pmid:30518755.

28. Kataoka C, Suzuki T, Kotani O, Iwata-Yoshikawa N, Nagata N, Ami Y, et al. The role of VP1 amino acid residue 145 of enterovirus 71 in viral fitness and pathogenesis in a cynomolgus monkey model. PLoS Pathog. 2015;11(7):e1005033. Epub 2015/07/17. pmid:26181772.

29. Nishimura Y, Lee H, Hafenstein S, Kataoka C, Wakita T, Bergelson JM, et al. Enterovirus 71 binding to PSGL-1 on leukocytes: VP1-145 acts as a molecular switch to control receptor interaction. PLoS Pathog. 2013;9(7):e1003511. Epub 2013/08/13. pmid:23935488.

30. Zaini Z, Phuektes P, McMinn P. A reverse genetic study of the adaptation of human enterovirus 71 to growth in Chinese hamster ovary cell cultures. Virus Res. 2012;165(2):151–6. Epub 2012/03/13. pmid:22406130.

31. Reddi HV, Lipton HL. Heparan sulfate mediates infection of high-neurovirulence Theiler's viruses. J Virol. 2002;76(16):8400–7. Epub 2002/07/23. pmid:12134043.

32. Lee E, Lobigs M. Mechanism of virulence attenuation of glycosaminoglycan-binding variants of Japanese Encephalitis Virus and Murray Valley Encephalitis Virus. J Virol 2002;76(10):4901–11. pmid:11967307

33. Lee E, Hall RA, Lobigs M. Common E protein determinants for attenuation of glycosaminoglycan-binding variants of Japanese encephalitis and West Nile viruses. J Virol. 2004;78(15):8271–80. Epub 2004/07/16. pmid:15254199.

34. Lee E, Lobigs M. E protein domain III determinants of yellow fever virus 17D vaccine strain enhance binding to glycosaminoglycans, impede virus spread, and attenuate virulence. J Virol. 2008;82(12):6024–33. Epub 2008/04/11. pmid:18400851.

35. Mandl CW, Kroschewski H, Allison SL, Kofler R, Holzmann H, Meixner T, et al. Adaptation of tick-borne encephalitis virus to BHK-21 cells results in the formation of multiple heparan sulfate binding sites in the envelope protein and attenuation in vivo. J Virol. 2001;75(12):5627–37. Epub 2001/05/18. pmid:11356970.

36. Byrnes AP, Griffin DE. Large-plaque mutants of Sindbis virus show reduced binding to heparan sulfate, heightened viremia, and slower clearance from the circulation. J Virol. 2000;74(2):644–51. Epub 2000/01/07. pmid:10623725.

37. Gardner CL, Ebel GD, Ryman KD, Klimstra WB. Heparan sulfate binding by natural eastern equine encephalitis viruses promotes neurovirulence. Proc Natl Acad Sci U S A. 2011;108(38):16026–31. Epub 2011/09/08. pmid:21896745.

38. Gardner CL, Choi-Nurvitadhi J, Sun C, Bayer A, Hritz J, Ryman KD, et al. Natural variation in the heparan sulfate binding domain of the eastern equine encephalitis virus E2 glycoprotein alters interactions with cell surfaces and virulence in mice. J Virol. 2013;87(15):8582–90. Epub 2013/05/31. pmid:23720725.

39. Tan CW, Sam IC, Lee VS, Wong HV, Chan YF. VP1 residues around the five-fold axis of enterovirus A71 mediate heparan sulfate interaction. Virology. 2017;501:79–87. Epub 2016/11/23. pmid:27875780.

40. Xu Y, Ma S, Zhu L, Huang Z, Chen L, Xu Y, et al. Clinically isolated enterovirus A71 subgenogroup C4 strain with lethal pathogenicity in 14-day-old mice and the application as an EV-A71 mouse infection model. Antiviral Res. 2017;137:67–75. Epub 2016/11/20. pmid:27864074.

41. Fujii K, Sudaka Y, Takashino A, Kobayashi K, Kataoka C, Suzuki T, et al. VP1 amino acid residue 145 of enterovirus 71 is a key residue for its receptor attachment and resistance to neutralizing antibody during cynomolgus monkey infection. J Virol. 2018. Epub 2018/06/01. pmid:29848582.

42. Meng T, Kwang J. Attenuation of human enterovirus 71 high-replication-fidelity variants in AG129 mice. J Virol. 2014;88(10):5803–15. Epub 2014/03/14. pmid:24623423.

43. Sadeghipour S, Bek EJ, McMinn PC. Ribavirin-resistant mutants of human enterovirus 71 express a high replication fidelity phenotype during growth in cell culture. J Virol. 2013;87(3):1759–69. Epub 2012/11/24. pmid:23175376.

44. Vignuzzi M, Wendt E, Andino R. Engineering attenuated virus vaccines by controlling replication fidelity. Nat Med. 2008;14(2):154–61. Epub 2008/02/05. pmid:18246077.

45. Acevedo A, Woodman A, Arnold JJ, Yeh MT, Evans D, Cameron CE, et al. Genetic recombination of poliovirus facilitates subversion of host barriers to infection. J bioRxiv 2018:273060.

46. Zaini Z, Phuektes P, McMinn P. Mouse adaptation of a sub-genogroup B5 strain of human enterovirus 71 is associated with a novel lysine to glutamic acid substitution at position 244 in protein VP1. Virus Res. 2012;167(1):86–96. Epub 2012/05/12. pmid:22575826.

47. Caine EA, Moncla LH, Ronderos MD, Friedrich TC, Osorio JE. A single mutation in the VP1 of enterovirus 71 Is responsible for increased virulence and neurotropism in adult interferon-deficient mice. J Virol 2016;90(19):8592–604. pmid:27440896

48. Pfeiffer JK, Kirkegaard K. Increased fidelity reduces poliovirus fitness and virulence under selective pressure in mice. PLoS Pathog. 2005;1(2):e11. pmid:16220146

49. Vignuzzi M, Stone JK, Arnold JJ, Cameron CE, Andino R. Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature. 2006;439(7074):344–8. Epub 2005/12/04. pmid:16327776.

50. Zhu W, Li J, Liang G. How does cellular heparan sulfate function in viral pathogenicity? Biomed Environ Sci. 2011;24(1):81–7. Epub 2011/03/29. pmid:21440844.

51. Liu J, Dong W, Quan X, Ma C, Qin C, Zhang L. Transgenic expression of human P-selectin glycoprotein ligand-1 is not sufficient for enterovirus 71 infection in mice. Arch Virol. 2012;157(3):539–43. Epub 2011/12/22. pmid:22187102.

52. Fujii K, Nagata N, Sato Y, Ong KC, Wong KT, Yamayoshi S, et al. Transgenic mouse model for the study of enterovirus 71 neuropathogenesis. Proc Natl Acad Sci U S A. 2013;110(36):14753–8. Epub 2013/08/21. pmid:23959904.

53. Lin YW, Yu SL, Shao HY, Lin HY, Liu CC, Hsiao KN, et al. Human SCARB2 transgenic mice as an infectious animal model for enterovirus 71. PLoS One. 2013;8(2):e57591. Epub 2013/03/02. pmid:23451246.

54. Yamayoshi S, Koike S. Identification of a human SCARB2 region that is important for enterovirus 71 binding and infection. J Virol 2011;85(10):4937–46. pmid:21389126

55. Miyamura K, Nishimura Y, Abo M, Wakita T, Shimizu H. Adaptive mutations in the genomes of enterovirus 71 strains following infection of mouse cells expressing human P-selectin glycoprotein ligand-1. J Gen Virol. 2011;92(Pt 2):287–91. pmid:20943886.

56. Victorio CB, Xu Y, Ng Q, Meng T, Chow VT, Chua KB. Cooperative effect of the VP1 amino acids 98E, 145A and 169F in the productive infection of mouse cell lines by enterovirus 71 (BS strain). Emerg Microbes Infect. 2016;5:e60. Epub 2016/06/23. pmid:27329847.

57. Chang CK, Wu SR, Chen YC, Lee KJ, Chung NH, Lu YJ, et al. Mutations in VP1 and 5'-UTR affect enterovirus 71 virulence. Sci Rep. 2018;8(1):6688. Epub 2018/04/29. pmid:29703921.

58. Chen YC, Yu CK, Wang YF, Liu CC, Su IJ, Lei HY. A murine oral enterovirus 71 infection model with central nervous system involvement. J Gen Virol. 2004;85(Pt 1):69–77. Epub 2004/01/14. pmid:14718621.

59. Bernard KA, Klimstra WB, Johnston RE. Mutations in the E2 glycoprotein of Venezuelan equine encephalitis virus confer heparan sulfate interaction, low morbidity, and rapid clearance from blood of mice. Virology. 2000;276(1):93–103. Epub 2000/10/07. pmid:11021998.

60. Kobayashi K, Sudaka Y, Takashino A, Imura A, Fujii K, Koike S. Amino acid variation at VP1-145 of enterovirus 71 determines attachment receptor usage and neurovirulence in human scavenger receptor B2 transgenic mice. J Virol. 2018. Epub 2018/06/01. pmid:29848584.

61. Wang X, Li S-H, Zhu L, Nian Q-G, Yuan S, Gao Q, et al. Near-atomic structure of Japanese encephalitis virus reveals critical determinants of virulence and stability. Nature Communications. 2017;8(1):14. pmid:28446752

62. Cheng HY, Huang YC, Yen TY, Hsia SH, Hsieh YC, Li CC, et al. The correlation between the presence of viremia and clinical severity in patients with enterovirus 71 infection: a multi-center cohort study. BMC Infect Dis. 2014;14:417. Epub 2014/07/30. pmid:25069383.

63. Feng M, Guo S, Fan S, Zeng X, Zhang Y, Liao Y, et al. The preferential infection of astrocytes by enterovirus 71 plays a key role in the viral neurogenic pathogenesis. Frontiers in cellular and infection microbiology. 2016;6:192-. pmid:28066727.

64. Tan SH, Wong KT, Ong KC. Enterovirus 71 can directly infect the brainstem via cranial nerves and infection can be ameliorated by passive immunization. J Neuropathol Exp Neurol. 2014;73(11):999–1008. pmid:25289894

65. Ong KC, Badmanathan M, Devi S, Leong KL, Cardosa MJ, Wong KT. Pathologic characterization of a murine model of human enterovirus 71 encephalomyelitis. J Neuropathol Exp Neurol. 2008;67(6):532–42. Epub 2008/06/04. pmid:18520772.

66. Racaniello VR. One hundred years of poliovirus pathogenesis. Virology. 2006;344(1):9–16. https://doi.org/10.1016/j.virol.2005.09.015. pmid:16364730

67. Lancaster KZ, Pfeiffer JK. Limited trafficking of a neurotropic virus through inefficient retrograde axonal transport and the type I interferon response. PLoS Pathog. 2010;6(3):e1000791. pmid:20221252

68. Chen CS, Yao YC, Lin SC, Lee YP, Wang YF, Wang JR, et al. Retrograde axonal transport: a major transmission route of enterovirus 71 in mice. J Virol. 2007;81(17):8996–9003. Epub 2007/06/15. pmid:17567704.

69. Couderc T, Delpeyroux F, Le Blay H, Blondel B. Mouse adaptation determinants of poliovirus type 1 enhance viral uncoating. J Virol. 1996;70(1):305–12. Epub 1996/01/01. pmid:8523541.

70. Rincon V, Rodriguez-Huete A, Lopez-Arguello S, Ibarra-Molero B, Sanchez-Ruiz JM, Harmsen MM, et al. Identification of the structural basis of thermal lability of a virus provides a rationale for improved vaccines. Structure. 2014;22(11):1560–70. Epub 2014/10/14. pmid:25308865.

71. Yuan S, Li G, Wang Y, Gao Q, Wang Y, Cui R, et al. Identification of positively charged residues in enterovirus 71 capsid protein VP1 essential for production of infectious particles. J Virol. 2016;90(2):741–52. Epub 2015/10/30. pmid:26512078.

72. Bordo D, Argos P. Suggestions for "safe" residue substitutions in site-directed mutagenesis. J Mol Biol. 1991;217(4):721–9. Epub 1991/02/20. pmid:2005621.

73. Cohen FE. The Statistics of Sequence Similarity Scores. NCBI; 2006.

74. Stern A, Yeh MT, Zinger T, Smith M, Wright C, Ling G, et al. The evolutionary pathway to virulence of an RNA virus. Cell. 2017;169(1):35–46.e19. Epub 2017/03/25. pmid:28340348.

75. Zhang YX, Huang YM, Li QJ, Li XY, Zhou YD, Guo F, et al. A highly conserved amino acid in VP1 regulates maturation of enterovirus 71. PLoS Pathog. 2017;13(9):e1006625. Epub 2017/09/25. pmid:28938017.

76. Tseligka ED, Sobo K, Stoppini L, Cagno V, Abdul F, Piuz I, et al. A VP1 mutation acquired during an enterovirus 71 disseminated infection confers heparan sulfate binding ability and modulates ex vivo tropism. PLoS Pathog. 2018;14(8):e1007190. pmid:30075025

77. Cordey S, Petty TJ, Schibler M, Martinez Y, Gerlach D, van Belle S, et al. Identification of site-specific adaptations conferring increased neural cell tropism during human enterovirus 71 infection. PLoS Pathog. 2012;8(7):e1002826. Epub 2012/08/23. pmid:22910880.

78. Huang SW, Huang YH, Tsai HP, Kuo PH, Wang SM, Liu CC, et al. A selective bottleneck shapes the evolutionary mutant spectra of enterovirus A71 during viral dissemination in humans. J Virol. 2017;91(23). Epub 2017/09/22. pmid:28931688.

79. Lum LC, Wong KT, Lam SK, Chua KB, Goh AY, Lim WL, et al. Fatal enterovirus 71 encephalomyelitis. J Pediatr. 1998;133(6):795–8. Epub 1998/12/08. pmid:9842048.

80. Gao LD, Hu SX, Zhang H, Luo KW, Liu YZ, Xu QH, et al. Correlation analysis of EV71 detection and case severity in hand, foot, and mouth disease in the Hunan Province of China. PLoS One. 2014;9(6):e100003. Epub 2014/06/19. pmid:24941257.

81. Sabanathan S, Tan LV, Thwaites L, Wills B, Qui PT, Rogier van Doorn H. Enterovirus 71 related severe hand, foot and mouth disease outbreaks in South-East Asia: current situation and ongoing challenges. J Epidemiol Commun H 2014;68(6):500–2. pmid:24652348

82. Xing W, Liao Q, Viboud C, Zhang J, Sun J, Wu JT, et al. Hand, foot, and mouth disease in China, 2008–12: an epidemiological study. Lancet Infect Dis. 2014;14(4):308–18. Epub 2014/02/04. pmid:24485991.

83. Zhang J, Sun J, Chang Z, Zhang W, Wang Z, Feng Z. Characterization of hand, foot, and mouth disease in China between 2008 and 2009. Biomed Environ Sci. 2011;24(3):214–21. Epub 2011/07/26. pmid:21784305.

84. Yang CH, Liang CT, Jiang ST, Chen KH, Yang CC, Cheng ML, et al. A novel murine model expressing a chimeric mSCARB2/hSCARB2 receptor is highly susceptible to oral infection with clinical isolates of enterovirus 71. J Virol. 2019;93(11):e00183–19. pmid:30894476

85. Liou AT, Wu SY, Liao CC, Chang YS, Chang CS, Shih C. A new animal model containing human SCARB2 and lacking stat-1 is highly susceptible to EV71. Sci Rep. 2016;6:31151. Epub 2016/08/09. pmid:27499235.

86. Tan CW, Tee HK, Lee MH, Sam IC, Chan YF. Enterovirus A71 DNA-launched infectious clone as a robust reverse genetic tool. PLoS One. 2016;11(9):e0162771. Epub 2016/09/13. pmid:27617744.

87. Kammerer U, Kapp M, Gassel AM, Richter T, Tank C, Dietl J, et al. A new rapid immunohistochemical staining technique using the EnVision antibody complex. J Histochem Cytochem. 2001;49(5):623–30. Epub 2001/04/17. pmid:11304800.

88. Momany FA, Rone R. Validation of the general purpose QUANTA 3.2/CHARMm force field. Comput. Chem. 1992;13(7):888–900.