ORIGINAL RESEARCH article
Molecular Epidemiology and Clinical Features of Enteroviruses-Associated Hand, Foot, and Mouth Disease and Herpangina Outbreak in Zunyi, China, 2019
- 1Department of Clinical Laboratory, Zunyi Medical University Third Affiliated Hospital, Zunyi, China
- 2Department of Pediatrics and Child Health, Zunyi Medical University Third Affiliated Hospital, Zunyi, China
- 3Department of Scientific Research Laboratory, Zunyi Medical University Third Affiliated Hospital, Zunyi, China
Background: Hand, foot and mouth disease (HFMD) and herpangina (HA), two of the most common childhood infectious diseases, are associated with enteroviruses (EVs) infection. The aim of this study was to identify the molecular epidemiology of enterovirus causing HFMD/HA in Zunyi, China, during 2019, and to describe the clinical features of the cases.
Methods: We collected the information on demographic and clinical characteristics, laboratory data of laboratory-confirmed EVs associated HFMD/HA cases in Zunyi Medical University Third Affiliated Hospital between March 1 and July 31, 2019. EV types were determined by either one-step real time RT-PCR or partial VP1 gene sequencing and sequence alignment. Phylogenetic analysis of CVA6, CVA2, and CVA5 were established based on the partial VP1 gene sequences by neighbor-joining method. Differences in clinical characteristics and laboratory results of the cases were compared among patients infected with the most prevalent EV types.
Results: From 1 March to 31 July 2019, 1,377 EVs associated HFMD/HA inpatients were confirmed. Of them, 4 (0.3%, 4/1,377) were EV-A71-associated cases, 84 (6.1%, 84/1,377) were CVA16-associated cases, and 1,289 (93.6%, 1,289/1,377) were non-EV-A71/CVA16-associated cases. Of the randomly selected 372 non-EV-A71/CVA16 cases, EV types have been successfully determined in 273 cases including 166 HFMD and 107 HA cases. For HFMD cases, the three most common types were CVA6 (80.7%, 134/166), CVA2 (5.4%, 9/166) and CVA5 (3.0%, 5/166); similarly, for HA cases, the three most prevalent serotypes were CVA6 (36.5%, 39/107), CVA2 (21.5%, 23/107) and CVA5 (18.7%, 20/107). Phylogenetic analysis showed that subclade D of CVA5, and subclade E of CVA6 and CVA2 were predominant in Zunyi during the outbreak in 2019. Compared with the cases caused by CVA16, the incidence of high fever and severe infection associated with CVA2, CVA5, and CVA6 was higher.
Conclusions: The recent HFMD/HA outbreak in Zunyi is due to a larger incidence of CVA6, CVA2, and CVA5. Novel diagnostic reagents and vaccines against these types would be important to monitor and control EV infections.
Hand, foot, and mouth disease (HFMD) and herpangina (HA) are the two major childhood infectious diseases caused by enteroviruses. HFMD/HA are recognized as self-limiting diseases, but may present with encephalitis, meningitis, myocarditis, acute flaccid paralysis, pulmonary edema in severe cases and can occasionally lead to death (1–4). HFMD/HA outbreaks occurred frequently worldwide, especially in the Asia–Pacific region. According to the most recent data released by Chinese Center for Disease Control and Prevention, ~2 million cases of HFMD were diagnosed in China, 2019, and 20 deaths were reported [quoted from (5)].
EVs belong to the family Picornaviridae, which comprises four human disease-associated species of Enterovirus A-D (6, 7). It has been reported that more than 20 EV types were capable to cause HFMD/HA (8). Of the EV types, enterovirus A71 (EV-A71) and coxsackievirus A16 (CVA16) has been recognized as the most predominant EV types associated with HFMD (9); while CVA2, CVA6, and CVA10 were the most prevalent types associated with HA (10, 11). Based on currently available data, it seems that CVA6 has replaced EV-A71 and CVA16 and became the most prevalent serotype associated with HFMD in regions of China (12–14). However, national epidemiological data on HFMD/HA associated EV types of recent 3 years remain scarce in China. We have observed that, between the year of 2012 and 2014, HFMD outbreaks were associated with EV other than EV-A71/CVA16, yet the subtype was not identified (15). Our limited knowledge on the epidemiology of EVs has hampered the development of effective diagnosis and prevention options. Therefore, the present study was conducted to illustrate the latest molecular epidemiology of enteroviruses causing HFMD/HA in Zunyi, China, and to compare the clinical characteristics between patients infected with the most common EV types.
We collected and reviewed the electronic records of all cases with laboratory-confirmed HFMD/HA in Zunyi Medical University Third Affiliated Hospital (The First People's Hospital of Zunyi) between March 1 and July 31, 2019. Patients who presenting with oral vesicular exanthema/ulcers on the anterior tonsillar pillars, soft palate, buccal mucosa, or uvula are preliminarily diagnosed as clinical cases of HA, and patients who presenting with oral vesicular exanthema/ulcers on the tongue or the buccal mucosa, and vesicular rashes over the palms, soles, buttocks or the trunk are preliminarily diagnosed as clinical cases of HFMD. Laboratory-confirmed cases of HFMD/HA were diagnosed if a clinical case was tested positive for enterovirus nucleic acids. Patients were classified as severe in case they had any neurological complications, or cardiopulmonary complications (16, 17).
Specimen, Data Collection, and Processing
Pharyngeal swabs were the only sample sources, and the aliquot of remaining clinical samples were stored at −80°C. The presence of enteroviruses and the types of EV-A71 and CVA16 were detected using real-time RT-PCR kits following the manufacturer's instructions (Da AN Gene Co. Ltd, Guangzhou, China). For identification of non-EV71/CVA16 types, 372 confirmed inpatient cases (216 HFMD cases and 156 HA cases) were randomly selected for further partial VP1 gene amplification and sequencing. The demographic information, clinical data, and laboratory findings of patients were retrospectively collected through the hospital health/laboratory information system.
Partial VP1 Gene Amplification and Enterovirus Typing
Serotypes other than EV-A71/CVA16 were identified with reverse transcription-nested PCR (nested RT-PCR) as previously described (18). Viral RNA was extracted from 140 μl clinical samples using the Viral RNA Kit (TINGEN, Beijing, China) according to the manufacturer's protocols. Once the viral RNA was obtained, cDNA was generated in a 20 μl volume reaction mixture immediately using the PrimeScript™ RT reagent kit (TAKARA, Japan) according to the manufacturer's instruction. The preparation of VP1 fragments includes two amplification processes. First, a 50 μl reaction mixture containing 10 μl cDNA was used in the first round PCR using primers 222 and 224, and the reaction was under these conditions: 94°C for 5 min; 40 cycles at 94°C for the 30 s, annealing at 42°C for 30 s, and extension at 72°C for 1 min; and a final extension step at 72°C for 5 min. Then, 3 μl amplification product from the first round PCR was added to a second PCR with the primers AN88 and AN89 under the amplification conditions as described for the first round PCR but with annealing temperature of 60°C. The amplified products were purified from 1.2% agarose gels and subjected to sequencing using ABI PRISM310 Genetic Analyzer. Finally, the Basic Local Alignment Search Tool (www.ncbi.nlm.nih.gov/blast) was used to classify the types of enteroviruses.
Sequence analysis of major serotypes of CVA6, CVA2, and CVA5 was performed using Molecular Evolutionary Genetics Analysis software (MEGA, version 6.0). Phylogenetic trees were built based on partial VP1 genes using the neighbor-joining method with MEGA 6.0 software. Validation of reconstructed evolution trees was supported statistically using the bootstrap with 1000 replicates.
Data were processed and analyzed using Microsoft Excel 2007 and IBM SPSS V.19.0 software (IBM Corp, Armonk, NY USA). All categorical variables were reported as frequency and percentage, and their differences among the groups were compared using the Pearson χ2 or Fisher's exact tests. Quantitative data were reported as median with interquartile range and were compared among the groups using the Kruskal–Wallis test. P-value < 0.05 was considered as statistically significant.
Overall Distribution of EVs
We have previously showed that a seasonal peak appears between March and July every year in Zunyi, China (15); we therefore conducted the epidemiological surveillance during the period. From 1 March to 31 July 2019, a total of 1,377 EVs associated HFMD/HA inpatient cases (823 HFMD and 554 HA cases) were enrolled. Of them, four were EV-A71 (0.3%, 4/1,377), 84 were CVA16 (6.1%, 84/1,377), and 1,289 cases (93.6%, 1,289/1,377) were associated with non-EV-A71/CVA16 infection (Figure 1).
Typing of Non-EV-A71/CVA16 EVs
To identify the latest epidemic pathogen of HFMD/HA, 372 samples identified as types other than EV-A71/CVA16 cases were randomly selected and subjected to typing by partial VP1 gene sequencing. A total of 14 EV types have been successfully identified in 273 samples (73.4%), including 166 HFMD cases and 107 HA cases. Collectively, CVA6 (63.4%, 173/273), CVA2 (11.7%, 32/273), CVA5 (9.2%, 25/273) and CVA4 (3.3%, 9/273) were the most prevalent EV types associated with HFMD/HA cases (Figure 2). For HFMD cases, the three most common types were CVA6 (80.7%, 134/166), CVA2 (5.4%, 9/166) and CVA5 (3.0%, 5/166); similarly, for HA cases, the three most prevalent serotypes were CVA6 (36.5%, 39/107), CVA2 (21.5%, 23/107) and CVA5 (18.7%, 20/107). The proportion of CVA2, CVA4, and CVA5 were higher in HA cases than in HFMD cases, whereas the proportion of CVA6 was higher in HFMD cases than in HA cases (Table 1).
Figure 2. Percentage of non-EV-A71/CVA16 serotypes in HFMD/HA patients. Left: HFMD (166 cases in total); Right: HA (107 cases in total).
Clinical Characteristics and Laboratory Results of Patients With Diverse Enterovirus Infections
To analyze the clinical features of EVs associated HFMD/HA, we included 84 patients who were infected with CVA16 as control at the same time. In all, clinical information and laboratory data were collected from 361 HFMD/HA cases. As shown in Table 2, there were 232 (64.3%) males and 129 (35.7%) females, with a male/female ratio of 1.8:1. The majority of patients were under 5 years old (n = 338, 93.6%), and about half of the patients were aged 1–3 years (54.6%, 197/361).
Clinical features of patients who were infected with CVA2, CVA5, CVA6 and CVA16 are shown in Table 3. Of the 361 cases, 60 (16.6%) were severe cases and 301 (83.4%) were mild HFMD/HA cases. The majority of the patients had initial clinical symptom of fever. Fever and fever higher than 39°C were significantly associated with enterovirus types. Compared to CVA16, CVA2, CVA5, and CVA6 were more likely to cause fever and high fever. In addition, compared with other EVs, a lower incidence of neurological complications was observed in CVA16 infected patients. In this study, all severe cases had neurological signs, including startled reaction (48.3%), twitching (31.7%), hand and foot trembling, and vomiting. Startled reaction was observed more frequently in patients who were infected CVA5 than CVA6. CVA5-infected cases more frequently developed severe HFMD/HA. Compared to CVA6, CVA5 and CVA2, CVA16 was less likely to cause severe HFMD/HA in this cohort. All patients were cured and discharged from the hospital.
Laboratory data showed the neutrophil percentage was significantly higher in CVA5-associated cases than in CVA6- and CVA16-associated cases, and the lymphocyte percentage was correspondingly significantly lower in CVA5-associated cases than in CVA6- and CVA16-associated cases (Table 4). Elevated serum C-reactive protein (>10 mg/L) was more frequently observed in patients infected with CVA2 and CVA5 than patients infected with CVA16 (Table 4). Moreover, levels of CRP were higher than in patients infected with CVA5 than CV-A6 (Figure 3).
Figure 3. Levels of C-reactive protein (CRP) in cases with CVA6, CVA2, CVA5, and CVA16 (*P < 0.05; **P < 0.01; ***P < 0.001).
Since CVA6, CVA2, and CVA5 were the main EV types circulating in Zunyi, we analyzed the evolutionary relatedness with their corresponding reference strains sequences available from GenBank. According to the phylogenetic relations, CVA6 evolutional dendrogram was separated into A–E. The nucleotide identities of CVA6 sequences in this study were 91–96.3%, and all CVA6 strains clustered in to subclade E. Phylogenetic analysis showed CVA6 circulating in Zunyi was in a close genetic relationship to the strains identified in Tianjin, Shanghai, Beijing, and Guangzhou, China (Figure 4). VP1 based evolutional dendrogram of CVA5 separated at the viruses phylogenetically into clades A-D. The present results showed that the nucleotide identities of CVA5 sequences in this study were 84.4–95.7%, and CVA5 strains belonged to subclade D, which were in a close genetic relationship to Chinese strains identified in other cities or provinces, such as Yunnan, Jiangxi, and Shengzhen, and they also showed high sequence identity to two strains from Japan and Australian (Figure 5). An evolutional dendrogram of CVA2 separated into A–E. In this study, the nucleotide identities of CVA2 sequences were 83.1–97.3%. Like CVA2 strains identified in other cities of China, CVA2 strains identified in Zunyi city segregated into subclade E (Figure 6).
Figure 4. Phylogenetic analysis of CVA6 based on partial VP1 gene sequence. Strains identified from HFMD are marked by circle, and those identified from HA are marked by triangle.
Figure 5. Phylogenetic analysis of CVA5 based on partial VP1 gene sequence. Strains identified from HFMD are marked by circle, and those identified from HA are marked by triangle.
Figure 6. Phylogenetic analysis of CVA2 based on partial VP1 gene sequence. Strains identified from HFMD are marked by circle, and those identified from HA are marked by triangle.
Here we provide epidemiologic information of EV associated HFMD/HA in Zunyi, 2019. Although EV-A71 and CVA16 have long time been the main EV types associated with HFMD (19), it seems that CVA6 has replaced the two types and became the most predominant EV type associated with HFMD/HA in some regions in recent years (14, 20, 21). Our work provides important evidence for the requirement of further EV surveillance, diagnostic reagent and vaccine development.
In this study, in addition to CVA6, CVA2, and CVA5 were the most prevalent types causing HFMD/HA. Although they have occurred recently in Europe in 2016 and in the United States in 2001, few studies reported CVA2 and CVA5 as the main pathogens responsible for HFMD compared with EV-A71, CVA6, CVA10, and CVA4 (22–24). Given the high prevalence of CVA6, CVA2, and CVA5 in patients with HFMD/HA, sensitive and reliable typing tests for them may be crucial for daily clinical diagnosis.
Unlike HFMD, the hazards of HA disease may be underestimated because it isn't a notifiable disease in China and few studies analyzed the etiology of HA. The present study showed that CVA6 (36.4%), CVA2 (21.5%), and CVA5 (18.7%) were the most three prevalent types associated with HA. CVA10, CVA4, CVA5, and CVA2 have aslo been reported as the major types associated with HA in other regions, including France, Hangzhou, Korea, and Thailand (8, 23, 25, 26). Together these data suggest that HA outbreaks caused by CVA2 and CVA5 should be monitored in future. It is notable that only a very limited number of HFMD/HA cases were caused by EV-A71 in this study. Previous findings showed the number of HFMD caused by EV-A71 was reduced in Sichuan, China, which was thought to result from the implementation of the EV-A71 vaccination (27). In Zunyi, it remains unclear how much the EV-A71 vaccination has affected the local prevalence of EVs, since the EV-A71 vaccines were only available in March 2019, and few children were vaccinated with the vaccines in this region.
While CVA16-associated case is generally thought to self-limited, a fraction of patients also develops severe case (28–30). In agreement with previous findings, the majority of CVA16-associated cases presented with mild symptom and only a small number of cases were severe in this study. Of note, when compared to cases infected with other EV types, the incidence of fever caused by CVA16 was relatively lower and levels of white blood cell (WBC) and CRP were also lower.
To our knowledge, clinical features of CVA5-associated HFMD/HA cases were seldom reported. In the present study, it was found that, compared to CVA16, CVA5 was associated with higher incidence of severe cases, and patients with CVA5 infection seemed to have a higher WBC number and serum CRP level. Although all cases have rashes, the rashes caused by CVA5 were usually confined to one or two sites. Further multicenter studies are needed to confirm this finding.
It has been reported that CVA6 is associated with higher incidence of severe HFMD cases than other EV types and patients infected with CVA6 commonly present with high fever, onychomadesis, pigmentation, and vesicular rashes (31). In addition, it has been shown that CVA6 was more likely to cause atypical clinical symptoms than other EV types in children and adults (32). Indeed, we also observed onychomadesis and atypical vesicular rashes in some cases associated with CVA6 infection, however, this information was missing in the majority of the patients, and their differences were not able to analyze in the present study.
Phylogenetic analysis of CVA6, CVA5, and CVA2 strains showed that strains identified in this region were close to the corresponding strains from other regions of China, which suggest that the prevalent strains in this region were widespread circulating in China. The strains identified from HFMD and HA cases did not form distinctive clusters, suggesting that even the same EV type can cause different clinical signs and symptoms. Further studies are needed to determine the underlying mechanisms involved in the differences in clinical symptoms between HFMD and HA caused by the same EV types.
The present study has some limitations. First, the cases of HFMD included in this study were from one city. Second, some EV types were not determined because the partial VP1 fragment could not be amplified using the primers and some sequences were not determined by sequencing using the degenerate primers. Third, EVs co-infection with other viruses or bacteria was not determined which may influence accuracy on the description of clinical features and laboratory results for a specific EV type.
The results of our study indicate that CVA6, and CVA2 and CVA5 were the dominant serotypes causing HFMD/HA outbreak in Zunyi, China. Future studies in different regions are required to confirm the results of this study. This finding is of great importance as it strongly indicates that the currently routine clinical practice of identifying EV-A71 and CVA16 is not enough and diagnostic kits for CVA6, CVA5 and CVA2 are required. In addition, novel EV vaccines need to be able to cover the emerging new EV types.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Sequencing data was downloaded from the GenBank database (URL) under the accession numbers U05876, AB777928, KP143074, KP143078, KP143073, KP143075, KP143076, JQ364886, AB090799, AB114106, AB114103, AY421764, AF081297, KF412903, AY919527, JN203517, AY919579, JN203516, KF734951, KJ848305, KP717645, KF647890, KJ774102, KT186273, AY421763, MF589268, MF589266, KP006020, MH716179, MT409180, LC481443, KC879500, KC879505, JF773168, JF773161, HQ728261, MT081376, MH111030, LC555996, KP006020, FR797989, AB162720, KC879493, KC879509, KC879510, JX139774, KC879541, JN034200, JX538054, JX538144, JF773151, JF773154, JN896767, JN896771, KM816569, KM816571, KC879523, JN168998, JN034233, KU560973, KU677983, KF780565, LC038050, JX473465, KP289357, KP006003, AB794105, LC013390, JQ968964, JX088582, JX569050, AY421760.
This study was in compliance with the Helsinki Declaration and was approved by the Human Research Ethics Committee of Zunyi Medical University Third Affiliated Hospital (The First People's Hospital of Zunyi). Informed consents were obtained from parents or guardians before sample collection.
WZ and KW: conceptualization and design. YA: writing. KW: funding acquisition, and supervision. YA, WZ, JW, JZ, MS, SY, CD, XL, DW, PT, XC, and HZ: experiments. YA, WZ, JZ, MS, SY, HZ, and JW: data collection and analysis. All authors contributed to the article and approved the submitted version.
This work was supported by the special fund for Science and Technology Cooperation between Guizhou Province and Zunyi City (2015-38), R&D project provided by The First People's Hospital of Zunyi (Yuankezi 2020-5), 15851 elite project provided by Zunyi City in 2016, Innovation Group Project provided by Education Department of Guizhou Province (QiankeheKYzi 019), and Cooperation Project provided by Zunyi city and The First People's Hospital of Zunyi (HZzi-2020-2).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
2. Cai K, Wang Y, Guo Z, Yu H, Li H, Zhang L, et al. Clinical characteristics and managements of severe hand, foot and mouth disease caused by enterovirus A71 and coxsackievirus A16 in Shanghai, China. BMC Infect Dis. (2019) 19:285. doi: 10.1186/s12879-019-3878-6
3. Hoang MTV, Nguyen TA, Tran TT, Vu TTH, Le NTN, Nguyen THN, et al. Clinical and aetiological study of hand, foot and mouth disease in southern Vietnam, 2013-2015: inpatients and outpatients. Int J Infect Dis. (2019) 80:1–9. doi: 10.1016/j.ijid.2018.12.004
4. Qiu J, Yan H, Cheng N, Lu X, Hu X, Liang L, et al. The clinical and epidemiological study of children with hand, foot, and mouth disease in hunan, China from 2013 to 2017. Sci Rep. (2019) 9:11662. doi: 10.1038/s41598-019-48259-1
5. Chinese Center for Disease Control and Prevention. Available online at: http://www.nhc.gov.cn/jkj/s3578/202004/b1519e1bc1a944fc8ec176db600f68d1.shtml (accessed April 2020).
6. Oberste MS, Maher K, Kilpatrick DR, Flemister MR, Brown BA, Pallansch MA. Typing of human enteroviruses by partial sequencing of VP1. J Clin Microbiol. (1999) 37:1288–93. doi: 10.1128/JCM.37.5.1288-1293.1999
8. Li W, Li C, Liu L, Liu X, Shang S, Mao H, et al. Molecular epidemiology of enterovirus from children with herpangina or hand, foot, and mouth disease in Hangzhou, 2016. Arch Virol. (2019) 164:2565–71. doi: 10.1007/s00705-019-04356-0
9. Hu YQ, Xie GC, Li DD, Pang LL, Xie J, Wang P, et al. Prevalence of coxsackievirus A6 and enterovirus 71 in hand, foot and mouth disease in Nanjing, China in 2013. Pediatr Infect Dis J. (2015) 34:951–7. doi: 10.1097/INF.0000000000000794
10. Zeng HR, Zheng HY, Liu L, Tan XH, Sun LM, Fang L, et al. Etiological characteristics of herpangina cases in Guangzhou city in Guangdong province, 2015. (in Chinese). Chin J Exp Clin Virol. (2017) 31:409–13. doi: 10.3760/cma.j.issn.1003-9279.2017.05.007
11. Li JM, LI L, Gao L, Dong XC, Liu XH, et al. Characteristics of herpangina epidemic and related environmental impact factors in Tianjin (in Chinese). J Environ Health. (2018) 35:326–9. doi: 10.16241/j.cnki.1001-5914.2018.04.012
12. Tan X, Li L, Zhang B, Jorba J, Su X, Ji T, et al. Molecular epidemiology of coxsackievirus A6 associated with outbreaks of hand, foot, and mouth disease in Tianjin, China, in 2013. Arch Virol. (2015) 160:1097–104. doi: 10.1007/s00705-015-2340-3
13. Weng Y, Chen W, He W, Huang M, Zhu Y, Yan Y. Serotyping and genetic characterization of hand, foot, and mouth disease (HFMD)-associated enteroviruses of No-EV71 and Non-CVA16 circulating in Fujian, China, 2011–2015. Med Sci Monit. (2017) 23:2508–18. doi: 10.12659/MSM.901364
15. Zhang W, Huang B, She C, Liu Y, Tong H, Wang F, et al. An epidemic analysis of hand, foot, and mouth disease in Zunyi, China between 2012 and 2014. Saudi Med J. (2015) 36:593–8. doi: 10.15537/smj.2015.5.10859
16. Subspecialty Group of Infectious Diseases t.S.o.P.C.M.A., and National Center for Quality Control of Infectious, D. [Expert consensus on diagnosis and treatment of herpangina. (version 2019)]. Zhonghua Er Ke Za Zhi. (2019) 57:177−80. doi: 10.3760/cma.j.issn.0578-1310.2019.03.004
17. Li XW, Ni X, Qian SY, Wang Q, Jiang RM, Xu WB, et al. Chinese guidelines for the diagnosis and treatment of hand, foot and mouth disease (2018 edition). World J Pediatr. (2018) 14:437–47. doi: 10.1007/s12519-018-0189-8
18. Nix WA, Oberste MS, Pallansch MA. Sensitive, seminested PCR amplification of VP1 sequences for direct identification of all enterovirus serotypes from original clinical specimens. J Clin Microbiol. (2006) 44:2698–704. doi: 10.1128/JCM.00542-06
19. Lu S, Ang LW, Tay J, Phoon MC, Hsu JP, Cutter J, et al. Seroepidemiology of coxsackievirus A6, coxsackievirus A16, and enterovirus 71 infections among children and adolescents in Singapore, 2008-2010. PLoS ONE. (2015) 10:e0127999. doi: 10.1371/journal.pone.0127999
21. Bian L, Wang Y, Yao X, Mao Q, Xu M, Liang Z. Coxsackievirus A6: a new emerging pathogen causing hand, foot and mouth disease outbreaks worldwide. Exp Rev Anti Infect Ther. (2015) 13:1061–71. doi: 10.1586/14787210.2015.1058156
22. Yang Q, Zhang Y, Yan D, Zhu S, Wang D, Ji T, et al. Two genotypes of coxsackievirus A2 associated with hand, foot, and mouth disease circulating in China since 2008. PLoS ONE. (2016) 11:e0169021. doi: 10.1371/journal.pone.0169021
23. Chansaenroj J, Auphimai C, Puenpa J, Mauleekoonphairoj J, Wanlapakorn N, Vuthitanachot V, et al. High prevalence of coxsackievirus A2 in children with herpangina in Thailand in 2015. Virus Dis. (2017) 28:111–4. doi: 10.1007/s13337-017-0366-8
24. Gao L, Zou G, Liao Q, Zhou Y, Liu F, Dai B, et al. Spectrum of enterovirus serotypes causing uncomplicated hand, foot, and mouth disease and enteroviral diagnostic yield of different clinical samples. Clin Infect Dis. (2018) 67:1729–35. doi: 10.1093/cid/ciy341
25. Mirand A, Henquell C, Archimbaud C, Ughetto S, Antona D, Bailly JL, et al. Outbreak of hand, foot and mouth disease/herpangina associated with coxsackievirus A6 and A10 infections in 2010, France: a large citywide, prospective observational study. Clin Microbiol Infect. (2012) 18:E110–8. doi: 10.1111/j.1469-0691.2012.03789.x
26. Park K, Lee B, Baek K, Cheon D, Yeo S, Park J, et al. Enteroviruses isolated from herpangina and hand-foot-and-mouth disease in Korean children. Virol J. (2012) 9:205. doi: 10.1186/1743-422X-9-205
27. Peng D, Ma Y, Liu Y, Lv Q, Yin F. Epidemiological and aetiological characteristics of hand, foot, and mouth disease in Sichuan Province, China, 2011-2017. Sci Rep. (2020) 10:6117. doi: 10.1038/s41598-020-63274-3
31. Yang X, Li Y, Zhang C, Zhan W, Xie J, Hu S, et al. Clinical features and phylogenetic analysis of severe hand-foot-and-mouth disease caused by Coxsackievirus A6. Infect Genet Evol. (2020) 77:104054. doi: 10.1016/j.meegid.2019.104054
32. Ramirez-Fort MK, Downing C, Doan HQ, Benoist F, Oberste MS, Khan F, et al. Coxsackievirus A6 associated hand, foot and mouth disease in adults: clinical presentation and review of the literature. J Clin Virol. (2014) 60:381–6. doi: 10.1016/j.jcv.2014.04.023
Keywords: molecular epidemiology, clinical features, enterovirus types, herpangina, hand foot and mouth disease
Citation: Ai Y, Zhang W, Wu J, Zhang J, Shen M, Yao S, Deng C, Li X, Wu D, Tian P, Cheng X, Zha H and Wu K (2021) Molecular Epidemiology and Clinical Features of Enteroviruses-Associated Hand, Foot, and Mouth Disease and Herpangina Outbreak in Zunyi, China, 2019. Front. Med. 8:656699. doi: 10.3389/fmed.2021.656699
Received: 21 January 2021; Accepted: 31 March 2021;
Published: 26 April 2021.
Edited by:Sara Manti, University of Catania, Italy
Reviewed by:Nadia Yehya Soudani, Lebanese University, Lebanon
Aiping Wu, Suzhou Institute of Systems Medicine (ISM), China
Copyright © 2021 Ai, Zhang, Wu, Zhang, Shen, Yao, Deng, Li, Wu, Tian, Cheng, Zha and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kaifeng Wu, firstname.lastname@example.org
†These authors have contributed equally to this work