Search

Community-associated, methicillin-resistant Staphylococcus aureus (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (R e ) and sustained transmission (R e > 1) coincided with spread of progenitor methicillin-susceptible S. aureus (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated S. aureus lineages from America, Asia, Australasia, and Europe. Surges in R e were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers.

Background Nasopharyngeal colonisation with nontypeable Haemophilus influenzae (NTHi) is associated with development of infections including pneumonia and otitis media. The 10-valent pneumococcal conjugate vaccine (PCV10) uses NTHi Protein D (PD) as a carrier. Papua New Guinean children have exceptionally early and dense NTHi carriage, and high rates of NTHi-associated disease. Vaccination with PCV10 could potentially reduce NTHi carriage and disease in this population by inducing a NTHi PD immune response. Methods Serum and nasopharyngeal swabs were collected from 101 Papua New Guinean children at 1, 4, 9, 10, 23 and 24 months of age. Children received PCV10 (n = 55) or PCV13 (not containing NTHi PD) (n = 46) at 1, 2 and 3 months of age. NTHi carriage density was measured in swabs by qPCR. Serum PD-IgG levels were measured by bead-based immunoassay. Results Papua New Guinean children did naturally develop PD-IgG antibodies whose levels were increased at 4 months of age with PCV10 vaccination at 1–2-3 months. Despite this, most children were colonised with NTHi by 4 months of age (~95%) regardless of being vaccinated with PCV10 or PCV13, and PCV10 had no impact on NTHi carriage density. Conclusion Early vaccination of infants with PCV10 elicited a robust PD antibody response but this had no impact on NTHi carriage. Trial registration ClinicalTrials.gov CTN NCT01619462.

Streptococcus pneumoniae is a key contributor to childhood morbidity and mortality in Papua New Guinea (PNG). For the first time, whole genome sequencing of 174 isolates has enabled detailed characterisation of diverse S. pneumoniae causing invasive disease in young children in PNG, 1989-2014. This study captures the baseline S. pneumoniae population prior to the introduction of 13-valent pneumococcal conjugate vaccine (PCV13) into the national childhood immunisation programme in 2014. Relationships amongst lineages, serotypes and antimicrobial resistance traits were characterised, and the population was viewed in the context of a global collection of isolates. The analyses highlighted adiverse S. pneumoniae population associated with invasive disease in PNG, with 45 unique Global Pneumococcal Sequence Clusters (GPSCs) observed amongst the 174 isolates reflecting multiple lineages observed in PNG that have not been identified in other geographic locations. The majority of isolates were from children with meningitis, of which 52% (n=72) expressed non-PCV13 serotypes. Over a third of isolates were predicted to be resistant to at least one antimicrobial. PCV13 serotype isolates had 10.1 times the odds of being multidrug-resistant (MDR) compared to non-vaccine serotype isolates, and no isolates with GPSCs unique to PNG were MDR. Serotype 2 was the most commonly identified serotype; we identified a highly clonal cluster of serotype 2 isolates unique to PNG, and a distinct second cluster indicative of long-distance transmission. Ongoing surveillance, including whole-genome sequencing, is needed to ascertain the impact of the national PCV13 programme upon the S. pneumoniae population, including serotype replacement and antimicrobial resistance traits.

Nanopore sequencing and phylodynamic modeling have been used to reconstruct the transmission dynamics of viral epidemics, but their application to bacterial pathogens has remained challenging. Cost-effective bacterial genome sequencing and variant calling on nanopore platforms would greatly enhance surveillance and outbreak response in communities without access to sequencing infrastructure. Here, we adapt random forest models for single nucleotide polymorphism (SNP) polishing developed by Sanderson and colleagues (2020. High precision Neisseria gonorrhoeae variant and antimicrobial resistance calling from metagenomic nanopore sequencing. Genome Res. 30(9):1354–1363) to estimate divergence and effective reproduction numbers (Re) of two methicillin-resistant Staphylococcus aureus (MRSA) outbreaks from remote communities in Far North Queensland and Papua New Guinea (PNG; n = 159). Successive barcoded panels of S. aureus isolates (2 × 12 per MinION) sequenced at low coverage (>5× to 10×) provided sufficient data to accurately infer genotypes with high recall when compared with Illumina references. Random forest models achieved high resolution on ST93 outbreak sequence types (>90% accuracy and precision) and enabled phylodynamic inference of epidemiological parameters using birth–death skyline models. Our method reproduced phylogenetic topology, origin of the outbreaks, and indications of epidemic growth (Re > 1). Nextflow pipelines implement SNP polisher training, evaluation, and outbreak alignments, enabling reconstruction of within-lineage transmission dynamics for infection control of bacterial disease outbreaks on portable nanopore platforms. Our study shows that nanopore technology can be used for bacterial outbreak reconstruction at competitive costs, providing opportunities for infection control in hospitals and communities without access to sequencing infrastructure, such as in remote northern Australia and PNG.

BackgroundTyphoid fever, a systemic infection caused by Salmonella enterica serovar Typhi, remains a considerable public health threat in impoverished regions within many low- and middle-income settings. However, we still lack a detailed understanding of the emergence, population structure, molecular mechanisms of antimicrobial resistance (AMR), and transmission dynamics of S. Typhi across many settings, particularly throughout the Asia-Pacific islands. Here we present a comprehensive whole genome sequence (WGS) based overview of S. Typhi populations circulating in Papua New Guinea (PNG) over 30 years.Principle findingsBioinformatic analysis of 86 S. Typhi isolates collected between 1980-2010 demonstrated that the population structure of PNG is dominated by a single genotype (2.1.7) that appears to have emerged in the Indonesian archipelago in the mid-twentieth century with very limited evidence of inter-country transmission. Genotypic and phenotypic data demonstrated that the PNG S. Typhi population appears to be susceptible to former first line drugs for treating typhoid fever (chloramphenicol, ampicillin and co-trimoxazole), as well as fluoroquinolones, third generation cephalosporins, and macrolides. PNG genotype 2.1.7 was genetically conserved, with very few deletions, and no evidence of plasmid or prophage acquisition. Genetic variation among this population was attributed to either single point mutations, or homologous recombination adjacent to repetitive ribosomal RNA operons.SignificanceAntimicrobials remain an effective option for the treatment of typhoid fever in PNG, along with other intervention strategies including improvements to water, sanitation and hygiene (WaSH) related infrastructure and potentially the introduction of Vi-conjugate vaccines. However, continued genomic surveillance is warranted to monitor for the emergence of AMR within local populations, or the introduction of AMR associated genotypes of S. Typhi in this setting.Author SummaryTyphoid fever, caused by Salmonella enterica serovar Typhi, is a systemic infection common to many low- to middle-income settings. While the population structure of S. Typhi has been genetically characterised using whole genome sequencing in many endemic countries throughout Sub-Saharan Africa and South Asia, we are still lacking a detailed understanding for many regions including those among the Asia-Pacific islands. Genomic surveillance of isolates spanning 30 years demonstrated a population structure of S. Typhi in Papua New Guinea (PNG) dominated by a single genotype (2.1.7) that emerged in the mid-twentieth century, is genetically homogeneous, and sensitive to a wide range of antibiotics commonly used in the treatment of typhoid. There was little evidence of inter-country transmission and the setting appeared free of S. Typhi genotypes commonly associated with AMR e.g. H58 (genotype 4.3.1). These data suggest that former first line drugs (chloramphenicol, ampicillin and co-trimoxazole), fluoroquinolones, third generation cephalosporins and macrolides all remain viable options for controlling typhoid in addition to the introduction of Vi-conjugate vaccines and improvements to water, sanitation and hygiene (WaSH) related infrastructure. Routine molecular surveillance is necessary to monitor for introduced or emerging AMR to inform treatment guidelines and intervention strategies.

BackgroundDevelopment of vaccines to prevent disease and death from Streptococcus pneumoniae, and nontypeable Haemophilus influenzae (NTHi), the main pathogens that cause otitis media, pneumonia, meningitis and sepsis, are a global priority. Children living in low and lower-middle income settings are at the highest risk of contracting and dying from these diseases. Improved vaccines with broader coverage are required. Data on the natural development of antibodies to putative vaccine antigens, especially in high-risk settings, can inform the rational selection of the best antigens for vaccine development.MethodsSerum IgG titres to four pneumococcal proteins (PspA1, PspA2, CbpA, and Ply) and five NTHi antigens (P4, P6, OMP26, rsPilA and ChimV4) were measured in sera collected from 101 Papua New Guinean children at 1, 4, 9, 10, 23 and 24 months of age using multiplexed bead-based immunoassays. Carriage density of S. pneumoniae and H. influenzae were assessed by quantitative PCR on genomic DNA extracted from nasopharyngeal swabs using species-specific primers and probes. All data were log-transformed for analysis using Student’s unpaired t-tests with geometric mean titre (GMT) or density (GMD) calculated with 95% confidence intervals (CI).ResultsSerum -pneumococcal protein-specific IgG titres followed a “U” shaped pattern, with a decrease in presumably maternally-derived IgG titres between 1 and 4 months of age and returning to similar levels as those measured at 1 month of age by 24 months of age. In contrast, NTHi protein-specific IgG titres steadily increased with age. There was no correlation between antibody titres and carriage density for either pathogen.ConclusionThis longitudinal study indicates that the waning of maternally- derived antibodies that is usually observed in infants, after infants does not occur for NTHi antigens in Papua New Guinean infants. Whether NTHi antigen IgG can be transferred maternally remains to be determined. Vaccines that are designed to specifically increase the presence of protective NTHi antibodies in the first few months of life may be most effective in reducing NTHi disease.Clinical Trial Registrationhttps://clinicaltrials.gov/, identifier NCT01619462.

Nanopore sequencing and phylodynamic modelling have been used to reconstruct the transmission dynamics of viral epidemics, but their application to bacterial pathogens has remained challenging. Here, we implement Random Forest models for single nucleotide polymorphism (SNP) polishing to estimate divergence and effective reproduction numbers (Re) of two community-associated, methicillin-resistant Staphylococcus aureus (MRSA) outbreaks in remote Far North Queensland and Papua New Guinea (n = 159). Successive bar-coded panels of S. aureus isolates (2 × 12 per MinION) sequenced at low-coverage (> 5x - 10x) provided sufficient data to accurately infer assembly genotypes with high recall when compared with Illumina references. De novo SNP calling with Clair was followed by SNP polishing using intra- and inter-species models trained on Snippy reference calls. Models achieved sufficient resolution on ST93 outbreak sequence types (> 70 - 90% accuracy and precision) for phylodynamic modelling from lineage-wide hybrid alignments and birth-death skyline models in BEAST2. Our method reproduced phylogenetic topology, geographical source of the outbreaks, and indications of sustained transmission (Re > 1). We provide Nextflow pipelines that implement SNP polisher training, evaluation, and outbreak alignments, enabling reconstruction of within-lineage transmission dynamics for infection control of bacterial disease outbreaks using nanopore sequencing.

Community-associated, methicillin-resistantStaphylococcus aureus(MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 years. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole genome sequencing data for the Australian sequence type (ST) 93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (Re) and sustained transmission (Re> 1) coincided with spread of progenitor methicillin-susceptibleS. aureus(MSSA) in remote northern Australia, dissemination of the ST93-MRSA-IV geno-type into population centers on the Australian East Coast, and sub-sequent importation into the highlands of Papua New Guinea and Far North Queensland. Analysis of a ST772-MRSA-V cluster in Pakistan suggests that sustained transmission in the community following importation of resistant genotypes may be more common than previously thought. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associatedS. aureuslineages from America, Asia, Australasia and Europe. Surges in Rewere observed at the divergence of antibiotic resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed amongst drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA and MSSA lineages in the late 20th century was driven by a combination of antibiotic resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers.

Shigella is a common cause of diarrhoea in Papua New Guinea (PNG) and other Oceania countries. However, little is known about the strains causing infection. Archived Shigella isolates (n = 72) were obtained from research laboratories in PNG and reference laboratories in Australia. Shigella virulence genes were detected by PCR, and antimicrobial susceptibility was determined by disk diffusion. The ipaH virulence gene was present in all 72 isolates. The prevalence of other virulence genes was variable, with ial, invE, ipaBCD, sen/ospD3 and virF present in 60% of isolates and set1A and set1B genes present in 42% of isolates. Most S. flexneri isolates contained genes encoding enterotoxin 1 and/or enterotoxin 2. Resistance to antibiotics was common, with 51/72 isolates resistant to 2-4 antimicrobials. A greater proportion of bacteria isolated since 2010 (relative to pre-2010 isolates) were resistant to commonly used antibiotics such as ampicillin, chloramphenicol, tetracycline, and trimethoprim-sulfamethoxazole; suggesting that antimicrobial resistance (AMR) in Shigella is increasing over time in the Oceania region. There is a need for improved knowledge regarding Shigella circulation in the Oceania region and further monitoring of AMR patterns.

Background Bacterial meningitis remains an important infection globally, with the greatest burden in children in low-income settings, including Papua New Guinea (PNG). We present serotype, antimicrobial susceptibility and outcome data from paediatric meningitis patients prior to introduction of Haemophilus influenzae type b (Hib) and pneumococcal conjugate vaccines (PCVs) in PNG, providing a baseline for evaluation of immunisation programs. Methods Cerebrospinal fluid (CSF) was collected from children admitted to Goroka General Hospital with suspected meningitis between 1996 and 2005. Culture and sensitivity was conducted, and pneumococci and H. influenzae were serotyped. Laboratory findings were linked to clinical outcomes. Results We enrolled 1884 children. A recognised pathogen was identified in 375 children (19.9 %). Streptococcus pneumoniae (n = 180) and Hib (n = 153) accounted for 88.8 % of pathogens isolated. 24 different pneumococcal serogroups were identified; non-PCV types 2, 24 and 46 accounted for 31.6 % of pneumococcal meningitis. 10- and 13-valent PCVs would cover 44.1 % and 45.4 % of pneumococcal meningitis respectively. Pneumococcal isolates were commonly resistant to penicillin (21.5 %) and 23 % of Hib isolates were simultaneously resistant to ampicillin, co-trimoxazole and chloramphenicol. The case fatality rate in patients with a recognised bacterial pathogen was 13.4 % compared to 8.5 % in culture-negative patients. Conclusions If implemented in routine expanded programme of immunisation (EPI) with high coverage, current PCVs could prevent almost half of pneumococcal meningitis cases. Given the diversity of circulating serotypes in PNG serotype replacement is of concern. Ongoing surveillance is imperative to monitor the impact of vaccines. In the longer term vaccines providing broader protection against pneumococcal meningitis will be needed. Electronic supplementary material The online version of this article (doi:10.1186/s12879-015-1197-0) contains supplementary material, which is available to authorized users.

Overview of microscopy results for cases of suspected meningitis in children in Goroka, Papua New Guinea. In this study, unusually high rates of S. aureus positive CSF was detected. Microscopy results support the notion that S. aureus is likely a contaminant in the majority of samples from which it was isolated. (DOCX 15Â kb)

Prevalence of serogroups and serotypes of S. pneumoniae isolated from children with meningitis admitted to Goroka General Hospital. (DOCX 22Â kb)

Comparison of number of polymorphonucleocytes in CSF samples positive for S. aureus isolation versus recognised pathogens (S. pneumoniae and H. influenzae), other pathogens (non-pneumococcus, non-Hi), probable contaminants and samples from which no bacteria were isolated. There is no evidence of significance difference between PMN numbers in CSF with S. aureus compared to CSF with no pathogens isolated. Analysis conducted using a non-parametric independent samples median test. (DOCX 16Â kb)