Phylogenetic evidence for inter-typic recombination in the emergence of human enterovirus 71 subgenotypes
© Yoke-Fun and AbuBakar; licensee BioMed Central Ltd. 2006
Received: 22 June 2006
Accepted: 30 August 2006
Published: 30 August 2006
Human enterovirus 71 (EV-71) is a common causative agent of hand, foot and mouth disease (HFMD). In recent years, the virus has caused several outbreaks with high numbers of deaths and severe neurological complications. Several new EV-71 subgenotypes were identified from these outbreaks. The mechanisms that contributed to the emergence of these subgenotypes are unknown.
Six EV-71 isolates from an outbreak in Malaysia, in 1997, were sequenced completely. These isolates were identified as EV-71 subgenotypes, B3, B4 and C2. A phylogenetic tree that correlated well with the present enterovirus classification scheme was established using these full genome sequences and all other available full genome sequences of EV-71 and human enterovirus A (HEV-A). Using the 5' UTR, P2 and P3 genomic regions, however, isolates of EV-71 subgenotypes B3 and C4 segregated away from other EV-71 subgenotypes into a cluster together with coxsackievirus A16 (CV-A16/G10) and EV-71 subgenotype C2 clustered with CV-A8. Results from the similarity plot analyses supported the clustering of these isolates with other HEV-A. In contrast, at the same genomic regions, a CV-A16 isolate, Tainan5079, clustered with EV-71. This suggests that amongst EV-71 and CV-A16, only the structural genes were conserved. The 3' end of the virus genome varied and consisted of sequences highly similar to various HEV-A viruses. Numerous recombination crossover breakpoints were identified within the non-structural genes of some of these newer EV-71 subgenotypes.
Phylogenetic evidence obtained from analyses of the full genome sequence supports the possible occurrence of inter-typic recombination involving EV-71 and various HEV-A, including CV-A16, the most common causal agent of HFMD. It is suggested that these recombination events played important roles in the emergence of the various EV-71 subgenotypes.
More than sixty different types of human enteroviruses have been identified . Species human enterovirus A (HEV-A) which include 11 members of the coxsackievirus A (CV-A) group; CV-A2, CV-A3, CV-A4, CV-A5, CV-A6, CV-A7, CV-A8, CV-A10, CV-12, CV-14, CV-A16 and human enterovirus 71 (EV-71) are associated with several human diseases [1, 2]. CV-A2, CV-A6, CV-A8 and CV-A10 are known to cause herpangina ; CV-A2, CV-A4, CV-A7 and CV-A10 cause aseptic meningitis  and CV-A5, CV-A10, CV-A16 and EV-71 usually cause hand, foot and mouth disease (HFMD) . Amongst these HEV-A viruses, only EV-71 is associated with major outbreaks involving infections with severe neurological complications and deaths. Recent outbreaks in Malaysia, Taiwan and Vietnam claimed more than 130 deaths [3, 4]. In these outbreaks, EV-71 was isolated from patients' specimens and post-mortem tissues, suggesting its direct involvement as the etiologic agent [5, 6]. Molecular epidemiological investigations performed using the partial gene sequences of all available EV-71 isolates suggested the presence of several EV-71 lineages co-circulating during most of the Asian outbreaks [3, 6–11]. EV-71 subgenotypes B3, B4, C1 and C2 were reported in Malaysia with subgenotype B3 being the most prevalent in the 1997 outbreak [3, 7–9]. The subgenotype B4 was present initially at a low level but became endemic and eventually replaced subgenotype B3 as the predominant virus in the subsequent outbreaks in Malaysia and the region [7, 9]. Neither the B3 nor B4 subgenotype, however, was associated with any fatal infections outside the Malaysia 1997 outbreak. Similarly, the C2 isolates were identified as the main causal agent of the HFMD-associated severe and fatal infections during the 1998 Taiwan outbreak [6, 10]. These isolates also caused severe infections in the Perth 1999 outbreak [7, 9] but no report of serious infection or death attributed to these isolates was recorded from the Malaysia 1997 outbreak [8, 12]. Recently, two new EV-71 subgenotypes, C3 and C4 were proposed to accommodate EV-71 strains isolated from cases of aseptic meningitis and paralysis in Korea and China, respectively [7, 13, 14]. There were no reports, however, that the subgenotype C4 isolates were associated with severe infections [13, 14]. The emergence of the various EV-71 subgenotypes with varying degree of clinical severity is puzzling. The viruses together with other HEV-A viruses are known to be present in the region and usually cause the typical mild HFMD [7–9, 15]. In the present study, the potential mechanisms that could play important role in the emergence of the different EV-71 subgenotypes are phylogenetically investigated.
Phylogenetic relationships between EV-71 and other HEV-A viruses
The complete genome sequences of all the EV-71 isolates were aligned with other available HEV-A complete genome sequences obtained from the Genbank. The total length of the genomes varied between 7395–7434 nucleotides and they encode for 2189–2201 amino acids. The differences in genome length were caused by the presence of deletions or insertions within the 5' UTR, capsid gene and 3' UTR sequences. The capsid genome region of EV-71 was longer (1 – 4 amino acids) or of the same length with other HEV-A viruses, with the exception of CV-A4 and CV-A6 which have eight extra amino acid insertions. No deletions and insertions were observed in the P2 and P3 genomic regions, hence, the gene length within these genomic regions remained constant for all the HEV-A viruses.
At the P2 genomic region, EV-71 BrCr clustered with CV-A2, CV-A3, CV-A6, CV-A10 and CV-A12, while the EV-71 subgenotype C2 isolates were phylogenetically closer to CV-A8 (Figure 2c). EV-71 isolates of subgenotype B2 and C4 segregated from all other EV-71 and formed a cluster together with CV-A4, CV-A5, CV-A14 and CV-A16/G10 (Figure 2c). CV-A16 Tainan5079 segregated away from the prototype CV-A16/G10 and grouped with EV-71 BrCr, CV-A2, CV-A3, CV-A6, CV-A10 and CV-A12 (100% bootstrap support, Figure 2c). At the P3 genomic region, the different EV-71 lineages formed four distinct clusters together with other HEV-A viruses. Isolates of subgenotype C2 were phylogenetically closer to CV-A8 and this was supported by a high bootstrap value (100%). Isolates of subgenotype B3 and C4 isolates clustered with CV-A4, CV-A14 and CV-A16 and support for the clustering was also significant (100% bootstrap). EV-71 isolates of subgenotype B2 and B4 were phylogenetically closer to CV-A5, although this did not receive high bootstrap support (65%). The genotype A isolate BrCr segregated from all the recent EV-71 subgenotypes and clustered with CV-A2, CV-A3, CV-A6, CV-A10 and CV-A12 (Figure 2d). At the P3 genomic region, isolate CV-A16 Tainan5079 segregated away from CV-A16/G10 and grouped with EV-71 BrCr, CV-A2, CV-A3, CV-A6, CV-A10 and CV-A12 (Figure 2d).
An unrooted phylogenetic tree drawn using only the 3' UTR showed clustering of EV-71 isolates of subgenotype of B3 and C4 with CV-A4, CV-A14 and CV-A16 (95% bootstrap support, not shown). No significant segregation (< 30% bootstrap support), however, was observed for the remaining isolates and this was perhaps due to the short sequence length of the 3' UTR. The clustering of isolates of subgenotype B3 and C4 with CV-A4, CV-A14 and CV-A16 at the 3' UTR genomic region was consistent with the previous clustering at the P2 and P3 genomic regions. Based on these results, it appeared that sequences of genes at the 3' half of the EV-71 genome contributed to the multiple and diverse EV-71 subgenotypes and these genes showed high similarity to different HEV-A viruses.
Inter-typic recombination between EV-71 and other HEV-A viruses
Support for inter-typic recombination involving CV-A16 and EV-71 were also obtained from the phylogenetic investigation of the CV-A16 Tainan5079. Several recombination breakpoints were revealed when the CV-A16 Tainan5079 whole genome sequence was queried against several potential parental genome sequences (Figure 4c). A major recombination breakpoint was identified between nucleotides 3616–3781 within the 2A gene (≥ 84% sequence similarity), hence switching the CV-A16 Tainan5079 genome sequences to that of the non-structural gene sequences of EV-71 BrCr (nucleotides 3781–7409). Two additional breakpoints were noted between nucleotides 168–193 and 513–584 of the 5' UTR sequences switching the genome sequence from CV-A16-like to that of EV-71 subgenotype C2 (Figure 4c). Acquisition of the EV-71 subgenotype C2 sequences also resulted in the CV-A16 Tainan5079 isolate acquiring the EV-71-like 5' UTR RNA stem loop structure (not shown). When the amino acid sequences of the non-structural genes were examined, at least 45 amino acids were found to be similar to that encoded by the ancestral EV-71 BrCr and not to CV-A16 G/10. Of these amino acids, 21 amino acid differences occurred within the 3D RNA polymerase gene alone, suggesting that the gene was very much EV-71-like than CV-A16 (not shown). These results revealed that the recently isolated CV-A16 had most of its genome 5' half including all the structural genes except the 5' UTR consisting of CV-A16-like sequences and the remaining 3' half of the genome consisting of the supposedly no longer circulating EV-71 genotype A BrCr sequence.
Results from the full genome sequence similarity plot analyses and phylogenetic investigations presented here support the possible occurrence of inter-typic recombination involving EV-71 and various HEV-A, including CV-A16, the common causal agent of HFMD. All the EV-71 isolates examined had conserved 5' half of the genome, which consisted of mainly the structural genes including VP1, VP2, VP3 and VP4. These structural gene sequences particularly VP1 and VP4 when used to construct the phylogenetic trees clustered all EV-71 isolates into a unique clade of its own, away from other HEV-A viruses. This justifies usage of these genes for the molecular identification of EV isolates [19, 20]. It is possible that the structural genes especially the virus capsid gene, carries genetic restrictions required to ensure virus receptor recognition, binding and entry into host cells, and hence remain highly conserved. CV-A16-like sequences were present within the genome 3' end of EV-71 subgenotype B3 and C4. Conversely, EV-71-like sequences were also found within the genome 3' end of CV-A16 Tainan5079. These highlight the possible occurrence of recombination events in the wild between the two common causal agents of HFMD. This possibility is not difficult to envisage considering that both EV-71 and CV-A16 viruses are known to co-circulate within the same population and cause a common disease. It is particularly true in parts of the world where the level of hygiene practices and sanitation are lacking or less developed, as in these circumstances the possibility of co-infection is expected to be higher.
Results from the present study also suggest that inter-typic recombination events occurred between EV-71, CV-A16 and other HEV-A viruses. The incongruent phylogeny patterns observed suggest the possible occurrence of recombination in EV-71 in a similar manner to those previously reported for a number of enteroviruses [21–27]. It is reported that most of the current circulating HEV-B viruses are actually recombinants relative to its prototype strains isolated in the 1950s [23, 27]. Recombination amongst HEV-A viruses has also been described in CV-A6, CV-A10 and CV-A12  and a novel group of HEV-A consisting of EV-76, EV-89, EV-90 and EV-91 . The P1 genomic region of EV-70, a virus that causes acute haemorrhagic conjunctivitis is highly similar to EV and rhinovirus. In addition, the P2 and P3 genomic regions are phylogenetically closer to CV-B and swine vesicular disease virus . Another EV, CV-B5, has its major capsid proteins that are most closely related to swine vesicular disease virus, with the 5' and 3' UTR and the P3 genomic region that are highly similar to the corresponding regions of other CV-B viruses . Other example of EV with highly diverse genome sequence origin includes the prototype EV-73 isolated in California in 1988. This EV has P1 genomic region that is highly similar to echovirus 6, echovirus 29, echovirus 21 and echovirus 30, with the P2 and P3 genomic regions that share high similarity to CV-B3 . Together, these findings reiterate the suggestion that recombination is one of the potential mechanisms by which genetic diversity, is generated in EV-71 and other HEV-A viruses especially within the non-structural genes where recombination hot spots are mainly located [21, 23, 32–35]. While recombination could occur within the structural gene region, it perhaps yielded nonviable or unstable progeny viruses that did not survive and this could be due to generation of less fit variants with deleterious mutations . Selection of a better fit virus is likely to have contributed towards the demise of the recombinant EV-71 B3 virus once found prevalent in Peninsula Malaysia and Sarawak during the 1997 outbreak, Singapore and Japan in 1997 [7, 8] and Perth in 1999 . This virus was shown to replicate less efficiently in comparison to EV-71 B4 viruses , hence it could be the reason why it was rapidly replaced with EV-71 B4 viruses which remained endemic in the region for at least up to 2002, causing infections in Sarawak, Singapore, Taiwan and Japan [7, 8]. The acquisition of a segment of genes from CV-A16 by EV-71 B3 virus [21, 37] could have rendered the virus less fit to adapt to a new environment especially the host immunity.
Phylogenetic evidence of possible inter-typic recombination amongst HEV-A viruses including between the two causal agents of HFMD, EV-71 and CV-A16, were obtained from analyses using the virus isolates' full genome sequences. The finding highlights the potential role of recombination in the emergence of the newer EV-71 lineages including those perhaps with varied virulence and clinical manifestations. Active molecular surveillance of EV-71 and CV-A16 would help to ensure rapid identification of emergence of any potentially highly virulent EV-71 or CV-A16 strains. This is particularly pertinent in regions where recurring HFMD outbreaks involving the different HEV-A have occurred.
RNA isolation and sequencing
EV-71 and HEV-A sequences used for sequence analysis and construction of phylogenetic trees.
Paralytic illness +
Not known +
Bulbar poliomyelitis +
Sore throat, nuchal rigidity, muscle weakness +
Flies in polio community +
Oligonucleotide primers for cDNA synthesis and PCR.
Nucleotide sequence (5' to 3')
Nucleotide position *
TAA AAC AGC YKK KGG KTK KYA CC
TGM ACR TGR ATG CAR AAC C
CCD GAT GTR YTR ACW GAA ACC
ATC ACA TCT GCC ACY CTR TCY
GGR GCR CCY AAY ACW GCY TAY ATW WTR GC
CYT GRG CRG TRG TWG AHG AHA CKA RB
CAR GTY TCH GTB CCR TTY ATG TCA CC
GCT GTY TTB GMY TTR AYC CAV GC
YAA CTC WMD SAG AAA RCA CTA YCC
GWS GAK GCA ATS ACA AAC TTD GAK G
CCW GAY CAY TTT GAY GGR TAY MAR C
AHC CRT ACT GVV CVA CRT CHY C
GTT GAG CGH CAC YTS AAY AGA GY
GYA RRT CYA ACC CAT ACT TRT CC
CAT YAA GAA RAR RGA CAT YYT BG
TTG CTA TTC TGG TTA TAA CAA ATT TAC CC
Overlapping DNA sequences with at least 85% sequence homology and a minimum of 20 overlaps were assembled into contigs to generate consensus sequences using Sequencher version 4.0.5 (Gene Codes Corporation, USA). The consensus sequences were aligned against other EV-71 and HEV-A complete genome sequences retrieved from the GenBank (Table 1) using Clustal X version 1.81 . All completely sequenced EV-71 isolates available in the GenBank and the six isolates obtained from the present study were used for construction of phylogenetic trees. The genetic distance was determined by a pairwise estimation of the sequences percent divergence. Positions with gaps were included, transition and transversion ratio was fixed at 10 and corrections for multiple substitutions were made. These options were chosen to take into account the ambiguous part of the alignments and to correct for more than one substitutions happening at many sites that may underestimate the actual genetic distances . Phylogenetic trees presented here were constructed using the neighbour-joining method . All the unrooted trees were displayed using TreeView version 1.66 . The strength of the phylogenetic trees was estimated by bootstrap analyses using 1000 random samplings. A bootstrap value of ≥ 70% indicates a strong support for the tree topology . Amino acid sequences were examined after stripping the 5' UTR and 3' UTR sequences and consensus sequences for each EV-71 genotype were established.
Potential recombination amongst the isolates was determined as previously described . Initially, the neighbour-joining trees were constructed using the various gene regions of the complete genome. The trees were manually inspected to determine if there was any shifting of the tree topology in comparison to that constructed using the whole genome sequence. The putative recombinant isolate, determined from the shifting of the tree topology, was then queried against all the existing HEV-A virus genomes. Initially, all the 31 manually corrected and aligned complete nucleotide sequences were used. Isolates that shared similar patterns were grouped together to generate consensus sequences. Recombination was identified when conflicting genome sequence profiles appeared suggesting acquisition of sequences from a different parental genotype. Phylogenetic trees were constructed for each of the putative recombinant sequences using the neighbour-joining tree algorithm and maximum likelihood distance model and support for the tree topology was determined by bootscanning analyses that utilized the bootstrapping procedures with 100 resamplings [42, 43]. Bootstrap values of > 70% were used to indicate robust support for the tree topology . The crossover breakpoints were identified when χ2 values were at the maximum .
For the analysis of EV-71 and HEV-A viruses, all gaps were stripped and the sequences were corrected for multiple substitutions, transition to tranversion ratio of ten was used and 50% consensus files were used to exclude the poorly conserved sites. A sliding window of 1000 nucleotides was moved in increments of 20 nucleotides at a time and the resulting similarity values were plotted along the virus genome. Transition and transversion ratio of 10 was chosen based on that used for HEV-A viruses to avoid underestimation of actual genetic distances . Both the similarity plot and bootscanning analysis sliding windows were based on that used for HEV-B to increase the phylogenetic signals and reduce the noise of unrelated phylogeny . However, for the analysis amongst EV-71, the sequences were not corrected for multiple substitutions, and a transition to tranversion ratio of two and a sliding window of 400 nucleotides was used to map the possible recombination breakpoints .
- HFMD Hand:
foot and mouth disease
This study is funded in parts by a grant from the Ministry of Science, Technology and Innovation, Malaysia # 06-02-09-001-BTK/TD/002.
- Poyry T, Kinnunen L, Hyypia T, Brown B, Horsnell C, Hovi T, Stanway G: Genetic and phylogenetic clustering of enteroviruses. J Gen Virol. 1996, 77: 1699-1717.View ArticlePubMedGoogle Scholar
- Melnick JL: Enteroviruses: polioviruses, coxsackieviruses, echoviruses and newer enteroviruses. Fields Virology. Edited by: Fields BN, Knipe DM, Howley PM, Chanock RM, Melnick JL, Moath TP, Roizman B, Straus SE. 1996, Philadelphia: Lippincott-Raven Publishers, 655-711. 3Google Scholar
- AbuBakar S, Chee HY, Al-Kobaisi MF, Xiaoshan J, Chua KB, Lam SK: Identification of enterovirus 71 isolates from an outbreak of hand, foot and mouth disease (HFMD) with fatal cases of encephalomyelitis in Malaysia. Virus Res. 1999, 61: 1-9. 10.1016/S0168-1702(99)00019-2.View ArticlePubMedGoogle Scholar
- Wu TN, Tsai SF, Li SF, Lee TF, Huang TM, Wang ML, Hsu KH, Shen CY: Sentinel surveillance for enterovirus 71, Taiwan, 1998. Emerg Infect Dis. 1999, 5: 458-460.PubMed CentralView ArticlePubMedGoogle Scholar
- Lum LCS, Wong KT, Lam SK, Chua KB, Goh AYT, Lim WL, Ong BB, Paul G, AbuBakar S, Lambert M: Fatal enterovirus 71 encephalomyelitis. J Pediatr. 1998, 133: 795-798. 10.1016/S0022-3476(98)70155-6.View ArticlePubMedGoogle Scholar
- Yan JJ, Wang JR, Liu CC, Yang HB, Su IJ: An outbreak of enterovirus 71 in Taiwan 1998: a comprehensive pathological, virological, and molecular study on a case of fulminant encephalitis. J Clin Virol. 2000, 17: 13-22. 10.1016/S1386-6532(00)00067-6.View ArticlePubMedGoogle Scholar
- Cardosa MJ, Perera D, Brown BA, Cheon D, Chan HM, Chan KP, Cho P, McMinn P: 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.View ArticlePubMedGoogle Scholar
- Herrero LJ, Lee CSM, Hurrelbrink RJ, Chua BH, Chua KB, McMinn PC: Molecular epidemiology of enterovirus 71 in peninsular Malaysia, 1997–2000. Arch Virol. 2003, 148: 1369-1385. 10.1007/s00705-003-0100-2.View ArticlePubMedGoogle Scholar
- McMinn P, Lindsay K, Perera D, Chan HM, Chan KP, Cardosa MJ: Phylogenetic analysis of enterovirus 71 strains isolated during linked epidemics in Malaysia, Singapore, and Western Australia. J Virol. 2003, 75: 7732-7738. 10.1128/JVI.75.16.7732-7738.2001.View ArticleGoogle Scholar
- Chu PY, Lin KH, Hwang KP, Chou LC, Wang CF, Shih SR, Wang JR, Shimada Y, Ishiko H: Molecular epidemiology of enterovirus 71 in Taiwan. Arch Virol. 2001, 146: 589-600. 10.1007/s007050170164.View ArticlePubMedGoogle Scholar
- AbuBakar S, Chee HY, Shafee N, Chua KB, Lam SK: Molecular detection of enteroviruses from an outbreak of hand, foot and mouth disease in Malaysia in 1997. Scand J Infect Dis. 1999, 31: 331-335. 10.1080/00365549950163734.View ArticlePubMedGoogle Scholar
- AbuBakar S, Chan YF, Lam SK: Outbreaks of enterovirus 71 infection. N Engl J Med. 2000, 342: 355-356. 10.1056/NEJM200002033420513.View ArticlePubMedGoogle Scholar
- Shimizu H, Utama A, Onnimala N, Li C, Li-Bi Z, Yu-Jie M, Pongsuwanna Y, Miyamura T: Molecular epidemiology of enterovirus 71 infection in the Western Pacific region. Pediatr Int. 2004, 46: 231-235. 10.1046/j.1442-200x.2004.01868.x.View ArticlePubMedGoogle Scholar
- Mizuta K, Abiko C, Murata T, Matsuzaki Y, Itagaki T, Sanjoh K, Sakamoto M, Hongo S, Murayama S, Hayasaka K: 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. 10.1128/JCM.43.12.6171-6175.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Infectious agent surveillance report (IASR): Virus isolation/detection from hand, foot and mouth disease cases, 2002-2006 (Infectious Agent Surveillance Report: Data based on the reports received before September 6, 2006 from public health institutes). [https://hasseidoko.mhlw.go.jp/Byogentai/Pdf/data37e.pdf]
- Brown BA, Oberste MS, Alexander JP, Kennett ML, Pallansch MA: Molecular epidemiology and evolution of enterovirus 71 strains isolated from 1970 to 1998. J Virol. 1999, 73: 9969-9975.PubMed CentralPubMedGoogle Scholar
- Singh S, Chow VTK, Chan KP, Ling AE, Poh CL: RT-PCR, nucleotide, amino acid phylogenetic analyses of enterovirus type 71 strains from Asia. J Virol Methods. 2000, 88: 193-204. 10.1016/S0166-0934(00)00185-3.View ArticlePubMedGoogle Scholar
- Oberste MS, Penaranda S, Maher K, Pallansch MA: Complete genome sequences of all members of the species Human enterovirus A. J Gen Virol. 2004, 85: 1597-1607. 10.1099/vir.0.79789-0.View ArticlePubMedGoogle Scholar
- Oberste MS, Maher K, Kilpatrick DR, Pallansch MA: Molecular evolution of human enteroviruses: Correlation of serotype with VP1 sequence and application to picornavirus classification. J Virol. 1999, 73: 1941-1948.PubMed CentralPubMedGoogle Scholar
- Oberste MS, Maher K, Kilpatrick DR, Flemister MR, Brown BA, Pallansch MA: Typing of human enteroviruses by partial sequencing of VP1. J Clin Microbiol. 1999, 37: 1288-1293.PubMed CentralPubMedGoogle Scholar
- Chan YF, AbuBakar S: Recombinant human enterovirus 71 in hand, foot and mouth disease patients. Emerg Infect Dis. 2004, 10: 1468-1470.View ArticlePubMedGoogle Scholar
- Lindberg AM, Andersson P, Savolainen C, Mulders MN, Hovi T: Evolution of the genome of Human enterovirus B: incongruence between phylogenies of the VP1 and 3CD regions indicates frequent recombination within the species. J Gen Virol. 2003, 84: 1223-1235. 10.1099/vir.0.18971-0.View ArticlePubMedGoogle Scholar
- Lukashev AN, Lashkevich VA, Ivanova OE, Koroleva GA, Hinkkanen AE, Ilonen J: Recombination in circulating enteroviruses. J Virol. 2003, 77: 10423-10431. 10.1128/JVI.77.19.10423-10431.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Kew O, Morris-Glasgow V, Landaverde M, Burns C, Shaw J, Garib Z, Andre J, Blackman E, Freeman CJ, Jorba J, Sutter R, Tambini G, Venczel L, Pedreira C, Laender F, Shimizu H, Yoneyama Y, Miyamura T, van Der Avoort H, Obertse MS, Kilpatrick D, Cochi S, Pallansch MA, de Quadros C: Outbreak of poliomyelitis in Hispaniola associated with circulating type 1 vaccine-derived poliovirus. Science. 2002, 296: 356-359. 10.1126/science.1068284.View ArticlePubMedGoogle Scholar
- Santti J, Hyypia T, Kinnunen L, Salminen M: Evidence of recombination among enteroviruses. J Virol. 1999, 73: 8741-8749.PubMed CentralPubMedGoogle Scholar
- Oberste MS, Maher K, Pallansch MA: Evidence for frequent recombination within species human enterovirus B based on complete genomic sequences of all thirty-seven serotypes. J Virol. 2004, 78: 855-867. 10.1128/JVI.78.2.855-867.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Oberste MS, Penaranda S, Pallansch MA: RNA recombination plays a major role in genomic change during circulation of coxsackie B viruses. J Virol. 2004, 78: 2948-2955. 10.1128/JVI.78.6.2948-2955.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Oberste MS, Maher K, Michele SM, Belliot G, Uddin M, Pallansch MA: Enterovirus 76, 89, 90 and 91 represent a novel group within the species human enterovirus A. J Gen Virol. 2005, 86: 445-451. 10.1099/vir.0.80475-0.View ArticlePubMedGoogle Scholar
- Ryan MD, Jenkins O, Hughes PJ, Brown A, Knowles NJ, Booth D, Minor PD, Almond JW: The complete nucleotide sequence of enterovirus type 70: relationships with other members of the picornaviridae. J Gen Virol. 1990, 71: 2291-2299.View ArticlePubMedGoogle Scholar
- Zhang G, Wilsden G, Knowles NJ, McCauley JW: Complete nucleotide sequence of a coxsackievirus B5 virus and its relationship to swine vesicular disease virus. J Gen Virol. 1993, 74: 845-853.View ArticlePubMedGoogle Scholar
- Norder H, Bjerregaard L, Magnius LO: Open reading frame sequence of an Asian enterovirus 73 strain reveals that the prototype from California is recombinant. J Gen Virol. 2002, 83: 1721-1728.View ArticlePubMedGoogle Scholar
- Andersson P, Edman K, Lindberg AM: Molecular analysis of the echovirus 18 prototype: evidence of interserotypic recombination with echovirus 9. Virus Res. 2002, 85: 71-83. 10.1016/S0168-1702(02)00019-9.View ArticlePubMedGoogle Scholar
- Lukashev AN, Lashkevich VA, Koroleva GA, Ilonen J, Hinkkanen AE: Recombination in uveitis-causing enterovirus strains. J Gen Virol. 2004, 85: 463-470. 10.1099/vir.0.19469-0.View ArticlePubMedGoogle Scholar
- Oprisan G, Combiescu M, Guillot S, Caro V, Combiescu A, Delpeyroux F, Crainic R: Natural genetic recombination between co-circulating heterotypic enteroviruses. J Gen Virol. 2002, 83: 2193-2200.View ArticlePubMedGoogle Scholar
- Simmonds P, Welch J: Frequency and dynamics of recombination within different species of human enteroviruses. J Virol. 2006, 80: 483-493. 10.1128/JVI.80.1.483-493.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- Stanway G, Hughes PJ, Westrop GD, Evans DM, Dunn G, Minor PD, Schild GC, Almond JW: Construction of poliovirus intertypic recombinants by the use of cDNA. J Virol. 1986, 57: 1187-1190.PubMed CentralPubMedGoogle Scholar
- Chan YF, AbuBakar S: Human enterovirus 71 subgenotype B3 lacks coxsackievirus A16-like neurovirulence in mice infection. Virology Journal. 2005, 2: 74-10.1186/1743-422X-2-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nuc Acids Res. 1997, 24: 4876-4882. 10.1093/nar/25.24.4876.View ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Page RD: TreeView: An application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999, 73: 152-160.PubMed CentralPubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Robertson DL, Hahn BH, Sharp PM: Recombination in AIDS viruses. J Mol Evol. 1995, 20: 249-259. 10.1007/BF00163230.View ArticleGoogle Scholar
- Singh S, Poh CL, Chow VTK: Complete sequence analyses of enterovirus 71 strains from fatal and non-fatal cases of the hand, foot and mouth disease outbreak in Singapore (2000). Microbiol Immunol. 2002, 46: 801-808.View ArticlePubMedGoogle Scholar
- Brown BA, Pallansch MA: Complete nucleotide sequence of enterovirus 71 is distinct from poliovirus. Virus Res. 1995, 39: 195-205. 10.1016/0168-1702(95)00087-9.View ArticlePubMedGoogle Scholar
- Shih SR, Ho MS, Lin KH, Wu SL, Chen YT, Wu CN, Lin TY, Chang LY, Tsao KC, Ning HC, Chang PY, Jung SM, Hsueh C, Chang KS: Genetic analysis of enterovirus 71 isolated from fatal and non-fatal cases of hand, foot and mouth disease during an epidemic in Taiwan, 1998. Virus Res. 2001, 68: 127-136. 10.1016/S0168-1702(00)00162-3.View ArticleGoogle Scholar
- Yang JY, Sun HC, Tseng TC, Chen HY: Complete nucleotide sequence of enterovirus 71 isolated from a patient with hand-foot-mouth disease in Taiwan. 1999, UnpublishedGoogle Scholar
- Yan JJ, Su IJ, Chen PF, Liu CC, Yu CK, Wang JR: Complete genome analysis of enterovirus 71 isolated from an outbreak in Taiwan and rapid identification of enterovirus 71 and coxsackievirus A16 by RT-PCR. J Med Virol. 2001, 65: 331-339. 10.1002/jmv.2038.View ArticlePubMedGoogle Scholar
- Yang F, Jin Q, He Y, Li L, Hou Y: The complete genome of enterovirus 71 China strain. Sci China Ser C. 2001, 44: 178-183.View ArticleGoogle Scholar
- Zhou S, Li L, He Y, Jin Q: The complete genome sequencing of EV71 SHZH03 isolated in China. 2003, UnpublishedGoogle Scholar
- Poyry T, Hyypia T, Horsnell C, Kinnunen L, Hovi T, Stanway G: Molecular analysis of coxsackievirus A16 reveals a new genetic group of enteroviruses. Virology. 1994, 202: 982-987. 10.1006/viro.1994.1423.View ArticlePubMedGoogle Scholar
This article is published under license to 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.