Skip to main content

Characterisation of Staphylococci species from neonatal blood cultures in low- and middle-income countries

Abstract

Background

In low- and middle-income countries (LMIC) Staphylococcus aureus is regarded as one of the leading bacterial causes of neonatal sepsis, however there is limited knowledge on the species diversity and antimicrobial resistance caused by Gram-positive bacteria (GPB).

Methods

We characterised GPB isolates from neonatal blood cultures from LMICs in Africa (Ethiopia, Nigeria, Rwanda, and South Africa) and South-Asia (Bangladesh and Pakistan) between 2015–2017. We determined minimum inhibitory concentrations and performed whole genome sequencing (WGS) on Staphylococci isolates recovered and clinical data collected related to the onset of sepsis and the outcome of the neonate up to 60 days of age.

Results

From the isolates recovered from blood cultures, Staphylococci species were most frequently identified. Out of 100 S. aureus isolates sequenced, 18 different sequence types (ST) were found which unveiled two small epidemiological clusters caused by methicillin resistant S. aureus (MRSA) in Pakistan (ST8) and South Africa (ST5), both with high mortality (n = 6/17). One-third of S. aureus was MRSA, with methicillin resistance also detected in Staphylococcus epidermidis, Staphylococcus haemolyticus and Mammaliicoccus sciuri. Through additional WGS analysis we report a cluster of M. sciuri in Pakistan identified between July-November 2017.

Conclusions

In total we identified 14 different GPB bacterial species, however Staphylococci was dominant. These findings highlight the need of a prospective genomic epidemiology study to comprehensively assess the true burden of GPB neonatal sepsis focusing specifically on mechanisms of resistance and virulence across species and in relation to neonatal outcome.

Peer Review reports

Background

Infections are a major cause of neonatal mortality [1, 2] with over 200,000 deaths reported each year [3]. Globally there is a large variation in the incidence of neonatal sepsis, although the burden is greatest in low- and middle-income countries (LMICs) [4, 5]. Despite Staphylococcus aureus being frequently reported as one of the top three bacterial causes of neonatal sepsis in LMICs along with GNB including Klebsiella pneumoniae and Escherichia coli [2, 6,7,8,9], surveillance data particularly using whole genome sequencing (WGS) is scarce. Understanding the true burden of coagulase negative Staphylococci (CoNS) in neonatal sepsis is far more complex. Despite multiple studies over the past decade suggesting that Staphylococcus epidermidis and Staphylococcus haemolyticus are important causes of nosocomial and bloodstream infections in neonates [10,11,12,13], there is often little capacity to implement routine surveillance, and these species are in many cases still considered contaminants, as they commonly reside on skin.

Studies have shown that neonates born at a very low birth weight (VLBW, < 1,500 g) may have a greater risk of contracting CoNS [11, 14]. Over the past two decades the number of specialised neonatal units in LMICs has fortunately increased, further increasing survival of VLBW infants. However, the likelihood of neonates exposed to nosocomial acquired CoNS is therefore higher, due to increased hospital stay for VLBW infants. Microbial investigation of CoNS isolated from blood cultures in countries with limited resources is often problematic due to a lack of simple and cost-effective methods of identification [9, 15]. Even if laboratories are equipped with manual or automated biochemical tests, these are laborious leading to subjective and inaccurate results; the tests add to the laboratory costs, and this may be difficult to justify given the often uncertain significance of the isolates. Although studies have reported cases of Staphylococcus epidermidis and Staphylococcus haemolyticus neonatal sepsis [9, 12, 13] and a PubMed search of “Staphylococci neonatal sepsis” produced over 300 hits (the majority of studies conducted in high-income countries) there is still a lack of data revealing the genomic diversity and resistome and virulome of staphylococcal species in LMICs.

BARNARDS is a network of 12 clinical sites distributed among four African (Ethiopia, Nigeria, Rwanda, South Africa) and three South-Asian (Bangladesh, India, Pakistan) countries aiming to assess the burden of antimicrobial resistance (AMR) in neonates. Herein we report data on sepsis caused by Gram-positive bacteria (GPB). We focus on the reporting of genomic diversity, genomic traits including antimicrobial resistance and virulence genes, and the phenotypic susceptibility profiles against 14 antibiotics for S. aureus isolates. Where possible, we performed genomics analysis on CoNS isolates to disclose AMR genotypes and whether mobile genetic elements were present. Finally, using WGS we report on an epidemiological cluster of neonatal sepsis caused by Mammaliicoccus sciuri, previously known as Staphylococcus sciuri in Islamabad, Pakistan between July-November 2017.

Methods

Study design

Ethical approval was obtained prior to the start of the study (Additional file 1: Table S1). We have used the term ‘neonate’ for all enrolments, this includes up to 60 days of age. Upon clinical signs of sepsis, a blood culture (using automated blood culture systems either BD BACTEC™ or BacT/ALERT® bioMérieux) was collected and processed using automated blood culture systems. Early onset sepsis (EOS) was defined as sepsis occurring < 72 h and late onset sepsis (LOS) > 72 h after birth. Neonates were followed up until 60 days. Clinical data collected included onset and outcome of sepsis. Due to high reporting of GPB sepsis and in particular S. aureus during a BARNARDS network meeting in 2017 (a program originally designed with an emphasis on carbapenem resistance in GNB neonatal sepsis) we retrospectively recovered GPB isolates for genomics characterisation in addition to GNB [16]. Unfortunately, at this point, 970/1266 GPB isolates, predominantly non-S. aureus were not stored or recoverable. Blood isolates were stored on charcoal swabs (Deltalab) for transport under UN3373 regulations to Cardiff University, UK.

Antimicrobial susceptibility testing

In Cardiff, GPB were plated onto Columbia Blood Agar supplemented with 5% sterile sheep/horse blood and incubated aerobically at 37 °C for 24–48 h. Isolates were identified using Microflex LT MALDI-TOF MS (Bruker Daltonik, GmbH, UK) with α-Cyano-4-hydroxycinnamic acid (HCCA) matrix (Sigma Aldrich) and stored in TS/72 beads (Technical Service Consultants, UK) at -80 °C. Minimum inhibitory concentrations (MICs) were determined by in-house agar dilution for 14 antibiotics (Additional file 1: Table S2) using Staphylococcus aureus ATCC 29,213, Escherichia coli ATCC 25,922 and Pseudomonas aeruginosa ATCC 2785 for quality control and interpreted according to the EUCAST v11 guidelines [17]. The MIC50 and MIC90 for each antibiotic was determined for S. aureus. Each bacterial isolate was tested once per antibiotic concentration to generate a single data point following validation using appropriate control strains listed above.

Whole genome sequencing and bioinformatics analysis

Assembly and annotation

WGS, assembly and annotation were performed as previously described [16], and individual accession numbers are listed in Additional file 1: Table S3. On average, isolates were paired-end sequenced to > 40 × coverage using Nextera XT v2 on an Illumina Miseq generating 300 bp reads. A summary of the bioinformatics programs and parameters used in this study can be found in Additional file 1: Fig. S1. Genome coverage was assessed using Qualimap (v2.2.1) [18]. ABRicate (v0.9.7) [19] (> 98% coverage and identity) was used to detect antimicrobial resistance genes (ARG) [20] and virulence factors (VF) [21]. Species identification was confirmed using BLAST v2.7.1 [22]. In silico MLST was determined for S. aureus, S. epidermidis and S. haemolyticus using mlst (v2.17.6) [23]. Novel alleles and profiles were submitted using PubMLST [24]. Staphylococcus aureus specific staphylococcal protein A (spa) gene typing, IS256 screening, and SCCmec cassette determination was performed using CGE [25] in combination with the following databases; spa typer (v1.0) [26], MGE (v1.0.2) [27] and SCCmecFinder (1.2) [28].

Comparative genomics

S. aureus (n = 251) from a European collection by Aanensen et al., 2016 [29] were incorporated into a species wide analysis using Prokka (v1.14.0) [30], Roary (v3.12.0) [31] and FastTree (v2.1.11) [32]. For all core genome SNP analyses performed (depicted in Additional file 1: Fig. S1), a representative reference of each dominant ST was selected from the genomes submitted to NCBI marked as ‘complete’, for variant calling using Snippy (v4.0.5) [33], Gubbins (v2.3.4) [34] snp-sites (v2.5.1)[35] and Raxml-ng (v0.9.0) [36] and IQ-tree (v2.0) where specified [37]. SNP pairwise distances were generated using snp-dists (v0.6) [38]. Phylogenetic trees were mid-rooted and annotated using iTOL (v5.7) [39]. In total n = 2,327 S. epidermidis and n = 310 S. haemolyticus genomes were downloaded from NCBI on 28th October 2020. Strains of the same ST (assessed using mlst v.2.17.6 [23]) of those found within this study were retained for comparative analysis. Metadata on whether these were reported as clinical, commensal, or other was recorded (Additional file 1: Table S4). SNP phylogeny workflows were performed as described above.

WGS analysis of Mammaliicoccus sciuri

An isolate belonging to the cluster was selected for complementary long read sequencing to produce a high-quality reference. gDNA was extracted and quantified as described [16], and SPRI beads (Mag-Bind TotalPure, Omega) were used to concentrate and purify at a 1:1 ratio with a final elution volume of 15µL to achieve an optimal range between 40-60 ng/µL. gDNA was quantified with the dsDNA BR assay kit using the qubit 4.0 fluorometer (Thermofisher). Genomic libraries were prepared using the Rapid Barcoding Kit (SQK-RBK004; Oxford Nanopore), sequenced on a R9.4 flow cell using a MinION (Oxford Nanopore) and basecalling was performed using Guppy within MinKnow. Long reads were demultiplexed using Porechop (v0.2.4) [40] and assembled with corresponding short reads generated from the Illumina MiSeq using Unicycler (v0.4.7) [41]. Mash dist (mash v2.2) [42] was performed on all M. sciuri genomes available in NCBI to determine genomic distance estimation to that of the cluster strain; the ten isolates with the greatest genomic similarity were incorporated into core SNP analysis as outlined previously. Metadata from NCBI was annotated within the phylogeny (Additional file 1: Table S5).

Statistical analysis

Statistical associations between clinical data (onset; EOS v LOS and outcome of sepsis; reported deceased v non reported deceased), phenotypic (MRSA) and genomic traits (ARG/ST/VF) were explored for S. aureus using Chi-Square (χ2) tests, comparing means T test and the Log Rank test SPSS (version 25.0.0.1). Kaplan–Meier survival curves were produced using the outcome follow up data (up to 60 days, or the last follow up record was used) in SPSS (version 25.0.0.1). The presence/absence of S. aureus VFs were analysed according to clinical onset and outcome of sepsis. Firstly, VFs associated with sepsis [43] were counted to provide a tailored VF score and a Kaplan–Meier survival curve was produced using the presence of Panton-Valentine Leukocidin (PVL), chosen as a frequently used clinical marker of S. aureus virulence. In addition, the frequency of all VFs detected > 98% identity and coverage from VFDB was counted against the clinical onset and outcome of sepsis to look for trends between EOS/LOS and alive/deceased. Statistical significance was taken at P value ≤ 0.05, and the false discovery rate was controlled using the Benjamini–Hochberg procedure [44]. A presence/absence matrix of all VFs detected per isolate (output from ABRicate vfdb) was generated using splist2presabs in the fuzzysim package in R 3.6.2.

Results

Characteristics of Gram-positive blood cultures reported during BARNARDS

In total, 1,266 GPB sepsis isolates were recorded from 1,239 neonates compared to 1,038 GNB sepsis isolates recorded from 1,007 neonates (Additional file 1: Table S6, Additional file 1: Fig. S2 [16]). From 296 GPB isolates available for identification, we identified 14 different species (Additional file 1: Fig. S3) and of these, 271 isolates were Staphylococci with 101 Staphylococcus aureus and 170 CoNS. Twenty-six isolates were identified as Staphylococci sciuri, but this species has since been reclassified to Mammaliicoccus sciuri. Overall, the onset of sepsis was 36% (n = 62) EOS, 54% (n = 71) LOS and for 10% (n = 17) cases the onset was undetermined (Additional file 1: Fig. S4). For S. aureus, the onset was 40% (n = 40) EOS, 47% (n = 47) LOS and undetermined for 13% (n = 14, Additional file 1: Fig. S4). Both S. epidermidis and S. haemolyticus was more common in LOS whereas M. sciuri was more common in EOS (Additional file 1: Fig. S4), however this data is indicative only.

This study specifically reports the antimicrobial susceptibility profiles of 144 isolates (n = 100 S. aureus and n = 44 CoNS) and WGS of 130 Staphylococci/Mammaliicocci sepsis isolates, (n = 100 S. aureus and n = 30 CoNS) which we were able to recover for microbiology characterisation after storage at − 80 °C.

Staphylococcus aureus and CoNS phylogeny

S. aureus for WGS characterisation included isolates predominantly from Pakistan (n = 42), Nigeria (n = 26), and South Africa (n = 16), along with nine from Bangladesh, two from Ethiopia and one from Rwanda, (n = 100 in total). In total, we detected 18 different sequence types (ST) (Fig. 1a) including the assignment of four novel STs; ST5426, ST5494, ST5945 and ST5496. We detected seven S. aureus clonal complex ST22 (CC22) from Bangladesh (n = 5), Pakistan (n = 1) and Rwanda (n = 1) and three ST5426 (within CC22) from Pakistan. Within BARNARDS, the commonest clonal complexes were CC5 (n = 36; 36%, consisting of ST6 [n = 18], ST5 [n = 17], ST650 [n = 1]) and CC8 (n = 13; 13%, all isolates belonged to ST8). Comparative genomic analyses including a previously known European collection, revealed ST5 isolates were distributed along 16 countries 45 (Fig. 1b). Interestingly, ST152 (n = 15) which is not assigned to a CC was the third most common ST identified. These isolates collectively represented 64% of the dataset and were isolated across the sampling period (Fig. 1a, Fig. 3a). All seventeen ST5 (Nigeria n = 6; South Africa, n = 11) were found carrying mecA (n = 12; 71%; Nigeria n = 3, South Africa n = 9; Fig. 2a) and the Nigerian and South African isolates grouped on different clades (Additional file 1: Fig. S5). Nine ST5 isolates collected from distinct South African neonates between July 2016 and July 2017 were within 10 pairwise SNPs (Supplementary Fig. 5), all were MRSA spa type t045, contained IS256 and virulence factors hlg, scn, spa and hly. Mortality following biological sepsis (MfBS) was reported for four out of these nine neonates. Additionally, 12 ST152 were isolated from African sites (Ethiopia n = 1, Nigeria n = 10 and South Africa n = 1, Fig. 1a, Additional file 1: Fig. S6) and were MSSA.

Fig. 1
figure 1

Core genome phylogenetic analysis of Staphylococcus aureus including a global contextual analysis. a Detailed core genome characterisation of 100 Staphylococcus aureus isolates (BARNARDS) using Roary (v3.12.0) and Fasttree (v2.1.11). Isolate labels are coloured according to clinical site of origin. The in-silico sequence type (ST) is shown outside of the isolate code (leaf). Presence of mecA, and whether the isolate was classified phenotypically as MRSA (as inferred from oxacillin MIC > 2 mg/l) is denoted by a filled triangle and/or circle respectively. Presence of IS256 is denoted by a filled rectangle. b Core genome characterisation of 351 Staphylococcus aureus isolates, incorporating a European collection [29] using Roary (v3.12.0) and Fasttree (v2.1.11). Coloured ranges in blue represent a S. aureus from the BARNARDS collection. Branch labels are coloured according to country of origin. Symbol represents source of isolate

Fig. 2
figure 2

Core genome phylogenetic analysis of coagulase negative Staphylococcus blood culture isolates. Core genome phylogeny of CoNS isolates displaying key genomic traits for comparison using Roary (v3.12.0) and Fasttree (v2.1.11). Isolate labels are coloured according to clinical site of origin. Clades are highlighted according to species. The in-silico sequence type (ST) is shown outside of the isolate code (leaf). Presence of mecA, and whether the isolate was classified phenotypically as MRSA (as inferred from oxacillin MIC > 2 mg/l) is denoted by a filled triangle and/or circle respectively. Presence of IS256 is denoted by a filled rectangle

ST6 (CC5) and ST8 (CC8) were found in South Asia (ST6 Pakistan n = 17/18, Bangladesh n = 1/18 and ST8 Pakistan n = 13, Fig. 1b, Supplementary Figs. 7–8). In the species-wide phylogeny, ST8 was found in 11 European countries, most notably Belgium and France (Fig. 1b). Within BARNARDS we saw three major ST8 clades, of which two were predominantly MRSA. (Fig. 1a, Supplementary Fig. 8). All ST8 spa type t064 genomes (n = 4/13) contained IS256 and genes coding for VFs associated with leukocyte targeting (Hlg) and toxin production (Hla). All 18 ST6 isolates were spa type t304, regardless of the clinical site of origin or whether they were MRSA (Additional file 1: Fig. S8). Twelve of the 17 ST6 S. aureus from Pakistan collected between August 2016-July 2017 grouped within 10 SNPs. Overall, we found a significant difference between neonatal survival and ST type (χ2, P 0.016) and a Kaplan survival curve analysis displaying patient survival for the four most dominant STs (ST5, ST6, ST8, ST152) shows ST5 and ST8 S. aureus being more likely to result in a fatal outcome (Fig. 3b). Of note, two neonates with a ST5 S. aureus sepsis from South Africa and two neonates with a ST8 S. aureus sepsis from Pakistan were also culture positive for at least one GNB isolate from the same blood sample and were excluded from the survival curve analysis.

Fig. 3
figure 3

Temporal frequency and survival curve data for Staphylococcus aureus blood culture isolates. a Stacked bar graph to show the temporal frequency of S. aureus isolates recovered from blood cultures during the BARNARDS sampling, per month. The bar graph is coloured according to the dominant STs, with all other STs being grouped as ‘Other’. b A Kaplan Meier survival plot comparing the four dominant ST groups (with all other STs detected grouped together as a single group called ‘other’) against the age at outcome for the neonate and up to 60 days. Findings are suggestive as the data presented is from a sample size (ST152 n = 15, ST other n = 37, ST5 n = 17, ST6 n = 18 and ST8 n = 13; overall comparison P 0.041)

A total of 30 CoNS were characterised by WGS (Fig. 2), S. epidermidis (n = 8), M. sciuri (n = 14) and S. haemolyticus (n = 8). The majority of S. epidermidis and S. haemolyticus were from Nigeria (NK n = 6, NN n = 7, NW n = 3), whereas n = 13/14 M. sciuri were from Asia, 12 being from one site in Pakistan (PP) and one from Bangladesh (BC). A single M. sciuri was recovered for WGS from Nigeria (NW). We found five STs from eight S. haemolyticus isolates including a novel assignment ST65 (Fig. 2). The most prominent ST (n = 3) found from all three sites in Nigeria was ST30, and all contained IS256. In n = 6/8 S. haemolyticus isolates IS256 was detected (Fig. 2). When screening public databases, we detected five ST30 strains, of which one was a sepsis case from India (Additional file 1: Fig. S9). ST30 isolates within this study were at least 150 SNPs different, and further distinct to the Indian sepsis isolate deposited within NCBI. We detected one ST1 and one ST3 S. haemolyticus. There are currently 44 ST1 and 49 ST3 S. haemolyticus submitted to NCBI and associated metadata reveals that n = 20/44 ST1 and n = 21/49 ST3 caused sepsis, largely within Europe, however there were sporadic entries from Russia, Brazil, India, and USA (Additional file 1: Figs. S10 and S11). ST56 was only submitted once into NCBI (up to October 2020), and this originated from a non-clinical sample from Germany in 2016. Here we report two ST56 isolates from blood cultures collected in Nigeria, both were collected from the same clinical site (NK), both were resistant to methicillin (Fig. 2) and had 707 pairwise SNPs.

We found that all 8 S. epidermidis isolates belonged to different STs, including the assignment of three novel STs (ST888-890), of which two (ST889 and ST890) were methicillin resistant (Fig. 2). IS256 was found in the ST2 S. epidermidis isolate from NN, Nigeria, however, was not found in any other isolates, although we present a limited dataset with available genomics data. ST490 was only detected once in the NCBI collection and was collected from a human sample. ST439 was detected twice, and an isolate from 2013 caused neonatal sepsis in the UK (ASM96671v1). Similarly, there were rare reports of ST210 within the NCBI collection (n = 5), and these were UK sepsis isolates from 2012 and 2017 (n = 3) and from human origin from the USA. ST2 was found in greater numbers in NCBI (n = 222), and the phylogeny separates into two main clades. The ST2 we report from Nigeria can be seen within the smaller clade alongside isolates from the environment, human commensal samples and those causing sepsis (Additional file 1: Fig. S12).

Antimicrobial resistance

In S. aureus, resistance rates were 53% (n = 49) to azithromycin, 32% (n = 32 to methicillin (oxacillin), 29% (n = 29) to tobramycin and 18% (n = 18) for gentamicin, (n = 100 isolates, except azithromycin n = 93; Table 1). Vancomycin and linezolid were the most active agents (0% resistance). The MIC90 of amikacin, tigecycline, minocycline, rifampicin, vancomycin and linezolid were all below the respective susceptible breakpoints (Table 1). Thirty-three were methicillin resistant S. aureus (MRSA), however the mecA gene was detected in 43% (n = 43). We also identified the reverse discrepancy, with three isolates negative for mecA (and mecC), but were classified MRSA by phenotypic methods. Typing of the staphylococcal cassette chromosome mec (SCCmec), revealed SCCmec_type_IVa(2B) was the most frequent (n = 30/43). All MRSA from Bangladesh carried SCCmec_type_IVa(2B) whereas all MRSA from South Africa carried SCCmec_type_IVa(1B). Although a 20% reduction in patient survival at 60 days was seen for neonates with MRSA compared to MSSA (Fig. 4a), this was not significant (Log Rank test, P 0.061), which may be due to the small sample size (n = 33). We found no significant difference between the onset of sepsis and MRSA (χ2, P 0.935 for phenotypic method and P 0.135 when mecA was used to denote MRSA).

Table 1 Antibiotic susceptibility testing for Gram-positive isolates recovered from blood cultures
Fig. 4
figure 4

Neonatal survival curve data for Staphylococcus aureus and pathogenicity markers: MRSA and PVL. Kaplan Meier survival plot comparing a MRSA v MSSA and time to neonatal outcome, and b The presence of virulence factor PVL and time to neonatal outcome censored at the last available follow up appointment

Although the MIC values indicate almost all S. aureus isolates were sensitive to both minocycline (Table 1; 93% [n = 93] sensitive) & tigecycline (100% [n = 100] sensitive), 25 isolates carried at least one tetracycline ARG (Fig. 2; n = 18 tet(K), n = 7 tet(M) and n = 3 tet(L). Additionally, aminoglycoside resistance genes were found in all S. aureus, and up to seven genes, including aac(3), aadD, sat4a and str, were found in a cluster of isolates from South Africa (Fig. 1). The insertion element IS256 is associated with the aminoglycoside resistance-mediating transposon Tn4001, and IS256 was found in 13/100 isolates. Macrolide ARG (ermA, ermC, InuA, mphC and msrA) were predominantly found in South Asian isolates, including ST1, ST6, ST8, ST22 and ST121 S. aureus. We only detected macrolide ARG in ST5 isolates from Africa, where msrA genes were detected in isolates from Nigeria, and ermA genes were detected in ST5 from South Africa (Fig. 1a).

For CoNS, isolates’ resistance rates were 50% (n = 22) to tobramycin and 41% (n = 18) to gentamicin (Table 2). Fourteen isolates carried blaZ, 10 of which also carried at least one aminoglycoside ARG (aadD and aph(3) genes; Fig. 2). We also detected IS256 in seven isolates (six S. haemolyticus and one S. epidermidis). Resistance to methicillin in CoNS (MR-CoNS) was seen in n = 11/18 isolates of S. epidermidis, n = 9/11 of S. haemolyticus and n = 2/13 of M. sciuri (Table 1, Fig. 2). mecA was detected in an additional four isolates, however all were sensitive to oxacillin. All M. sciuri isolates had the mecA1 gene, however this did not confer resistance to ß-lactam antibiotics, as reported previously 46. We detected two copies of mecA (mecA1_2 and mecA_9) in one isolate from Bangladesh, ultimately rendering this isolate resistant to oxacillin with an MIC > 8 mg/l (Fig. 2, Table 1). Furthermore, a single M. sciuri isolate from Pakistan was resistant to methicillin, tigecycline and vancomycin. Aside from these two isolates, the remaining 13 M. sciuri tested were sensitive to all other antibiotics. We also found a single vgaALC and catPC gene in two S. haemolyticus isolates, both of which also carried mecA, and aminoglycoside resistance genes (aadD and aph(3)-III).

Table 2 In silico Staphylococcus aureus virulence factors analysed for both onset and outcome of sepsis

Virulence factors

From 100 S. aureus isolates, 5494 genes putatively related to virulence (> 98% coverage and identity) were detected in total. All VF hits were summarised according to the gene presence/absence against sepsis onset (EOS v LOS) and outcome (alive v deceased) extending on analysis reported in [46]. For example, for the gene geh, we saw 32/40 (80%) presence in EOS as opposed to 24/46 (52.2%) in LOS. For the six genes that were significant before FDR correction (geh, cap8H, sak, cap8I, cap8K and scn), all had 13% or greater presence for EOS than for LOS (Additional file 1: Table S7). Similarly, lukD was present in 16/18 (88.9%) of S. aureus isolates recovered from deceased neonates, but only in 50/82 (61%) of neonates recorded as alive at the latest follow up appointment/60 days. For the eight genes that were significant before FDR correction, (lukD, sbi, aur, cap8B, scn, esaE, essC, hlgB), all bar one had at least 19% greater gene presence in the deceased than the survivors (Additional file 1: Table S7). The exception being scn where we saw 14/18 (77.8%) in the deceased with 77/82 (93.9%) in the survivors. This analysis was performed on the presence/absence of single genes, and further work is needed to investigate the role of multiple VF in onset/outcome of sepsis.

We also grouped VF related to bloodstream invasion/infection (n = 656) according to their biological functions as described by Powers et al., 2014 43 (Table 2, Additional file 1: Fig. S13). hlgA/B/C and lukD genes were found in 100% (n = 100) and 66% (n = 66) of S. aureus respectively, and PVL was found in 30% (n = 30). In neonates reported as deceased (n = 17), we found that the presence of PVL was found to significantly reduce the time to death (comparing means test, P 0.028), with the mean time to mortality at eight days compared to 31 days for PVL negative S. aureus. Similarly, the Kaplan Meier survival curve suggests the presence of PVL may result in a rapid MfBS, and this was more pronounced in the first 30 days of life, however, was not statistically significant (Log Rank test, P 0.251, Fig. 4b). In addition to virulence associated genes targeting leukocytes, we detected genes that contribute to inhibiting host complement systems (chp, and scn) and genes that contribute to S. aureus aggregation (Table 2, coa n = 8 isolates, cflB n = 4, clfA n = 16 and vWbp n = 17).

We found eight virulence associated genes from seven CoNS isolates. Two hits of the antiphagocytosis genes cap8G and cap8E within the same contig were detected in a S. haemolyticus isolate, whereas the remaining genes were detected within S. epidermidis. VF from the biofilm related ica operon (icaA) were detected in three S. epidermidis isolates, of which one contained IS256 and was MRSA (from Nigeria National Hospital, Abuja (NN)), whereas esxA, a gene associated with bacterial survival (type VII secretion system), were found in three S. epidermidis isolates (NK n = 2, NN n = 1). VFs present in VFDB were not detected within M. sciuri isolates.

M. sciuri cluster

Twenty-two M. sciuri were collected from Islamabad, Pakistan, between July-November 2017 (Fig. 5a). The blood culture date was not available for eight and we inferred approximate dates based on the consecutive enrolment numbers for the study (Fig. 5a). Where the onset of sepsis was known (n = 13), all were EOS suggesting the same strain may have spread within the hospital environment, perhaps related to the delivery area/delivery practices as, commonly reported, hospital contamination is more often linked to LOS [5, 47]. WGS data was available for 11 isolates which were all within four pairwise SNPs and were assigned to ST82, suggesting these represented the same/phylogenetically very close strains. The hybrid assembly of short and long read sequences produced a single contig for the bacterial chromosome at 2,932,999 bp (Additional file 1: Fig. S14). From genome estimation screening using mash (v2.2), the 10 genetically closest genomes to the Pakistani M. sciuri were > 1,400 SNPs away. The isolate that shares the same branch point in the collective phylogeny (Fig. 5b), was an M. sciuri blood culture isolate collected from Bangladesh during this study and was assigned to ST81. Similarly, the single M. sciuri isolate available for WGS characterisation from Nigeria was assigned to ST80. The MLST scheme for M. sciuri was first made available within pubmlst.org on the 1st of February 2021. Of the M. sciuri genomes that had the greatest Jaccard similarity index and were included into the phylogenetic analysis, three were clinical samples from Europe (Denmark and Czech Republic) and these were distinct throughout the tree suggesting different lineages causing human infection.

Fig. 5
figure 5

a Analysis of Mammaliicoccus sciuri recovered from blood cultures and a WGS contextual analysis. A timeline of M. sciuri neonatal sepsis in Pakistan, indicating which were available for whole genome characterisation. The blocks represent an individual case, and are colour coded according to the clinical outcome of sepsis. b A core genome SNP tree for all PP-BB isolates with WGS data available, performed using snippy-gubbins and Raxml (please refer to the text for details on methods) c Comparative phylogenetic tree of all S. sciuri with WGS data available in BARNARDS including the 10 genetically closest strains of M. sciuri from NCBI

Discussion

When a neonate is presenting with clinical signs of sepsis, empirical antibiotic treatment will be prescribed. Often, in LMICs hospitals there are not sufficient resources to perform necessary laboratory diagnostics, or a confirmatory blood culture. There is currently little to no data extending beyond an initial species identification (often via Gram staining and biochemical tests) and an antibiogram when Staphylococci/GPB are reported from neonatal blood cultures. We performed this retrospective study to gain an insight into the diversity of GPB isolated from blood cultures, and we highlight the need for future genomic and epidemiological studies to uncover the true burden of GPB neonatal sepsis in LMICs. A recent review by Garcia et al., 2022 describe a “Paucity of [care] bundles and elements aimed at HCAI detection (4% of 295 elements)” and “barriers to strengthening laboratory capacity”[48] when addressing matters of infection control in neonatal care in LMICs, emphasising the need for surveillance studies also mirrored elsewhere [49].

S. aureus is often the most frequently identified Gram-positive cause of neonatal sepsis in LMICs [2, 7, 50]. One of the most common STs found in this study was ST152, and although ST152 is not associated to a particular CC, scanning the list on the PubMLST database isolate collection [24], we noticed that of the 40 isolates with this ST, 12 were collected from blood cultures and all were from Africa (n = 9 The Gambia and n = 3 Mozambique). Herein, we found ST152 S. aureus neonatal sepsis in Ethiopia, Nigeria, and South Africa. Additionally, we detected a cluster of ST5 (CC5) in South Africa and previous published work from the same hospital suggests that ST5 is a successful and problematic lineage [51]. Many countries have reported CC5 MRSA from routine hospital molecular surveillance causing nosocomial infections further suggesting its global importance [51,52,53,54,55,56]. Frequently, however, nosocomial epidemiological clusters or outbreaks go undetected perhaps due to limited resources and/or protocols in place to detect outbreaks, outlining the importance of genomic surveillance studies to tackle nosocomial infection and provide insight into transmission dynamics of bacteria within neonatal units. Similar to all ST152 S. aureus being found only in Africa, during this study, all ST8 isolates were from neonates in South Asia. The numerous clades within the S. aureus ST8 group (as inferred from SNP analysis of the ST8 reads) suggest multiple introductions of ST8 causing neonatal sepsis in Pakistan. CC8 is well established and has commonly been reported from clinical infections worldwide, including South Asia [57, 58]. Importantly, we found that infection with ST8 and ST5 S. aureus had lower survival in the Kaplan Meir curve compared to other ST groups suggesting that certain lineages of S. aureus may be more pathogenic and problematic in neonatal sepsis in LMICs. Previous work on S. aureus bacteraemia also suggests that clonal differences may be a contributing indicator of infection outcome [59], and through preliminary data analysis we indicate that certain STs, ARGs, MGEs and VFs may be used as indicators of neonatal outcome.

Although S. aureus neonatal sepsis isolates in this study remained susceptible to vancomycin, similar to many previous reports [8, 60], we report an MRSA rate of 33% slightly higher than reported in a single site study in Uganda [61], and similar to a study in Nigeria [62], emphasising MRSA as an important source of infection and further highlighting the burden of AMR in neonatal sepsis in LMICs. Importantly, CoNS isolates often carry methicillin resistant genes [9, 13, 15], and through WGS and MIC determination, we identified methicillin resistance in one-third of isolates analysed, highest in S. epidermidis (n = 6/8) and S. haemolyticus (n = 7/8). Meric et al. [10] showed the importance of horizontal gene transfer events to increase the likelihood of divergent lineages causing infection due to the acquisition of genetic traits (virulence genes and ARG, e.g. SCCmecA) potentially conferring niche adaptations associated with pathogenicity. Kozitskaya et al.[63] also reports a significant association between the presence of IS256 in clinical S. epidermidis, perhaps suggesting that isolates with increased antimicrobial resistance are successfully colonising hospital environments. Although our WGS dataset for CoNS is too small for such comparisons, and found relatively few IS256 in S. epidermidis, we have shown concordance with previous studies and published genomics data, indicating that further surveillance studies are needed to characterise GPB isolates recovered from neonatal blood cultures.

The biggest question raised upon the isolation of CoNS from blood cultures is whether this confirms biological sepsis or represents contamination from transient colonisation of the skin. The ANISA study classified all CoNS as “definite contaminants”[64] and CoNS were likewise excluded from Okomo et al.’s.,[2] meta-analysis, many neonatal studies report CoNS as a leading cause of neonatal sepsis [65,66,67]. Our results must be interpreted with this in mind. A second blood culture is common practice following the identification of CoNS, however this remains challenging in neonates, as often a small volume of blood was collected (volume not recorded), and it was not possible for sites to consistently perform duplicate blood cultures following the local species identification of a CoNS. Speciation of CoNS can be difficult in laboratories where biochemical and conventional microbiology techniques are routinely employed [9], and long-term improvements in access to microbiology diagnostics in LMICs will be vital to understand the true burden of Staphylococci neonatal sepsis in these sites. Without this capacity, the cluster of M. sciuri described here could not be recognised and managed, as one clear example of the value of diagnostic capacity. M. sciuri have a broad range of habitats including animals, humans and the environment [45]. During a recent phylogenomic comparison of Staphylococci, several CoNS species were reclassified including Staphylococcus sciuri [68]. Although infection caused by M. sciuri has largely been reported in animals, there are case reports of M. sciuri causing peritonitis [69] and two distinct cases of catheter associated sepsis in Turkey [70] and Japan [71]. The clinical and microbiology team in Pakistan were routinely reporting and storing all GPB isolates from blood cultures, in contrast to other sites in South Asia which is likely to account for the steep difference in GPB neonatal sepsis, and the majority (n = 21, 81%) of M. sciuri isolates from Pakistan (Fig. 5). We did however identify five M. sciuri isolates from other clinical sites in Bangladesh, Ethiopia, and Nigeria. From these, two isolates (from Bangladesh and Nigeria) were shown to be sitting on separate phylogenetic branches through genomics analysis, (Fig. 5b). Furthermore, the identification of M. sciuri from 5/12 clinical sites highlights the potential for uncommon pathogens to be present in clusters or cause outbreaks.

There are limitations to this study. BARNARDS was established to determine the aetiology of GNB sepsis, therefore there was heterogeneity in the reporting and storage of GPB isolates among the clinical sites. Although many clinical sites that participated in BARNARDS identified GPB species from blood cultures as part of their routine standard operating procedures and guided by local practices, not all sites had the capacity and resources to identify bacteria to the species level. The retrospective recovery of Gram-positive isolates at local sites was greatest in Nigeria and Pakistan, therefore caution must be taken in any interpretation of results. There was discrepancy between the degree of reporting GPB neonatal sepsis within the clinical sites therefore we make no reference to prevalence or specific country analyses within this dataset. All statistical analyses between clinical data and microbiology findings are exploratory and hypothesis forming. Group B Streptococci (GBS), considered the dominant cause of sepsis in high-income countries [47, 72], was not found/reported, and this has also been observed in other studies [1, 5, 6, 8, 60, 73]. With limited microbiology resources in many LMICs however (GBS require additional media and growth requirements), septic neonates may be deceased prior to a blood culture and Vera et al., suggest GBS are perhaps underreported due to these infections presenting very early into life [47]. Additionally, due to the nature of our retrospective recovery, the loss of fastidious organisms at the local sites and/or transportation to the UK might have hindered the recovery of GBS and further studies are needed to determine the incidence of Staphylococci and GBS neonatal sepsis.

Conclusion

S. aureus is a leading cause of neonatal sepsis in LMICs, and we found one third of S. aureus to be MRSA. Future studies should assess the use of genomic characterisation to predict whether certain lineages of Staphylococci, and GPB species have greater pathogenic potential in neonatal sepsis. Unveiling a cluster of M. sciuri in Pakistan further demonstrates the importance of GPB in neonatal sepsis. There is an urgent need for microbiological and genomic surveillance studies with strict infection control practices to categorically evaluate the true burden of GPB in neonatal sepsis in LMICs, laying emphasis on whether associations exist between certain species, lineages and virulence factors and the clinical onset or outcome of sepsis.

Availability of data and materials

Sequences have been deposited in ENA with the project number: PRJEB40908, with individual accessions reported in the supplementary dataset (Additional file 1: Table S3). The M. sciuri genome used as a reference for SNP calling has been deposited in Genbank with the accession number CP071138.

Abbreviations

AMR:

Antimicrobial resistance

ARG:

Antibiotic resistance genes

ATCC:

American Type Culture Collection

BARNARDS:

Burden of Antibiotic Resistance in Neonates from Developing Societies

CC:

Clonal complex

CoNS:

Coagulase negative Staphylococci

EOS:

Early onset sepsis

EUCAST:

European Committee on Antimicrobial Susceptibility Testing

GNB:

Gram negative bacteria

GPB:

Gram positive bacteria

IS:

Insertion sequence

LMIC:

Low- and- middle income country

LOS:

Late onset sepsis

MfBs:

Mortality following biological sepsis

MGE:

Mobile genetic element

MIC:

Minimum inhibitory concentration

MLST:

Multi locus sequence typing

MRSA:

Methicillin resistant Staphylococcus aureus

PVL:

Panton-Valentine Leukocidin

SNP:

Single nucleotide polymorphism

ST:

Sequence type

VF:

Virulence factors

VFDB:

Virulence factor database

VLBW:

Very low birth weight

WGS:

Whole genome sequencing

References

  1. Hamer DH, Darmstadt GL, Carlin JB, Zaidi AKM, Bs MB. Etiology of bacteremia in young infants in six countries. Pediatr Infect Dis J. 2015;34(1):1–8.

    Article  Google Scholar 

  2. Okomo U, Akpalu ENK, le Doare K, Roca A, Cousens S, Jarde A, et al. Articles Aetiology of invasive bacterial infection and antimicrobial resistance in neonates in sub-Saharan Africa : a systematic review and meta-analysis in line with the STROBE-NI reporting guidelines. Lancet Infect Dis. 2019;3099(19):1–16. https://doi.org/10.1016/S1473-3099(19)30414-1.

    Article  Google Scholar 

  3. Naghavi M, Abajobir AA, Abbafati C, Abbas KM, Abd-Allah F, Abera SF, et al. Global, regional, and national age-sex specifc mortality for 264 causes of death, 1980–2016: a systematic analysis for the Global Burden of Disease Study 2016. Lancet. 2017;390(10100):1151–210.

    Article  Google Scholar 

  4. Fleischmann C, Reichert F, Cassini A, Horner R, Harder T, Markwart R, et al. Global incidence and mortality of neonatal sepsis : a systematic review and meta- analysis. Arch Dis Child. 2021;106:1–8.

    Article  Google Scholar 

  5. Popescu CR, Cavanagh MMM, Tembo B, Lufesi N, Goldfarb DM, Kissoon N, et al. Neonatal sepsis in low-income countries: epidemiology, diagnosis and prevention. Expert Rev Anti Infect Ther. 2020. https://doi.org/10.1080/14787210.2020.1732818.

    Article  PubMed  Google Scholar 

  6. Dramowski A, Madide A, Bekker A. Neonatal nosocomial bloodstream infections at a referral hospital in a middle-income country: burden, pathogens, antimicrobial resistance and mortality. Paediatr Int Child Health. 2015;35(3):265–72.

    PubMed  Article  Google Scholar 

  7. Medugu N, Iregbu K, Tam Pying I, Obaro S. Aetiology of neonatal sepsis in Nigeria, and relevance of Group b streptococcus: a systematic review. PLoS ONE. 2018. https://doi.org/10.1371/journal.pone.0200350.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Sankar MJ. Neonatal sepsis in South Asia: huge burden and spiralling antimicrobial resistance. BMJ. 2019;1(364): k5314.

    Google Scholar 

  9. Asante J, Amoako DG, Abia ALK, Somboro AM, Govinden U, Bester LA, et al. Review of clinically and epidemiologically relevant coagulase-negative staphylococci in Africa. Microb Drug Resist. 2020;26(8):951–70.

    CAS  PubMed  Article  Google Scholar 

  10. Méric G, Mageiros L, Pensar J, Laabei M, Yahara K, Pascoe B, et al. Disease-associated genotypes of the commensal skin bacterium Staphylococcus epidermidis. Nat Commun. 2018;9(1):1–11.

    Article  CAS  Google Scholar 

  11. Mary Healy C, Palazzi DL, Edwards MS, Campbell JR, Baker CJ. Features of invasive staphylococcal disease in neonates. Pediatrics. 2004;114(4):953–61.

    PubMed  Article  Google Scholar 

  12. Foka A, Chini V, Petinaki E, Kolonitsiou F, Anastassiou D, Dimitracopoulos G, et al. Clonality of slime-producing methicillin- resistant coagulase-negative staphylococci disseminated in the neonatal intensive care unit of a university hospital. Clin Microbiol Infect. 2006;12:1230–3.

    CAS  PubMed  Article  Google Scholar 

  13. Cavanagh JP, Hjerde E, Holden MTG, Kahlke T, Klingenberg C, Flægstad T, et al. Whole-genome sequencing reveals clonal expansion of multiresistant Staphylococcus haemolyticus in European hospitals. J Antimicrob Chemother. 2014;69(11):2920–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  14. Hornik CP, Fort P, Clark RH, Watt K, Benjamin DK, Smith PB, et al. Early and late onset sepsis in very-low-birth-weight infants from a large group of neonatal intensive care units. Early Hum Dev. 2012;88(SUPPL 2):1–13.

    Google Scholar 

  15. Sah S, Bordoloi P, Vijaya D, Amarnath SK, Sheela Devi C, Indumathi VA, et al. Simple and economical method for identification and speciation of Staphylococcus epidermidis and other coagulase negative Staphylococci and its validation by molecular methods. J Microbiol Methods. 2018;149(May):106–19.

    CAS  PubMed  Article  Google Scholar 

  16. Sands K, Carvalho MJ, Portal E, Thomson K, Dyer C, Akpulu C, et al. Characterization of antimicrobial-resistant Gram-negative bacteria that cause neonatal sepsis in seven low- and middle-income countries. Nat Microbiol. 2021;6(April):512–23.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  17. EUCAST. EUCAST V11: The European Committee on Antimicrobial Susceptibility Testing. Breakpoint tables for interpretation of MICs and zone diameters. Version 11.0, 2021. 2021.

  18. Okonechnikov K, Conesa A, García-Alcalde F. Qualimap 2: Advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics. 2016;32(2):292–4.

    CAS  PubMed  Google Scholar 

  19. Seemann T. ABRicate [Internet]. https://github.com/tseemann/abricate

  20. Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V, et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J Antimicrob Chemother. 2020;75(12):3491–500.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: Hierarchical and refined dataset for big data analysis—10 years on. Nucleic Acids Res. 2016;44(D1):D694–7.

    CAS  PubMed  Article  Google Scholar 

  22. Mcginnis S, Madden TL. BLAST : at the core of a powerful and diverse set of sequence analysis tools. Nucleic Acids Res. 2004;32:20–5.

    Article  CAS  Google Scholar 

  23. Seemann T. mlst, Github [Internet]. 2019. https://github.com/tseemann/mlst

  24. Jolley KA, Bray JE, Maiden MCJ. Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications [version 1; referees: 2 approved]. Wellcome Open Res. 2018;3:1–20.

    Article  CAS  Google Scholar 

  25. Thomsen MCF, Ahrenfeldt J, Cisneros JLB, Jurtz V, Larsen MV, Hasman H, et al. A bacterial analysis platform: an integrated system for analysing bacterial whole genome sequencing data for clinical diagnostics and surveillance. PLoS ONE. 2016;11(6):1–14.

    Article  CAS  Google Scholar 

  26. Bartels MD, Petersen A, Worning P, Nielsen JB, Larner-Svensson H, Johansen HK, et al. Comparing whole-genome sequencing with sanger sequencing for spa typing of methicillin-resistant Staphylococcus aureus. J Clin Microbiol. 2014;52(12):4305–8.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  27. Johansson MHK, Bortolaia V, Tansirichaiya S, Aarestrup FM, Roberts AP, Petersen TN. Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder. J Antimicrob Chemother. 2021;76(1):101–9.

    CAS  PubMed  Article  Google Scholar 

  28. International Working Group on the Classification of Staphylococcal Cassette Chromosome Elements (IWG-SCC). Classification of Staphylococcal Cassette Chromosome mec (SCC mec): Guidelines for Reporting Novel SCC mec Elements. J Antimicrob Chemother. 2009;53(12):4961–7.

    Article  CAS  Google Scholar 

  29. Aanensen DM, Feil EJ, Holden MTG, Dordel J, Yeats CA, Fedosejev A. Whole-genome sequencing for routine pathogen surveillance in public health: a population snapshot of invasive Staphylococcus aureus in Europe. MBio. 2016;7(3):1–15.

    Article  Google Scholar 

  30. Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30(14):2068–9.

    CAS  PubMed  Article  Google Scholar 

  31. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Sequence analysis Roary: rapid large-scale prokaryote pan genome analysis. Bioinformatics. 2015;31(July):3691–3.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  32. Price MN, Dehal PS, Arkin AP. FastTree 2—approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010. https://doi.org/10.1371/journal.pone.0009490.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Seemann T. snippy: fast bacterial variant calling from NGS reads. [snippy] 2015. https://github.com/tseemann/snippy.

  34. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43(3): e15.

    PubMed  Article  CAS  Google Scholar 

  35. Page AJ, Taylor B, Delaney AJ, Soares J, Seemann T, Keane JA, et al. SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. Microb Genom. 2016;2(4): e000056.

    PubMed  PubMed Central  Google Scholar 

  36. Kozlov AM, Darriba D, Flouri T, Morel B, Stamatakis A. RAxML-NG: A fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics. 2019;35(21):4453–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  37. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, Von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  38. Seemann T. snp-dists [Internet]. https://github.com/tseemann/snp-dists. Accessed Dec 2018.

  39. Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  40. Wick RR. Porechop v0.2.4 [Internet]. https://github.com/rrwick/Porechop. Accessed Nov 2018.

  41. Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13(6):1–22.

    Article  CAS  Google Scholar 

  42. Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH, Koren S, et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol. 2016. https://doi.org/10.1186/s13059-016-0997-x.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Powers ME, Wardenburg JB. Igniting the fire: Staphylococcus aureus virulence factors in the pathogenesis of sepsis. PLoS Pathog. 2014;10(2):10–3.

    Article  CAS  Google Scholar 

  44. Url S, Society RS. Controlling the false discovery rate: a practical and powerful approach to multiple testing Yoav Benjamini; Yosef Hochberg Journal of the Royal Statistical Society. Series B (Methodological), Vol. 57, No. 1. (1995), J R Statist Soc B 1995;57(1):289–300.

  45. Nemeghaire S, Argudín MA, Feßler AT, Hauschild T, Schwarz S, Butaye P. The ecological importance of the Staphylococcus sciuri species group as a reservoir for resistance and virulence genes. Vet Microbiol. 2014;171(3–4):342–56.

    PubMed  Article  Google Scholar 

  46. Thomson K, Dyer C, Liu F, Sands K, Carvalho MJ, Barrell M, et al. Evaluating the roles of antibiotic resistance, drug target attainment, bacterial pathogenicity, virulence and antibiotic access and affordability in affecting outcomes in neonatal sepsis: an international microbiology and drug evaluation prospective study. Lancet Infect Dis. 2021. https://doi.org/10.1016/S1473-3099(21)00050-5.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Zea-Vera A, Ochoa TJ. Challenges in the diagnosis and management of neonatal sepsis. J Trop Pediatr. 2015;61(1):1–13.

    PubMed  PubMed Central  Article  Google Scholar 

  48. Molina García A, Cross JH, Fitchett EJA, Kawaza K, Okomo U, Spotswood NE, et al. Infection prevention and care bundles addressing health care-associated infections in neonatal care in low-middle income countries: a scoping review. EClinicalMedicine. 2022;44:101259.

    PubMed  PubMed Central  Article  Google Scholar 

  49. Fitzgerald FC, Zingg W, Chimhini G, Chimhuya S, Wittmann S, Brotherton H, et al. the impact of interventions to prevent neonatal healthcare-associated infections in low- and middle-income countries: a systematic review. Pediatr Infect Dis J. 2022;41(3):S26-35.

    PubMed  PubMed Central  Article  Google Scholar 

  50. Ullah O, Khan A, Ambreen A, Ahmad I, Akhtar T, Grandapor AJ, et al. Antibiotic sensitivity pattern of bacterial isolates of neonatal septicemia in Peshawar, Pakistan. Arch Iran Med. 2016;19(12):866–9.

    PubMed  Google Scholar 

  51. Jansen van Rensburg MJ, Eliya Madikane V, Whitelaw A, Chachage M, Haffejee S, Gay Elisha B. The dominant methicillin-resistant Staphylococcus aureus clone from hospitals in Cape Town has an unusual genotype: ST612. Clin Microbiol Infect. 2011;17(5):785–92.

  52. Williamson DA, Roberts SA, Ritchie SR, Coombs GW, Fraser JD, Heffernan H. Clinical and molecular epidemiology of methicillin-resistant Staphylococcus aureus in New Zealand: rapid emergence of sequence type 5 (ST5)-SCCmec-IV as the dominant community-associated MRSA Clone. PLoS ONE. 2013;8(4):1–7.

    Article  CAS  Google Scholar 

  53. Challagundla L, Reyes J, Rafiqullah I, Sordelli DO, Echaniz-Aviles G, Velazquez-Meza ME, et al. Phylogenomic classification and the evolution of Clonal complex 5 methicillin-resistant Staphylococcus aureus in the Western Hemisphere. Front Microbiol. 2018;9:1–14.

    Article  Google Scholar 

  54. Teixeira MM, Araújo MC, Silva-Carvalho MC, Beltrame CO, Oliveira CCHB, Figueiredo AMS, et al. Emergence of clonal complex 5 (CC5) methicillin-resistant Staphylococcus aureus (MRSA) isolates susceptible to trimethoprim-sulfamethoxazole in a Brazilian hssospital. Braz J Med Biol Res. 2012;45(7):637–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  55. Kawaguchiya M, Urushibara N, Yamamoto D, Yamashita T, Shinagawa M, Watanabe N, et al. Characterization of PVL/ACME-positive methicillin-resistant staphylococcus aureus (Genotypes ST8-MRSA-IV and ST5-MRSA-II) isolated from a university hospital in Japan. Microb Drug Resist. 2013;19(1):48–56.

    CAS  PubMed  Article  Google Scholar 

  56. Kawaguchiya M, Urushibara N, Ghosh S, Kuwahara O, Morimoto S, Ito M, et al. Genetic diversity of emerging Panton-Valentine leukocidine/arginine catabolic mobile element (ACME)-positive ST8 SCCmec-IVa meticillinresistant Staphylococcus aureus (MRSA) strains and ACME-positive CC5 (ST5/ST764) MRSA strains in northern Japan. J Med Microbiol. 2013;62:1852–63.

    PubMed  Article  Google Scholar 

  57. Gorwitz RJ. Community-associated methicillin-resistant Staphylococcus aureus: epidemiology and update. Pediatric Infect Dis J. 2008;27(10):925–6.

    Article  Google Scholar 

  58. Strauß L, Stegger M, Akpaka PE, Alabi A, Breurec S, Coombs G, et al. Origin, evolution, and global transmission of community-acquired Staphylococcus aureus ST8. Proc Natl Acad Sci U S A. 2017;114(49):E10596–604.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  59. Recker M, Laabei M, Toleman MS, Reuter S, Saunderson RB, Blane B, et al. Clonal differences in Staphylococcus aureus bacteraemia-associated mortality. Nat Microbiol. 2017;2(10):1381–8.

    CAS  PubMed  Article  Google Scholar 

  60. Okomo U, Akpalu ENK, Le DK, Roca A, Cousens S, Jarde A, et al. Articles Aetiology of invasive bacterial infection and antimicrobial resistance in neonates in sub-Saharan Africa: a systematic review and meta-analysis in line with the STROBE-NI reporting guidelines. Lancet Infect Dis. 2019;3099(19):1–16.

    Google Scholar 

  61. Tumuhamye J, Sommerfelt H, Bwanga F, Ndeezi G, Mukunya D, Napyo A, et al. Neonatal sepsis at Mulago national referral hospital in Uganda: etiology, antimicrobial resistance, associated factors and case fatality risk. PLoS ONE. 2020. https://doi.org/10.1371/journal.pone.0237085.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Akindolire AE, Tongo O, Dada-Adegbola H, Akinyinka O. Etiology of early onset septicemia among neonates at the university college hospital, Ibadan, Nigeria. J Infect Dev Ctries. 2016;10(12):1338–44.

    CAS  PubMed  Article  Google Scholar 

  63. Kozitskaya S, Cho SH, Dietrich K, Marre R, Naber K, Ziebuhr W. The bacterial insertion sequence element is256 occurs preferentially in nosocomial Staphylococcus epidermidis isolates: association with biofilm formation and resistance to aminoglycosides. Infect Immun. 2004;72(2):1210–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  64. Hossain B, Weber MW, Hamer DH, Hibberd PL, Ahmed A, Marzan M, et al. Classification of blood culture isolates into contaminants and pathogens on the basis of clinical and laboratory data. Pediatr Infect Dis J. 2016;35(5):S52–4.

    PubMed  Article  Google Scholar 

  65. Ozkan H, Cetinkaya M, Koksal N, Celebi S, Hacimustafaoglu M. Culture-proven neonatal sepsis in preterm infants in a neonatal intensive care unit over a 7 year period: coagulase-negative Staphylococcus as the predominant pathogen. Pediatr Int. 2014;56(1):60–6.

    PubMed  Article  Google Scholar 

  66. Lamba M, Sharma R, Sharma D, Choudhary M, Maheshwari RK. Bacteriological spectrum and antimicrobial susceptibility pattern of neonatal septicaemia in a tertiary care hospital of North India. J Matern Fetal Neonatal Med. 2016;29(24):3993–8.

    CAS  PubMed  Google Scholar 

  67. Kohli-Kochhar R, Omuse G, Revathi G. A ten-year review of neonatal bloodstream infections in a tertiary private hospital in Kenya. J Infect Dev Ctries. 2011. https://doi.org/10.3855/jidc.1674.

    Article  PubMed  Google Scholar 

  68. Madhaiyan M, Wirth J, Saravanan V. Phylogenomic analyses of the Staphylococcaceae family suggest the reclassification of five species within the genus Staphylococcus as heterotypic synonyms, the promotion of five subspecies to novel species, the taxonomic reassignment of five Staphylococcu. Int J Syst Evol Microbiol. 2020. https://doi.org/10.1099/ijsem.0.004498.

    Article  PubMed  Google Scholar 

  69. Meservey A, Sullivan A, Wu C, Lantos PM. Staphylococcus sciuri peritonitis in a patient on peritoneal dialysis. Zoonoses Public Health. 2020;67(1):93–5.

    PubMed  Article  Google Scholar 

  70. Koçoğlu E, Karabay O. Catheter associated Staphylococcus sciuri sepsis in a patient with acute myeloid leukemia. Mikrobiyol Bul. 2006;40(4):397–400.

    PubMed  Google Scholar 

  71. Horri T, Suzuki Y, Kimura T, Kanno T, Maekawa M. Intravenous catheter-related septic shock caused by Staphylococcus sciuri and Escherichia vulneris. Scand J Infect Dis. 2001;33(12):930–2.

    Article  Google Scholar 

  72. Shabayek S, Ferrieri P, Spellerberg B. Group B streptococcal colonization in African countries: prevalence, capsular serotypes, and molecular sequence types. Pathogens. 2021;10:1606. https://doi.org/10.3390/pathogens10121606.

    Article  PubMed  PubMed Central  Google Scholar 

  73. Investigators of the Delhi Neonatal Infection Study (DeNIS) collaboration. Characterisation and antimicrobial resistance of sepsis pathogens in neonates born in tertiary care centres in Delhi, India: a cohort study. Lancet Glob Health. 2016;4(10):e752–60. https://doi.org/10.1016/S2214-109X(16)30148-6.

    Article  Google Scholar 

Download references

Acknowledgements

We would firstly like to acknowledge all the mothers’ and neonates that enrolled into the BARNARDS study, and all hospital and research staff at each of the participating clinical sites. We thank and acknowledge the hard work of the team at Liofilchem®, Roseto, Italy, specifically F. Brocco and M. Collett. We would like to acknowledge Specialist Antimicrobial Chemotherapy Unit (SACU) at Public Health Wales for support with bacterial isolate identification using MALDI-TOF MS. We thank Maksuda Islam for providing access to the clinical sites and epidemiology data in Bangladesh. We would like to acknowledge Rubina Kamran the Microbiologist from PIMS who sadly passed away in 2018. We would like to thank Julian Parkhill for advice on SNP and phylogenetic analysis. We would like to thank the team at Bill and Melinda Gates Foundation namely Padmini Srikantiah, Rasa Izadnegahdar, Keith Klugman and Sara Vernam. We would also like to acknowledge Wales Gene Park and ARCCA for their continued bioinformatics support. Bioinformatics analysis was undertaken using the supercomputing facilities at Cardiff University operated by Advanced Research Computing at Cardiff (ARCCA) on behalf of the Cardiff Supercomputing Facility and the HPC Wales and Supercomputing Wales (SCW) projects. We acknowledge the support of the latter, which is part-funded by the European Regional Development Fund (ERDF) via the Welsh Government. The Centre for Trials Research is funded by Health and Care Research Wales and Cancer Research UK. We thank the team of curators for the databases hosted on PubMLST https://pubmlst.org/databases/.

Funding

BARNARDS was funded by two awards $4.28M (OPP1119772) and $849K (OP1191522) from the Bill and Melinda Gates Foundation at 500 Fifth Avenue North-Seattle, WA 98109. The funders had no role in the design of the study, data collection, analysis, interpretation of data or manuscript preparation.

Author information

Authors and Affiliations

Authors

Consortia

Contributions

KS guided the analysis, wrote the manuscript and produced the Figures. KS performed the bioinformatics analysis supported by RA, MJC, JP and BS. KS, EP, and JM performed whole genome sequencing experiments. KT performed the minimum inhibitory concentration (MIC) experiments. KT and MJC analysed the MIC dataset. WJW and MJC performed the statistical analysis. KS, MJC, EP, KT, CD, CA, AF, TH, JM, and MN performed microbiology culture and sample processing in Cardiff University. CD and RM provided the managed epidemiological dataset. GC, DB, SS, SB, RKN, SN, KI, FM, SU, RZ, HS, AM, JBM, AR, LG, SM, ANHB and AW collected and processed blood cultures and collected clinical data at the clinical sites. BS and JM guided the analysis and writing of the manuscript. TRW, KS, and MJC, RM, GC, SB, KI, RZ, JBM and SM designed the BARNARDS study. All authors reviewed and commented on the manuscript.

Corresponding author

Correspondence to Kirsty Sands.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained at each of the seven participating countries (Additional file 1: Table S1). Bangladesh: Ethical Review Committee, Bangladesh Institute of Child Health (BICH-ERC-4/3/2015), Ethiopia: Boston Children's Hospital (IRB-P00023058), India: Institutional Ethics Committee, National Institute of Cholera and Enteric Diseases and Institute of Post Graduate Medical Education and Research, IPGME&R Research Oversight Committee (A-I/2016-IEC and Inst/IEC/2016/508), Nigeria: Kano State Hospitals Management Board (8/10/1437AH), Health Research Ethics Committee (HREC) and National Hospital, Abuja (NHA/EC/017/2015), Pakistan: Shaheed Zulfiqar Ali Bhutto Medical University, Pakistan Institute of Medical Sciences (PIMS) Islamabad (Ref No NA, signed letter from Prof. Tabish Hazir). Rwanda: Republic of Rwanda, National Ethics Committee (No342/RNEC/2015), South Africa: Stellenbosch University and Tygerberg Hospital, Research projects, Western Cape Government (N15/07/063). In local languages, research nurses provided mothers with study information and collected consent for mother and/or neonatal enrolment. Informed consent was obtained in writing unless this was not possible (due to literacy barriers), and oral consent was collected from the mothers by trained researchers. Verbal consent was approved for use as informed consent if written consent was not possible by each committee and this was documented by the participant signing/marking the consent form.

Consent for publication

Not applicable.

Competing interests

The authors have no competing interests to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

Additional tables and figures.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sands, K., Carvalho, M.J., Spiller, O.B. et al. Characterisation of Staphylococci species from neonatal blood cultures in low- and middle-income countries. BMC Infect Dis 22, 593 (2022). https://doi.org/10.1186/s12879-022-07541-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-022-07541-w

Keywords

  • LMIC
  • Staphylococci
  • Neonatal sepsis
  • Genomics
  • Mammaliicocci
  • Mortality
  • Early onset
  • Late onset
  • Epidemiology