Overlapping Streptococcus pyogenes and Streptococcus dysgalactiae subspecies equisimilis household transmission and mobile genetic element exchange.

Streptococcus dysgalactiae subspecies equisimilis (SDSE) and Streptococcus pyogenes share skin and throat niches with extensive genomic homology and horizontal gene transfer (HGT) possibly underlying shared disease phenotypes. It is unknown if cross-species transmission interaction occurs. Here, we conduct a genomic analysis of a longitudinal household survey in remote Australian First Nations communities for patterns of cross-species transmission interaction and HGT. Collected from 4547 person-consultations, we analyse 294 SDSE and 315 S. pyogenes genomes. We find SDSE and S. pyogenes transmission intersects extensively among households and show that patterns of co-occurrence and transmission links are consistent with independent transmission without inter-species interference. We identify at least one of three near-identical cross-species mobile genetic elements (MGEs) carrying antimicrobial resistance or streptodornase virulence genes in 55 (19%) SDSE and 23 (7%) S. pyogenes isolates. These findings demonstrate co-circulation of both pathogens and HGT in communities with a high burden of streptococcal disease, supporting a need to integrate SDSE and S. pyogenes surveillance and control efforts.

Evaluating the role of asymptomatic throat carriage of Streptococcus pyogenes in impetigo transmission in remote Aboriginal communities in Northern Territory, Australia: a retrospective genomic analysis.

BACKGROUND: Streptococcus pyogenes, or group A Streptococcus (GAS), infections contribute to a high burden of disease in Aboriginal Australians, causing skin infections and immune sequelae such as rheumatic heart disease. Controlling skin infections in these populations has proven difficult, with transmission dynamics being poorly understood. We aimed to identify the relative contributions of impetigo and asymptomatic throat carriage to GAS transmission. METHODS: In this genomic analysis, we retrospectively applied whole genome sequencing to GAS isolates that were collected as part of an impetigo surveillance longitudinal household survey conducted in three remote Aboriginal communities in the Northern Territory of Australia between Aug 6, 2003, and June 22, 2005. We included GAS isolates from all throats and impetigo lesions of people living in two of the previously studied communities. We classified isolates into genomic lineages based on pairwise shared core genomes of more than 99% with five or fewer single nucleotide polymorphisms. We used a household network analysis of epidemiologically and genomically linked lineages to quantify the transmission of GAS within and between households. FINDINGS: We included 320 GAS isolates in our analysis: 203 (63%) from asymptomatic throat swabs and 117 (37%) from impetigo lesions. Among 64 genomic lineages (encompassing 39 emm types) we identified 264 transmission links (involving 93% of isolates), for which the probable source was asymptomatic throat carriage in 166 (63%) and impetigo lesions in 98 (37%). Links originating from impetigo cases were more frequent between households than within households. Households were infected with GAS for a mean of 57 days (SD 39 days), and once cleared, reinfected 62 days (SD 40 days) later. Increased household size and community presence of GAS and scabies were associated with slower clearance of GAS. INTERPRETATION: In communities with high prevalence of endemic GAS-associated skin infection, asymptomatic throat carriage is a GAS reservoir. Public health interventions such as vaccination or community infection control programmes aimed at interrupting transmission of GAS might need to include consideration of asymptomatic throat carriage. FUNDING: Australian National Health and Medical Research Council.

minSNPs: an R package for the derivation of resolution-optimised SNP sets from microbial genomic data.

Here, we present the R package, minSNPs. This is a re-development of a previously described Java application named Minimum SNPs. MinSNPs assembles resolution-optimised sets of single nucleotide polymorphisms (SNPs) from sequence alignments such as genome-wide orthologous SNP matrices. MinSNPs can derive sets of SNPs optimised for discriminating any user-defined combination of sequences from all others. Alternatively, SNP sets may be optimised to determine all sequences from all other sequences, i.e., to maximise diversity. MinSNPs encompasses functions that facilitate rapid and flexible SNP mining, and clear and comprehensive presentation of the results. The minSNPs' running time scales in a linear fashion with input data volume and the numbers of SNPs and SNPs sets specified in the output. MinSNPs was tested using a previously reported orthologous SNP matrix of Staphylococcus aureus and an orthologous SNP matrix of 3,279 genomes with 164,335 SNPs assembled from four S. aureus short read genomic data sets. MinSNPs was shown to be effective for deriving discriminatory SNP sets for potential surveillance targets and in identifying SNP sets optimised to discriminate isolates from different clonal complexes. MinSNPs was also tested with a large Plasmodium vivax orthologous SNP matrix. A set of five SNPs was derived that reliably indicated the country of origin within three south-east Asian countries. In summary, we report the capacity to assemble comprehensive SNP matrices that effectively capture microbial genomic diversity, and to rapidly and flexibly mine these entities for optimised marker sets.

Host-dependent resistance of Group A Streptococcus to sulfamethoxazole mediated by a horizontally-acquired reduced folate transporter.

Described antimicrobial resistance mechanisms enable bacteria to avoid the direct effects of antibiotics and can be monitored by in vitro susceptibility testing and genetic methods. Here we describe a mechanism of sulfamethoxazole resistance that requires a host metabolite for activity. Using a combination of in vitro evolution and metabolic rescue experiments, we identify an energy-coupling factor (ECF) transporter S component gene (thfT) that enables Group A Streptococcus to acquire extracellular reduced folate compounds. ThfT likely expands the substrate specificity of an endogenous ECF transporter to acquire reduced folate compounds directly from the host, thereby bypassing the inhibition of folate biosynthesis by sulfamethoxazole. As such, ThfT is a functional equivalent of eukaryotic folate uptake pathways that confers very high levels of resistance to sulfamethoxazole, yet remains undetectable when Group A Streptococcus is grown in the absence of reduced folates. Our study highlights the need to understand how antibiotic susceptibility of pathogens might function during infections to identify additional mechanisms of resistance and reduce ineffective antibiotic use and treatment failures, which in turn further contribute to the spread of antimicrobial resistance genes amongst bacterial pathogens.

Identification and Discrimination of Chlamydia trachomatis Ocular and Urogenital Strains and Major Phylogenetic Lineages by CtGEM Typing, A Double-Locus Genotyping Method.

CtGEM typing was developed to subdivide the bacterial species Chlamydia trachomatis on the basis of genome phylogeny and anatomical tropism. The rationale was facilitation of surveillance for ocular strains, although the method is applicable to essentially any C. trachomatis surveillance application that does not require high resolution. CtGEM is a double-locus genotyping method. The loci included in the assay were identified by computerized analysis of 65 complete genomes for resolution optimized sets of single nucleotide polymorphisms (SNPs). From this, two PCR amplifiable fragments were defined. One, rg1, is within a hypothetical gene annotated as Jali-1891 within the C. trachomatis B_Jali20 genome. The other, ofr, is within the ompA gene which encodes the major outer membrane protein. Variation in rg1 is conferred by two SNPs defining four haplotypes that exhibit concordance with genome phylogeny. Variation within ofr is more complex and allows for inference of ompA genotype, either to the level of single genotype, or group of closely related genotypes. Two CtGEM formats were developed. One is based on interrogation of the two loci by high resolution melting analysis (HRMA), and the other based on analysis of the loci by Sanger sequencing. The genotypes defined identify known ocular genotypes, discriminate known ocular genotypes from each other, discriminate the major phylogenetic lineages of the species, and discriminate all ompA genotypes with the exception of closely related variants within the genotypes H, I, J cluster. The Sanger sequencing format provides slightly more resolution that the HRMA format with respect to ompA genotype. An unusual aspect of this method is that all possible combinations of rg1 haplotype, and inferred ompA genotype(s) have been given CtGEM typing numbers. This includes types that at this time have not been shown to exist.

Author Correction: Atlas of group A streptococcal vaccine candidates compiled using large-scale comparative genomics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Concerns for efficacy of a 30-valent M-protein-based Streptococcus pyogenes vaccine in regions with high rates of rheumatic heart disease.

The prevalence of rheumatic heart disease (RHD) in the Aboriginal population of the Australian Northern Territory is high, and Streptococcus pyogenes skin infections likely contribute to this. A promising candidate S. pyogenes "30mer" vaccine is composed of 30 pharyngitis associated type-specific antigens from the S. pyogenes M protein. Cross opsonisation experiments suggest that 30mer vaccine protection may extend to non-cognate emm types. A new "emm cluster" scheme for classifying M protein is based on the full-length coding sequence, and correlates with functional and immunological properties, and anatomical tropism. Twenty-seven years of research in the Northern Territory has yielded 1810 S. pyogenes isolates with clinical and emm type data. The primary aim was to analyse these data with reference to the emm cluster scheme and cross opsonisation information, to inform estimation of 30mer vaccine efficacy in the Northern Territory. The isolates encompass 101 emm types. Variants of cluster A-C were enriched in throat isolates, and variants of emm cluster D enriched in skin isolates. Throat isolates were enriched for 30mer vaccine cognate emm types in comparison with skin isolates of which only 25% were vaccine emm types. While cross opsonisation data indicates potential for enhancing 30mer vaccine coverage, more than one third of skin isolates were within 38 emm types untested for cross opsonisation. Emm cluster D variants, in particular emm cluster D4, were not only all non-cognate with the vaccine, but were abundant and diverse, and less likely to be cross-opsonisation positive than other emm clusters. Long term persistence of many emm types in the study area was revealed. It was concluded that the 30mer vaccine efficacy in the Northern Territory will likely require both cross protection, and additional measures to elicit immunity against variants of emm cluster D.

Atlas of group A streptococcal vaccine candidates compiled using large-scale comparative genomics.

Group A Streptococcus (GAS; Streptococcus pyogenes) is a bacterial pathogen for which a commercial vaccine for humans is not available. Employing the advantages of high-throughput DNA sequencing technology to vaccine design, we have analyzed 2,083 globally sampled GAS genomes. The global GAS population structure reveals extensive genomic heterogeneity driven by homologous recombination and overlaid with high levels of accessory gene plasticity. We identified the existence of more than 290 clinically associated genomic phylogroups across 22 countries, highlighting challenges in designing vaccines of global utility. To determine vaccine candidate coverage, we investigated all of the previously described GAS candidate antigens for gene carriage and gene sequence heterogeneity. Only 15 of 28 vaccine antigen candidates were found to have both low naturally occurring sequence variation and high (>99%) coverage across this diverse GAS population. This technological platform for vaccine coverage determination is equally applicable to prospective GAS vaccine antigens identified in future studies.

CtGEM typing: Discrimination of Chlamydia trachomatis ocular and urogenital strains and major evolutionary lineages by high resolution melting analysis of two amplified DNA fragments.

Chlamydia trachomatis infects the urogenital tract (UGT) and eyes. Anatomical tropism is correlated with variation in the major outer membrane protein encoded by ompA. Strains possessing the ocular ompA variants A, B, Ba and C are typically found within the phylogenetically coherent "classical ocular lineage". However, variants B, Ba and C have also been found within three distinct strains in Australia, all associated with ocular disease in children and outside the classical ocular lineage. CtGEM genotyping is a method for detecting and discriminating ocular strains and also the major phylogenetic lineages. The rationale was facilitation of surveillance to inform responses to C. trachomatis detection in UGT specimens from young children. CtGEM typing is based on high resolution melting analysis (HRMA) of two PCR amplified fragments with high combinatorial resolving power, as defined by computerised comparison of 65 whole genomes. One fragment is from the hypothetical gene defined by Jali-1891 in the C. trachomatis B_Jali20 genome, while the other is from ompA. Twenty combinatorial CtGEM types have been shown to exist, and these encompass unique genotypes for all known ocular strains, and also delineate the TI and T2 major phylogenetic lineages, identify LGV strains and provide additional resolution beyond this. CtGEM typing and Sanger sequencing were compared with 42 C. trachomatis positive clinical specimens, and there were no disjunctions. CtGEM typing is a highly efficient method designed and tested using large scale comparative genomics. It divides C. trachomatis into clinically and biologically meaningful groups, and may have broad application in surveillance.

Contaminated fingers: a potential cause of Chlamydia trachomatis-positive urine specimens.

OBJECTIVES: The detection of an STI agent in a urogenital tract (UGT) specimen from a young child is regarded as being indicative of sexual abuse. However, the probabilities of contamination events that could conceivably lead to STI positive specimens in the absence of sexual contact are unclear. The objective was to estimate the potential for fingers that have come in contact with Chlamydia trachomatis-positive urine to detectably contaminate C. trachomatis-negative urine. METHODS: The study design was based on self-experimentation. Dilutions of C. trachomatis elementary bodies (EBs) were prepared. A participant contacted an EB dilution then a urine surrogate specimen. The experiment was performed by three participants using three C. trachomatis isolates, of genotype E, F and B. Two surrogate urine contact methods were used to mimic contamination of a carer assisting with a child's urine collection. All EB dilutions and urine surrogate specimens were subjected to C. trachomatis assay and quantification in a real-time PCR-based diagnostic system. RESULTS: The amplimer crossing point (Cq) for EB dilutions was 10.0±1.6 less than for corresponding finger contacted urine specimens, which corresponds to ~10 µL of EB suspension transferred. This was largely independent of participant identity, C. trachomatis strain or EB dilution. Hand decontamination led to large reductions in EBs transferred, but transfer remained consistently detectable. Recent Cq data from C. trachomatis-positive clinical urine specimens were collated, and 20% clearly contained sufficient C. trachomatis to detectably contaminate another specimen by finger-mediated transfer, as in this experiment. CONCLUSIONS: This study directly demonstrated the potential for urine contaminated fingers to convert a C. trachomatis-negative urine specimen to C. trachomatis positive as a result of contact. Accordingly, procedures for urine specimen collection, particularly from children, need to be designed to prevent contamination.

Comprehensive global genome dynamics of Chlamydia trachomatis show ancient diversification followed by contemporary mixing and recent lineage expansion.

Chlamydia trachomatis is the world's most prevalent bacterial sexually transmitted infection and leading infectious cause of blindness, yet it is one of the least understood human pathogens, in part due to the difficulties of in vitro culturing and the lack of available tools for genetic manipulation. Genome sequencing has reinvigorated this field, shedding light on the contemporary history of this pathogen. Here, we analyze 563 full genomes, 455 of which are novel, to show that the history of the species comprises two phases, and conclude that the currently circulating lineages are the result of evolution in different genomic ecotypes. Temporal analysis indicates these lineages have recently expanded in the space of thousands of years, rather than the millions of years as previously thought, a finding that dramatically changes our understanding of this pathogen's history. Finally, at a time when almost every pathogen is becoming increasingly resistant to antimicrobials, we show that there is no evidence of circulating genomic resistance in C. trachomatis.

Primary health clinic toilet/bathroom surface swab sampling can indicate community profile of sexually transmitted infections.

BACKGROUND: The microbiome of built environment surfaces is impacted by the presence of humans. In this study, we tested the hypothesis that analysis of surface swabs from clinic toilet/bathroom yields results correlated with sexually transmitted infection (STI) notifications from corresponding human populations. We extended a previously reported study in which surfaces in toilet/bathroom facilities in primary health clinics in the Australian Northern Territory (NT) were swabbed then tested for nucleic acid from the STI agents Chlamydia trachomatis, Neisseria gonorrhoeae and Trichomonas vaginalis. This was in the context of assessing the potential for such nucleic acid to contaminate specimens collected in such facilities. STIs are notifiable in the NT, thus allowing comparison of swab and notification data. METHODS: An assumption in the design was that while absolute built environment loads of STI nucleic acids will be a function of patient traffic density and facility cleaning protocols, the relative loads of STI nucleic acids from different species will be largely unaffected by these processes. Another assumption was that the proportion of swabs testing positive for STIs provides a measure of surface contamination. Accordingly, "STI profiles" were calculated. These were the proportions that each of the three STIs of interest contributed to the summed STI positive swabs or notifications. Three comparisons were performed, using swab data from clinics in remote Indigenous communities, clinics in small-medium towns, and a single urban sexual health clinic. These data were compared with time and place-matched STI notifications. RESULTS: There were significant correlations between swab and notifications data for the both the remote Indigenous and regional data. For the remote Indigenous clinics the p values ranged from 0.041 to 0.0089, depending on data transformation and p value inference method. Further, the swab data appeared to strongly indicate known higher relative prevalence of gonorrhoeae in central Australia than in northern Australia. Similarly, the regional clinics yielded p values from 0.0088-0.0022. In contrast, swab and notifications data from the sexual health clinic were not correlated. DISCUSSION: Strong correlations between swab and notifications were observed. However, there was evidence for limitations of this approach. Despite the correlation observed with the regional clinics data, one clinic yielded zero positive swabs for C. trachomatis, although this STI constituted 25.1% of the corresponding notifications. This could be ascribed to stochastic effects. The lack of correlation observed for sexual health clinic data was also likely due to stochastic effects. It was concluded that toilet/bathroom surface swab sampling has considerable potential for public health surveillance. The approach may be applicable in situations other than primary health clinics, and for targets other than STIs.

Chlamydia trachomatis from Australian Aboriginal people with trachoma are polyphyletic composed of multiple distinctive lineages.

Chlamydia trachomatis causes sexually transmitted infections and the blinding disease trachoma. Current data on C. trachomatis phylogeny show that there is only a single trachoma-causing clade, which is distinct from the lineages causing urogenital tract (UGT) and lymphogranuloma venerum diseases. Here we report the whole-genome sequences of ocular C. trachomatis isolates obtained from young children with clinical signs of trachoma in a trachoma endemic region of northern Australia. The isolates form two lineages that fall outside the classical trachoma lineage, instead being placed within UGT clades of the C. trachomatis phylogenetic tree. The Australian trachoma isolates appear to be recombinants with UGT C. trachomatis genome backbones, in which loci that encode immunodominant surface proteins (ompA and pmpEFGH) have been replaced by those characteristic of classical ocular isolates. This suggests that ocular tropism and association with trachoma are functionally associated with some sequence variants of ompA and pmpEFGH.

Chlamydia trachomatis genotypes in a cross-sectional study of urogenital samples from remote Northern and Central Australia.

OBJECTIVES: The objective was to determine the frequency of trachoma genotypes of Chlamydia trachomatis-positive urogenital tract (UGT) specimens from remote areas of the Australian Northern Territory (NT). SETTING: The setting was analysis of remnants of C. trachomatis positive primarily UGT specimens obtained in the course of clinical practice. The specimens were obtained from two pathology service providers. PARTICIPANTS: From 3356 C. trachomatis specimens collected during May 2012-April 2013, 439 were selected for genotyping, with a focus on specimens from postpubescent patients, in remote Aboriginal communities where ocular trachoma is potentially present. PRIMARY AND SECONDARY OUTCOME MEASURES: The primary outcome measure was the proportion of successfully genotyped UGT specimens that were trachoma genotypes. The secondary outcome measures were the distribution of genotypes, and the frequencies of different classes of specimens able to be genotyped. RESULTS: Zero of 217 successfully genotyped UGT specimens yielded trachoma genotypes (95% CI for frequency=0-0.017). For UGT specimens, the genotypes were E (41%), F (22%), D (21%) and K (7%), with J, H and G and mixed genotypes each at 1-4%. Four of the five genotyped eye swabs yielded trachoma genotype Ba, and the other genotype J. Two hundred twenty-two specimens (50.6%) were successfully genotyped. Urine specimens were less likely to be typable than vaginal swabs (p<0.0001). CONCLUSIONS: Unlike in some other studies, in the remote NT, trachoma genotypes of C. trachomatis were not found circulating in UGT specimens from 2012 to 2013. Therefore, C. trachomatis genotypes in UGT specimens from young children can be informative as to whether the organism has been acquired through sexual contact. We suggest inclusion of C. trachomatis genotyping in guidelines examining the source of sexually transmitted infections in young children in areas where trachoma genotypes may continue to circulate, and continued surveillance of UGT C. trachomatis genotypes.

Whole genome sequencing to investigate a putative outbreak of the virulent community-associated methicillin-resistant Staphylococcus aureus ST93 clone in a remote Indigenous community.

We report two cases of severe pneumonia due to clone ST93 methicillin-resistant Staphylococcus aureus (MRSA) presenting from a remote Australian Indigenous community within a 2-week period, and the utilization of whole genome sequences to determine whether these were part of an outbreak. S. aureus was isolated from 12 of 92 nasal swabs collected from 25 community households (including the two index households); one isolate was ST93. Three of five skin lesion S. aureus isolates obtained at the community were ST93. Whole genome sequencing of the ST93 isolates from this study and a further 20 ST93 isolates from the same region suggested that recent transmission and progression to disease had not taken place. The proximity in time and space of the two severe pneumonia cases is probably a reflection of the high burden of disease due to ST93 MRSA in this population where skin infections and household crowding are common.

Novel staphylococcal species that form part of a Staphylococcus aureus-related complex: the non-pigmented Staphylococcus argenteus sp. nov. and the non-human primate-associated Staphylococcus schweitzeri sp. nov.

We define two novel species of the genus Staphylococcus that are phenotypically similar to and have near identical 16S rRNA gene sequences to Staphylococcus aureus. However, compared to S. aureus and each other, the two species, Staphylococcus argenteus sp. nov. (type strain MSHR1132(T) = DSM 28299(T) = SSI 89.005(T)) and Staphylococcus schweitzeri sp. nov. (type strain FSA084(T) = DSM 28300(T) = SSI 89.004(T)), demonstrate: 1) at a whole-genome level considerable phylogenetic distance, lack of admixture, average nucleotide identity <95 %, and inferred DNA-DNA hybridization <70 %; 2) different profiles as determined by MALDI-TOF MS; 3) a non-pigmented phenotype for S. argenteus sp. nov.; 4) S. schweitzeri sp. nov. is not detected by standard nucA PCR; 5) distinct peptidoglycan types compared to S. aureus; 6) a separate ecological niche for S. schweitzeri sp. nov.; and 7) a distinct clinical disease profile for S. argenteus sp. nov. compared to S. aureus.

Multisite Direct Determination of the Potential for Environmental Contamination of Urine Samples Used for Diagnosis of Sexually Transmitted Infections.

BACKGROUND: The detection of a sexually transmitted infection (STI) agent in a urine specimen from a young child is regarded as an indicator of sexual contact. False positives may conceivably arise from the transfer of environmental contaminants in clinic toilet or bathroom facilities into urine specimens. METHODS: The potential for contamination of urine specimens with environmental STI nucleic acid was tested empirically in the male and female toilets or bathrooms at 10 Northern Territory (Australia) clinics, on 7 separate occasions at each. At each of the 140 experiments, environmental contamination with Chlamydia trachomatis, Neisseria gonorrhoeae, and Trichomonas vaginalis nucleic acid contamination was determined by swabbing 10 locations, and urine collection was simulated 5 times, using a (1) synthetic urine surrogate and (2) a standardized finger contamination procedure. RESULTS: The most contaminated toilets and bathrooms were in remote Indigenous communities. No contamination was found in the Northern Territory Government Sexual Assault Referral Centre clinics, and intermediate levels of contamination were found in sexual health clinics and in clinics in regional urban centres. The frequency of surrogate urine sample contamination was low but non-zero. For example, 4 of 558 of the urine surrogate specimens from remote clinics were STI positive. CONCLUSIONS: This is by far the largest study addressing the potential environmental contamination of urine samples with STI agents. Positive STI tests arising from environmental contamination of urine specimens cannot be ruled out. The results emphasize that urine specimens from young children taken for STI testing should be obtained by trained staff in clean environments, and duplicate specimens should be obtained if possible.

Invasive Staphylococcus aureus Infections in Children in Tropical Northern Australia.

BACKGROUND: Despite a high burden of staphylococcal skin disease in children and high incidence of Staphylococcus aureus bacteremia in adult Indigenous populations in northern Australia, there are few studies describing incidence or clinical information of invasive S aureus (ISA) infections in children. METHODS: We conducted a retrospective review for all cases of S aureus bacteremia and sterile site infections, for children under 15 years, in northern Australia over a 4-year period (2007-2010). Cases were categorized as neonatal (<28 days) and pediatric (≥28 days). RESULTS: Forty-four cases (9 neonatal, 35 pediatric) were identified. The annual incidence of ISA was 27.9 cases per 100 000 population. Among pediatric cases, the annual incidence was significantly higher in the Indigenous (46.6) compared with the non-Indigenous (4.4) population (IRR: 10.6 [95% confidence interval, 3.8-41.4]). Pediatric infections were predominantly community-associated (86%). Clinical infection sites included osteoarticular (66%), pleuropulmonary (29%), and endocarditis (9%), and multifocal disease was common (20%). Eighty-three percent of pediatric cases presented with sepsis; 34% resulted in intensive care admission. Neonatal cases were all born prematurely; 89% were late-onset infections. Overall, 27% of infections were due to methicillin-resistant S aureus (MRSA). Compared with methicillin-sensitive S aureus (MSSA), there was no difference in severity or presentation in pediatric MRSA cases, but a higher proportion of MRSA cases were readmitted. CONCLUSIONS: The annual incidence of ISA infection in this study is among the highest described, largely due to a disproportionate burden in Indigenous children. Infections are frequently severe and infection with MRSA is common. Children presenting with suspected ISA in this region should be treated empirically for MRSA.