Email updates

Keep up to date with the latest news and content from Virology Journal and BioMed Central.

Open Access Research

Detection and whole genome sequence analysis of an enterovirus 68 cluster

Angela K Todd1*, Richard J Hall1, Jing Wang1, Mathew Peacey1, Sharla McTavish1, Christy J Rand2, Jo-Ann Stanton2, Susan Taylor3 and Q Sue Huang1

Author Affiliations

1 Institute of Environmental Science and Research Limited, National Centre for Biosecurity and Infectious Disease, 66 Ward Street, Wallaceville, Upper Hutt, 5018, New Zealand

2 Department of Anatomy, University of Otago, P. O. Box 913, Dunedin, 9054, New Zealand

3 Middlemore Hospital, Otahuhu, Auckland, 2025, New Zealand

For all author emails, please log on.

Virology Journal 2013, 10:103  doi:10.1186/1743-422X-10-103


The electronic version of this article is the complete one and can be found online at: http://www.virologyj.com/content/10/1/103


Received:2 October 2012
Accepted:20 March 2013
Published:2 April 2013

© 2013 Todd et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

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.

Objectives

To genetically characterise enterovirus 68 specimens identified in New Zealand in 2010.

Study design

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.

Results

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.

Conclusions

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.

Keywords:
Enterovirus 68; Respiratory tract infections; Molecular genome sequencing

Background

Enteroviruses belong to the Picornaviridae family and are among the most commonly identified aetiological agents of human disease [1]. 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 [3]. Since then, EV68 has been isolated rarely; only 26 strains have been identified over 36 years in the USA [2]. 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].

Objectives

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.

Study design

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.

RNA extractions

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 [17]. An EV68 RT-PCR assay was designed in-house using Primer3 v 0.4.0 [18] 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).

Phylogenetic analysis

Sequences were edited and assembled using Bionumerics 5.10 (Applied Maths, Belgium) and compared to known VP1 sequences in GenBank using BLAST [19]. 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 [20] (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 [21]. 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 [22] and FASTX-Toolkit [23] 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.

Results

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).

Table 1. Characteristics of fifteen patients in whom enterovirus 68 was isolated in New Zealand in 2010

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 [17]. 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.

thumbnailFigure 1. The evolutionary history was 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 MEGA520 (Center for Evolutionary Medicine and Informatics, Tempe, AZ, USA). 1000 bootstrap replicates were performed to determine the consensus tree presented in this figure, support for nodes present in greater than 75% of the trees are annotated. The viruses detected in this study are marked with an asterisk. The GenBank accession number and the strain designation is shown for each virus.

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.

Table 2. Pairwise comparison of nucleotide sequence identities of NZ-2010-541 with other EV68 strains (shown in %)

Discussion

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 [2]. 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 [5]. 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 [12]. 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 [12]. 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 [12]. 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.

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

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.

Acknowledgements

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.

Funding

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.

References

  1. 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. OpenURL

  2. Khetsuriani N, Lamonte-Fowlkes A, Oberst S, Pallansch MA: Enterovirus surveillance–United States, 1970–2005.

    MMWR Surveill Summ 2006, 55:1-20. PubMed Abstract | Publisher Full Text OpenURL

  3. Schieble JH, Fox VL, Lennette EH: A probable new human picornavirus associated with respiratory diseases.

    Am J Epidemiol 1967, 85:297-310. PubMed Abstract | Publisher Full Text OpenURL

  4. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. 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. PubMed Abstract | Publisher Full Text OpenURL

  6. 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. PubMed Abstract | Publisher Full Text OpenURL

  7. 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. PubMed Abstract | Publisher Full Text OpenURL

  8. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. 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. OpenURL

  10. 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. PubMed Abstract | Publisher Full Text OpenURL

  11. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. 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. PubMed Abstract | Publisher Full Text OpenURL

  13. 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. PubMed Abstract | Publisher Full Text OpenURL

  14. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. 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. PubMed Abstract | Publisher Full Text OpenURL

  17. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. 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. OpenURL

  19. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215:403-410. PubMed Abstract OpenURL

  20. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. 454 Sequencing System Software Manual, v 2.5p1:

    Part C - GS De novo assembler, GS reference mapper, SFF tools. 2010.

    In press

    OpenURL

  23. FASTX-Toolkit:

    FASTQ/A short-reads pre-processing tools.

    Available from: http://hannonlab.cshl.edu/fastx_toolkit/ webcite

    OpenURL

  24. 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. PubMed Abstract | Publisher Full Text OpenURL

  25. 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. PubMed Abstract | Publisher Full Text OpenURL

  26. Affairs SNZaMoPI: Health and Pacific peoples in New Zealand. Wellington, New Zealand: Statistics New Zealand and Ministry of Pacific Island Affairs; 2011. OpenURL

  27. Tukuitonga CR, Bell S, Robinson E: Hospial admission among Pacific children Auckland 1992–97.

    N Z Med J 2000, 113:358-361. PubMed Abstract OpenURL

  28. Gracey M, Gee V: Hospitalization of infants for infections in Western Australia, 1980–91.

    J Paediatr Child Health 1994, 30:502-505. PubMed Abstract | Publisher Full Text OpenURL

  29. 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. PubMed Abstract | Publisher Full Text OpenURL

  30. 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. PubMed Abstract | Publisher Full Text OpenURL