- Open Access
Human enterovirus 71 subgenotype B3 lacks coxsackievirus A16-like neurovirulence in mice infection
Virology Journal volume 2, Article number: 74 (2005)
At least three different EV-71 subgenotypes were identified from an outbreak in Malaysia in 1998. The subgenotypes C2 and B4 were associated with the severe and fatal infections, whereas the B3 virus was associated with mild to subclinical infections. The B3 virus genome sequences had ≥85% similarity at the 3' end to CV-A16. This offers opportunities to examine if there are characteristic similarities and differences in virulence between CV-A16, EV-71 B3 and EV-71 B4 and to determine if the presence of the CV-A16-liked genes in EV-71 B3 would also confer the virus with a CV-A16-liked neurovirulence in mice model infection.
Analysis of human enterovirus 71 (EV-71) subgenotype B3 genome sequences revealed that the 3D RNA polymerase and domain Z of the 3'-untranslating region RNA secondary structure had high similarity to CV-A16. Intracerebral inoculation of one-day old mice with the virus resulted in 16% of the mice showing swollen hind limbs and significantly lower weight gain in comparison to EV-71 B4-infected mice. None of the mice presented with hind leg paralysis typical in all the CV-A16 infected mice. CV-A16 genome sequences were amplified from the CV-A16-infected mice brain but no amplification was obtained from all the EV-71-inoculated mice suggesting that no replication had taken place in the suckling mice brain.
The findings presented here suggest that EV-71 B3 viruses had CV-A16-liked non-structural gene features at the 3'-end of the genome. Their presence could have affected virulence by affecting the mice general health but was insufficient to confer the EV-71 B3 virus a CV-A16-liked neurovirulence in mice model infection.
Enterovirus 71 (EV-71) was first described in 1969 during an outbreak with central nervous system complications in California . Since then, EV-71 infections have been associated with a number of outbreaks with wide clinical manifestations, ranging from mild hand, foot and mouth disease (HFMD) to severe neurological complications and deaths. These include outbreaks in Bulgaria , Hungary , Japan  and more recently Malaysia [5, 6], Taiwan  and Singapore . In the later three outbreaks, more than a hundred deaths in total were reported, elevating EV-71 infection as one of the most deadly virus infection to date amongst young children below the age of 3 years old in Asia. The sudden emergence of the deadly forms of EV-71 infection in Asia was puzzling, as the virus together with other human enterovirus A viruses especially coxsackievirus A5 (CV-A5), CV-A10 and CV-A16 have been noted to cause HFMD in the region for sometime . During the outbreak in Malaysia, at least three different EV-71 subgenotypes were identified. The subgenotypes C2 and B4 were associated with the severe and fatal infections, whereas, mild to subclinical infections were associated with the B3 viruses [10–12]. Unlike the earlier two subgenotypes, the B3 virus circulated for only a brief period during the outbreak and they have since not been isolated from patients from the later outbreaks [11, 12]. A recent study reported that the B3 virus genome sequences had ≥93% similarity to EV-71 at the 5' end whereas the P3 genome region and 3'UTR had ≥85% similarity to CV-A16 . CV-A16 is known to be the most common causative agent for the self-limiting HFMD. It is usually characterized by mild fever, oral ulcers and vesicular lesions on palms and soles and is not known to cause severe and fatal CNS infections. It is not presently understood why EV-71 infections tend to cause the more severe form of HFMD in comparison to CV-A16. The findings that EV-71 B3 viruses had high sequence similarity to CV-16 at the 3' end of the genome and that the viruses were not associated with the severe form of HFMD, offered opportunities to examine the potential roles of the respective genes in determining virulence. Hence, the present study was undertaken to examine if there are characteristic similarities and differences between CV-A16, EV-71 B3 and the more virulent EV-71 B4 virus and to determine if the presence of the CV-A16-liked genes in the EV-71 B3 virus genome would also confer the virus a CV-A16-liked neurovirulence in mice.
Results and Discussion
The consensus amino acid sequences of the two available EV-71 B3 virus genomes (SHA63 and SHA66) were compared to other available subgenotype B4 and CV-A16/G10 genome sequences from the Genbank. Several amino acids (His1775, Thr1947, Ile1806, Gln1825, Thr1928, Thr1947, Asn2099, Glu2114 and Gln2159) that were characteristic of the CV-A16/G10 were found in EV-71 B3 isolates. These amino acid differences occurred only within the 3D RNA polymerase gene, suggesting that this gene is very much CV-A16 than it is EV-71. Comparisons of the EV-71 B3 amino acid sequences against all other EV-71 and CV-A16 also revealed at least 12 amino acids (Asn1124, Arg1152, Ser1335, Ser1641, Tyr1799, Asp1822, Val1860, Ser1864, Val1997, Ala2039 Asp2101 and Leu2125) that were unique to the EV-71 B3 isolates. Eight of these amino acid differences occurred within the 3D RNA polymerase gene. Two of these unique mutations found were located between amino acids 176–348 genome region essential for RNA-protein interactions  (Fig. 1). Alignment of the EV-71 B3 (SHA66) and EV-71 B4 (UH1) isolates RNA polymerase against the three-dimensional crystal structure of poliovirus 1 Mahoney strain 3D RNA polymerase (PDB: 1RDR) was performed to locate these mutations. Of these eight mutations in EV-71 B3 virus, three were located within the finger subdomain and two were located at the palm motif suggesting that the EV-71 B3 virus amino acid substitutions were mainly located within the 3D RNA polymerase functional domains. The highly 'flexible' finger domain is involves in modulating substrate recognition and oligomerization of the polymerase for binding to nucleotides . In poliovirus, mutations within the 3D RNA polymerase located to the 3' end of the genome have been shown to affect neurovirulence [16, 17]. Hence, this highlights the potential importance of the 3D RNA polymerase in determining the virus neurovirulence. It was also found that in addition to the presence of CV-A16 or CV-A16-liked 3D RNA polymerase gene sequences, the EV-71 B3 viruses also shared a similar predicted 3' UTR secondary structures with CV-A16/G10 at domain Z (Fig. 2), a domain reported as important in determining cardiovirulence of CV-B3 . Mutations that affect the stem-and-loop structures have been shown earlier to abolish infectivity and virus RNA synthesis [19, 20]. The predicted domain Y known to form a tertiary RNA 'kissing' structure with domain X of the EV-71 B3 virus, however, differed from the EV-71 B4 and CV-A16/G10 (Fig. 2).
Inoculation of one day-old newborn mice showed that all mice inoculated with CV-A16 had the typical signs and symptoms of CV-A16 infections by day two post-inoculation. The mice were lethargic, had floppy tails, tremoring, uncoordinated movement and reduced average body weight in comparison to EV-71 B3- or EV-71 B4-inoculated mice (Fig. 3a,3g). Approximately 17% (4/24) of the mice had hind leg paralysis by day three post-inoculation and one died (Fig. 3b,3e,3f, Additional file: 1). By day four post-inoculation, all the CVA16-inoculated mice had developed hind leg paralysis and subsequently died (Fig. 3b,3e,3f). A 150 bp enterovirus genome sequence were amplified and sequenced from the total RNA of the brain of all the CV-A16-inoculated mice confirming the presence of CV-A16 in the mice brain (Fig. 4). Mice-inoculated with EV-71 B3 and EV-71 B4 viruses also had significantly reduced average body weight in comparison to the control mock-infected mice (Student's t-test, P < 0.05, Fig. 3d,3g). Mice inoculated with EV-71 B3 virus, however, had significantly reduced average body weight in comparison to those inoculated with the EV-71 B4 virus (Fig. 3g). These mice appeared lethargic and uncoordinated beginning on day two post-inoculation. Of these, 16% (4/25) developed swollen hind legs and one subsequently died on day five post-inoculation (Fig. 3c,3e,3f). There were no hind leg paralysis noted and the remaining surviving mice recovered, fed well and regained balance after day six post-inoculation. In contrast, about 20% (6/31) of the mice inoculated with EV-71 B4 virus developed swollen fore limbs or hind legs and of these, three died after day four post-inoculation (Fig. 3e,3f). After day eight post-inoculation, the B4-inoculated mice also recovered, became more active and fed well. Pairwise comparison of the clinical illness and survival probability between the virus-inoculated groups and control were significant suggesting that the three viruses, CV-A16, EV-71 B3 and EV-71 B4 viruses caused death in mice (log rank survival analysis, P < 0.05, Fig. 3e,3f) but only infection with CV-A16 lead to 100% mortality. In contrast to CV-A16 infection, no amplification of the enterovirus sequence was detected in the selected EV-71 B4- and EV-71 B3-inoculated mice brain, suggesting that EV-71 B3 and EV-71 B4 viruses perhaps did not replicate in the mice brain when introduced intracerebrally (Fig. 4). This may help to explain the absence of hind leg paralysis in all the EV-71-infected mice and the complete recovery of all the surviving mice. Death seen amongst these mice may have been caused by infection of other tissues as manifested in mice with swollen limbs and legs. Evidence suggesting that EV-71 strains isolated during the Bulgaria poliomyelitis-like epidemic had higher tropism for mouse muscle tissues than the brain tissues  and EV-71 neurovirulence mimicking human infection was achieved only from using a mouse-adapted virus strain but not the parental strain [21, 22] support the findings from the present study that EV-71 B3 and B4 did not infect the brain. The infection, however, manifests clinically in some mice as non-specific swollen limbs and legs. Hence it is possible that, though both EV-71 and CV-A16 viruses are closely related, different receptors are utilize for the respective virus entry into the different tissues and this could be mediated through the virus structural proteins. The mutations that occurred within the 3D RNA polymerase of the EV-71 B3 virus along with the presence of CV-A16-liked 3' UTR domain Z RNA secondary structure then could contribute to virulence but by themselves did not affect EV71 neurovirulence in mice as in contrast to CV-A16, the B3 virus lacks tropism for the mice brain. Since the major differences between the EV-71 B4 and EV-71 B3 viruses occurred at the 3' end of the genome, this support the view that the structural genes of EV-71 and CV-A16 determined tissue tropisms.
Results from the present study, also did not support the possibility that acquisition of CV-A16-liked genome sequences alone is sufficient to confer the EV-71 B3 virus a CVA16-liked neurovirulence in mice. The significant mice weight gain differences noted between mice infected with EV-71 B3 and EV-71 B4 viruses, with the later performing much better, however, suggested that EV-71 B3 virus infection somehow did affect mice general health. As weight gain differences are the only biological parameter that differentiate between the B3 and B4 viruses, it does appears that EV-71 B3 affected mice more than the EV-71 B4 virus. It is also worth noting that in contrast to infection in mice, CV-A16 infection in human in general does not result in severe infection as oppose to EV-71, particularly the EV-71 B4 virus infection. In parallel manner, the EV-71 B3 viruses, while they affected mice, they did not cause severe or fatal infection in humans. These implied that the EV-71 B3 virus is truly different and as its genome suggested, it has to some extent features of both EV-71 and CV-A16 infection in-vivo.
Results from the present study suggest that EV-71 B3 virus had CV-A16-liked non-structural gene 3D RNA polymerase and 3' UTR features at the 3' end of the genome. Their presence affected virulence differently from infection with EV-71 B4 and CV-A16 by affecting the mice general health. The presence of the CV-A16-liked genes, however, was insufficient to markedly influence the neurovirulence properties of EV-71 B3 virus in mice.
Materials and methods
Two EV-71 isolates identified from the 1997 HFMD outbreak in Malaysia were used. The subgenotype B3 isolate, SHA66 (EMBL: AJ238457) was isolated from a HFMD patient presented with mild infection [6, 23]. The subgenotype B4 isolate, UH1 (EMBL: AJ238455) on the other hand, was isolated from the brain of a patient who died of EV-71-associated neurogenic pulmonary edema [5, 6, 24]. The CV-A16 isolate used was previously isolated from a HFMD patient seen at the University Malaya Medical Centre. This CV-A16 isolate was identified and characterized using monoclonal antibody staining (Chemicon Cat #3323, California, USA) and amplification of partial 5' UTR gene (data not shown).
Amino acid sequence analysis
Amino acid sequences were examined after stripping the 5' UTR and 3' UTR sequences and consensus sequences of EV-71 B3 and EV-71 B4 viruses were aligned and manually edited using GeneDoc software . The previously published three-dimensional crystal structure of the 3D RNA polymerase was downloaded as template for the alignment. Using the WHAT IF program , domains that represent the conserved regions, loops, insertion or deletions were manually visualized to generate a structural alignment.
RNA secondary structure prediction
The 3' UTR RNA secondary structure was predicted using Zuker optimal and suboptimal minimal free energy folding algorithms, as implemented in RNA Structure version 3.71 software . Part of the poly A tract was incorporated into the sequences.
Determination of virulence in mice
A total of 24, 25 and 31 one-day old newborn ICR mice were inoculated intracerebrally with either CV-A16 or SHA66 (B3 virus) or UH1 (B4 virus) virus inoculum. The virus inoculum with infectivity of ~1 × 103 p.f.u. was injected in a volume of 10–20 μl into the mice brain. The mice were closely monitored for any clinical symptoms, paralysis and death and the weight of each surviving mouse was recorded daily up to day 11 post-inoculation. Another litter with at least 10 one-day old newborn mice was injected with comparable growth medium and used as controls. At selected intervals post-infection, some of the mice were sacrificed and the brain tissues were harvested for total RNA using the TRI Reagent™ (Molecular Research Centre, Inc., Cincinnati, USA) following the manufacturer's recommended protocols. The RT-PCR amplification for the detection of enterovirus sequence was performed using 1 μg of RNA. Access RT-PCR kit (Promega, USA) and primer pairs, EntabF (5'-TCC TCC GGC CCC TGA ATG CGG CTA AT-3'; nucleotide positions 449–474, based on MS87 strain, Genbank: U22522) and EVRR (5'-AAT TGT CAC CAT AAG CAG GC-3'; nucleotide positions 586–606) were used. Reverse transcription was performed at 42°C for an hour followed by amplification steps; 95°C-30 seconds, 55°C-30 seconds and 72°C-30 seconds for 30 cycles and finally with 5 minutes extension at 72°C using the PTC thermal cycler (MJ Research, Massachusetts, USA). When no amplicon was obtained, the number of cycle was increased to 40. Alternatively, a second step PCR using similar parameters was performed using ten-fold diluted RT-PCR product as template. The amplified DNA fragments were electrophoresed using 2% agarose gel in 0.5 × tris-acetate EDTA buffer (0.02 M Tris base, 0.5 mM EDTA pH 8.0, 0.057% glacial acetic acid) and sequence confirmation was made by sequencing the DNA fragment.
Student's t-test was used to evaluate if the differences in weight between the virus-inoculated mice and control mice was significant. Wilcoxon signed rank test was used to compare the survival and paralysis probability between the virus-inoculated mice and control mice. All statistical analyses were implemented using SPSS for Windows version 11.5 (SPSS Inc., Illinois, USA). All tests were two-sided and P < 0.05 was considered as statistically significant.
- HFMD Hand:
foot and mouth disease
Schmidt NJ, Lennette EH, Ho HH: An apparently new enterovirus isolated from patients with disease of the central nervous system. J Infect Dis 1974, 129: 304-309.
Chumakov M, Voroshilova M, Shindarov L, Lavrova I, Gracheva L, Koroleva G, Vasilenko S, Brodvarova I, Nikolova M, Gyurova S, Gacheva M, Mitov G, Ninov N, Tsylka E, Robinson I, Frolova M, Bashkirtsev V, Martiyanova L, Rodin V: Enterovirus 71 isolated from cases of epidemic poliomyelitis-like disease in Bulgaria. Arch Virol 1979, 60: 329-340. 10.1007/BF01317504
Nagy G, Takatsy S, Kukan E, Mihaly I, Domok I: Virological diagnosis of enterovirus type 71 infections: experiences gained during an epidemic of acute CNS diseases in Hungary in 1978. Arch Virol 1982, 71: 217-227. 10.1007/BF01314873
Hagiwara A, Yoneyama T, Takami S, Hashimoto I: Genetic and phenotypic characteristics of enterovirus 71 isolates from patients with encephalitis and with hand, foot and mouth disease. Arch Virol 1984, 79: 273-283. 10.1007/BF01310816
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.
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
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.
Singh S, Chow VTK, Chan KP, Ling AE, Poh CL: RT-PCR, nucleotide, amino acid and phylogenetic analyses of enterovirus type 71 strains from Asia. J Virol Meth 2000, 88: 193-204. 10.1016/S0166-0934(00)00185-3
Infectious agent surveillance report (IASR). Isolation and detection report of viruses by years 1982–July 2005: enterovirus 1 & enterovirus 2[http://idsc.nih.go.jp/iasr/virus/virus-kk.html]
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 2001, 75: 7732-7738. 10.1128/JVI.75.16.7732-7738.2001
Cardosa MJ, Perera D, Brown BA, Cheon D, Chan HM, Chan KP, Cho H, 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.
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
Chan YF, AbuBakar S: Recombinant human enterovirus 71 in hand, foot and mouth disease patients. Emerg Infect Dis 2004, 10: 1468-1470.
Cornell CT, Semler BL: Subdomain specific functions of the RNA polymerase region of poliovirus 3CD polypeptide. Virology 2002, 298: 200-213. 10.1006/viro.2002.1481
Hansen JL, Long AM, Schultz SC: Structure of the RNA-dependent RNA polymerase of poliovirus. Structure 1997, 5: 1109-1122. 10.1016/S0969-2126(97)00261-X
Christodoulou C, Colbere-Garapin F, Macadam A, Taffs LF, Marsden S, Minor P, Horaud F: Mapping of mutations associated with neurovirulence in monkeys infected with Sabin 1 poliovirus revertants selected at high temperature. J Virol 1990, 64: 4922-4929.
Tardy-Panit M, Blondel B, Martin A, Tekaia F, Horaud F, Delpeyroux F: A mutation in the RNA polymerase of poliovirus type 1 contributes to attenuation in mice. J Virol 1993, 67: 4630-4638.
Merkle I, van Ooij MJ, van Kuppeveld FJ, Glaudemans DH, Galama JM, Henke A, Zell R, Melchers WJ: Biological significance of a human enterovirus B-specific RNA elements in the 3' nontranslated region. J Virol 2002, 76: 9900-9909. 10.1128/JVI.76.19.9900-9909.2002
Melchers WJ, Hoenderop JG, Bruins Slot HJ, Pleij CW, Pilipenko EV, Agol VI, Galama JM: Kissing of the two predominant hairpin loops in the coxsackie B virus 3' untranslated region is the essential structural feature of the origin of the replication required for the negative-strand RNA synthesis. J Virol 1997, 71: 686-696.
Mirmomeni MH, Hughes PJ, Stanway G: An RNA tertiary structure in the 3' untranslated region of enteroviruses is necessary for efficient replication. J Virol 1997, 71: 2363-2370.
Chen YC, Yu CK, Wang YF, Liu CC, Su IJ, Lei HY: A murine oral enterovirus 71 infection model with central nervous system involvement. J Gen Virol 2004, 85: 69-77. 10.1099/vir.0.19423-0
Wang YF, Chou CT, Lei HY, Liu CC, Wang SM, Yan JJ, Su IJ, Wang JR, Yeh TM, Chen SH, Yu CK: A mouse-adapted enterovirus 71 strain causes neurological disease in mice after oral infection. J Virol 2004, 78: 7916-7924. 10.1128/JVI.78.15.7916-7924.2004
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
AbuBakar S, Shafee N, Chee HY: Induction of apoptosis in neurogenic pulmonary edema. Med J Malaysia 1999, 54: 402-403.
Nicholas KB, Nicholas HB Jr, Deerfield DW: GeneDoc: analysis and visualization of genetic variation. EMBNEW NEWS 1997, 4: 14.
Vriend G: WHAT IF: A molecular modeling and drug design program. J Mol Graph 1990, 8: 52-56. 10.1016/0263-7855(90)80070-V
Matthews DH, Sabina J, Zuker M, Turner DM: Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J Mol Biol 1999, 288: 911-940. 10.1006/jmbi.1999.2700
This study is funded in parts by grants from the Ministry of Science, Technology and Innovation, Malaysia # 06-02-09-001-BTK/TD/002.
The author(s) declare that they have no competing interests.
The corresponding author, Sazaly AbuBakar is the principal investigator of the study; is involved in the design, supervision, data analyses and writing of the report. Chan Y-F performed all the virological investigations, nucleotide sequencing and analyses of data. All authors were involved in the preparation of this "Research Article" and figures.