Skip to main content

Genetic characterization of a Coxsackie A9 virus associated with aseptic meningitis in Alberta, Canada in 2010



An unusually high incidence of aseptic meningitis caused by enteroviruses was noted in Alberta, Canada between March and October 2010. Sequence based typing was performed on the enterovirus positive samples to gain a better understanding of the molecular characteristics of the Coxsackie A9 (CVA-9) strain responsible for most cases in this outbreak.


Molecular typing was performed by amplification and sequencing of the VP2 region. The genomic sequence of one of the 2010 outbreak isolates was compared to a CVA-9 isolate from 2003 and the prototype sequence to study genetic drift and recombination.


Of the 4323 samples tested, 213 were positive for enteroviruses (4.93%). The majority of the positives were detected in CSF samples (n = 157, 73.71%) and 81.94% of the sequenced isolates were typed as CVA-9. The sequenced CVA-9 positives were predominantly (94.16%) detected in patients ranging in age from 15 to 29 years and the peak months for detection were between March and October. Full genome sequence comparisons revealed that the CVA-9 viruses isolated in Alberta in 2003 and 2010 were highly homologous to the prototype CVA-9 in the structural VP1, VP2 and VP3 regions but divergent in the VP4, non-structural and non-coding regions.


The increase in cases of aseptic meningitis was associated with enterovirus CVA-9. Sequence divergence between the prototype strain of CVA-9 and the Alberta isolates suggests genetic drifting and/or recombination events, however the sequence was conserved in the antigenic regions determined by the VP1, VP2 and VP3 genes. These results suggest that the increase in CVA-9 cases likely did not result from the emergence of a radically different immune escape mutant.


The enterovirus genome is comprised of a single open reading frame flanked by the 5’ and 3’ untranslated regions (UTRs), and the encoded polyprotein is cleaved to produce the structural and nonstructural proteins [1]. Although most infections are asymptomatic, non-polio enteroviruses are the most common infectious cause of aseptic meningitis [1]. Several outbreaks resulting from different enterovirus serotypes have been described including but not restricted to Coxsackievirus A9 (CVA-9) [2], EV71 [3, 4], EV68 [5], Coxsackievirus A24 [6], echovirus 18 [7], echovirus 30 [8, 9], and Coxsackievirus A16 [10].

Non-polio enteroviruses were traditionally classified into serotypes, based on a neutralization assay and included 64 classical isolates consisting of Coxsackieviruses A, Coxsackieviruses B, echoviruses and several numbered enteroviruses [11]. Genetic studies led to a new classification scheme, grouping the enteroviruses into species A, B, C and D, although serotype determination is still widely used, including for epidemiological purposes. Typing methods based on sequencing have been previously described in the VP1 [1214], VP2 [15, 16], and VP4 [17] regions. Molecular typing methods based on the P1 genomic region sequences, have been found to yield equivalent results to neutralizing assays, which they have essentially replaced [1, 14]. Molecular typing has also led to the characterization of new enterovirus serotypes [4, 1820] and provides information on recombination between enteroviruses. This is a common phenomenon which almost always occurs between viruses belonging to the same species [1, 21, 22]. Serotype determination remains important for molecular epidemiology, in part because the capsid contains the neutralization epitopes, and a rise or fall in incidence of a serotype can be linked to serotype-specific humoral immunity within the population [21]. The capsid also determines the cellular receptors used by the virus [23] and is thus an important determinant of viral pathogenesis, although non-structural proteins also contribute to the pathogenesis of the virus [24].

An unusually high incidence of aseptic meningitis was noted in Alberta, Canada between March and October 2010 [25]. A high proportion of CSF samples submitted to the Provincial Laboratory For Public Health of Alberta were positive for enteroviruses and serotyping by molecular methods revealed that majority of these enteroviruses belonged to the CVA-9 serotype. The primary goal of this study was to provide a full genetic characterization of the CVA-9 isolates responsible for this outbreak. Secondary goals were to comment on the methods of molecular serotyping for the diagnostic laboratory and on the use of long RT-PCR as a convenient method to obtain the near full-length genomic sequence of enteroviruses, especially when recombination events are involved in the genesis of the viral genome. The genome sequence of these isolates was determined and compared to the sequence of a CVA-9 isolate from Alberta in 2003 and to the prototype CVA-9 sequence (strain: Griggs; D00627) [26]. Serotyping by sequencing within the VP2 region was found to be more reliable than within the VP4 region. A comparison of the different genes revealed a higher nucleotide and amino acid conservation in the structural regions, and analysis of the sequence of the non-structural region pointed to recombination events in the genesis of the 2010 isolate.


Detection and serotyping

Of the 4323 samples tested, 213 were positive for enteroviruses (4.93%). The majority of the positives were CSF samples (n = 157, 73.71%) followed by specimens of respiratory origin (n = 17, 7.98%), stools (n = 10, 4.70%) and plasma (n = 8, 3.76%). The rest of this analysis concentrates on CSF samples. A total of 72 positive CSF samples were randomly selected to represent different age groups, geographic locations and disease severity for serotype determination based on phylogenetic clustering of the partial sequence of the VP2 region with the prototype sequences. More than 85% of samples could be typed without the need for culture or nested PCR (data not shown). This RT-PCR detects both enteroviruses and rhinoviruses, and allows for differentiation based on the length of the amplicon.

The most common serotype detected was CVA-9 (n = 59, 81.94%) and it was clearly the predominant serotype in the population. Closer examination of age distribution revealed that the majority (n = 56, 94.92%) of the CVA-9 positives were from patients ranging in age from 15 to 29 years, only three (5.08%) CVA-9 positives were detected in children less than a year old. The other frequently detected serotypes were CVB-3 (n = 3) and CVB-4 (n = 2). Figure 1A and 1B show the distribution of the different serotypes detected from infants less than one year old and from patients older than one year. A phylogenetic tree including representatives of the different enterovirus serotypes detected and some representative prototype strains in the VP2 region is indicated in Figure 2. The VP2 sequence from all the sequenced isolates was identical, showing that the outbreak was associated with a clonal or nearly clonal CVA-9 population.

Figure 1

Distribution of the different serotypes detected in the sequenced samples. A indicates serotype distribution in samples from patients less than one year and B from patients over one year old.

Figure 2

Phylogenetic tree including representatives of the different serotypes detected and corresponding prototype sequences (in italics) . The alignment used for the phylogenetic tree included 173 bases of the VP2 region from the clinical specimens, the same region from the prototype strains is also included. The Genbank accession number for the prototype strains is included in legend for Table 2. Patient samples are identified using a laboratory number followed by the enterovirus type.

Age and gender distribution

Of the 59 samples positive for CVA-9, roughly equal numbers were detected in males (n = 30) and females (n = 29). Figure 3 shows the age distribution and percentage of CSF samples positive for CVA-9 within the different age groups tested. The data indicates that the majority of the CVA-9 positives were detected in patients between 15 to 30 years old.

Figure 3

Age distribution of patients positive for enterovirus CVA-9. The bars represent the total number of samples positive for CVA-9 and the broken line indicates the percentage of samples positive for CVA-9 among the serotyped samples in each age group. The solid line indicates the percentage of samples positive for all enteroviruses in each age group based on the number of CSF samples tested within that age group.


Samples from January 2010 to December 2010 were used in this study; peak months for the detection of CVA-9 positives in CSF samples were between March and October.

Comparison of near full length genomes

The genomic sequences of the CVA-9_2003 and CVA-9_2010 have been deposited in GenBank under accession number JQ837913 and JQ837914, respectively. Comparison of the nucleotide and amino acid sequences in each of the gene segments from CVA-9_2003 and CVA-9_2010 to the prototype CVA-9 sequence are shown in Table 1. Nucleotide identity in the 5’ untranslated region (5’ UTR) between CVA-9_2003 and CVA-9_2010 was higher (96.5%) than with the prototype CVA-9 sequence which was around 85%. The nucleotide and amino acid identity was high (around 95%) between CVA-9_2003 and CVA-9_2010 samples in the structural genes (VP4, VP2, VP3 and VP1). The nucleotide identity in the structural genes for both these samples was lower (around 80%) when compared to the prototype CVA-9 sequence but amino acid identity was high (around 90%) indicating that the nucleotide changes consisted mostly of synonymous mutations. CVA-9_2010 showed some sequence diversity from the CVA-9_2003 isolate in the 3’ non-structural domain. This domain encodes the 2A protease, 2B protein, 2C helicase, 3A protein, 3B protein (VPg), 3C protease and 3D protein (polymerase). The nucleotide sequence identity in the 2A protease and 2B protein regions was 96.15% and 93.60% respectively; identity in the 2C helicase, 3A protein and 3D protein was around 80%. The lowest nucleotide identity was observed in the 3B protein at 74.63% and 3C protease at 77.78%, however amino acid conservation in these genes was higher as indicated in Table 1. Comparison of these gene segments to the prototype CVA-9 strain showed a higher percent of amino acid identity as compared to nucleotide identity.

Table 1 Percent identity comparison of nucleotide and amino acid sequences for isolates from 2003 (JQ837913), 2010 (JQ837914) and prototype CVA-9 sequence (D00627)

Phylogenetic analysis comparing the nucleotide sequences of all the species B prototype sequences to CVA-9_2003 and CVA-9_2010 shows that in the VP2 region, these sequences cluster with the prototype CVA-9 sequence with a high bootstrap value; a similar pattern was observed for the VP1 and VP3 regions. However, in the VP4 region, the closest neighbour to CVA-9_2003 and CVA-9_2010 sequences is enterovirus 101 (EV-101) and there is significant divergence from the prototype CVA-9 sequence at the nucleotide level (data not shown). Table 2 lists the nearest neighbour based on phylogenetic analysis for the structural and non-structural gene segments, compared to the gene sequences from prototype strains. As indicated in the table the regions 2A, 2B, and 3B from the CVA-9_2003 and CVA-9_2010 cluster together, however, the 2C, 3A, 3C and 3D regions from CVA-9_2003 and CVA-9_2010 cluster with different enterovirus species B prototype sequences indicating base pair changes or recombination in these regions. The probable mosaic nature of the genome of CVA-9_2003 and CVA-9_2010 was further explored by performing Simplot analysis. Figure 4A compares the sequence of CVA-9_2003 to that of the enteroviruses identified in Table 2 as the most homologous in each gene. All the viruses are highly homologous within the 5’NC, and the higher score of EV-74 may or may not reflect the origin of the 5’UTR segment of the Alberta isolates; nevertheless, sequencing the 5’UTR would not allow for accurate serotyping. It shows that in the VP2-VP3-VP1 region the CVA-9_2003 is most homologous to the prototype CVA-9; it is speculated that the differences are attributable to genetic drifting over the years. It is interesting to note that in the VP4 region the highest homology is with EV-101, suggesting that recombination with EV-101 (or a very similar virus) might have occurred. This suggests that relying only on the VP4 sequence would not allow for accurate serotyping. Simplot analysis also shows that in the non-structural coding region, recombination with other viruses of species B has occurred, although the dominance of individual viruses is not clear enough to attribute the origin of the segments with certainty. Similarly, Figure 4B compares the sequence of CVA-9_2010 to the enteroviruses identified in Table 2. Similar comments can be made regarding the sequence of CVA-9_2010, although in the non-structural domain, identifying CVB-6, EV-86 and EV-100 as contributing to the recombination events appears more convincing. Tables 1 and 2 show that there is a high nucleotide identity between CVA-9_2003 and CVA-9_2010 in the region upstream of that coding for protein 2B. It can be speculated that CVA-9_2010 arose from recombination events between CVA-9_2003, CVB-6, EV-86 and EV-100. This hypothesis is examined by Simplot analysis in Figure 4C. This analysis suggests that CVA-9 (2010 and 2003) are highly homologous from the 5’NC to the 2B domain; recombination with other viruses clearly occurred beyond this point, with contributions from EV-86 and EV-100, and likely from CVB-6.

Table 2 Closest neighbour to the 2003 and 2010 samples based on phylogenetic analysis
Figure 4

Simplot analysis of full genome sequences of CVA-9_2003, CVA-9_2010 and related species B sequences. Genomic regions are indicated at the bottom of each panel based on numbering of the prototype CVA-9 (D00627) 5’UTR:1–722; VP4:723–929; VP2:930–1712; VP3:1713–2426; VP1:2427–3332; 2A:3333–3773; 2B:3774–4069; 2C:4070–5057; 3A:5058–5323; 3B:5324–5390; 3C:5391–5938; 3D:5939–7324; 3’UTR:7325 - >7403. The query sequences used are CVA-9_2003 in panel 4A and CVA-9_2010 in panels 4B and 4C. Arrows in panels 4A and 4B indicate recombination with EV-101 in the VP4 region.


Recombination between enteroviruses is quite frequent, and occurs typically within the same species [1, 27]. It was initially thought that the 5’UTR and P1 regions move independently during recombination, and the frequency of recombination was higher in the non-structural coding region [22]. This lead to the conclusion that serotyping by molecular methods involving sequencing of the P1 region should be equivalent to sequencing segments of this region (VP1,2,3 or 4). Several studies have confirmed the finding that sequencing of the VP1 or VP2 regions yields similar typing results compared to sequencing the entire P1 region [1, 14, 28]. Typing using the VP4 region has yielded conflicting results including studies reporting both consistent and inconsistent results with classical serotyping [28, 29]; Perera et al. have noted that the VP4 sequence is unreliable for serotyping of viruses in species B and C [29]. This inconsistency seems to result from recombination events involving the VP4 region [22]. Intuitively, one may still continue to regard the segment coding for the external domains of the capsid to recombine as a block since the VP4 domain remains completely inside the capsid and is inaccessible to antibodies [1].

The amplicon used for serotype determination includes a segment from the 5’UTR, the VP4 region and a segment of the VP2 region; but given the possibility of recombinations outside the VP2-VP3-VP1 region, only the sequence within VP2 was used for the typing of samples in this study. It was noted that if serotype assignment is made based on Blast comparisons with the GenBank database sequences, there is a potential for inaccurate serotype assignment since contemporary isolates submitted into GenBank may have been serotyped based on the VP4 sequence, which can reassort independently of the VP2-VP3-VP1 segment. In this study, the serotype assignment was achieved by constructing a phylogenetic tree with the partial VP2 sequence obtained by sequencing the amplicon and the corresponding VP2 sequences of the original prototype strains [1]. In spite of the genetic drift evident in the sequences of contemporary viruses, this method provided unambiguous serotype assignment. Instances of recombination within the VP2-VP3-VP1 segment have only been rarely reported. Al-Hello and colleagues reported on an isolate that was identified as echovirus 11 by genetic analysis, but could be neutralized by antisera specific for echovirus 11 or CVA-9. A peptide scan analysis confirmed the presence of epitopes recognized by both antisera, whereas a Simplot analysis failed to reveal a recombination event within the capsid [30]. Should reports of isolates with similar properties become more frequent, the very concept of enterovirus serotype might become untenable. Notwithstanding these considerations, at this point in time it remains valuable to serotype isolates, because the serotype is determined by the neutralizing epitopes, and remains important for the molecular epidemiology of these viruses. In addition, the capsid determines the choice of the cellular receptor, consequently the serotype is an important factor to determine the pathogenicity of an isolate [23]; however, the remainder of the genome also contributes to pathogenicity [24] and complete genome sequencing is required for full characterisation of an isolate.

During the 2010 outbreak, a different serotype distribution was noted among infants (< 1year) and older individuals, this difference in distribution of serotypes between age groups has been noted before [31, 32]. The serotypes found among infants included in our study have been typically associated with neonatal disease, for example, Coxsackie B3 and B4, echovirus 9, and EV-71. It was interesting to observe one isolate with the VP2 region most similar to the newly described EV-97 [20]. Among the rest of the population, CVA-9 was clearly the predominant serotype. Phylogenetic analysis of the VP2 region from the CVA-9 isolates showed that all the isolates had the same sequences and therefore originated from the same strain.

Comparison of the capsid regions showed the CVA-9_2003 and CVA-9_2010 isolates to be highly homologous and this was corroborated by Simplot anaysis; consequently it is speculated that the sudden increase in detection of CVA-9 in CSF samples is attributable to an increased number of individuals susceptible to CVA-9 in the population rather than to the emergence of an immune escape mutant in 2010. For an infectious agent to cause an outbreak, it is necessary that the herd immunity of the population drops below a threshold that is determined by the basic reproduction number R0[33]. For CVA-9 to re-enter the population of Alberta in 2010, either the new isolate was an immune escape mutant against which the population had no prior immunity, or the herd immunity had declined below the threshold. Because of the very high homology in the capsid sequences between the 2003 and 2010 isolates, it does not seem likely that the 2010 isolate was an immune escape mutant.

A comparison with the sequence of the prototype CVA-9 shows that the two Alberta isolates inherited their VP2-VP3-VP1 segment from the prototype, but the VP4 segment shows sequence divergence. For all the other protein coding segments, phylogenetic trees of the nucleotide sequence were constructed to compare the Alberta strains with the prototype strains within species B, and Simplot analysis were performed. Two important observations can be made: firstly, the two Alberta isolates are highly homologous in most regions, with few changes that can be largely attributed to sequence drift; however, recombination has likely occurred in the 2C, 3A, 3C, and 3D regions. Secondly, the Alberta isolates are very different from the prototype CVA-9 in the VP4 segment (supporting the concept that VP4 can be a site for recombination). The 5’NC, structural and non-structural regions appear to be components of a mosaic genome modified from the prototype CVA-9 by recombination and drifting. Similar observations have been made on other contemporary CVA-9 strains [34], and indeed in other contemporary enteroviruses within species B [22, 27]. How much the pathogenicity of these CVA-9 viruses has been modified by the non-structural genes, compared to the prototype strain, is unclear [22].


In summary the sudden increase in cases of aseptic meningitis in Alberta in 2010 was associated with enterovirus CVA-9. The capsid region was highly homologous to the capsid of a 2003 isolate, suggesting that the infections were not the result of the emergence of an immune escape mutant. We thus speculate that the increase in the number of infections may have resulted from a decline in the level of herd immunity against this virus to a level where the virus was able to penetrate the population. When compared to the prototype strain of CVA-9, the Alberta isolates displayed signs of multiple recombination events in addition to genetic drifting.


Screening for enteroviruses

Specimens submitted to the Provincial Laboratory for Public Health (ProvLab) for enterovirus testing from January 1 to December 31, 2010 were included in this study. A total of 4323 samples from patients ranging in age from one day old to 97 years of age were tested; the most common specimen types tested included CSF (n = 2687, 62.16%), plasma (n = 497, 11.50%), nasopharyngeal and throat swabs (n = 213, 4.93%) and stools (n = 103, 2.40%). Viral RNA was extracted using the easyMAG® automated extractor (BioMérieux, Durham, NC, USA) and the extracted nucleic acid was screened for enteroviruses using a previously published NASBA assay [35].

Serotyping of enteroviruses

Enteroviruses were serotyped using one-step RT-PCR to amplify the partial 5’untranslated region (293 bp), the VP4 region (207 bp) and partial VP2 region (250 bp) using previously described primers 1-EV/RV and 2-EV/RV from RNA extracted directly from the specimen [36]. Sequencing was performed without the need for nested amplification. Amplification was performed using the One-step RT-PCR kit from Qiagen (Ontario, Canada) using 10 μl of 5X buffer, 10 μl of Q solution, 2 μl of 10 mM dNTPs, 2 μl of enzyme, 0.125 μl of RNaseOUT (Life Technologies, Ontario, Canada), 0.8 μM of primers and 5 μl of template nucleic acid in a total volume of 50 μl. The reverse transcription step was performed at 50°C for 30 mins, followed by enzyme activation at 95°C for 15 mins. The amplification protocol included 45 cycles of denaturation at 95°C for 30 seconds, followed by annealing at 55°C for 1.5 minutes and amplification at 72°C for 60 seconds. A final extension step was performed for 10 minutes at 72°C followed by cooling. Amplified products were sequenced in both directions on the ABI PRISM 3130-Avant Genetic Analyzer (Applied Biosystems (ABI), Foster City, CA).

Complete genome sequencing of CVA-9

Two enterovirus positive samples from 2003 (CVA-9_2003) and 2010 (CVA-9_2010) with a high viral load as estimated by the Ct value, that exhibited strong growth in primary Rhesus Monkey Kidney cells from Diagnostic Hybrids (Ohio, USA) were used. The near-complete genome of these viruses was amplified by long RT-PCR as previously described [37]. The amplicons were sequenced by genome walking and contig assembly was performed using seqscape v2.5 (ABI). The genomic sequences of the CVA-9_2003 and CVA-9_2010 have been deposited in GenBank under accession numbers JQ837913 and JQ837914, respectively.

Sequence analysis

Sequences were aligned using ClustalX (Version 1.81, March 2000; and phylogenetic analysis was conducted using Treecon [38]. Distance estimation was performed using the Jukes and Cantor distance correction (1969), topology inference was performed using the neighbour joining method, and bootstraping was done using 1000 replicates and the tree was re-rooted at the internode. Simplot analysis were performed using alignments done with CLustalX and the Simplot for Windows v3.5.1 program [39].



Cerebrospinal fluid


Coxsackie virus A


Coxsackie virus B




Untranslated region


Nucleic acid sequence based amplification


Reverse-transcription polymerase chain reaction


Applied biosystems


Coxsackie virus A9


  1. 1.

    Pallansch M, Roos R: Enteroviruses: Polioviruses, Coxsackieviruses, Echoviruses, and Newer Enteroviruses. In Fields Virology. 5th edition. Edited by: Knipe DM, Howley PM. Philadelphia, PA: Lippincott Williams Wilkins; 2011:839-894.

    Google Scholar 

  2. 2.

    Cui A, Yu D, Zhu Z, Meng L, Li H, Liu J, Liu G, Mao N, Xu W: An outbreak of aseptic meningitis caused by coxsackievirus A9 in Gansu, the People’s Republic of China. Virol J 2010, 7: 72. 10.1186/1743-422X-7-72

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Solomon T, Lewthwaite P, Perera D, Cardosa MJ, McMinn P, Ooi MH: Virology, epidemiology, pathogenesis, and control of enterovirus 71. Lancet Infect Dis 2010, 10: 778-790. 10.1016/S1473-3099(10)70194-8

    PubMed  Article  Google Scholar 

  4. 4.

    Stanway G, Brown F, Christian P: Picornaviridae. In Virus Taxonomy-classification and nomenclature of viruses. 8th report of the International Committee on the Taxonomy of Viruses. Edited by: Fauquet CM, Mayo MA, Maniloff J. Amsterdam, The Netherlands: Elsevier Academic Press; 2005:757-778.

    Google Scholar 

  5. 5.

    Centers for Disease Control and Prevention: Clusters of acute respiratory illness associated with human enterovirus 68--Asia, Europe, and United States, 2008–2010. MMWR Morb Mortal Wkly Rep 2011, 60: 1301-1304.

    Google Scholar 

  6. 6.

    Tavares FN, Campos RM, Burlandy FM, Fontella R, de Melo MM, da Costa EV, da Silva EE: Molecular characterization and phylogenetic study of coxsackievirus A24v causing outbreaks of acute hemorrhagic conjunctivitis (AHC) in Brazil. PLoS One 2011, 6: e23206. 10.1371/journal.pone.0023206

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  7. 7.

    Tsai HP, Huang SW, Wu FL, Kuo PH, Wang SM, Liu CC, Su IJ, Wang JR: An echovirus 18-associated outbreak of aseptic meningitis in Taiwan: epidemiology and diagnostic and genetic aspects. J Med Microbiol 2011, 60: 1360-1365. 10.1099/jmm.0.027698-0

    PubMed  Article  Google Scholar 

  8. 8.

    Perevoscikovs J, Brila A, Firstova L, Komarova T, Lucenko I, Osmjana J, Savrasova L, Singarjova I, Storozenko J, Voloscuka N, Zamjatina N: Ongoing outbreak of aseptic meningitis in South-Eastern Latvia, June - August 2010. Euro Surveill 2010,15(32):9-11.

    Google Scholar 

  9. 9.

    Savolainen-Kopra C, Paananen A, Blomqvist S, Klemola P, Simonen ML, Lappalainen M, Vuorinen T, Kuusi M, Lemey P, Roivainen M: A large Finnish echovirus 30 outbreak was preceded by silent circulation of the same genotype. Virus Genes 2011, 42: 28-36. 10.1007/s11262-010-0536-x

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Wu PC, Huang LM, Kao CL, Fan TY, Cheng AL, Chang LY: An outbreak of coxsackievirus A16 infection: comparison with other enteroviruses in a preschool in Taipei. J Microbiol Immunol Infect 2010, 43: 271-277. 10.1016/S1684-1182(10)60043-6

    PubMed  Article  Google Scholar 

  11. 11.

    Minor P, Brown F, Domingo E: Virus Taxonomy. Classification and Nomenclature of Viruses. Sixth report of the International Committee on Taxonomy of Viruses. Edited by: Murphy FA, Fauquet CM, Bishop DHL. Vienna, Austria: Springer-Verlag; 1995:329-336.

    Google Scholar 

  12. 12.

    Caro V, Guillot S, Delpeyroux F, Crainic R: Molecular strategy for ‘serotyping’ of human enteroviruses. J Gen Virol 2001, 82: 79-91.

    PubMed  CAS  Article  Google Scholar 

  13. 13.

    Norder H, Bjerregaard L, Magnius LO: Homotypic echoviruses share aminoterminal VP1 sequence homology applicable for typing. J Med Virol 2001, 63: 35-44. 10.1002/1096-9071(200101)63:1<35::AID-JMV1005>3.0.CO;2-Q

    PubMed  CAS  Article  Google Scholar 

  14. 14.

    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  CAS  PubMed Central  Google Scholar 

  15. 15.

    Nasri D, Bouslama L, Pillet S, Bourlet T, Aouni M, Pozzetto B: Basic rationale, current methods and future directions for molecular typing of human enterovirus. Expert Rev Mol Diagn 2007, 7: 419-434. 10.1586/14737159.7.4.419

    PubMed  CAS  Article  Google Scholar 

  16. 16.

    Oberste MS, Maher K, Pallansch MA: Molecular phylogeny of all human enterovirus serotypes based on comparison of sequences at the 5’ end of the region encoding VP2. Virus Res 1998, 58: 35-43. 10.1016/S0168-1702(98)00101-4

    PubMed  CAS  Article  Google Scholar 

  17. 17.

    Ishiko H, Shimada Y, Yonaha M, Hashimoto O, Hayashi A, Sakae K, Takeda N: Molecular diagnosis of human enteroviruses by phylogeny-based classification by use of the VP4 sequence. J Infect Dis 2002, 185: 744-754. 10.1086/339298

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    King AMQ, Brown F, Christian P, Hovi T, Hyypia T, Knowles NJ, Lemon SM, Minor PD, Palmenberg AC, Skern T: Picornaviridae. In Virus Taxonomy. Seventh report of International Commitee for Taxonomy of Viruses. Edited by: van Regenmortel MHV, Fauquet CM, Bishop DHL. New York-San Diego, USA: Academic Press; 2000:657-673.

    Google Scholar 

  19. 19.

    Oberste MS, Michele SM, Maher K, Schnurr D, Cisterna D, Junttila N, Uddin M, Chomel JJ, Lau CS, Ridha W, al-Busaidy S, Norder H, Magnius LO, Pallansch MA: Molecular identification and characterization of two proposed new enterovirus serotypes, EV74 and EV75. J Gen Virol 2004, 85: 3205-3212. 10.1099/vir.0.80148-0

    PubMed  CAS  Article  Google Scholar 

  20. 20.

    Oberste MS, Maher K, Nix WA, Michele SM, Uddin M, Schnurr D, al-Busaidy S, Akoua-Koffi C, Pallansch MA: Molecular identification of 13 new enterovirus types, EV79-88, EV97, and EV100-101, members of the species human enterovirus B. Virus Res 2007, 128: 34-42. 10.1016/j.virusres.2007.04.001

    PubMed  CAS  Article  Google Scholar 

  21. 21.

    Lukashev AN: Role of recombination in evolution of enteroviruses. Rev Med Virol 2005, 15: 157-167. 10.1002/rmv.457

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Lukashev AN, Lashkevich VA, Ivanova OE, Koroleva GA, Hinkkanen AE, Ilonen J: Recombination in circulating human enterovirus B: independent evolution of structural and non-structural genome regions. J Gen Virol 2005, 86: 3281-3290. 10.1099/vir.0.81264-0

    PubMed  CAS  Article  Google Scholar 

  23. 23.

    Domingo E, Martin V, Perales C, Escarmis C: Coxsackieviruses and quasispecies theory: evolution of enteroviruses. Curr Top Microbiol Immunol 2008, 323: 3-32. 10.1007/978-3-540-75546-3_1

    PubMed  CAS  Google Scholar 

  24. 24.

    Racaniello VR: Picornaviridae: The Viruses and Their Replication. In Fields Virology. 5th edition. Edited by: Knipe DM, Howley PM. Philadelphia, PA: Lippincott Williams Wilkins; 2011:795-838.

    Google Scholar 

  25. 25.

    Alberta Health Services. Increased incidence of Enteroviruses causing Aseptic Meningitis this spring and summer.

  26. 26.

    Chang KH, Auvinen P, Hyypia T, Stanway G: The nucleotide sequence of coxsackievirus A9; implications for receptor binding and enterovirus classification. J Gen Virol 1989,70(Pt 12):3269-3280.

    PubMed  CAS  Article  Google Scholar 

  27. 27.

    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  CAS  PubMed Central  Article  Google Scholar 

  28. 28.

    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: 649-657. 10.1002/jmv.21652

    PubMed  CAS  Article  Google Scholar 

  29. 29.

    Kubo H, Iritani N, Seto Y: Molecular classification of enteroviruses not identified by neutralization tests. Emerg Infect Dis 2002, 8: 298-304. 10.3201/eid0803.010200

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  30. 30.

    Al-Hello H, Paananen A, Eskelinen M, Ylipaasto P, Hovi T, Salmela K, Lukashev AN, Bobegamage S, Roivainen M: An enterovirus strain isolated from diabetic child belongs to a genetic subcluster of echovirus 11, but is also neutralised with monotypic antisera to coxsackievirus A9. J Gen Virol 2008, 89: 1949-1959. 10.1099/vir.0.83474-0

    PubMed  CAS  Article  Google Scholar 

  31. 31.

    Abzug MJ: Presentation, diagnosis, and management of enterovirus infections in neonates. Paediatr Drugs 2004, 6: 1-10. 10.2165/00148581-200406010-00001

    PubMed  Article  Google Scholar 

  32. 32.

    Romero JR: Enteroviruses. In Clinical Virology. 3rd edition. Edited by: Richman DD, Whitley RJ, Hayden FG. Washington: ASM Press; 2009.

    Google Scholar 

  33. 33.

    Anderson RM, May RM: Infectious Diseases of Humans. Great Britain: Oxford Science Publications; 1991.

    Google Scholar 

  34. 34.

    Santti J, Harvala H, Kinnunen L, Hyypia T: Molecular epidemiology and evolution of coxsackievirus A9. J Gen Virol 2000, 81: 1361-1372.

    PubMed  CAS  Article  Google Scholar 

  35. 35.

    Fox JD, Han S, Samuelson A, Zhang Y, Neale ML, Westmoreland D: Development and evaluation of nucleic acid sequence based amplification (NASBA) for diagnosis of enterovirus infections using the NucliSens Basic Kit. J Clin Virol 2002, 24: 117-130. 10.1016/S1386-6532(01)00241-4

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Coiras MT, Aguilar JC, Garcia ML, Casas I, Perez-Brena P: Simultaneous detection of fourteen respiratory viruses in clinical specimens by two multiplex reverse transcription nested-PCR assays. J Med Virol 2004, 72: 484-495. 10.1002/jmv.20008

    PubMed  CAS  Article  Google Scholar 

  37. 37.

    Martino TA, Tellier R, Petric M, Irwin DM, Afshar A, Liu PP: The complete consensus sequence of coxsackievirus B6 and generation of infectious clones by long RT-PCR. Virus Res 1999, 64: 77-86. 10.1016/S0168-1702(99)00081-7

    PubMed  CAS  Article  Google Scholar 

  38. 38.

    Van de PY, De WR: Construction of evolutionary distance trees with TREECON for windows: accounting for variation in nucleotide substitution rate among sites. Comput Appl Biosci 1997, 13: 227-230.

    Google Scholar 

  39. 39.

    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  CAS  PubMed Central  Google Scholar 

Download references


This work was supported by the Provincial Laboratory for Public Health of Alberta and the University of Calgary.

Author information



Corresponding author

Correspondence to Raymond Tellier.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

EC carried out the molecular studies and participated in the sequence alignment. KP, SW and RT conceived of the study, participated in its design, coordination, and analysis. KP wrote the manuscript, SW and RT edited the manuscript. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Pabbaraju, K., Wong, S., Chan, E.N.Y. et al. Genetic characterization of a Coxsackie A9 virus associated with aseptic meningitis in Alberta, Canada in 2010. Virol J 10, 93 (2013).

Download citation


  • Human enteroviruses
  • Coxsackie A9
  • Aseptic meningitis
  • Serotyping