Pathological Characteristics of Echovirus 30 Infection in a Mouse Model.

Echovirus 30 (E30), a member of species B enterovirus, is associated with outbreaks of aseptic meningitis and has become a global health emergency. However, the pathogenesis of E30 remains poorly understood due to the lack of appropriate animal models. In this study, we established a mouse infection model to explore the pathogenicity of E30. The 2-day-old IFNAR-/- mice infected with E30 strain WZ16 showed lethargy and paralysis, and some died. Obvious pathological changes were observed in the skeletal muscle, brain tissue, and other tissues, with the highest viral load in the skeletal muscles. Transcriptome analysis of brain and skeletal muscle tissues from infected mice showed that significant differentially expressed genes were enriched in complement response and neuropathy-related pathways. Using immunofluorescence assay, we found that the viral double-stranded RNA (dsRNA) was detected in the mouse brain region and could infect human glioma (U251) cells. These results indicated that E30 affects the nervous system, and they provide a theoretical basis for understanding its pathogenesis. IMPORTANCE Echovirus 30 (E30) infection causes a wide spectrum of diseases with mild symptoms, such as hand, foot, and mouth disease (HFMD), acute flaccid paralysis, and aseptic meningitis and other diseases, especially one of the most common pathogens causing aseptic meningitis outbreaks. We established a novel mouse model of E30 infection by inoculating neonatal mice with clinical isolates of E30 and observed the pathological changes induced by E30. Using the E30 infection model, we found complement responses and neuropathy-related genes in the mice tissues at the transcriptome level. Moreover, we found that the viral dsRNA localized in the mouse brain and could replicate in human glioma cell line U251 rather than in the neuroblastoma cell line, SK-N-SH.

Whole-genome analysis of coxsackievirus B3 reflects its genetic diversity in China and worldwide.

BACKGROUND: Coxsackievirus B3 (CVB3) has emerged as an active pathogen in myocarditis, aseptic meningitis, hand, foot, and mouth disease (HFMD), and pancreatitis, and is a heavy burden on public health. However, CVB3 has not been systematically analyzed with regard to whole-genome diversity and recombination. Therefore, this study was undertaken to systematically examine the genetic characteristics of CVB3 based on its whole genome. METHODS: We combined CVB3 isolates from our national HFMD surveillance and global sequences retrieved from GenBank. Phylogenetic analysis was performed to examine the whole genome variety and recombination forms of CVB3 in China and worldwide. RESULTS: Phylogenetic analysis showed that CVB3 strains isolated worldwide could be classified into clusters A-E based on the sequence of the entire VP1 region. The predominant CVB3 strains in China belonged to cluster D, whereas cluster E CVB3 might be circulated globally compared to other clusters. The average nucleotide substitution rate in the P1 region of CVB3 was 4.82 × 10-3 substitutions/site/year. Myocarditis was more common with cluster A. Clusters C and D presented more cases of acute flaccid paralysis, and cluster D may be more likely to cause HFMD. Multiple recombination events were detected among CVB3 variants, and there were twenty-three recombinant lineages of CVB3 circulating worldwide. CONCLUSIONS: Overall, this study provides full-length genomic sequences of CVB3 isolates with a wide geographic distribution over a long-term time scale in China, which will be helpful for understanding the evolution of this pathogen. Simultaneously, continuous surveillance of CVB3 is indispensable to determine its genetic diversity in China as well as worldwide.

Coxsackievirus B4: an underestimated pathogen associated with a hand, foot, and mouth disease outbreak.

In order to discover the causes of a coxsackievirus B4 (CV-B4)-associated hand, foot, and mouth disease (HFMD) outbreak and to study the evolutionary characteristics of the virus, we sequenced isolates obtained during an outbreak for comparative analysis with previously sequenced strains. Phylogenetic and evolutionary dynamics analysis was performed to examine the genetic characteristics of CV-B4 in China and worldwide. Phylogenetic analysis showed that CV-B4 originated from a common ancestor in Shandong. CV-B4 strains isolated worldwide could be classified into genotypes A-E based on the sequence of the VP1 region. All CV-B4 strains in China belonged to genotype E. The global population diversity of CV-B4 fluctuated substantially over time, and CV-B4 isolated in China accounted for a significant increase in the diversity of CV-B4. The average nucleotide substitution rate in VP1 of Chinese CV-B4 (5.20 × 10-3 substitutions/site/year) was slightly higher than that of global CV-B4 (4.82 × 10-3 substitutions/site/year). This study is the first to investigate the evolutionary dynamics of CV-B4 and its association with an HFMD outbreak. These findings explain both the 2011 outbreak and the global increase in CV-B4 diversity. In addition to improving our understanding of a major outbreak, these findings provide a basis for the development of surveillance strategies.

Two Coxsackievirus B3 outbreaks associated with hand, foot, and mouth disease in China and the evolutionary history worldwide.

BACKGROUND: Coxsackievirus B3 (CV-B3) is usually associated with aseptic meningitis and myocarditis; however, the association between CV-B3 and hand, foot, and mouth disease (HFMD) has not been clearly demonstrated, and the phylogenetic dynamics and transmission history of CV-B3 have not been well summarized. METHOD: Two HFMD outbreaks caused by CV-B3 were described in Hebei Province in 2012 and in Shandong Province in 2016 in China. To analyze the epidemiological features of two CV-B3 outbreaks, a retrospective analysis was conducted. All clinical specimens from CV-B3 outbreaks were collected and disposed according to the standard procedures supported by the WHO Global Poliovirus Specialized Laboratory. EV genotyping and phylogenetic analysis were performed to illustrate the genetic characteristics of CV-B3 in China and worldwide. RESULTS: Two transmissible lineages (lineage 2 and 3) were observed in Northern China, which acted as an important "reservoir" for the transmission of CV-B3. Sporadic exporting and importing of cases were observed in almost all regions. In addition, the global sequences of CV-B3 showed a tendency of geographic-specific clustering, indicating that geographic-driven adaptation plays a major role in the diversification and evolution of CV-B3. CONCLUSIONS: Overall, our study indicated that CV-B3 is a causative agent of HFMD outbreak and revealed the phylogenetic dynamics of CV-B3 worldwide, as well as provided an insight on CV-B3 outbreaks for effective intervention and countermeasures.

Synergetic association between coxsackievirus A16 genotype evolution and recombinant form shifts.

Coxsackievirus A16 (CVA16) is a major pathogen that causes hand, foot, and mouth disease (HFMD). The recombination form (RF) shifts and global transmission dynamics of CVA16 remain unknown. In this retrospective study, global sequences of CVA16 were retrieved from the GenBank database and analyzed using comprehensive phylogenetic inference, RF surveys, and population structure. A total of 1,663 sequences were collected, forming a 442-sequences dataset for VP1 coding region analysis and a 345-sequences dataset for RF identification. Based on the VP1 coding region used for serotyping, three genotypes (A, B, and D), two subgenotypes of genotype B (B1 and B2), and three clusters of subgenotype B1 (B1a, B1b, and B1c) were identified. Cluster B1b has dominated the global epidemics, B2 disappeared in 2000, and D is an emerging genotype dating back to August 2002. Globally, four oscillation phases of CVA16 evolution, with a peak in 2013, and three migration pathways were identified. Europe, China, and Japan have served as the seeds for the global transmission of CVA16. Based on the 3D coding region of the RFs, five clusters of RFs (RF-A to -E) were identified. The shift in RFs from RF-B and RF-C to RF-D was accompanied by a change in genotype from B2 to B1a and B1c and then to B1b. In conclusion, the evolution and population dynamics of CVA16, especially the coevolution of 3D and VP1 genes, revealed that genotype evolution and RF replacement were synergistic rather than stochastic.

Evolution, recombination and geographic spreading of global Coxsackievirus A6.

BACKGROUND: The increasing incidence of hand, foot, and mouth disease (HFMD) associated with Coxsackievirus A6 (CVA6) has become a very significant public health problem. The aim of this study is to investigate the recombination, geographic transmission, and evolutionary characteristics of the global CVA6. METHODS: From 2019 to 2022, 73 full-length CVA6 sequences were obtained from HFMD patients in China and analyzed in combination with 1032 published whole genome sequences. Based on this dataset, the phylogenetic features, recombinant diversity, Bayesian phylodynamic characteristics, and key amino acid variations in CVA6 were analyzed. RESULTS: The four genotypes of CVA6, A, D, E, and F, are divided into 24 recombinant forms (RFs, RF-A - RF-X) based on differences in the P3 coding region. The eastern China region plays a key role in the dissemination of CVA6 in China. VP1-137 and VP1-138 are located in the DE loop on the surface of the CVA6 VP1 protein, with the former being a highly variable site and the latter having more non-synonymous substitutions. CONCLUSIONS: Based on whole genome sequences, this study contributes to the CVA6 monitoring, early warning, and the pathogenic mechanism by studying recombination diversity, geographical transmission characteristics, and the variation of important amino acid sites.

Genomic Epidemiology and Transmission Dynamics of Global Coxsackievirus B4.

The aim of this study was to determine the global genetic diversity and transmission dynamics of coxsackievirus B4 (CVB4) and to propose future directions for disease surveillance. Next-generation sequencing was performed to obtain the complete genome sequence of CVB4, and the genetic diversity and transmission dynamics of CVB4 worldwide were analyzed using bioinformatics methods such as phylogenetic analysis, evolutionary dynamics, and phylogeographic analysis. Forty complete genomes of CVB4 were identified from asymptomatic infected individuals and hand, foot, and mouth disease (HFMD) patients. Frequent recombination between CVB4 and EV-B multiple serotypes in the 3Dpol region was found and formed 12 recombinant patterns (A-L). Among these, the CVB4 isolated from asymptomatic infected persons and HFMD patients belonged to lineages H and I, respectively. Transmission dynamics analysis based on the VP1 region revealed that CVB4 epidemics in countries outside China were dominated by the D genotype, whereas the E genotype was dominant in China, and both genotypes evolved at a rate of > 6.50 × 10-3 substitutions/site/year. CVB4 spreads through the population unseen, with the risk of disease outbreaks persisting as susceptible individuals accumulate. Our findings add to publicly available CVB4 genomic sequence data and deepen our understanding of CVB4 molecular epidemiology.

Epidemiology of hand, foot, and mouth disease and the genetic characteristics of Coxsackievirus A16 in Taiyuan, Shanxi, China from 2010 to 2021.

Hand, foot, and mouth disease (HFMD) is a common childhood infectious disease caused by human enteroviruses (EV). This study aimed to describe the epidemiological features of HFMD and the genetic characteristics of Coxsackievirus A16 (CVA16) in Taiyuan, Shanxi, China, from 2010 to 2021. Descriptive epidemiological methods were used to analyze the time and population distribution of HFMD and the genetic characteristics of CVA16. Except being affected by the COVID-19 epidemic in 2020, HFMD epidemics were sporadic from January to March each year, and began to increase in April, with a major epidemic peak from May to August, which declined in September, followed by a secondary peak from October to December. The prevalence of EV infection was the highest in children aged one to five years (84.42%), whereas its incidence was very low in children under one year of age (5.48%). Enterovirus nucleic acid was detected by real-time reverse transcription polymerase chain reaction in 6641 clinical specimens collected from patients with HFMD from 2010 to 2021, and 4236 EV-positive specimens were detected, including 988 enterovirus A71 (EV-A71), 1488 CVA16, and 1760 other enteroviruses. CVA16 remains prevalent and has co-circulated with other EVs in Taiyuan from 2010 to 2021. A phylogenetic tree constructed based on the VP1 region showed that all CVA16 strains belonged to two different clades of the B1 genotype, B1a and B1b. They showed a nucleotide similarity of 86.5-100%, and an amino acid similarity of 96.9-100%. Overall, these findings add to the global genetic resources of CVA16, demonstrate the epidemiological characteristics of HFMD as well as the genetic features of CVA16 in Taiyuan City during 2010-2021, and provide supporting evidence for the prevention and control of HFMD.

Molecular Characteristics and Genetic Evolution of Echovirus 33 in Mainland of China.

Echovirus, a member of the Enterovirus B (EV-B) family, has led to numerous outbreaks and pandemics, causing a broad spectrum of diseases. Based on the national hand, foot, and mouth disease (HFMD) surveillance system, seven strains of echovirus 33 (E33) were isolated from Mainland of China between 2010 and 2018. The whole genomes of these strains were isolated and sequenced, and phylogenetic trees were constructed based on the gene sequences in different regions of the EV-B prototype strains. It was found that E33 may be recombined in the P2 and P3 regions. Five genotypes (A-E) were defined based on the entire VP1 region of E33, of which the C gene subtype was the dominant gene subtype at present. Recombinant analysis showed that genotype C strains likely recombined with EV-B80, EV-B85, E13, and CVA9 in the P2 and P3 regions, while genotype E had the possibility of recombination with CVB3, E3, E6, and E4. Results of Bayesian analysis indicated that E33 may have appeared around 1955 (95% confidence interval: 1945-1959), with a high evolutionary rate of 1.11 × 10-2 substitution/site/year (95% highest posterior density (HPD): 8.17 × 10-3 to 1.4 × 10-2 substitution/site/year). According to spatial transmission route analysis, two significant transmission routes were identified: from Australia to India and from Oman to Thailand, which the E33 strain in Mainland of China likely introduced from Mexico and India. In conclusion, our study fills the gaps in the evolutionary analysis of E33 and can provide important data for enterovirus surveillance.

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 3Dpol 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 3Dpol, 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 3Dpol, 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.

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.

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.

Genetic recombination in fast-spreading coxsackievirus A6 variants: a potential role in evolution and pathogenicity.

Hand, foot, and mouth disease (HFMD) is a common global epidemic. From 2008 onwards, many HFMD outbreaks caused by coxsackievirus A6 (CV-A6) have been reported worldwide. Since 2013, with a dramatically increasing number of CV-A6-related HFMD cases, CV-A6 has become the predominant HFMD pathogen in mainland China. Phylogenetic analysis based on the VP1 capsid gene revealed that subtype D3 dominated the CV-A6 outbreaks. Here, we performed a large-scale (near) full-length genetic analysis of global and Chinese CV-A6 variants, including 158 newly sequenced samples collected extensively in mainland China between 2010 and 2018. During the global transmission of subtype D3 of CV-A6, the noncapsid gene continued recombining, giving rise to a series of viable recombinant hybrids designated evolutionary lineages, and each lineage displayed internal consistency in both genetic and epidemiological features. The emergence of lineage -A since 2005 has triggered CV-A6 outbreaks worldwide, with a rate of evolution estimated at 4.17 × 10-3 substitutions site- 1 year-1 based on a large number of monophyletic open reading frame sequences, and created a series of lineages chronologically through varied noncapsid recombination events. In mainland China, lineage -A has generated another two novel widespread lineages (-J and -L) through recombination within the enterovirus A gene pool, with robust estimates of occurrence time. Lineage -A, -J, and -L infections presented dissimilar clinical manifestations, indicating that the conservation of the CV-A6 capsid gene resulted in high transmissibility, but the lineage-specific noncapsid gene might influence pathogenicity. Potentially important amino acid substitutions were further predicted among CV-A6 variants. The evolutionary phenomenon of noncapsid polymorphism within the same subtype observed in CV-A6 was uncommon in other leading HFMD pathogens; such frequent recombination happened in fast-spreading CV-A6, indicating that the recovery of deleterious genomes may still be ongoing within CV-A6 quasispecies. CV-A6-related HFMD outbreaks have caused a significant public health burden and pose a great threat to children's health; therefore, further surveillance is greatly needed to understand the full genetic diversity of CV-A6 in mainland China.

Genomic epidemiology of coxsackievirus A16 in mainland of China, 2000-18.

Hand, foot, and mouth disease (HFMD), which is a frequently reported and concerning disease worldwide, is a severe burden on societies globally, especially in the countries of East and Southeast Asia. Coxsackievirus A16 (CV-A16) is one of the most important causes of HFMD and a severe threat to human health, especially in children under 5 years of age. To investigate the epidemiological characteristics, spread dynamics, recombinant forms (RFs), and other features of CV-A16, we leveraged the continuous surveillance data of CV-A16-related HFMD cases collected over an 18-year period. With the advent of the EV-A71 vaccine since 2016, which targeted the EV-A71-related HFMD cases, EV-A71-related HFMD cases decreased dramatically, whereas the CV-A16-related HFMD cases showed an upward trend from 2017 to October 2019. The CV-A16 strains observed in this study were genetically related and widely distributed in the mainland of China. Our results show that three clusters (B1a-B1c) existed in the mainland of China and that the cluster of B1b dominates the diffusion of CV-A16 in China. We found that eastern China played a decisive role in seeding the diffusion of CV-A16 in China, with a more complex and variant transmission trend. Although EV-A71 vaccine was launched in China in 2016, it did not affect the genetic diversity of CV-A16, and its genetic diversity did not decline, which confirmed the epidemiological surveillance trend of CV-A16. Two discontinuous clusters (2000-13 and 2014-18) were observed in the full-length genome and arranged along the time gradient, which revealed the reason why the relative genetic diversity of CV-A16 increased and experienced more complex fluctuation model after 2014. In addition, the switch from RFs B (RF-B) and RF-C co-circulation to RF-D contributes to the prevalence of B1b cluster in China after 2008. The correlation between genotype and RFs partially explained the current prevalence of B1b. This study provides unprecedented full-length genomic sequences of CV-A16 in China, with a wider geographic distribution and a long-term time scale. The study presents valuable information about CV-A16, aimed at developing effective control strategies, as well as a call for a more robust surveillance system, especially in the Asia-Pacific region.