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

Clinical features of hand, foot and mouth disease caused by Coxsackievirus A6 in Xi'an, China, 2013-2019: A multicenter observational study

PURPOSE

To investigate the clinical features of hand, foot, and mouth disease (HFMD) caused by coxsackievirus A6 (CVA6) and this work may help early diagnose of atypical HFMD.

METHODS

From January 2013 to December 2019, a total of 7,208 patients with a clinical diagnosis of HFMD in Xi'an Children's Hospital, Xi'an Central Hospital, and Xi'an Jiaotong University Second Affiliated Hospital, were included in this observational study. The clinical data, specimens and follow-up results were collected. Real-time RT‒PCR was performed for the detection and typing of enterovirus nucleic acids.

RESULTS

Of the 7,208 clinically diagnosed HFMD patients, 5,622 were positive for enterovirus nucleic acids, and the positive proportions of CVA6, enterovirus 71 (EV-A71), coxsackievirus A16 (CVA16), and other enteroviruses were 31.0% (1,742/5,622), 27.0% (1,518/5,622), 35.0% (1,968/5,622), and 7.0% (394/5,622), respectively. Based on the etiology, patients were divided into CVA6 group, EV-A71group, and CVA16 group. The mean age at onset was significantly higher in the CVA6 group (4.62±2.13 years) than in the EV-A71 group and CVA16 group (3.45±2.25 years and 3.35±2.13 years, respectively; both P < 0.05). The male/female ratio was 1.45 (1,031/711) in the CVA6 group and was not significantly different from the other two groups. The incidence of fever was significantly higher in the CVA6 group [82.5% (1,437/1,742)] than in the EV-A71 group [51.3% (779/1,518)] and the CVA16 group [45.9% (903/1,968)] (P < 0.05). In the CVA6 group, the rashes were more frequently on the trunk and elbows/knees and were significantly different from the other two groups (P < 0.05). The number of patients with two or more rash morphologies was significantly higher in the CVA6 group than in the other two groups (P < 0.05). The incidence of bullous rash in the CVA6 group [20.2%; n = 352] was higher than in the EV-A71 group [0.33%; n = 5] and CVA16 group [0.66%; n = 13] (P < 0.05). The incidence of neurological complications was significantly higher in the EV-A71 group [52.1% (791/1,518)] than in the CVA16 group [5.1% (100/1,968)] and the CVA6 group [0.8% (14/1,742)] (P < 0.05). In the follow-up period, 160 patients (9.2%) with CVA6 HFMD experienced onychomadesis, but no onychomadesis was observed in the EV-A71 and CVA16 groups. The average WBC count was significantly higher in the CVA6 group than in the CVA16 group (P < 0.05). The number of patients with increased CRP was significantly larger in the CVA6 group than in the CVA16 group but was significantly smaller than that in the EV-A71 group (P < 0.05).

CONCLUSIONS

CVA6 has become one of the main pathogens of HFMD in the Xi'an area during 2013-2019. The main clinical manifestations were slightly different from those of HFMD caused by EV-A71 or CVA16, with a higher frequency of fever, diverse morphologies and diffuse distribution of rashes, fewer neurological complications and some onychomadesis.

A multiplex one-step fluorescence quantitative differential diagnosis method for severe hand, foot and mouth disease caused by coxsackievirus A16

Hand foot and mouth disease (HFMD) is a common childhood infectious disease which is caused by human enterovirus. The objective of this study was to develop a rapid, sensitive, and accurate method for detecting severe HFMD caused by coxsackievirus A16 (CV-A16). A closed-tube sensitive multiplex one-step reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) was applied to detect CV-A16 in the early stage of severe HFMD. This assay targeted the CV-A16 structure protein VP1 to distinguish CV-A16 from other coxsackieviruses The 5'UTR region of enteric viruses was used for detecting the enterovirus and ribonuclease P (RNaseP) was adopted as the internal reference gene. The multiplex MGB probe assay system was used to detect PCR amplicons with different fluorescence reporters in the same system. The limit of detection (LOD) of the RT-qPCR assay for the CV-A16 VP1 gene was 125.893 copies/μl, for the 5' UTR was 50.1187 copies/μl and for the RNaseP gene was 158.49 copies/μl. Furthermore, specificity analysis showed that the multiplex RT-PCR had no cross-reactivity with the influenza virus, herpangina virus and SARS-COV-2. In correlation analysis, the sensitivity of the multiplex RT-qPCR assay for CV-A16 detection was 100 % (288/288) and the specificity of the multiplex RT-qPCR assay was 99.94 % (3395/3397). The overall agreement between the multiplex RT-qPCR and the results of clinical diagnosis was 99.95 % (3683/3685) and kappa value was 0.996 (p<0.001). The entire procedure, from specimen processing to result reporting, could be completed within 1.5 hours. The one-step multiplex RT-qPCR assay for detecting CV-A16 developed in this study is a good laboratory diagnostic tool for rapid and reliable distinguished detection of CV-A16, especially for severe HFMD patients at an early stage in the disease with low virus load of CV-A16.

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

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

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

Introduction

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

Methods

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

Results

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

Discussion

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

Molecular evolutionary dynamics of enterovirus A71, coxsackievirus A16 and coxsackievirus A6 causing hand, foot and mouth disease in Thailand, 2000-2022

Hand, foot and mouth disease (HFMD) is a public health threat worldwide, particularly in the Asia-Pacific region. Enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16), and CVA6 are the major pathogens causing HFMD outbreaks in several countries, including Thailand. We retrieved 385 VP1 nucleotide sequences, comprising 228 EV-A71, 33 CVA16, and 124 CVA6, deposited in the databases between 2000 and 2022 for molecular evolutionary characterization using Bayesian phylogeny. All EV-A71 identified belonged to genotype B, subgenotypes B4, and B5, and to genotype C, subgenotypes C1, C2, C4a, C4b, and C5. The analyzes demonstrated these viruses' co-circulation and subgenotypic changes throughout the past two decades. The CVA16 was grouped in genotype B1, predominantly subgenotype B1a, and the CVA6 was grouped in subgenotype D3, clades 1-4. The tMRCA of EV-A71 genotypes B and C, CVA16 B1, and CVA6 D3 dated 1993.79, 1982.62, 1995.86, and 2007.31, respectively, suggesting that the viruses were likely introduced and cryptically circulated in Thailand before the HFMD cases were recognized. We demonstrated these viruses' fluctuation and cyclical pattern throughout the two decades of observation. This study provided insight into evolutionary dynamics concerning molecular epidemiology and supported the selection of current genotype-matched vaccines, vaccine development, and implementation.

Highly sensitive and rapid identification of coxsackievirus A16 based on reverse transcription multiple cross displacement amplification combined with nanoparticle-based lateral flow biosensor assay

Introduction

One of the main pathogens responsible for human hand, foot, and mouth disease (HFMD), coxsackievirus A16, has put young children’s health at danger, especially in countries in the Asia-Pacific region. Early quick identification is essential for the avoidance and control of the disorder since there are no vaccinations or antiviral medications available to prevent and manage CVA16 infection.

Methods

Here, we describe the creation of an easy, speedy, and accurate CVA16 infection detection approach using lateral flow biosensors (LFB) and reverse transcriptionmultiple cross displacement amplification (RT-MCDA). A group of 10 primers was developed for the RT-MCDA system in order to amplify the genes in an isothermal amplification device while targeting the highly conserved region of the CVA16 VP1 gene. Then, without requiring any extra tools, RT-MCDA amplification reaction products might well be detected by visual detection reagent (VDR) and LFB.

Results

The outcomes showed that 64°C within 40 min was the ideal reaction setting for the CVA16-MCDA test. Target sequences with &lt;40 copies might be found using the CVA16-MCDA. There was no cross-reaction among CVA16 strains and other strains. The findings demonstrated that the CVA16-MCDA test could promptly and successfully identify all of the CVA16-positive (46/220) samples identified by the traditional real-time quantitative polymerase chain reaction (qRT-PCR) assays for 220 clinical anal swab samples. The whole process, such as the processing of the sample (15 min), the MCDA reaction (40 min), and the documenting of the results (2 min), could be finished in 1 h.

Conclusion

The CVA16-MCDA-LFB assay, which targeted the VP1 gene, was an efficient, simple, and highly specific examination that might be used extensively in rural regions’ basic healthcare institutions and point-of-care settings.

Epidemiological and etiological characteristics of mild hand, foot and mouth disease in children under 7 years old, Nanjing, China, 2010-2019

Background

Mild hand, foot and mouth disease (HFMD) cases make up a relatively high proportion of HFMD while have often been overlooked. This study aimed to investigate the epidemiological and etiological characteristics of mild HFMD in Nanjing.

Methods

Data on mild HFMD cases, during 2010–2019 in Nanjing, were collected from the China Information System for Disease Control and Prevention. This study mainly focused on mild cases aged < 7 years. Descriptive analysis was used to summarize epidemiological and etiological characteristics of mild cases. Flexible spatial scan statistic was used to detect spatial clusters of mild cases.

Results

A total of 175,339 mild cases aged < 7 years were reported, accounting for 94.4% of all mild cases. There was a higher average annual incidence of mild HFMD in children aged < 7 years (4,428 cases/100,000) compared with children aged ≥ 7 years (14 cases/100,000, P < 0.001), and especially children aged 1-year-old (7,908 cases/100,000). Mild cases showed semi-annual peaks of activity, including a major peak (April to July) and a minor peak (September to November). The average annual incidence was higher in males (5,040 cases/100,000) than females (3,755 cases/100,000). Based on the cumulative reported cases, the most likely cluster was detected, including Yuhuatai District, Jiangning District, Jiangbei new Area, and Pukou District. The annual distribution of enterovirus serotypes showed a significant difference. During 2010–2016, Enterovirus 71 (EV71), Coxsackievirus A16 (Cox A16), and other non-EV71/Cox A16 EVs, accounted for 29.1%, 34.6%, 36.3% of all the enterovirus test positive cases, respectively. Moreover, during 2017–2019, Cox A6, Cox A16, EV71, and other non-EV71/Cox A16/Cox A6 EVs, accounted for 47.3%, 32.5%, 10.7%, 9.5%, respectively.

Conclusions

Children under 7 years old are at higher risk of mild HFMD. Regions with high risk are mainly concentrated in the areas surrounding central urban areas. Cox A16 and Cox A6 became the dominant serotypes and they alternated or were co-epidemic. Our findings could provide valuable information for improving the regional surveillance, prevention and control strategies of HFMD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13690-022-00974-4.

Multiple molecular characteristics of circulating enterovirus types among pediatric hand, foot and mouth disease patients after EV71 vaccination campaign in Wuxi, China

The molecular properties of the circulating causative agents of hand, foot and mouth disease (HFMD) in Wuxi remain unclear, posing diagnostic and prevention challenges. Additionally, in several regions of mainland China, the EV71 immunisation drastically reduced related cases and altered the HFMD pathogen spectrum, while the precise situation in Wuxi remained unknown. To address these issues, paediatric HFMD cases diagnosed in the clinic were enrolled and anal swabs were acquired in the spring of 2019. The 5′-UTR and VP1 genes were interpreted using RT-nPCR with degenerate primers to confirm their genotypes. Following that, the entire genome sequences of each viral type were recovered, allowing for the interpretation of several molecular properties. A total of 249 clinically confirmed HFMD cases had their anal swabs taken for viral identification, from which the genome sequences of seven genotypes were recovered. Coxsackievirus A16 is the most prevalent type, followed by Coxsackievirus A6, A10, A2, A4, A5 and Echovirus 11, all of which were genetically determined for the first time in Wuxi. Phylogenetic and recombination analyses were used to evaluate their evolutionary relationships with other strains found in other regions. Noticeably, a CVA16 subtype, responsible for a large proportion of the observed cases, phylogenetically clustered within clade B1a along with some strains from other countries, was the first one to be reported in China. Furthermore, some recombination events were inferred from strains detected in sporadic cases, particularly the recombination between CVA2 and CVA5 strains. Our investigation elucidated the multiple molecular characteristics of the HFMD causal enterovirus strains in Wuxi, underlining the potential hazards associated with these circulating viral types in the population and aiding in future surveillance and prevention of this disease.

Coxsackievirus A16 in Southern Vietnam

Background: Hand, Foot and Mouth Disease (HFMD) is a major public health concern in the Asia-Pacific region. Most recent HFMD outbreaks have been caused by enterovirus A71 (EV-A71), coxsackievirus A16 (CVA16), CVA10, and CVA6. There has been no report regarding the epidemiology and genetic diversity of CVA16 in Vietnam. Such knowledge is critical to inform the development of intervention strategies.

Materials and Methods: From 2011 to 2017, clinical samples were collected from in- and outpatients enrolled in a HFMD research program conducted at three referral hospitals in Ho Chi Minh City (HCMC), Vietnam. Throat or rectal swabs positive for CVA16 with sufficient viral load were selected for whole genome sequencing and evolutionary analysis.

Results: Throughout the study period, 320 CVA16 positive samples were collected from 2808 HFMD patients (11.4%). 59.4% of patients were male. The median age was 20.8 months (IQR, 14.96–31.41). Patients resided in HCMC (55.3%), Mekong Delta (22.2%), and South East Vietnam (22.5%). 10% of CVA16 infected patients had moderately severe or severe HFMD. CVA16 positive samples from 153 patients were selected for whole genome sequencing, and 66 complete genomes were obtained. Phylogenetic analysis demonstrated that Vietnamese CVA16 strains belong to a single genogroup B1a that clusters together with isolates from China, Japan, Thailand, Malaysia, France and Australia. The CVA16 strains of the present study were circulating in Vietnam some 4 years prior to its detection in HFMD cases.

Conclusion: We report for the first time on the molecular epidemiology of CVA16 in Vietnam. Unlike EV-A71, which showed frequent replacement between subgenogroups B5 and C4 every 2–3 years in Vietnam, CVA16 displays a less pronounced genetic alternation with only subgenogroup B1a circulating in Vietnam since 2011. Our collective findings emphasize the importance of active surveillance for viral circulation in HFMD endemic countries, critical to informing outbreak response and vaccine development.

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.

Epidemiology and Sequence-Based Evolutionary Analysis of Circulating Non-Polio Enteroviruses

Enteroviruses (EVs) are positive-sense RNA viruses, with over 50,000 nucleotide sequences publicly available. While most human infections are typically associated with mild respiratory symptoms, several different EV types have also been associated with severe human disease, especially acute flaccid paralysis (AFP), particularly with endemic members of the EV-B species and two pandemic types—EV-A71 and EV-D68—that appear to be responsible for recent widespread outbreaks. Here we review the recent literature on the prevalence, characteristics, and circulation dynamics of different enterovirus types and combine this with an analysis of the sequence coverage of different EV types in public databases (e.g., the Virus Pathogen Resource). This evaluation reveals temporal and geographic differences in EV circulation and sequence distribution, highlighting recent EV outbreaks and revealing gaps in sequence coverage. Phylogenetic analysis of the EV genus shows the relatedness of different EV types. Recombination analysis of the EV-A species provides evidence for recombination as a mechanism of genomic diversification. The absence of broadly protective vaccines and effective antivirals makes human enteroviruses important pathogens of public health concern.

Molecular epidemiology of enteroviruses associated with severe hand, foot and mouth disease in Shenzhen, China, 2014-2018

In this study, we investigated the epidemiology and molecular characteristics of enteroviruses associated with severe hand, foot and mouth disease (HFMD) in Shenzhen, China, during 2014-2018. A total of 137 fecal specimens from patients with severe HFMD were collected. Enterovirus (EV) types were determined using real-time reverse transcription polymerase chain reaction (RT-PCR), RT nested PCR, and sequencing. Sequences were analyzed using bioinformatics programs. Of 137 specimens tested, 97 (70.8%), 12 (8.8%), and 10 (7.3%) were positive for EV-A71, coxsackievirus A6 (CVA6), and CVA16, respectively. Other pathogens detected included CVA2 (2.9%, 4/137), CVA10 (2.9%, 4/137), CVA5 (0.7%, 1/137), echovirus 6 (E6) (0.7%, 1/137) and E18 (0.7%, 1/137). The most frequent complication in patients with proven EV infections was myoclonic jerk, followed by aseptic encephalitis, tachypnea, and vomiting. The frequencies of vomiting and abnormal eye movements were higher in EV-A71-infected patients than that in CVA6-infected or CVA16-infected patients. Molecular phylogeny based on the complete VP1 gene revealed no association between the subgenotype of the virus and disease severity. Nevertheless, 12 significant mutations that were likely to be associated with virulence or the clinical phenotype were observed in the 5’UTR, 2A^pro, 2C, 3A, 3D^pol and 3’UTR of CVA6. Eight significant mutations were observed in the 5’UTR, 2B, 3A, 3D^pol and 3’UTR of CVA16, and 10 significant mutations were observed in the 5’UTR, VP1, 3A and 3C^pro of CVA10. In conclusion, EV-A71 is still the main pathogen causing severe HFMD, although other EV types can also cause severe complications. Potential virulence or phenotype-associated sites were identified in the genomes of CVA6, CVA16, and CVA10.