The draft genomes and investigation of serotype distribution, antimicrobial resistance of group B Streptococcus strains isolated from urine in Suzhou, China

Background The group B Streptococcus (GBS) is a human commensal bacterium, which is capable of causing several infectious diseases in infants, and people with chronic diseases. GBS has been the most common cause of infections in urinary tract of the elders, but relatively few studies reported the urine-isolated GBS and their antimicrobial susceptibilities. Hence, we decided to investigate GBS specially isolated from urine in Suzhou, China. Methods 27 GBS samples were isolated from urine in Suzhou, China. The PCR and agarose gel electrophoresis were used to identify the serotype distribution. Susceptibility tests were based on MIC test and Kirby–Bauer test. Genome were sequenced via Illumina Hiseq platform and assembled by SPAdes. Genomes of five isolates were sequenced and submitted to NCBI genome database. The sequencing files in fastq format were submitted to NCBI SRA database. Results Five serotypes were identified. The resistant rates measured for tetracycline, erythromycin, clindamycin and fluoroquinolones were 74.1, 63.0, 44.4 and 48.1%, respectively. 18.5% of the isolates were nonsusceptible to nitrofurantoin. The resistance to tetracycline was mainly associated with the gene tetM. The erythromycin resistance was mainly associated with the genes ermB and mefE. The genes ermB and lnuB were the prevalent genes in cMLSB type. No known nitrofurantoin resistance gene was found in nitrofurantoin-nonsusceptible GBS. Conclusions Five serotypes were identified in our study. High rates of GBS isolates were resistant to tetracycline, erythromycin, clindamycin and fluoroquinolones. The genes ermB and lnuB occupied high rates in cMLSB phenotype. Electronic supplementary material The online version of this article (10.1186/s12941-018-0280-y) contains supplementary material, which is available to authorized users.


Background
Group B Streptococcus (GBS, i.e. Streptococcus agalactiae) is a gram-positive microbe with bacterial morphology of sphere and oval observed via microscope [1]. Tilapia mossambica and cow are susceptible to GBS, leading to the reduced milk production [2,3]. GBS can be divided into ten serotypes, Ia, Ib, and II-IX, based on variant capsular polysaccharide. The serotypes Ia, Ib, III and V are the predominant pathogens, primarily liable for human infection [4]. GBS has traditionally been considered as a perinatal pathogen, in relation to certain infectious diseases in pregnant or postpartum women encomapassing endometritis, septicemia, chorioamnionitis and urinary tract infection [5]. Furthermore, GBS has been proved causing severe diseases in infants, such as pneumonia and meningitis [6]. Infectious pattern in infants were separated into two ways: mother-to-child vertically transmitted infection and birth canal transmission [7]. Screening for GBS in pregnant women via swabs taken from the vagina, of which the colonized GBS were mostly harmless, following by antibiotic treatment to prevent the transmission to the babies. GBS can also cause asymptomatic bacteriuria, cystitis, pyelonephritis, urethritis and urosepsis in urinary tract [8]. In a series of studies, bacteremia and urinary tract infections caused by GBS were the most common diagnoses among people over the age of 70, with 13 of 33 cases (39.4%) [9]. Previous studies indicated that all serotypes of GBS are susceptible to penicillin and have high rates of resistance to tetracycline, clindamycin and erythromycin [10], but little attention was paid to molecular mechanisms underlying those antimicrobial susceptibilities of GBS in urinary tract. A total of 27 GBS samples were isolated from urine in this study to investigate the serotype distribution and antibiotic susceptibilities in Suzhou, China.

Samples collecting and typing
All the GBS samples were isolated from urine between November 2015 and March 2017 in Suzhou Hospital of Traditional Chinese Medicine, Suzhou, China. The isolates were cultured on Blood agar medium at 37 °C with 5% CO 2 for 24 h. The DNA was extracted via the Ezup Column Bacteria Genomic DNA Purification Kit. The polymerase chain reaction (PCR) and agarose gel electrophoresis were applied (the primer sequences are in Additional file 1: Table S1, the PCR amplification system is in Additional file 1: Table S2 and the procedure of PCR is in Additional file 1: Table S3) to determine the serotype according to the expectation of electrophoresis result (Additional file 1: Table S4).

Antimicrobial susceptibility test and drug resistance gene PCR
Susceptibility tests to tetracycline, erythromycin, clindamycin, fluoroquinolones, nitrofurantoin, vancomycin, linezolid, quinupristin-dalfopristin and tigecycline were based on the MIC test. Susceptibility to ceftriaxone and penicillin was based on the Kirby-Bauer test recommended by the clinical and laboratory standards institute (CLSI) guidelines [11]. The clindamycin-induced resistance was detected based on D test [12]. The PCR of drug resistance genes was applied to find 12 resistance genes for tetracycline, erythromycin and clindamycin (Additional file 1: Table S5).

Genome sequencing and bioinformatics analysis
Genomes of four isolates (No.8,No.11,No.17 and No.24) showing multiple resistances and one isolate (No.20) indicating susceptive to all antibiotic were sequenced. The DNA was extracted by QIAGEN DNeasy ® Blood & Tissue and sequenced via Illumina Hiseq platform. The genomes were assembled by SPAdes [13] and the k-mer value was set to 127. The scaffolds < 1 kb in length were dropped. The genome sequences were submitted to Antibiotic Resistance Genes Database (ARDB, http://ardb. cbcb.umd.edu/) [14] to query drug resistance genes. The genes were predicted by Prokka (version 1.12).
Two protein sequences of nitrofurantoin-resistant genes, nsfA and nsfB, were received from The Comprehensive Antibiotic Resistance Database (CARD, https :// card.mcmas ter.ca/) [15] and then aligned to No.8 genome and No.11 genome sequences via TBLASTN (0.05 e-value) in order to find matching sequences.
The assembled draft genomes of the five isolates were submitted to NCBI genome and SRA database (SAMN08287475, SAMN08287476, SAMN08287477, SAMN08287478 and SAMN0828747). The sequencing files were submitted to NCBI SRA database (SRP128497).
Six complete serotype III genomes were download from NCBI (NZ_CP022537.

Resistance of tetracycline, erythromycin and clindamycin
All the 27 samples were susceptible to ceftriaxone, penicillin, vancomycin, linezolid, quinupristin-dalfopristin and tigecycline, while the highest percentage of 74.1% (20/27, Table 1 and Additional file 1: Table S6) samples were resistant to tetracycline. In previous studies, only the PCR primers of genes involved in resistance to tetracycline, erythromycin and clindamycin has been reported. In our PCR test of antibiotic resistance genes, tetM was the dominating resistance gene to tetracycline with 90.0% (18/20) occurrence. Other genes resistant to tetracycline were tetO, tetK and tetL, and the rates were 15.0% (3/20), 5.0% (1/20) and 5.0% (1/20) respectively.

Table 1 PCR result of antibiotic resistance genes
The order was sorted based on Additional file 1: Table S6. Tet genes are resistant to tetracycline; erm genes are resistant to clindamycin, erythromycin and streptogramin_b; mef genes are resistant to erythromycin; lnuB gene is resistant to clindamycin. Erythromycin is a kind of macrolides; clindamycin a kind of lincosamides TET tetracycline, CLI clindamycin, E erythromycin, FQNS fluoroquinolones. cMLS B constitutive resistance to macrolides, lincosamides, and streptograminB, iMLS B inducible resistance to macrolides, lincosamides, and streptograminB, MS resistance to macrolides and susceptibility to lincosamides

Serotype
Sample tetM tetO tetK tetL ermA ermB ermC ermM ermTR mefA mefE lnuB Resistance summary D test  Table S6). We sequenced the genome of No.11 isolate and found another known erythromycin resistance gene, ermT ( Table 3), suggesting that potentially there were other known resistance genes in the remaining six isolates.

Multiple resistance and detection for novel antibiotic resistance genes
Isolates of No. 8,No.11,No.17 and No.24 were nonsusceptible to at least three antibiotics. As the control, No.20 isolate was susceptible to all antibiotics. These five isolates were sequenced. The draft genome sizes were from 2.06 to 2.20 MB, and the scaffold N50 s were from 107,763 to 468,051 bp ( Table 2). As expected, none known resistance genes were found in PCR test and draft genome of No.20 isolate.
The No.8 isolate was resistant to tetracycline, erythromycin, clindamycin, and intermediate to nitrofurantoin. Genes tetO (tetracycline resistant), ermB (erythromycin and clindamycin resistant), mefE (erythromycin resistant) and lnuB (clindamycin resistant) were found in the PCR test. Genes tetO (tetracycline resistant), mefA (erythromycin resistant) and lnuB (clindamycin resistant) were found in the draft genome. No nitrofurantoin resistance gene was found (Tables 1, 3). No.11 isolate was resistant to tetracycline, erythromycin, and intermediate to nitrofurantoin. Tetracycline resistance gene tetM was found in PCR test and in the draft genome, but no erythromycin resistance gene was found in PCR test. An erythromycin and clindamycin resistance gene ermB was found in the draft genome. However, no nitrofurantoin resistance gene was found. About 48.1% (13/27, Additional file 1: Table S6) of the isolates were resistant to fluoroquinolones. One of them, No.17 isolate, was resistant to tetracycline, erythromycin and fluoroquinolones.   [16]. About two-thirds of the length of nsfB protein were aligned to all of the five draft genomes by TBLASTN search (E-value threshold 0.05, Additional file 1: Table S7), in which the genome area contained a 'Ribosomal large subunit pseudouridine synthase B' f in No.8 isolate, a 'glyoxal reductase' in No.11, a "putative NAD(P)H nitroreductase YfkO" in No.17,No.20 and No.24 isolates. So that the phenotype of nitrofurantoin intermediation was not caused by nsfB in our cases.
Based on the assumption that the mechanism of the nitrofurantoin intermediation in two isolates was similar, six complete serotype III genomes download from NCBI were included and built orthologous list with genomes of the nitrofurantoin-intermediate samples to search for candidate genes, given that the nitrofurantoin intermediation could not be explained by reported nitrofurantoinresistant genes. Considering that a detailed information for nitrofurantoin resistance of these complete serotype III genomes is limited, we supposed that at most one of these genomes has a nitrofurantoin-resistant phenotype according to the reference resistance occurrence. Seventeen orthologous contain genes in genomes of No.8, No.11 isolates only or mixed with at most one complete genome downloaded from NCBI. There were 18 genes of isolate No.8 and 17 genes of isolate No.11 in these seventeen orthologous (Additional file 1: Table S8). Most of those genes (16 of 35) were annotated as "hypothetical proteins" since their biological function has still not been clarified. We also identified several genes associated with transposition (6 of 35) and bacteriolysis (2 of 35), which appear to be unessential for the nitrofurantoin intermediation. Therefore, the nitrofurantoin-resistant related genes probably exist in those "hypothetical proteins". There were 27 genes distributed in the orthologous which contain genes only in isolate No.8 or mixed with at most one NCBI genome. 60 genes were distributed in the orthologous which contain genes only in isolate No.11 or mixed with at most one NCBI genome. Totally 148 genes of No.8 and 53 genes of No.11 were not recorded in this orthologous list. These genes cannot be found in another orthologous list which built between No.8 genome and No.11 genome either, suggesting that these 148 genes in No.8 and 53 genes in No.11 are unique genes of their own.
To further investigate the nitrofurantoin intermediation related genes in isolates No.8 and No.11 genome, the orthologous list of isolates No. 8,No.11,No.17,No.20 and No.24 was built, with the other three genomes as the nitrofurantoin-susceptible control. A high proportion of genes were annotated as "hypothetical proteins" (10 of 28, Additional file 1: Table S9), which was consistent with the above-mentioned observation. On the basis of our findings, it can be predicted that there is a higher possibility for the intersection of genes identified in both two orthologous lists to be the candidate genes, which were classified as hypothetical protein, helix-turnhelix domain protein and recombinase. With the range of candidates narrowed down, the four genes annotated as hypothetical proteins (Additional file 1: Table S9) were more likely to be the nitrofurantoin-resistant related genes, providing valuable reference for future experimental study.

Discussion
Five serotypes were detected in our research, encompassing Ia, Ib, III, V and VI. The serotype III, V and Ib were of the highest proportion. The serotype III GBS was the most common type detected in our isolates, which was in line with studies from other countries [17,18].
All of the GBS isolates studied in this work were susceptible to penicillin, indicating that penicillin is still a considerable choice for treatment of GBS infection. There was high rate of resistance to tetracycline (74.1%) in the test. However, this rate was lower than that reported in Beijing city (83.9%, isolated from vagina or rectum) [19], Taiwan (97.2%, isolated from urine, genital tract, genital tract, wound or pus, blood and respiratory tract) [20] and North American (89.2%, isolated from vagina or rectum) [17]. Compared to Shanghai, a city close to Suzhou, the rate of resistance to erythromycin in Suzhou (63.0%) is higher than reported in Shanghai (51.9%), however the rate of resistance to clindamycin in Suzhou (44.4%) is lower than that of Shanghai (61.5%) [21]. Note that, 80% of the nitrofurantoin-nonsusceptible GBSs in Suzhou were serotype III isolates, much higher than the 29.6% proportion of serotype III occurrence in the studied samples.
The tetM gene was the prevalent resistance gene (18/20) in tetracycline resistance isolates. The other detected tetracycline-nonsusceptible related genes were tetO, tetK and tetL. Besides, ermB (8/17) was the prevalent gene related to erythromycin resistance, which had been reported in Korea (isolated from vagina) [22]. It has been reported that the predominant genes in cMLS B were ermB in Spain [23,24] and ermB and ermB-mef(A/F) in Shanghai city [21]. Interestingly, except ermB, we found lnuB also showed a high proportion (7/10) in cMLS B in Suzhou (70.0%), five isolates have both ermB and lnuB genes, three isolates have ermB gene only and two isolates have lnuB gene only. Although this high rate of lnuB had not been reported in previous studies.
We suspected that new resistance genes probably exist in the genomes of isolates No.8 and No.11 since the phenotype of nitrofurantoin intermediation in our cases was not caused by nsfA and nsfB, nitrofurantoin-resistant genes that were recorded in CARD database. If the nitrofurantoin resistance mechanism is similar in the two samples, the 35 genes in the orthologous which contain genes both in the genomes of isolates No.8 and No.11 only or mixed with at most one NCBI genome would be potential candidates for the nitrofurantoin resistance. We also built another orthologous list with the other three genomes of the five isolates as the nitrofurantoin-susceptible control to integrate our finding. There is a greater likelihood for the four genes annotated as hypothetical proteins in the intersection of the two orthologous lists to be the nitrofurantoin-resistant related genes.

Conclusions
We isolated 27 GBS samples from urine in Suzhou, China. Five serotypes were identified (Ia, Ib, III, V and VI). The serotype III was the primary type. All of the samples were susceptible to ceftriaxone, penicillin, vancomycin, linezolid, quinupristin-dalfopristin and tigecycline. The resistance rates to tetracycline, erythromycin, clindamycin and fluoroquinolones were 74.1, 63.0, 44.4 and 48.1% respectively. The resistance to tetracycline was mainly associated with the gene tetM (18/20). The erythromycin resistance was mainly associated with the genes ermB (8/17) and mefE (7/17). The genes ermB (8/10) and lnuB (7/10) occupied high rates in cMLS B phenotype. We found five samples were nonsusceptible to nitrofurantoin, and four of them were serotype III, among which two genomes of isolates do not contain any reported nitrofurantoin resistance genes. Five isolates genomes were sequenced and submitted to NCBI genomes database (SAMN08287475, SAMN08287476, SAMN08287477, SAMN08287478 and SAMN0828747) and SRA database (SRP128497).
Authors' contributions GY collected the GBS samples, carried out the serotyping and drug testing experiment and sequenced genomes. DX carried out the bioinformatics analysis and drafted the manuscript. LY and ZL helped to draft the manuscript. ZGP raised useful suggestions. ZY conceived of the study, participated in its design and coordination, helped to draft the manuscript. All authors read and approved the final manuscript.