Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998

Genetic characteristics and phylogenetic analysis of coxsackievirus A6 isolated in Linyi, China, 2022-2023

Hand, foot and mouth disease (HFMD) was one of the most common infectious disease in the past few decades. After 2013, coxsackievirus A6 (CVA6) has replaced enterovirus 71 (EV-A71) and coxsackievirus A16 (CVA16), becoming the predominant pathogen responsible for HFMD in many areas in China. The objective of this study is to investigate the genetic characteristics and molecular epidemiology of CVA6 in Linyi from 2022 to 2023. A total of 965 HFMD cases were enrolled in this study and analyses based on VP1 nucleotide sequences were performed to determine the evolutionary trajectory of CVA6. In 2022, 281 (281/386, 72.8%) were positive for enterovirus (EVs), and 217 (217/281, 77.2%) were CVA6 positive. In 2023, 398 (398/579, 68.7%) samples were positive for EVs, and 243 (243/398, 61.1%) were CVA6 positive. Six sequences were selected from each year for the homology analysis. The results showed that 12 strains isolated in Linyi were far from the prototype strain (AY421764) and the first CVA6 strain reported in China (JQ364886). Phylogenetic analysis showed that the CVA6 strains isolated in Linyi all belonged to D3 subgenotype. CVA6 is emerging as a common pathogen causing HFMD in Linyi, and continuous surveillance of HFMD etiological agents is necessary.

Inhibition of EV71 replication by an interferon-stimulated gene product L3HYPDH

Enterovirus 71 (EV71) is the common causative agent of hand-foot-mouth disease (HFMD). Despite evidence in mice model suggested that the interferon (IFN) signaling pathways play a role in defending against this virus, knowledge on the IFN-mediated antiviral response is still limited. Here we identified an IFN-stimulated gene (ISG) called L3HYPDH, whose expression inhibits EV71 replication. Mapping assay indicated that amino acids 61-120 and 295-354 are critical for its optimal antiviral activity. Mechanismly, L3HYPDH specifically inhibits protein translation mediated by EV71 internal ribosome entry site (IRES). Our data thus uncovered a new mechanism utilized by the host cell to restrict EV71 replication.

Phylogenetic characteristics and recombination analysis of echovirus 5 associated with severe acute respiratory infection in China

IMPORTANCE

This study is the first report of echovirus 5 (E5) associated with severe acute respiratory infection and obtained the first E5 whole-genome sequence in China. Combined with the sequences available in the GenBank database, the first genotyping, phylogenetic characteristics, recombination, and genetic evolutionary analysis of E5 was performed in this study. Our findings providing valuable information on global E5 molecular epidemiology.

Characterizing enterovirus C96 genome and phylodynamics analysis

Enterovirus C96 (EV-C96) is a recently discovered serotype belonging to enterovirus C species. It had been isolated from patients with acute flaccid paralysis, hand, foot, and mouth disease, diarrhea, healthy people, or environmental specimens. Despite increasing reports of the virus, the small number of full-length genomes available for EV-C96 has limited molecular epidemiological studies. In this study, newly collected rare EV-C96 strains in China from 1997 to 2020 were combined with sequences available in GenBank for comprehensive analyses. Sequence analysis revealed that the nucleotide sequence similarity of EV-C96 and the prototype strain (BAN00-10488) was 75%-81.8% and the amino acid sequence similarity was 85%-94.9%. EV-C96 had a high degree of genetic variation and could be divided into 15 genogroups. The mean evolutionary rate was 5.16 × 10-3 substitution/site/year, and the most recent common ancestor was dated to 1925. A recombination analysis revealed that EV-C96 may be a recombinant derived from other serotypes in the EV-C group in the nonstructural protein coding region. This comprehensive and integrated analysis of the whole genome sequence of EV-C96 provides valuable data for further studies on the molecular epidemiology of EV-C96 worldwide.

Surveillance, epidemiology, and impact of the coronavirus disease 2019 interventions on the incidence of enterovirus infections in Nanchang, China, 2010-2022

Introduction

Pathogen spectrum of Hand, foot and mouth disease (HFMD) has substantially changed in the past decade in China. Growing evidence has indicated that anti-COVID-19 nonpharmaceutical interventions (NPIs) can support control of various infectious diseases, including intestinal diseases.

Methods

In this study, HFMD cases were enrolled from sentinel hospitals of Nanchang, Jiangxi province, and enteroviruses were genotyped using specific real time RT-PCR. We systematically characterized the epidemiology of HFMD based on the continuous molecular surveillance and estimated the impact of COVID-19 intervention on HFMD incidence using seasonal autoregressive integrated moving average (ARIMA) models.

Results

A total of 10247 HFMD cases were included during 2010-2022, of which 6121 enterovirus (EV)-positive cases (59.7%) were identified by real-time RT-PCR. Over 80% cases were associated with EV-A71 and coxsackievirus A16 (CVA16) during 2010-2012, while the type distribution significantly changed as CVA6 emerged to be dominant, accounting for 22.6%-59.6% during 2013-2022. It was observed that the prevalence patterns of EV-A71 and CVA16 were similar and both of them peaked in the second quarter and then leveled off. However, CVA6 was generally prevalent around the fourth quarter, demonstrating a staggered prevalence during 2010-2019. During the COVID-19 epidemic, the seasonal HFMD epidemic peak was restrained, and the ARIMA analysis indicated that the COVID-19 intervention had mitigated EV transmission during the first COVID-19 outbreak in early 2020. In addition, bivariate Spearman's cross-correlation coefficients were estimated for the major types CVA6, CVA16 and EV-A71. Our analyses indicated the possible existence of correlations among CVA6, CVA16 and EV-A71 prevalence in the epidemiological level.

Discussion

Taken together, the type distribution of HFMD has substantially changed over the last decade and CVA6 and CVA16 are currently the most predominant types co-circulating in Nanchang. The anti-COVID-19 NPIs significantly reduced the incidence of EV infections.

Epidemiology of Hand, Foot, and Mouth Disease and Genetic Evolutionary Characteristics of Coxsackievirus A10 in Taiyuan City, Shanxi Province from 2016 to 2020

In recent years, the prevalence of hand, foot, and mouth disease (HFMD) caused by enteroviruses other than enterovirus A71 (EV-A71) and coxsackievirus A16 (CVA16) has gradually increased. The throat swab specimens of 2701 HFMD cases were tested, the VP1 regions of CVA10 RNA were amplified using RT-PCR, and phylogenetic analysis of CVA10 was performed. Children aged 1–5 years accounted for the majority (81.65%) and boys were more than girls. The positivity rates of EV-A71, CVA16, and other EVs were 15.22% (219/1439), 28.77% (414/1439), and 56.01% (806/1439), respectively. CVA10 is one of the important viruses of other EVs. A total of 52 CVA10 strains were used for phylogenetic analysis based on the VP1 region, 31 were from this study, and 21 were downloaded from GenBank. All CVA10 sequences could be assigned to seven genotypes (A, B, C, D, E, F, and G), and genotype C was further divided into C1 and C2 subtypes, only one belonged to subtype C1 and the remaining 30 belonged to C2 in this study. This study emphasized the importance of strengthening the surveillance of HFMD to understand the mechanisms of pathogen variation and evolution, and to provide a scientific basis for HFMD prevention, control, and vaccine development.

Current status of hand-foot-and-mouth disease

Hand-foot-and-mouth disease (HFMD) is a viral illness commonly seen in young children under 5 years of age, characterized by typical manifestations such as oral herpes and rashes on the hands and feet. These symptoms typically resolve spontaneously within a few days without complications. Over the past two decades, our understanding of HFMD has greatly improved and it has received significant attention. A variety of research studies, including epidemiological, animal, and in vitro studies, suggest that the disease may be associated with potentially fatal neurological complications. These findings reveal clinical, epidemiological, pathological, and etiological characteristics that are quite different from initial understandings of the illness. It is important to note that HFMD has been linked to severe cardiopulmonary complications, as well as severe neurological sequelae that can be observed during follow-up. At present, there is no specific pharmaceutical intervention for HFMD. An inactivated Enterovirus A71 (EV-A71) vaccine that has been approved by the China Food and Drug Administration (CFDA) has been shown to provide a high level of protection against EV-A71-related HFMD. However, the simultaneous circulation of multiple pathogens and the evolution of the molecular epidemiology of infectious agents make interventions based solely on a single agent comparatively inadequate. Enteroviruses are highly contagious and have a predilection for the nervous system, particularly in child populations, which contributes to the ongoing outbreak. Given the substantial impact of HFMD around the world, this Review synthesizes the current knowledge of the virology, epidemiology, pathogenesis, therapy, sequelae, and vaccine development of HFMD to improve clinical practices and public health efforts.

Antigenic mapping of enterovirus A71 from Taiwan and Southeast Asia

Enterovirus A71 (EV-A71) is a non-enveloped virus possessing 4 capsid proteins: VP1-VP4. The outermost capsid protein, VP1, plays roles in both antigenicity and virulence of the virus. The concept of generating other EV-A71 genotypes of reverse genetics (rg) viruses by replacing VP1 can be made possible with synthetic biotechnology, allowing us to redesign organisms, creating unavailable ones. To determine suitable vaccine candidates against EV-A71 infections, we combined synthetic biotechnology, rg-virus production and high-fidelity determinants to produce genetically stable viruses. With the use of antigenic cartography, we are able to view the antigenic distance among various points. We analyzed and generated various EV-A71 VP1 sequences from Taiwan and Southeast Asian (SEA) countries, which were then used to produce recombinant rg-viruses and the viral proteins were purified for immunization of mice and rabbits. Antisera against various EV-A71 genotypes were used in neutralization assays against various Taiwan and SEA EV-A71 genotypes. Based on neutralization data from mice and rabbit antisera, we found that antisera produced from several genotypes were able to effectively neutralize the various Taiwan and SEA EV-A71 genotypes. Additionally, comparing the antigenic maps produced from mouse, rabbit and human antisera against different EV-A71 genotypes, a difference in clustering was seen and the spacing between points also differed. Based on antigenic mapping and neutralizing activities, B4 7008-HF and C4 M79 may be good potential vaccine candidates against EV-A71.

Enterovirus A71 Promotes Exosome Secretion by the Nonstructural Protein 3A Interacting with Rab27a

Exosomes are small membrane-bound vesicles which are intraluminal vesicles (ILVs) secreted to the extracellular space after multivesicular bodies (MVBs) fuse with the plasma membrane. Although it is known that exosomes play a multitude of roles during viral infection, the mechanism that regulates their secretion during viral infection is unknown. Here, we found that enterovirus A71 (EV-A71) infection increased exosome secretion both in vivo and in vitro. Importantly, the expression of nonstructural protein 3A was sufficient to promote exosome secretion, while a mutation affecting the amino acid 18 position abrogated this effect, without changing the size of exosomes in vivo or in vitro. Transmission electron microscopy (TEM) analysis revealed that 3A decreases the number of MVBs and ILVs in vivo and in vitro, which suggested 3A may boost the fusion between MVBs and the plasma membrane. Furthermore, we demonstrated that an interaction between 3A and the small GTPase protein, Rab27a, protected Rab27a from ubiquitination, resulted in increasing exosome release. Data indicated a novel mechanism by which EV-A71 3A modifies exosome secretion during viral infection. IMPORTANCE Research has shown that viral infection impacts exosome secretion, but its regulation mechanisms remain poorly understood. Nonstructural protein 3A of EV-A71 interacts with many host factors and is involved in the remodeling of cellular membranes. In this investigation, we applied exogenous expression of 3A protein for exploring its regulation on exosome secretion and utilized immunoprecipitation combined with proteomics approaches to identify 3A-interacting factors. Our results demonstrate that 3A protein upregulates the release of the exosomes and that the 3A mutant strain of EV-A71 induce less exosome release compared with the EV-A71 wild type. Viral 3A protein interacts with the host factor Rab27a to prevent it from being ubiquitinated, which in turn improves exosome secretion both in vitro and in vivo. EV-A71 3A protein is a novel viral factor in the control of exosome production.

Genomic Epidemiology and Phylodynamic Analysis of Enterovirus A71 Reveal Its Transmission Dynamics in Asia

Enterovirus A71 (EV-A71) is one of the main pathogens causing hand, foot, and mouth disease (HFMD) outbreaks in Asian children under 5 years of age. In severe cases, it can cause neurological complications and be life-threatening. In this study, 200 newly sequenced EV-A71 whole-genome sequences were combined with 772 EV-A71 sequences from GenBank for large-scale analysis to investigate global EV-A71 epidemiology, phylogeny, and Bayesian phylodynamic characteristics. Based on the phylogenetic analysis of the EV-A71 3D^pol region, six new evolutionary lineages (lineages B, J, K, O, P, and Q) were found in this study, and the number of evolutionary lineages was expanded from 11 to 17. Temporal dynamics and recombination breakpoint analyses based on genotype C revealed that recombination of nonstructural protein-coding regions, including 3D^pol, is an important reason for the emergence of new lineages. The EV-A71 epidemic in the Asia-Pacific region is complex, and phylogeographic analysis found that Vietnam played a key role in the spread of subgenotypes B5 and C4. The origin of EV-A71 subgenotype C4 in China is East China, which is closely related to the prevalence of subgenotype C4 in the south and throughout China. Selection pressure analysis revealed that, in addition to VP1 amino acid residues VP1-98 and VP1-145, which are associated with EV-A71 pathogenicity, amino acid residues VP1-184 and VP1-249 were also positively selected, and their functions still need to be determined by biology and immunology. This study aimed to provide a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development through a comprehensive analysis. IMPORTANCE EV-A71 is one of the most important pathogens causing HFMD outbreaks; however, large-scale studies of EV-A71 genomic epidemiology are currently lacking. In this study, 200 new EV-A71 whole-genome sequences were determined. Combining these with 772 EV-A71 whole-genome sequences in the GenBank database, the evolutionary and transmission characteristics of global and Asian EV-A71 were analyzed. Six new evolutionary lineages were identified in this study. We also found that recombination in nonstructural protein-coding regions, including 3D^pol, is an important cause for the emergence of new lineages. The results provided a solid theoretical basis for EV-A71-related disease surveillance and prevention, antiviral research, and vaccine development.

Molecular Epidemiology and Evolution of Coxsackievirus A9

Nineteen CVA9 isolates were obtained between 2010 and 2019 from six provinces of mainland China, using the HFMD surveillance network established in China. Nucleotide sequencing revealed that the full-length VP1 of 19 CVA9 isolates was 906 bases encoding 302 amino acids. The combination of the thresholds of the phylogenetic tree and nucleotide divergence of different genotypes within the same serotype led to a value of 15-25%, and enabled CVA9 worldwide to be categorized into ten genotypes: A-J. The phylogenetic tree showed that the prototype strain was included in genotype A, and that the B, C, D, E, H, and J genotypes disappeared during virus evolution, whereas the F, I, and G genotypes showed co-circulation. Lineage G was the dominant genotype of CVA9 and included most of the strains from nine countries in Asia, North America, Oceania, and Europe. Most Chinese strains belonged to the G genotype, suggesting that the molecular epidemiology of China is consistent with that observed worldwide. The 165 partial VP1 strains (723 nt) showed a mean substitution rate of 3.27 × 10-3 substitution/site/year (95% HPD range 2.93-3.6 × 10-3), dating the tMRCA of CVA9 back to approximately 1922 (1911-1932). The spatiotemporal dynamics of CVA9 showed the spread of CVA9 obviously increased in recent years. Most CVA9 isolates originated in USA, but the epidemic areas of CVA9 are now concentrated in the Asia-Pacific region, European countries, and North America. Recombination analysis within the enterovirus B specie (59 serotypes) revealed eight recombination patterns in China at present, CVB4, CVB5, E30, CVB2, E11, HEV106, HEV85, and HEV75. E14, and E6 may act as recombinant donors in multiple regions. Comparison of temperature sensitivity revealed that temperature-insensitive strains have more amino acid substitutions in the RGD motif of the VP1 region, and the sites T283S, V284M, and R288K in the VP1 region may be related to the temperature tolerance of CVA9.

A novel subgenotype C6 Enterovirus A71 originating from the recombination between subgenotypes C4 and C2 strains in mainland China

Recombination plays important roles in the genetic diversity and evolution of Enterovirus A71 (EV-A71). The phylogenetics of EV-A71 in mainland China found that one strain DL71 formed a new subgenotype C6 with unknown origin. This study investigated the detailed genetic characteristics of the new variant. DL71 formed a distinct cluster within genotype C based on the genome and individual genes (5′UTR, VP4, VP1, 2A, 2B, 2C, 3D, and 3′UTR). The average genetic distances of the genome and individual genes (VP3, 2A, 2B, 2C, 3A, 3C, and 3D) between DL71 and reference strains were greater than 0.1. Nine recombination events involving smaller fragments along DL71 genome were detected. The strains Fuyang-0805a (C4) and Tainan/5746/98 (C2) were identified as the parental strains of DL71. In the non-recombination regions, DL71 had higher identities with Fuyang-0805a than Tainan/5746/98, and located in the cluster with C4 strains. However, in the recombination regions, DL71 had higher identities with Tainan/5746/98 than Fuyang-0805a, and located in the cluster with C2 strains. Thus, DL71 was a novel multiple inter-subgenotype recombinant derived from the dominant subgenotype C4 and the sporadic subgenotype C2 strains. Monitoring the emergence of new variants by the whole-genome sequencing remains essential for preventing disease outbreaks and developing new vaccines.

A Novel Attenuated Enterovirus A71 Mutant with VP1-V238A,K244R Exhibits Reduced Efficiency of Cell Entry/Exit and Augmented Binding Affinity to Sulfated Glycans

Enterovirus A71 (EV-A71) is one of the major etiological agents of hand, foot, and mouth disease (HFMD), and infection occasionally leads to fatal neurological complications in children. However, only inactivated whole-virus vaccines against EV-A71 are commercially available in Mainland China. Furthermore, the mechanisms underlying the infectivity and pathogenesis of EV-A71 remain to be better understood. By adaptation of an EV-A71 B5 strain in monkey Vero cells in the presence of brilliant black BN (E151), an anti-EV-A71 agent, a double mutant with VP1-V238A,K244R emerged whose infection was enhanced by E151. The growth of the reverse genetics (RG) mutant RG/B5-VP1-V238A,K244R (RG/B5-AR) was promoted by E151 in Vero cells but inhibited in other human and murine cells, while its parental wild type, RG/B5-wt, was strongly prevented by E151 from infection in all tested cells. In the absence of E151, RG/B5-AR exhibited defective cell entry/exit, resulting in reduced viral transmission and growth in vitro. It had augmented binding affinity to sulfated glycans, cells, and tissue/organs, which probably functioned as decoys to restrict viral dissemination and infection. RG/B5-AR was also attenuated, with a 355 times higher 50% lethal dose (LD₅₀) and a shorter timing of virus clearance than those of RG/B5-wt in suckling AG129 mice. However, it remained highly immunogenic in adult AG129 mice and protected their suckling mice from lethal EV-A71 challenges through maternal neutralizing antibodies. Overall, discovery of the attenuated mutant RG/B5-AR contributes to better understanding of virulence determinants of EV-A71 and to further development of novel vaccines against EV-A71.

IMPORTANCE Enterovirus A71 (EV-A71) is highly contagious in children and has been responsible for thousands of deaths in Asia-Pacific region since the 1990s. Unfortunately, the virulence determinants and pathogenesis of EV-A71 are not fully clear. We discovered that a novel EV-A71 mutant, VP1-V238A,K244R, showed growth attenuation with reduced efficiency of cell entry/exit. In the Vero cell line, which has been approved for manufacturing EV-A71 vaccines, the growth defects of the mutant were compensated by a food dye, brilliant black BN. The mutant also showed augmented binding affinity to sulfated glycans and other cellular components, which probably restricted viral infection and dissemination. Therefore, it was virulence attenuated in a mouse model but still retained its immunogenicity. Our findings suggest the mutant as a promising vaccine candidate against EV-A71 infection.

Characterization of Coxsackievirus A6 Strains Isolated From Children With Hand, Foot, and Mouth Disease

Coxsackievirus A6 (CVA6) is a key pathogen causing hand, foot and mouth disease (HFMD). However, there are currently no specific antiviral drugs or vaccines for treating infections caused by CVA6. In this study, human rhabdomyosarcoma (RD), African green monkey kidney (Vero), and human embryonic lung diploid fibroblast (KMB17) cells were used to isolate CVA6 from 327 anal swab and fecal samples obtained during HFMD monitoring between 2009 and 2017. The VP1 genes of the isolates were sequenced and genotyped, and the biological characteristics of the representative CVA6 strains were analyzed. A total of 37 CVA6 strains of the D3 gene subtypes were isolated from RD cells, all of which belonged to the epidemic strains in mainland China. Using the adaptive culture method, 10 KMB17 cell-adapted strains were obtained; however, no Vero cell-adapted strains were acquired. Among the KMB17 cell-adapted strains, only KYN-A1205 caused disease or partial death in suckling mice, and its virulence was stronger than its RD cell-adapted strain. The pathogenic KYN-A1205 strain caused strong tropism to the muscle tissue and led to pathological changes, including muscle necrosis and nuclear fragmentation in the forelimb and hindlimb. Sequence analysis demonstrated that the KYN-A1205 strain exhibited multiple amino acid mutations after KMB17 cell adaptation. Moreover, it showed strong pathogenicity, good immunogenicity and genetic stability, and could be used as an experimental CVA6 vaccine candidate.

Global spatiotemporal transmission patterns of human enterovirus 71 from 1963 to 2019

Enterovirus 71 (EV71) can cause large outbreaks of hand, foot, and mouth disease (HFMD) and severe neurological diseases, which is regarded as a major threat to public health, especially in Asia-Pacific regions. However, the global spatiotemporal spread of this virus has not been identified. In this study, we used large sequence datasets and a Bayesian phylogenetic approach to compare the molecular epidemiology and geographical spread patterns of different EV71 subgroups globally. The study found that subgroups of HFMD presented global spatiotemporal variation, subgroups B0, B1, and B2 have caused early infections in Europe and America, and then subgroups C1, C2, C3, and C4 replaced B0-B2 as the predominant genotypes, especially in Asia-Pacific countries. The dispersal patterns of genotype B and subgroup C4 showed the complicated routes in Asia and the source might in some Asian countries, while subgroups C1 and C2 displayed more strongly supported pathways globally, especially in Europe. This study found the predominant subgroup of EV71 and its global spatiotemporal transmission patterns, which may be beneficial to reveal the long-term global spatiotemporal transmission patterns of human EV71 and carry out the HFMD vaccine development.

Functional Insights into Silymarin as an Antiviral Agent against Enterovirus A71 (EV-A71)

Enterovirus A71 (EV-A71) is a major neurovirulent agent capable of causing severe hand, foot and mouth disease (HFMD) associated with neurological complications and death. Currently, no FDA-approved antiviral is available for the treatment of EV-A71 infections. The flavonoid silymarin was shown to exert virucidal effects, but the binding site on the capsid was unknown. In this study, the ligand interacting site of silymarin was determined in silico and validated in vitro. Moreover, the potential of EV-A71 to develop resistance against silymarin was further evaluated. Molecular docking of silymarin with the capsid of EV-A71 indicated that silymarin binds to viral protein 1 (VP1) of EV-A71, specifically at the GH loop of VP1. The in vitro binding of silymarin with VP1 of EV-A71 was validated using recombinant VP1 through ELISA competitive binding assay. Continuous passaging of EV-A71 in the presence of silymarin resulted in the emergence of a mutant carrying a substitution of isoleucine by threonine (I97T) at position 97 of the BC loop of EV-A71. The mutation was speculated to overcome the inhibitory effects of silymarin. This study provides functional insights into the underlying mechanism of EV-A71 inhibition by silymarin, but warrants further in vivo evaluation before being developed as a potential therapeutic agent.

First evidence of enterovirus A71 and echovirus 30 in Uruguay and genetic relationship with strains circulating in the South American region

Human enteroviruses (EVs) comprise more than 100 types of coxsackievirus, echovirus, poliovirus and numbered enteroviruses, which are mainly transmitted by the faecal-oral route leading to diverse diseases such as aseptic meningitis, encephalitis, and acute flaccid paralysis, among others. Since enteroviruses are excreted in faeces, wastewater-based epidemiology approaches are useful to describe EV diversity in a community. In Uruguay, knowledge about enteroviruses is extremely limited. This study assessed the diversity of enteroviruses through Illumina next-generation sequencing of VP1-amplicons obtained by RT-PCR directly applied to viral concentrates of 84 wastewater samples collected in Uruguay during 2011-2012 and 2017-2018. Fifty out of the 84 samples were positive for enteroviruses. There were detected 27 different types belonging to Enterovirus A species (CVA2-A6, A10, A16, EV-A71, A90), Enterovirus B species (CVA9, B1-B5, E1, E6, E11, E14, E21, E30) and Enterovirus C species (CVA1, A13, A19, A22, A24, EV-C99). Enterovirus A71 (EV-A71) and echovirus 30 (E30) strains were studied more in depth through phylogenetic analysis, together with some strains previously detected by us in Argentina. Results unveiled that EV-A71 sub-genogroup C2 circulates in both countries at least since 2011-2012, and that the C1-like emerging variant recently entered in Argentina. We also confirmed the circulation of echovirus 30 genotypes E and F in Argentina, and reported the detection of genotype E in Uruguay. To the best of our knowledge this is the first report of the EV-A71 C1-like emerging variant in South-America, and the first report of EV-A71 and E30 in Uruguay.

Cross-Antigenicity between EV71 Sub-Genotypes: Implications for Vaccine Efficacy

Enterovirus A-71 (EV71) is a global, highly contagkkious pathogen responsible for severe cases of hand-food-mouth-disease (HFMD). The use of vaccines eliciting cross neutralizing antibodies (NTAbs) against the different circulating EV71 sub-genotypes is important for preventing HFMD outbreaks. Here, we tested the cross-neutralizing activities induced by EV71 genotype/sub-genotype A, B0-B4, C1, C2, C4, and C5 viruses using rats. Differences were noted in the cross-neutralization of the 10 sub-genotypes tested but there were generally good levels of cross-neutralization except against genotype A virus, against which neutralization antibody titres (NTAb) where the lowest with NTAbs being the highest against sub-genotype B4. Moreover, NTAb responses induced by C4, B4, C1, and C2 viruses were homogenous, with values of maximum/minimum NTAb ratios (MAX/MIN) against all B and C viruses ranging between 4.0 and 6.0, whereas MAX/MIN values against B3 and A viruses were highly variable, 48.0 and 256.0, respectively. We then dissected the cross-neutralizing ability of sera from infants and children and rats immunized with C4 EV71 vaccines. Cross-neutralizing titers against the 10 sub-genotypes were good in both vaccinated infants and children and rats with the MAX/MIN ranging from 1.8–3.4 and 5.1–7.1, respectively, which were similar to those found in naturally infected patients (2.8). Therefore, we conclude that C4 EV71 vaccines can provide global protection to infants and children against HFMD caused by different sub-genotypes.

Safety and Immunogenicity of a Stable, Cold-Adapted, Temperature-Sensitive/Conditional Lethal Enterovirus A71 in Monkey Study

Enterovirus A71 (EV-A71) and coxsackievirus A16 (CA16) are major etiological agents of hand foot and mouth disease (HFMD) in children, which may result in fatal neurological complications. The development of safe, cost effective vaccines against HFMD, especially for use in developing countries, is still a top public health priority. We have successfully generated a stable, cold-adapted, temperature sensitive/conditional lethal EV-A71 through adaptive culturing in Vero cells at incrementally lower cultivation temperatures. An additional 40 passages at an incubation temperature of 28 °C, and a temperature reversion study at an incubation temperature of 37 °C and 39.5 °C, reveals the virus’s phenotypic and genetic stability at the predefined culture conditions. Six unique mutations (two in noncoding regions and four in nonstructural protein-coding genes) in combination may have contributed to its stable phenotype and inability to fully revert to its original wild phenotype. The safety and immunogenicity of this stable, cold-adapted, temperature sensitive/conditional lethal EV-A71 was performed in six monkeys. None of the inoculated monkeys developed any obvious clinical illness except one which developed a transient spike of fever. No gross postmortem lesion or abnormal histological finding was noted for all monkeys at autopsy. No virus was reisolated although EV-A71 specific RNA was detected in serum samples collected on both day 4 and day 8 postinoculation. Only EV-A71 RNA and viral antigen were detected in the spleen homogenate and peripheral blood mononuclear cells, respectively, collected on day 4. The two remaining monkeys developed good humoral immune response on day 14 and day 30 post-inoculation.

Molecular epidemiological characteristics of echovirus 6 in mainland China: extensive circulation of genotype F from 2007 to 2018

Echovirus 6 (E6) is associated with various clinical diseases and is frequently detected in environmental sewage. Despite its high prevalence in humans and the environment, little is known about its molecular phylogeography in mainland China. In this study, 114 of 21,539 (0.53%) clinical specimens from hand, foot, and mouth disease (HFMD) cases collected between 2007 and 2018 were positive for E6. The complete VP1 sequences of 87 representative E6 strains, including 24 strains from this study, were used to investigate the evolutionary genetic characteristics and geographical spread of E6 strains. Phylogenetic analysis based on VP1 nucleotide sequence divergence showed that, globally, E6 strains can be grouped into six genotypes, designated A to F. Chinese E6 strains collected between 1988 and 2018 were found to belong to genotypes C, E, and F, with genotype F being predominant from 2007 to 2018. There was no significant difference in the geographical distribution of each genotype. The evolutionary rate of E6 was estimated to be 3.631 × 10^-3 substitutions site^-1 year^-1 (95% highest posterior density [HPD]: 3.2406 × 10^-3-4.031 × 10^-3 substitutions site^-1 year^-1) by Bayesian MCMC analysis. The most recent common ancestor of the E6 genotypes was traced back to 1863, whereas their common ancestor in China was traced back to around 1962. A small genetic shift was detected in the Chinese E6 population size in 2009 according to Bayesian skyline analysis, which indicated that there might have been an epidemic around that year.

Natural intertypic and intratypic recombinants of enterovirus 71 from mainland China during 2009-2018: a complete genome analysis

Surveillance of recombinant enterovirus 71 (EV71) and subgenotype replacement is vital for preventing and controlling hand, foot, and mouth disease (HFMD) outbreaks. Despite this, data on recombinant variants and phylogeny of circulating EV71 strains in mainland China are limited. In this study, recombinant variants of EV71 were identified in mainland China from 2009 to 2018. Phylogenetic analysis indicated that except for individual strains (CQ2014-86/CQ/CHN/2014 and EV71/Xiamen/2009 (B5)), almost all of the EV71 strains in mainland China belonged to the subgenotype C4a. Analysing complete genome sequences of 196 EV71 isolates, 3 intertypic recombination strains (VR1432, 30-2/2015/BJ, and Guangdong-2009) and 5 intratypic recombination strains (EV71/P1034/2013, VR1432, Henan-ZMD/CHN/2012, Hubei-WH/CHN/2012, and EV71/P868/2013/China) were identified among naturally circulating EV71. The breakpoints of these recombinant strains were located within the P1, P2, and P3 encoding regions. Notably, a double recombinant (VR1432) resulting from recombination between EV71 subgenotype C4a and C4b strain SHZH98 and a CA8 strain Donovan was identified. This study reports these specific intertypic and intratypic recombination events for the first time highlighting the importance of genetic recombination in the emergence of new enterovirus variants.

Pan-Enterovirus Amplicon-Based High-Throughput Sequencing Detects the Complete Capsid of a EVA71 Genotype C1 Variant via Wastewater-Based Epidemiology in Arizona

We describe the complete capsid of a genotype C1-like Enterovirus A71 variant recovered from wastewater in a neighborhood in the greater Tempe, Arizona area (Southwest United States) in May 2020 using a pan-enterovirus amplicon-based high-throughput sequencing strategy. The variant seems to have been circulating for over two years, but its sequence has not been documented in that period. As the SARS-CoV-2 pandemic has resulted in changes in health-seeking behavior and overwhelmed pathogen diagnostics, our findings highlight the importance of wastewater-based epidemiology (WBE ) as an early warning system for virus surveillance.

Global Spread of the B5 Subgenotype EV-A71 and the Phylogeographical Analysis of Chinese Migration Events

The subgenotype B5 of EV-A71 is a widely circulating subgenotype that frequently spreads across the globe. Several outbreaks have occurred in nations, such as Malaysia, Thailand, Vietnam, and Japan. Appearing first in Taiwan, China, the subgenotype has been frequently reported in mainland of China even though no outbreaks have been reported so far. The current study reconstructed the migration of the B5 subgenotype of EV-A71 in China via phylogeographical analysis. Furthermore, we investigated its population dynamics in order to draw more credible inferences. Following a dataset cleanup of B5 subgenotype of EV-A71, we detected earlier B5 subgenotypes of EV-A71 sequences that had been circulating in Malaysia and Singapore since the year 2000, which was before the 2003 outbreak that occurred in Sarawak. The Bayesian inference indicated that the most recent common ancestor of B5 subgenotype EV-A71 appeared in September, 1994 (1994.75). With respect to the overall prevalence, geographical reconstruction revealed that the B5 subgenotype EV-A71 originated singly from single-source cluster and subsequently developed several active lineages. Based on a large amount of data that was accumulated, we conclude that the appearance of the B5 subgenotype of EV-A71 in mainland of China was mainly due to multiple migrations from different origins.

A Polysaccharide Purified From Ganoderma lucidum Acts as a Potent Mucosal Adjuvant That Promotes Protective Immunity Against the Lethal Challenge With Enterovirus A71

Enterovirus A71 (EV-A71), the pathogen responsible for the seasonal hand-foot-and-mouth epidemics, can cause significant mortality in infants and young children. The vaccine against EV-A71 could potentially prevent virus-induced neurological complications and mortalities occurring due to the high risk of poliomyelitis-like paralysis and fatal encephalitis. It is known that polysaccharide purified from Ganoderma lucidum (PS-G) can effectively modulate immune function. Here, we used PS-G as an adjuvant with the EV-A71 mucosal vaccine and studied its effects. Our data showed that PS-G-adjuvanted EV-A71 generated significantly better IgA and IgG in the serum, saliva, nasal wash, bronchoalveolar lavage fluid (BALF), and feces. More importantly, these antibodies could neutralize the infectivity of EV-A71 (C2 genotype) and cross-neutralize the B4, B5, and C4 genotypes of EV-A71. In addition, more EV-A71-specific IgA- and IgG- secreting cells were observed with the used of a combination of EV-A71 and PS-G. Furthermore, T-cell proliferative responses and IFN-γ and IL-17 secretions levels were notably increased in splenocytes when the EV-A71 vaccine contained PS-G. We also found that levels of IFN-γ and IL-17 released in Peyer's patch cells were significantly increased in EV-A71, after it was combined with PS-G. We further demonstrated that both CD4⁺ and CD8⁺ T cells were able to generate IFN-γ and IL-17 in the spleen. An easy-accessed model of hybrid hSCARB2⁺/⁺/stat-1^-/^- mice was used for EV-A71 infection and pathogenesis. We infected the mouse model with EV-A71, which was premixed with mouse sera immunized with the EV-A71 vaccine with or without PS-G. Indeed, in the EV-A71 + PS-G group, the levels of VP1-specific RNA sequences in the brain, spinal cord, and muscle decreased significantly. Finally, hSCARB2-Tg mice immunized via the intranasal route with the PS-G-adjuvanted EV-A71 vaccine resisted a subsequent lethal oral EV-A71 challenge. Taken together, these results demonstrated that PS-G could potentially be used as an adjuvant for the EV-A71 mucosal vaccine.

Genetic Diversity Analysis of Coxsackievirus A8 Circulating in China and Worldwide Reveals a Highly Divergent Genotype

Coxsackievirus A8 (CV-A8) is one of the pathogens associated with hand, foot and mouth disease (HFMD) and herpangina (HA), occasionally leading to severe neurological disorders such as acute flaccid paralysis (AFP). Only one study aimed at CV-A8 has been published to date, and only 12 whole-genome sequences are publicly available. In this study, complete genome sequences from 11 CV-A8 strains isolated from HFMD patients in extensive regions from China between 2013 and 2018 were determined, and all sequences from GenBank were retrieved. A phylogenetic analysis based on a total of 34 complete VP1 sequences of CV-A8 revealed five genotypes: A, B, C, D and E. The newly emerging genotype E presented a highly phylogenetic divergence compared with the other genotypes and was composed of the majority of the strains sequenced in this study. Markov chain Monte Carlo (MCMC) analysis revealed that genotype E has been evolving for nearly a century and somehow arose in approximately 2010. The Bayesian skyline plot showed that the population size of CV-A8 has experienced three dynamic fluctuations since 2001. Amino acid residues of VP1100N, 103Y, 240T and 241V, which were embedded in the potential capsid loops of genotype E, might enhance genotype E adaption to the human hosts. The CV-A8 whole genomes displayed significant intra-genotypic genetic diversity in the non-capsid region, and a total of six recombinant lineages were detected. The Chinese viruses from genotype E might have emerged recently from recombining with European CV-A6 strains. CV-A8 is a less important HFMD pathogen, and the capsid gene diversity and non-capsid recombination variety observed in CV-A8 strains indicated that the constant generation of deleterious genomes and a constant selection pressure against these deleterious mutations is still ongoing within CV-A8 quasispecies. It is possible that CV-A8 could become an important pathogen in the HFMD spectrum in the future. Further surveillance of CV-A8 is greatly needed.

Emergence of genotype C1 Enterovirus A71 and its link with antigenic variation of virus in Taiwan

An outbreak of the hand-foot-mouth disease with severe neurological cases, mainly caused by the genotype C1 enterovirus A71 (EV-A71), occurred in Taiwan between 2018 and early 2019. In the recent decade, the most dominant EV-A71 genotypes in Taiwan were B5 and C4 but changed to C1 in 2018. Antibody-mediated immunity plays a key role in limiting the EV-A71 illness in humans. However, the level of neutralizing activities against genotype C1 virus by human polyclonal and monoclonal antibodies (MAbs) remains largely unclear. In the study, we demonstrated that that 39% (9 in 23) of post-infection sera from the genotype B5- or C4-infected patients in 2014–2017 exhibit reduced titers with the 2018–2019 genotype C1 viruses than with the earlier B5 and C4 viruses tested. This finding with polyclonal sera is confirmed with human MAbs derived from genotype B5 virus-infected individuals. The 2018–2019 genotype C1 virus is resistant to the majority of canyon-targeting human MAbs, which may be associated with the residue change near or at the bottom of the canyon region on the viral capsid. The remaining three antibodies (16-2-11B, 16-3-4D, and 17-1-12A), which target VP1 S241 on the 5-fold vertex, VP3 E81 on the 3-fold plateau and VP2 D84 on the 2-fold plateau of genotype C1 viral capsid, respectively, retained neutralizing activities with variable potencies. These neutralizing antibodies were also found to be protective against a lethal challenge of the 2018–2019 genotype C1 virus in an hSCARB2-transgenic mice model. These results indicate that the EV-A71-specific antibody response may consist of a fraction of poorly neutralizing antibodies against 2018–2019 genotype C1 viruses among a subset of previously infected individuals. Epitope mapping of protective antibodies that recognize the emerging genotype C1 virus has implications for anti-EV-A71 MAbs and the vaccine field.

EV-A71 is a cause of hand-foot-mouth disease, epidemics of which still regularly occur around the globe. Given that EV-A71 immune protection from the disease correlates with neutralizing antibody responses, but the responses in humans prior to an outbreak are still poorly understood. An outbreak of hand-foot-mouth disease among children emerged in Taiwan from 2018 to 2019, and genotype C1 EV-A71 caused most of the cases. Here, we characterized EV-A71-neutralizing antibody profiles in details at both the serological and monoclonal levels and showed that antibodies generated by humans prior to the emergence of genotype C1 EV-A71 less effectively neutralize C1 compared to the prior circulating genotypes, which implies the presence of antigenic variation in the EV-A71 genotypes. We further identified and mapped critical neutralizing epitopes of 2018–2019 genotype C1 EV-A71 on the top and margin of the viral capsid pentamer and demonstrated the in vivo protective effect of human monoclonal antibodies, which highlight the properties of human antibody-neutralizing sites on EV-A71 and the potential of human antibodies as antiviral agents.

Evaluation of the diagnostic performance and its associated factors of a commercial anti-EV-A71 IgM-capture ELISA kit in hospitalized children with clinical diagnostic HFMD

OBJECTIVES

Enterovirus A71 (EV-A71) is the main pathogen of severe hand, foot, and mouth disease (HFMD). Commercial enzyme-linked immunosorbent assays (ELISAs) are widely used in Chinese hospitals for the rapid diagnosis of acute EV-A71 infections. We present an evaluation of the diagnostic performance of a commercial anti-EV-A71 IgM-capture ELISA kit.

METHODS

A prospective, hospital-based HFMD cohort was established in Henan Children's Hospital (February 2017 - February 2018). Stool and blood specimens were collected from 1413 participants for diagnosing EVA71 by quantitative Real-Time Reverse Transcription-Polymerase Chain Reaction (qRT-PCR) and anti-EV-A71 ELISA.

RESULTS

Detection yields of EV-A71 IgM increased from 6.5 % (95 % CI:3.3 %-11.4 %) at 0∼24 h, to 42 % (95 % CI:28.3 %-57.8) at 120∼144 h from onset to sampling, and stabilized at ∼40 % after 144 h. With increased time from onset to sampling, the sensitivity of the commercial ELISA increased from 0.54 (95 % CI:0.25-0.81) to 0.74 (95 % CI:0.43-0.66), while specificity decreased from 0.97 (95 % CI:0.93-0.99) to 0.80 (95 % CI:0.69-0.89), and PPV decreased from 0.96 (95 % CI:0.92-0.99) to 0.84 (95 % CI:0.73-0.92). Multivariate analysis found age, EV-A71 vaccination, previous HFMD/Herpangina infection, disease severity, infection during peak EV-A71 season, and sampling time after symptom onset were significantly associated with the diagnostic performance of this anti-EV-A71 IgM-capture ELISA.

CONCLUSION

Achieving satisfactory specificity and sensitivity scores, this commercial anti-EV-A71 IgM-capture ELISA kit is suitable for clinical EV-A71 diagnosis, particularly in resource-poor areas. However, clinicians should interpret results in the context of patient history and epidemiological setting.

Sub-genotype change and recombination of coxsackievirus A6s may be the cause of it being the predominant pathogen for HFMD in children in Beijing, as revealed by analysis of complete genome sequences

OBJECTIVES

To investigate why coxsackievirus A6 (CVA6) has replaced enterovirus A71 (EV71) and coxsackievirus A16 (CVA16), which used to be the most predominant etiological agents, for hand, foot and mouth disease (HFMD) in children in Beijing, China.

METHODS

Sixty-four CVA6-positive samples were identified from 2010 to 2016 and selected for whole genome sequence amplification and analysis.

RESULTS

It was demonstrated that the whole genome sequences of CVA6s in this study were 7432-7435 nucleotides in length, and the different lengths were only in the 5'UTR region. The phylogenetic tree analysis of the full-length VP1 region of CVA6s indicated that the prevalent CVA6s in Beijing changed from the previous D2 sub-genotype to the D3 sub-genotype in 2013. In this study, two recombinant forms (RFs)- RF-C and RF-D - of CVA6 mainly appeared in 2010 and 2011. Since 2013, three recombinant CVA6 variants - RF-A, J and L - have been prevalent in children with HFMD in Beijing. The recombination region of RF-J was located at the 2C region, while RF-L had a new recombination point in the 3D region. The recombination of prevalent CVA6s in Beijing from 2013 to 2016 occurred within non-capsid regions of the genome, especially the P2 and P3 regions.

CONCLUSIONS

The sub-genotype change and recombination of CVA6s indicated from this study may explain why CVA6 has become the predominant pathogen causing HFMD since 2013.

Characterization of Plaque Variants and the Involvement of Quasi-Species in a Population of EV-A71

Viral plaque morphologies in human cell lines are markers for growth capability and they have been used to assess the viral fitness and selection of attenuated mutants for live-attenuated vaccine development. In this study, we investigate whether the naturally occurring plaque size variation reflects the virulence of the variants of EV-A71. Variants of two different plaque sizes (big and small) from EV-A71 sub-genotype B4 strain 41 were characterized. The plaque variants displayed different in vitro growth kinetics compared to the parental wild type. The plaque variants showed specific mutations being present in each variant strain. The big plaque variants showed four mutations I97L, N104S, S246P and N282D in the VP1 while the small plaque variants showed I97T, N237T and T292A in the VP1. No other mutations were detected in the whole genome of the two variants. The variants showed stable homogenous small plaques and big plaques, respectively, when re-infected in rhabdomyosarcoma (RD) and Vero cells. The parental strain showed faster growth kinetics and had higher viral RNA copy number than both the big and small plaque variants. Homology modelling shows that both plaque variants have differences in the structure of the VP1 protein due to the presence of unique spontaneous mutations found in each plaque variant This study suggests that the EV-A71 sub-genotype B4 strain 41 has at least two variants with different plaque morphologies. These differences were likely due to the presence of spontaneous mutations that are unique to each of the plaque variants. The ability to maintain the respective plaque morphology upon passaging indicates the presence of quasi-species in the parental population.

A Large-Scale Outbreak of Echovirus 30 in Gansu Province of China in 2015 and Its Phylodynamic Characterization

Background

Echovirus 30 (E-30) has been investigated and reported worldwide and is closely associated with several infectious diseases, including encephalitis; myocarditis; and hand, foot, and mouth disease. Although many E-30 outbreaks associated with encephalitis have been reported around the world, it was not reported in northwest China until 2015.

Methods

The clinical samples, including the feces, serum, throat swabs, and cerebrospinal fluid, were collected for this study and were analyzed for diagnosis. E-30 was isolated and processed according to the standard procedures. The epidemiological and phylogenetic analysis were performed to indicate the characteristics of E-30 outbreaks and phylodynamics of E-30 in China.

Results

The E-30 outbreaks affected nine towns of Gansu Province in 2015, starting at a school of Nancha town and spreading to other towns within 1 month. The epidemiological features showed that children aged 6–15 years were more susceptible to E-30 infection. The genotypes B and C cocirculated in the world, whereas the latter dominated the circulation of E-30 in China. The genome sequences of this outbreak present 99.3–100% similarity among these strains, indicating a genetic-linked aggregate outbreak of E-30 in this study. Two larger genetic diversity expansions and three small fluctuations of E-30 were observed from 1987 to 2016 in China, which revealed the oscillating patterns of E-30 in China. In addition, the coastal provinces of China, such as Zhejiang, Fujian, and Shandong, were initially infected, followed by other parts of the country. The E-30 strains isolated from mainland of China may have originated from Taiwan of China in the last century.

Conclusion

The highly similar E-30 genomes in this outbreak showed an aggregate outbreak of E-30, with nine towns affected. Our results suggested that, although the genetic diversity of E-30 oscillates, the dominant lineages of E-30 in China has complex genetic transmission. The coastal provinces played an important role in E-30 spread, which implied further development of effective countermeasures. This study provides a further insight into the E-30 outbreak and transmission and illustrates the importance of valuable surveillance in the future.

Newly emerged enterovirus-A71 C4 sublineage may be more virulent than B5 in the 2015-2016 hand-foot-and-mouth disease outbreak in northern Vietnam

Enterovirus-A71 (EV-A71) is a common cause of hand-foot-and-mouth disease (HFMD) and, rarely, causes severe neurological disease. This study aimed to elucidate the epidemiological and genetic characteristics and virulence of EV-A71 strains isolated from children diagnosed with HFMD. Rectal and throat swabs were collected from 488 children with HFMD in Hanoi, Vietnam, in 2015-2016. From 391 EV-positive patients, 15 EVs, including coxsackievirus A6 (CV-A6; 47.1%) and EV-A71 (32.5%, n = 127), were identified. Of the 127 EV-A71 strains, 117 (92.1%) were the B5 subgenotype and 10 (7.9%) were the C4 subgenotype. A whole-genome analysis of EV-A71 strains showed that seven of the eight C4a strains isolated in 2016 formed a new lineage, including two possible recombinants between EV-A71 C4 and CV-A8. The proportion of inpatients among C4-infected children was higher than among B5-infected children (80.0% vs. 27.4%; P = 0.002). The virulence of EV-A71 strains was examined in human scavenger receptor class B2 (hSCARB2)-transgenic mice, and EV-A71 C4 strains exhibited higher mortality than B5 strains (80.0% vs. 30.0%, P = 0.0001). Thus, a new EV-A71 C4a-lineage, including two possible recombinants between EV-A71 C4 and CV-A8, appeared in 2016 in Vietnam. The EV-A71 C4 subgenotype may be more virulent than the B5 subgenotype.

Recombination Events and Conserved Nature of Receptor Binding Motifs in Coxsackievirus A9 Isolates

Coxsackievirus A9 (CVA9) is an enterically transmitted enterovirus and one of the most pathogenic type among human enteroviruses. CVA9 isolates use a distinctive RGD (Arg-Gly-Asp) motif within VP1 capsid protein that defines its ability to bind to integrin receptor(s) for cellular entry. To investigate CVA9 evolution and pathogenicity, genetic relationships and recombination events were analyzed between 54 novel clinical isolates of CVA9, as well as 21 previously published full length CVA9 sequences from GenBank. Samples were investigated by partial sequencing of the novel VP1 and 3Dpol genes, as well as including the corresponding areas from GenBank sequences. Phylogenetic analyses were combined with clinical data in a further attempt to analyze whether sequence evolution reflects CVA9 pathogenicity in the phylogenies. Furthermore, VP1 gene was also analyzed for receptor binding sites including the RGD motif and the putative heparan sulfate (HS) site. Analysis of the 559-nucleotide-long VP1 sequences identified six clades. Although most of the strains within each clade showed geographical clustering, the grouping pattern of the isolates in the analysis of the VP1 gene was strikingly different from grouping of 3Dpol, which suggests that recombination events may have occurred in the region encoding the nonstructural proteins. Inclusion of clinical data did not provide any evidence of symptom based phylogenetic clustering of CVA9 isolates. Amino acid sequence analysis of the VP1 polypeptide demonstrated that the RGD motif was fully conserved among the isolates while the putative HS binding site was only found in one isolate. These data suggest that integrin binding is essential for virus tropism, but do not explain the symptom repertoire.

Establishment of Asia-Pacific Network for Enterovirus Surveillance

Enteroviruses (EV), the major pathogens of hand, foot, and mouth disease (HFMD) and herpangina, affect millions of children each year. Most human enteroviruses cause self-limited infections except polioviruses, enterovirus A71 (EV-A71), enterovirus D68 (EV-D68), and several echoviruses (Echo) and coxsackieviruses (CV). Especially, EV-A71 has repeatedly caused large-scale outbreaks in the Asia-Pacific region since 1997. Some Asian countries have experienced cyclical outbreaks of severe EV-A71 infections and initiated development of EV-A71 vaccines. Five EV-A71 vaccine candidates have been clinically evaluated and three of them were approved for marketing in China. However, none of the China-approved products seek marketing approval in other countries. This situation supports a role for collaboration among Asian countries to facilitate clinical trials and licensure of EV-A71 vaccines. Additionally, enterovirus D68 outbreaks have been reported in the US and Taiwan currently and caused severe complications and deaths. Hence, an Asia-Pacific Network for Enterovirus Surveillance (APNES) has been established to estimate disease burden, understand virus evolution, and facilitate vaccine development through harmonizing laboratory diagnosis and data collection. Founded in 2017, the APNES is comprised of internationally recognized experts in the field of enterovirus in Asian countries working to raise awareness of this potentially fatal and debilitating disease. This article demonstrated the summaries of the first expert meeting, 2017 International Workshop on Enterovirus Surveillance and Vaccine Development, held by APNES in Taipei, Taiwan, March 2017.

Enterovirus A71: virulence, antigenicity, and genetic evolution over the years

As a neurotropic virus, enterovirus A71 (EV-A71) emerge and remerge in the Asia-Pacific region since the 1990s, and has continuously been a threat to global public health, especially in children. Annually, EV-A71 results in hand-foot-and-mouth disease (HFMD) and occasionally causes severe neurological disease. Here we reviewed the global epidemiology and genotypic evolution of EV-A71 since 1997. The natural selection, mutation and recombination events observed in the genetic evolution were described. In addition, we have updated the antigenicity and virulence determinants that are known to date. Understanding EV-A71 epidemiology, genetic evolution, antigenicity, and virulence determinants can expand our insights of EV-A71 pathogenesis, which may benefit us in the future.

The History of Enterovirus A71 Outbreaks and Molecular Epidemiology in the Asia-Pacific Region

Enterovirus A71 (EV-A71) is one of the common causative pathogens for hand foot and mouth disease (HFMD) affecting young children. HFMD outbreak can result in a substantial pediatric hospitalization and burden the healthcare services, especially in less-developed countries. Since the initial epidemic of predominantly EV-A71 in California in 1969, the high prevalence of HFMD in the Asia-pacific region and elsewhere around the world represents a significant morbidity in this age group. With the advent of rapid and accurate diagnostic tools, there has been a dramatic increase in the number of laboratory-confirmed EV-A71 infection over the past two decades. The population, cultural, and socioeconomic diversity among countries in the Asia-Pacific region all influence the transmission and morbidity associated with HFMD. This review summarizes the current state of epidemiology of EV-A71 in Asia-Pacific countries based on the most recent epidemiological data and available information on the prevalence and disease burden. This knowledge is important in guiding the prevention, control and future research on vaccine development of this highly contagious disease of significant socioeconomic implications in public health.

Severe enterovirus A71 associated hand, foot and mouth disease, Vietnam, 2018: preliminary report of an impending outbreak

Since January 2018, over 53,000 hospitalisations and six deaths due to hand, foot and mouth disease (HFMD) have occurred across Vietnam with most cases from September onward. In a large tertiary referral hospital, Ho Chi Minh City, enterovirus A71 subgenogroup C4 was predominant, while B5 was only sporadically detected. The re-emergence of C4 after causing a severe HFMD outbreak with > 200 deaths in 2011–12 among susceptible young children raises concern of another impending severe outbreak.

Multiple genotypes of Echovirus 11 circulated in mainland China between 1994 and 2017

Echovirus 11 (E-11) is one of the most frequently isolated enteroviruses causing meningitis and other diseases such as hand, foot, and mouth disease (HFMD) and acute flaccid paralysis (AFP). Fifty-nine newly determined E-11 VP1 sequences from the China AFP and HFMD surveillance network and 500 E-11 VP1 sequences obtained from the GenBank database, which were associated with 12 categories of diseases, were screened for phylogenetic analysis. Based on the standard method of genotype classification, E-11 strains circulated worldwide were reclassified into six genotypes as A, B, C, D, E, and F, in which genotype F is newly divided, and genotypes A and C are further divided into A1–5 and C1–4 by this research, whereas genotype D was still divided into D1–5 as in a previous study of Oberste et al. Sub-genotype A1 was the predominant sub-genotype in mainland China between 2008–2017, whereas sub-genotype D5 was the predominant sub-genotype circulated outside China from 1998–2014. However, genotype and sub-genotype spectra showed statistical significance among AFP and HFMD cases (χ² = 60.86, P < 0.001), suggesting that different genotypes might have a tendency to cause different diseases. Strengthening the surveillance of E-11 might provide further information about pathogenic evolution or specific nucleotide mutation associated with different clinical diseases.

Chemical components from the twigs of Caesalpinia latisiliqua and their antiviral activity

Three new compounds, (3S)-dihydrobonducellin 8-O-β-D-glucopyranoside (1), 3',5'-dimethoxy-jezonolid (2), and latisilinoid (3), along with 16 known compounds, were isolated from the twigs of Caesalpinia latisiliqua (Leguminosae). The known compounds were identified as flavonoids, stilbenes, and phenolics as determined by extensive spectroscopic methods, including 1D and 2D NMR. All the isolated compounds were evaluated for their antiviral activity in HRV1B-, CVB3-, and EV71-infected cells. Among the tested compounds, three flavonoids (4-6) and two stilbenes (12 and 14) exhibited significant antiviral activity. This is the first phytochemical investigation of C. latisiliqua twigs.

Epidemiological study on hand, foot and mouth disease in Tongzhou District, Beijing, 2013-2017

Objective

To study the epidemiological characteristics of hand, foot and mouth disease (HFMD) in Tongzhou District, Beijing between 2013 and 2017.

Methods

Data on HFMD infections from 1 January 2013 to 31 December 2017 were collected from the Notifiable Infectious Diseases Reporting Information System and analysed. Serotyping of enteroviruses from samples from patients with HFMD was undertaken using reverse transcription–polymerase chain reaction.

Results

A total of 15 341 patients with HFMD were reported and 32 patients (0.2%) were classified as having severe HFMD. The annual mean incidence rate of HFMD was 219.3/100 000 of the general population. The incidence and case-severity rates of HFMD generally decreased between 2013 and 2017. In the floating migrant population, the incidence and cases-severity rates of HFMD were significantly higher than in the local population. The peak incidence and severity-case rates were at 2 years of age and > 90% of patients were ≤5 years. Enterovirus A71 and Coxsackievirus A16 were the predominant pathogens in 2013–2017.

Conclusions

During the 5-year period 2013–2017, the incidence rate and case-severity rate of HFMD generally decreased in Tongzhou District, Beijing. The floating migrant population and children ≤5 years of age were at the highest risk of HFMD.

Replication and cytokine profiles of different subgenotypes of enterovirus 71 isolated from Thai patients in peripheral blood mononuclear cells

Enterovirus 71 (EV71) and coxsackievirus A16 (CA16) are common causative agents of mild and self-limiting symptoms of childhood hand, foot, and mouth disease (HFMD). However, some EV71-infected HFMD patients can develop severe neurological and/or fatal cardiopulmonary complications. In Thailand, HFMD associated with the EV71 subgenotypes C4a and B5 were reported to be associated with diverse outcomes. However, variations in enterovirus subgenotypes and virulence factors have not been fully elucidated; this study elucidated these variations in peripheral blood mononuclear cells (PBMCs) exposed to different subgenotypes of isolated enteroviruses for 24 and 48 h. Following infection, viral titers were determined by plaque assay. Infected cells and intracellular cytokines were quantified using flow cytometry, and multiplex assay was used to examine cytokine release. All isolated subgenotypes showed replication capability in PBMCs; specifically, the replication titer of EV71 C4a tended to be higher than titers of EV71 B5 and CA16. Additionally, the infectivity of EV71 B5 was higher in monocytes than in lymphocytes. Compared with EV71 B5, EV71 C4a and CA16 had greater ability to induce intra- and extracellular cytokine responses. These findings provide new insights into variations in cellular immune responses to different EV71 subgenotypes isolated from Thai patients, which should be considered for the development of vaccines and therapeutic agents.

Molecular characterization of enterovirus-A71 in children with acute flaccid paralysis in the Philippines

Background

Several inactivated enterovirus-A71 (EV-A71) vaccines are currently licensed in China; however, the development of additional EV-A71 vaccines is ongoing, necessitating extensive analysis of the molecular epidemiology of the virus worldwide. Until 2012, laboratory confirmation of EV-A71 for hand, foot, and mouth disease (HFMD) and other associated diseases had not occurred in the Philippines. Because EV-A71 has been linked with cases of acute flaccid paralysis (AFP), AFP surveillance is one strategy for documenting its possible circulation in the country. To expand current knowledge on EV-A71, molecular epidemiologic analysis and genetic characterization of EV-A71 isolates were performed in this study.

Methods

A retrospective study was performed to identify and characterize nonpolio enteroviruses (NPEVs) associated with AFP in the Philippines, and nine samples were found to be EV-A71–positive. Following characterization of these EV-A71 isolates, the complete viral protein 1 (VP1) gene was targeted for phylogenetic analysis.

Results

Nine EV-A71 isolates detected in 2000 (n = 2), 2002 (n = 4), 2005 (n = 2), and 2010 (n = 1) were characterized using molecular methods. Genomic regions spanning the complete VP1 region were amplified and sequenced using specific primers. Phylogenetic analysis of the full-length VP1 region identified all nine EV-A71 Philippine isolates as belonging to the genogroup C lineage, specifically the C2 cluster. The result indicated a genetic linkage with several strains isolated in Japan and Taiwan, suggesting that strains in the C2 cluster identified in the Asia-Pacific region were circulating in the Philippines.

Conclusion

The study presents the genetic analysis of EV-A71 in the Philippines. Despite some limitations, the study provides additional genetic data on the circulating EV-A71 strains in the Asia-Pacific region, in which information on EV-A71 molecular epidemiology is incomplete. Considering that EV-A71 has a significant public health impact in the region, knowledge of its circulation in each country is important, especially for formulating vaccines covering a wide variety of strains.

Molecular characteristic analysis for the VP1 region of coxsackievirus A6 strains isolated in Jiujiang area, China, from 2012 to 2013

BACKGROUND

Enterovirus 71 (EV-A71) and Coxsackievirus A16 (CV-A16) are the most common causative agents causing hand, foot, and mouth disease (HFMD). However, coxsackievirus A6 (CV-A6), previously largely ignored, became the predominant pathogen in China in 2012. The objective of this study is to investigate the genetic characteristics and molecular epidemiology of HFMD caused by CV-A6 to guide the diagnosis and treatment of the disease, as well as disease prevention.

MATERIAL AND METHODS

A total of 138 suspected HFMD cases were enrolled in this study and analyses based on complete VP1 nucleotide sequences were performed to determine the evolutionary trajectory of emerging CV-A6.

RESULTS

Among 138 samples in Jiujiang, 125 (90.58%) were positive for enterovirus, the most frequently presented serotypes were CV-A6 (77, 61.60%), CV-A16 (28, 22.40%), EV-A71 (6, 4.80%) and untyped enteroviruses (14, 11.20%). Seventy-seven CV-A6 positive specimens were analyzed for the complete VP1 sequences by sequencing and 36 representative isolates were selected to perform nucleotide sequence similarity analysis. The results showed that 36 strains isolated from HFMD patients were clustered closely to the mainland China and were far from prototype strain CV-A6/Gdula (AY421764) and other international subtypes. Moreover, phylogenetic analysis of the VP1 gene revealed that 36 circulating strains were not significantly concentrated in one branch, but were widely distributed in each branch.

CONCLUSIONS

Continuous surveillance of HFMD etiological agents other than EV-A71 and CV-A16 is necessary. CV-A6 is emerging as the most common pathogen causing HFMD. Closely monitoring the magnitude and trend of CV-A6 epidemic and the trend of pathogenic spectrum changes can provide scientific basis for this disease prevention and control to the department of disease control.

A decade of sustained selection pressure on two surface sites of the VP1 protein of Enterovirus A71 suggests that immune evasion may be an indirect driver for virulence

Enterovirus A71 (EV-A71) is an emerging pathogen in the Enterovirus A species group. EV-A71 causes hand, foot and mouth disease (HFMD), with virulent variants exhibiting polio-like acute flaccid paralysis and other central nervous system manifestations. We analysed all enterovirus A71 complete genomes with collection dates from 2008 to mid-2018. All sub-genotypes exhibit a strong molecular clock with omega (dN/dS) suggesting strong purifying selection. In sub-genotypes B5 and C4, positive selection can be detected at two surface sites on the VP1 protein, also detected in positive selection studies performed prior to 2008. Toggling of a limited repertoire of amino acids at these positively selected residues over the last decade suggests that EV-A71 may be undergoing a sustained frequency-dependent selection process for immune evasion, raising issues for vaccine development. These same sites have also been previously implicated in virus-host binding and strain-associated severity of HFMD, suggesting that immune evasion may be an indirect driver for virulence (154 words).

Emerging Coxsackievirus A6 Causing Hand, Foot and Mouth Disease, Vietnam

Hand, foot and mouth disease (HFMD) is a major public health issue in Asia and has global pandemic potential. Coxsackievirus A6 (CV-A6) was detected in 514/2,230 (23%) of HFMD patients admitted to 3 major hospitals in southern Vietnam during 2011–2015. Of these patients, 93 (18%) had severe HFMD. Phylogenetic analysis of 98 genome sequences revealed they belonged to cluster A and had been circulating in Vietnam for 2 years before emergence. CV-A6 movement among localities within Vietnam occurred frequently, whereas viral movement across international borders appeared rare. Skyline plots identified fluctuations in the relative genetic diversity of CV-A6 corresponding to large CV-A6–associated HFMD outbreaks worldwide. These data show that CV-A6 is an emerging pathogen and emphasize the necessity of active surveillance and understanding the mechanisms that shape the pathogen evolution and emergence, which is essential for development and implementation of intervention strategies.

Surveillance for severe hand, foot, and mouth disease from 2009 to 2015 in Jiangsu province: epidemiology, etiology, and disease burden

Background

Severe hand, foot, and mouth disease (HFMD) is a common childhood illness caused by various enteroviruses. The disease has imposed increased burden on children younger than 5 years old. We aimed to determine the epidemiology, CNS complication, and etiology among severe HFMD patients, in Jiangsu, China.

Methods

Epidemiological, clinical, and laboratory data of severe HFMD cases were extracted from 2009 to 2015. The CNS complication, annually severe illness rates, mortality rates, severity-PICU admission rates, severity-hospitalization rates, and so on were analyzed to assess the disease burden of severe HFMD. All analyses were stratified by time, region, population, CNS involvement and serotypes. The VP1 gene from EV-A71, CV-A16, CV-A6, CV-A10 and other enteroviruses isolates was amplified. Phylogenetic analysis was performed using MEGA5.0.

Results

Seven thousand nine hundred ninety-four severe HFMD cases were reported, of them, 7224 cases were inpatients, 611 were PICU inpatients, and 68 were fatal. The average severe illness rate, mortality rate, severity−fatality rate, severity-PICU admission rate, and severity-hospitalization rate were 14.54, 0.12,8506, 76,430, and 903,700 per 1 million, respectively. The severe illness rate was the highest in the 12–23 months age group, and the greatest mortality rate was in the 6–11 months age group. Geographical difference in severe illness rate and mortality were found. Patients infected with EV-A71 were at a higher proportion in different CNS involvement even death. EV-A71, CV-A16 and other enteroviruses accounted for 79.14, 6.49, and 14.47%, respectively. A total of 14 non-EV-A71/ CV-A16 genotypes including CV-A2, CV-A4, CV-A 6, CV-A9, CV-A10, CV-B1, CV-B2, CV-B3, CV-B4, CV-B5, E-6, E-7, E-18, and EV-C96 were identified. Phylogentic analyses demonstrated that EV-A71 strains belonged to subgenotype C4a, while CV-A16 strains belonged to sub-genotype B1a and sub-genotype B1b of genotype B1. CV-A6 strains were assigned to genogroup F, and CV-A10 strains belonged to genogroup D.

Conclusions

Future mitigation policies should take into account the age, region heterogeneities, CNS conditions and serotype of disease. Additional a more rigorous study between the mild and severe HFMD should be warranted to elucidate the difference epidemiology, pathogen spectrum and immunity patterns and to optimize interventions in the following study.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3659-7) contains supplementary material, which is available to authorized users.

Longitudinal study on enterovirus A71 and coxsackievirus A16 genotype/subgenotype replacements in hand, foot and mouth disease patients in Thailand, 2000-2017

BACKGROUND

Enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) are the major causative agents of hand, foot and mouth disease (HFMD) worldwide, particularly in the Asia-Pacific region. Several strains have emerged, circulated, and faded out over time in recent decades. This study investigated the EV-A71 and CV-A16 circulating strains and replacement of genotypes/subgenotypes in Thailand during the years 2000-2017.

METHODS

The complete VP1 regions of 92 enteroviruses obtained from 90 HFMD patients, one asymptomatic adult contact case, and one encephalitic case were sequenced and investigated for serotypes, genotypes, and subgenotypes using a phylogenetic analysis.

RESULTS

The 92 enterovirus isolates were identified as 67 (72.8%) EV-A71 strains comprising subgenotypes B4, B5, C1, C2, C4a, C4b and C5, and 25 (27.2%) CV-A16 strains comprising subgenotypes B1a and B1b. Genotypic/subgenotypic replacements were evidenced during the study period. EV-A71 B5 and C4a have been the major circulating strains in Thailand for more than a decade, and CV-A16 B1a has been circulating for almost two decades.

CONCLUSIONS

This study provides chronological data on the molecular epidemiology of EV-A71 and CV-A16 subgenotypes in Thailand. Subgenotypic replacement frequently occurred with EV-A71, but not CV-A16. Monitoring for viral genetic and subgenotypic changes is important for molecular diagnosis, vaccine selection, and vaccine development.

Complete genome analysis demonstrates multiple introductions of enterovirus 71 and coxsackievirus A16 recombinant strains into Thailand during the past decade

Hand, foot, and mouth disease (HFMD) caused by enteroviruses remains a public health threat, particularly in the Asia-Pacific region during the past two decades. Moreover, the introduction of multiple subgenotypes and the emergence of recombinant viruses is of epidemiological importance. Based on either the full genome or VP1 sequences, 32 enteroviruses (30 from HFMD patients, 1 from an encephalitic patient, and 1 from an asymptomatic contact case) isolated in Thailand between 2006 and 2014 were identified as 25 enterovirus 71 (EV71) isolates (comprising 20 B5, 1 C2, 2 C4a, and 2 C4b subgenotypes) and 7 coxsackievirus A16 (CA16) isolates (comprising 6 B1a and 1 B1b subgenotypes). The EV71 subgenotype C4b was introduced into Thailand for the first time in 2006 and was replaced by subgenotype C4a strains in 2009. Phylogenetic, similarity plot and bootscan analyses of the complete viral genomes identified 12 recombinant viruses among the 32 viral isolates. Only one EV71-B5 isolate out of 20 was a recombinant virus with one region of intratypic or intertypic recombination, while all four EV71-C4 isolates were recombinant viruses having undergone double recombination, and all seven CA16 isolates were recombinant viruses. The recombination breakpoints of these recombinants are located solely within the P2 and P3 regions. Surveillance for circulating strains and subgenotype replacement are important with respect to molecular epidemiology and the selection of the upcoming EV71 vaccine. In addition, the clinical importance of recombinant viruses needs to be further explored.

Severity of enterovirus A71 infection in a human SCARB2 knock-in mouse model is dependent on infectious strain and route

Enterovirus A71 (EV-A71) is a major etiological agent of human hand, foot and mouth disease, and it can cause severe neurological complications. Although several genotypes of EV-A71 strains are prevalent in different regions of the world, the genotype C4 has circulated in mainland China for more than 20 years. The pathogenicity of different EV-A71 clinical isolates varies and needs to be explored. In this study, hSCARB2 knock-in mice (N = 181) with a wide range of ages were tested for their susceptibility to two EV-A71 strains with the subgenotypes C4 and C2, and two infection routes (intracranial and venous) were compared. The clinical manifestations and pathology and their relationship to the measured viral loads in different tissues were monitored. We observed that 3 weeks is a crucial age, as mice younger than 3-week-old that were infected became extremely ill. However, mice older than 3 weeks displayed diverse clinical symptoms. Significant differences were observed in the pathogenicity of the two strains with respect to clinical signs, disease incidence, survival rate, and body weight change. We concluded that hSCARB2 knock-in mice are a sensitive model for investigating the clinical outcomes resulting from infection by different EV-A71 strains. The intracranial infection model appears to be suitable for evaluating EV-A71 neurovirulence, whereas the venous infection model is appropriate for studying the pathogenicity of EV-A71.

Hand, foot and mouth disease and herpangina caused by enterovirus A71 infections: a review of enterovirus A71 molecular epidemiology, pathogenesis, and current vaccine development

Enterovirus A71 (EV-A71) infections are one of the main etiological agents of hand, foot and mouth disease (HFMD) and herpangina worldwide. EV-A71 infection is a life-threatening communicable disease and there is an urgent global need for the development of vaccines for its prevention and control. The morbidity rate of EV-A71 infection differs between countries. The pathogen’s genetic lineages are undergoing rapid evolutionary changes. An association between the occurrence of EV-A71 infection and the circulation of different genetic strains of EV-A71 virus has been identified around the world. In this review, we present and discuss the molecular epidemiology and pathogenesis of the human disease caused by EV-A71 infection, as well as current prospects for the development of an EV-A71 vaccine.

Molecular epidemiology analysis of enterovirus 71 strains isolated in Dak Lak, Vietnam, 2011-2016

Human enterovirus 71 (EV71) is the major etiologic agent of hand, foot, and mouth disease (HFMD). EV71 outbreaks have been reported in Dak Lak in recent years, however, the genotypes/subgenotypes information and phylogeny of circulating EV71 strains are limited. The objectives of this study were to determine the genotypes/subgenotypes and investigate the phylogeny of EV71 isolates in Dak Lak over a 6-year period. Viruses were isolated from clinical samples from patients with HFMD. In total, 43 EV71 isolates circulated in Dak Lak during 2011-2016 were used for the phylogenetic analysis using complete VP1 gene. The phylogenetic analysis of the VP1 gene revealed that two major genotypes, B and C, were found. Among the 43 EV71 strains, 29 belonged to subgenotype C4, 2 belonged to subgenotype C5, and 12 belonged to subgenotype B5. Of these, the subgenotype C4 was predominant in 2011-2013 and this was later replaced by the subgenotype B5 in 2014. The subgenotype B5 was dominant between 2014 and 2015, and then C4 recirculated in 2016. Our study also indicated that the subgenotypes C4 and B5 emerged into Dak Lak were closely related to variants causing epidemics of HFMD in the southern and central region of Vietnam and Thailand. Sequence analysis showed that nine amino acid mutations were detected in the VP1 region. Our results identified two significant amino acid substitutions (D31N and E145G/Q) associated with enhancing EV71 virulence.