Virol JVirology Journal1743-422XBioMed Central2095902029756441743-422X-7-27710.1186/1743-422X-7-277ResearchGenetic diversity and C2-like subgenogroup strains of enterovirus 71, Taiwan, 2008HuangYuan-Pin1yphuang@cdc.gov.twLinTsuey-Li1tllin@cdc.gov.twHsuLi-Ching1liching@cdc.gov.twChenYu-Ju1yjchen@cdc.gov.twTsengYin-Hsin1hyork@cdc.gov.twHsuChiu-Chu1jenn@cdc.gov.twFanWen-Bin1fanndino@cdc.gov.twYangJyh-Yuan1jyyang@cdc.gov.twChangFeng-Yee1fychang@cdc.gov.twWuHo-Sheng12wuhs@cdc.gov.twResearch and Diagnostic Center, Centers for Disease Control, Department of Health, Taipei, Taiwan, R.O.CSchool of Medical Laboratory Science and Biotechnology, Taipei Medical University, Taipei, Taiwan, R.O.C201020102010727727718201020102010Copyright ©2010 Huang et al; licensee BioMed Central Ltd.2010Huang et al; licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.Background

Human enterovirus 71 (EV-71) is known of having caused numerous outbreaks of hand-foot-mouth disease, and other clinical manifestations globally. In 2008, 989 EV-71 strains were isolated in Taiwan.

Results

In this study, the genetic and antigenic properties of these strains were analyzed and the genetic diversity of EV-71 subgenogroups surfacing in Taiwan was depicted, which includes 3 previously reported subgenogroups of C5, B5, and C4, and one C2-like subgenogroup. Based on the phylogenetic analyses using their complete genome nucleotide sequences and neutralization tests, the C2-like subgenogroup forms a genetically distinct cluster from other subgenogroups, and the antisera show a maximum of 128-fold decrease of neutralization titer against this subgenogroup. In addition, the subgenogroup C4 isolates of 2008 were found quite similar genetically to the Chinese strains that caused outbreaks in recent years and thus they should be carefully watched.

Conclusions

Other than to be the first report describing the existence of C2-like subgenogroup of EV-71 in Taiwan, this article also foresees a potential of subgenogroup C4 outbreaks in Taiwan in the near future.

Background

Belonging to the genus Enterovirus of the family Picornaviridae, human enterovirus 71 (EV-71) is one of the most causative pathogens infecting humans and may cause outbreaks of hand-foot-mouth disease (HFMD), herpangina, and severe neurological symptoms, especially in young children [1]. There are over one hundred serotypes identified in the genus Enterovirus [2], which was originally classified into polioviruses, coxsackievirus A, coxsackievirus B, and echoviruses on the basis of differences in cell tropism, infectivity, antigenicity, and pathogenicity [1]. In recent years, the genus Enterovirus was re-classified into ten species, Human enterovirus A, Human enterovirus B, Human enterovirus C, Human enterovirus D, Simian enterovirus A, Bovine enterovirus, Porcine enterovirus B, Human rhinovirus A, Human rhinovirus B, and Human rhinovirus C based on the molecular characteristics. Former Coxsackievirus A2 (CV-A2), CV-A3, CV-A4, CV-A5, CV-A6, CV-A7, CV-A8, CV-A10, CV-A12, CV-A14, CV-A16, EV-71, EV-76, EV-89, EV-90, EV-91, EV-92, Simian enteroviruses SV19, SV43, SV46, and A13 are now members of Human enterovirus A [3-5].

The positive-stranded RNA genome of EV-71 possesses approximately 7,500 nucleotides and includes three genomic regions designated P1, P2, and P3. P1 region encodes four structural capsid proteins (VP4, VP2, VP3, and VP1), while P2 and P3 encodes seven nonstructural proteins (2A, 2B, 2C, 3A, 3B, 3C, and 3D). The nonstructural proteins are involved in polyprotein processing, and the capsid proteins, especially VP1, contain many neutralization antigenic sites and correspond to the virus serotyping [6]. In previous studies, the N-terminal portion of the VP1 capsid protein (composed of 297 amino acids) was likely to contain a major antigenic region and had important neutralizing antibody determinants [7,8]. But in another study, two synthetic peptides containing the C-terminal part of the VP1 protein (amino acid 163-177 and 208-222) were capable of eliciting neutralizing antibodies against EV-71 [9]. In addition, three regions on the VP1 protein (amino acid 66-77, 145-159, and 247-261) were identified to be capable of inducing human EV-71-specific CD4+ T-cell proliferation [10]. However, the accurate locations of neutralizing epitopes are still uncertain. Recombination found in the same serotype (intratypic) or in the different serotype (intertypic) and point mutation events result in the evolution of EV. Multiple strains circulating at the same area may increase the possibility of recombination, and many recombinants have been observed in EV [11-13].

EV-71 is genetically divided into three genogroups, A, B, and C, on the basis of the VP1 sequences analyses [14]. Genogroups B and C are each further divided into five subgenogroups, designated as B1-B5 and C1-C5, while genogroup A contains only one strain, the prototype strain BrCr [15,16]. In addition, some uncommon subgenogroups were also identified. For instance, isolates of subgenogroups B0 were first observed in The Netherlands in 1963 [17], and those of subgenogroup C0 were observed in Japan in 1978 [18,19]. One Indian isolate in 2001 was genetically distinct from all other EV-71 strains and designated as genotype D [20].

Since EV-71 was first isolated in California in 1969, many EV-71 outbreaks have been reported worldwide, for instance, several outbreaks took place in the USA, Japan, and other countries in the 1970s (subgenogroup B1), in Hong Kong, Australia, and the USA in the 1980s (subgenogroups B1, B2, and C1), and especially in the Asian Pacific region in recent years [21,22]. Subgenogroup B3 was described in Sarawak, Singapore, and Australia in 1997, 1998, and 1999, respectively, while subgenogroup C4 was identified on Mainland China in 1998. After that, EV-71 epidemics of subgenogroup B4 were reported in Singapore, Sarawak, and Sydney, and those of subgenogroup C3 were described in Korea in 2000 [15]. Subgenogroup B5 was identified in Sarawak, Japan, and Singapore in the last decade and subgenogroup C5 in southern Vietnam in 2005 [16]. Since one subgenogroup could be found from different countries in the same or different period, to predict the epidemiological pattern of EV-71 infections is not easy. For example, subgenogroup C1 was first described in the United States in 1986 [14], but caused several outbreaks in Germany, Australia, the United Kingdom and other countries [23-25]. On the other hand, one subgenogroup could be identified in the same area during a long period; for instance, subgenogroup C4 showed up repeated on Mainland China from 1998 to 2008 [26].

In Taiwan, a large outbreak was reported in 1998, followed by two lesser outbreaks in 2000 and 2001, and one more in 2008 [27-29]. Based on a study covering 8-years, the incidence of mild cases of HFMD/herpangina was reported as 0.8 to 19.9 cases per sentinel physician per week. The seasonal incidence varied, but usually peaked in the summer [30]. Over the past several years, co-circulation patterns of various genetic subgenogroups were frequently observed in Taiwan. Back in 1998 for instance, the subgenogroup C2 was found to be the major one with subgenogroups B4 and C4 as two minors. Afterwards, the subgenogroup B4 was singled out as the major cause of the outbreaks with C4 as a minor in 2002, and then subgenogroup B5 became the major one with a minor C5 from 2006 to 2008 [21,31]. In such a situation, it is expected that the possibility of recombination between various subgenogroups of EV-71 increases. Therefore, we persistently analyzed all EV-71 isolates collected by our surveillance system, and tried to find out if any isolates were genetically distinct from those EV-71 strains isolated from earlier outbreaks by phylogenetic analyses and neutralization tests.

ResultsEpidemiological results

According to our laboratory surveillance data, EV-71 viruses of various subgenogroups were isolated from 989 patients in Taiwan in 2008. They were 413 females, 564 males, and 12 with gender not specified, and no significant differences were observed in gender distribution (p > 0.05). Among these patients with age ranging from 1 week to 38 years old, most (810/989, 81.9%) were under 5, including 342 girls, 460 boys and 8 with missing data of gender. EV-71 infections were reported throughout the year with a peak in the summer, roughly between May and July.

Basic Local Alignment Search Tool (BLAST) result

Four subgenotypes of EV-71, including 980 subgenogroup B5 isolates, 6 subgenogroup C4 isolates, 1 subgenogroup C5 isolate, and 2 subgenogroup C2-like isolates, were identified according to the BLAST results of partial VP1 region nucleotide sequences (Figure 1). All isolates showed extremely high identities with their respective reference strains (>97%), except the two C2-like isolates (<93%). The genotyping of the subgenogroup C2-like isolates were thus further confirmed by phylogenetic analysis. These two isolates, 2008-07776 and 2008-00643, were collected in Taipei County in May and August, respectively.

Different subgenogroups of 989 enterovirus 71 (EV-71) isolates in Taiwan in 2008 according to the BLAST results. The subgenogroup was determined by BLAST analysis of partial VP1 region nucleotide sequences. There were 980 subgenogroup B5 isolates, 6 subgenogroup C4 isolates, 1 subgenogroup C5 isolate, and 2 subgenogroup C2-like isolates identified according to the BLAST analysis.

Phylogenetic analysis and recombination analysis

After the BLAST process, four subgenogroup B5 and four subgenogroup C4 isolates randomly chosen, along with the only one subgenogroup C5, and two subgenogroup C2-like isolates, were used in a phylogenetic analysis on partial VP1 gene nucleotide sequence (Figure 2). The B5 and C5 isolates turned out to be genetically similar to the Taiwan strains isolated in 2007, while the C4 isolates tested were close to those China strains isolated in 2008-2009. Besides, the C2-like isolates were located in genogroup C, but not within any known subgenogroup.

Phylogenetic analysis of enterovirus 71 strains based on partial VP1 gene sequence (nucleotide position 16-418). Phylogenetic analysis was performed based on partial VP1 gene nucleotide sequences of reference strains from the GenBank and 11 representative isolates chosen from 989 sequenced isolates from Taiwan in 2008. The phylogenetic tree was constructed by the neighbor-joining method with MEGA version 4 software, and the reliabilities indicated at the branch nodes were evaluated using 1,000 bootstrap replications. Only values of over 70% were shown. The prototype coxsackievirus A16 (CA16) G-10 strain was used as an out-group.

Due to the uncertain genotyping on partial VP1 gene, with no more than 93% of the nucleotide identity between the C2-like isolates and the reference strains, each gene region of these two isolates was further sequenced and recombination analyses conducted. The nucleotide and amino acid identities between EV-71 subgenogroups were presented in Table 1. No amino acid changes were observed for C2-like isolates in the two regions of VP1 protein which were capable of eliciting neutralizing antibodies (amino acids 163-177 and 208-222). Moreover, there were no unique changes in three regions of VP1 protein, which were capable of inducing human EV-71-specific CD4+ T-cell proliferation (amino acids 66-77, 145-159, and 247-261). The phylogenetic analysis results showed that these 2 subgenogroup C2-like isolates formed a distinct cluster within genogroup C based on P1 and P2 region nucleotide sequences (Figure 3, panels A-B), and within genogroup B based on P3 region nucleotide sequences (Figure 3, panel C). The phylogenetic trees of each gene sequences were shown in Additional File 1.

Percent identity (%) of nucleotide and amino acid sequences in different gene fragment between subgenogroup C2-like and other subgenogroups of enterovirus 71*

Sub- genogroupGene

5'-UTRVP4VP2VP3VP12A2B2C3A3B3C3DComplete
Ant83.184.080.8-81.183.0-83.383.879.176.078.377.9-78.272.776.576.879.3-79.4
aa10097.697.994.996.692.996.697.690.993.491.795.3
B1nt84.2-84.478.784.5-84.679.6-79.883.5-83.776.8-78.072.3-73.082.4-82.581.3-81.781.879.5-80.577.7-78.380.2-80.6
aa98.5-10097.697.596.9-97.394.0-94.691.9-92.995.1-95.496.595.496.793.2-94.395.8-95.9
B2nt85.1-85.280.184.180.3-80.583.3-83.679.173.482.678.6-79.080.378.578.3-78.480.3-80.5
aa10097.697.197.394.692.994.597.695.496.793.9-94.195.8
B3nt84.1-84.682.6-83.082.4-82.579.2-79.682.9-83.378.8-79.176.783.4-83.683.7-84.490.985.4-85.684.7-85.282.4-82.6
aa10097.6-98.096.696.995.394.997.5-97.898.895.498.396.5-97.497.0-97.3
B4nt83.9-84.583.083.3-83.480.3-80.483.3-83.680.675.082.3-82.679.0-79.883.378.6-78.878.2-78.580.4-80.6
aa10098.097.197.695.392.996.395.3-96.510095.694.5-94.896.2-96.3
B5nt79.8-83.482.6-83.581.7-82.880.9-82.383.6-84.580.6-81.773.0-73.782.1-82.879.4-81.084.8-86.377.7-78.177.5-78.180.3-80.5
aa10097.2-98.097.197.9-98.394.0-95.492.995.7-96.397.695.4-10095.0-96.194.1-94.596.1-96.3
C1nt81.6-82.288.488.5-88.988.2-89.688.4-89.085.5-86.284.8-85.179.7-79.975.9-77.177.2-78.775.2-75.579.4-79.982.4-83.0
aa98.5-10099.2-99.699.1-99.599.696.6-97.394.996.6-96.991.895.493.4-93.994.3-94.596.7-96.9
C2nt81.2-83.391.7-94.693.7-94.493.8-95.191.9-93.090.8-92.690.9-91.578.7-80.477.1-79.075.7-77.273.2-75.078.9-80.084.0-85.4
aa10099.2-10010098.6-99.698.691.9-93.991.1-96.987.2-93.095.487.4-93.993.2-94.395.1-96.9
C3nt83.2-83.587.990.4-90.590.4-91.088.1-88.686.8-87.188.578.8-79.077.1-77.975.775.4-75.778.5-78.982.9-83.2
aa10099.610099.696.6-97.393.997.293.095.492.3-92.893.9-94.196.8
C4nt83.1-84.188.4-90.888.7-90.689.2-90.486.1-87.982.6-84.073.7-75.482.4-83.580.2-81.383.3-86.383.0-83.781.8-83.683.8-84.4
aa97.1-10099.6-10098.7-99.598.6-99.695.3-96.092.9-94.997.2-98.497.6-98.895.496.1-97.295.4-96.597.4-97.9
C5nt82.5-82.989.3-90.388.0-88.787.1-87.487.886.2-86.882.4-83.178.3-79.374.0-74.880.3-81.876.3-76.577.3-77.481.4-82.5
aa10099.610098.6-99.696.091.9-93.997.293.0-94.195.493.4-93.993.2-93.996.6-96.7

*Subgenogroup A: BrCr-CA-70 (GenBank accession no. U22521), B1: 236-TW86 (FJ357379) and 244-TW86 (FJ357381), B2: 26M/AUS/4/99 (EU364841), B3: SAR/SHA66 (AM396586) and 26M/AUS/4/99 (EU364841), B4: 5865/SIN/000009 (AF316321) and 5666/SIN/002209 (AF352027), B5: S19841-SAR-03 (DQ341363), 2007-08747 (EU527985) and 2009-03531 (HM622390), C1: 804/NO/03 (DQ452074) and 1M-AUS-12-00 (DQ341361), C2: 1245a/98/TW (AF176044), ENT/PM/SHA71 (AM396585), Tainan/5746/98 (AF304457), TW/2086/98 (AF119796) and 6F/AUS/6/99 (DQ381846), C3: 06/KOR/00 (DQ341355) and 03/KOR/00 (DQ341355), C4: 1235/04/TW (DQ133459), BJ08/Z004/3 (FJ606447), 1/SHENZHEN/08 (FJ607334), Shanghai/036/2009 (FJ713137) and SHZH03-CHN (AY465356), C5: E2005125-TW (EF063152) and 2007-07364 (EU527983), C2-like: 2008-00643 (HM622391) and 2008-07776 (HM622392).

Phylogenetic analyses of enterovirus 71 strains. The phylogenetic tree was constructed by the neighbor-joining method with MEGA version 4 software, and the reliabilities indicated at the branch nodes were evaluated using 1,000 bootstrap replications. Only values of over 70% were shown. The prototype coxsackievirus A16 (CA16) G-10 strain was used as an out-group. The tree was drawn on the basis of the P1 region nucleotide sequences (A), the P2 region nucleotide sequences (B), and the P3 region nucleotide sequences (C).

One suspected recombination event was shown in the similarity plot and bootscan analyses between subgenogroup C2 and subgenogroup B3 of EV-71 (p < 0.01) (Figure 4).

Bootscan analyses of enterovirus 71 nucleotide sequences. The subgenogroup C2-like strain 2008-00643 was queried against other subgenogroups of enterovirus 71 using SimPlot, version 3.5.1, in a sliding window of 400 nucleotides with a 20 nucleotides step.

Preparation of anti-enterovirus rabbit serum, and neutralization test

Anti-EV-71 rabbit sera against three subgenogroups (C2, C5, and B5) of EV-71 virus, with 100 cell culture infective dose (CCID50) viruses per 50 μl for immunization, were used for neutralization test. Table 2 shows the neutralization antibody titers against different subgenogroups of EV-71. Based on the data against their homo-subgenogroup viruses, antisera C2, C5, and B5 showed a 2- to 16-fold decrease in titers against their hetero-subgenogroups. However, the result of neutralization antibody titers of the same antisera against the C2-like subgenogroup showed an obvious difference (p < 0.05), with an 8- to 128-fold decrease compared to those of their homo-subgenogroup.

Neutralization antibody titers of rabbits antisera against enterovirus 71 (EV-71) from different subgenogroups

Antisera no.Subgenogroup of immunogenEV-71 strain

97111207 (C2) *E2004104 (C4) *E2006125 (C5) *E2002042 (B4) *E2007599 (B5) *C2-like
18,1921,0241,0244,0964,096256
2C232,7688,1924,096131,072131,0722,048
332,76832,76816,38465,53665,5364,096

432,768262,144131,07265,53665,5364,096
532,768524,288131,072262,144262,1444,096
6C565,536524,288262,144262,144262,1448,192

78,192262,14432,76832,76832,7682,048
816,3848,19232,768262,144131,0722,048
9B54,0964,0964,09632,7688,192128
1032,768131,072131,072262,144131,0721,024

* Statistically significant difference in log10-transformed data when compared to subgenogroup C2-like group (p < 0.05).

In addition, there were 11 pairs of serum samples used for neutralization test in this study, including acute-phase serum (3-7 days post infection) and recovery-phase serum (15-39 days post infection) (Table 3). Sera obtained from the patients with EV-71 infection belonging to subgenogroups B4, C4, C5, and B5 showed a maximum of 16-fold decrease in neutralization titers against hetero-subgenogroups of EV-71 as compared to the ones against their homo-subgenogroup. On the contrary, sera showed a maximum of 128-fold decrease against the C2-like subgenogroup. Taken together, these results indicated a divergence of antigenic relationship between the subgenogroup C2-like and other subgenogroups.

Serum neutralization antibody titers against different subgenogroups of enterovirus 71 (EV-71)

Antisera no.Subgenogroup of EV-71 infectionSampling period (days post infection)*EV-71 strain

97111207 (C2) *E2004104 (C4) *E2006125 (C5) *E2002042 (B4) *E2007599 (B5) *C2-like
1B4AP(6)1,0241,0241,0241,0242,04816
RP(32)1,0241,0241,0241,0242,04864
2C4AP(5)641,0242561,02451232
RP(13)5121,0241,0241,0244,09664
3C5AP(5)1,0242,0482,0481,0241,02464
RP(16)1,0241,0244,0962,048512256
4C5AP(4)2565121,0242,0482,04816
RP(16)2561,0241,0241,0241,02432
5C4AP(4)1281,0245122,0482,0488
RP(17)5124,0962,04816,38416,384128
6B5AP(6)2565125125121,02432
RP(27)2,0484,0962,0484,0964,096256
7B5AP(7)6432326464<8
RP(22)1281282565121,02464
8B5AP(3)8321664128<8
RP(15)1,0248,1924,0964,0964,09664
9B5AP(4)326464643216
RP(39)12812825651225664
10B5AP(4)1,0245125125121,02464
RP(23)5122562561,02451264
11B5AP(4)1285122562561,02432
RP(16)2,0488,1924,0964,0968,192256

* Statistically significant difference in log10-transformed data when compared to subgenogroup C2-like group (p < 0.05). AP: acute phase; RP: recovery phase.

Discussion

Enterovirus infections, especially EV-71, were associated with HFMD, herpangina, and neurological diseases and very common in the West Pacific region where Taiwan locates. There has been about two thousands isolates in Taiwan reported by the surveillance program each year since 2001 [16,32]. Moreover, emergence of new EV-71 subgenogroups was reported continuously. Because the VP1 gene is highly related to host neutralization antibodies and viral virulence, determining the genogroup of EV-71 is generally based on the VP1 gene sequence [17], and three genotypes were recognized accordingly [14]. A combination of VP1 and 3D gene sequences was proposed to be used for initial genotyping [19]. However, only a few studies about the antigenic variances of EV-71 have been reported [29,33].

In this study, we reported a genetic and antigenic diversity of EV-71 subgenogroups in Taiwan in 2008, including 3 previously reported subgenogroup C5, B5, C4, and one C2-like subgenogroup. The surveillance results of EV-71 molecular epidemiology in Taiwan was quite different from those in other counties, for example, genogroup C was the only one spotted in the United Kingdom from 1998 to 2006 and in Germany from 1997 to 2007 [24,25].

EV-71 of subgenogroup C5 was first isolated in southern Vietnam in 2005 and caused an outbreak with neurological disease and high prevalence [16]. According to our surveillance data in 2008, the isolate of subgenogroup C5 was identified in July, and this subgenogroup was still recognized in 2009 (unpublished data). Although these subgenogroup C5 strains were in low numbers and did not result in outbreaks in Taiwan in recent years [21], a previous report of EV-71 showed that the genogroups which caused outbreaks were usually in circulation 2 to 5 years before the onset of the outbreaks [29]. Hence we could not exclude the possibility of an outbreak caused by subgenogroup C5 strains in the subsequent years. Subgenogroup B5 strains were isolated in Taiwan in 2003 and 2007, and became the dominant genogroup in outbreaks in 2008. The antigenic variation of subgenogroup B5 strains had been discussed previously [21,29], and B1/B4, B5, and C2/C4 were divided into different groups in the antigenic map. But in another study, subgenogroup B5 was proposed to be redesignated as B4 based on the genetic analysis of complete genome nucleotide sequences [19]. More studies are needed to explain the inconsistent results between antigenic and genetic typing.

Subgenogroup C4 circulated and evolved in neighbouring countries in recent years chronologically, especially in China. There were two clusters of subgenogroup C4 strains in China from 1998 to 2008, C4b (from 1998 to 2004) and C4a (from 2003 to 2008), and the Shandong C4a strains were further divided into three lineages [26]. In Taiwan, subgenogroup C4 was first isolated in 1998 (as C4b cluster in China), and then caused outbreaks from 2004 to 2005 (as C4a cluster in China) [31]. According to the sequence analyses in this study, we identified several C4 isolates which were correlated well with C4 strains in China in 2008-2009, but not correlated with those isolated in Taiwan in 2004-2005, indicating that the virus was supposed to be transmitted from China (Figure 2). This subgenogroup caused several outbreaks in China over the last four years [26,34], but not in Taiwan, which was possibly due to herd immunity related to the subgenogroup C4 epidemic in Taiwan from 2004 to 2005. However, we still detected several subgenogroup C4 isolates in 2009 (unpublished data), and an increase of severe cases in early 2010, indicating that a potential of subgenogroup C4 outbreak in 2010 was expected, and to maintain a comprehensive surveillance system for enteroviruses seems to be a must.

Inter-genogroup, inter-subgenogroup and intra-subgenogroup average divergences of EV-71 complete genome nucleotide sequences were 17-22%, 10-14% and 1-10%, respectively [19]. However, further evidence is needed to designate the subgenogroup C2-like as a new subgenogroup. On the other hand, the lower neutralization antibody titers of subgenogroup C2-like (with a maximum of 128-fold decrease) indicated the antigenic differences with other subgenogroups (Table 2, Table 3). In previous study, a close antigenic relationship among the EV-71 isolates belonging to genogroups B and C was reported. The neutralization titers of the antisera for different genogroups of EV-71 ranged from 512 to >1,024, while the titers of the antisera for homologous EV-71 isolates were >1,024 [33]. The antigenic diversity of subgenogroup C2-like viruses displayed in this study may result in the inefficiency of herd immunity, and cause concerns on vaccine development for EV-71, e.g., monovalent or polyvalent vaccine. In addition, to further clarify the divergences, more researches using EV-71 monoclonal antibodies are needed for identification of neutralization epitopes.

The subgenogroup C2-like was supposed to be a recombinant originated from subgenogroup C2 and B3 based on a bootscan analysis. In addition, the subgenogroup C2-like viruses were isolated from different patients in different month, demonstrating that this subgenogroup was not a single case but circulated for a period of time. In Taiwan, subgenogroup C2 strains were only observed in 1998 [35], but subgenogroup B3 strains were never reported before. It is difficult to trace the actual spread route due to the recently more frequent international travel and fluxes of laborers. However, each gene region of the subgenogroup C2-like was 73.2-95.1% identical to that of other subgenogroups (Table 1), so it is supposed probably that the ancestors of this subgenogroup were imported into Taiwan before 2008, experienced recombination events, and then evolved into a unique subgenogroup. For enteroviruses, recombination was most reported to occur in the nonstructural protein region [36], while few reports demonstrated recombination in the structural capsid protein region [37]. The putative recombination breakpoint at 2B gene in this study was not reported yet. Other breakpoints at the 3'-termini of the 2A and 3C regions [38], 3D and 3'UTR regions [39] were identified in previous reports. It was speculated that the higher degree of similarity in nonstructural protein region may favor the occurrence of recombination. However, variants with recombination or deletion mutations, especially in structural protein region, may not survive or replicate less efficiently [13,40]. The subgenogroup C2-like strains showed lower CCID50 than other subgenogroups (data not shown), and it may explain why this subgenogroup did not cause outbreaks in 2008. Another possibility was that the prevalence of subgenogroup C2-like might be underestimated due to asymptomatic infections or mild illness despite a surveillance system had been set up.

Conclusions

In summary, firstly, we described a genetic and antigenic diversity of EV-71 subgenogroups in Taiwan in 2008, including 3 previously reported subgenogroups C5, B5, and C4, and one C2-like subgenogroup. Secondly, the subgenogroup C4 isolates in 2008 were genetically similar to the Chinese strains causing outbreaks in recent years, so we need to closely monitor if these subgenogroup C4 outbreaks happen or not in Taiwan in the next few years. Thirdly, due to the diversity of phylogeny, rapid changing of subgenogroups, and the potential of severe and fatal outbreaks on their way, it is a must to monitor the recombination events as well as antigenic and genetic evolution of EV-71 very attentively and carefully.

MethodsVirus isolation and identification

EV-71 viruses used in this study were collected by the surveillance systems under Centers for Disease Control, Taiwan (Taiwan CDC). These 989 strains were isolated from throat swabs, stools, sera, or cerebrospinal fluid specimens taken from patients with HFMD, herpangina, and other symptoms related to enterovirus infection. Virus isolation was carried out using rhabdomyosarcoma (RD), human diploid fibroblast (MRC-5), African green monkey kidney (Vero), human lung carcinoma (A549), monkey kidney (LLC-MK2), or human epidermoid carcinoma (HEp-2) cell lines until cytopathic effects (CPE) were observed. The isolates were then identified by immunofluorescence assay (IFA) using an EV-71 commercial monoclonal antibody (Light Diagnostic, USA). The CCID50 of the virus was calculated by the Reed and Muench method [41].

RT-PCR and Sequencing

Viral RNA was extracted according to the manufactory protocol from 140 μl of culture supernatant by QIAamp Viral RNA Mini Kit (Qiagen, Santa Clara, CA). One-step RT-PCR of VP1 gene was performed for all 989 EV-71 isolates with primer set 159/162 [14], and full-length RT-PCR was performed for two isolates (2008-07776 and 2008-00643) as described previously [13]. The products were confirmed by agarose electrophoresis and were stained with ethidium bromide. DNA was sequenced in both directions using BigDye Terminator Ready Reaction Cycle Sequencing Kit and an automated sequencer ABI 3730 (Applied Biosystems, Foster City, CA, USA).

Sequence analysis and recombination analysis

Identification and subtyping was carried out by sequence comparisons with reference EV sequences in GenBank using the BLAST [42] and confirmed by phylogenetic analysis. The DNA sequences were assembled and then aligned with reference sequences using the Clustal W program by BioEdit (version 7.0.9.0) software [43]. Phylogenetic trees were constructed using the neighbor-joining method by MEGA version 4 software [44] with 1,000 replications of bootstrap analyses. The prototype coxsackievirus A16 (CA16) G-10 strain was used as an out-group. Detection of recombination events among the subgenogroups of EV-71 viruses using the full genome sequence was determined by similarity plot and bootscan analyses using SimPlot, version 3.5.1 [45] as previously described [21,46]. The nucleotide identity was calculated using the Kimura 2-parameter method with a transition-transversion ratio of 10 [47] and a sliding window of 400 nucleotides in 20 nucleotide steps. The recombination breakpoints were determined by the maximization of χ2 analysis [48], and the p values for the resulting informative sites were calculated using the χ2 test.

Preparation of different subgenogroups EV-71 virus as immunogen for rabbits immunization

Three ancient EV-71 strains of subgenogroups C2, C5 and B5 (AFP98111207, E2006125, and E2007599, respectively) in Taiwan were selected for antiserum preparation. These strains were propagated in RD cells, and the CCID50 was determined before animal inoculations. Anti-enterovirus rabbit serum was prepared as described preciously [49]. Briefly, New Zealand White rabbits were immunized intravenously with 5 ml of UV-inactivated virus stock (>108 CCID50/ml) of above three subgenogroups of EV-71. The animals were subsequently boosted four times with the same dose at a 2-day interval, except with a double dose (10 ml) at the final boosting on day 42, and the sera were tested for neutralization antibodies on day 49.

Determination of neutralization antibody titers

Rabbit antisera and pairs of serum samples collected during the acute-phase and recovery-phase from patients with EV-71 infection were examined for neutralization antibodies. All sample determinations were performed in duplicate. Sera were first inactivated at 56ºC for 30 min, and then diluted two-fold serially in DMEM from 1:8 to 1:1,024. One-hundred CCID50 viruses (50 μl) were added to the well contained above serially diluted antiserum, and the mixtures were then incubated in a CO2 incubator at 36ºC for 60 min. Later, 100 μl of RD cell suspension containing approximately 3 × 104 cells was added to each well, and the CPE was recorded during the next 4 days. The neutralization end-point titer is defined as the highest dilution fold at which 50% of cells showing complete inhibition of CPE formation.

Statistical analysis

Differences between proportions were tested using the χ2 test. The neutralization antibody titers were compared between the subgenogroup C2-like group and other subgenogroup groups by using Student's t-test with log10-transformed data. The p value < 0.05 is taken to indicate statistically significance.

Nucleotide sequence accession numbers

The nucleotide sequences newly determined in this study have been submitted to the GenBank under the accession no. HM622381 to HM622392.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

YPH, TLL drafted the manuscript. YPH, WBF performed sequence and data analysis. TLL, YHT, CCH performed virus isolation, viral identification and neutralization test. LCH, YJC collected epidemiological information and edited the manuscript. JYY, FYC provided consultation and editing of the manuscript. HSW revised the manuscript. All authors read and approved the final manuscript.

Supplementary MaterialAdditional file 1

Phylogenetic analysis of enterovirus 71. The phylogenetic tree was constructed by the neighbor-joining method with MEGA version 4 software, and the reliabilities indicated at the branch nodes were evaluated using 1,000 bootstrap replications. Only values of over 70% were shown. The prototype coxsackievirus A16 (CA16) G-10 strain was used as an out-group. The tree was drawn based on the 5'UTR (A), VP4 (B), VP2 (C), VP3 (D), VP1 (E), 2A (F), 2B (G), 2C (H), 3A (I), 3B (J), 3C (K), and 3D (L) region nucleotide sequences.

Click here for file

Acknowledgements

We would like to thank the chiefs of Taiwan CDC Contracted Virology Laboratories for their cooperation to make this study possible. They are Chuan-Liang Kao, Jang-Jih Lu, Yu-Jiun Chan, Kuo-Chien Tsao, Ming-Jer Ding, Mu-Chin Shih, Chi-Ho Chan, Jen-Shiou Lin, Jen-Ren Wang, Kuei-Hsiang Lin, Yung-Ching Liu, Hock-Liew Eng, and Li-Kuang Chen.

This study was supported financially by research grants from Taiwan CDC and National Research Program for Genomic Medicine.

MelnickJLField BN, DMK, Howley PM, Channock RM, Melnick JL, Monath TP, Roizman B, Straus SEEnteroviruses: polioviruses, coxackieviruses, echoviruses, and newer enterovirusesFields Virology19963Philadelphia, Pa: Lippincott-Raven Publishers655712ObersteMSMaherKNixWAMicheleSMUddinMSchnurrDal-BusaidySAkoua-KoffiCPallanschMAMolecular identification of 13 new enterovirus types, EV79-88, EV97, and EV100-101, members of the species Human Enterovirus BVirus Res2007128344210.1016/j.virusres.2007.04.00117485128Enterovirus update pagehttp://www.picornaviridae.com/enterovirus/enterovirus.htmStanwayGBrownFChristianPFauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LAPicornaviridaeVirus taxonomy: classification and nomenclature of viruses Eighth report of the International Committee on the Taxonomy of Viruses2005Amsterdam, The Netherlands: Elsevier Academic Press757778New enterovirus specieshttp://www.picornastudygroup.com/types/enterovirus/ev_species.htmNixWAObersteMSPallanschMASensitive, seminested PCR amplification of VP1 sequences for direct identification of all enterovirus serotypes from original clinical specimensJ Clin Microbiol2006442698270410.1128/JCM.00542-0616891480TanCSCardosaMJHigh-titred neutralizing antibodies to human enterovirus 71 preferentially bind to the N-terminal portion of the capsid protein VP1Arch Virol20071521069107310.1007/s00705-007-0941-117318736SivasamughamLACardosaMJTanWSYusoffKRecombinant Newcastle Disease virus capsids displaying enterovirus 71 VP1 fragment induce a strong immune response in rabbitsJ Med Virol2006781096110410.1002/jmv.2066816789020FooDGAlonsoSPhoonMCRamachandranNPChowVTPohCLIdentification of neutralizing linear epitopes from the VP1 capsid protein of Enterovirus 71 using synthetic peptidesVirus Res2007125616810.1016/j.virusres.2006.12.00517222936FooDGMacaryPAAlonsoSPohCLIdentification of human CD4 T-cell epitopes on the VP1 capsid protein of enterovirus 71Viral Immunol20082121522410.1089/vim.2007.008918419256Rakoto-AndrianariveloMGumedeNJegouicSBalanantJAndriamamonjySNRabemanantsoaSBirminghamMRandriamanalinaBNkolomoniLVenterMReemergence of recombinant vaccine-derived poliovirus outbreak in MadagascarJ Infect Dis20081971427143510.1086/58769418419577BouslamaLNasriDCholletLBelguithKBourletTAouniMPozzettoBPilletSNatural recombination event within the capsid genomic region leading to a chimeric strain of human enterovirus BJ Virol2007818944895210.1128/JVI.00180-0717537864ChanYFAbuBakerSPhylogenetic evidence for inter-typic recombination in the emergence of human enterovirus 71 subgenotypesBMC Microbiol200667410.1186/1471-2180-6-7416939656BrownBAObersteMSAlexanderJPJrKennettMLPallanschMAMolecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998J Virol1999739969997510559310CardosaMJPereraDBrownBACheonDChanHMChanKPChoHMcMinnPMolecular epidemiology of human enterovirus 71 strains and recent outbreaks in the Asia-Pacific region: comparative analysis of the VP1 and VP4 genesEmerg Infect Dis2003946146812702227TuPVThaoNTPereraDHuuTKTienNTThuongTCHowOMCardosaMJMcMinnPCEpidemiologic and virologic investigation of hand, foot, and mouth disease, southern Vietnam, 2005Emerg Infect Dis2007131733174118217559van der SandenSKoopmansMUsluGvan der AvoortHEpidemiology of enterovirus 71 in the Netherlands, 1963 to 2008J Clin Microbiol2009472826283310.1128/JCM.00507-0919625480MunemuraTSaikusaMKawakamiCShimizuHOsetoMHagiwaraAKimuraHMiyamuraTGenetic diversity of enterovirus 71 isolated from cases of hand, foot and mouth disease in Yokohama City between 1982 and 2000Arch Virol200314825326310.1007/s00705-002-0920-512556991ChanYFSamICAbuBakarSPhylogenetic designation of enterovirus 71 genotypes and subgenotypes using complete genome sequencesInfect Genet Evol20101040441210.1016/j.meegid.2009.05.01019465162DeshpandeJMNadkarniSSFrancisPPEnterovirus 71 isolated from a case of acute flaccid paralysis in India represents a new genotypeCurr Sci20038413501353HuangYPLinTLKuoCYLinMWYaoCYLiaoHWHsuLCYangCFYangJYChenPJWuHSThe circulation of subgenogroups B5 and C5 of enterovirus 71 in Taiwan from 2006 to 2007Virus Res200813720621210.1016/j.virusres.2008.07.01518706461ShimizuHUtamaAOnnimalaNLiCLi-BiZYu-JieMPongsuwannaYMiyamuraTMolecular epidemiology of enterovirus 71 infection in the Western Pacific RegionPediatr Int20044623123510.1046/j.1442-200x.2004.01868.x15056257McMinnPLindsayKPereraDChanHMChanKPCardosaMJPhylogenetic analysis of enterovirus 71 strains isolated during linked epidemics in Malaysia, Singapore, and Western AustraliaJ Virol2001757732773810.1128/JVI.75.16.7732-7738.200111462047BibleJMIturriza-GomaraMMegsonBBrownDPantelidisPEarlPBendigJTongCYMolecular epidemiology of human enterovirus 71 in the United Kingdom from 1998 to 2006J Clin Microbiol2008463192320010.1128/JCM.00628-0818650362DiedrichSWeinbrechtASchreierESeroprevalence and molecular epidemiology of enterovirus 71 in GermanyArch Virol20091541139114210.1007/s00705-009-0413-x19506798ZhangYTanXJWangHYYanDMZhuSLWangDYJiFWangXJGaoYJChenLAn outbreak of hand, foot, and mouth disease associated with subgenotype C4 of human enterovirus 71 in Shandong, ChinaJ Clin Virol20094426226710.1016/j.jcv.2009.02.00219269888HoMChenERHsuKHTwuSJChenKTTsaiSFWangJRShihSRAn epidemic of enterovirus 71 infection in Taiwan. Taiwan Enterovirus Epidemic Working GroupN Engl J Med199934192993510.1056/NEJM19990923341130110498487LinTYTwuSJHoMSChangLYLeeCYEnterovirus 71 outbreaks, Taiwan: occurrence and recognitionEmerg Infect Dis2003929129312643822HuangSWHsuYWSmithDJKiangDTsaiHPLinKHWangSMLiuCCSuIJWangJRReemergence of enterovirus 71 in 2008 in Taiwan: dynamics of genetic and antigenic evolution from 1998 to 2008J Clin Microbiol2009473653366210.1128/JCM.00630-0919776232ChenSCChangHLYanTRChengYTChenKTAn eight-year study of epidemiologic features of enterovirus 71 infection in TaiwanAm J Trop Med Hyg20077718819117620652LinKHHwangKPKeGMWangCFKeLYHsuYTTungYCChuPYChenBHChenHLEvolution of EV71 genogroup in Taiwan from 1998 to 2005: an emerging of subgenogroup C4 of EV71J Med Virol20067825426210.1002/jmv.2053416372302TsengFCHuangHCChiCYLinTLLiuCCJianJWHsuLCWuHSYangJYChangYWEpidemiological survey of enterovirus infections occurring in Taiwan between 2000 and 2005: analysis of sentinel physician surveillance dataJ Med Virol2007791850186010.1002/jmv.2100617935170KungSHWangSFHuangCWHsuCCLiuHFYangJYGenetic and antigenic analyses of enterovirus 71 isolates in Taiwan during 1998-2005Clin Microbiol Infect20071378278710.1111/j.1469-0691.2007.01745.x17488328YangFRenLXiongZLiJXiaoYZhaoRHeYBuGZhouSWangJQiJEnterovirus 71 outbreak in the People's Republic of China in 2008J Clin Microbiol2009472351235210.1128/JCM.00563-0919439545ShihSRHoMSLinKHWuSLChenYTWuCNLinTYChangLYTsaoKCNingHCGenetic analysis of enterovirus 71 isolated from fatal and non-fatal cases of hand, foot and mouth disease during an epidemic in Taiwan, 1998Virus Res20006812713610.1016/S0168-1702(00)00162-310958984OprisanGCombiescuMGuillotSCaroVCombiescuADelpeyrouxFCrainicRNatural genetic recombination between co-circulating heterotypic enterovirusesJ Gen Virol2002832193220012185273BlomqvistSBruuALStenvikMHoviTCharacterization of a recombinant type 3/type 2 poliovirus isolated from a healthy vaccinee and containing a chimeric capsid protein VP1J Gen Virol20038457358010.1099/vir.0.18708-012604808HuangSCHsuYWWangHCHuangSWKiangDTsaiHPWangSMLiuCCLinKHSuIJWangJRAppearance of intratypic recombination of enterovirus 71 in Taiwan from 2002 to 2005Virus Res200813125025910.1016/j.virusres.2007.10.00218036697ChanYFAbuBakarSHuman enterovirus 71 subgenotype B3 lacks coxsackievirus A16-like neurovirulence in mice infectionVirol J200527410.1186/1743-422X-2-7416122396StanwayGHughesPJWestropGDEvansDMDunnGMinorPDSchildGCAlmondJWConstruction of poliovirus intertypic recombinants by use of cDNAJ Virol198657118711903005620ReedLJMuenchHA simple method of estimating fifty percent endpointsAm J Hyg193827493497AltschulSFGishWMillerWMyersEWLipmanDJBasic local alignment search toolJ Mol Biol19902154034102231712HallTABioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NTNucleic Acids Symp Ser1999419598TamuraKDudleyJNeiMKumarSMEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0Mol Biol Evol2007241596159910.1093/molbev/msm09217488738LoleKSBollingerRCParanjapeRSGadkariDKulkarniSSNovakNGIngersollRSheppardHWRaySCFull-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombinationJ Virol1999731521609847317ObersteMSPenarandaSMaherKPallanschMAComplete genome sequences of all members of the species Human enterovirus AJ Gen Virol2004851597160710.1099/vir.0.79789-015166444SalminenMOCarrJKBurkeDSMcCutchanFEIdentification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanningAIDS Res Hum Retroviruses1995111423142510.1089/aid.1995.11.14238573403RobertsonDLHahnBHSharpPMRecombination in AIDS virusesJ Mol Evol19954024925910.1007/BF001632307723052LinTLLiYSHuangCWHsuCCWuHSTsengTCYangCFRapid and highly sensitive coxsackievirus a indirect immunofluorescence assay typing kit for enterovirus serotypingJ Clin Microbiol20084678578810.1128/JCM.01114-0718032614