# Molecular epidemiology of human enterovirus 71 at the origin of an epidemic of fatal hand, foot and mouth disease cases in Cambodia

Veasna Duong, Channa Mey, Marc Eloit, Huachen Zhu, Lucie Danet, Zhong Huang, Gang Zou, Arnaud Tarantola, Justine J. Cheval, Philippe Pérot, et al.

► To cite this version:

Veasna Duong, Channa Mey, Marc Eloit, Huachen Zhu, Lucie Danet, et al.. Molecular epidemiology of human enterovirus 71 at the origin of an epidemic of fatal hand, foot and mouth disease cases in Cambodia. Emerging microbes & infections, 2016, 5 (9), pp.e104. 10.1038/emi.2016.101. pasteur-01739359

HAL Id: pasteur-01739359

https://pasteur.hal.science/pasteur-01739359

Submitted on 21 Mar 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ORIGINAL ARTICLE

# Molecular epidemiology of human enterovirus 71 at the origin of an epidemic of fatal hand, foot and mouth disease cases in Cambodia

Veasna Duong<fnref n="1" />, Channa Mey<fnref n="1" />, Marc Eloit², Huachen Zhu³, Lucie Danet<fnref n="1" />, Zhong Huang⁴, Gang Zou⁴, Arnaud Tarantola<fnref n="1" />, Justine Cheval⁵, Philippe Perot², Denis Laurent⁶, Beat Richner⁶, Santy Ky⁶, Sothy Heng⁶, Sok Touch⁷, Ly Sovann⁷, Rogier van Doorn⁸, Thanh Tan Tran⁸, Jeremy J Farrar⁸, David E Wentworth⁹, Suman R Das⁹, Timothy B Stockwell⁹, Jean-Claude Manuguerra², Francis Delpeyroux²,¹⁰, Yi Guan³, Ralf Altmeyer⁴ and Philippe Buchy<fnref n="1" />,¹¹

Human enterovirus 71 (EV-A71) causes hand, foot and mouth disease (HFMD). EV-A71 circulates in many countries and has caused large epidemics, especially in the Asia-Pacific region, since 1997. In April 2012, an undiagnosed fatal disease with neurological involvement and respiratory distress occurred in young children admitted to the Kantha Bopha Children's Hospital in Phnom Penh, Cambodia. Most died within a day of hospital admission, causing public panic and international concern. In this study, we describe the enterovirus (EV) genotypes that were isolated during the outbreak in 2012 and the following year. From June 2012 to November 2013, 312 specimens were collected from hospitalized and ambulatory patients and tested by generic EV and specific EV-A71 reverse transcription PCR. EV-A71 was detected in 208 clinical specimens while other EVs were found in 32 patients. The VP1 gene and/or the complete genome were generated. Our phylogenetic sequencing analysis demonstrated that 80 EV-A71 strains belonged to the C4a subgenotype and 3 EV-A71 strains belonged to the B5 genotype. Furthermore, some lineages of EV-A71 were found to have appeared in Cambodia following separate introductions from neighboring countries. Nineteen EV A (CV-A6 and CV-A16), 9 EV B (EV-B83, CV-B3, CV-B2, CV-A9, E-31, E-2 and EV-B80) and 4 EV C (EV-C116, EV-C96, CV-A20 and Vaccine-related PV-3) strains were also detected. We found no molecular markers of disease severity. We report here that EV-A71 genotype C4 was the main etiological agent of a large outbreak of HFMD and particularly of severe forms associated with central nervous system infections. The role played by other EVs in the epidemic could not be clearly established.

*Emerging Microbes & Infections* (2016) **5**, e104; doi:10.1038/emi.2016.101; published online 21 September 2016

**Keywords:** Cambodia; epidemic; hand foot and mouth disease; human enterovirus 71; severe disease

## INTRODUCTION

Enteroviruses (EVs) belong to the genus *Enterovirus* within the family *Picornaviridae*, which are small non-enveloped, single-stranded positive-sense RNA viruses with a genome of 7500 bases.<fnref n="1" /> Human enteroviruses (HEV) are classified into three main species: EV-A, EV-B, EV-C, and other less common EV D-J species.<fnref n="1" />,² Human enterovirus 71 (EV-A71) belongs to the EV-A species, which along with other EVs can cause hand, foot and mouth disease (HFMD).³

EV-A71 was first isolated in California in 1969,⁴ although it was suggested that the virus may have circulated in the Netherlands as early as 1963.⁵ EV-A71 is present in most countries. Large-scale epidemics, however, have been observed in Asia-Pacific since 1997.

EV-A71 has become endemic in this region, with regular emergences of new genetic lineages.⁶

Sequencing of the P1 coding region for capsid proteins is well correlated with the previous definition of serotype. The highly variable sequences of the EV-A71 VP1 gene have been used to group EV-A71 isolates into four genotypes (A, B, C and D).⁷,⁸ EV strains of the same genotype have >92% nucleotide sequence identity, whereas isolates with a diversity of >15% are considered as different genotypes.⁹ The A group consists of viruses isolated in the United States of America in 1969, including the prototype BrCr strain.⁴ Genotype B has five subgenotypes and has been circulating in Asia, particularly in Malaysia and Singapore, since the 1970s.¹⁰ Subgenotype B5 was reported during

<footnote n="1">¹Pasteur Institute in Cambodia, Phnom Penh 12000, Cambodia; ²Pasteur Institute, Paris 75724, France; ³State Key Laboratory of Emerging Infectious Diseases, School of Public Health, The University of Hong Kong, Hong Kong, China; ⁴Institut Pasteur in Shanghai, Shanghai 200025, China; ⁵PathoQuest, Paris 75015, France; ⁶Kantha Bopha Hospital, Phnom Penh 12000, Cambodia; ⁷Ministry of Health, Phnom Penh 12000, Cambodia; ⁸Oxford University Clinical Research Unit, Hospital for Tropical Diseases, Ho Chi Minh P1Q5, Vietnam; ⁹J. Craig Venter Institute, Rockville, MD 292037, USA; ¹⁰National Institute for Health and Medical Research, INSERM U994, Paris 75000, France and ¹¹GlaxoSmithKline Vaccines R&D, Singapore 189720, Singapore</footnote>
recent outbreaks in Thailand,¹¹,¹² Taiwan¹³,¹⁴ and Vietnam.¹⁵ Geno-
type C also comprises five subgenotypes. C1 emerged in the
mid-1980s¹⁰ and replaced genotype B as the predominant genotype
in Europe and America.¹⁰ Subgenotype C2 emerged in 1995 in
Australia and the United States of America.¹⁰ The C3 subgenotype
emerged in Korea and Japan.¹⁶,¹⁷ C4 has circulated predominantly in
China since 2000 and is present in many other Asian countries.¹⁸⁻²¹
The C5 subgenotype was reported in Japan, Vietnam and Taiwan.²⁰,²²
Other less frequent EV-A71 genotypes (D–J) have been isolated in
India, Africa and Madagascar.²³

In recent years, numerous HFMD outbreaks were reported in Asian countries, including China, Vietnam, Singapore, South Korea and Thailand.²⁴ In April 2012, an undiagnosed fatal disease with neurological involvement and respiratory failure was detected at Kantha Bopha Children's Hospital in Phnom Penh, Cambodia. Most of the children who were affected were aged below three years and mainly originated from the southern and central parts of the country.²⁵ Most of them died (56/61 cases), usually within a day after admission to the hospital.²⁵ Etiological and epidemiological investigations were initiated in June 2012 after the first series of clinical samples were submitted to the Virology Unit at Institut Pasteur in Cambodia (IPC), Phnom Penh, Cambodia. EV-A71 followed by other non-EV-A71 EVs were the microorganisms that were most commonly detected in these patients.²⁶ The present study aims to describe the EV genotypes that were isolated during the outbreak in 2012 and 2013.

MATERIALS AND METHODS

Sample collection

From June 2012 to November 2013, clinical specimens such as throat, nasopharyngeal and rectal swabs, cerebrospinal fluid (CSF), bronch-oalveolar lavage and sera were collected from hospitalized patients and outpatients at Kantha Bopha Hospital in Phnom Penh and Jayavarman VII Hospital in Siem Reap (Northwest Cambodia).

Patients who were admitted to the two hospitals with signs and symptoms of EV infection (i.e., all patients hospitalized with severe disease as well as some patients randomly selected in the outpatient ward presenting with mild disease) were classified into the following three main clinical categories: HFMD, HFMD with central nervous system involvement (CNSI) alone and HFMD with cardiopulmonary failure (CPF). All patients in the CPF group also experienced CNSI.²⁷

The patient demographic data, clinical symptoms, diagnostic test
results and medical imaging findings were collected by the clinical
team. Because the initial series of samples used in this study were
collected for diagnostic purposes and characterization of the micro-
organism involved in an outbreak setting, approval was not sought
from the National Ethics Committee. Once the causative agent was
identified, patients with encephalitis were included in an ongoing
study on the etiology of pediatric encephalitis; this study was approved
by the National Ethics Committee (approval no. 107NECHR).
Additional authorization was obtained to include the outpatients
presenting with HFMD (approval no. 0059NECHR). For all patients,
a signed informed consent was obtained from their guardians before
inclusion.

EV-A71 diagnosis

Diagnostic reverse transcription PCR (RT-PCR). RNA was extracted from 140 µL of rectal, throat and nasopharyngeal swabs collected in viral transport medium, serum and CSF using the QIAamp Viral RNA Mini Kit (Qiagen1, Hilden, Germany) according to the manufacturer's recommendations.

A specific real-time RT-PCR assay was used for the detection of EV-A71 as described elsewhere.²⁸ Briefly, five µL of RNA was reverse-transcribed and amplified using the SuperScript III Platinum One-Step qRT-PCR Kit (Invitrogen, Carlsbad, CA, USA) in a mixture of 10 µL of 2x Reaction Mix with ROX, 10 µM forward primer (EV71-VP1-634F: 5'-GGA GAA CAC AAR CAR GAG AAA GA-3'), 10 µM reverse primer (EV71-VP1-743R: 5'-ACT AAA GGG TAC TTG GAC TTV GA-3'), 10 µM of EV-A71-specific probe (EV71-VP1-TaqMan: 5'-FAM-TGA TGG GCA CGT TCT CAG TGC G-BHQ1-3') and one µL of SuperScript III RT/Platinum Taq Mix. Molecular-grade water was added to the reaction to obtain a final volume of 20 µL. Thermocycling was conducted as follows: reverse transcription at 50 °C for 30 min, then holding at 95 °C for 2 min, followed by 40 cycles of 95 °C for 15 s, 55 °C for 30 s and 72 °C for 20 s.

The samples that tested negative for EV-A71 were further investi-
gated using pan-EV real-time RT-PCR adapted from Beld et al.,²⁹
using a few modifications for the detection of other HEVs. Specifically,
the RT step was performed prior to TaqMan real-time PCR. The
extracted RNAs were reverse-transcribed using random hexamers. The
RT was conducted by adding seven µL of RNA into a mixture
containing 1500 ng of random hexamers and 10 mM of each
deoxynucleoside triphosphate. Molecular-grade water was added to
obtain the final volume of 12 µL. After incubation for five min at 60 °C
and one min on ice, seven µL of the second mix (containing 4 µL of
5 × buffers, 0.1 M of DTT, 40 U of RNaseOut and five U of Super-
Script III Reverse Transcriptase (Invitrogen) and 1.4 µL of molecular-
grade water) was added to the first mixture. The second mixture was
incubated at 25 °C for 10 min, at 50 °C for one h and at 72 °C for
15 min.

Virus isolation. Patient's samples were inoculated onto Vero E6 cells. Briefly, 150 µL of CSF, 20 µL of sera or 100 µL of rectal, nasopharyngeal and throat swab samples collected in viral transport medium were added to 80% confluent Vero E6 cells in a 12-well plate. After incubating for 1 h at 37 °C with 5% CO₂, 1.7 mL of maintenance medium (Dulbecco's Modified Eagle Medium+2% of heat inactivated fetal calf serum and 1% of Penicillin–Streptomycin (Sigma-Aldrich, Saint Louis, MO, USA)) was added to each well. The plate was then incubated for seven days with the same incubation conditions. The presence of the virus in the supernatant was revealed by the diagnostic real-time PCR method described above.

Sequencing of EV-A71

Sequencing of the VP1 gene. RNA from the PCR-positive cell culture supernatants was extracted using the QIAamp Viral RNA Mini Kt (Qiagen, Hilden, Germany) and the VP1 gene was sequenced using the primers described by Nix et al.³⁰ DNA sequencing was performed in a commercial facility using an ABI 3730XL Analyzer (96-capillary) using the ABI Prism BigDye Terminator Cycle Sequencing Kit according to the manufacturer's instructions (Macrogen, Seoul, South Korea).

Whole-genome sequencing. Whole-genome sequencing of the high-
titrated cell culture supernatants was performed using high-
throughput sequencing. The culture supernatants of two isolates were
clarified by centrifugation at 3500–4000 r for 10 min at 4 °C to remove
the cell debris. They were then ultracentrifuged at 50 000 r.p.m. to
collect the viral pellet. The pellet was resuspended in 140 µL of 10 mM
Tris-Cl buffer (pH 7.5) and treated with RNase, followed by the
addition of 560 µL of prepared AVL Buffer (Qiagen) not containing
carrier RNAs to cleave the viral particles. RNA extraction was
conducted using the QIAamp Viral RNA Mini Kit (Qiagen) and DNAs were removed using the RNase-Free DNase Set (Qiagen). Double-stranded cDNAs were synthesized using the cDNA Synthesis System (Roche Applied Science, Penzberg, Germany) and random amplification of the viral cDNAs was then carried out as previously described.³¹ Whole-genomic sequencing was conducted on a Genome Sequencer Junior (Roche Applied Science), and the original reads were assembled into contigs using Lasergene, version 9.0 (http://www.dnastar.com).

For sample preparation, sequencing and bioinformatic analyses were conducted as described in Gagnieur et al.,³² except that the human reads were filtered against NCBI build 37.1/assembly hg19. Each sample was sequenced with a depth ranging from 67 to 233 million reads each.

*Whole-genome sequencing by Sanger method.* The total RNA was extracted from biological samples or cell culture suspension using the QIAamp Viral RNA Mini Kit (Qiagen) as described above. cDNA was synthesized from 3 µL of RNA using two µL of each nucleotide triphosphate (Eurobio, Mannheim, Germany), 10 pmol of each EV-A71 primer (Supplementary Table S1) and 7.5 µL of sterilized distilled water. The mixture was incubated at 70 °C for 5 min and then stored on ice for 2 min. Each tube was then incubated with 1 µL of SuperScript II Reverse Transcriptase (Invitrogen, Life Technologies) and 0.5 µL of RNasin (Promega, Madison, WI, USA) in a final volume of 20 µL for 50 min at 42 °C, then at 95 °C for 5 min and at 4 °C for 1 min.

5 µL of cDNA was then amplified by sequencing PCR using 10 pmol of paired primers for each EV-A71 fragment (Supplementary Table S1), one µL of AmpliTaq DNA polymerase (Applied Biosystems, Carlsbad, CA, USA; Life Technologies), 10 nmol of each nucleotide triphosphate and 10 × taq buffer for a total of 5 µL, in a final volume of 50 µL with H₂O added as necessary. The cycling conditions included initial denaturation at 94 °C for 3 min, 35 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, extension at 72 °C for 1 min and final extension at 72 °C for 5 min. The PCR amplicons were sequenced by the Sanger method at a commercial facility (Macrogen).

All partial VP1 and full-genome sequences generated in this study were submitted to the GenBank sequence database (NCBI, Bethesda, MD, USA) with accession numbers as follows: KX197414–KX197422 (EV-A71 VP1), KX197423–KX197454 (Non EV-A71, VP1) and KX197455–KX197465 (EV-A71 full genome) (Supplementary Table S2).

## Sequence analysis and phylogenetic analysis

Sequences generated from the PCR products obtained for each strain were analyzed and assembled using the CLC Main Workbench 5.5 package (CLC bio A/S, Aarhus, Denmark). The EV-A71 reference strains consist of 242 sequences of the VP1 gene and 457 full-genome sequences available in GenBank belonging to two genotypes: B1–B5 and C1–C5. For the phylogenetic analysis of the other EV isolates, 108 VP1 sequences of reference strains, including all four EV species (A, B, C and D), were used.

A Parechovirus (JX575746) sequence was used as an out-group in the phylogenetic analysis. Multiple sequence alignment of Cambodian strains with reference strains available in GenBank was conducted using Muscle,³³ which is available in the Seaview version 4.2.5 package.³⁴ Phylogenetic analyses were performed using the maximum likelihood method using the GTR model suggested by Jmodeltest,³⁵ which is available in MEGA 5.2³⁶ with 1000 bootstrap re-sampling.

## Recombination analysis

We used version 3 of the Recombinant Detection Program (RDP 4.56)³⁷ to screen for recombination in our set of EV-A71 full-genome sequence

data. This program implements several methods for the identification of recombinant sequences and recombination breakpoints. Accordingly, we used the RDP,³⁸ LARD,³⁹ Bootscanning,⁴⁰ Maxchi,⁴¹ Chimaera,⁴² GeneConv,⁴³ 3Seq⁴⁴ and Sis-scan⁴⁵ methods, using the default settings in each case. We only considered recombination events that were identified by at least three methods. Statistical significance was set to the *P* < 0.01 level after considering the Bonferroni correction for multiple comparisons as implemented in RDP.

## Statistical analyses

The means of continuous variables and the proportions of the categorical variables were compared using Student's *t*-test, Chi² test and Fisher's exact test as appropriate. Statistical significance was set to 0.05. All statistical analyses were carried out using the Stata 12 statistical software for Windows (StataCorp LP, College Station, TX, USA).

# RESULTS

## Patient characteristics and EV genotyping

Between June 2012 and November 2013, 142 cases of HFMD, 20 cases with CNSI and 150 cases of encephalitis with CPF were recruited by the two participating hospitals for inclusion in the present study (Table 1). Among the suspected cases, EV-A71 was detected by either RT-PCR and/or virus isolation in 96 (67.6%) of the HFMD patients, 11 (55.0%) patients with CNSI and 101 (67.3%) patients with CNSI associated with CPF. Additionally, co-infection with EV-A71 and non-EV-A71 was detected in five cases. This was immediately reported to the health authorities. The differences in EV-A71 prevalence by disease severity group were not statistically significant (*P* = 0.389).

The children included in this study were relatively young in all three clinical groups, with mean ages of 2.18 (95% confidence interval (CI): 1.78–2.58), 2.51 (95% CI: 1.04–4.03) and 2.84 (95% CI: 2.28–3.40) years in HFMD, CNSI and CPF cases, respectively (Table 1). There was no significant association between the age of the children and the severity of disease (*P* = 0.92). However, there were significantly more girls in CPF cases than in HFMD (*P* = 0.023) except in ‘pure’ encephalitic syndromes (*P* = 0.33).

The case fatality rate (CFR) was higher among CPF compared with CNSI cases (63.3% vs 15.0%, *P* < 0.001). No fatalities were reported in HFMD patients.

Laboratory-confirmed EV-A71 and other EV cases originated from 22 of the 24 Cambodian provinces. Most cases, however, came from the northwestern (78/240) and southern (98/240) parts of the country. This pattern corresponds both to the density of the population and to the catchment area of the two participating hospitals (Supplementary Figure S1).

Non-EV-A71 viruses were detected in 37 cases (11.9%). Among these, five were co-infections with EV-A71. Additionally, five EVs detected by generic RT-PCR were unsubtypeable (VP1 sequencing attempt failed) (Table 1). Of all the non-EV-A71 viruses, EV-A was the most frequent (59.5%) followed by EV-B (28.0%) and EV-C (12.5%). The mean age of non-EV-A71 cases was two years (n = 36; 95% CI: 1.54–2.77; range 0.25–7). The sex ratio was 0.76 males per female.

## Phylogenetic analysis

The phylogenetic analysis of the VP1 gene sequence of EV-A71 viruses was based on 80 sequences from Cambodia. This included nine VP1 sequences in addition to 11 VP1 sequences that were from fully sequenced strains and that were generated in the present study by either Sanger or high-throughput sequencing methods
<table><caption>Table 1 Patient characteristics and symptoms/clinical signs</caption><thead><tr><th></th><th>Total (n=312)</th><th>HFMD (n=142)</th><th>Central nervous system involvement (n=20)</th><th>Cardiopulmonary failure (n=150)</th></tr></thead><tbody><tr><td><strong>Age (years)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mean (95% CI)</td><td>2.51</td><td>2.18<sup>a</sup> (1.78-2.58)</td><td>2.51 (1.04-4.03)</td><td>2.84<sup>b</sup> (2.28-3.40)</td></tr><tr><td>Median</td><td>2</td><td>2</td><td>1.7</td><td>2</td></tr><tr><td>Min.-max.</td><td>0.16-27</td><td>0.75-27</td><td>0.91-14</td><td>0.16-27</td></tr><tr><td>Sex (% female)</td><td>134 (42.8%)</td><td>52 (36.6%)</td><td>8 (40.0%)</td><td>74<sup>c</sup> (50.3%)</td></tr><tr><td><strong>Hospitalization</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>No<sup>d</sup></td><td>62 (20.0%)</td><td>62 (43.4%)</td><td>0</td><td>0</td></tr><tr><td>Hospitalized</td><td>250 (80.0%)</td><td>80 (56.6%)</td><td>20 (100%)</td><td>150 (100%)</td></tr><tr><td><strong>Clinical symptoms</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Papulovesicular rash</td><td>130 (42.4%)</td><td>103 (72.0%)</td><td>5 (35.7%)</td><td>22 (14.7%)</td></tr><tr><td>Fatal cases</td><td>98 (31.4%)</td><td>0</td><td>3 (15.0%)</td><td>95 (63.3%)</td></tr><tr><td><strong>Laboratory diagnostic tests (RT-PCR and/or virus isolation)</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>EV-A71</td><td>208 (66.6%)</td><td>96 (67.6%)</td><td>11 (55.0%)</td><td>101 (67.3%)</td></tr><tr><td>Other than EV-A71</td><td>32 (10.3%)</td><td>18 (12.7%)</td><td>5 (25.0%)</td><td>9 (6.0%)</td></tr><tr><td>Subtyped</td><td>27 (73%)</td><td>15</td><td>5</td><td>7</td></tr><tr><td>Unsubtypable</td><td>5 (13.5%)</td><td>3</td><td>0</td><td>2</td></tr><tr><td>Co-infection with EV-A71 and non-EV-A71</td><td>5 (1.6%)</td><td>1 (0.7%)</td><td>0</td><td>4 (2.7%)</td></tr><tr><td>Negative</td><td>67 (21.5%)</td><td>27 (19.0%)</td><td>4 (20.0%)</td><td>36 (24.0%)</td></tr><tr><td><strong>Sequencing</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>EV-A71</td><td></td><td></td><td></td><td></td></tr><tr><td>VP1</td><td>9</td><td>3 (33.3%)</td><td>4 (44.5%)</td><td>2 (22.2%)</td></tr><tr><td>Full genome</td><td>11</td><td>4 (36.5%)</td><td>2 (18.0%)</td><td>5 (45.5%)</td></tr><tr><td>Other than EV-A71</td><td></td><td></td><td></td><td></td></tr><tr><td>VP1</td><td>32</td><td>16 (50.0%)</td><td>5 (15.6%)</td><td>11 (35.4%)</td></tr><tr><td>Enterovirus A</td><td>18 (56.3%)</td><td>14 CV-A6, 1 CV-A16</td><td></td><td>2 CV-A6, 1 CV-A16</td></tr><tr><td>Enterovirus B</td><td>7 (21.9%)</td><td>0</td><td>2 CV-B3, 1 E-2 and 1 EV-B80</td><td>1 EV-B83, 1 CV-B2, 1 E-31 and 1 CV-A9</td></tr><tr><td>Enterovirus C</td><td>2 (6.3%)</td><td>0</td><td>1 CV-A20</td><td>1 EV-C116</td></tr><tr><td>Co-infection with EV-A71</td><td>5 (15.5%)</td><td>1 CV-A6 (A)</td><td>0</td><td>2 CV-B2 (B), 1 EV-C96 (C) and 1 PV-3 (C)</td></tr></tbody></table>

Abbreviations: enterovirus A, (A); enterovirus B, (B); enterovirus C, (C); confidence interval, CI; coxsackievirus A, CVA; coxsackievirus B, CVB; echovirus, E; enterovirus, EV; human enterovirus 71, EV-A71; hand, foot and mouth disease, HFMD; poliovirus, P; reverse transcription PCR, RT-PCR; poliovirus E vaccine-related, VRPV.

(Table 1 and Supplementary Table S1). Additionally, 60 VP1 sequences were from full genomes of Cambodian strains generated in collaboration with the J. Craig Venter Institute and already available in GenBank. Phylogenetic analysis showed that all the strains that were isolated in Cambodia belonged to the C4a subgenotype except for three strains that clustered closely with the B5 subgenotype (Figure 1). A phylogenetic analysis of 70 full genomes of Cambodian EV-A71 strains was also conducted (Supplementary Figure S2). However, it did not demonstrate any significant difference compared with the VP1 gene analysis except for minor position changes, all within the same clusters. Consequently, only VP1 phylogeny will be discussed in this article. All the 77 Cambodian C4a strains clustered in three different lineages. Lineage 1 consisted of one strain that was detected in 2012 in Banteay Meanchey Province (Northwest region, close to the border with Thailand) that clustered closely to a Thai strain that was isolated in 2008. Five strains in lineage 2 were detected between 2012 and 2013 and were closely related to viruses originating from China between 2010 and 2012. Lineage 3 consisted of 76 Cambodian strains that were

detected between 2012 and 2013 from different provinces across Cambodia that clustered with viruses that were isolated in Vietnam in 2011 and 2012. The three Cambodian subgenotype B5 strains were detected in 2012 in Banteay Meanchey Province (northwest region, close to the border with Thailand). Of the three B5 strains, two clustered with strains that were isolated in Thailand in 2011. One strain clustered with both Vietnamese strains that were isolated during the 2011–2013 epidemics as well as with Thai strains that were detected in 2012 (Figure 1).

We did not observe any association between the precise phylogeny of the Cambodian strains and the severity of disease or the geographical distribution of cases (data not shown).

There was 97.8% nucleotide identity between the VP1 and/or full-genome sequences of the Cambodian strains of lineage 1 and Thai strain FJ556875 THA08. Nucleotide identity of lineage 2 and the sequences of Chinese viruses isolated from 2010 to 2012 ranged from 97.4% to 99.1%. Nucleotide identity of lineage 3 and the sequences of the Vietnamese viruses that were isolated in 2011 ranged from 97.4%

<footnote><sup>a</sup>n=141.</footnote>

<footnote><sup>b</sup>n=146.</footnote>

<footnote><sup>c</sup>n=147.</footnote>

<footnote><sup>d</sup>Patients recruited at the outpatient department.</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1** Maximum likelihood phylogenetic tree of 269 EV-A71 VP1 sequences. The branches are color-coded according to the location of sample collection (Cambodia = red, Vietnam = blue, Thailand = purple and China = green). For better clarity, taxon names are not shown, and sequences from other locations beside the four countries are not color-coded. The sequences isolated from Cambodian patients with different disease categories of severity are marked by white triangles for hand, foot and mouth disease, white rectangles for central nervous system involvement and black circles for cardiopulmonary failure. The tree was built using the maximum likelihood method based on the GTR+G4 model. The robustness of nodes was assessed with 1000 bootstrap replicates. Bootstrap values <70 are not shown.</figcaption></figure>

to 99.1%. The percentage of homology between the sequences of the Cambodian B5 subgenotype strains and the corresponding strains of Vietnamese and Thai origin with which they clustered ranged from 98.7% to 99.6% and from 98.7% to 99.0%, respectively.

Phylogenetic analysis of 32 non-EV-A71 strains demonstrated that the Cambodian sequences belonged to three EV species (Figure 2), including 19 strains belonging to EV-A (17 Coxsackievirus (CV) A6 and 2 CV-A16), nine strains belonging to EV-B (1 EV-B3, 2 CV-B3, 2 CV-B2, 1 CV-A9, 1 Echovirus (E)-31, 1 E-2 and 1 EV-B80) and four strains belonging to EV-C (one each of EV-C116, EV-C96, CV-A20 and poliovirus 3). The poliovirus 3 strain had a homology of >99% with a homologous Sabin vaccine strain and was classified as poliovirus vaccine-related. Of all non-EV-A71 viruses, CV-A6 was detected significantly more frequently in mild HFMD cases than in severe cases (*P*<0.001), while EV-B and EV-C were only found in the

most severe cases (*P*<0.001; Table 1). Phylogenetic analysis of the CV-A6 sequences displayed close genetic relations to CV-A6 strains that circulated in Thailand and China (Figure 2).

## Mutations and recombination

We did not observe any of the known mutations that are associated with a specific clinical group. A recombination event was detected by five methods (RDP, Genecov, MaxChi, Chimaera and SiScan) in the RDP software with *P*-value <0.001 in one Cambodian EV-A71 strain (KP308450) collected in 2012. The recombination analysis detected two breakpoints at nucleotide positions 390 and 3664. The major parent was HM622391 (detected in Taiwan in 2008) and the minor parent was EF373576 (Taiwan, 2004) from the same genotype. No other recombination events were detected in the 2012–2013 Cambodian strains.
<figure><img src="image_3.png" /><figcaption>**Figure 2** Maximum likelihood (ML) phylogenetic tree of 140 VP1 sequences from all four enterovirus species (A, B, C and D). The tree was built using the ML method based on the GTR+G4 model. The robustness of nodes was assessed with 1000 bootstrap replicates. Bootstrap values <70 are not shown. The sequences isolated from Cambodian patients with different disease categories of severity are marked by white triangles for hand, foot and mouth disease, white rectangles for central nervous system involvement and black circles for cardiopulmonary failure.</figcaption></figure>

DISCUSSION

This work describes the molecular characteristics of the EV-A71
viruses causing the first reported outbreak of EV infection with a high
apparent number of severe outcomes in Cambodia between 2012 and
2013. The outbreak of HFMD in early July 2012 was marked by a very
high CFR (92%) among patients at the intensive care units of Kantha
Bopha Children's Hospitals. This caused panic among the public,
which was abundantly reported by international media and commen-
ted in the scientific community.²⁵,⁴⁶,⁴⁷ This study demonstrates that
EV-A71 was the main etiological agent of the Cambodian outbreak.
The circulation of various non-EV-A71 EV types was also detected.
However, EVs in general are extremely common, especially in
children, and cause a spectrum of usually mild infectious diseases
that are rarely associated with epidemics. Their detection, essentially
by rectal swabs, was therefore most likely not directly related to the
outbreak of severe disease. Most patients originated from the North-
eastern and Southern regions of Cambodia in which the two large
pediatric hospitals participating in the study were located. These
regions may represent a bias and limitation of geographical recruit-
ment, especially as these hospitals provide care free of charge. Similar
to other reports, young children were the most affected age group in
our study with a median age of two years.²⁰,⁴⁸,⁴⁹ We also observed that
the proportion of girls was significantly higher in the clinically severe
group. In a meta-analysis of 19 separate studies, however, gender was
not a risk factor, whereas young age was a risk factor.⁵⁰ Indeed,
children aged below two years were often found to be at a higher risk

of developing severe complications.⁵¹⁻⁵³ We did not observe any
difference in median age between the three clinical groups in the
present study. The association of other factors with poor outcomes
will be assessed in another article.

The phylogenetic analysis of 80 Cambodian strains demonstrated that the viruses clustered in subgenotypes C4a and B5. Subgenotype C4a has been reported to cause large outbreaks in Asia, including in China in 2012,⁵⁴ Vietnam in 2011²⁸ and Thailand in 2008–2009.⁴⁸ Other outbreaks were reported with co-circulation of different genotypes or subgenotypes and other EVs of species such as in Vietnam in 2011 with subgenotypes EV-A71 C4 and C5,²⁸ Thailand with EV-A71, CV-A16 and CV-A10⁴⁸ and Singapore in which EV-A71 co-circulated with CV-A6, CV-A10, CV-A16 and CV-A4.⁴⁹ Although subgenotype B5 was detected in Cambodia in a province bordering Thailand in 2012, this genotype was not detected subsequently. Genotype B5 virus was possibly introduced into Cambodia from neighboring Thailand or Vietnam where this virus circulated previously and caused outbreaks of HFMD in 2012¹¹,¹²,⁵⁵ if not earlier, according to serological studies conducted in Thailand.¹² In 2012, Vietnam experienced a genotype switch from subgenogroup C4 to B5. Although the switch to B5 was associated with small outbreaks in Vietnam as in Sarawak (Malaysia) in 2003 and 2006,⁵⁶ Taiwan experienced a large outbreak of subgenotype B5 strains in 2008.⁵⁷

Interestingly, three distinct lineages were identified in subgenotype
C4. Lineage 1 grouped closely to a Thai strain, lineage 2 was closely
related to Chinese strains and lineage 3 clustered with Vietnamese
strains. Cambodian B5 strains clustered with strains that were detected previously or concurrently in Thailand and Vietnam. This observation suggests several separate introductions of EV-A71 viruses from neighboring countries into Cambodia, in view of the very high percentage of nucleotide identity (>97%) between the Cambodian strains and those isolated in China, Vietnam and Thailand. The circulation of viruses between neighboring countries is often observed with EVs and has been regularly documented in the past.⁵⁸⁻⁶¹ Several genotypes (e.g., C1, C4, C5, B5) of EV-A71 and some non-EV-A71 viruses (e.g., CV-A16) circulated in Vietnam since 2005 and in Thailand since 2008,¹²,⁴⁸ with C5 being the predominant EV-A71 subgenotype.²⁰,²⁸ The reason why HFMD was never detected or was detected but never reported in Cambodia before 2012 remains unclear. One of the best hypotheses is that the disease was misdiagnosed by clinicians. Indeed, our recent retrospective serological study using a EV-A71 strain isolated during the 2012 Cambodian outbreak to detect specific neutralizing antibodies showed that EV-A71 was widely circulating in the country at least during the preceding decade, with peaks in 2001, 2005, 2007 and 2010.⁶² These peaks coincided with outbreaks that occurred in Vietnam in 2005²⁸ and 2011²⁰ as well as in China and Taiwan in 2004–2005, 2008 and 2010–2012.¹⁴,⁵³,⁵⁴,⁶³

During the 2012 outbreak, many severe cases experiencing CPF did not present the typical rash usually observed in classic HFMD. Indeed, a papulovesicular rash was observed only in 14.7% of the CPF cases vs 72% in mild HFMD cases. Additionally, most admitted patients were critical and died within a few hours after admission. This patient outcome did not provide enough time for clinicians to conduct extensive diagnostic tests. The high CFR (92%) drew the attention of clinicians while many mild HFMD cases likely went unnoticed during the first outbreak period.

The high CFR was observed at the beginning of the outbreak in 2012 and then decreased slightly to an average of 63.3% among the patients presenting with CPF. The CFR in Cambodia was extremely high compared with that observed in other countries such as in China, Taiwan and Vietnam, where the CFR ranged from 0.14% to 34% depending on year, age group and gender.⁵²,⁶⁴,⁶⁵ Subgenotype C4a, which was detected in most Cambodian cases, was already suspected to both cause more severe disease (CPF) and be associated with higher lethality.⁵⁰,⁵²,⁶⁶⁻⁶⁸ We did not find any significant molecular difference, however, between EV-A71 subgenotype C4a sequences isolated from mild HFMD and severe (CNSI and CPF) cases. Similar observations were also reported by other authors.⁵²,⁶⁹ There may be other potential explanations for the severity of EV-A71 infections reported in Cambodia such as malnutrition, vitamin deficits, co-infection with other infectious diseases or delay in hospital admission combined with a lack of specific intensive care equipment (such as extracorporeal membrane oxygenation) or intravenous immunoglobulins.⁵⁰,⁷⁰⁻⁷³ Host genetics may also have a role. Data to explore these hypotheses were not collected during the outbreak in 2012 and 2013. Further studies to explore the possible role of these factors would facilitate understanding the unique situation that Cambodia experienced in 2012. Co-infections involving several EVs were quite common in HFMD as well as in severe EV-A71 infections.⁷¹,⁷⁴ Four out of the five co-infections detected in this study, including the poliovirus vaccine-related EV, were associated with severe disease.

Recombinations are relatively common in EVs.⁵⁵,⁷⁵,⁷⁶ A single case of recombination was detected in Cambodia; interestingly, the parent strain appeared to originate from Taiwan. Studies from Vietnam⁵⁵ and

Thailand⁷⁷ also reported a low number of recombination events in their EV-A71 strains.

In addition to EV-A71, several other EVs were detected in the various clinical forms of the disease, including CV-A6, which was the most common of these non-EV-A71 viruses. In our study, CV-A6 was only associated with mild disease, while EV-A71 was detected in similar proportions among the three clinical forms. Our CV-A6 findings are consistent with other reports from the region.⁴⁹,⁷⁸⁻⁸¹ CV-A6 progressively became the most common virus detected among mild HFMD cases, apparently replacing other non-EV-A71 viruses, including CV-A16 in China, Thailand, Malaysia and Singapore.¹¹,⁴⁹,⁸²⁻⁸⁷ CV-A6 seems now to be associated with many outbreaks worldwide.⁸⁸ Unfortunately, no recent data regarding the circulation of CV-A6 in Vietnam is available. EV-B and EV-C species were only detected in severe (CNSI and CPF) cases. We cannot, however, confirm the association between these species and disease severity because none of these viruses were detected in CSF except for one case with CV-B3. In all other cases, the viruses were detected in other clinical specimens, especially in rectal swabs (17/32).

In this study, we report the identification of the origin of a large outbreak of EV-A71 and particularly of severe forms of central nervous system infections associated with fatal CPF. The main causative agent was EV-A71 genotype C4. The role of non-EV-A71 EVs such as CV-A6 in HFMD outbreaks will be studied further for potential consideration as a component of new vaccines to prevent HFMD. Physician training, improved reporting from Cambodian hospitals and cross-border surveillance should be implemented to deliver timely adequate care and curb the outbreak. Further study of the risk factors for severe HFMD, including host genetics, nutritional status, co-infection and comorbidity, should be investigated in order to better understand the pathogenesis of this disease.

## ACKNOWLEDGEMENTS

We thank the Institut Pasteur in Shanghai (Ke Lan, Peijun Ren, Qianqian Zhu, Shanshan Xu, Jia Liu, Lili Wang, Qingwei Liu), Institut Pasteur in Cambodia (Vincent Deubel), Institut Pasteur in Paris (Anna-Maria Burgiere, Erika Muth, Marc Jouan, Magali Herrant, Jerome Salomon, Marie-Line Joffret and Romain Volle) and the J. Craig Venter Institute (Rebecca Halpin, Dana Busan, Wei Wang and Xudong Lin), as well as the patients and the health-care teams from Kantha Bopha Hospital. This study was financially supported by the Outbreak investigation funds of the International Division of the Institut Pasteur. The sequencing work performed by the J. Craig Venter Institute was funded in whole or part with federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services under contract number HHSN272200900007C. The opinions expressed in this article are the authors' own and do not reflect the view of the Centers for Disease Control and Prevention, the Department of Health and Human Services or the United States Government.

1. Pallansch M, Roos R. Enteroviruses: polioviruses, Coxsackieviruses, echoviruses, and newer enteroviruses. In: Knipe DM, Howley PM (eds). *Field Virology*, 5th edn. Vol. 1. Lippincott Williams & Wilkins: Philadelphia, USA. 2007, 840–893.

2. Adams MJ, Lefkowitz EJ, King AM *et al*. Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses. *Arch Virol* 2015; **160**: 1837–1850.

3. Cabrerizo M, Tarrago D, Munoz-Almagro C *et al*. Molecular epidemiology of enterovirus 71, coxsackievirus A16 and A6 associated with hand, foot and mouth disease in Spain. *Clin Microbiol Infect* 2013; **20**: 150–156.

4. Schmidt NJ, Lennette EH, Ho HH. An apparently new enterovirus isolated from patients with disease of the central nervous system. *J Infect Dis* 1974; **129**: 304–309.

5. van der Sanden S, Koopmans M, Uslu G, van der Avoort H. Dutch Working Group for Clinical V. Epidemiology of enterovirus 71 in the Netherlands, 1963 to 2008. *J Clin Microbiol* 2009; **47**: 2826–2833.

6. Yip CC, Lau SK, Woo PC, Yuen KY. Human enterovirus 71 epidemics: what's next? *Emerg Health Threats J* 2013; **6**: 19780.
7 Brown BA, Oberste MS, Alexander Jr JP, Kennett ML, Pallansch MA. Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998. *J Virol* 1999; **73**: 9969-9975.

8 Nasri D, Bouslama L, Pillet S et al. Basic rationale, current methods and future directions for molecular typing of human enterovirus. *Expert Rev Mol Diagn* 2007; **7**: 419–434.

9 McMinn PC. Recent advances in the molecular epidemiology and control of human enterovirus 71 infection. *Curr Opin Virol* 2012; **2**: 199–205.

10 Solomon T, Lewthwaite P, Perera D et al. Virology, epidemiology, pathogenesis, and control of enterovirus 71. *Lancet Infect Dis* 2010; **10**: 778–790.

11 Puenpa J, Mauleekoonphairoj J, Linsuwanon P et al. Prevalence and characterization of enterovirus infections among pediatric patients with hand foot mouth disease, herpangina and influenza like illness in Thailand, 2012. *PLoS One* 2014; **9**: e98888.

12 Linsuwanon P, Puenpa J, Huang SW et al. Epidemiology and seroepidemiology of human enterovirus 71 among Thai populations. *J Biomed Sci* 2014; **21**: 16.

13 Luo ST, Chiang PS, Chung WY et al. Reemergence of enterovirus 71 epidemic in northern taiwan, 2012. *PLoS One* 2015; **10**: e0116322.

14 Wu WH, Kuo TC, Lin YT et al. Molecular epidemiology of enterovirus 71 infection in the central region of Taiwan from 2002 to 2012. *PLoS One* 2013; **8**: e83711.

15 Tan, le V, Tuyen NT, Thanh TT et al. A generic assay for whole-genome amplification and deep sequencing of enterovirus A71. *J Virol Methods* 2015; **215-216**: 30–36.

16 Cardosa MJ, Perera D, Brown BA et al. Molecular epidemiology of human enterovirus 71 strains and recent outbreaks in the Asia-Pacific region: comparative analysis of the VP1 and VP4 genes. *Emerg Infect Dis* 2003; **9**: 461–468.

17 Iwai M, Masaki A, Hasegawa S et al. Genetic changes of coxsackievirus A16 and enterovirus 71 isolated from hand, foot, and mouth disease patients in Toyama, Japan between 1981 and 2007. *Jpn J Infect Dis* 2009; **62**: 254–259.

18 He SJ, Han JF, Ding XX, Wang YD, Qin CF. Characterization of enterovirus 71 and coxsackievirus A16 isolated in hand, foot, and mouth disease patients in Guang-dong, 2010. *Int J Infect Dis* 2013; **17**: e1025-e1030.

19 Puenpa J, Theamboonlers A, Korkong S et al. Molecular characterization and complete genome analysis of human enterovirus 71 and coxsackievirus A16 from children with hand, foot and mouth disease in Thailand during 2008-2011. *Arch Virol* 2011; **156**: 2007-2013.

20 Tu PV, Thao NT, Perera D et al. Epidemiologic and virologic investigation of hand, foot, and mouth disease, southern Vietnam, 2005. *Emerg Infect Dis* 2007; **13**: 1733-1741.

21 Nguyen VH, Sibounheuang B, Phommasone K et al. First isolation and genomic characterization of enterovirus A71 and coxsackievirus A16 from hand and mouth disease patients in the Lao PDR. *New Microbes New Infect* 2014; **2**: 170-172.

22 Thoa, le PK, Chiang PS, Khanh TH et al. Genetic and antigenic characterization of enterovirus 71 in Ho Chi Minh City, Vietnam, 2011. *PLoS One* 2013; **8**: e69895.

23 Bessaud M, Razafindratsimyresy R, Nougairede A et al. Molecular comparison and evolutionary analyses of VP1 nucleotide sequences of new African human enterovirus 71 isolates reveal a wide genetic diversity. *PLoS One* 2014; **9**: e90624.

24 Wu JS, Zhao N, Pan H et al. Patterns of polymorphism and divergence in the VP1 gene of enterovirus 71 circulating in the Asia-Pacific region between 1994 and 2013. *J Virol Methods* 2013; **193**: 713–728.

25 ProMED-mail. *Undiagnosed Illness, Fatal, Child-Cambodia* (03): EV71. 20120708.1193960. International Society for Infectious Diseases: Brookline, MA, USA; 2012. Available at http://www.promedmail.org.

26 Buchy P. *Human Enterovirus 71–Cambodia: Genotype C4*. 20120726.1216290. International Society for Infectious Diseases: Brookline, MA, USA; 2012. Available at http://www.promedmail.org.

27 World Health Organization. *A Guide to Clinical Management and Public Health Response for Hand, Foot and Mouth Disease (HFMD)*. WHO: Geneva, Switzerland; 2011. Available at http://www.wprowho.int/publications/docs/Guidancefortheclinical-managementofHFMD.pdf (accessed 2 January 2016).

28 Khanh TH, Sabanathan S, Thanh TT et al. Enterovirus 71-associated hand, foot, and mouth disease, Southern Vietnam, 2011. *Emerg Infect Dis* 2012; **18**: 2002–2005.

29 Beld M, Minnaar R, Weel J et al. Highly sensitive assay for detection of enterovirus in clinical specimens by reverse transcription-PCR with an armored RNA internal control. *J Clin Microbiol* 2004; **42**: 3059–3064.

30 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–2704.

31 Wang D, Coscoy L, Zylberberg M et al. Microarray-based detection and genotyping of viral pathogens. *Proc Natl Acad Sci USA* 2002; **99**: 15687–15692.

32 Gagnieur L, Cheval J, Gratigny M et al. Unbiased analysis by high throughput sequencing of the viral diversity in fetal bovine serum and trypsin used in cell culture. *Biologicals* 2014; **42**: 145-152.

33 Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. *BMC Bioinformatics* 2004; **5**: 113.

34 Gouy M, Guindon S, Gascuel O. SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. *Mol Biol Evol* 2010; **27**: 221–224.

35 Posada D. Selection of models of DNA evolution with *jModelTest*. *Methods Mol Biol* 2009; **537**: 93–112.

36 Tamura K, Peterson D, Peterson N et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. *Mol Biol Evol* 2011; **28**: 2731–2739.

37 Martin DP, Lemey P, Lott M et al. RDP3: a flexible and fast computer program for analyzing recombination. *Bioinformatics* 2010; **26**: 2462–2463.

38 Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. *Bioinformatics* 2000; **16**: 562–563.

39 Holmes EC, Worobey M, Rambaut A. Phylogenetic evidence for recombination in dengue virus. *Mol Biol Evol* 1999; **16**: 405–409.

40 Salminen MO, Carr JK, Burke DS, McCutchan FE. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. *AIDS Res Hum Retroviruses* 1995; **11**: 1423–1425.

41 Smith JM. Analyzing the mosaic structure of genes. *J Mol Evol* 1992; **34**: 126-129.

42 Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. *Proc Natl Acad Sci USA* 2001; **98**: 13757-13762.

43 Padidam M, Sawyer S, Fauquet CM. Possible emergence of new geminiviruses by frequent recombination. *Virology* 1999; **265**: 218-225.

44 Boni MF, Posada D, Feldman MW. An exact nonparametric method for inferring mosaic structure in sequence triplets. *Genetics* 2007; **176**: 1035-1047.

45 Gibbs MJ, Armstrong JS, Gibbs AJ. Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. *Bioinformatics* 2000; **16**: 573-582.

46 Dussart P, Baril L, Petit L et al. Clinical and virological study of dengue cases and the members of their households: the multinational DENFRAME Project. *PLoS Negl Trop Dis* 2012; **6**: e1482.

47 Seiff A. Cambodia unravels cause of mystery illness. *Lancet* 2012; **380**: 206.

48 Chatproedprai S, Theanboonlers A, Korkong S et al. Clinical and molecular characterization of hand-foot-and-mouth disease in Thailand, 2008-2009. *Jpn J Infect Dis* 2010; **63**: 229-233.

49 Wu Y, Yeo A, Phoon MC et al. The largest outbreak of hand; foot and mouth disease in Singapore in 2008: the role of enterovirus 71 and coxsackievirus A strains. *Int J Infect Dis* 2010; **14**: e1076-e1081.

50 Fang Y, Wang S, Zhang L et al. Risk factors of severe hand, foot and mouth disease: a meta-analysis. *Scand J Infect Dis* 2014; **46**: 515-522.

51 Xu Z, Shen H, Wang Z et al. The population-based health effect of hand, foot and mouth disease in children in Shanghai. *Pediatr Infect Dis J* 2014; **33**: 448-452.

52 Liu SL, Pan H, Liu P et al. Comparative epidemiology and virology of fatal and nonfatal cases of hand, foot and mouth disease in mainland China from 2008 to 2014. *Rev Med Virol* 2015; **25**: 115–128.

53 Xing W, Liao Q, Viboud C et al. Hand, foot, and mouth disease in China, 2008-12: an epidemiological study. *Lancet Infect Dis* 2014; **14**: 308-318.

54 Ni H, Yi B, Yin J et al. Epidemiological and etiological characteristics of hand, foot, and mouth disease in Ningbo, China, 2008-2011. *J Clin Virol* 2012; **54**: 342-348.

55 Geoghegan JL, Tan le V, Kuhnert D et al. Phylodynamics of enterovirus A71-associated hand, foot, and mouth disease in Viet Nam. *J Virol* 2015; **89**: 8871-8879.

56 Chu KB, Kasri AR. Hand foot and mouth disease due to enterovirus 71 in Malaysia. *Viro Sin* 2011; **26**: 221-228.

57 Huang SW, Hsu YW, Smith DJ et al. Reemergence of enterovirus 71 in 2008 in taiwan: dynamics of genetic and antigenic evolution from 1998 to 2008. *J Clin Microbiol* 2009; **47**: 3653-3662.

58 Mizuta K, Abiko C, Murata T et al. Frequent importation of enterovirus 71 from surrounding countries into the local community of Yamagata, Japan, between 1998 and 2003. *J Clin Microbiol* 2005; **43**: 6171-6175.

59 Wang HY, Tsao KC, Hsieh CH et al. Inferring nonneutral evolution from contrasting patterns of polymorphisms and divergences in different protein coding regions of enterovirus 71 circulating in Taiwan during 1998-2003. *BMC Evol Biol* 2010; **10**: 294.

60 Kung SH, Wang SF, Huang CW et al. Genetic and antigenic analyses of enterovirus 71 isolates in Taiwan during 1998-2005. *Clin Microbiol Infect* 2007; **13**: 782-787.

61 AbuBakar S, Sam IC, Yusof J et al. Enterovirus 71 outbreak, Brunei. *Emerg Infect Dis* 2009; **15**: 79-82.

62 Horwood PF, Andronico A, Tarantola A et al. Seroepidemiology of Human Enterovirus 71 Infection among Children, Cambodia. *Emerg Infect Dis* 2016; **22**: 92-95.

63 Yang F, Ren L, Xiong Z et al. Enterovirus 71 outbreak in the People's Republic of China in 2008. *J Clin Microbiol* 2009; **47**: 2351-2352.

64 Nguyen NT, Pham HV, Hoang CQ et al. Epidemiological and clinical characteristics of children who died from hand, foot and mouth disease in Vietnam, 2011. *BMC Infect Dis* 2014; **14**: 341.

65 Chen KT, Chang HL, Wang ST, Cheng YT, Yang JY. Epidemiologic features of hand-foot-mouth disease and herpangina caused by enterovirus 71 in Taiwan, 1998-2005. *Pediatrics* 2007; **120**: e244-e252.

66 Gao LD, Hu SX, Zhang H et al. Correlation analysis of EV71 detection and case severity in hand, foot, and mouth disease in the Hunan Province of China. *PLoS One* 2014; **9**: e100003.

67 Yang T, Xu G, Dong H, Ye M, He T. A case-control study of risk factors for severe hand-foot-mouth disease among children in Ningbo, China, 2010-2011. *Eur J Pediatr* 2012; **171**: 1359-1364.

68 Xu W, Liu CF, Yan L et al. Distribution of enteroviruses in hospitalized children with hand, foot and mouth disease and relationship between pathogens and nervous system complications. *Viro J* 2012; **9**: 8.

69 Shimizu H, Utama A, Yoshii K et al. Enterovirus 71 from fatal and nonfatal cases of hand, foot and mouth disease epidemics in Malaysia, Japan and Taiwan in 1997-1998. *Jpn J Infect Dis* 1999; **52**: 12-15.

70 Chen S, Yang Y, Yan X et al. Influence of vitamin A status on the antiviral immunity of children with hand, foot and mouth disease. *Clin Nutr* 2012; **31**: 543-548.

71 Han JF, Zhang Y, Hou PQ et al. Human enterovirus co-infection in severe HFMD patients in China. *J Clin Virol* 2014; **61**: 621-622.

72 Chen SM, Du JW, Jin YM et al. Risk factors for severe hand-foot-mouth disease in children in Hainan, China, 2011-2012. *Asia Pac J Public Health* 2015; **27**: 715-722.

73 Chea S, Cheng YB, Chokephaibulkit K et al. Workshop on use of intravenous immunoglobulin in hand, foot and mouth disease in Southeast Asia. *Emerg Infect Dis* 2015; **21**: doi: 10.3201/eid2101140992.
74 Liu LJ, Xu HM, Li XJ *et al*. Co-detection in the pathogenesis of severe hand-foot-mouth disease. *Arch Virol* 2012; **157**: 2219-2222.

75 Zhang Y, Tan X, Cui A *et al*. Complete genome analysis of the C4 subgenotype strains of enterovirus 71: predominant recombination C4 viruses persistently circulating in China for 14 years. *PLoS One* 2013; **8**: e56341.

76 Yip CC, Lau SK, Zhou B *et al*. Emergence of enterovirus 71 "double-recombinant" strains belonging to a novel genotype D originating from southern China: first evidence for combination of intratypic and intertypic recombination events in EV71. *Arch Virol* 2010; **155**: 1413-1424.

77 Mauleekoonphairoj J, Vongpunsawad S, Puenpa J, Korkong S, Poovorawan Y. Complete genome sequence analysis of enterovirus 71 isolated from children with hand, foot, and mouth disease in Thailand, 2012-2014. *Virus Genes* 2015; **51**: 290-293.

78 Xu M, Su L, Cao L *et al*. Enterovirus genotypes causing hand foot and mouth disease in Shanghai, China: a molecular epidemiological analysis. *BMC Infect Dis* 2013; **13**: 489.

79 Park SH, Choi SS, Oh SA *et al*. Detection and characterization of enterovirus associated with herpangina and hand, foot, and mouth disease in Seoul, Korea. *Clin Lab* 2011; **57**: 959-967.

80 Zeng H, Lu J, Zheng H *et al*. The epidemiological study of coxsackievirus A6 revealing hand, foot and mouth disease epidemic patterns in Guangdong, China. *Sci Rep* 2015; **5**: 10550.

81 Huang Y, Zhou Y, Lu H *et al*. Characterization of severe hand, foot, and mouth disease in Shenzhen, China, 2009-2013. *J Med Virol* 2015; **87**: 1471-1479.

82 He YQ, Chen L, Xu WB *et al*. Emergence, circulation, and spatiotemporal phylogenetic analysis of coxsackievirus a6- and coxsackievirus a10-associated hand, foot, and mouth disease infections from 2008 to 2012 in Shenzhen, China. *J Clin Microbiol* 2013; **51**: 3560-3566.

83 Di B, Zhang Y, Xie H *et al*. Circulation of Coxsackievirus A6 in hand-foot-mouth disease in Guangzhou, 2010-2012. *Virol J* 2014; **11**: 157.

84 Puenpa J, Chieochansin T, Linsuwanon P *et al*. Hand, foot, and mouth disease caused by coxsackievirus A6, Thailand, 2012. *Emerg Infect Dis* 2013; **19**: 641-643.

85 Chen W, Weng YW, He WX *et al*. [Molecular epidemiology of HFMD-associated pathogen coxsackievirus A6 in Fujian Province, 2011–2013]. *Bing Du Xue Bao* 2014; **30**: 624-629 (in Chinese).

86 Xu M, Su L, Cao L *et al*. Genotypes of the enterovirus causing hand foot and mouth disease in Shanghai, China, 2012–2013. *PLoS One* 2015; **10**: e0138514.

87 Ling BP, Jalilian FA, Harmal NS, Yubbu P, Sekawi Z. Detection and characterization of viruses causing hand, foot and mouth disease from children in Seri Kembangan, Malaysia. *Trop Biomed* 2014; **31**: 654-662.

88 Bian L, Wang Y, Yao X *et al*. Coxsackievirus A6: a new emerging pathogen causing hand, foot and mouth disease outbreaks worldwide. *Expert Rev Anti Infect Ther* 2015; **13**: 1061-1071.

<figure><img src="image_4.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2016

Supplementary Information for this article can be found on the *Emerging Microbes & Infections* website (http://www.nature.com/emi)