Geographic consistency in dominant, non-typeable Haemophilus influenzae genotypes colonising four distinct Australian paediatric groups: a cohort study
- Heidi C. Smith-Vaughan1, 2Email author,
- Jemima Beissbarth1,
- Jacinta Bowman3,
- Kim M. Hare1,
- Erin P. Price1,
- Janessa Pickering4, 5,
- Deborah Lehmann4,
- Anne B. Chang1, 6,
- Peter S. Morris1, 7,
- Robyn L. Marsh†1 and
- Amanda J. Leach†1
© The Author(s) 2016
Received: 14 December 2015
Accepted: 8 August 2016
Published: 18 August 2016
Non-typeable Haemophilus influenzae (NTHi)-associated ear and respiratory diseases (including pneumonia) represent a major health burden in many parts of the world. NTHi strains retrieved from the upper airways commonly reflect those found in the lower airways. Despite growing genomic and genotyping data on NTHi, there remains a limited understanding of global and regional NTHi population structures. The aim of this study was to determine whether nasopharyngeal carriage in four Australian paediatric groups at varying risk of NTHi colonisation was dominated by the same NTHi genotypes. Genotyping data generated by PCR-ribotyping were evaluated for 3070 NTHi isolates colonising the nasopharynges of Aboriginal and non-Aboriginal children enrolled in four longitudinal studies in three separate urban and remote regions of Australia. Several NTHi PCR-ribotypes dominated in nasopharyngeal carriage, irrespective of study setting. Principal coordinates analysis confirmed a cluster of common PCR-ribotypes among all cohorts. In conclusion, we identified dominant PCR-ribotypes common to geographically disparate Australian paediatric populations. Future genomic analyses will shed further light on the precise factors underlying the dominance of certain NTHi strains in nasopharyngeal carriage.
KeywordsNon-typeable Haemophilus influenzae Genotyping Carriage Paediatric PCR-ribotyping
Diseases associated with non-typeable Haemophilus influenzae (NTHi) represent a major health burden worldwide . This Gram-negative bacterium is commonly associated with adult community-acquired pneumonia (CAP), particularly among those with underlying respiratory disease . Because young children generally cannot expectorate, the aetiology of paediatric CAP is poorly understood and a role for NTHi remains controversial. NTHi is also associated with chronic lung diseases; for example, it is the most commonly detected pathogen in chronic obstructive pulmonary disease, a disease identified by the World Health Organization as the third leading cause of mortality globally . NTHi is the most common bacterium in the lower airways of adults and children with non-cystic fibrosis bronchiectasis [4, 5], being cultured from the sputa of up to 70 % of adults , and 78 % of bronchoalveolar lavage specimens from Australian Indigenous children with bronchiectasis . NTHi is also a major otitis media pathogen and has been detected in 89 % of ear discharge specimens from children with suppurating ears in remote Australian Indigenous communities .
Genetic and genomic studies have shown that NTHi is a genetically diverse organism that frequently undergoes lateral gene transfer [7, 8], a trait that has hindered our understanding of the population biology of this organism. As an example of this diversity, the H. influenzae multilocus sequence typing (MLST) database (http://pubmlst.org/hinfluenzae/) currently lists 1369 NTHi isolates with MLST data from across the globe, which are represented by 966 distinct sequence types (STs). Similarly, our studies of NTHi nasopharyngeal carriage isolates have identified 130 PCR-ribotypes (PRTs) to date. Despite this high level of diversity, several studies have identified a population structure among NTHi [7, 9]. Of note, MLST and whole genome analyses show some correlation between population structure and virulence potential, although no geographic clustering was detected [7, 9]. While these data contribute to our understanding of NTHi populations, large genomic studies are necessary to improve our understanding of geographic patterns, temporal changes (including outbreaks), adaptive evolutionary pressures such as those driven by antibiotic or vaccine selection, and tissue tropism. Whole genome sequencing is ideal for exploration of bacterial population structure, and that is the direction for future studies; however, in the current study we used PRT. In addition to its cost effectiveness, PRT has the advantage of interrogating 16S rDNA, which is an evolutionarily stable marker across bacterial species. This slow rate of evolution potentially provides a more robust signal for detecting strain relatedness than other genotyping methods that are based on non-ribosomal loci. We compared PRT data for 3070 nasopharyngeal carriage NTHi isolates obtained from four paediatric groups from three distinct geographic regions of Australia. Our aim was to determine whether the same NTHi PRTs dominated in nasopharyngeal carriage in these geographically distinct populations.
This was a retrospective cohort study that analysed NTHi genotyping data generated previously in completed studies.
Participants and setting
Description of the 4 Australian paediatric non-typeable Haemophilus influenzae (NTHi)-carriage cohorts used in this study
Study design, setting and participant ethnicity
Year of data
No. enrolled children
(NTHi carriage rate)
Child eligibility age & other key criteriab
(No. NTHi positive swabs)
No. unique NTHi typeda
(Total NTHi typed)
(No. per 100 swabs)
Cluster RCT of a hygiene intervention in 20 Darwin (NT) child care centres.
90 % were non-Aboriginal 
0–4 yrs, attending 3 d/wk
Fortnightly for 6 m
RCT of azithromycin vs amoxicillin for AOM in 16 remote NT communities. All were Aboriginal children 
(316 with swab data)
0.5–6 yrs with AOM
Days 0 and 6–11. Additionally, Day 12–21 if perforation.
1 wk to 24 m
Ages 1–3, 6–8 wks, then months 4, 6, 12, 18, 24.
1 wk to 24 m
Ages 1–3, 6–8 wks, then months 4, 6, 12, 18, 24.
Laboratory and statistical analyses
For this study, we used genotyping data generated by PRT, a method that interrogates genetic polymorphisms in the H. influenzae ribosomal operons using restriction digestion . Briefly, approximately 6 kb ribosomal operons (16S–5S) were amplified by long-PCR, and HaeIII restriction fragments of the resulting product were separated by agarose gel electrophoresis. PRTs were assigned based on restriction fragment lengths (0.4 to 1.2 kb) following alignment to a 2-log DNA ladder (New England BioLabs, United States) and an in-house standard. To aid pattern recognition, the images were divided into 16 predetermined sections.
Results and discussion
In the Northern Territory paediatric cohorts (Study 1 and Study 2), the mean NTHi carriage rates across all swabs were approximately 50 % (unpublished) for Study 1 and 85 % (269/316) at baseline for Study 2  (Table 1). In the Kalgoorlie studies, NTHi carriage rates among the Aboriginal (Study 3) and non-Aboriginal (Study 4) children were 36.3 % (183/504) and 9.4 % (98/1045), respectively . In Study 1, 84 different PRTs were identified among 2179 NTHi isolates. In Study 2, 73 different PRTs were identified among 551 NTHi isolates. In Study 3, 65 different PRTs were detected among 231 NTHi isolates, and 37 PRTs were detected among 109 isolates in Study 4.
Principal coordinates analysis (PCoA) was used to visualise the PRT distribution among the four study cohorts. This analysis showed a dominant cluster of PRTs common to all studies (Fig. 1), thereby supporting the absence of geographic or population-specific clustering of NTHi PRTs. The highest number of PRTs was detected in Study 1. Higher dispersion of PRTs across PCO1 and PCO2 for Studies 1–2 (compared to Studies 3–4) reflects the increased number of PRTs unique to Studies 1–2. For Study 1, this finding likely represents the large sample (2179 NTHi isolates). However, for Study 2, this finding likely reflects the larger number of strains dominating in this population, which may result from the high NTHi carriage rates (85 %) in this cohort or its broader geographic sampling (16 remote communities). This finding could also reflect a difference in disease state (all children had a diagnosis of AOM), though we previously showed that NTHi from paired nasopharyngeal and ear discharge swabs from children with AOM with perforation were generally the same PRT (unpublished observations).
Another carriage study design question relates to the number of isolates selected for analysis from each swab. Multiple NTHi genotypes during carriage have previously been observed and modelled ; however, the value of assessing multiple strain carriage in different populations is not well understood. In Studies 1–2 where additional colonies were selected only if morphologically distinct, multiple PRT carriage was detected in 8 and 17 % of swabs, respectively. For Studies 3–4 where 2 colonies were selected as per protocol, including any that were morphologically distinct, multiple PRT carriage was detected in 20 and 7 % of swabs, respectively. Although our data cannot inform the most appropriate number of colonies to select per specimen for measuring multiple PRT carriage, nor the method of selection, it is likely that these methods would vary by population. Indeed, the value of typing additional colonies appears to be greater for populations with high carriage rates.
In conclusion, analysis of NTHi nasopharyngeal carriage PRTs in four cohorts with varying NTHi colonisation rates allowed us to identify dominant genotypes common to four Australian paediatric populations. More than 3000 isolates from three distinct urban and remote geographic regions in Australia (Kalgoorlie in southern Western Australia, the urban Darwin region of the Northern Territory, and several remote communities across the Northern Territory) were examined. This finding provides an opportunity to study the factors conferring an advantage to the dominant PRTs in NTHi nasopharyngeal carriage, thereby advancing our understanding of the constraints that support NTHi population structure. Continued efforts to understand NTHi population dynamics will improve the identification of useful vaccine candidates for targeting universal or dominant strains, and may inform models of the potential efficacy of NTHi vaccines in current and future production pipelines. Further characterisation of carriage- and disease-related (acute [e.g. pneumonia] and chronic [e.g. bronchiectasis]) NTHi populations by whole genome sequencing or other genotyping means is necessary to determine the relevance of carriage versus disease genotype frequency, and may shed light on the factors supporting dominance of certain NTHi genotypes. Although the Australian experience may not necessarily apply elsewhere, whole genome sequence data from Australian and international NTHi isolates have (to date) demonstrated a level of core genome similarity between global and Australian NTHi populations . Thus, we expect our findings to be applicable to other study settings.
We thank the families who participated in these studies for making this research possible. We also thank the many coordinators, nurses and microbiologists who contributed to the original studies, which provided the specimens used in this particular study.
The original studies or specimen analyses were funded through the Australian National Health and Medical Research Council (NHMRC grant numbers 212044, 436023, 100009, 1023781), and Studies 3–4 through two Healthway grants (grant numbers 6028 and 10564). The Darwin component of the paper is associated with an NHMRC Centre for Research Excellence in Otitis Media of Aboriginal and Torres Strait Islander Children (grant number 1078557) and an NHMRC Centre for Research Excellence in Lung Health of Aboriginal and Torres Strait Islander Children (grant number 1040830). RM is supported by NHMRC Frank Fenner Early Career Fellowship 1034703. DL was funded through NHMRC program (grant number 353514). AC is supported by NHMRC Practitioner Fellowship 545216. AL is supported by NHMRC Elizabeth Blackburn Research Fellowship 1020561. KH is supported by NHMRC Early Career Fellowship 1072870. The funders had no role in study design, collection and analysis of data, decision to publish, or writing of the manuscript.
Availability of data and materials
Data and materials that are proprietary to Telethon Kids Institute and Menzies School of Health Research under the custodianship of the study’s principal investigator can be requested from the Custodian or relevant Custodial Committees; and is subject to permission from the relevant Human Research Ethics Committees and existing contractual and third party agreements.
All the authors met ICMJE authorship criteria. JB, JP, DL, KH, PM, EP, AC, AL contributed to data acquisition. HS-V, JB, RM conducted the data analysis and interpretation. HS-V wrote the first draft of the manuscript; JB, JP, DL, KH, PM, AL, EP, AC critically reviewed the manuscript. All authors advised on intellectual content of the manuscript. All authors agree with the manuscript’s results and conclusions. All authors approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Ethical approvals for Studies 1 (#99/68) and 2 (#02–44) were obtained from the Human Research Ethics Committees of the Menzies School of Health Research and the Northern Territory Department of Health, and from the Central Australian Human Research Ethics Committee. Study 2 was registered with the Australian Clinical Trials Registry (ACTRN 12609000691246). Ethical approval to conduct Studies 3 and 4 was provided by the Western Australian Aboriginal Health Information and Ethics Committee (WAAHIEC), the Northern Goldfields Health Service and Nursing Education Ethics Committee in Kalgoorlie, Princess Margaret Hospital Ethics Committee, and the Confidentiality of Health Information Committee of the Health Department of Western Australia. Written informed consent was obtained from families to enrol their child in these studies.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Van Eldere J, Slack MP, Ladhani S, Cripps AW. Non-typeable Haemophilus influenzae, an under-recognised pathogen. Lancet Infect Dis. 2014;14:1281–92. PMID:25012226, http://dx.doi.org/10.1016/S1473-3099(14)70734–0.View ArticlePubMedGoogle Scholar
- Slack MP. A review of the role of Haemophilus influenzae in community-acquired pneumonia. Pneumonia. 2015;29:26–43. http://dx.doi.org/10.15172/pneu.2015.6/520.View ArticleGoogle Scholar
- WHO. The Top 10 Causes of Death. 2014. Available from: http://www.who.int/mediacentre/factsheets/fs310/en/. Accessed 5 Jan 2015.Google Scholar
- Hare KM, Grimwood K, Leach AJ, Smith-Vaughan H, Torzillo PJ, Morris PS, et al. Respiratory bacterial pathogens in the nasopharynx and lower airways of Australian indigenous children with bronchiectasis. J Pediatr. 2010;157:1001–5. PMID:20656297, http://dx.doi.org/10.1016/j.jpeds.2010.06.002.View ArticlePubMedGoogle Scholar
- King PT, Holdsworth SR, Freezer NJ, Villanueva E, Holmes PW. Microbiologic follow-up study in adult bronchiectasis. Respir Med. 2007;101:1633–8. PMID:17467966, http://dx.doi.org/10.1016/j.rmed.2007.03.009.View ArticlePubMedGoogle Scholar
- Smith-Vaughan HC, Binks MJ, Marsh RL, Kaestli M, Ward L, Hare KM, et al. Dominance of Haemophilus influenzae in ear discharge from Indigenous Australian children with acute otitis media with tympanic membrane perforation. BMC Ear Nose Throat Disord. 2013;13:12–3. PMID:24099576, http://dx.doi.org/10.1186/1472-6815-13-12.View ArticlePubMedPubMed CentralGoogle Scholar
- De Chiara M, Hood D, Muzzi A, Pickard DJ, Perkins T, Pizza M, et al. Genome sequencing of disease and carriage isolates of nontypeable Haemophilus influenzae identifies discrete population structure. Proc Natl Acad Sci USA. 2014;111:5439–44. PMID:24706866, http://dx.doi.org/10.1073/pnas.1403353111.View ArticlePubMedPubMed CentralGoogle Scholar
- Price EP, Sarovich DS, Nosworthy E, Beissbarth J, Marsh RL, Pickering J, et al. Haemophilus influenzae: using comparative genomics to accurately identify a highly recombinogenic human pathogen. BMC Genomics. 2015;16:641. PMID:26311542, http://dx.doi.org/10.1186/s12864-015-1857-x.View ArticlePubMedPubMed CentralGoogle Scholar
- LaCross NC, Marrs CF, Gilsdorf JR. Population structure in nontypeable Haemophilus influenzae. Infect Genet Evol. 2013;14:125–36. PMID:23266487 http://dx.doi.org/10.1016/j.meegid.2012.11.023.View ArticlePubMedGoogle Scholar
- Bailie R, Morris P, Leach A, Wilson C, Dunne K, Si D, et al. The effect of additional training in recommended hygiene practices on rates of bacterial cross-infection and respiratory illness in Australian child care centres: A randomised controlled trial. Epidemiology. 2003;14:S15. http://dx.doi.org/10.1097/00001648-200309001-00009.View ArticleGoogle Scholar
- Morris PS, Gadil G, McCallum GB, Wilson CA, Smith-Vaughan HC, Torzillo P, et al. Single-dose azithromycin versus seven days of amoxycillin in the treatment of acute otitis media in Aboriginal children (AATAAC): a double blind, randomised controlled trial. Med J Aust. 2010;192:24–9. PMID:20047544.PubMedGoogle Scholar
- Watson K, Carville K, Bowman J, Jacoby P, Riley TV, Leach AJ, et al. Kalgoorlie Otitis Media Research Project Team. Upper respiratory tract bacterial carriage in Aboriginal and non-Aboriginal children in a semi-arid area of Western Australia. Pediatr Infect Dis J. 2006;25:782–90. PMID:16940834, http://dx.doi.org/10.1097/01.inf.0000232705.49634.68.View ArticlePubMedGoogle Scholar
- Smith-Vaughan HC, Sriprakash KS, Mathews JD, Kemp DJ. Long PCR-ribotyping of nontypeable Haemophilus influenzae. J Clin Microbiol. 1995;33:1192–5. PMID:7542264.PubMedPubMed CentralGoogle Scholar
- Pickering J, Smith-Vaughan H, Beissbarth J, Bowman JM, Wiertsema S, Riley TV, et al. Diversity of nontypeable Haemophilus influenzae strains colonizing Australian Aboriginal and non-Aboriginal children. J Clin Microbiol. 2014;52:1352–7. PMID:24501028, http://dx.doi.org/10.1128/JCM.03448-13.View ArticlePubMedPubMed CentralGoogle Scholar
- Binks MJ, Temple B, Kirkham LA, Wiertsema SP, Dunne EM, Richmond PC, et al. Molecular surveillance of true nontypeable Haemophilus influenzae: an evaluation of PCR screening assays. PLoS ONE. 2012;7:e34083. PMID:22470516, http://dx.doi.org/10.1371/journal.pone.0034083.View ArticlePubMedPubMed CentralGoogle Scholar
- Beissbarth J, Binks M, Scarff K, Marsh R, Nosworthy E, Chang AB, Leach AJ, Smith-Vaughan HC. Nontypeable Haemophilus influenzae and H. haemolyticus identification in the Northern Territory of Australia. Presented at: OMOZ, 2014; Melbourne, Australia.Google Scholar
- Clarke KG, R. PRIMER v6: User Manual/Tutorial. Plymouth: PRIMER E2006Google Scholar
- Smith-Vaughan HC, McBroom J, Mathews JD. Modelling of endemic carriage of Haemophilus influenzae in Aboriginal infants in Northern Australia. FEMS Immunol Med Microbiol. 2001;31:137–43. PMID:11549421, http://dx.doi.org/10.1111/j.1574-695X.2001.tb00510.x.View ArticlePubMedGoogle Scholar