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

The impact of EV71 vaccination program on hand, foot and mouth disease in Zhejiang Province, China: A negative control study

Objective

To estimate the potential causal impact of Enterovirus A71 (EV71) vaccination program on the reduction of EV71-infected hand, foot, and mouth disease (HFMD) in Zhejiang Province.

Methods

We utilized the longitudinal surveillance dataset of HFMD and EV71 vaccination in Zhejiang Province during 2010-2019. We estimated vaccine efficacy using a Bayesian structured time series (BSTS) model, and employed a negative control outcome (NCO) model to detect unmeasured confounding and reveal potential causal association.

Results

We estimated that 20,132 EV71 cases (95% CI: 16,733, 23,532) were prevented by vaccination program during 2017-2019, corresponding to a reduction of 29% (95% CI: 24%, 34%). The effectiveness of vaccination increased annually, with reductions of 11% (95% CI: 6%, 16%) in 2017 and 66% (95% CI: 61%, 71%) in 2019. Children under 5 years old obtained greater benefits compared to those over 5 years. Cities with higher vaccination coverage experienced a sharper EV71 reduction compared to those with lower coverage. The NCO model detected no confounding factors in the association between vaccination and EV71 cases reduction.

Conclusions

This study suggested a potential causal effect of the EV71 vaccination, highlighting the importance of achieving higher vaccine coverage to control the HFMD.

An enterovirus A71 virus-like particle with replaced loops confers partial cross-protection in mice

Enterovirus A71 (EV-A71), coxsackievirus A16 (CV-A16), and CV-A10 belong to the main prevailing types causing hand-foot-and-mouth disease. Since EV-A71 monovalent vaccine does not confer cross-protection, developing a multivalent vaccine is essential. In this study, a trivalent chimeric virus-like particle of EV-A71 (EV-A71-VLPCHI3) was constructed based on EV-A71-VLP backbone by replacing the corresponding surface loops with CV-A16 VP1 G-H, CV-A10 VP1 B-C and E-F loops, which are critical for immunogenic neutralization. The baculovirus-insect cell expression system was employed for EV-A71-VLPCHI3 production. EV-A71-VLPCHI3 was purified by sucrose density gradient and observed by transmission electron microscopy. The immunogenicity and protective efficacy of EV-A71-VLPCHI3 were evaluated in mice. Our results revealed that EV-A71-VLPCHI3 had a similar morphology to inactivated EV-A71 particles and could induce specific IgG antibodies against EV-A71, CV-A16 and CV-A10 in mice. More importantly, EV-A71-VLPCHI3 enhanced cross-reactive protection against CV-A16 and CV-A10, by 20 % and 40 %, compared to inactivated EV-A71 counterparts, respectively. In conclusion, the successful construction of EV-A71-VLPCHI3 suggested that loop-dependent heterologous protection could be transferred by loops replacement on the surface of viral capsid. This strategy may also supplement the development of multivalent vaccines against other infectious viral diseases.

The neuropathological mechanism of EV-A71 infection attributes to inflammatory pryoptosis and viral replication via activating the hsa_circ_0045431/ hsa_miR_584/NLRP3 regulatory axis

Neuropathological damage has been considered to be the main cause of death from EV-A71 infection, but the underlying mechanism has not been elucidated. Pyroptosis, a new form of inflammatory programmed cell death, has been verified to be involved in the pathogenesis of various viruses. circRNAs are a novel type of endogenous noncoding RNA gaining research interest in recent years, especially their special roles in the process of virus infection. Thus, in this study, we combined EV-A71, pyroptosis and circRNA to find a breakthrough in the pathogenesis of EV-A71 infection. Firstly, whether EV-A71 infection leaded to pyroptosis formation was examined by a series detection of cell death, cell viability, LDH release, caspase 1 activity, the expression levels of pyroptosis-related molecules and the concentrations of IL-1β and IL-18. Secondly, high-throughput sequencing of circRNAs was carried out to excavate the circRNA-miRNA-mRNA regulatory axis which might be associated with pyroptosis formation. Finally, the gain- and loss-of-functional experiments were further conducted to identify their functions. Our results showed that EV-A71 infection caused pyroptosis formation in SH-SY5Y cells. The circRNA sequencing analyzed the differentially expressed circRNAs and their possible functions. It was found that the hsa_circ_0045431/hsa_miR_584/NLRP3 regulatory axis might be involved in pyroptosis formation during EV-A71 infection. Then, hsa_circ_0045431 sponged hsa_miR_584 and hsa_miR_584 directly targeted NLRP3 were validated by IF, dual-luciferase, qRT-PCR and WB assays. Functional experiments were performed to further uncover that the up-regulation of hsa_circ_0045431 and NLRP3 promoted the inflammatory pyroptosis and viral replication, while the up-regulation of hsa_miR_584 suppressed the inflammatory pyroptosis and viral replication, and vice versa. Collectively, our study demystified that EV-A71 infection induced pyroptosis formation by activating hsa_circ_0045431/hsa_miR_584/NLRP3 regulatory axis, which could further effect viral replication. These findings provided novel insights into the pathogenesis of EV-A71 infection, and meanwhile revealed that the hsa_circ_0045431/ hsa_miR_584/NLRP3 regulatory axis can serve as a potential biological therapeutic target for EV-A71 infection.

HMGB1 is involved in viral replication and the inflammatory response in coxsackievirus A16-infected 16HBE cells via proteomic analysis and identification

Coxsackievirus A16 (CV-A16) is still an important pathogen that causes hand, foot and mouth disease (HFMD) in young children and infants worldwide. Previous studies indicated that CV-A16 infection is usually mild or self-limiting, but it was also found that CV-A16 infection can trigger severe neurological complications and even death. However, there are currently no vaccines or antiviral compounds available to either prevent or treat CV-A16 infection. Therefore, investigation of the virus‒host interaction and identification of host proteins that play a crucial regulatory role in the pathogenesis of CV-A16 infection may provide a novel strategy to develop antiviral drugs. Here, to increase our understanding of the interaction of CV-A16 with the host cell, we analyzed changes in the proteome of 16HBE cells in response to CV-A16 using tandem mass tag (TMT) in combination with LC‒MS/MS. There were 6615 proteins quantified, and 172 proteins showed a significant alteration during CV-A16 infection. These differentially regulated proteins were involved in fundamental biological processes and signaling pathways, including metabolic processes, cytokine‒cytokine receptor interactions, B-cell receptor signaling pathways, and neuroactive ligand‒receptor interactions. Further bioinformatics analysis revealed the characteristics of the protein domains and subcellular localization of these differentially expressed proteins. Then, to validate the proteomics data, 3 randomly selected proteins exhibited consistent changes in protein expression with the TMT results using Western blotting and immunofluorescence methods. Finally, among these differentially regulated proteins, we primarily focused on HMGB1 based on its potential effects on viral replication and virus infection-induced inflammatory responses. It was demonstrated that overexpression of HMGB1 could decrease viral replication and upregulate the release of inflammatory cytokines, but deletion of HMGB1 increased viral replication and downregulated the release of inflammatory cytokines. In conclusion, the results from this study have helped further elucidate the potential molecular pathogenesis of CV-A16 based on numerous protein changes and the functions of HMGB1 Found to be involved in the processes of viral replication and inflammatory response, which may facilitate the development of new antiviral therapies as well as innovative diagnostic methods.

Tandem mass tag (TMT) labeling-based quantitative proteomic analysis reveals the cellular protein characteristics of 16HBE cells infected with coxsackievirus A10 and the potential effect of HMGB1 on viral replication

Coxsackievirus A10 (CV-A10) is recognized as one of the most important pathogens associated with hand, foot, and mouth disease (HFMD) in young children under 5 years of age worldwide, and it can lead to fatal neurological complications. However, available commercial vaccines fail to protect against CV-A10. Therefore, there is an urgent need to study new protein targets of CV-A10 and develop novel vaccine-based therapeutic strategies. Advances in proteomics in recent years have enabled a comprehensive understanding of host pathogen interactions. Here, to study CV-A10-host interactions, a global quantitative proteomic analysis was conducted to investigate the molecular characteristics of host cell proteins and identify key host proteins involved in CV-A10 infection. Using tandem mass tagging (TMT)-based mass spectrometry, a total of 6615 host proteins were quantified, with 293 proteins being differentially regulated. To ensure the validity and reliability of the proteomics data, three randomly selected proteins were verified by Western blot analysis, and the results were consistent with the TMT results. Further functional analysis showed that the upregulated and downregulated proteins were associated with diverse biological activities and signaling pathways, such as metabolic processes, biosynthetic processes, the AMPK signaling pathway, the neurotrophin signaling pathway, the MAPK signaling pathway, and the GABAergic synaptic signaling. Moreover, subsequent bioinformatics analysis demonstrated that these differentially expressed proteins contained distinct domains, were localized in different subcellular components, and generated a complex network. Finally, high-mobility group box 1 (HMGB1) might be a key host factor involved in CV-A10 replication. In summary, our findings provide comprehensive insights into the proteomic profile during CV-A10 infection, deepen our understanding of the relationship between CV-A10 and host cells, and establish a proteomic signature for this viral infection. Moreover, the observed effect of HMGB1 on CV-A10 replication suggests that it might be a potential therapeutic target treatment of CV-A10 infection.

MicroRNA expression profile of human umbilical vein endothelial cells in response to coxsackievirus A10 infection reveals a potential role of miR-143-3p in maintaining the integrity of the blood-brain barrier

Coxsackievirus A10 (CV-A10) has been one of the main etiologies of hand, foot, and mouth disease (HFMD) epidemics in recent years and can cause mild to severe illness and even death. Most of these severe and fatal cases were closely associated with neurological impairments, but the potential mechanism of neuropathological injury triggered by CV-A10 infection has not been elucidated. MicroRNAs (miRNAs), implicated in the regulation of gene expression in a post-transcriptional manner, play a vital role in the pathogenesis of various central nervous system (CNS) diseases; therefore, they serve as diagnostic biomarkers and are emerging as novel therapeutic targets for CNS injuries. To gain insights into the CV-A10-induced regulation of host miRNA-processing machinery, we employed high-throughput sequencing to identify differentially expressed miRNAs in CV-A10-infected human umbilical vein endothelial cells (HUVECs) and further analyzed the potential functions of these miRNAs during CV-A10 infection. The results showed that CV-A10 infection could induce 189 and 302 significantly differentially expressed miRNAs in HUVECs at 24 and 72 hpi, respectively, compared with the uninfected control. Moreover, the expression of four selected miRNAs and their relevant mRNAs was determined to verify the sequencing data by quantitative reverse transcription-polymerase chain reaction (RT-qPCR) methods. After that, gene target prediction and functional annotation revealed that the targets of these dysregulated miRNAs were mostly enriched in cell proliferation, signal transduction, cAMP signalling pathway, cellular response to interleukin-6, ventral spinal cord interneuron differentiation, negative regulation of glial cell differentiation, neuron migration, positive regulation of neuron projection development, etc., which were primarily involved in the processes of basic physiology, host immunity, and neurological impairments and further reflected vital regulatory roles of miRNA in viral pathogenicity. Finally, the construction of a miRNA-regulated network also suggested that the complex regulatory mechanisms mediated by miRNAs might be involved in viral pathogenesis and virus-host interactions during CV-A10 infection. Furthermore, among these dysregulated miRNAs, miR-143-3p was demonstrated to be involved in the maintenance of blood-brain barrier (BBB) integrity.

Effect of proteins isolated from Brazilian snakes on enterovirus A71 replication cycle: An approach against hand, foot and mouth disease

Enterovirus A71 (EVA71) belongs to the Picornaviridae family and is the main etiological agent of hand, foot, and mouth disease (HFMD). There is no approved antiviral against EVA71, and therefore the search for novel anti-EVA71 therapeutics is essential. In this context, the antiviral activity of proteins isolated from snake venoms has been reported against a range of viruses. Here, the proteins CM10 and CM14 isolated from Bothrops moojeni, and Crotamin and PLA2CB isolated from Crotalus durissus terrificus were investigated for their antiviral activity against EVA71 infection. CM14 and Crotamin possessed a selective index (SI) of 170.8 and 120.4, respectively, while CM10 and PLA2CB had an SI of 67.4 and 12.5, respectively. CM14 inhibited all steps of viral replication (protective effect: 76 %; virucidal: 99 %; and post-entry: 99 %). Similarly, Crotamin inhibited up to 99 % of three steps. In contrast, CM10 and PLA2CB impaired one or two steps of EVA71 replication, respectively. Further dose-response assays using increasing titres of EVA71 were performed and CM14 and Crotamin retained functionality with high concentrations of EVA71 (up to 1000 TCID50). These data demonstrate that proteins isolated from snake venom are potent inhibitors of EVA71 and could be used as scaffolds for future development of novel antivirals.

Global quantitative proteomic analysis profiles of host protein expression in response to Enterovirus A71 infection in bronchial epithelial cells based on tandem mass tag (TMT) peptide labeling coupled with LC-MS/MS uncovers the key role of proteasome in virus replication

Enterovirus A71 (EV-A71) is a neurotropic human pathogen which mainly caused hand, foot and mouth disease (HFMD) mostly in children under 5 years-old. Generally, EV-A71-associated HFMD is a relatively self-limiting febrile disease, but there will still be a small percentage of patients with rapid disease progression and severe neurological complications. To date, the underlying mechanism of EV-A71 inducing pathological injury of central nervous system (CNS) remains largely unclear. It has been investigated and discussed the changes of mRNA, miRNA and circRNA expression profile during infection by EV-A71 in our previous studies. However, these studies were only analyzed at the RNA level, not at the protein level. It's the protein levels that ultimately do the work in the body. Here, to address this, we performed a tandem mass tag (TMT) peptide labeling coupled with LC-MS/MS approach to quantitatively identify cellular proteome changes at 24 h post-infection (hpi) in EV-A71-infected 16HBE cells. In total, 6615 proteins were identified by using TMT coupled with LC-MS/MS in this study. In the EV-A71- and mock-infected groups, 210 differentially expressed proteins were found, including 86 upregulated and 124 downregulated proteins, at 24 hpi. To ensure the validity and reliability of the proteomics data, 3 randomly selected proteins were verified by Western blot and Immunofluorescence analysis, and the results were consistent with the TMT results. Subsequently, functional enrichment analysis indicated that the up-regulated and down-regulated proteins were individually involved in various biological processes and signaling pathways, including metabolic process, AMPK signaling pathway, Neurotrophin signaling pathway, Viral myocarditis, GABAergic synapse, and so on. Moreover, among these enriched functional analysis, the "Proteasome" pathway was up-regulated, which has caught our attention. Inhibition of proteasome was found to obviously suppress the EV-A71 replication. Finally, further in-depth analysis revealed that these differentially expressed proteins contained distinct domains and localized in different subcellular components. Taken together, our data provided a comprehensive view of host cell response to EV-A71 and identified host proteins may lead to better understanding of the pathogenic mechanisms and host responses to EV-A71 infection, and also to the identification of new therapeutic targets for EV-A71 infection.

Advances and Breakthroughs in IRES-Directed Translation and Replication of Picornaviruses

Viruses lack the properties to replicate independently due to the limited resources encoded in their genome; therefore, they hijack the host cell machinery to replicate and survive. Picornaviruses get the prerequisite for effective protein synthesis through specific sequences known as internal ribosome entry sites (IRESs). In the past 2 decades, significant progress has been made in identifying different types of IRESs in picornaviruses. This review will discuss the past and current findings related to the five different types of IRESs and various internal ribosome entry site trans-acting factors (ITAFs) that either promote or suppress picornavirus translation and replication. Some IRESs are inefficient and thus require ITAFs. To achieve their full efficiency, they recruit various ITAFs, which enable them to translate more effectively and efficiently, except type IV IRES, which does not require any ITAFs. Although there are two kinds of ITAFs, one promotes viral IRES-dependent translation, and the second type restricts. Picornaviruses IRESs are classified into five types based on their use of sequence, ITAFs, and initiation factors. Some ITAFs regulate IRES activity by localizing to the viral replication factories in the cytoplasm. Also, some drugs, chemicals, and herbal extracts also regulate viral IRES-dependent translation and replication. Altogether, this review will elaborate on our understanding of the past and recent advancements in the IRES-dependent translation and replication of picornaviruses. IMPORTANCE The family Picornaviridae is divided into 68 genera and 158 species. The viruses belonging to this family range from public health importance, such as poliovirus, enterovirus A71, and hepatitis A virus, to animal viruses of great economic importance, such as foot-and-mouth disease virus. The genomes of picornaviruses contain 5' untranslated regions (5' UTRs), which possess crucial and highly structured stem-loops known as IRESs. IRES assemble the ribosomes and facilitate the cap-independent translation. Virus-host interaction is a hot spot for researchers, which warrants deep insight into understanding viral pathogenesis better and discovering new tools and ways for viral restriction to improve human and animal health. The cap-independent translation in the majority of picornaviruses is modulated by ITAFs, which bind to various IRES regions to initiate the translation. The discoveries of ITAFs substantially contributed to understanding viral replication behavior and enhanced our knowledge about virus-host interaction more effectively than ever before. This review discussed the various types of IRESs found in Picornaviridae, past and present discoveries regarding ITAFs, and their mechanism of action. The herbal extracts, drugs, and chemicals, which indicated their importance in controlling viruses, were also summarized. In addition, we discussed the movement of ITAFs from the nucleus to viral replication factories. We believe this review will stimulate researchers to search for more novel ITAFs, drugs, herbal extracts, and chemicals, enhancing the understanding of virus-host interaction.

Etiological Analysis of Viral Encephalitis in Children in Zhejiang Province from 2018 to 2019

Objective: To investigate the common pathogens of viral encephalitis (VE) in children, and to provide guidance for the empirical diagnosis and treatment of patients with VE. Methods: A total of 227 cerebrospinal fluid (CSF) samples were collected from pediatric patients with VE in Zhejiang province from January 2018 to December 2019. The samples were tested using multiplex and singleplex Reverse Transcription-Polymerase Chain Reaction (RT-PCR) with primers specific to enterovirus (EV), varicella-zoster virus (VZV), mumps virus (MuV), cytomegalovirus (CMV), herpes simplex virus type 1 (HSV-1)/type 2 (HSV-2), Epstein–Barr virus (EBV), and human herpesvirus 6 (HHV-6). The data of the two analyses were compared and then verified using Sanger sequencing. Results: Of the 227 CSF samples, 90 were shown to be positive for multiplex RT-PCR with a positivity rate of 39.65% and a 95% confidence interval (33.2%, 46.1%). EV was the most common cause of VE, followed by EBV, HHV-6, MuV, CMV, VZV, and HSV-1. Most included cases occurred in summer, accounting for 49.78% of all cases. For EV, EBV, and HSV-2, multiplex RT-PCR showed a positivity rate of 34.36%, which was not statistically different from that of 30.4% shown by singleplex RT-PCR. The sequences of EV, EBV, VZV, MuV, CMV, HSV-1, HHV-6, and HSV-2 were confirmed by sequencing the PCR products obtained from multiplex and singleplex PCR. Conclusions: In children, VE is more prevalent in the summer than in other seasons in Zhejiang province, and EV may be the most common causative pathogen.

Prediction of hand, foot, and mouth disease epidemics in Japan using a long short-term memory approach

Hand, foot, and mouth disease (HFMD) is a common febrile illness caused by enteroviruses in the Picornaviridae family. The major symptoms of HFMD are fever and a vesicular rash on the hand, foot, or oral mucosa. Acute meningitis and encephalitis are observed in rare cases. HFMD epidemics occur annually in Japan, usually in the summer season. Relatively large-scale outbreaks have occurred every two years since 2011. In this study, the epidemic patterns of HFMD in Japan are predicted four weeks in advance using a deep learning method. The time-series data were analyzed by a long short-term memory (LSTM) approach called a Recurrent Neural Network. The LSTM model was trained on the numbers of weekly HFMD cases in each prefecture. These data are reported in the Infectious Diseases Weekly Report, which compiles the national surveillance data from web sites at the National Institute of Infectious Diseases, Japan, under the Infectious Diseases Control Law. Consequently, our trained LSTM model distinguishes between relatively large-scale and small-scale epidemics. The trained model predicted the HFMD epidemics in 2018 and 2019, indicating that the LSTM approach can estimate the future epidemic patterns of HFMD in Japan.

Comprehensive profiling and characterization of cellular microRNAs in response to coxsackievirus A10 infection in bronchial epithelial cells

Coxsackievirus A10 (CV-A10), the causative agent of hand, foot, and mouth disease (HFMD), caused a series of outbreaks in recent years and often leads to neurological impairment, but a clear understanding of the disease pathogenesis and host response remains elusive. Cellular microRNAs (miRNAs), a large family of non-coding RNA molecules, have been reported to be key regulators in viral pathogenesis and virus-host interactions. However, the role of host cellular miRNAs defensing against CV-A10 infection is still obscure. To address this issue, we systematically analyzed miRNA expression profiles in CV-A10-infected 16HBE cells by high-throughput sequencing methods in this study. It allowed us to successfully identify 312 and 278 miRNAs with differential expression at 12 h and 24 h post-CV-A10 infection, respectively. Among these, 4 miRNAs and their target genes were analyzed by RT-qPCR, which confirmed the sequencing data. Gene target prediction and enrichment analysis revealed that the predicted targets of these miRNAs were significantly enriched in numerous cellular processes, especially in regulation of basic physical process, host immune response and neurological impairment. And the integrated network was built to further indicate the regulatory roles of miRNAs in host-CV-A10 interactions. Consequently, our findings could provide a beneficial basis for further studies on the regulatory roles of miRNAs relevant to the host immune responses and neuropathogenesis caused by CV-A10 infection.

Immunogenicity and safety of the inactivated enterovirus 71 vaccine administered concomitantly with the measles-rubella vaccine in infants aged 8 months in China: A noninferiority randomized controlled trial

BACKGROUND

To evaluate the immunogenicity and safety of simultaneous administration of the enterovirus 71 (EV71) vaccine with the measles and rubella (MR) combined vaccine.

METHODS

In this phase 4, randomized, open-label and noninferiority study, a total of 680 infants aged 8 months were enrolled and assigned to the simultaneous administration group (infants received the first dose of EV71 vaccine and MR vaccine on Day 0, and the second dose of EV71 vaccine on Day 28), or the separate administration groups (EV71 group: infants received two doses of EV71 vaccine on Day 0 and Day 28, respectively; MR group: infants received MR vaccine on Day 0). Blood sample was obtained on Day 0 and Day 56 to measure antibody responses to each of the antigens in terms of antibody titer or concentration, respectively. Local and systemic adverse reactions (ARs) and other adverse events (AEs) following each dose were monitored and compared among groups.

RESULTS

After vaccination, simultaneous administration group showed similar seroconversion rates of antibody against EV71(97.9%), measles (97.4%), and rubella (94.3%) compared to EV71 group (99.6% for anti-EV71) or MR group (98.4% for anti-measles and 98.9% for anti-rubella, respectively). Noninferiority was demonstrated for all antibodies as the lower limits of two-sided 97.5% confidence intervals (CIs) of the difference in seroconversion rates between simultaneous administration group and separate administration groups were above the predefined margin of -10%. Additionally, the adverse reaction rates were comparable among groups (54.4% in the simultaneous group versus 43.9% in the MR group versus 52.6% in the EV71 group).

CONCLUSION

Antibody responses induced by simultaneous administration of EV71 vaccine with MR vaccine were robust and noninferior to those by single administration alone. Like the previous findings by single administration alone, simultaneous administration demonstrated comparable reactogenicity and safety profiles.

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.

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.

circRNA expression patterns and circRNA-miRNA-mRNA networks during CV-A16 infection of SH-SY5Y cells

Coxsackievirus A16 (CV-A16) has caused worldwide epidemics of hand, foot, and mouth disease (HFMD) in infants and preschool children. Circular RNAs (circRNAs), a class of noncoding RNA molecules, participate in the progression of viral infectious diseases. Although the function of circRNAs has been a heavily researched topic, their role in CV-A16 infection is still unclear. In this study, the viral effects of CV-A16 on the cellular circRNA transcriptome were investigated using next-generation sequencing technology. The results showed that a total of 8726, 8611, and 6826 circRNAs were identified at 0, 12, and 24 h postinfection, respectively. Moreover, it was found that 1769 and 1192 circRNAs were differentially expressed in at 12 and 24 h postinfection, respectively. The common differentially expressed circRNAs were used for functional annotation analysis, and it was found that the parent genes of differentially expressed circRNAs might be associated with the viral infection process, especially the “Immune system process” in GO analysis and the “Inflammation mediated by chemokine and cytokine signaling pathway” in KEGG analysis. Subsequently, circRNA-miRNA-mRNA regulatory networks were constructed, and the hsa_circ_0004447/hsa-miR-942-5p/MMP2, hsa_circ_0078617/hsa-miR-6780b-5p/MMP2 and hsa_circ_0078617/hsa-miR-5196-5p/MMP2 regulatory axes were identified by enrichment analysis as important networks during the progression of CV-A16 infection. Finally, six dysregulated circRNAs were selected for validation and were verified to be consistent with the sequencing results. Considering all of these results, to the best of our knowledge, this study is the first to present a comprehensive overview of circRNAs induced by CV-A16 infection, and this research demonstrated that a network of enriched circRNAs and circRNA-associated competitive endogenous RNAs (ceRNAs) is involved in the regulation of CV-A16 infection, thereby helping to elucidate the mechanisms underlying CV-A16-host interactions.

Comprehensive analysis of the circRNA expression profile and circRNA-miRNA-mRNA network in the pathogenesis of EV-A71 infection

Enterovirus A71 (EV-A71) is an important emerging virus posing a threat to children under five years old. Circular RNAs (circRNAs), a novel class of endogenous RNAs, have been recognized to play important roles in the onset and development of viral diseases. However, it has not been determined which specific circRNAs are involved in the pathological mechanisms of EV-A71 infection. In this study, RNA-sequencing (RNA-seq) was conducted to characterize differentially expressed circRNAs during the process of EV-A71 infection. Overall, 8726, 10405 and 4710 circRNAs were detected in the control, EV-A71-12 h and EV-A71-24 h groups, respectively, of which 1851 and 951 circRNAs were differentially expressed in the EV-A71-12 h and EV-A71-24 h groups versus the control group. The overlapping circRNAs in the EV-A71-infected groups were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, which further demonstrated that the host genes of these circRNAs were principally implicated in activities associated with the progression of viral infection, such as immune system process, Wnt signaling pathway, etc. Additionally, qRT-PCR detection showed that six selected circRNAs were identical to the sequencing data. To excavate the key circRNAs in EV-A71 infection, we comprehensively evaluated and integrated the relationship of circRNA/miRNA/mRNA, and eventually screened 2 key circRNA regulatory axes, namely hsa_circ_0017115/hsa-miR-150-5p/EGR1 and hsa_circ_0005060/hsa-miR-4685-5p/MMP2. In summary, our findings not only provide the first comprehensive expression and functional profile of circRNAs in response to EV-A71 infection, but also offer a novel direction to elucidate the molecular mechanism underlying viral pathogenesis and the cellular immune response in host-EV-A71 interactions.

Enterovirus A71 causing meningoencephalitis and acute flaccid myelitis in a patient receiving rituximab

We present the case of a young woman being treated with rituximab for rheumatoid arthritis who developed a severe enteroviral meningoencephalitis and acute flaccid myelitis (AFM). Cerebrospinal fluid (CSF) and stool reverse transcription-polymerase chain reaction (RT-PCR) testing confirmed the diagnosis and additional sequencing studies performed at the CDC further characterized the enterovirus as enterovirus A71 (EV-A71). After treatment with intravenous immunoglobulin (IVIg) and fluoxetine (based on previous reports of possible efficacy) the patient experienced a remarkable improvement over time. This case highlights the importance of considering enteroviral infection in patients treated with rituximab, depicts a possible clinical course of enteroviral meningoencephalitis and AFM, and illustrates the importance of testing multiple sites for enterovirus infection (CSF, stool, nasopharyngeal swab, blood). Here we present the case with a brief review of the literature pertaining to EV-A71.

Enterovirus A71 Vaccines

Enterovirus A71 (EV-A71) is a major causative agent of hand, foot, and mouth disease (HFMD) and herpangina. Moreover, EV-A71 infection can lead to neurological complications and death. Vaccination is the most efficient way to control virus infection. There are currently three inactivated, whole EV-A71 vaccines licensed by the China NMPA (National Medical Products Administration). Several other types of vaccines, such as virus-like particles and recombinant VP1 (capsid protein), are also under development. In this review, we discuss recent advances in the development of EV-A71 vaccines.

The epidemiological risk factors of hand, foot, mouth disease among children in Singapore: A retrospective case-control study

The incidence of hand, foot, and mouth disease (HFMD) is increasing over the years despite current prevention and control policies in Singapore. A retrospective case-control study was conducted among parents whose children attended childcare centres in Singapore to assess the epidemiological risk factors associated with HFMD among children below 7 years old. Parents of 363 children with HFMD (as cases) and 362 children without HFMD (as controls) were enrolled from 22 childcare centres. Data of potential risk factors were collected through a standardised self-administered questionnaire from parents which include demographics and hygiene practices. Multivariate analysis were adjusted for age group, parent’s education level, mother's age, HFMD-infected siblings, and preschool admission period. Child’s age between 1.5 and 4.9 years, child who had been in childcare for more than 1.9years, having HFMD-infected siblings, two or more children in a family, higher educated parents, parents who had HFMD episode previously, wash toys with soap once every two to three weeks, sanitise toys once every two to three weeks, out-sourced cleaner in childcare centre, no domestic helper at home and more than 22 children in a classroom were independent risk factors of HFMD. These evidence provide crucial implications to guide more effective prevention and control of HFMD in Singapore.

Global profiling of the alternative splicing landscape reveals transcriptomic diversity during the early phase of enterovirus 71 infection

The alteration of host cell splicing is a major strategy favouring viral replication; however, the interaction between human tonsillar epithelial cells (HTECs) and enterovirus 71 (EV71) has not been fully elucidated. Here, a total of 201 differentially expressed genes (DEGs) and 3266 novel genes with coding potential were identified. A total of 3479 skipped exons (SEs), 515 alternative 3' splice sites (A3SSs), 391 alternative 5' splice sites (A5SSs), 531 mutually exclusive exons (MXEs) and 825 retained introns (RIs) were identified as significantly altered alternative splicing (AS) events. The enriched DEGs were mainly related to the cell cycle, spliceosome, and Toll-like receptor (TLR) signalling pathways. Finally, the replication of EV71 was significantly inhibited by TLR2 heterodimers. Our findings suggest that AS events induced by EV71 increase the transcriptomic diversity of HTECs in response to EV71 infection. Additionally, TLR2 heterodimers have the potential to protect HTECs against EV71.

Co-circulation of coxsackieviruses A-6, A-10, and A-16 causes hand, foot, and mouth disease in Guangzhou city, China

Background

Hand, foot, and mouth disease (HFMD) is a common infectious disease occurring in children under 5 years of age worldwide, and Enterovirus A71 (EV-A71) and Coxsackievirus A16 (CVA-16) are identified as the predominant pathogens. In recent years, Coxsackievirus A6 (CVA-6) and Coxsackievirus A10 (CVA-10) have played more and more important role in a series of HFMD outbreaks. This study aimed to understand the epidemic characteristics associated with HFMD outbreak in Guangzhou, 2018.

Methods

The clinical and laboratory data of 1220 enterovirus-associated HFMD patients in 2018 were analysed in this study. Molecular diagnostic methods were performed to identify its serotypes. Phylogenetic analyses were depicted based on the complete VP1 gene.

Results

There were 21 enterovirus serotypes detected in Guangzhou in 2018. Three serotypes of enterovirus, CVA-6 (364/1220, 29.8%), CVA-10 (305/1220, 25.0%), and CVA-16 (397/1220, 32.5%), were identified as the causative pathogens and accounted for 87.3% among all 1220 HFMD patients. In different seasons, CVA-6 was the predominant pathogen of HFMD during autumn, and CVA-10 as well as CVA-16 were more prevalent in summer. Patients infected by CVA-6, CVA-10 or CVA-16 showed similar clinical features and laboratory characteristics, and the ratios of severe HFMD were 5.8, 5.9, and 1.5% in the three serotypes. Phylogenetic analyses of VP1 sequences showed that the CVA-6, CVA-10, and CVA-16 sequences belonged to the sub-genogroup E2, genogroup E, and genogroup B1, respectively.

Conclusions

CVA-6, CVA-10, and CVA-16 were the predominant and co-circulated serotypes in Guangzhou China, 2018, which should be the new target for prevention and control of HFMD. Our findings provide useful information for diagnosis, treatment, and prevention of HFMD.

Transcriptome sequencing analysis of SH-SY5Y cells infected with EV71 reveals the potential neuropathic mechanisms

Enterovirus A71 (EV71) remains the most common causative agent of hand, foot, and mouth disease (HFMD), and the neurological complications induced by EV71 are usually the leading cause of death in children with HFMD. However, the mechanism of nervous system changes caused by EV71 infection is still unclear. Therefore, in the current study, EV71 was inoculated into the human neuroblastoma cell line SH-SY5Y and subsequent transcriptome sequencing was used to examine the alterations of the transcriptome in infected SH-SY5Y cells. It is expected to determine the underlying mechanism of neurological diseases in response to EV71 infection. As a result, a total of 82,406,974, 112,410,808 and 87,780,371 clean reads were found in the control, EV71-12 h and EV71-24 h groups, respectively. Moreover, 160 and 745 differentially expressed genes were identified in the EV71-12 h and EV71-24 h groups, respectively, as compared to the control group. Next, to further explore the pathogenic mechanism triggered by EV71 infection, we mainly focused on the common differentially expressed genes at different time points of EV71 infection. And it was discovered that there were 95 common differentially expressed genes, which were used to conduct GO and pathway analysis. GO enrichment analysis demarcated related biological processes, molecular functions and cellular components, and KEGG pathway analysis enabled annotations of metabolic pathways and revealed interactions among the significantly enriched pathways. The results showed that the enriched GO term "Nervous system development" and enriched pathway "CCKR signaling map" might be important contributors to EV71-induced neuropathological mechanisms. In addition, we also screened 10 up- and down-regulated non-protein coding genes with significantly different expression in our transcriptome profiling, which suggested that these abnormally regulated non-protein-encoding genes might also play important roles in the pathogenesis of EV71 infection. Eventually, RT-qPCR technology was adopted to validate the transcriptome sequencing data and the experiment demonstrated that the RT-qPCR and transcriptome sequencing results were basically consistent. In summary, this is the first transcriptome analysis of SH-SY5Y cells in response to EV71 infection and provides valuable cues for further exploring the mechanism of nervous system changes caused by EV71 infection.

A review and meta-analysis of the epidemiology and clinical presentation of coxsackievirus A6 causing hand-foot-mouth disease in China and global implications

Coxsackievirus A6 (CV-A6) has been associated with increasingly occurred sporadic hand-foot-mouth disease (HFMD) cases and outbreak events in many countries. In order to understand epidemiological characteristics of CV-A6, we collected the information describing HFMD caused by CV-A6 to describe the detection rate, severe rate and onychomadesis rate, which is defined as one or more nails defluvium, caused by CV-A6 from 2007 to 2017. The results showed that there was an outbreak of CV-A6 every other year, and overall trend of the epidemic of CA6-associated HFMD was increasing in China. The detection rate of CV-A6 in other countries was 32.0% (95% CI: 25.0%~40.0%) before 2013 and 28.0% (95% CI: 20.0%~36.0%) after 2013, respectively. Although the severe rate of HFMD caused by CV-A6 was low (0.10%, 95% CI: 0.01%~0.20%), CV-A6 can cause a high incidence of onychomadesis (28.0%, 95%CI: 21.9%-34.3%). Thus, it would be worthwhile to research and develop an effective multivalent vaccine for CV-A6 to achieve a more powerful prevention of HMFD.

Interplays between Enterovirus A71 and the innate immune system

Enterovirus A71 (EV-A71) is a growing threat to public health, particularly in the Asia-Pacific region. EV-A71 infection is most prevalent in infants and children and causes a wide spectrum of clinical complications, including hand-foot-and-mouth disease (HFMD), pulmonary and neurological disorders. The pathogenesis of EV-A71 infection is poorly understood at present. It is likely that viral factors and host immunity, and their interplay, affect the pathogenesis and outcome of EV-A71 infection. The mammalian innate immune system forms the first layer of defense against viral infections and triggers activation of adaptive immunity leading to full protection. In this review, we discuss recent advances in our understanding of the interaction between EV-A71 and the innate immune system. We discuss the role of pattern-recognition receptors (PRRs), including Toll-like receptors (TLRs), RIG-I-like receptors (RLRs), and inflammasomes, in the detection of EV-A71 infection and induction of antiviral immunity. As a counteraction, EV-A71 viral proteins target multiple innate immune pathways to facilitate viral replication in host cells. These novel insights at the virus-host interphase may support the future development of vaccines and therapeutics against EV-A71 infection.