Detection and whole genome sequence analysis of an enterovirus 68 cluster
Virology Journal volume 10, Article number: 103 (2013)
Enteroviruses are a common cause of human disease and are associated with a wide range of clinical manifestations. Enterovirus 68 is rarely detected yet was reported in many countries in 2010. Here enterovirus 68 was identified for the first time in New Zealand in 2010 and was detected in a further fourteen specimens over a six month period.
To genetically characterise enterovirus 68 specimens identified in New Zealand in 2010.
The genome sequence of a New Zealand representative enterovirus 68 isolate was obtained. Ten clinical specimens were analysed by sequencing the VP1 region of the enterovirus 68 genome.
Based on sequence analysis of the VP1 region and the full genome of one representative isolate, the New Zealand enterovirus 68 isolates clustered with contemporary enterovirus 68 viruses and do not show any clear distinguishing genetic diversity when compared to other strains. All fifteen specimens showed high similarity with enterovirus 68 by VP1 sequencing. The majority of New Zealand patients suffered from bronchiolitis, were less than two years of age and were of Pacific Island or Maori descent.
We document the rare occurrence of an enterovirus 68 cluster in New Zealand in 2010. These viruses shared similarity with other clusters of enterovirus 68 that occurred globally in 2010. A greater awareness in enterovirus 68 infection may help detect this virus with increased frequency and enable us to better understand the role this strain plays in disease and the reasons behind this global emergence in 2010.
Enteroviruses belong to the Picornaviridae family and are among the most commonly identified aetiological agents of human disease . There are approximately 100 enterovirus serotypes which cause a range of clinical manifestations; from asymptomatic infections to more serious illnesses such as aseptic meningitis, myocarditis and acute flaccid paralysis [1, 2]. Enterovirus 68 (EV68) is a member of the Human enterovirus D species and was first isolated in California, USA in 1962 from children who were hospitalised with lower respiratory tract infections . Since then, EV68 has been isolated rarely; only 26 strains have been identified over 36 years in the USA . EV68 is unique among enteroviruses in that it has a lower than optimum growth temperature and is acid sensitive [4–6]. As such, it shares characteristics with human rhinovirus [4, 5]. It is further unique in that it is almost exclusively associated with respiratory disease [4, 5]. Recently EV68 has been isolated with increased frequency. Its isolation has been reported in Germany, the Philippines, Thailand, Italy, Japan, the United States, the United Kingdom and the Netherlands with the majority of these reports occurring in 2010 [7–15]. To date, only four full genome sequences of EV68 have been published; that of the prototype Fermon strain, the French 37–99 strain and two EV68 strains that were circulating in Japan in 2010 [5, 9, 16].
In the current study we describe fifteen cases of EV68 isolated from samples taken from March to August 2010 in New Zealand. All cases were initially identified by partial VP1 sequencing. Due to limiting sample volumes complete VP1 sequencing was performed on only ten of the fifteen EV68 samples confirming EV68 infection. Additionally, characterisation of the full genome sequence of a representative New Zealand EV68 isolate was achieved by Roche 454 sequencing.
Patients and specimen collection
The National Poliovirus and Enterovirus Identification Reference Laboratory at the Institute of Environmental Science and Research Limited, National Centre for Biosecurity and Infectious Disease routinely receives untyped enterovirus clinical specimens or cell culture isolates from four major hospitals (based in Auckland, Waikato, Wellington and Christchurch) as part of the New Zealand enterovirus surveillance network.
Viruses and cells
Human rhabdomyosarcoma (RD) cells (passage 242–256) were propagated in 10% Hanks Minimal Essential Medium (Gibco, Life Technologies, Carlsbad, CA, USA) supplemented with 10% (v/v) foetal bovine serum (HyClone, New Zealand), 7.5% (v/v) sodium bicarbonate, 1% (v/v) 1 M hepes and antibiotics.
Viral nucleic acid was extracted (400 μL of clinical specimen or cell culture isolate) using the Zymo ZR Viral RNA Kit™ (Zymo Research Corporation, Irvine, CA, USA) as per the manufacturer’s instructions. Nucleic acid extracts (20 μL) were stored at −80°C until required.
Partial and full VP1 RT-PCR and sequencing
Initially, specimens were characterised by amplification and sequencing of a 375 bp partial VP1 region as described previously . An EV68 RT-PCR assay was designed in-house using Primer3 v 0.4.0  in order to specifically amplify and sequence the entire VP1 region. The forward primer was designated EV68-VP1-Forward (5’-GCA-GCC-TAT-CAG-GTG-GAG-AG-3’) and the reverse primer was designated EV68-VP1-Reverse (5’-TGC- TCA-TGT-ATG-GCA-TGG-TT-3’) resulting in a product that was 805 bp in length. One-step RT-PCR assays were performed using the SuperScript III Platinum One-Step System (Invitrogen, Life Technologies). Each 50 μL reaction contained 5 μL of RNA, 25 μL of 2X Reaction Mix, 0.5 μL of reverse transcriptase-Taq polymerase enzyme mix, 20 units of RNase inhibitor(Invitrogen) and 0.2 μM of primers. Following an initial 50 minute reverse transcription at 50°C and 15 minute denaturation step at 94°C, a three-step cycling procedure of denaturation at 94°C for 30 seconds, annealing at 55°C for 30 seconds and extension at 72°C for 1 minute over 40 cycles was used. A final extension was performed at 72°C for 5 minutes. Each assay included negative and positive extraction controls and RNase-free reagents and procedures were used to minimise contamination.
Samples were run on an E-Gel® EX 2% precast agarose gel (Invitrogen). Bands of interest were purified using the Qiagen QIAquick PCR Purification Kit (Qiagen Incorporated, Valencia, CA, USA) and used as templates for sequencing on an ABI 3130XL automated DNA sequence using ABI Big Dye v3.1 technology (Applied Biosystems, Life Technologies).
Sequences were edited and assembled using Bionumerics 5.10 (Applied Maths, Belgium) and compared to known VP1 sequences in GenBank using BLAST . Samples were determined to be EV68 strains if they were ≥ 75% similar to the prototype EV68 Fermon strain (GenBank accession number: AY426531).
The phylogenetic relationships were inferred using the Neighbour-Joining method (Tamura 3-parameter model with gamma distributed rates among sites) on nucleotide sequence data from the VP1 region, incorporating 651 base pairs. All positions containing gaps and missing data were eliminated. Evolutionary analyses were conducted in MEGA5  (Center for Evolutionary Medicine and Informatics, Tempe, AZ, USA). 1000 bootstrap replicates were performed to determine the consensus tree.
Enterovirus genome sequencing
Full-length genomic sequence of EV68 was obtained by sequencing on a Roche GS FLX instrument. Passage 1 viral culture was firstly purified to enrich for virus particles. A 1 mL volume of viral culture was subjected to three cycles of freeze-thaw lysis between −80°C and 37°C and then to centrifugation at 3000 × g for 15 minutes. To remove cellular RNA, 400 μL of supernatant was treated with RNAse A (Ambion, Life Technologies) at a final concentration of 10 μg.mL -1 at 37°C for one hour. RNA was extracted from this fraction using the Ambion MagMax Viral RNA kit (Ambion) with elution into 40 μL of EB buffer. The RNA concentration was not high enough for direct application to the Roche GS FLX sequencer so a random amplification method for low yield samples was used, as has been described previously . 500 ng of this material was used to make a barcoded, GS FLX Rapid library (Roche, Mannheim, Germany) according to the manufacturer’s instructions except that the sample was not nebulized. Completed library length was between 500 bp and 900 bp and the concentration was 8.35 × 109 molecules/μL. Two emulsion PCR reactions were performed using a library concentration of 4 counts per bead and the enriched beads combined to give a total of 137,200. These were loaded onto a 1/16th region of a pico titre plate and sequenced using the XLR70 Titanium chemistry (Roche). The utilities sffinfo  and FASTX-Toolkit  were used to remove adaptor and primer sequences and reads filtered based upon data quality. This also increased BLAST search performance because redundant sequence information was not submitted. Human EV68 sequence was assembled using a reference strain, Fermon (GenBank accession number: AY426531). Assembly and annotation was achieved using Geneious 5.5.6 (Biomatters Ltd., New Zealand), on enterovirus hits obtained from the BLASTN search. Areas of assembly with less than 5X coverage were resequenced by the sanger-sequencing methodology using custom designed PCR primers.
Patients and clinical data
In 2010, EV68 was detected in a total of fifteen samples. EV68 positive patients ranged in age from 1 month to 48 years, median age 11 years and 0 months. The majority of patients (n = 11; 73.3%) were less than 2 years of age. Patients testing positive for EV68 were predominantly male (n = 10; 66.7%) and most were identified as being of Pacific Island or Maori descent (n = 12; 80.0%) (Table 1). The first EV68 positive sample was taken in March of 2010. The peak incidence was observed in June of 2010 (n = 7; 46.7%) and this coincides with the winter period in New Zealand, typical of respiratory disease. Most of the samples originated in the Auckland region (n = 13; 86.7%) with the remaining two samples coming from the Waikato region. Respiratory samples included throat swabs (n = 1; 6.7%), nasopharyngeal swabs (NPS) (n = 4; 26.7%), nasopharyngeal aspirates (NPA) (n = 9; 60.0%) and a tracheal aspirate (n = 1; 6.7%). Patients presented with bronchiolitis (n = 6; 40.0%), asthma (n = 2; 13.3%) or other symptoms including cough, stridor, coryza and wheezing (Table 1).
VP1 Phylogenetic analysis
Based upon sequencing of the VP1 gene segment of ten New Zealand clinical isolates, the viruses grouped in one distinct phylogenetic cluster (Figure 1). Enteroviruses were judged to be the same serotype if they were ≥ 75% similar to the prototype Fermon strain (GenBank accession number: AY426531) but had < 70% sequence similarity to any other reference strain . The ten New Zealand EV68 isolates had an average of 88.6% similarity to the VP1 of the prototype Fermon strain and only 66.5-68.8% similarity to that of enterovirus 70 (GenBank accession number: DQ201177) and 94 (GenBank accession number: DQ916376) (other species D enteroviruses), thus confirming these viruses as EV68. The New Zealand isolates were 98.0-99.5% identical to each other. The New Zealand EV68 VP1 nucleotide sequences were compared with the sequences of EV68 from the Philippines, Japan and the Netherlands. The VP1 phylogenetic tree demonstrated that the New Zealand EV68 isolates clustered most closely with EV68 isolates that were detected in the Netherlands in 2010 but were not sufficiently distinct from any other EV68 sequence so as to identify them as a genetically distinct group.
Enterovirus whole genome sequencing
We describe the fifth published complete genome of an EV68 virus with the exception of 32 nucleotides in the 5’ UTR which could not be determined either by high throughput or sanger sequencing. The Fermon (GenBank accession number: AY426531), 37–99 (GenBank accession number: EF107098) and two Japanese strains (GenBank accession numbers: AB601882 and AB601883) are the only EV68 strains for which complete genome sequences are available. High-throughput sequencing of NZ-2010-541 (passage 1, viral culture) on 1/16th of a Roche GS FLX Titanium plate produced 41,349 pass filter reads. 40,091 reads were identified as human enterovirus by a BLASTN search of the GenBank non-redundant nucleotide database. 26,523 high quality reads were used to assemble a genome sequence for NZ-2010-541 of 7,320 bp (GenBank accession number: JX070222). The closest complete genome sequence identified was that of the strain 37–99 (GenBank accession number: EF107098) which showed 93.40% nucleotide sequence identity overall (Table 2). The genome sequence of NZ-2010-541 is unremarkable when compared to the other published genomes and no evidence for recombination was revealed.
We have reported an emergence of EV68 infections in New Zealand in 2010. EV68 was detected in a total of fifteen specimens between March and August 2010. We also describe a full length genome sequence of a representative EV68 isolate, which contributes to the four published full genome sequences of EV68; that of the prototype Fermon, that of the French 37-99 strain and two EV68 strains that were circulating in Japan in 2010 [5, 9, 16]. Until recently, EV68 was only isolated rarely . However in the past couple of years EV68 has been isolated with increased frequency with its isolation being reported in various countries worldwide beginning in 2009 [7–13]. The isolation of enterovirus 68 in New Zealand in 2010 is consistent with this trend. EV68 is thought to be associated with more severe lower respiratory tract infections such as bronchiolitis and pneumonia as well as causing exacerbations of asthma . New Zealand EV68 patients primarily presented with bronchiolitis, cough, coryza, wheezing and the exacerbation of asthma. Eight patients were admitted to hospital although two of those patients had underlying illness. The high number of EV68 hospitalisations within a short period of time suggests that EV68 could be highly transmissible in the community and capable of causing more severe disease particularly in young children.
The majority (80.0%) of EV68 strains identified in 2010 in New Zealand were from patients of Pacific Island or Maori descent. Of note was that the other 20.0% of cases specified their ethnicity as other suggesting ethnicity other than New Zealand European, Maori, Pacific Island or Asian. Despite the over-representation of theses ethnicities being disproportionately affected by EV68 this dataset is not large enough to conclude EV68 infection is more likely to affect New Zealand Pacific Island and Maori populations but does warrant enhanced surveillance or further study.
An alternative explanation for this is the bias of hospital users, as Pacific Island or Maori ethnic groups are usually over-represented in hospital admissions for acute respiratory tract infections [24–27]. A number of factors such as nutrition and living conditions, exposure to cigarette smoke and lack of early medical care due to cost may lead to the higher hospitalisation rates among these communities [24, 26, 27]. New Zealand minority groups are not alone in experiencing higher hospitalisation rates. Internationally, minority groups show higher hospital admission rates for respiratory infections than those in the majority culture [28–30]. Phylogenetic analysis of the variable VP1 region of the New Zealand EV68 isolates showed that our clinical isolates grouped in one distinct phylogenetic cluster. Comparing the New Zealand EV68 VP1 nucleotide sequences with those from the Philippines, Japan and the Netherlands showed that the New Zealand isolates clustered with EV68 isolates that were detected in the Netherlands in 2010. EV68 isolates from the Netherlands clustered into two distinct phylogenetic groups; those detected in 2009 that were circulating in the community and those detected in 2010 from EV68 positive patients that required hospitalisation . Additional analysis of the EV68 isolates from the Netherlands demonstrated an amino acid deletion at residue 242 in the VP1 protein in the 2010 cluster compared with the 2009 cluster. This deletion was also observed in the New Zealand isolate NZ-2010-541 that was sequenced in this study . This deletion was common to both New Zealand and Netherland hospitalised patients but not to Netherland patients treated in the community and this could indicate increased virulence leading to more severe infection . Further studies are required to understand the role of this deletion in relation to virulence and disease severity and confirm this hypothesis.
It was expected that a full coverage of the NZ-2010-541 genome was not achieved by high throughput sequencing, despite our anticipation that 40,000 sequence reads with an average read length of approximately 500 base pairs should provide an average coverage of approximately 2700X. We noted 17 gaps in the high throughput sequence data with less than 5X coverage which ranged in size from 3 to 180 bp, and custom sanger sequencing primers were designed to fill these gaps.
In summary, this report describes the occurrence of a rare EV68 infection associated with lower respiratory tract infection in young children in New Zealand. The national poliovirus/enterovirus surveillance system is an essential and effective tool in identifying emerging and re-emerging enterovirus outbreaks and has the potential to inform upon public health interventions. An increased awareness of EV68 may help to detect this virus with greater frequency thus providing valuable information to healthcare providers. Generating useful phylogenetic and epidemiological data may help to reveal the reason for the increased activity of this virus in 2010.
Pallansch M, Roos R: Enteroviruses, polioviruses, echoviruses, and newer enteroviruses. In Fields virology. Edited by: Fields BN, Knipe DM, Howley PM, Griffin DE, Lamb RA, Martin MA. Philadelphia: Wolters Kluwer Health/Lippincott, Williams & Wilkins; 2007.
Khetsuriani N, Lamonte-Fowlkes A, Oberst S, Pallansch MA: Enterovirus surveillance–United States, 1970–2005. MMWR Surveill Summ 2006, 55: 1-20.
Schieble JH, Fox VL, Lennette EH: A probable new human picornavirus associated with respiratory diseases. Am J Epidemiol 1967, 85: 297-310.
Blomqvist S, Savolainen C, Raman L, Roivainen M, Hovi T: Human rhinovirus 87 and enterovirus 68 represent a unique serotype with rhinovirus and enterovirus features. J Clin Microbiol 2002, 40: 4218-4223. 10.1128/JCM.40.11.4218-4223.2002
Oberste MS, Maher K, Schnurr D, Flemister MR, Lovchik JC, Peters H: Enterovirus 68 is associated with respiratory illness and shares biological features with both the enteroviruses and the rhinoviruses. J Gen Virol 2004, 85: 2577-2584. 10.1099/vir.0.79925-0
Savolainen C, Blomqvist S, Mulders MN, Hovi T: Genetic clustering of all 102 human rhinovirus prototype strains: serotype 87 is close to human enterovirus 70. J Gen Virol 2002, 83: 333-340.
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.
Imamura T, Fuji N, Suzuki A, Tamaki R, Saito M, Aniceto R: Enterovirus 68 among children with severe acute respiratory infection, the Philippines. Emerg Infect Dis 2011, 17: 1430-1435.
Kaida A, Kubo H, Sekiguchi J, Kohdera U, Togawa M, Shiomi M: Enterovirus 68 in children with acute respiratory tract infections, Osaka. Japan Emerg Infect Dis 2011, 17: 1494-1497.
Meijer A, van der Sanden S, Snijders BE, Jaramillo-Gutierrez G, Bont L, van der Ent CK: Emergence and epidemic occurrence of enterovirus 68 respiratory infections in The Netherlands in 2010. Virology 2012, 423: 49-57. 10.1016/j.virol.2011.11.021
Piralla A, Baldanti F, Gerna G: Phylogenetic patterns of human respiratory picornavirus species, including the newly identified group C rhinoviruses, during a 1-year surveillance of a hospitalized patient population in Italy. J Clin Microbiol 2011, 49: 373-376. 10.1128/JCM.01814-10
Rahamat-Langendoen J, Riezebos-Brilman A, Borger R, van der Heide R, Brandenburg A, Scholvinck E: Upsurge of human enterovirus 68 infections in patients with severe respiratory tract infections. J Clin Virol 2011, 52: 103-106. 10.1016/j.jcv.2011.06.019
Renwick N, Schweiger B, Kapoor V, Liu Z, Villari J, Bullmann R: A recently identified rhinovirus genotype is associated with severe respiratory-tract infection in children in Germany. J Infect Dis 2007, 196: 1754-1760. 10.1086/524312
Linsuwanon P, Puenpa J, Suwannakarn K, Auksornkitti V, Vichiwattana P, Korkong S: Molecular epidemiology and evolution of human enterovirus serotype 68 in Thailand, 2006–2011. PLoS One 2012, 7: e35190. 10.1371/journal.pone.0035190
Lauinger IL, Bible JM, Halligan EP, Aarons EJ, Macmahon E, Tong CY: Lineages, sub-lineages and variants of enterovirus 68 in recent outbreaks. PLoS One 2012, 7: e36005. 10.1371/journal.pone.0036005
Smura TP, Junttila N, Blomqvist S, Norder H, Kaijalainen S, Paananen A: Enterovirus 94, a proposed new serotype in human enterovirus species D. J Gen Virol 2007, 88: 849-858. 10.1099/vir.0.82510-0
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: 2698-2704. 10.1128/JCM.00542-06
Rozen S, Skaletsky HJ: Primer3 on the www for general users and for biologist programmers. In Bioinformatics methods and protocols: methods in molecular biology. Edited by: Krawetz S, Misener S. Totowa, NJ: Humana Press; 2000:365-386.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215: 403-410.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011, 28: 2731-2739. 10.1093/molbev/msr121
Palacios G, Quan PL, Jabado OJ, Conlan S, Hirschberg DL, Liu Y: Panmicrobial oligonucleotide array for diagnosis of infectious diseases. Emerg Infect Dis 2007, 13: 73-81. 10.3201/eid1301.060837
454 Sequencing System Software Manual, v 2.5p1: Part C - GS De novo assembler, GS reference mapper, SFF tools. 2010. In press
FASTX-Toolkit: FASTQ/A short-reads pre-processing tools. Available from: http://hannonlab.cshl.edu/fastx_toolkit/
Buchanan L, Malcolm J: The challenge of providing child health care in the indigenous population of New Zealand. J Paediatr Child Health 2010, 46: 471-474. 10.1111/j.1440-1754.2010.01838.x
Grant CC, Wall CR, Gibbons MJ, Morton SM, Santosham M, Black RE: Child nutrition and lower respiratory tract disease burden in New Zealand: a global context for a national perspective. J Paediatr Child Health 2011, 47: 497-504. 10.1111/j.1440-1754.2010.01868.x
Affairs SNZaMoPI: Health and Pacific peoples in New Zealand. Wellington, New Zealand: Statistics New Zealand and Ministry of Pacific Island Affairs; 2011.
Tukuitonga CR, Bell S, Robinson E: Hospial admission among Pacific children Auckland 1992–97. N Z Med J 2000, 113: 358-361.
Gracey M, Gee V: Hospitalization of infants for infections in Western Australia, 1980–91. J Paediatr Child Health 1994, 30: 502-505. 10.1111/j.1440-1754.1994.tb00721.x
Moore H, Burgner D, Carville K, Jacoby P, Richmond P, Lehmann D: Diverging trends for lower respiratory infections in non-Aboriginal and Aboriginal children. J Paediatr Child Health 2007, 43: 451-457. 10.1111/j.1440-1754.2007.01110.x
Williams P, Gracey M, Smith P: Hospitalization of aboriginal and non-aboriginal patients for respiratory tract diseases in Western Australia, 1988–1993. Int J Epidemiol 1997, 26: 797-805. 10.1093/ije/26.4.797
We wish to acknowledge the help of the technical staff in the National Poliovirus and Enterovirus Reference Laboratory at the Institute of Environmental Science and Research Limited. As well, the staff at Middlemore Hospital and Waikato Hospitals for the acquisition of specimens and for providing us with the necessary clinical details.
The national poliovirus/enterovirus surveillance is funded by the Ministry of Health who kindly permitted the use of relevant data for publication. The Special Pathogens program is funded by the Institute of Environmental Science and Research Limited core fund provided by the Ministry for Business, Innovation and Employment.
The authors declare that they have no competing interest.
AKT participated in the design and co-ordination of the study, carried out the molecular PCR assays and wrote the manuscript. RJH participated in the design of the study, prepared a representative enterovirus strain for whole genome sequencing and assisted with sequence alignment. JW performed sequence alignment and phylogenetic analysis. MP participated in the design of the study and assisted with sequence alignment. SMcT performed initial molecular PCR assays. CJR carried out enterovirus sequencing. J-AS carried out enterovirus sequencing and assisted with sequence alignment and phylogenetic analysis. ST participated in study design and assisted with assay development. QSH participated in study design and contributed to writing of the final manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Todd, A.K., Hall, R.J., Wang, J. et al. Detection and whole genome sequence analysis of an enterovirus 68 cluster. Virol J 10, 103 (2013). https://doi.org/10.1186/1743-422X-10-103