Whole genome sequence based capsular typing and antimicrobial resistance prediction of Group B streptococcal isolates from colonized pregnant women in Nigeria

Background Streptococcus agalactiae (Group B Streptococcus, GBS) is one of the major bacterial pathogens responsible for neonatal sepsis. Whole genome sequencing has, in recent years, emerged as a reliable tool for capsular typing and antimicrobial resistance prediction. This study characterized vaginal and rectal isolates of Group B Streptococcus obtained from pregnant women in Port Harcourt, Nigeria using a whole-genome sequence-based approach. Results Capsular types Ia, Ib, II, III, IV and V were detected among the 43 isolates sequenced. Twelve sequence types (STs) were identified, with ST19 (n = 9, 27.3 %) and ST486 (n = 5, 15.2 %) the most frequent among non-duplicated isolates. Of the alpha-like proteins (alp) identified, Alp1 was the most prevalent in 11 (33.3 %) isolates. Macrolide and lincosamide resistance determinants were present in 15 (45.5 %) isolates; ermB was detected in 1 (3 %), ermTR in 7 (21.2 %) isolates, lnu gene was detected in 6 (18.2 %) and mef was identified in 3 (9.1 %) isolates. Resistance of GBS to erythromycin and clindamycin (predicted from presence of erm or mef genes) was found to be 30.3 % and 24.2 %, respectively. All isolates were predicted resistant to tetracycline with only the tetM gene identified. Fluoroquinolone-resistance conferring substitutions in gyrA + parC were detected in 9 (27.3 %) isolates and chloramphenicol resistance was predicted from presence of aac6’-aph2 gene in 11 (33.3 %). Conclusions The data available from the whole genome sequencing of these isolates offers a small but insightful description of common serotypes and resistance features within colonizing GBS in Nigeria. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07929-z.


Introduction
Streptococcus agalactiae, also known as Group B Streptococcus (GBS) is increasingly being recognized in Africa as one of the major pathogens responsible for neonatal sepsis [1][2][3]. Vertical transmission from a rectovaginally colonized mother is the main predisposing factor for the development of early onset neonatal GBS sepsis [4,5]. Some countries, including the United States and Australia, have policies regarding intrapartum antibiotic prophylaxis (IAP) to prevent early-onset disease (EOD) with resulting benefits [6,7]. However, many low and middle income countries are yet to develop or implement such policies [8]. Whilst vaccine development is still underway [9,10], IAP remains one of the most effective tools for the prevention of EOD [7,11]. The use of IAP may inadvertently result in the increase of antimicrobial resistance (AMR) [8]. For this reason, determination of antimicrobial resistance patterns is important, especially in settings such as Nigeria where there is high incidence of indiscriminate antibiotic use [12]. Whole genome sequencing (WGS) has emerged as a reliable and cost-effective method for capsular typing and for accurate predictions of minimum inhibitory concentrations for antimicrobials in GBS [13][14][15].
This study analysed whole genome sequences of colonizing Group B Streptococcus isolates from vaginal and rectal swabs obtained from a previously described study in pregnant women in Port Harcourt, Nigeria [16]. Data from this study will add to existing evidence to guide vaccine efforts as well as contribute to local recommendations for intrapartum antibiotic prophylaxis.

Bacterial isolates
A total of 46 GBS isolates, obtained from separate vaginal and rectal samples of 185 pregnant women in their third trimester, collected retrospectively between March and July 2018 were used in this study. Samples were collected from women receiving antenatal care at the University of Port Harcourt Teaching Hospital in Port Harcourt Nigeria. Ethical approval was received from the hospital's ethics board and each participant signed a consent form prior to enrolment. Details of the study design and sample collection are described previously by Bob-Manuel et al. [16]. Briefly, the rectovaginal samples were incubated overnight in LIM RambaQUICK™ (CHROMagar, France) and subcultured onto CHROMagar™ StrepB (CHROMagar, France). Streptococcus agalactiae were presumptively identified as small, mauve-coloured colonies on CHROMagar StrepB. Isolates were stored in tryptic soy broth containing 5 % glycerol at -70°C. The isolates were shipped on Columbia agar with chocolated horse blood (Oxoid, UK) to the Streptococcus Laboratory, Centers for Disease Control and Prevention (CDC), for whole genome sequencing.

Whole genome sequencing and strain feature prediction
At CDC, GBS isolates were cultured on Trypticase soy agar supplemented with 5 % sheep blood prior to DNA extraction from single colonies as previously described [13]. Following shearing of samples using a Covaris M220 ultrasonicator (Covaris, Inc, MA, USA) genomic libraries were constructed using sparQ DNA Library Kit (Quantabio, MA USA) [13,14]. Short read sequences were generated on Illumina Miseq platform with MiSeq v2 500 cycle kit (Illumina Inc, USA) and assembled de novo using Velve-tOptimiser [17]. Serotypes, multilocus sequence typing (MLST), antibiotic resistance determinants, and predicted minimum inhibitory concentrations (MICs) were determined using a previously validated GBS bioinformatics pipeline [13]. Isolates flagged for phenotypic testing (those with lnu gene and ermTR + lnu) were tested by Etest (Biomerieux USA) using Mueller-Hinton agar with 5 % sheep blood (Becton Dickinson, USA), according to manufacturer's instructions. Interpretation of susceptible, intermediate and resistant from MIC WGS-prediction and Etest were determined using the standard Clinical Laboratory Standards Institute (CLSI) breakpoints for beta-hemolytic streptococci [18]. Additional strain features, including surface protein genes encoding the hypervirulent GBS adhesin (hvga), serine-rich repeat (srr) proteins, alpha protein family (alpha, Rib, Alp2/3, Alp1), and pilus proteins predictive of pilus islands (designated here as PI1, PI2A, and PI2B) were extracted from genomic data (https:// github.com/BenJamesMetcalf) [14]. Isolate identifiers, pipeline features, antimicrobial MIC predictions and interpretations, assembly metrics and NCBI genome accession numbers are listed in Supplementary Table 1. Contingency tables and a χ 2 test (or Fisher's exact test) were used to determine the significance of associations (at α = 0.05).

Results
The proportion of vaginal/rectal samples positive for GBS was 16.8 % (31/185 women) (sTable 1); 12 women had both vaginal and rectal swabs positive for GBS with 10 women having identical serotypes and therefore counted only once in further analyses. Whole genome sequencing revealed 3 of the isolates to be non-GBS. Six of ten known capsular serotypes (Ia, Ib, II, III, IV and V) were detected within this sampling, with serotype V the most frequent (n = 11/33, 33.3 %) ( Table 1). Multilocus sequence typing (MLST) identified a total of 12 STs among the 33 strains with ST19 (n = 9, 27.3 %) and ST486 (n = 5, 15.2 %) most frequently occurring. Eight (88.9 %) of the ST19 sequence strains were serotype V isolates while serotype III accounted for 1 (11.1 %) isolate. Only 1 serotype III isolate (3 %) had the highly virulent ST17 clonal type.

Discussion
This study reports the molecular characterization of colonizing strains of Group B Streptococcus in the vagina and rectum of pregnant women. Serotype V was the most frequent serotype in this population representing 33.3 % of the population.
Sequence types ST19, 486, 182, 24 and 2 accounted for 69.7 % of the total isolates. This diversity is different from other reports that have reported ST19, ST23 and ST1 more commonly [19]. The highly virulent ST17 responsible for most invasive disease, especially late onset disease in neonates [20], was rare in this study (3 %) population which may suggest a better prognosis for neonates who may become colonized at birth in this environment. While ST182 was reported as a frequent clonal type by Medugu et al. in another local study in Abuja, Nigeria [2], ST486 has not previously been documented as a common sequence type among GBS isolates in Nigeria. The pubMLST [21]  Alp1 (Epsilon) (33.3 %) was the most common of all the alpha-like proteins among the isolates in this study, similar to findings reported in Egypt (27 %) [23]. Other studies have documented Rib as the most frequently detected Alp gene in other geographic regions like China and Iran (37.5 and 53 % respectively) [24,25]. The serine-rich repeat protein Srr-1 was identified in almost all isolates, however only a single isolate had the Srr-2 known to be associated with greater binding affinity and higher morbidity [26]. The same isolate was positive for the hypervirulent GBS adhesin (HvgA), another surface protein demonstrated to increase adhesion and virulence particularly in the ST17 lineage [27].
Similar to other studies from Nigeria [2,28], all isolates were determined to be susceptible to beta-lactam antibiotics including 3rd generation cephalosporins. This is quite reassuring, suggesting that penicillin and ampicillin remain as viable options for intrapartum antibiotic prophylaxis in this region despite concerns for antimicrobial resistance. For women allergic to penicillin and at high risk of anaphylaxis, guidelines recommend clindamycin as an alternative option for prophylaxis [29]. In this study however, clindamycin and erythromycin resistance were relatively high at 24.2 and 30.3 %, respectively. Data from an earlier study in Nigeria reported erythromycin and clindamycin resistance rates of 6.5 % [28] suggesting increased resistance rates within the country. However, all isolates were susceptible to vancomycin in this study, so patients with a high risk for penicillin anaphylaxis who are colonized with a clindamycin resistant GBS strain, could receive vancomycin as the next best option [29].
Although other antibiotic classes such as fluoroquinolones, chloramphenicol, and tetracycline are not recommended for intrapartum prophylaxis, determining their susceptibility is important for tracking of resistance among GBS in general and as useful options for other invasive group B streptococcal infections in other age groups. All isolates were predicted as tetracycline resistant which is consistent with several other reports documenting high levels of resistance to this antibiotic [28,30,31]. Resistance to chloramphenicol and levofloxacin were also high at 25.6 %. Similarly, while these antibiotics may not be recommended for neonatal sepsis, they may be useful options for other invasive group B streptococcal infections in other age groups.
The only tetracycline gene present in this population was tetM. Of all the tetracycline resistance determinants, the tetM has the most spread geographically among GBS [28,30,31]. The presence of the double gyrA + parC mutants conferring full fluoroquinolone-resistance was high and were found predominantly in serotype V/ST-19 clone. The most frequently reported fluoroquinolone non-susceptible serotype/ST combinations are III/ST-19 and V/ST-19 [32,33] supporting the data from this study. Occurrence of either of these clonal complexes is likely a reflection of the circulating clones in those geographical areas in which they are found.
This study has demonstrated some salient characteristics of the prevalent colonizing strains of GBS in Nigeria, thus, contributing to the body of knowledge concerning GBS, as well as providing AMR data if IAP is considered in the future.