Letter to the Editor

Hand, foot, and mouth disease outbreak caused by coxsackievirus A6, China, 2013

To investigate the evolution and genetic characteristics of Coxsackievirus A6 (CVA6) which acted as the predominant pathogen of hand, foot and mouth disease (HFMD) in children in Beijing, China, 2017–2019.

Throat swab specimens were collected for general Enterovirus (EV), enterovirus A71 (EV-A71) and CVA16 detection by Real-time PCR. These general EV-positive samples were identified by semi-nested RT-PCR method and sequencing. The CVA6 VP1 gene and genome sequences were amplified and sequenced. The phylogenetic, variation and recombination analyses were performed.

A total of 1721 HFMD patients were enrolled in this study, with the male to female ratio of 1.62:1. The majority of cases were less than five years, which accounted for 73.50%. The overall detection rate of EV was 88.32% (1520/1721). A total of 8 EV types were identified, including CVA6 (55.86%), CVA16 (26.32%), EV-A71 (2.24%), CVA10 (2.04%), CVA4 (1.05%), CVA5 (0.59%), CVA2 (0.33%), and CVA8 (0.07%), while 175 (11.51%) EV were untyped. The main pathogen of HFMD was CVA6 from 2017 to 2018, while CVA6 and CVA16 were the main causative pathogens in 2019. The nucleotide and amino acid sequence identities of the 120 CVA6 complete VP1 gene sequences in this study were 91.2%–100.0% and 97.7%–100.0%, respectively. Compared with the prototype strain (Gdula) of CVA6, the nucleotide and amino acid sequence identities were 81.7%–84% and 94.7%–96.3%, respectively. The phylogenetic tree indicated that all 120 CVA6 sequences belonged to sub-genotype D3, while 119 CVA6 sequences belonged to evolutionary branch D3a, except one from 2017 belonged to D3b. Recombination analysis based on the complete genome sequences showed that potential multiple recombination may have occurred in 2B and 3D protein coding regions with EV-A114.

The main pathogens of HFMD were CVA6 and CVA16 in Beijing, China, 2017–2019; while these CVA6, as recombination strains, belonged to the D3a evolutionary branch.

Hand, foot, and mouth disease (HFMD) is an important public health problem. A monovalent EV-A71 vaccine was launched in China in 2016. Previous studies showed that inactivated monovalent EV-A71 vaccines were highly efficient against HFMD associated with EV-A71 but not against HFMD with other etiologies, leading to a hypothesis that the introduction of EV-A71 vaccines might change the pathogen spectrum and epidemiological trend of HFMD. In this study, we described for the first time the changing epidemiological characteristics of HFMD after the launch of the EV-A71 vaccine.

We extracted individual-based epidemiological data on HFMD cases reported to the Chinese Center for Disease Control and Prevention between January 2013 and December 2019. We described the changing epidemiological characteristics of HFMD before and after vaccine launch according to the distribution of diseases characteristics (demographic, temporal, and geographical) and evaluated the potential changes in risk factors of severe patients. All analyses were stratified by the phase before and after vaccine launch, and by enterovirus serotype.

During 2013−2019, 15,316,710 probable cases of HFMD were reported. Of these, 787,197 (5·1%) were laboratory confirmed and 76,982 (0·5%) were severe. After the launch of the EV-A71 vaccine, the median age of HFMD patients infected with EV-A71 increased from 2·24 years (IQR:1·43, 3·56) to 2·81 years (IQR:1·58, 4·01). The proportion of patients less than 3 years of age decreased while the proportion of patients 3−5 years of age increased. There was a large decrease (60·7%) in the proportion of severe cases as well as a decline (28·3%) in HFMD patients infected with EV-A71. After the launch of the EV-A71 vaccine, the severe illness rate and mortality rate of HFMD patients in all age groups has decreased sharply, 62·20% and 83·78% respectively. The timing of the HFMD epidemic peak was delayed (1-2 months) . After the launch of EV-A71 vaccine, the risk of becoming a severe case for EV-A71 serotype was decreased, whereas that risk was instead increased for CV-A16 (from 0·17 (95% CI:0·16, 0·18) to 0·23 (95% CI:0·21, 0·25)) and other enterovirus compared to EV-A71 (from 0·38 (95% CI:0·37, 0·39) to 0·58 (95% CI:0·56, 0·61)). The longer the time from onset to diagnosis, the higher was the risk of being a severe case, but the effect size was decreased.

The introduction of the EV-A71 vaccine has effectively reduced the proportion of severe HFMD cases and mortality, but changes to the dominant serotypes should be closely monitored. Development of multivalent vaccines to avoid an increased case burden due to other enteroviruses is greatly needed.

This research was supported by the National Natural Science Foundation of China (81973102, 81773487), Public Health Talents Training Program of Shanghai Municipality (GWV-10.2-XD21), the 5th Three-year Action Program of Shanghai Municipality for Strengthening the Construction of Public Health System (GWV-10.1-XK05), the Major Project of Scientific and Technical Winter Olympics from National Key Research and Development Program of China (2021YFF0306000), 13th Five-Year National Science and Technology Major Project for Infectious Diseases (2018ZX10725-509) and Key projects of the PLA logistics Scientific research Program (BHJ17J013).

Repeated outbreaks of Hand, foot and mouth disease (HFMD) infections have been observed in recent decades and dominated by various enteroviral serotypes. In particular, enterovirus 71 (EV-A71), coxsackievirus A16 (CV-A16) and coxsackievirus A6 (CV-A6) dominated the prevalence of HFMD infections alternatively in recent years with various outbreak sizes in Baoji, a city of Shaanxi Province in Northwest China. Estimating the reproduction number for various enteroviruses serotypes in northwest China (north temperate zone) and identification of cyclicity of HFMD infections are therefore an issue of great importance for future epidemics prediction and control. The basic/effective reproduction numbers for EV-A71, CV-A16 and CV-A6 were estimated based on daily new cases in 2010, 2011 and 2018, respectively, in which the corresponding pathogen dominated the epidemic. Two different methods based on serial interval were adopted and the basic reproduction number were estimated to be in the range of (1.33, 1.46) for CV-A16, (1.20, 1.29) for EV-A71, and (1.38, 1.59) for CV-A6, respectively. The estimated daily effective reproduction numbers significantly fluctuated before June or after July but varied mildly in (0.5,2) in around June to July for three serotypes. The weekly effective reproduction number for HFMD was estimated based on weekly new cases from year 2010 to 2018, and in most years it peaked in the range of (1.6,2.0) in February to March as well as in the range of (1.0,1.2) in September to October. The wavelet analysis based on the time series of HFMD cases from 2008 to 2018 showed obvious annual and semi-annual cyclicity, while the inter-annual cycles are infeasible. In this study we found that CV-A6 shows the greatest transmission ability among these three pathogens while EV-A71 exhibits the weakest ability of transmission, and moreover, the estimated values of basic reproduction number in northwest China are lower than those in Singapore, Hongkong and Guangdong, which may be due to different climatic circumstances.