Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Analysis of Enterovirus 68 Strains from the 2014 North American Outbreak Reveals a New Clade, Indicating Viral Evolution

  • Juan Du,

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Baisong Zheng,

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Wenwen Zheng,

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Peng Li,

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Jian Kang,

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Jingwei Hou,

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Richard Markham,

    Affiliation Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

  • Ke Zhao ,

    paraake@icloud.com (KZ); xfyu@tju.edu.cn (XFY)

    Affiliation Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China

  • Xiao-Fang Yu

    paraake@icloud.com (KZ); xfyu@tju.edu.cn (XFY)

    Affiliations Institute of Virology and AIDS Research, First Hospital of Jilin University, Changchun, Jilin, China, Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of America

Abstract

Enterovirus 68 (EVD68) causes respiratory illness, mostly in children. Despite a reported low-level of transmission, the occurrence of several recent outbreaks worldwide including the 2014 outbreak in North America has raised concerns regarding the pathogenesis and evolution of EVD68. To elucidate the phylogenetic features of EVD68 and possible causes for the 2014 outbreak, 216 EVD68 strain sequences were retrieved from Genbank, including 22 from the 2014 outbreak. Several geographic and genotypic origins were established for these 22 strains, 19 of which were classified as Clade B. Of these 19 strains, 17 exhibited subsequent clustering and variation in protein residues involved in host-receptor interaction and/or viral antigenicity. Approximately 18 inter-clade variations were detected in VP1, which led to the identification of a new Clade D in EVD68 strains. The classification of this new clade was also verified by the re-construction of a Neighbor-Joining tree during the phylogenetic analysis. In addition, our results indicate that members of Clade B containing highly specific alterations in VP1 protein residues were the foremost contributors to the 2014 outbreak in the US. Altered host-receptor interaction and/or host immune recognition may explain the evolution of EVD68 as well as the global emergence and ongoing adaptation of this virus.

Introduction

Human Enterovirus 68 is a member of the genus Enterovirus that belongs to the family Picornaviridae. Human enteroviruses can be classified into four different groups (HEV-A to HEV-D) based on serological and phylogenetic analyses [1]. Enterovirus 68 has been classified as a member of group D and thus designated as EVD68. In contrast to other common enteroviruses such as Enterovirus 71 (EVA71) and Coxsackievirus A16 (CV-A16), both of which are categorized as HEV-A and are responsible for the widely spread hand, foot, and mouth disease (HFMD) [19], EVD68 causes respiratory illness with mild (sore throat, cough, and fever) to serious (including dyspnea or even respiratory failure) symptoms [1019].

Similar to other enteroviruses, EVD68 is a non-enveloped virus containing a ~7.4-kb RNA genome coding for one polyprotein. Subsequent to translation by the host machinery, this polyprotein is self-digested by viral proteins 2A and 3C yielding four structural proteins (VP1, VP2, VP3, and VP4) and seven non-structural proteins (2A, 2B, 2C, 3A, 3B, 3C, and an RNA-dependent RNA polymerase). Although all structural and non-structural proteins are important for the replication of EVD68, viral proteins VP1-4 participate directly in the assembly of new EVD68 virions as well as in host-receptor recognition during further rounds of infection [20].

EVD68 infections were first detected in 1962 [21] and were considered rare until the recent worldwide occurrences of clusters and outbreaks [15, 2231]. These outbreaks most likely indicate that EVD68 has continued to evolve over the past decades, resulting in increased viral infectivity, replication potency, and/or host immune evasion thus facilitating viral transmission. Previous analysis indicated that the most recent common ancestor (TMRCA) of EVD68 existed around 1960 [26], indicating that EVD68 is a relatively new enterovirus, compared to Enterovirus 71 (TMRCA arose in 1931) [4] or Coxsackievirus A16 (first isolated in 1951) [32], and that EVD68 continues to evolve and adapt to its host.

The recent EVD68 outbreaks indicate that EVD68 is well adjusted to transmission across hosts; the most recent outbreak occurred in North America in 2014 (the 2014 outbreak); over a thousand children in the US and another hundred in Canada were infected, with a minimum death toll of 14 [33]. Although similar outbreaks have been reported, as well as the involvement of EVD68 viral sequences [34, 35], there were no obvious cause(s) for the occurrence of the 2014 outbreak. Thus, phylogenetic analysis of EVD68 strains from the 2014 outbreak could provide up-to-date information regarding the evolution status of the virus.

In this study, we analyzed the VP1 region of 216 EVD68 strains by phylogenetic methods, including 22 from the 2014 outbreak. We determined that strains of Clade B EVD68 that were not detected in US in the past contributed most to the recent outbreak in 2014. Further analyses of VP1 variation indicated that altered viral antigenicity was among the potential causes of the 2014 outbreak. Meanwhile, we also detected a new clade, namely Clade D, based on several criteria for EVD68 clade classification. Our study therefore provides explanations for the occurrence of the 2014 outbreak and reveals the most recent evolutionary status of EVD68, which are essential for future control and prevention of EVD68 transmission.

Materials and Methods

Sequences

Initially, we intended to use the full-length VP1 region of EVD68 (nucleotides 2389–3315, corresponding to the EVD68 Fermon sequence [accession number, AY426531]) for this study, since it best correlates with the antigenicity of enteroviruses and has been widely used for enteroviral subtyping. However, in order to include additional EVD68 sequences without losing the number of sites required for the subsequent phylogenetic analyses, we decided to use a shorter region (2476–3129). Based on these criteria, 216 EVD68 sequences, including 22 from strains isolated during the 2014 outbreak, were retrieved from GenBank (S1 Table). However, full-length VP1 region of fewer (117 strains) EVD68 sequences were tested during some of the analyses, which is clearly indicated.

Bayesian analysis

Bayesian evolutionary analysis was performed using procedures similar to previously described methods [26]. Briefly, the retrieved EVD68 sequences were first aligned with MEGA5 (version 5.2.2) [36]. The evolutionary relationships of the viral sequences, along with TMRCA, were then inferred using the Bayesian Markov chain Monte Carlo (MCMC) method implemented with the BEAST package (version 1.8.1) [37]. The uncorrelated log-normal relaxed clock was used with an Extended Bayesian Skyline Plot and both HKY and SRD06 models for nucleotide substitution. The MCMC chains represent 200-million-iterations and the parameters were subsampled every 10,000 iterations.

The collected data were first evaluated with Tracer (version 1.6) to calculate the effective sample size (ESS); an ESS >300 was required for each of the parameters in order to confirm sampling efficiency. Next, TreeAnnotator (version 1.8.1) was used to generate a Bayesian evolutionary tree, removing 10% burn-in and summarizing the maximum clade credibility trees. Finally, the tree was visualized and modified with FigTree (version 1.4.2).

Phylogenetic analyses

To confirm the clade classification, a Neighbor-Joining tree was re-constructed with full-length VP1 sequences. In total, 117 EVD68 VP1 sequences were used and first aligned with MEGA5 [36]. The Neighbor-Joining tree was then re-constructed by using the Kimura-2-parameter model. Additional parameters were included during the re-construction to correlate with those for the Bayesian evolutionary analysis above: Gamma Distributed rates were applied with a parameter of 4, and only the 1st and 2nd positions of each codon was used. The re-constructed tree was then subjected to a phylogeny test based on the bootstrap method for 1,000 replicates.

To reveal possible evolutionary selection of each codon of EVD68 VP1, strains from each clade were aligned and the values of dN-dS for each codon were calculated with MEGA5 [36].

Protein remodeling

The structure of EVD68 VP1 was based on the published structure of the EVD68 virion (PDB: 4WM8) [20]. Remodeling and color manipulation were performed using PyMOL software (http://www.pymol.org, version 1.7.6.0).

In addition, the structure of Echovirus 7 (PDB, 2X5I) [38] was used as the template to reconstruct the structure of EVD68 with additional intact loop regions. The structure was proposed with MODELLER software (version 9.14) [39] using CA/AFP/v14T04344 (accession number, KM892502) as the subject. Remodeling, structural alignment, and color manipulation were performed using PyMOL.

Results

To gain a general understanding of the evolutionary status and viral relationships, a Bayesian phylogenetic tree was reconstructed using BEAST (version 1.8.1) [37] based on 216 published EVD68 VP1 sequences, including 22 sequences from the 2014 outbreak (Fig 1). Addition of the 2014 outbreak strains did not significantly alter the estimated year of TMRCA (1959.09, with a 95% highest probability density [HPD] of 1950.39–1961.34), compared to previously published studies [26]. Interestingly, the 2014 outbreak strains were located in several different clusters (Fig 1), with two strains being classified as the original Clade A, one as Clade C, and the remaining 19 as Clade B. The clustering of the 2014 strains with strains isolated outside of the US also suggested that the strains that caused the 2014 outbreak had multiple origins.

thumbnail
Fig 1. Bayesian MCMC analysis revealed Clade B members as the key contributors to the 2014 outbreak and demonstrated the emergence of a new Clade D from the original Clade.

A. Reported EVD68 sequences retrieved from GenBank are labeled in different colors in this tree based on sample origin, with the 2014 outbreak strains in red. The Bayesian phylogenetic tree was reconstructed with the partial VP1 coding region (2,435–3,128) according to a previously described method [26]. The original Clade A (ori A) classification is shown in addition to the four clades classified in this study. The scale bar below the tree indicates the time of strain sampling. Triangle markers indicate nodes for new clade classification. Specific strain labeling: Hashtag, EVD68 Fermon strain; 1, CA/AFP/11-1767; 2, US/KY/14-18953; 3, NYC strains; 4, CA/RESP/10-786; 5, MD02-1 and MD02-2. The year of the most recent common ancestor (TMRCA) and the Bayesian posterior probability (BPP) for each clade are shown in the table.

https://doi.org/10.1371/journal.pone.0144208.g001

Of the 19 Clade B strains that predominately contributed (86.36%) to the last known EVD68 outbreak in North America in 2014, 17 belonged to a subclade that separated from other strains around the year 2007 (Fig 1). The clustering of this subclade was confirmed with a Bayesian Posterior Probability (BPP) value of 1. Moreover, among these 17 strains, CA/AFP/11-1767 (strain number 1 in Fig 1) was a relatively early strain, indicating that the remaining 16 might be evolutionary progeny. A more detailed examination of the VP1 protein sequences indicated that CA/AFP/11-1767 contains two residues at positions 98 and 290 that differ from those found in most B1 strains (Fig 2A, position 99 is indicated instead of 98 since the latter was missing from the published structure). Interestingly, based on published structural information for EVD68 [20], position 290 occupies a surface position that is mostly likely linked with the antigenicity of VP1 (Fig 2B). Position 98 is located on a loop region surrounding the icosahedral five-fold axes (or the “canyon” region) responsible for host-receptor interaction (Fig 2B). The changes to both residues altered the length of the side chain; the change at position 98 also changed the polarity of the side chain (Fig 2A).

thumbnail
Fig 2. Specific residue variations detected in the EVD68 VP1 protein from the 2014 outbreak strains.

(A) VP1 protein alignment based on the Clade B subcluster. The alignment was performed with MEGA5 [36]. Hyphens indicate residues identical to those in CA/AFP/11-1767; “X” indicates unidentified amino acid due to unidentified nucleotides in the submitted CA/AFP/v12T04950 sequence. Numbers indicate protein positions. Sequence CQ5585 is missing the fragment surrounding position 290. The subtree for the subclustered Clade B strains is extracted from the tree in Fig 1. (B) Structures of viral proteins VP1 (green), VP2 (purple), and VP3 (cyan) in the surface model. The remodeling was based on the published EVD68 structure (PDB: 4WM8) [20]. Positions 218 and 290, located on the surface of VP1, are labeled in red. Position 99 is shown in red instead of 98 since the latter is missing in the published EVD68 structure. (C) Position 194, is located beneath the viral surface and is essential for β-sheet formation, which may be important for supporting EVD68 virion “canyon” formation. Left panel: the surface model of VP1 rotated 120° from the model shown in Fig 2B; right panel: a cartoon model of the β-sheet in which residue 194 interacts with residue 183.

https://doi.org/10.1371/journal.pone.0144208.g002

Similar residue substitutions occurred at positions 194 and 218 in most of the remaining B1 strains (12 of 16; 75.00%; Fig 2A). In contrast to positions 98 and 290 that are located on the central surface region, position 218 is located on the side of VP1 and is possibly involved in VP1-VP1 interactions and hence in “canyon” formation (Fig 2B). Position 194 is not located on the surface of EVD68 but rather faces inward (Fig 2C). Despite the moderate decrease in the length of the side chain, the alteration from methionine to isoleucine at position 194 is highly conserved in these 12 B1 strains. Furthermore, the cartoon model indicated that although position 194 was very close to VP2 (Fig 2C), it was not involved in VP1-VP2 interactions. Rather, it participated in the formation of a β-sheet connected to loops that contributed to “canyon” formation (Fig 2C). Importantly, the alterations at positions 194 and 218 seem to co-exist in addition to the changes at positions 98 and 290 (Fig 2A), suggesting a possible synergistic effect on the binding of VP1 to the cellular receptor. Therefore, the 2014 outbreak was very likely due to modification of the receptor interaction affecting local cell tropism; the inability of previously established anti-EVD68 antibodies to detect the altered residue at position 290 might have contributed to the occurrence of the outbreak.

Although other EVD68 clades contributed slightly to the 2014 outbreak, one of the three non-Clade B strains, US/KY/14-18953 (strain number 2 in Fig 1), piqued our interest. This strain was first determined as belonging to Clade A based on a previous classification [26], however it was evolutionarily distant from other Clade A strains isolated in the US, such as the NYC strains from 2009 (orange color in Clade A, Fig 1) and CA/RESP/10-786 from 2014 (strain number 4 in Fig 1). Moreover, unlike previous reconstructions of EVD68 phylogenetic trees [26, 30], we detected an additional cluster with a BPP value of 1 that included strain US/KY/14-18953 (Fig 1). Most of the remaining Clade A sequences formed another cluster (Fig 1) originating in 2001; this new clade also had a BPP of 1.

Full-length VP1 is used most frequently to determined the subtype for other enteroviruses such as EVA71 and CA16. To avoid any bias introduced to the Bayesian analysis above based on only a partial VP1 sequence, a different approach based on full-length VP1 was necessary to confirm the separation from Clade A. Out of 226 strains used in this study, 117 of them contained full-length VP1 and these full-length VP1 sequenes were then subjected to a Neighbor-Joining (NJ) tree re-construction with MEGA5 [36] (Fig 3). Interestingly, a more significant separation from Clade A was observed, and the cluster containing USA/KY/14-18953 was phylogenetically distant from the remaining Clade A strains (Fig 3A and 3B). On the other hand, the NJ tree presented clustering of Clade A, B, and C, similar to that detected in the Bayesian analysis (Fig 3B, 3C and 3D), confirming the validity of the NJ tree re-construction. In addition, an amino acid comparison of the full-length VP1 regions showed that the difference between the two Clade A clusters were similar to those between the other clades (Table 1), indicating that US/KY/14-18953 belongs to a new clade (Clade D, Fig 1). Clade D is distinct from Clade A and emerged during recent instances of EVD68 transmission.

thumbnail
Fig 3. Phylogenetic analysis of full-length VP1 region of EVD68 strains.

117 full-length EVD68 VP1 sequences were used for the re-construction of the Neighbor-Joining tree (A). Similar clade formation was detected compared to that in Fig 1, and each clade was shown as a subtree (B, C, and D, for Clade A, B, and C respectively). The tree was re-constructed with the help of MEGA5 [36]. Only bootstrap values >70 were shown. For a clearer presentation of the relationship between strains, only the topology was shown. Circle symbols indicate EVD68 strain from the 2014 outbreak, while square symbols indicate compressed EVD68 clades.

https://doi.org/10.1371/journal.pone.0144208.g003

thumbnail
Table 1. Amino acid variation in the VP1 residues of different EV68 clades.

https://doi.org/10.1371/journal.pone.0144208.t001

To further confirm the accuracy of our new clade classification, the primary protein sequences of the VP1 region were compared in detail; these sequences have been used in the phylogenetic classification of EVD68, which correlated with the serological classification. In total, 18 positions with amino acids that differed among the clades were detected (Table 2). A purifying selection was observed for each of these 18 positions within at least each of two clades, where a dN-dS value ≤0 was observed (Fig 4). Therefore, the pattern for both intra-clade conservation and inter-clade variation in Table 2 validated the novel clade classification.

thumbnail
Fig 4. Evolutionary selection on the 18 clade-specific positions.

After sequence alignment, the dN-dS values for each clade were calculated for these clade-specific positions with the help of MEGA5 [36]. Values >0 indicate positive selection, while values <0 indicate purifying selection. Each codon with dN-dS = 0 in this study shared conserved nucleotides and thus was considered as purifyingly selected. Clade D was excluded from this analysis because only two strains contained full-length VP1.

https://doi.org/10.1371/journal.pone.0144208.g004

thumbnail
Table 2. Specific amino acids for EV68 clade classification.

https://doi.org/10.1371/journal.pone.0144208.t002

An interesting question arising from the novel EVD68 clade classification was why these specific 18 positions varied. To better understand the impact of these positions on the EVD68 virion, we attempted to apply these variations to previously published EVD68 structures [20]. However, many of these positions (positions 2, 92, 141, 143, 144, 145, and 148) could not be located on the model (Fig 5A). Therefore, we proposed another structural model of EVD68 VP1 based on an additional published enterovirus structure (PDB: 2X5I) [38]. Structural alignment using the PyMol software suggested that the proposed model was similar to the original EVD68 structure, despite minor differences at relaxed regions (Fig 5B). Using both structures, we were able to locate 17 positions listed in Table 2. Interestingly, 11 of these positions are located at the VP1 “canyon” formation region (Fig 5A and 5C); all 11 residues are situated on four loops, including the BC-loop and DE-loop (Fig 5C) surrounding the icosahedral five-fold axes (i.e., the “canyons”), which are considered as the most important sites for both receptor interaction [20] and neutralizing immunogenicity [40]. In addition, positions 297 and 305 are located on the surface of EVD68, and therefore, might contribute to virion antigenicity. Variations in these positions would benefit EVD68 during adaptation of local host tropism, further confirming the new classification proposed in this study.

thumbnail
Fig 5. Host-receptor binding and/or viral antigenicity correlate with EVD68 clade classification.

(A). Surface models of EVD68 VP1 (green), VP2 (purple), and VP3 (cyan), with clade-specific residues labeled in red. The remodeling was based on the published EVD68 structure (PDB: 4WM8) [20]. Positions 140 and 149 are labeled in red to indicate the approximate positions of residues 141, 143, 144, 145, and 148 that are missing in the published structure. (B). Structural alignment of VP1 from the published EVD68 structure (green) or the proposed structure (blue) based on the viral structure of Echovirus 7 (PDB: 2X5I) [38]. The alignment was performed with PyMol. (C). Surface model of the proposed EVD68 VP1 structure (blue), with clade-specific residues indicated in red. The panel on the right provides a detailed cartoon model showing clade-specific residues (in red) located on the loops forming the “canyon” responsible for host-receptor interaction. The BC-loop and DE-loop, thought to be critical in both host-receptor binding and viral antigenicity [40], are labeled with orange frames.

https://doi.org/10.1371/journal.pone.0144208.g005

Discussion

Despite their size and dependence on host cells, viruses continue to evolve, through mechanisms such as ligand adaptation, immune evasion, host machinery hijacking, and/or viral enzyme improvement, in order to enhance their replication within the host. In EVD68, the former two mechanisms require variation of viral capsid proteins. The results of our present study analyzing the VP1 region of EVD68 strains from the 2014 outbreak revealed a recently derived novel EVD68 clade, Clade D, which contains unique residues at certain positions, including ones responsible for receptor and antibody recognition [40]. In addition, we propose that the TMRCA for Clade C be moved to the year 2003.6, since no strains similar to MD02-1 and MD02-2 (number 5 in Fig 1) were isolated in the US. This new clade classification provides a better interpretation of the Bayesian phylogenetic tree, since many reported EVD68 strains isolated prior to 2003 (such MD02-1 and MD02-2) stand alone instead of causing clusters or outbreaks. Sequence analysis of these unclassified strains revealed the occurrence of reverse mutations (data not shown), indicating that these strains were adapting to the host and were probably evolutionary intermediates of EVD68. Interestingly, this classification determined similar generation time points for each of the clades (2005.0, 2004.6, 2003.6, and 2005.5 for Clade A, B, C, and D, respectively), possibly indicating the existence of a yet-unknown event that occurred around 2003–2005 during the emergence of global EVD68 transmission.

Although many of the early EVD68 strains such as the Fermon reference strain were isolated in the US, no clear cluster or outbreak was detected in the US, except for one that occurred in New York City in 2009 involving Clade A strains [26]. In contrast, although all the 2014 outbreak strains clustered with strains from other countries, there was no clear evidence suggesting that they were epidemiologically linked. In particular, transmission of the Clade B strains, which were key contributors to the 2014 outbreak, most likely originated with US strain CA/AFP/11-1767 and then spread through Asia and/or Europe and back to the US where they caused the outbreak. It is also possible that CA/AFP/11-1767 continued to evolve in the US, until it was maximally adapted to the local host. Therefore, it is difficult to determine whether the strains from the 2014 outbreak were of domestic or foreign origin.

Similar to other enteroviruses, EVD68 VP1 is mostly involved in viral binding to host receptors as well as in host immune recognition [20]. Viral evolution involving VP1 is obviously persistent, since variations at four positions contributed to the severe outbreak of EVD68 in the US. These changes are most likely essential, since no previous Clade B infections were detected in the US, although many Clade B cases had been reported in Europe and Asia prior to 2014. Indeed, based on the published EVD68 structure [20], these four positions are associated with host-receptor recognition and/or antigenicity. Alterations at these positions could promote cell entry as well as enhance immune evasion. In addition to these four intra-clade variations, we identified 18 inter-clade variations in VP1, most of which were also located on the viral surface of EVD68 and were associated with “canyon” formation, supporting both the phylogenetic and serological classification of EVD68.

Our findings have important implications in terms of treatment and vaccine development for EVD68 or enteroviruses in general. In theory, preventing viral entry by blocking the “canyon”-receptor interaction would constitute the best way to cure or, at the least, control enteroviral infection. Previous studies have confirmed the role of VP1 variation as a major factor affecting enteroviral immunogenicity [9, 41]. In the case of EVD68, the anti-rhinovirus drug pleconaril has already been shown to interact with VP1 and suppress viral pathogenesis in HeLa cells [20]. The published viral inhibitors of EVA71, one of the major causative pathogens of HFMD, were also designed to interact with the “canyon” region [42]. In contrast to these previous findings, the data obtained in this study suggests that the VP1 protein is capable of changing under selective pressure; a concept applicable to EVD68 and enteroviruses in general. Moreover, the identification of the Clade B subclade indicates that such changes can occur very easily, since the four protein variations occurred within just a few years. In addition, we have previously observed that a single amino acid variation in EVA71 VP1 leads to significant differences in the host immune response [9]. Such occurrences in EVD68 (and probably other enteroviruses as well) poses challenges in the development of antivirals and/or vaccines, the success of which may require the development of novel strategies, such as the design of antiviral drugs that suppress the function of viral protease 3C [43].

In summary, we determined that Clade B EVD68 strains had the greatest impact on the 2014 outbreak in North America. Both sequence and structural analyses indicate that viral adaptation to the local host environment and/or viral evasion from host immune recognition most likely contributed to the spread of Clade B EVD68 strains in the US. In addition, a new Clade D was classified based on Bayesian evolutionary analysis, phylogenetic analysis, as well as differences in the VP1 protein. Whether this new Clade may contribute to future EVD68 transmission is an interesting question. Furthermore, detailed structural information indicated that variation in “canyon” formation in VP1 is the probable cause of the generation of different EVD68 clades and altered viral antigenicity. Multiple analyses of the VP1 protein suggest that the recently generated EVD68 strains continue to evolve and adapt to the human host. Therefore, further efforts are required to monitor EVD68 transmission and evolution in order to better control and prevent further damage caused by this emerging pathogen.

Supporting Information

S1 Table. EVD68 strains used in this study.

https://doi.org/10.1371/journal.pone.0144208.s001

(DOCX)

Acknowledgments

We wish to thank Drs. Jingliang Li and Junliang Chang for interesting discussions. In addition we thank Editage for the provided editing services.

Author Contributions

Conceived and designed the experiments: KZ XFY. Performed the experiments: JD BZ WZ PL JK JH KZ. Analyzed the data: JD BZ RM KZ XFY. Contributed reagents/materials/analysis tools: JD BZ KZ XFY. Wrote the paper: JD RM KZ XFY.

References

  1. 1. Zhao K, Han X, Wang G, Hu W, Zhang W, Yu XF. Circulating coxsackievirus A16 identified as recombinant type A human enterovirus, China. Emerg Infect Dis. 2011;17(8):1537–1540. pmid:21801645; PubMed Central PMCID: PMC3381541.
  2. 2. McWilliam Leitch EC, Cabrerizo M, Cardosa J, Harvala H, Ivanova OE, Koike S, et al. The association of recombination events in the founding and emergence of subgenogroup evolutionary lineages of human enterovirus 71. J Virol. 2012;86(5):2676–2685. pmid:22205739; PubMed Central PMCID: PMC3302253.
  3. 3. Wang X, Zhu C, Bao W, Zhao K, Niu J, Yu XF, et al. Characterization of full-length enterovirus 71 strains from severe and mild disease patients in northeastern China. PLoS One. 2012;7(3):e32405. pmid:22479324; PubMed Central PMCID: PMC3315555.
  4. 4. Bessaud M, Razafindratsimandresy R, Nougairede A, Joffret ML, Deshpande JM, Dubot-Peres A, et al. Molecular comparison and evolutionary analyses of VP1 nucleotide sequences of new African human enterovirus 71 isolates reveal a wide genetic diversity. PLoS One. 2014;9(3):e90624. pmid:24598878; PubMed Central PMCID: PMC3944068.
  5. 5. Li J, Chang J, Liu X, Yang J, Guo H, Wei W, et al. Protection from lethal challenge in a neonatal mouse model by circulating recombinant form coxsackievirus A16 vaccine candidates. J Gen Virol. 2014;95(Pt 5):1083–1093. pmid:24496826; PubMed Central PMCID: PMC3983757.
  6. 6. Li Z, Yu J, Liu L, Wei Z, Ehrlich ES, Liu G, et al. Coxsackievirus A16 infection induces neural cell and non-neural cell apoptosis in vitro. PLoS One. 2014;9(10):e111174. pmid:25350381; PubMed Central PMCID: PMC4211689.
  7. 7. Wei W, Guo H, Li J, Ren S, Wei Z, Bao W, et al. Circulating HFMD-associated coxsackievirus A16 is genetically and phenotypically distinct from the prototype CV-A16. PLoS One. 2014;9(4):e94746. pmid:24736564; PubMed Central PMCID: PMC3988102.
  8. 8. Chang J, Li J, Liu X, Liu G, Yang J, Wei W, et al. Broad protection with an inactivated vaccine against primary-isolated lethal enterovirus 71 infection in newborn mice. BMC Microbiol. 2015;15:139. pmid:26169371; PubMed Central PMCID: PMC4501189.
  9. 9. Chang J, Li J, Wei W, Liu X, Liu G, Yang J, et al. Determinants of EV71 immunogenicity and protection against lethal challenge in a mouse model. Immunol Res. 2015;62(3):306–15. pmid:26025091.
  10. 10. Oberste MS, Maher K, Schnurr D, Flemister MR, Lovchik JC, Peters H, et al. Enterovirus 68 is associated with respiratory illness and shares biological features with both the enteroviruses and the rhinoviruses. J Gen Virol. 2004;85(Pt 9):2577–2584. pmid:15302951.
  11. 11. Hasegawa S, Hirano R, Okamoto-Nakagawa R, Ichiyama T, Shirabe K. Enterovirus 68 infection in children with asthma attacks: virus-induced asthma in Japanese children. Allergy. 2011;66(12):1618–20. pmid:21958204.
  12. 12. Kaida A, Kubo H, Sekiguchi J, Kohdera U, Togawa M, Shiomi M, et al. Enterovirus 68 in children with acute respiratory tract infections, Osaka, Japan. Emerg Infect Dis. 2011;17(8):1494–1497. pmid:21801632; PubMed Central PMCID: PMC3381549.
  13. 13. Kreuter JD, Barnes A, McCarthy JE, Schwartzman JD, Oberste MS, Rhodes CH, et al. A fatal central nervous system enterovirus 68 infection. Arch Pathol Lab Med. 2011;135(6):793–796. pmid:21631275.
  14. 14. Rahamat-Langendoen J, Riezebos-Brilman A, Borger R, van der Heide R, Brandenburg A, Scholvinck E, et al. Upsurge of human enterovirus 68 infections in patients with severe respiratory tract infections. J Clin Virol. 2011;52(2):103–106. pmid:21802981.
  15. 15. Ikeda T, Mizuta K, Abiko C, Aoki Y, Itagaki T, Katsushima F, et al. Acute respiratory infections due to enterovirus 68 in Yamagata, Japan between 2005 and 2010. Microbiol Immunol. 2012;56(2):139–143. pmid:22309616.
  16. 16. Jacobson LM, Redd JT, Schneider E, Lu X, Chern SW, Oberste MS, et al. Outbreak of lower respiratory tract illness associated with human enterovirus 68 among American Indian children. Pediatr Infect Dis J. 2012;31(3):309–312. pmid:22315004.
  17. 17. Xiang Z, Gonzalez R, Wang Z, Ren L, Xiao Y, Li J, et al. Coxsackievirus A21, enterovirus 68, and acute respiratory tract infection, China. Emerg Infect Dis. 2012;18(5):821–824. pmid:22516379; PubMed Central PMCID: PMC3358056.
  18. 18. Renois F, Bouin A, Andreoletti L. Enterovirus 68 in pediatric patients hospitalized for acute airway diseases. J Clin Microbiol. 2013;51(2):640–643. pmid:23224095; PubMed Central PMCID: PMC3553905.
  19. 19. Imamura T, Suzuki A, Lupisan S, Kamigaki T, Okamoto M, Roy CN, et al. Detection of enterovirus 68 in serum from pediatric patients with pneumonia and their clinical outcomes. Influenza Other Respir Viruses. 2014;8(1):21–24. pmid:24209770; PubMed Central PMCID: PMC4177794.
  20. 20. Liu Y, Sheng J, Fokine A, Meng G, Shin WH, Long F, et al. Structure and inhibition of EV-D68, a virus that causes respiratory illness in children. Science. 2015;347(6217):71–74. pmid:25554786; PubMed Central PMCID: PMC4307789.
  21. 21. Schieble JH, Fox VL, Lennette EH. A probable new human picornavirus associated with respiratory diseases. Am J Epidemiol. 1967;85(2):297–310. pmid:4960233.
  22. 22. Lauinger IL, Bible JM, Halligan EP, Aarons EJ, MacMahon E, Tong CY. Lineages, sub-lineages and variants of enterovirus 68 in recent outbreaks. PLoS One. 2012;7(4):e36005. pmid:22536453; PubMed Central PMCID: PMC3335014.
  23. 23. Linsuwanon P, Puenpa J, Suwannakarn K, Auksornkitti V, Vichiwattana P, Korkong S, et al. Molecular epidemiology and evolution of human enterovirus serotype 68 in Thailand, 2006–2011. PLoS One. 2012;7(5):e35190. pmid:22586446; PubMed Central PMCID: PMC3346751.
  24. 24. Meijer A, van der Sanden S, Snijders BE, Jaramillo-Gutierrez G, Bont L, van der Ent CK, et al. Emergence and epidemic occurrence of enterovirus 68 respiratory infections in The Netherlands in 2010. Virology. 2012;423(1):49–57. pmid:22177700.
  25. 25. Piralla A, Lilleri D, Sarasini A, Marchi A, Zecca M, Stronati M, et al. Human rhinovirus and human respiratory enterovirus (EV68 and EV104) infections in hospitalized patients in Italy, 2008–2009. Diagn Microbiol Infect Dis. 2012;73(2):162–167. pmid:22494556.
  26. 26. Tokarz R, Firth C, Madhi SA, Howie SR, Wu W, Sall AA, et al. Worldwide emergence of multiple clades of enterovirus 68. J Gen Virol. 2012;93(Pt 9):1952–1958. pmid:22694903; PubMed Central PMCID: PMC3542132.
  27. 27. Jaramillo-Gutierrez G, Benschop KS, Claas EC, de Jong AS, van Loon AM, Pas SD, et al. September through October 2010 multi-centre study in the Netherlands examining laboratory ability to detect enterovirus 68, an emerging respiratory pathogen. J Virol Methods. 2013;190(1–2):53–62. pmid:23458694.
  28. 28. Lu QB, Wo Y, Wang HY, Wei MT, Zhang L, Yang H, et al. Detection of enterovirus 68 as one of the commonest types of enterovirus found in patients with acute respiratory tract infection in China. J Med Microbiol. 2014;63(Pt 3):408–414. pmid:24324030.
  29. 29. Meijer A, Benschop KS, Donker GA, van der Avoort HG. Continued seasonal circulation of enterovirus D68 in the Netherlands, 2011–2014. Euro Surveill. 2014;19(42). pmid:25358039.
  30. 30. Opanda SM, Wamunyokoli F, Khamadi S, Coldren R, Bulimo WD. Genetic diversity of human enterovirus 68 strains isolated in Kenya using the hypervariable 3'-end of VP1 gene. PLoS One. 2014;9(7):e102866. pmid:25054861; PubMed Central PMCID: PMC4108357.
  31. 31. Piralla A, Girello A, Grignani M, Gozalo-Marguello M, Marchi A, Marseglia G, et al. Phylogenetic characterization of enterovirus 68 strains in patients with respiratory syndromes in Italy. J Med Virol. 2014;86(9):1590–1593. pmid:24155220.
  32. 32. Poyry T, Hyypia T, Horsnell C, Kinnunen L, Hovi T, Stanway G. Molecular analysis of coxsackievirus A16 reveals a new genetic group of enteroviruses. Virology. 1994;202(2):982–987. pmid:8030260.
  33. 33. Khan F. Enterovirus D68: acute respiratory illness and the 2014 outbreak. Emerg Med Clin North Am. 2015;33(2):e19–32. pmid:26065305.
  34. 34. Brown BA, Nix WA, Sheth M, Frace M, Oberste MS. Seven Strains of Enterovirus D68 Detected in the United States during the 2014 Severe Respiratory Disease Outbreak. Genome Announc. 2014;2(6). pmid:25414503; PubMed Central PMCID: PMC4239358.
  35. 35. Midgley CM, Jackson MA, Selvarangan R, Turabelidze G, Obringer E, Johnson D, et al. Severe respiratory illness associated with enterovirus D68—Missouri and Illinois, 2014. MMWR Morb Mortal Wkly Rep. 2014;63(36):798–799. pmid:25211545.
  36. 36. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–2739. pmid:21546353; PubMed Central PMCID: PMC3203626.
  37. 37. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–1973. pmid:22367748; PubMed Central PMCID: PMC3408070.
  38. 38. Plevka P, Hafenstein S, Harris KG, Cifuente JO, Zhang Y, Bowman VD, et al. Interaction of decay-accelerating factor with echovirus 7. J Virol. 2010;84(24):12665–12674. pmid:20881044; PubMed Central PMCID: PMC3004353.
  39. 39. Fiser A, Sali A. Modeller: generation and refinement of homology-based protein structure models. Methods Enzymol. 2003;374:461–491. pmid:14696385.
  40. 40. Imamura T, Okamoto M, Nakakita S, Suzuki A, Saito M, Tamaki R, et al. Antigenic and receptor binding properties of enterovirus 68. J Virol. 2014;88(5):2374–2384. pmid:24371050; PubMed Central PMCID: PMC3958110.
  41. 41. Li J, Liu G, Liu X, Yang J, Chang J, Zhang W, et al. Optimization and Characterization of Candidate Strain for Coxsackievirus A16 Inactivated Vaccine. Viruses. 2015;7(7):3891–3909. pmid:26193302; PubMed Central PMCID: PMC4517132.
  42. 42. De Colibus L, Wang X, Spyrou JA, Kelly J, Ren J, Grimes J, et al. More-powerful virus inhibitors from structure-based analysis of HEV71 capsid-binding molecules. Nat Struct Mol Biol. 2014;21(3):282–288. pmid:24509833; PubMed Central PMCID: PMC4530014.
  43. 43. Tan J, George S, Kusov Y, Perbandt M, Anemuller S, Mesters JR, et al. 3C protease of enterovirus 68: structure-based design of Michael acceptor inhibitors and their broad-spectrum antiviral effects against picornaviruses. J Virol. 2013;87(8):4339–4351. pmid:23388726; PubMed Central PMCID: PMC3624371.