- Open Access
Precise genotyping and recombination detection of Enterovirus
- Chieh-Hua Lin†1, 2,
- Yu-Bin Wang†3, 4,
- Shu-Hwa Chen†3,
- Chao Agnes Hsiung1, 2Email author and
- Chung-Yen Lin1, 3, 5Email author
© Lin et al. 2015
- Published: 9 December 2015
Enteroviruses (EV) with different genotypes cause diverse infectious diseases in humans and mammals. A correct EV typing result is crucial for effective medical treatment and disease control; however, the emergence of novel viral strains has impaired the performance of available diagnostic tools. Here, we present a web-based tool, named EVIDENCE (EnteroVirus In DEep conception, http://symbiont.iis.sinica.edu.tw/evidence), for EV genotyping and recombination detection. We introduce the idea of using mixed-ranking scores to evaluate the fitness of prototypes based on relatedness and on the genome regions of interest. Using phylogenetic methods, the most possible genotype is determined based on the closest neighbor among the selected references. To detect possible recombination events, EVIDENCE calculates the sequence distance and phylogenetic relationship among sequences of all sliding windows scanning over the whole genome. Detected recombination events are plotted in an interactive figure for viewing of fine details. In addition, all EV sequences available in GenBank were collected and revised using the latest classification and nomenclature of EV in EVIDENCE. These sequences are built into the database and are retrieved in an indexed catalog, or can be searched for by keywords or by sequence similarity. EVIDENCE is the first web-based tool containing pipelines for genotyping and recombination detection, with updated, built-in, and complete reference sequences to improve sensitivity and specificity. The use of EVIDENCE can accelerate genotype identification, aiding clinical diagnosis and enhancing our understanding of EV evolution.
The Enterovirus (EV) genus (family Picornaviridae) contains twelve species, including Enterovirus A to H and J, and Rhinovirus A to C. These viruses cause a wide range of diseases in humans and mammals. The single-stranded RNA genome of EV contains a single open reading frame (ORF) flanked by 5' and 3' untranslated regions (UTRs). The ORF encodes a polyprotein, which is further processed into 11 proteins: VP1-4 (structural proteins), and 2A-2C and 3A-3D (non-structural proteins) . The genetic diversity of EVs arises from the accumulation of single-base changes during viral propagation, as well as from recombination events that cause genome segments to be swapped between or within EV genotypes. To date, 308 Enterovirus genotypes have been reported (http://www.picornaviridae.com/enterovirus/enterovirus.htm, on 2015/04), and the number is rising.
Different enterovirus genotypes cause different clinical symptoms . Classical serotyping methods, such as serum neutralizing test and immunofluorescent assay, are not sufficient to specify all genotypes. For example, Tsao et al. (2010) reported that 15~30% of EV isolates failed to be serotyped in Taiwan . To overcome this problem, many clinicians have turned to sequence-based molecular typing methods, which assign viral genotypes based on nucleotide sequences; such techniques are more successful at resolving EV isolates to the corresponding genotype, and also provide rapid diagnosis . The VP1 capsid-coding region has been suggested to be the most suitable region for EV genome genotyping [4, 5]. In addition, the 5'UTR [6, 7], VP2 [8, 9], VP4 [10, 11] and 3D [11, 12] regions, as well as combinations of more than two regions, including the 5'UTR and VP4/VP2 , the 5'UTR and VP1, and VP1 and 3D , have been evaluated for their usefulness for improving the sensitivity and specificity of diagnosis. However, incongruent results may be obtained from different typing methods based on either single or multiple coding regions of the genome [16–19].
At present, there are two EV genotyping tools: enterovirus genotyping tool (version 0.1; National Institute of Public Health and the Environment (RIVM), the Netherlands)  and the genotyping tool of the NCBI . Both of these resolve genotypes on the VP1 region, and disregard the rest of the EV genome. This approach limits the ability to distinguish between strains that originated from recombination events. Moreover, EV genotype reference sequences are never updated in these libraries.
A fast, highly sensitive, and specific molecular typing tool is essential for clinical diagnosis and medical treatment. In this study, we developed a web tool, EVIDENCE (EnteroVirus In DEep coNCEption), a workbench for phylogenetic-based genotyping and recombination detection in EV genomes. Up-to-date EV classification data, nomenclature, and GenBank accession numbers for each genotype's prototype sequence were collected from the Picornaviridae Study Group website at http://www.picornaviridae.com/, and these were combined with sequences collected from the NCBI to build the genotyping reference set (GTRefSet). Phylogenetic inference was used to resolve the best-fit genotype of novel EV sequences using single or multiple genomic regions of interest. For detection of recombination events, the closeness between the suspected recombinant and reference sequences was measured as bootscanning supports by the phylogenetic method, and as sequence similarity by the distance method. The pipeline design enables users to seamlessly run recombination analyses with guidance for the choice of references. Furthermore, we revised the EV sequences in GenBank to standardize the nomenclature and to clarify genotype assignments. The collected sequences were built into the database, and can be retrieved in an indexed catalog or be searched for by keyword or sequence similarity.
EVIDENCE is the first web-based tool providing pipelines for genotyping and recombination detection based on both sequence context and phylogenetic inference. Furthermore, EVIDENCE uses the most complete and regularly updated reference sequences to maintain high sensitivity and specificity, thereby accelerating genotype identification in clinical diagnosis and enhancing our understanding of EV evolution. EVIDENCE is available at http://symbiont.iis.sinica.edu.tw/evidence.
Materials and methods
Reference sequence sets
Nucleotide sequences of EV prototype strains (Additional file 1: Table S1) listed in the Picornaviridae Study Group website (http://www.picornaviridae.com/) were fetched from GenBank database. If the complete genome of a prototype strain was not available, we collected sequences of all the other available regions instead. For a genotype without a reference prototype assignment, the longest and/or earliest reported sequence was selected as the genotype's representative reference. Finally, 396 nucleotide sequences are selected to build the three hundred and eight prototype models (the genotyping reference set, GTRefSet). Furthermore, we expanded the prototype GTRefSet to the extended reference set (ExRefSet, Additional file 1: Table S2) by including the sequences that is highly similar (sequences identity >75%)  to the representative prototype sequence of the same genotype. GTRefSet was the core reference set for phylogenetic analysis and recombination analysis, and ExRefSet was used for automatic re-assignment of EV sequence genotypes.
Re-classification of all EV sequences in GenBank
Let rank(braw qx ),rank(cov qx ),rank(identity qx ) be the ascending sorting rank of BLAST raw score, alignment coverage, and identity between q and x in Cq, respectively. We build a set of the rank sums rS Cq where
Finally, the best mixed-ranking score (BMS) for the query q is defined as the largest rank sum r in rS Cq ,
Sequences are assigned to the genotype/species of their BMS references if the BMS is only reached by one reference sequence or multiple references with consistent genotype/prototype assignment. If the BMS of a query sequence is reached by two or more different genotypes in one species, the nomenclature assignment of this sequence is set to the species level. Query sequences are assigned as "unclassified" species with "unclassified" genotype if a BMS from a highly related sequence is not detected in ExRefSet (e.g., Cqis an empty set), or multiple references reach the BMS but are incongruent at the species level.
The rationales of using three parameters in the mixed-ranking score function are described below. The expectation value (E-value) given by BLAST is often used to indicate the significance of an alignment, and is often inferred to the homology/similarity relationship of the hit to the query sequences. Calculation of E-value is affected by the content of searching database and the length of the matching segment. For example, a small E-value may be granted to a short region in high sequence similarity, leading to false positive results if we use an E-value cut-off for selecting sequence matches. In contrast, the BLAST raw score is directly derived from alignment segments can eliminate this artifact . Thus, we use BLAST raw score instead. The BLAST raw score is also depended on the scoring parameters being used (i.e., reward for matching base and penalty for mismatching base/gap) [28, 29]. In our case, the scoring parameters are using the defaults in discontiguous megablast (dc-megablast), says, match = 2, mismatch = −3, gap open = −5, and gap extension = −2. This empirical setting is optimized for catching highly similar sequences with gap allowance. Then, we calculated the BLAST raw score for the query sequence q to the reference x in the reference set, and granted the rank score Braw qx to each q-x pair by the sorting order of the score.
Furthermore, we adopted the percentage of identity and the coverage of the alignment in the algorithm to address the importance of the conservation of the base components and the overall alignments of any two sequences in comparison. These two indexes are necessary for identifying the closest homologous sequences. The numeric values of the identity and coverage were transformed to ranked scores, in which we can sum up for the importance evaluated by three different scoring scheme in a normalized scale. Therefore, the mixed-ranking score is based on a measure of sequence similarity, adding weights on the importance of the quality of alignment.
Phylogenetic-based genotyping procedure
Executing the phylogenetic analysisrequires one query sequence and at least two reference strain prototype sequences. The genome regions of interest are cropped from the query sequence and concatenated for phylogenetic analysis. Multiple sequence alignments of the concatenated query and reference fragments are performed by Clustal Omega version 1.2.0 . The phylogenetic inference was done in by PhyML 3.0  (model GTR + G + I) with user defined bootstrap iterations. The circular phylogenetic tree topology is generated by jsPhyloSVG library .
Recombination event detection
We implemented Bootscan  and SimPlot  methods in EVIDENCE to detect possible recombination events and graphically present the results. Briefly, multiple sequence alignments between the query and reference genotype sequences (built into GTRefSet or custom uploaded) were performed using Clustal Omeg. The distance and similarity between sequences are calculated using DNADIST in PHYLIP (package version 3.5c)  and phylogenetic inference is analyzed using the Neighbor-joining (NJ) method with the Kimura two-parameters substitution model. Parameters, including sliding window size, step size, bootstrapping iterations, and the use of trimA1 (version 1.2)  to remove gaps, are adjustable. For each sliding window, the bootstrap value of the reference that was the first neighbor clustered to the query was used to derive the percentage of bootstrap support. Finally, profiles of query-reference closeness, evaluated by similarity or percentage of bootstrap support of each query-reference pair in each sliding window segment, were plotted along the EV genome. A clear crossing-over of two query-reference profiles with sharp slopes suggests a swap of best-fit reference and the presence of a nearby recombination breaking point.
EVIDENCE (http://symbiont.iis.sinica.edu.tw/evidence/) is constructed on an open-source Linux Arch (version), Nginx (version 0.12.4), and SQLAlchemy and SQLite relational database (version 220.127.116.11) structure. Graphical visualization was provided using Canvas and SVG library. Scripts for joining software packages to seamless pipelines were written in Perl and Python. The whole system is run in a virtual machine (CPUs of 2.27GHz, 8 cores, 16 GB RAM) located in the Institute of Information Science, Academia Sinica, Taiwan.
Results and discussion
The usage of EVIDENCE
EVIDENCE is a searchable database for updated and re-classified EV sequences and a workbench for EV genotyping and recombination detection (Figure 2).
A total of 54,790 up-to-date EV sequences were curated and built into an indexed category of virus species, genotype, and epidemiological annotations (host, continent, country, and the reported year). This database can be browsed through the hierarchical structure or searched by keywords (Figure 2A, panel 1) or sequence similarity (e.g., BLASTN, TBLASTN, or TBLASTX) (Figure 2A, panel 2). For both sequence search functions, the basic features extracted from GenBank records and reclassified genotypes are shown in table format. The search results, including the brief table of reported entries and the FASTA file of hits, are made available for download.
EVIDENCE matches query sequences to references in GTRefSet using the BLAST algorithm (Figure 2B, panel 1). To perform genotyping, the query nucleotide sequence in FASTA format is pasted or uploaded through the sequence input interface, and then submitted by clicking on 'scan region'. Scores for the query to each EV genomic region of each prototype are calculated and displayed in the results table. The table functions like a flexible input interface for the next phylogenetic analysis step. Clicking on the table column (the genome region) or table row (sequence title) will select the region or the references, respectively. The table is decreasingly sorted by the mixed-ranking score of a single column, or by the sum of mixed-ranking scores of all the selected columns (Figure 2B). The sorting order implies the relative fitness of each query-reference pair with respect to the region(s) of interest.
The phylogenetic analysis step calculates the relatedness between the query and the selected references. A tree topology is generated by PhyML with adjustable bootstrap iterations (default: 100 iterations). Phylogenetic analysis outputs, including the tree topology as a newick file and a png file, multiple sequence alignment, and the selected reference sequences, are made available for download (Figure 2B).
The basic principle for detecting potential recombination events is to segment the whole EV genome into small overlapping segments, in order to identify swaps of the most similar reference (plot similarity) or changes in the relatively closest neighbor (bootscanning) in the successive segments; similarity is derived from the sequence distance for each query-reference pair, and sequence neighbors are determined based on the percentage of bootstrap iterations supporting the reference as the closest neighbor to the query.
The recombination analysis takes a single sequence query, accompanied with three or more references to give adequate estimations. The reference strains are selected from those built into GTRefSet and/or from uploaded 'custom references'. The parameters for phylogenetic inference include sliding window (default: 200), step size (default: 20), bootstrapping (default: 100 iterations), and whether or not gaps in the multiple sequence alignment are trimmed (default: no trimming); all of these parameters are adjustable (Figure 2C, panel 1).
Plotting of the similarity or bootstrapping results is optional. The EV genome diagram is plotted on top to help visualize the location of recombination events. Dynamic figures are used to present the data and allow the user to zoom into/out of the plot. The value of each sliding window on the plot is shown via mouseover events, and the user may zoom in on a specific region of the plot by cropping the region through mouse dragging (Figure 2C, panel 2). The bootscanning/similarity plot, multiple alignment table, and bootstrap/similarity value table can be downloaded.
The choice of reference strains for recombination detection is crucial. Using inappropriate reference strains to identify the recombination region may eliminate the significance of the bootscanning result, and increase the noise of recombination breakpoint determination. In EVIDENCE, users can start the analysis from genotyping. References in GTRefSet are evaluated for the fitness of each genome region, which can help users select appropriate regions. The data, including the query sequence and at least three selected references, are redirected by clicking "push data to recombination detection" after the phylogenetic analysis step.
It is worth noting that using too many reference sequences in a recombination detection analysis may return an insignificant bootscanning plot or a messy similarity plot. We suggest that analysis should begin with less than ten references, and then the non-informative ones should be removed to improve the resolution.
Statistics of GTRefSet and ExtRefSet
Updating the classification of sequences using a BLAST-based method
To date, 54,790 EV sequences are deposited in GenBank, but the genotype assignment of these sequences in GenBank may be inaccurate. One possible reason for any inaccuracies is that the genotype of a sequence is assigned by individual researchers upon submission. However, differences in the typing methods used may lead to different conclusions. As mentioned above, serological assays and molecular typing may generate conflicting results if the examined strain has undergone a recombination event that changed the viral genome and the correlation of epitopes with the reference virus strain; such inconsistency is independent of the specificity of the antibody used to resolve strains at the genotype level. A second possible reason is that the nomenclature and the virus classification system are not synchronized between reference databases. Virus classification for NCBI taxonomy is based on that of the International Nucleotide Sequence Database Collaboration (INSDC) . However, classification of viruses is less stable than that of other organism kingdoms due to the rapidly evolving nature of these infectious life forms. Irregular updates have led to the inconsistent nomenclature and taxonomic classification used by GenBank and the International Committee on Virus Taxonomy (ICTV). Moreover, the updated viral classification scheme has not been applied to sequences previously deposited in the database.
Assigned by EVIDENCE
Sequences with Congruent Assignment (% *)
unclassified Enterovirus unclassified Rhinovirus
Region2 (sequence number)
5'UTR (1), 2C (2), 3D (8)
5'UTR (1), 3D (1)
2C (2), 3D (11)
5'UTR (69), 2C (2), 3D (3)
2C (2), 3D (5)
2C (1), 3D (5)
5'UTR (85), VP4-VP2 (2), 2C (3), 3D (4)
VP1 (1), 2C (1), 3D (1)
5'UTR (28), 2BC (2), 3D (24)
5'UTR (10), VP4 (3), VP1 (3), 3D (6)
Number (%) of Genotype Discrepancies
RIVM Genotyping Tool
(49 only typed to species level)
(4 seqs typed to species level, 2 mistyped)
(4 seqs typed to species level, 2 mistyped)
(1 seq typed to species level, 1 mistyped)
(1 seq typed to species level)
(5 successfully typed)
(1 seq typed to species level, 4 successfully typed)
A case study
Here, we demonstrate a genotyping and recombination pipeline using two coxsackievirus A16 (CV-A16) strains with distinct pathogenesis.
Coxsackievirus strains CV-A16 GD09/24 (GenBank accession KC117317) and GD09/119 (GenBank accession KC117318), exhibiting differing levels of clinical virulence, were isolated in Guangdong, China, in 2009 . The authors performed phylogenetic analysis with 28 CV-A16 homologous strains and one EV-A71 prototype strain to assign GD09/24 and GD09/119 to the CV-A16 genotype. To detect recombination events, the authors compared two novel CA-A16 strains with two EV-A71 strains and one CV-A16 prototype strain. Bootscanning results indicated that GD09/24 and GD09/119 underwent homologous recombination with EV-A71 in the P2 and P3 regions.
This demonstration shows that the genotyping and recombination pipeline in EVIDENCE can provide suitable candidates as references for recombination detection. Additionally, users can download all output files, and perform analyses using different reference sequences or genomic region(s) with a user-friendly interface.
Classical EV typing is largely dependent on serotyping methods. VP1 has been the subject of extensive research on account of the neutralization potency of its antiserum [23, 58, 72]. Typing specificity can be improved by using a panel of antibodies against VP1 and other viral proteins [73–77]. Thus, genotypes may be assigned through observing the response of several antibodies raised from epitopes in different genome regions. If a novel EV strain emerged from a recombinant event that joined epitopes of different parent strains, the serological phenotype may fail to reflect clinical virulence. As more EV sequences are reported, it is increasingly apparent that recombination occurs frequently within inter- or intra- genotypes. Moreover, each genomic region is subject to distinct selective pressures, and thus their evolution is independent of one another [78, 79]. Increased genetic diversity often leads to phenotypic variation, which is problematic for clinical therapy.
EVIDENCE can be used to perform EV typing based on sequence context. This tool disassociates reference prototypes into functional components of the virus genome, and performs analysis in a modularized manner. However, the correlation of individual genome regions with genome virulence remains unclear. We hope that EVIDENCE can be used to address this question and provide insights into EV evolution, as well as facilitate the diagnosis of clinical specimens to ensure appropriate treatment.
The authors wish to thank Dr. Min-Shi Lee and Dr. Pao-Yang Chen for valuable discussions on epidemiology of enterovirus and model of molecular tying, respectively. We also thank the editor and anonymous reviewers for their helpful advice.
The research was funded by Ministry of Science and Technology (MOST), Taiwan, for financially supporting this research through MOST 104-2319-B-400-002 to CAH, MOST 103-2311-B-001-033-MY3 to CYL, MOST 101-2321-B-001-043-MY2 and MOST 102-2811-B-001-046 to SHC.
The publication costs for this article were funded by Ministry of Science and Technology (MOST), Taiwan.
This article has been published as part of BMC Genomics Volume 16 Supplement 12, 2015: Joint 26th Genome Informatics Workshop and 14th International Conference on Bioinformatics: Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S12.
- Tapparel C, Siegrist F, Petty TJ, Kaiser L: Picornavirus and enterovirus diversity with associated human diseases. Infect Genet Evol. 2013, 14: 282-293.PubMedView ArticleGoogle Scholar
- Tsao KC, Huang CG, Huang YL, Chen FC, Huang PN, Huang YC, et al: Epidemiologic features and virus isolation of enteroviruses in Northern Taiwan during 2000-2008. J Virol Methods. 2010, 165 (2): 330-332.PubMedView ArticleGoogle Scholar
- Caro V, Guillot S, Delpeyroux F, Crainic R: Molecular strategy for 'serotyping' of human enteroviruses. J Gen Virol. 2001, 82 (Pt 1): 79-91.PubMedView ArticleGoogle Scholar
- 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 (8): 2698-2704.PubMedView ArticleGoogle 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 (5): 1288-1293.PubMedGoogle Scholar
- Zhou F, Kong F, Wang B, McPhie K, Gilbert GL, Dwyer DE: Molecular characterization of enterovirus 71 and coxsackievirus A16 using the 5' untranslated region and VP1 region. J Med Microbiol. 2011, 60 (Pt 3): 349-358.PubMedView ArticleGoogle Scholar
- Zhou F, Wang Q, Sintchenko V, Gilbert GL, O'Sullivan MV, Iredell JR, Dwyer DE: Use of the 5' Untranslated Region and VP1 Region to Examine the Molecular Diversity in Enterovirus B Species. J Med Microbiol. 2014, 63 (Pt 10): 1339-1355.PubMedView ArticleGoogle Scholar
- Ibrahim W, Boukhadra N, Nasri-Zoghlami D, Berthelot P, Omar S, Bourlet T, et al: Partial sequencing of the VP2 capsid gene for direct enterovirus genotyping in clinical specimens. Clin Microbiol Infect. 2014, 20 (9): O558-O565.PubMedView ArticleGoogle Scholar
- Ibrahim W, Ouerdani I, Pillet S, Aouni M, Pozzetto B, Harrath R: Direct typing of human enteroviruses from wastewater samples. J Virol Methods. 2014, 207: 215-219.PubMedView ArticleGoogle Scholar
- Perera D, Shimizu H, Yoshida H, Tu PV, Ishiko H, McMinn PC, Cardosa MJ: A comparison of the VP1, VP2, and VP4 regions for molecular typing of human enteroviruses. J Med Virol. 2010, 82 (4): 649-657.PubMedView ArticleGoogle Scholar
- Chan YF, Wee KL, Chiam CW, Khor CS, Chan SY, Amalina WM, Sam IC: Comparative genetic analysis of VP4, VP1 and 3D gene regions of enterovirus 71 and coxsackievirus A16 circulating in Malaysia between 1997-2008. Trop Biomed. 2012, 29 (3): 451-466.PubMedGoogle Scholar
- Casas I, Palacios GF, Trallero G, Cisterna D, Freire MC, Tenorio A: Molecular characterization of human enteroviruses in clinical samples: comparison between VP2, VP1, and RNA polymerase regions using RT nested PCR assays and direct sequencing of products. J Med Virol. 2001, 65 (1): 138-148.PubMedView ArticleGoogle Scholar
- Bochkov YA, Grindle K, Vang F, Evans MD, Gern JE: Improved molecular typing assay for rhinovirus species A, B, and C. J Clin Microbiol. 2014, 52 (7): 2461-2471.PubMedView ArticleGoogle Scholar
- Thoelen I, Moes E, Lemey P, Mostmans S, Wollants E, Lindberg AM, et al: Analysis of the serotype and genotype correlation of VP1 and the 5' noncoding region in an epidemiological survey of the human enterovirus B species. J Clin Microbiol. 2004, 42 (3): 963-971.PubMedView ArticleGoogle Scholar
- Chan YF, Sam IC, AbuBakar S: Phylogenetic designation of enterovirus 71 genotypes and subgenotypes using complete genome sequences. Infect Genet Evol. 2010, 10 (3): 404-412.PubMedView ArticleGoogle 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. Journal of General Virology. 2003, 84 (Pt 5): 1223-1235.PubMedView ArticleGoogle Scholar
- Hu L, Zhang Y, Hong M, Zhu S, Yan D, Wang D, et al: Phylogenetic evidence for multiple intertypic recombinations in enterovirus B81 strains isolated in Tibet, China. Sci Rep. 2014, 4: 6035-PubMedGoogle Scholar
- Tang J, Yoshida H, Ding Z, Tao Z, Zhang J, Tian B, et al: Molecular Epidemiology and Recombination of Human Enteroviruses from AFP surveillance in Yunnan, China from 2006 to 2010. Sci Rep. 2014, 4: 6058-PubMedView ArticleGoogle Scholar
- Yip CC, Lau SK, Lo JY, Chan KH, Woo PC, Yuen KY: Genetic characterization of EV71 isolates from 2004 to 2010 reveals predominance and persistent circulation of the newly proposed genotype D and recent emergence of a distinct lineage of subgenotype C2 in Hong Kong. Virol J. 2013, 10: 222-PubMedView ArticleGoogle Scholar
- Kroneman A, Vennema H, Deforche K, van der Avoort H, Penaranda S, Oberste MS, et al: An automated genotyping tool for enteroviruses and noroviruses. Journal of Clinical Virology. 2011, 51 (2): 121-125.PubMedView ArticleGoogle Scholar
- Rozanov M, Plikat U, Chappey C, Kochergin A, Tatusova T: A web-based genotyping resource for viral sequences. Nucleic Acids Res. 2004, 32 (Web Server issue): W654-W659.PubMedView ArticleGoogle Scholar
- King AM, Adams MJ, Lefkowitz EJ, Carstens EB: Virus taxonomy: classification and nomenclature of viruses: Ninth Report of the International Committee on Taxonomy of Viruses. Elsevier. 2012, 9:Google Scholar
- Oberste MS, Pallansch MA: Enterovirus molecular detection and typing. Reviews in Medical Microbiology. 2005, 16 (4): 163-171.View ArticleGoogle Scholar
- Adams MJ, Lefkowitz EJ, King AM, Carstens EB: Ratification vote on taxonomic proposals to the International Committee on Taxonomy of Viruses (2014). Archives of virology. 2014, 159 (10): 2831-2841.PubMedView ArticleGoogle Scholar
- Bukreyev AA, Chandran K, Dolnik O, Dye JM, Ebihara H, Leroy EM, et al: Discussions and decisions of the 2012-2014 International Committee on Taxonomy of Viruses (ICTV) Filoviridae Study Group, January 2012-June 2013. Archives of Virology. 2014, 159 (4): 821-830.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.PubMedView ArticleGoogle Scholar
- Rasko DA, Myers GS, Ravel J: Visualization of comparative genomic analyses by BLAST score ratio. BMC Bioinformatics. 2005, 6: 2-PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- Korf I, Yandell M, Bedell J: BLAST. 2003, O'Reilly Media, IncorporatedGoogle Scholar
- Blom N, Hansen J, Blaas D, Brunak S: Cleavage site analysis in picornaviral polyproteins: discovering cellular targets by neural networks. Protein Sci. 1996, 5 (11): 2203-2216.PubMedView ArticleGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011, 7: 539-PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59 (3): 307-321.PubMedView ArticleGoogle 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 (11): 1423-1425.PubMedView ArticleGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999, 73 (1): 152-160.PubMedGoogle Scholar
- Felsenstein J: PHYLIP 3.5 (phylogeny inference package). Department of Genetics, University of Washington, Seattle. 1993Google Scholar
- Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T: trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009, 25 (15): 1972-1973.PubMedView ArticleGoogle Scholar
- Federhen S: The NCBI Taxonomy database. Nucleic Acids Res. 2012, 40 (Database issue): D136-D143.PubMedView ArticleGoogle Scholar
- Tu PV, Thao NT, Perera D, Huu TK, Tien NT, Thuong TC, et al: Epidemiologic and virologic investigation of hand, foot, and mouth disease, southern Vietnam, 2005. Emerg Infecti Dis. 2007, 13 (11): 1733-1741.View ArticleGoogle Scholar
- Sanders SA, Herrero LJ, McPhie K, Chow SS, Craig ME, Dwyer DE, et al: Molecular epidemiology of enterovirus 71 over two decades in an Australian urban community. Arch Virol. 2006, 151 (5): 1003-1013.PubMedView ArticleGoogle Scholar
- Nasri D, Bouslama L, Omar S, Saoudin H, Bourlet T, Aouni M, et al: Typing of human enterovirus by partial sequencing of VP2. J Clin Microbiol. 2007, 45 (8): 2370-2379.PubMedView ArticleGoogle Scholar
- Aw TG, Gin KY: Environmental surveillance and molecular characterization of human enteric viruses in tropical urban wastewaters. Journal of Applied Microbiology. 2010, 109 (2): 716-730.PubMedGoogle Scholar
- Zhong JY, Zhu B, Hua L, Wang CB, Kuang L, Xie JH, Chen Y: [Complete genomic sequence analysis on human enterovirus 71 strain in Guangzhou, in 2008 and 2010]. Zhonghua Liu Xing Bing Xue Za Zhi. 2011, 32 (7): 700-704.PubMedGoogle Scholar
- Bessaud M, Pillet S, Ibrahim W, Joffret ML, Pozzetto B, Delpeyroux F, Gouandjika-Vasilache I: Molecular characterization of human enteroviruses in the Central African Republic: uncovering wide diversity and identification of a new human enterovirus A71 genogroup. J Clin Microbiol. 2012, 50 (5): 1650-1658.PubMedView ArticleGoogle Scholar
- Sun Q, Zhang Y, Zhu S, Tian H, Huang G, Cui H, et al: Transmission of human enterovirus 85 recombinants containing new unknown serotype HEV-B donor sequences in Xinjiang Uighur autonomous region, China. PLoS One. 2013, 8 (1): e55480-PubMedView ArticleGoogle Scholar
- Lukashev AN, Shumilina EY, Belalov IS, Ivanova OE, Eremeeva TP, Reznik VI, et al: Recombination strategies and evolutionary dynamics of the Human enterovirus A global gene pool. J Gen Virol. 2014, 95 (Pt 4): 868-873.PubMedView ArticleGoogle Scholar
- Ni H, Yi B, Yin J, Fang T, He T, Du Y, et al: Epidemiological and etiological characteristics of hand, foot, and mouth disease in Ningbo, China, 2008-2011. J Clin Virol. 2012, 54 (4): 342-348.PubMedView ArticleGoogle Scholar
- Oberste MS, Maher K, Michele SM, Belliot G, Uddin M, Pallansch MA: Enteroviruses 76, 89, 90 and 91 represent a novel group within the species Human enterovirus A. J Gen Virol. 2005, 86 (Pt 2): 445-451.PubMedView ArticleGoogle Scholar
- Laitinen OH, Honkanen H, Pakkanen O, Oikarinen S, Hankaniemi MM, Huhtala H, et al: Coxsackievirus B1 is associated with induction of beta-cell autoimmunity that portends type 1 diabetes. Diabetes. 2014, 63 (2): 446-455.PubMedView ArticleGoogle Scholar
- Chiang PS, Huang ML, Luo ST, Lin TY, Tsao KC, Lee MS: Comparing molecular methods for early detection and serotyping of enteroviruses in throat swabs of pediatric patients. PLoS One. 2012, 7 (10): e48269-PubMedView ArticleGoogle Scholar
- Lu QB, Zhang XA, Wo Y, Xu HM, Li XJ, Wang XJ, et al: Circulation of Coxsackievirus A10 and A6 in hand-foot-mouth disease in China, 2009-2011. PLoS One. 2012, 7 (12): e52073-PubMedView ArticleGoogle 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. Scandinavian journal of infectious diseases. 1999, 31 (4): 331-335.PubMedView ArticleGoogle Scholar
- Craig ME, Howard NJ, Silink M, Rawlinson WD: Reduced frequency of HLA DRB1*03-DQB1*02 in children with type 1 diabetes associated with enterovirus RNA. J Infect Dis. 2003, 187 (10): 1562-1570.PubMedView ArticleGoogle Scholar
- Garcia J, Espejo V, Nelson M, Sovero M, Villaran MV, Gomez J, et al: Human rhinoviruses and enteroviruses in influenza-like illness in Latin America. Virol J. 2013, 10: 305-PubMedView ArticleGoogle Scholar
- Sadeuh-Mba SA, Bessaud M, Joffret ML, Endegue Zanga MC, Balanant J, Mpoudi Ngole E, et al: Characterization of Enteroviruses from non-human primates in cameroon revealed virus types widespread in humans along with candidate new types and species. PLoS Negl Trop Dis. 2014, 8 (7): e3052-PubMedView ArticleGoogle 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 (2): 855-867.PubMedView ArticleGoogle Scholar
- Siafakas N, Papaventsis D, Levidiotou-Stefanou S, Vamvakopoulos NC, Markoulatos P: Classification and structure of echovirus 5'-UTR sequences. Virus Genes. 2005, 31 (3): 293-306.PubMedView ArticleGoogle Scholar
- Oberste MS, Maher K, Kilpatrick DR, Pallansch MA: Molecular evolution of the human enteroviruses: correlation of serotype with VP1 sequence and application to picornavirus classification. J Virol. 1999, 73 (3): 1941-1948.PubMedGoogle Scholar
- Oberste MS, Maher K, Flemister MR, Marchetti G, Kilpatrick DR, Pallansch MA: Comparison of classic and molecular approaches for the identification of untypeable enteroviruses. J Clin Microbiol. 2000, 38 (3): 1170-1174.PubMedGoogle Scholar
- Wang JR, Tsai HP, Huang SW, Kuo PH, Kiang D, Liu CC: Laboratory diagnosis and genetic analysis of an echovirus 30-associated outbreak of aseptic meningitis in Taiwan in 2001. J Clin Microbiol. 2002, 40 (12): 4439-4444.PubMedView ArticleGoogle Scholar
- Chen Y, Li C, He D, Cheng T, Ge S, Shih JW, et al: Antigenic analysis of divergent genotypes human Enterovirus 71 viruses by a panel of neutralizing monoclonal antibodies: current genotyping of EV71 does not reflect their antigenicity. Vaccine. 2013, 31 (2): 425-430.PubMedView ArticleGoogle Scholar
- Han JF, Yu N, Pan YX, He SJ, Xu LJ, Cao RY, et al: Phenotypic and genomic characterization of human coxsackievirus A16 strains with distinct virulence in mice. Virus Research. 2014, 179: 212-219.PubMedView ArticleGoogle Scholar
- Junttila N, Leveque N, Magnius LO, Kabue JP, Muyembe-Tamfum JJ, Maslin J, et al: Complete coding regions of the prototypes enterovirus B93 and C95: phylogenetic analyses of the P1 and P3 regions of EV-B and EV-C strains. J Med Virol. 2015, 87 (3): 485-497.PubMedView ArticleGoogle Scholar
- Riquet FB, Blanchard C, Jegouic S, Balanant J, Guillot S, Vibet MA, et al: Impact of exogenous sequences on the characteristics of an epidemic type 2 recombinant vaccine-derived poliovirus. J Virol. 2008, 82 (17): 8927-8932.PubMedView ArticleGoogle Scholar
- Oberste MS: Comparative genomics of the coxsackie B viruses and related enteroviruses. Current topics in microbiology and immunology. 2008, 323: 33-47.PubMedGoogle Scholar
- Arita M, Shimizu H, Nagata N, Ami Y, Suzaki Y, Sata T, et al: Temperature-sensitive mutants of enterovirus 71 show attenuation in cynomolgus monkeys. J Gen Virol. 2005, 86 (Pt 5): 1391-1401.PubMedView ArticleGoogle Scholar
- Weinzierl AO, Rudolf D, Maurer D, Wernet D, Rammensee HG, Stevanovic S, Klingel K: Identification of HLA-A*01- and HLA-A*02-restricted CD8+ T-cell epitopes shared among group B enteroviruses. J Gen Virol. 2008, 89 (Pt 9): 2090-2097.PubMedView ArticleGoogle Scholar
- Mao Q, Wang Y, Yao X, Bian L, Wu X, Xu M, Liang Z: Coxsackievirus A16: epidemiology, diagnosis, and vaccine. Hum Vaccin Immunother. 2014, 10 (2): 360-367.PubMedView ArticleGoogle Scholar
- Zhang C, Zhu R, Yang Y, Chi Y, Yin J, Tang X, et al: Phylogenetic analysis of the major causative agents of hand, foot and mouth disease in Suzhou city, Jiangsu province, China, in 2012-2013. Emerg Microbes Infect. 2015, 4 (2): e12-PubMedView ArticleGoogle Scholar
- Tan X, Li L, Zhang B, Jorba J, Su X, Ji T, et al: Molecular epidemiology of coxsackievirus A6 associated with outbreaks of hand, foot, and mouth disease in Tianjin, China, in 2013. Arch Virol. 2015, 160 (4): 1097-1104.PubMedView ArticleGoogle Scholar
- Xu M, Su L, Cao L, Zhong H, Dong N, Xu J: Enterovirus genotypes causing hand foot and mouth disease in Shanghai, China: a molecular epidemiological analysis. BMC Infectious Diseases. 2013, 13: 489-PubMedView ArticleGoogle Scholar
- Foo DG, Alonso S, Phoon MC, Ramachandran NP, Chow VT, Poh CL: Identification of neutralizing linear epitopes from the VP1 capsid protein of Enterovirus 71 using synthetic peptides. Virus Res. 2007, 125 (1): 61-68.PubMedView ArticleGoogle Scholar
- Kiener TK, Jia Q, Lim XF, He F, Meng T, Chow VT, Kwang J: Characterization and specificity of the linear epitope of the enterovirus 71 VP2 protein. Virol J. 2012, 9: 55-PubMedView ArticleGoogle Scholar
- Buttinelli G, Donati V, Ruggeri FM, Joki-Korpela P, Hyypia T, Fiore L: Antigenic sites of coxsackie A9 virus inducing neutralizing monoclonal antibodies protective in mice. Virology. 2003, 312 (1): 74-83.PubMedView ArticleGoogle Scholar
- Abed Y, Wolf D, Dagan R, Boivin G: Development of a serological assay based on a synthetic peptide selected from the VP0 capsid protein for detection of human parechoviruses. J Clin Microbiol. 2007, 45 (6): 2037-2039.PubMedView ArticleGoogle Scholar
- Dang S, Gao N, Li Y, Li M, Wang X, Jia X, et al: Dominant CD4-dependent RNA-dependent RNA polymerase-specific T-cell responses in children acutely infected with human enterovirus 71 and healthy adult controls. Immunology. 2014, 142 (1): 89-100.PubMedView ArticleGoogle Scholar
- Varela-Calvino R, Skowera A, Arif S, Peakman M: Identification of a naturally processed cytotoxic CD8 T-cell epitope of coxsackievirus B4, presented by HLA-A2.1 and located in the PEVKEK region of the P2C nonstructural protein. J Virol. 2004, 78 (24): 13399-13408.PubMedView ArticleGoogle Scholar
- Lukashev AN: Role of recombination in evolution of enteroviruses. Reviews in medical virology. 2005, 15 (3): 157-167.PubMedView ArticleGoogle Scholar
- Lukashev AN, Lashkevich VA, Ivanova OE, Koroleva GA, Hinkkanen AE, Ilonen J: Recombination in circulating enteroviruses. J Virol. 2003, 77 (19): 10423-10431.PubMedView ArticleGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.