Molecular evidence of persistent epidemic and evolution of subgenotype B1 coxsackievirus A16-associated hand, foot, and mouth disease in China. J Clin Microbiol. 2010;48:619-622. [PMID: 20018819 DOI: 10.1128/JCM.02338-09]

Epidemiological and Genetic Characterization of Coxsackievirus A6-Associated Hand, Foot, and Mouth Disease in Gwangju, South Korea, in 2022

Coxsackievirus A6 (CV-A6) has emerged as the predominant causative agent of hand, foot, and mouth disease (HFMD) in young children. Since the declaration of coronavirus disease 2019 (COVID-19) as a global pandemic, the incidence of infectious diseases, including HFMD, has decreased markedly. When social mitigation was relaxed during the COVID-19 pandemic in 2022, the re-emergence of HFMD was observed in Gwangju, South Korea, and seasonal characteristics of the disease appeared to have changed. To investigate the molecular characteristics of enterovirus (EV) associated with HFMD during 2022, 277 specimens were collected. Children aged younger than 5 years accounted for the majority of affected individuals. EV detection and genotyping were performed using real-time RT-PCR and nested RT-PCR followed by sequence analysis. The EV detection rate was found to be 82.3%, and the main genotype identified was CV-A6. Sixteen CV-A6 samples were selected for whole genome sequencing. According to phylogenetic analysis, all CV-A6 strains from this study belonged to the sub-genotype D3 clade based on VP1 sequences. Analysis of 3D polymerase phylogeny showed that only the recombinant RF-A group was identified. In conclusion, circulating EV types should be continuously monitored to understand pathogen emergence and evolution during the post-pandemic era.

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.

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.

Metagenomic Surveillance of Viral Gastroenteritis in a Public Health Setting

Norovirus is the primary cause of viral gastroenteritis (GE). To investigate norovirus epidemiology, there is a need for whole-genome sequencing and reference sets consisting of complete genomes. To investigate the potential of shotgun metagenomic sequencing on the Illumina platform for whole-genome sequencing, 71 reverse transcriptase quantitative PCR (RT-qPCR) norovirus positive-feces (threshold cycle [CT], <30) samples from norovirus surveillance within The Netherlands were subjected to metagenomic sequencing. Data were analyzed through an in-house next-generation sequencing (NGS) analysis workflow. Additionally, we assessed the potential of metagenomic sequencing for the surveillance of off-target viruses that are of importance for public health, e.g., sapovirus, rotavirus A, enterovirus, parechovirus, aichivirus, adenovirus, and bocaparvovirus. A total of 60 complete and 10 partial norovirus genomes were generated, representing 7 genogroup I capsid genotypes and 12 genogroup II capsid genotypes. In addition to the norovirus genomes, the metagenomic approach yielded partial or complete genomes of other viruses for 39% of samples from children and 6.7% of samples from adults, including adenovirus 41 (N = 1); aichivirus 1 (N = 1); coxsackievirus A2 (N = 2), A4 (N = 2), A5 (N = 1), and A16 (N = 1); bocaparvovirus 1 (N = 1) and 3 (N = 1); human parechovirus 1 (N = 2) and 3 (N = 1); Rotavirus A (N = 1); and a sapovirus GI.7 (N = 1). The sapovirus GI.7 was initially not detected through RT-qPCR and warranted an update of the primer and probe set. Metagenomic sequencing on the Illumina platform robustly determines complete norovirus genomes and may be used to broaden gastroenteritis surveillance by capturing off-target enteric viruses. IMPORTANCE Viral gastroenteritis results in significant morbidity and mortality in vulnerable individuals and is primarily caused by norovirus. To investigate norovirus epidemiology, there is a need for whole-genome sequencing and reference sets consisting of full genomes. Using surveillance samples sent to the Dutch National Institute for Public Health and the Environment (RIVM), we compared metagenomics against conventional techniques, such as RT-qPCR and Sanger-sequencing, with norovirus as the target pathogen. We determined that metagenomics is a robust method to generate complete norovirus genomes, in parallel to many off-target pathogenic enteric virus genomes, thereby broadening our surveillance efforts. Moreover, we detected a sapovirus that was not detected by our validated gastroenteritis RT-qPCR panel, which exemplifies the strength of metagenomics. Our study shows that metagenomics can be used for public health gastroenteritis surveillance, the generation of reference-sets for molecular epidemiology, and how it compares to current surveillance strategies.

Phylogenetics and phylogeographic characteristics of coxsackievirus A16 in hand foot and mouth disease and herpangina cases collected in Beijing, China from 2019 to 2021

Coxsackievirus A16 (CV-A16) is a significant pathogen responsible for causing hand foot and mouth disease (HFMD) and herpangina (HA). This study aimed to investigate the recent evolution and spread of CV-A16 by monitoring HFMD and HA cases in 29 hospitals across 16 districts in Beijing from 2019 to 2021. The first five cases of HFMD and the first five cases of HA each month in each hospital were included in the study. Real-time reverse transcription polymerase chain reaction was used to identify CV-A16, CV-A6, and EV-A71. From each district, two to four CV-A16 positive samples with a relatively long sampling time interval every month were selected for sequencing. A total of 3344 HFMD cases and 2704 HA cases were enrolled in this study, with 76.0% (2541/3344) of HFMD and 45.4% (1227/2704) of HA cases confirmed to be infected by enterovirus. Among the EV-positive samples, CV-A16 virus was detected in 33.61% (854/2541) of HFMD cases and 13.4% (165/1227) of HA cases, with the predominant cluster being B1a. Both B1a and B1b had a co-circulation of local and imported strains, with different origin time (1993 vs. 1995), different global distribution (14 countries vs. 10 countries), and different transmission centers but mainly distributed in the southern and eastern regions of Beijing. Strengthening surveillance of HFMD in southern and eastern regions will improve the prevention and control efficiency of enterovirus infections.

Human SCARB2 Acts as a Cellular Associator for Helping Coxsackieviruses A10 Infection

Coxsackievirus A10 (CVA10) causes hand, foot, and mouth disease (HFMD) and herpangina, which can result in severe neurological symptoms in children. CVA10 does not use the common enterovirus 71 (EV71) receptor, human SCARB2 (hSCARB2, scavenger receptor class B, member 2), for infection but instead uses another receptor, such as KREMEN1. Our research has shown that CVA10 can infect and replicate in mouse cells expressing human SCARB2 (3T3-SCARB2) but not in the parental NIH3T3 cells, which do not express hSCARB2 for CVA10 entry. Knocking down endogenous hSCARB2 and KREMEN1 with specific siRNAs inhibited CVA10 infection in human cells. Co-immunoprecipitation confirmed that VP1, a main capsid protein where virus receptors for attaching to the host cells, could physically interact with hSCARB2 and KREMEN1 during CVA10 infection. It is the efficient virus replication following virus attachment to its cellular receptor. It resulted in severe limb paralysis and a high mortality rate in 12-day-old transgenic mice challenged with CVA10 but not in wild-type mice of the same age. Massive amounts of CVA10 accumulated in the muscles, spinal cords, and brains of the transgenic mice. Formalin inactivated CVA10 vaccine-induced protective immunity against lethal CVA10 challenge and reduced the severity of disease and tissue viral loads. This is the first report to show that hSCARB2 serves as an associate to aid CVA10 infection. hSCARB2-transgenic mice could be useful in evaluating anti-CVA10 medications and studying the pathogenesis induced by CVA10.

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.

Recent advances in anti-coxsackievirus A16 viral drug research

Hand, foot and mouth disease, a childhood disorder caused by enteroviruses, is intermittently endemic in the Asia-Pacific region and endangers the lives of many infants and young children. Coxsackievirus A16 (CV-A16) is one of the major pathogens causing hand, foot, and mouth disease on occasion, resulting in catastrophic neurological sequelae and patient death. Currently, no clinical interventions are available that completely block the CV-A16 infection. Therefore, research on anti-CV-A16 treatment continues to be a significant focus of interest. This report provides a detailed background on and an introduction to CV-A16; a description of the viral gene and protein structures and a summary of the current advances in pharmaceutical targets, drug research and other related areas.

Genetic characterization and molecular epidemiology of Coxsackievirus A12 from mainland China during 2010-2019

Coxsackievirus A12 (CVA12) is an enterovirus that has been isolated in many countries in recent years. However, studies on CVA12 are limited, and its effective population size, evolutionary dynamics and recombination patterns have not been clarified now. In this study, we described the phylogenetic characteristics of 16 CVA12 strains isolated from pediatric HFMD patients in mainland China from 2010 to 2019. Comparison of the nucleotide sequences and amino acid sequences with the CVA12 prototype strain revealed that the 16 CVA12 strains are identical in 78.8-79% and 94-94.2%, respectively. A phylodynamic analysis based on the 16 full-length VP1 sequences from this study and 21 sequences obtained from GenBank revealed a mean substitution rate of 6.61 × 10^-3 substitutions/site/year (95% HPD: 5.16-8.20 × 10^-3), dating the time to most recent common ancestor (tMRCA) of CVA12 back to 1946 (95% HPD: 1942-1947). The Bayesian skyline plot showed that the effective population size has experienced twice dynamic fluctuations since 2007. Phylogeographic analysis identified two significant migration pathways, indicating the existence of cross-provincial transmission of CVA12 in mainland China. Recombination analysis revealed two recombination patterns between 16 CVA12 strains and other EV-A, suggesting that there may be extensive genetic exchange between CVA12 and other enteroviruses. In summary, a total of 16 full-length CVA12 strains were reported in this study, providing valuable references for further studies of CVA12 worldwide.

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.

Genetic and Cross Neutralization Analyses of Coxsackievirus A16 Circulating in Taiwan from 1998 to 2021 Suggest Dominant Genotype B1 can Serve as Vaccine Candidate

Coxsackievirus A16 (CVA16) is well known for causing hand-foot-and-mouth disease (HFMD) and outbreaks were frequently reported in Taiwan in the past twenty years. The epidemiology and genetic variations of CVA16 in Taiwan from 1998 to 2021 were analyzed in this study. CVA16 infections usually occurred in early summer and early winter, and showed increased incidence in 1998, 2000-2003, 2005, 2007-2008, and 2010 in Taiwan. Little or no CVA16 was detected from 2017 to 2021. CVA16 infection was prevalent in patients between 1 to 3 years old. A total of 69 isolates were sequenced. Phylogenetic analysis based on the VP1 region showed that CVA16 subgenotype B1 was dominantly isolated in Taiwan from 1998 to 2019, and B2 was identified only from isolates collected in 1999 and 2000. There was a high frequency of synonymous mutations in the amino acid sequences of the VP1 region among CVA16 isolates, with the exception of position 145 which showed positive selection. The recombination analysis of the whole genome of CVA16 isolates indicated that the 5'-untranslated region and the non-structural protein region of CVA16 subgenotype B1 were recombined with Coxsackievirus A4 (CVA4) and enterovirus A71 (EVA71) genotype A, respectively. The recombination pattern of subgenotype B2 was similar to B1, however, the 3D region was similar to EVA71 genotype B. Cross-neutralization among CVA16 showed that mouse antisera from various subgenotypes viruses can cross-neutralize different genotype with high neutralizing antibody titers. These results suggest that the dominant CVA16 genotype B1 can serve as a vaccine candidate for CVA16.

Development and evaluation of an inactivated Coxsackievirus A16 vaccine in gerbils

Coxsackievirus A16 (CVA16) is one of the major pathogens responsible for human hand, foot, and mouth disease (HFMD), which has threatened the health of young children, particularly in Asia-Pacific nations. Vaccination is an effective strategy for protecting children from CVA16 infection. However, there is currently no licensed CVA16 vaccine for use in humans. In this study, we isolated a high-growth CVA16 virus strain in MRC-5 cells and developed an MRC-5-adapted vaccine candidate strain termed CVA16-393 via two rounds of plaque purification. The CVA16-393 strain was grouped into the B1b subgenotype and grew to a titre of over 10⁷ TCID₅₀/ml in MRC-5 cells. The VP1 gene region of this strain, which contains the major neutralizing epitopes, displayed high stability during serial passages. The inactivated whole-virus vaccine produced by the CVA16-393 strain induced an effective neutralizing antibody response in Meriones unguiculatus (gerbils) after two doses of intraperitoneal inoculation. One week after the booster immunization, the geometric mean titres of the neutralizing antibodies for the 10246, 40812TXT, 11203SD, TJ-224 and CA16-194 strains from different regions of China were 137.8, 97.8, 113.4, 64.1 and 122.3, respectively. A CVA16 vaccine dose above 25 U was also able to provide 100% cross-protection against lethal challenges with these five clinical strains in gerbils. Immunization at a one-week interval could maintain a high level of neutralizing antibody titres for at least 8 weeks. Thus, the vaccine produced by this CVA16-393 strain might be promising.

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.

The functional motions and related key residues behind the uncoating of coxsackievirus A16

The coxsackievirus A16 (CVA16) is a highly contagious virus that causes the hand, foot, and mouth disease, which seriously threatens the health of children. At present, there are still no available antiviral drugs or effective treatments against the infection of CVA16, and thus it is of great significance to develop anti-CVA16 vaccines. However, the intrinsic uncoating property of the capsid may destroy the neutralizing epitopes and influence its immunogenicity, which hinders the vaccine developments. In the present work, the functional-quantity-based elastic network model analysis method developed by our group was extended to combine with group theory to investigate the uncoating motions of the CVA16 capsid, and then the functionally key residues controlling the uncoating motions were identified by our functional-quantity-based perturbation method. Several motion modes encoded in the topological structure of the capsid were revealed to be responsible for the uncoating of CVA16 particle. These modes predominantly contribute to the fluctuation of the gyration radius of the capsid. Then, by using the perturbation method, four clusters of key sites involved in the uncoating motions were identified, whose perturbations induce significant changes in the fluctuation of the gyration radius. These key residues are mainly located at the 2-fold channels, the quasi 3-fold channels, the bottom of the canyons, and the inter-subunit interfaces around the 3-fold axes. Our studies are helpful for better understanding the uncoating mechanism of the CVA16 capsid and provide potential target sites to prevent the uncoating motions, which is valuable for the vaccine design against CVA16.

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.

Epidemiological and molecular characteristics of circulating CVA16, CVA6 strains and genotype distribution in hand, foot and mouth disease cases in 2017 to 2018 from Western India

Hand, Foot, and Mouth disease (HFMD) is a mild exanthematous and febrile disease occurs in children aged ≤10 years old. The present study highlights clinical, epidemiological characteristics, distribution of enterovirus (EV) types, and sub genotypes in HFMD cases reported during 2017 to 2018 in Western India. A total of 93 clinical samples collected from 68 HFMD cases were included. The presence of EV-RNA was determined by 5'UTR based nested reverse transcription polymerase chain reaction followed by molecular typing, sub genotyping by VP1/2A junction or VP1, full VP1 gene amplification, and phylogenetic analysis. The study reports 80.64% (75/93) EV positivity and 94.66% (71/75) typing rate, with a predominant circulation of CVA16 and CVA6 strains. Sequence analysis revealed the presence of coxsackievirus (CV)A16 (57.7%), CVA6 (40.8%), and Echo1 (1.4%) strains. EV infections were predominantly observed in children aged 1 to 3 years old (43.9%). Although cases were reported throughout the year, peaked in July (15.8%) and August (24.6%) months and persisted till September (19.3%). All the CVA16 and CVA6 positive strains were genotyped using full VP1 gene amplification. All CVA16 Indian strains (n = 41) were clustered with rarely reported B1c sub genotype and CVA6 strains (n = 29) with E2 sub-lineage. The study highlights the genetic characteristics of circulating CVA16, CVA6, and Echo1 strains in HFMD cases from Western India. The emergence of CVA16 B1c genotype and sub-lineage E2 of CVA6 strains and their constant circulation further demands systemic surveillance studies on HFMD from different parts of India to facilitate the rapid diagnosis of CVA16 and CVA6 strains using the molecular and serological based approach and for intervention strategies.

Molecular surveillance of coxsackievirus A16 in southern China, 2008-2019

A national surveillance system on hand, foot, and mouth disease (HFMD) was launched in 2008 in China. Since then, millions of HFMD cases have been reported each year, with enterovirus A71 (EV-A71), coxsackievirus A16 (CV-A16), and coxsackievirus A6 (CV-A6) as the major causative pathogens. Long-term surveillance of viral infection rates and genetic changes is essential for understanding the disease epidemiology pattern. Here, we analyzed molecular surveillance data on CV-A16 covering a period of 12 years (2008-2019) in Guangdong, China, one of the regions reporting the largest number of HFMD cases. Full VP1 sequences of 456 strains were determined to examine the genetic diversity and changes in the distribution of CV-A16 variants. Our study revealed an irregular pattern of CV-A16 infections in Guangdong. Different from the cyclic epidemics observed in some Asia-Pacific regions, there was a continuously high CV-A16 infection rate from 2008 to 2014, and after a period of lower epidemic activity in 2015-2017, an upsurge of CV-A16 infection was observed in 2018-2019. Cocirculation of subgenotypes B1a and B1b was observed, but while subgenotype B1a was predominant from 2008 to 2012, it appears to have been replaced by B1b, which has circulated as the predominant subgenotype since 2013. Phylogenetic analysis showed that most of the circulating CV-A16 strains are endemic, with occasional transmission between neighboring regions. The re-emergence of B1a in 2016-2019 in Guangdong was likely the result of introduction(s) from Southeast Asia. These results highlight the importance of continuous molecular surveillance from different areas, which will improve our understanding of the origin of the epidemic and facilitate the development of strategies for HFMD disease control.

Molecular epidemiology of coxsackievirus A16 circulating in children in Beijing, China from 2010 to 2019

BACKGROUND

Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot and mouth disease (HFMD). This study aimed to investigate the molecular epidemiology and evolutionary characteristics of CVA16.

METHODS

Throat swabs were collected from children with HFMD and suspected HFMD during 2010-2019. Enteroviruses (EVs) were detected and typed by real-time reverse transcription-polymerase chain reaction (RT-PCR) and RT-PCR. The genotype, evolutionary rate, the most recent common ancestor, population dynamics and selection pressure of CVA16 were analyzed based on viral protein gene (VP1) by bioinformatics software.

RESULTS

A total of 4709 throat swabs were screened. EVs were detected in 3180 samples and 814 were CVA16 positive. More than 81% of CVA16-positive children were under 5 years old. The prevalence of CVA16 showed obvious periodic fluctuations with a high level during 2010-2012 followed by an apparent decline during 2013-2017. However, the activities of CVA16 increased gradually during 2018-2019. All the Beijing CVA16 strains belonged to sub-genotype B1, and B1b was the dominant strain. One B1c strain was detected in Beijing for the first time in 2016. The estimated mean evolutionary rate of VP1 gene was 4.49 × 10^-3 substitution/site/year. Methionine gradually fixed at site-23 of VP1 since 2012. Two sites were detected under episodic positive selection, one of which (site-223) located in neutralizing linear epitope PEP71.

CONCLUSIONS

The dominant strains of CVA16 belonged to clade B1b and evolved in a fast evolutionary rate during 2010-2019 in Beijing. To provide more favorable data for HFMD prevention and control, it is necessary to keep attention on molecular epidemiological and evolutionary characteristics of CVA16.

Epidemical and etiological study on hand, foot and mouth disease following EV-A71 vaccination in Xiangyang, China

Coxsackievirus A6 (CV-A6) and Coxsackievirus A10 (CV-A10) have been emerging as the prevailing serotypes and overtaking Enterovirus A71 (EV-A71) and Coxsackievirus A16 (CV-A16) in most areas as main pathogens of hand, foot and mouth disease (HFMD) in China since 2013. To investigate whole etiological spectrum following EV-A71 vaccination of approximate 40,000 infants and young children in Xiangyang, enteroviruses were serotyped in 4415 HFMD cases from October 2016 to December 2017 using Real Time and conventional PCR and cell cultures. Of the typeable 3201 specimen, CV-A6 was the predominant serotype followed by CV-A16, CV-A10, CV-A5, CV-A2 and EV-A71 with proportions of 59.54%, 15.31%, 11.56%, 4.56%, 3.78% and 3.03%, respectively. Other 12 minor serotypes were also detected. The results demonstrated that six major serotypes of enteroviruses were co-circulating, including newly emerged CV-A2 and CV-A5. A dramatic decrease of EV-A71 cases was observed, whereas the total cases remained high. Multivalent vaccines against major serotypes are urgently needed for control of HFMD.

Coxsackieviruses A6 and A16 associated with hand, foot, and mouth disease in Vietnam, 2008-2017: Essential information for rational vaccine design

Development of multivalent hand, foot, and mouth disease (HFMD) vaccines against enterovirus A71 (EV-A71) and several non-EV-A71 enteroviruses is needed for this life-threatening disease with a huge economic burden in Asia-Pacific countries. Comprehensive studies on the molecular epidemiology and genetic and antigenic characterization of major causative enteroviruses will provide information for rational vaccine design. Compared with molecular studies on EV-A71, that for non-EV-A71 enteroviruses remain few and limited in Vietnam. Therefore, we conducted a 10-year study on the circulation and genetic characterization of coxsackievirus A16 (CV-A16) and CV-A6 isolated from patients with HFMD in Northern Vietnam between 2008 and 2017. Enteroviruses were detected in 2228 of 3212 enrolled patients. Of the 42 serotypes assigned, 28.4% and 22.4% accounted for CV-A6 and CV-A16, being the second and the third dominant serotypes after EV-A71 (31.7%), respectively. The circulation of CV-A16 and CV-A6 showed a wide geographic distribution and distinct periodicity. Phylogenetic analyses revealed that the majority of Vietnamese CV-A6 and CV-A16 strains were located within the largest sub-genotypes or sub-genogroups. These comprised strains isolated from patients with HFMD worldwide during the past decade and the Vietnamese strains have been evolving in a manner similar to the strains circulating worldwide. Amino acid sequences of the putative functional loops on VP1 and other VPs among Vietnamese CV-A6 and CV-A16 isolates were highly conserved. Moreover, the functional loop patterns of VP1 were similar to the dominant patterns found worldwide, except for the T164K substitution on the EF loop in Vietnamese CV-A16. The findings suggest that the development of a universal HFMD vaccine, at least in Vietnam, must target CV-A6 and CV-A16 as two of the three major HFMD-causing serotypes. Vietnamese isolates or their genome sequences can be considered for rational vaccine design.

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.

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 Variation of Multiple Serotypes of Enteroviruses Associated with Hand, Foot and Mouth Disease in Southern China

Enteroviruses (EVs) species A are a major public health issue in the Asia–Pacific region and cause frequent epidemics of hand, foot and mouth disease (HFMD) in China. Mild infections are common in children; however, HFMD can also cause severe illness that affects the central nervous system. To molecularly characterize EVs, a prospective HFMD virological surveillance program was performed in China between 2013 and 2016. Throat swabs, rectal swabs and stool samples were collected from suspected HFMD patients at participating hospitals. EVs were detected using generic real-time and nested reverse transcription-polymerase chain reactions (RT-PCRs). Then, the complete VP1 regions of enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16) and CVA6 were sequenced to analyze amino acid changes and construct a viral molecular phylogeny. Of the 2836 enrolled HFMD patients, 2,517 (89%) were EV positive. The most frequently detected EVs were CVA16 (32.5%, 819), CVA6 (31.2%, 785), and EV-A71 (20.4%, 514). The subgenogroups CVA16_B1b, CVA6_D3a and EV-A71_C4a were predominant in China and recombination was not observed in the VP1 region. Sequence analysis revealed amino acid variations at the 30, 29 and 44 positions in the VP1 region of EV-A71, CVA16 and CVA6 (compared to the respective prototype strains BrCr, G10 and Gdula), respectively. Furthermore, in 21 of 24 (87.5%) identified EV-A71 samples, a known amino acid substitution (D31N) that may enhance neurovirulence was detected. Our study provides insights about the genetic characteristics of common HFMD-associated EVs. However, the emergence and virulence of the described mutations require further investigation.

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⁻³ substitutions site^-1 year⁻¹ 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.

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.

Pathogenic characteristics of hand, foot and mouth disease in Shaanxi Province, China, 2010-2016

Hand, foot, and mouth disease (HFMD) is a common childhood illness caused by enteroviruses. We analyzed the pathogenic characteristics of HFMD in Shaanxi province, China, during 2010-2016. Clinical samples were collected from HFMD cases. Real-time PCR and RT-PCR were used to identify the enterovirus(EVs) serotypes. Viral RNA sequences were amplified using RT-PCR and compared by phylogenetic analysis. Descriptive epidemiological methods were used to analyze. A total of 16,832 HFMD positive cases were confirmed in the laboratory. EV-A71 and CV-A16 were the main pathogens in 2010. EV-A71 was the dominant pathogen in the periods of 2011 to 2012 and 2014, 2016. In 2013 and 2015, other EVs increased greatly, in which CV-A6 was the predominant pathogen. EV-A71 was more frequently detected in deaths and severe cases. Phylogenetic analysis revealed that EV-A71 belonged to the C4a evolution branch of C4 sub-genotype and CV-A16 belonged to the B1a or B1b evolution branch of B1 sub-genotype, whereas CV-A6 strains were assigned to D2 or D3 sub-genotype. The pathogen spectrum of HFMD has changed in 7 years, and the major serotypes EV-A71, CV- A16 and CV- A6 alternated or co-circulated. Long-term surveillance and research of EVs should be strengthened for the prevention and control of HFMD.

Emerging recombination of the C2 sub-genotype of HFMD-associated CV-A4 is persistently and extensively circulating in China

Sporadic outbreaks caused by coxsackievirus A4 (CV-A4) have been reported worldwide. To further elucidate the detailed genetic characteristics and evolutionary recombination events of CV-A4, virus samples from nationwide hand, foot and mouth disease (HFMD) surveillance, encompassing 27 out of the 31 provinces in China, were investigated. Comprehensive and systematic phylogenetic analyses were performed by using 29 complete genomes, 142 complete CV-A4 VP1 sequences. Four genotypes (A, B, C and D) and five sub-genotypes (C1-C5) were re-identified based on the complete VP1 sequences. C2 is the predominant sub-genotype of CV-A4 associated with HFMD and has evolved into 3 clusters. Cluster 1 is a major cluster that has been persistently and extensively circulating in China since 2006 and has been associated with all severe cases. All the sequences showed high homology with the CV-A4 prototype in the P1 region, while higher identities with CV-A5, CV-14 and CV-16 in the P2 and P3 regions. Recombination analysis revealed that C2 had two specific genetic recombination patterns with other EV-A prototypes in the 5'-UTR and 3D region compared with C5. These recombination patterns might be associated with the increased transmissibility of C2 viruses, which were obtained due to their persistent and extensive circulation in populations.

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.

Automated cell-based luminescence assay for profiling antiviral compound activity against enteroviruses

We describe the development, optimisation, and validation of an automated, cell-based and high-throughput screening assay using existing luminescence-based ATPlite reagents for identifying antiviral compounds that inhibit enterovirus replication. Antiviral efficacy was determined by measuring the ATP levels in cells that were protected from the viral cytopathic effect (CPE) by the antiviral compounds pleconaril and rupintrivir. CPE-based assay conditions were optimised at a cell density of 5000 cells/well and a viral infection dose of 100 CCID₅₀ in 384-well plates. The assay exhibited excellent robustness, with Z'-factor values between 0.75 and 0.82, coefficients of variation between 0.33% and 1.45%, and signal-to-background ratios ranging from 6.92 to 22.6 when testing three enterovirus A71 isolates circulating in China. The assay was also suitable for screening other picornaviruses, such as poliovirus, coxsackievirus, echovirus, and parechovirus.

Near-Complete Genome Sequences of 12 Coxsackievirus Group A Strains from Hand, Foot, and Mouth Disease and Herpangina Cases with Different Clinical Symptoms

Coxsackievirus group A (CV-A) strains are important pathogens of hand, foot, and mouth disease and herpangina. We report here the near-complete genome sequences of 12 CV-A strains isolated from infants and children with different clinical diseases. The presented data will be very useful for future genome-based epidemiological studies.

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.

An Environmental Surveillance in Uruguay Reveals the Presence of Highly Divergent Types of Human Enterovirus Species C and a High Frequency of Species A and B Types

Information about Human Enterovirus circulation in Uruguay is scarce. The aim of this study was to generate the first description about their circulation in the country through the study of sewage samples collected before and after the switch from Oral Poliovirus Vaccine to Inactivated Poliovirus Vaccine. Viruses were concentrated by an adsorption-elution to a negatively charged membrane, and real-time quantitative PCR and qualitative PCR methods were used to detect, quantify, and characterize enteroviruses. Positive samples were inoculated in RD cells and two passages were performed. Additionally, RD+ samples were subsequently passed onto L20B cells. Human Enteroviruses were detected in 67.6% of the samples, with concentrations between 4.9 and 6.6 Log₁₀ genomic copies per liter. 10% of positive samples replicated in RD cells, of which none in L20B cells. Molecular characterization of Human Enterovirus strains directly detected from sewage sample concentrates allowed the identification of highly divergent members of species C such as Enterovirus C99 and Coxsackievirus A13, as well as the frequent detection of species A and B members (particularly Coxsackievirus A16 and Echovirus 6, respectively). Other detected types were Coxsackievirus A2, A22, B1, B5, Echovirus 5, and 9. The characterization of viruses isolated in cell culture revealed the presence of Echovirus 6 and Coxsackievirus B3. Despite the absence of poliovirus, a wide circulation of different enterovirus types was evidenced in Uruguayan sewage samples, highlighting that the local populations are exposed to different kinds of diseases originated by several human enterovirus.

Molecular surveillance of coxsackievirus A16 reveals the emergence of a new clade in mainland China

Coxsackievirus A16 (CV-A16) of the genotypes B1a and B1b have co-circulated in mainland China in the past decades. From 2013 to 2017, a total of 3,008 specimens from 3,008 patients with mild hand, foot, and mouth disease were collected in the present study. Viral RNA was tested for CV-A16 by a real-time RT-PCR method, and complete VP1 sequences and full-length genome sequences of CV-A16 strains from this study were determined by RT-PCR and sequencing. Sequences were analyzed using a series of bioinformatics programs. The detection rate for CV-A16 was 4.1%, 25.9%, 10.6%, 28.1% and 12.9% in 2013, 2014, 2015, 2016 and 2017, respectively. Overall, the detection rate for CV-A16 was 16.5% (497/3008) in this 5-year period in Shenzhen, China. One hundred forty-two (142/155, 91.6%) of the 155 genotype B1 strains in the study belonged to subgenotype B1b, and 13 (13/155, 8.4%) strains belonged to subgenotype B1a. Two strains (CVA16/Shenzhen174/CHN/2017 and CVA16/Shenzhen189/CHN/2017) could not be assigned to a known genotype. Phylogenetic analysis of these two strains and other Chinese CV-A16 strains indicated that these two CV-A16 strains clustered independently in a novel clade whose members differed by 8.4%-11.8%, 8.4%-12.1%, and 14.6%-14.8% in their nucleotide sequences from those of Chinese B1a, B1b, and genotype D strains, respectively. Phylogenetic analysis of global CV-A16 strains further indicated that the two novel CV-A16 strains from this study grouped in a previously uncharacterized clade, which was designated as the subgenogroup B3 in present study. Meanwhile, phylogenetic reconstruction revealed two other new genotypes, B1d and B4, which included a Malaysian strain and two American strains, respectively. The complete genome sequences of the two novel CV-A16 strains showed the highest nucleotide sequence identity of 92.3% to the Malaysian strain PM-15765-00 from 2000. Comparative analysis of amino acid sequences of the two novel CV-A16 strains and their relatives suggested that variations in the nonstructural proteins may play an important role in the evolution of modern CV-A16.

Molecular characterization of Coxsackievirus A16 strains isolated from children with severe hand, foot, and mouth disease in Yunnan, Southwest China, during 2009-2015

Coxsackievirus A16 (CV-A16) commonly causes mild symptoms, but severe diseases, such as aseptic meningitis, encephalitis, and even fatal cases, have been reported. Thirteen CV-A16 strains were isolated from patients with severe hand, foot, and mouth disease in Yunnan, Southwest China, from 2009 to 2015. Subgenotype B1a and B1b of CV-A16 were predominantly circulating the region with B1b the predominant strain in recent years. The mean rate of nucleotide substitution based on the VP1 gene sequence was 4.545 × 10 ^-3  substitution per site per year from 2009 to 2015. These results may help in understanding the genetic diversity of CV-A16 and develop a CV-A16 vaccine.

Persistent circulation of genotype D coxsackievirus A2 in mainland of China since 2008

Coxsackievirus A2 (CV-A2) has emerged as an important etiological agent in the hand, foot, and mouth disease and herpangina pathogen spectrum because of its high global prevalence. In the present study, we investigated the evolutionary dynamics of CV-A2 circulating in China. We analyzed a total of 163 entire VP1 sequences of CV-A2, including 74 sequences generated from the present study and 89 sequences collected from the GenBank database. Phylogenetic analysis based on the entire VP1 nucleotide sequences confirmed the persistent circulation of the predominant genotype D in mainland of China since 2008. Cluster analysis grouped the sequences into two distinct clusters, clusters 1 and 2, with most grouped under cluster 2. After 2012, cluster 1 was gradually replaced by cluster 2. Results of Bayesian Markov chain Monte Carlo analysis suggested that multiple lineages of genotype D were transmitted in mainland of China at an estimated evolutionary rate of 6.32×10⁻³ substitutions per site per year, which is consistent with the global evolutionary rate of CV-A2 (5.82×10⁻³ substitutions per site per year). Continuous transmission and evolution of CV-A2 resulted in the genetic polymorphism.

Epidemiological characteristics and influential factors of hand, foot, and mouth disease reinfection in Wuxi, China, 2008-2016

Background

Hand, foot, and mouth disease (HFMD) is a viral disease caused by human enteroviruses. Although HFMD reinfection is common, studies investigating this phenomenon are insufficient.

Methods

The present study focused on HFMD reinfection in Wuxi from 2008 to 2016 using surveillance system data.

Results

Of 107,677 cases included in the study, 6470 cases were classified as reinfections. The overall reinfection rate was 6.01% (6.37% male and 5.48% female patients), which decreased with increasing age (χ² = 1125.477, p < 0.001). The rate was 6.17 and 5.79% in urban and rural areas, respectively, and 7.83 and 5.98% of the cases were severe and mild, respectively. Multivariate logistic regression analysis showed that male sex, younger age, residence in an urban area, and severe disease were risk factors for HFMD reinfection. The case-severity rate in secondary infection cases was lower than that in non-reinfection cases (odds ratio 0.675, 95% confidence interval 0.526–0.866).

Conclusions

Boys younger than 4 years of age living in urban areas were more prone to reinfection. Specific health education and intervention should be developed to protect these susceptible populations.

The emerging sub-genotype C2 of CoxsackievirusA10 Associated with Hand, Foot and Mouth Disease extensively circulating in mainland of China

Coxsackievirus A10 (CV-A10) associated with Hand, foot, and mouth disease (HFMD) cases emerged increasingly in recent years. In this study, the samples from nation-wide HFMD surveillance, including 27 out of 31 provinces in China were investigated, and the continuous and extensive virological surveillance, covered 13 years, were conducted to provide a comprehensive molecular characterization analysis of CV-A10. 855 CV-A10 viruses (33 severe cases included), were isolated from HFMD children patients during 2009 to 2016 in China. 164 representative sequences from these viruses, together with 117 CV-A10 sequences downloaded from GenBank based on entire VP1 were recruited in this study. Two new genotypes (F and G) and two sub-genotypes (C1 and C2) were identified. Among 264 Chinese sequences, 9 of them were genotype B, 8 of them were C1, and the other (247) were C2, the predominant sub-genotype in China since 2012. Chinese C2 viruses showed obvious temporal characteristics and can be divided into 3 clusters (cluster 1~3). Cluster 3 viruses was circulating extensively during 2014 and 2016 with more severe cases. It is very necessary and important to continuously conduct the extensive virological surveillance for CV-A10, and further evolutionary studies will provide more evidence on its evolution and virulence.

Genetic characteristics of the P1 coding region of Coxsackievirus A16 associated with hand, foot, and mouth disease in China

Coxsackievirus A16 (CVA16) is one of the major etiological agents of hand, foot, and mouth disease (HFMD) in young children. To investigate the genetic characteristics of the P1 coding region gene of CVA16 associated with HFMD in China, we included the sequences of CVA16 specimens obtained from outbreak investigations and sporadic HFMD cases between 1998 and 2014 in China from GenBank, we genotyped the CVA16 sequences and analyzed P1 coding region sequences that encode structural proteins with bioinformatics software. CVA16 was classified into genotypes A and B1 based on the VP1 gene; the B1b and B1a subgenotypes were the major CVA16 strains and predominated in the coastal areas of China. Four strains were found to show inter- and intra-typic recombination in the P1 region. The amino acid identities of VP1, VP2, VP3, and VP4 proteins in all Chinese CVA16 strains were 88.2-100%, 83.0-100%, 87.6-100%, and 72.4-100%, respectively. A total of 251 amino acid substitution sites were detected in the structural proteins encoded by the P1 coding region gene. The amino acid sequences of the P1 coding region in Chinese CVA16 strains were highly conserved, although some amino acid mutations occurred with high frequency: VP1-T11A (10%), N14S (14%), L23M/V (11%), T98M (16%), V107A (14%), N102D (6.1%), E145V (8.8%), N218D (10%), E241K (22%), T248A/I (6.8%); VP2-I217V (22%), T226A (38%); VP3-N141S/G (5.4%), and N240D (15%). The genetic characteristics of CVA16 in the P1 coding region gene may provide a basis for developing a CVA16 vaccine and preventing and controlling HFMD in China.

Epidemiological and serological surveillance of hand-foot-and-mouth disease in Shanghai, China, 2012-2016

Aside from enterovirus 71 (EV71) and coxsackie virus A16 (CV-A16), viruses that are known to cause hand-foot-and-mouth disease (HFMD), epidemiological profiles of other enteroviruses that induce HFMD are limited. We collected 9949 laboratory surveillance HFMD cases and 1230 serum samples from infants and children in Shanghai from 2012–2016. Since 2013, CV-A6 has displaced EV71 and CV-A16 to become the predominant serotype. Interestingly, novel epidemiological patterns in EV71 and CV-A16 infections were observed, with one large peak in both 2012 and 2014, followed by two smaller peaks in the respective following years (2013 and 2015). Through sequencing, we found that C4a, B1b, D-Cluster-1 and B constituted the major subgenotypes of EV71, CV-A16, CV-A6 and CV-A10, respectively. Among healthy individuals, 50.49% and 54.23% had positive neutralising antibodies (NtAbs) against EV71 and CV-A16, respectively, indicating that EV71 and CV-A16 silent infections were common. These populations may be an important potential source of infection. The overall seropositive rate of EV71 NtAbs showed a fluctuating, markedly downward trend, indicating the potential risk of a future EV71 epidemic. High CV-A16 NtAb seroprevalence corroborated a documented CV-A16 ‘silent’ epidemic. Children aged 1–5 years had the lowest EV71 NtAb seropositive rate, whereas those aged 1–2 years exhibited the lowest CV-A16 NtAb seropositive rate. This is the first comprehensive investigation of the epidemiology and aetiology, as well as the seroprevalence, of HFMD in Shanghai between 2012 and 2016. This study provides the latest insights into developing a more efficient HMFD vaccination programme.

A Selective Bottleneck Shapes the Evolutionary Mutant Spectra of Enterovirus A71 during Viral Dissemination in Humans

RNA viruses accumulate mutations to rapidly adapt to environmental changes. Enterovirus A71 (EV-A71) causes various clinical manifestations with occasional severe neurological complications. However, the mechanism by which EV-A71 evolves within the human body is unclear. Utilizing deep sequencing and haplotype analyses of viruses from various tissues of an autopsy patient, we sought to define the evolutionary pathway by which enterovirus A71 evolves fitness for invading the central nervous system in humans. Broad mutant spectra with divergent mutations were observed at the initial infection sites in the respiratory and digestive systems. After viral invasion, we identified a haplotype switch and dominant haplotype, with glycine at VP1 residue 31 (VP1-31G) in viral particles disseminated into the integumentary and central nervous systems. In vitro viral growth and fitness analyses indicated that VP1-31G conferred growth and a fitness advantage in human neuronal cells, whereas VP1-31D conferred enhanced replication in human colorectal cells. A higher proportion of VP1-31G was also found among fatal cases, suggesting that it may facilitate central nervous system infection in humans. Our data provide the first glimpse of EV-A71 quasispecies from oral tissues to the central nervous system within humans, showing broad implications for the surveillance and pathogenesis of this reemerging viral pathogen.IMPORTANCE EV-A71 continues to be a worldwide burden to public health. Although EV-A71 is the major etiological agent of hand, foot, and mouth disease, it can also cause neurological pulmonary edema, encephalitis, and even death, especially in children. Understanding selection processes enabling dissemination and accurately estimating EV-A71 diversity during invasion in humans are critical for applications in viral pathogenesis and vaccine studies. Here, we define a selection bottleneck appearing in respiratory and digestive tissues. Glycine substitution at VP1 residue 31 helps viruses break through the bottleneck and invade the central nervous system. This substitution is also advantageous for replication in neuronal cells in vitro Considering that fatal cases contain enhanced glycine substitution at VP1-31, we suggest that the increased prevalence of VP1-31G may alter viral tropism and aid central nervous system invasion. Our findings provide new insights into a dynamic mutant spectral switch active during acute viral infection with emerging viral pathogens.

Large-scale HFMD epidemics caused by Coxsackievirus A16 in Bangalore, India during 2013 and 2015

Hand foot and mouth disease (HFMD) is a relatively unreported disease in India. This study was undertaken to characterize the enterovirus type/s associated with two unexpectedly-massive epidemics that occurred in Bangalore, India in 2013 and 2015. Stool samples of 229 children with HFMD living in Northern and Southern areas of Bangalore were tested by RT-PCR; 189 (82.5%) were enterovirus positive. The Indian CV-A16 strains exhibited 98-99% sequence identity with those reported in France and China in the 5' untranslated region. BLAST and phylogenetic analyses of complete genomes of representative Indian isolates revealed that the 2015 epidemic was predominated by an inter-species recombinant between CV-A16 and coxsackievirus B5. The 2013 epidemic was primarily caused by nonrecombinant strains. The CV-A16 strains circulated in India since 2007 and phylogeographic analyses indicated imported cases in France and China. In conclusion, CV-A16-associated HFMD epidemics should be recognized as an emerging public health problem in India.

Phylogeography of Coxsackievirus A16 Reveals Global Transmission Pathways and Recent Emergence and Spread of a Recombinant Genogroup

Coxsackievirus A16 (CV-A16; Picornaviridae) is an enterovirus (EV) type associated with hand, foot, and mouth disease (HFMD) in children. To investigate the spatial spread of CV-A16, we used viral sequence data sampled during a prospective sentinel surveillance of HFMD in France (2010 to 2014) and phylogenetic reconstruction. A data set of 168 VP1 sequences was assembled with 416 publicly available sequences of various geographic origins. The CV-A16 sequences reported were assigned to two clades, genogroup B and a previously uncharacterized clade D. The time origins of clades B and D were assessed in 1978 (1973 to 1981) and 2004 (2001 to 2007), respectively. The shape of the global CV-A16 phylogeny indicated worldwide cocirculation of genetically distinct virus lineages over time and across geographic regions. Phylogenetic tree topologies and Bayes factor analysis indicated virus migration. Virus transportation events in clade B within Europe and Asia and between countries of the two geographic regions were assessed. The sustained transmission of clade D viruses over 4 years was analyzed at the township level in France and traced back to Peru in South America. Comparative genomics provided evidence of recombination between CV-A16 clades B and D and suggested an intertype recombinant origin for clade D. Time-resolved phylogenies and HFMD surveillance data indicated that CV-A16 persistence is sustained by continuing virus migration at different geographic scales, from community transmission to virus transportation between distant countries. The results showed a significant impact of virus movements on the epidemiological dynamics of HFMD that could have implications for disease prevention.IMPORTANCE Coxsackievirus A16 is one of the most prevalent enterovirus types in hand, foot, and mouth disease outbreaks reported in Southeast Asia. This study is based on epidemiological and viral data on HFMD caused by CV-A16 in a European country. The phylogeographic data complemented the syndromic surveillance with virus migration patterns between geographic regions in France. The results show how viral evolutionary dynamics and global virus spread interact to shape the worldwide pattern of an EV disease. CV-A16 transmission is driven by movements of infected individuals at different geographic levels: within a country (local dynamics), between neighboring countries (regional dynamics), and between distant countries (transcontinental dynamics). The results are consistent with our earlier data on EV-A71 and confirm the epidemiological interconnection of Asia and Europe with regard to EV infections.

Complete Genome Sequences of Four Coxsackievirus A16 Strains Isolated from Four Children with Severe Hand, Foot, and Mouth Disease

Here, we report the complete genome sequences of four coxsackievirus A16 strains isolated from four children with severe hand, foot, and mouth disease. Three of them were assigned to subgenotype B1b based on phylogenetic analysis of the VP1 gene, and the other one belonged to subgenotype B1a.

Persistent circulation of Coxsackievirus A6 of genotype D3 in mainland of China between 2008 and 2015

A total of 807 entire VP1 sequences of Coxsackievirus A6 (CV-A6) from mainland of China from 1992 to 2015, including 520 in this study and 287 from the GenBank database, were analysed to provide a basic framework of molecular epidemiological characteristics of CV-A6 in China. Sixty-five VP1 sequences including 46 representative CV-A6 isolates from 807 Chinese strains and 19 international strains from GenBank were used for describing the genotypes and sub-genotypes. The results revealed that CV-A6 strains can be categorised into 4 genotypes designated as A, B, C, and D according to previous data and can be further subdivided into B1–B2, C1–C2, and D1–D3 sub-genotypes. D3 is the predominant sub-genotype that circulated in recent years in mainland of China and represents 734 of 807 Chinese isolates. Sixty-six strains belong to D2, whereas B1 and C1 comprise a single strain each, and five AFP strains formed B2. Sub-genotype D3 first circulated in 2008 and has become the predominant sub-genotype since 2009 and then reached a peak in 2013, while D2 was mostly undetectable in the past years. These data revealed different transmission stages of CV-A6 in mainland of China and that sub-genotype D3 may have stronger transmission ability.

Prevalence of enteroviruses in healthy populations and excretion of pathogens in patients with hand, foot, and mouth disease in a highly endemic area of southwest China

Etiological carriers and the excretion of the pathogens causing hand, foot, and mouth disease (HFMD) in healthy persons, patients, and asymptomatic persons infected with HFMD as ongoing infection sources may play an important role in perpetuating and spreading epidemics of HFMD. The aims of this study were to determine the carrier status of EV-A71 and CV-A16 in healthy populations, as well as the duration of EV-A71 and CV-A16 shedding in the stools of HFMD patients in an epidemic area of southwest China. A cross-sectional study and a follow-up study were conducted in three HFMD endemic counties of Yunnan Province. Six hundred sixty-seven healthy subjects were recruited to participate in the cross-sectional study, and two stool specimens were collected from each subject. Among the healthy subjects, 90 (13.5%) tested positive for viral isolation, but neither EV-A71 nor CV-A16 was detected in healthy individuals. Of the 150 patients with probable HFMD, 55.3% (83/150) tested positive for viral isolation with presented serotypes such as EV-A71 (51.81%, 43/83), CV-A16 (32.53%, 27/83), other EVs (13.25%, 11/83), and mixed EV-A71 and CV-A16 (2.41%, 2/83). The longest duration of EV-A71 and CV-A16 shedding in stool specimens from patients with HFMD was >46 days after onset. The positive rate of EV-A71 in the stool specimens of confirmed patients dropped to 50% by the end of the third week, and the same occurred with CV-A16 by the end of approximately the seventh week after onset. Although carriers of major causative agents of HFMD in healthy populations are fewer in number, the prolonged shedding of pathogens in patients with HFMD may serve as an important factor in perpetuating and spreading HFMD epidemics.

Human picornaviruses associated with neurological diseases and their neutralization by antibodies

Picornaviruses are the most commonly encountered infectious agents in mankind. They typically cause mild infections of the gastrointestinal or respiratory tract, but sometimes also invade the central nervous system. There, they can cause severe diseases with long-term sequelae and even be lethal. The most infamous picornavirus is poliovirus, for which significant epidemics of poliomyelitis were reported from the end of the nineteenth century. A successful vaccination campaign has brought poliovirus close to eradication, but neurological diseases caused by other picornaviruses have increasingly been reported since the late 1990s. In this review we focus on enterovirus 71, coxsackievirus A16, enterovirus 68 and human parechovirus 3, which have recently drawn attention because of their links to severe neurological diseases. We discuss the clinical relevance of these viruses and the primary role of humoral immunity in controlling them, and summarize current knowledge on the neutralization of such viruses by antibodies.