Next Article in Journal
How to Inhibit Nuclear Factor-Kappa B Signaling: Lessons from Poxviruses
Next Article in Special Issue
Changes in Group B Streptococcus Colonization among Pregnant Women before and after the Onset of the COVID-19 Pandemic in Brazil
Previous Article in Journal
Tetraplex Fluorescent Microbead-Based Immunoassay for the Serodiagnosis of Newcastle Disease Virus and Avian Influenza Viruses in Poultry Sera
Previous Article in Special Issue
Group B Streptococcus among Pregnant Women and Neonates in Saudi Arabia: A Systemic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Epidemiology of Group B Streptococci in Lithuania Identifies Multi-Drug Resistant Clones and Sporadic ST1 Serotypes Ia and Ib

by
Jonah Rodgus
1,2,3,†,
Ruta Prakapaite
4,5,†,‡,
Panagiotis Mitsidis
1,2,
Ramune Grigaleviciute
6,‡,
Rita Planciuniene
7,‡,
Povilas Kavaliauskas
4,8,9,‡ and
Elita Jauneikaite
1,3,*,‡
1
Department of Infectious Disease Epidemiology, School of Public Health, Imperial College London, London W2 1PG, UK
2
MRC Centre for Molecular Bacteriology and Infection, Department of Infectious Disease, Imperial College London, London SW7 2AZ, UK
3
NIHR Health Protection Research Unit in Healthcare Associated Infections and Antimicrobial Resistance, Department of Infectious Disease, Imperial College London, Hammersmith Hospital, London W12 0NN, UK
4
Institute of Infectious Diseases and Pathogenic Microbiology, 59116 Prienai, Lithuania
5
Biozentrum, University of Basel, 4056 Basel, Switzerland
6
Biological Research Centre, Lithuanian University of Health Sciences, 47181 Kaunas, Lithuania
7
Institute of Microbiology and Virology, Lithuanian University of Health Sciences, 47181 Kaunas, Lithuania
8
Division of Infectious Diseases, Department of Medicine, Weill Cornell Medicine of Cornell University, New York, NY 10065, USA
9
Department of Microbiology and Immunology, University of Maryland School of Medicine, Baltimore, MD 21201, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
These authors are members of Surveillance and Molecular Epidemiology of Maternal Infections Consortium.
Pathogens 2022, 11(9), 1060; https://doi.org/10.3390/pathogens11091060
Submission received: 15 August 2022 / Revised: 11 September 2022 / Accepted: 13 September 2022 / Published: 17 September 2022
(This article belongs to the Special Issue Group B-Streptococcus (GBS))

Abstract

:
Streptococcus agalactiae (Group B Streptococcus, GBS) is a leading cause of neonatal infections. Yet, detailed assessment of the genotypic and phenotypic factors associated with GBS carriage, mother-to-baby transmission, and GBS infection in neonates and adults is lacking. Understanding the distribution of GBS genotypes, including the predominance of different serotypes, antimicrobial resistance (AMR) genes, and virulence factors, is likely to help to prevent GBS diseases, as well as inform estimates of the efficacy of future GBS vaccines. To this end, we set out to characterise GBS isolates collected from pregnant and non-pregnant women in Kaunas region in Lithuania. Whole genome sequences of 42 GBS isolates were analysed to determine multi-locus sequence typing (MLST), the presence of acquired AMR and surface protein genes, and the phylogenetic relatedness of isolates. We identified serotypes Ia (42.9%, 18/42), III (33.3%, 14/42), V (21.4%, 9/42), and a single isolate of serotype Ib. Genomic analyses revealed high diversity among the isolates, with 18 sequence types (STs) identified, including three novel STs. 85.7% (36/42) of isolates carried at least one AMR gene: tetM or tetO (35/42), ermB or lsaC (8/42) and ant6-Ia and aph3-III (2/42). This study represents the first genomic analysis of GBS isolated from women in Lithuania and contributes to an improved understanding of the global spread of GBS genotypes and phenotypes, laying the foundations for future GBS surveillance in Lithuania.

1. Introduction

Group B Streptococcus (GBS) is a β-haemolytic gram-positive bacterium associated with the colonisation of mucous membranes in the human body. A commensal in the gastrointestinal and lower rectogenital tracts of up to 36% of pregnant women in Europe [1,2], GBS can be transmitted from mother to neonate, with invasive GBS infections a leading cause of neonatal pneumonia, septicemia, and meningitis worldwide [3]. At present, the administration of intrapartum antimicrobials represents the only clinical intervention available to prevent GBS diseases in neonates [4]. However, antimicrobial resistance (AMR) and multi-drug resistance (MDR) conveyed by the presence of increasingly common genomic determinants (both genes and mutations) undermine the future viability of this intervention [5]. The presence of tetracycline resistance genes (tetM and tetO), for instance, is almost ubiquitous in GBS, with over 85% of isolates carrying either tetM or tetO [6,7]. Moreover, recent studies on GBS have reported an increase in macrolide, fluoroquinolone, and aminoglycoside resistance, driven by the genes ermB, lnuC, mefA, and msrD [5], mutations in gyrA and parC, and the genes aac(6’)-aph(2”), ant6-Ia, and aph3-III [5,6], respectively. Mutations in penicillin-binding-proteins (PBP), namely PBP1a, PBP2a, PBP2b, and PBP2x, that reduce the susceptibility of GBS to penicillin have been reported [8,9,10,11]. Mutations in pbp1a and pbp2x, specifically, have also been associated with reduced susceptibility to beta-lactams [12,13]. Given the increasing presence of these AMR determinants, there is a need for routine surveillance of GBS isolates and to better characterise circulating antimicrobial resistant GBS genotypes and phenotypes.
Monitoring the incidence of GBS serotypes, AMR determinants, and disease-causing GBS populations is likely to inform estimates of GBS vaccine efficacy. Of the ten known GBS serotypes (Ia, Ib, II-IX), just five serotypes, Ia, Ib, II, III, and V, account for 93-99% of cases of neonatal and adult infections globally [14]. As such, the main GBS vaccines in development include a trivalent vaccine covering serotypes Ia, Ib, and III [15]; a hexavalent vaccine covering serotypes Ia, Ib, and II-V [16,17], and a multivalent adjuvanted protein vaccine (NCT03807245). The latter of these vaccines exploits two fusion proteins constituting different combinations of two of the four most common Alp-family proteins, Alpha-C (also known as bca) and Rib, or Alp1 and Alp2/3 (https://www.minervax.com/product-pipeline/ accessed on 9 September 2022). Phase I clinical trials have shown high levels of antibody response in study participants, though further work is needed to investigate long term immunity and the potential interaction with pre-existing immunity to GBS [18,19].
Despite recognition of the need for a GBS vaccine more widely, there exists a distinct lack of detailed studies focused on the epidemiology of GBS carriage and infections in pregnant women, their neonates, and non-pregnant adults in Lithuania. Previous studies of GBS carriage in Lithuania have reported maternal and neonatal GBS carriage rates between 2006–2007 of 15.3% and 6.4% [20], respectively, with increased maternal GBS colonisation rates reported in later years (18.6% in 2012 and 22.9% in 2014 [21]), and mortality in neonates due to GBS infection having increased from 5.2% in 1999 to 18.1% in 2012 [22]. Here, we report the serotypes and genotypes of clones found among 42 clinical GBS isolates from pregnant and non-pregnant women in Lithuania, with findings laying the foundations for future surveillance studies in this country.

2. Results

2.1. Overview and Genomic Characterisation of Lithuanian GBS Isolates

Of the 42 GBS isolates, 83% (35/42) were retrieved from vaginal swab samples, 14% (6/42) from urine samples, and one from a blood sample of an individual with confirmed bacteraemia (Table S1). Most isolates (85.7%, 36/42) were from women reported to be pregnant at the time of sample collection, with 32y the median age for this cohort (Table S1). The isolates comprised four serotypes, Ia, Ib, III, and V. Serotype Ia was most prevalent, representing 42.9% (18/42) of isolates (including the single blood isolate, Table S1). Serotypes III, V, and Ib represented 33.3% (14/42), 21.4% (9/42), and 2.4% (1/42) of isolates, respectively. We identified 18 distinct sequence types (STs) using multi-locus sequence typing (MLST), including three novel STs (ST1865-1867). The most common STs were ST17, ST1, ST88, and ST144 (Table S1). STs clustered into five major clonal complexes (CCs), CC23 (n = 14), CC17 (n = 9), CC1 (n = 7), CC19 (n = 5), and CC452 (n = 4), each of which broadly corresponded to a major clade from the phylogenetic tree constructed for all isolates (Figure 1).
We identified one serotype Ib ST1 isolate and two serotype Ia isolates from CC1 - one ST1 and one ST1387 (a single locus variant of ST1). As such, further investigations were carried out to identify any potential capsular switching events amongst the CC1 isolates (constituting mainly ST1 genotypes). Compared to other publicly available serotype Ia ST1 isolates (n = 8), our CC1 isolates of serotypes Ia (n = 2), V (n = 4), and Ib (n = 1) appeared to have participated in a capsular switching event. Indeed, all serotype Ia isolates, except two (LTS018 and SRR7283561), were genetically similar to other serotype V isolates, as well as the serotype Ib ST1 isolate, with less than 90SNPs (recombination-free core SNPs) present (Figure S1a). Moreover, recombination borders were very similar for the Ib isolate (LTS066) and two Ia isolates (LTS045 and LTS018), strongly indicative of a capsular switching event between serotypes V and Ib, and serotypes V and Ia (Figure S1b).

2.2. Antimicrobial Resistance Profiles and Antimicrobial Resistance Associated Mobile Genetic Elements

Overall, 64% (27/42) of isolates carried only one tetracycline resistance gene, either tetM (n = 25) or tetO (n = 2). Two isolates (LTS011 and LTS031) carried one tetracycline resistance gene (tetM or tetO), the erythromycin resistance gene, ermB, as well as the aminoglycoside resistance genes, ant6-Ia and aph3-III (Figure 1). Nine isolates carried the macrolide resistance genes ermA, ermB, and lsaC (Figure 1). Four of the five isolates carrying ermB had confirmed resistance to erythromycin and clindamycin (Figure 1, Table S1). Six serotype Ia ST88 isolates did not carry any acquired AMR genes. All isolates were susceptible to penicillin and, as expected, did not carry any mutations in the genes pbp1A or pbp2X, previously associated with resistance to penicillin (Table S1). Additionally, we identified no mutations attributed to fluoroquinolone resistance in the quinolone resistance-determining regions, gyrA and parC [23].
Acquired AMR genes are often carried by mobile genetic elements (MGEs). Hence, we set out to investigate the presence of MGEs in our GBS isolates. The most common transposon, Tn916, known to carry tetM, was detected in 14 isolates. Furthermore, Tn3872 (Figure S2), carrying tetM and ermB, was detected in two ST1 isolates (LTS039 and LTS066). In one of the two multi-drug resistant isolates (LTS011, serotype III, ST17), we detected a putative integrative conjugative element (ICE, named here as LTS_ICESag37), downstream of the recombination hotspot gene, rplL, with this ICE carrying the acquired AMR genes, tetO, ant6-Ia, aph3-III, and ermB. Reconstruction of LTS_ICESag37 showed that, at 120Kbp, it was much longer than the 79Kbp ICESag37 reference sequence described previously [24]. This was due to the potential insertion of a region containing the phage assembly and recombination genes present in the Sag153 reference genome. However, LTS_ICESag37 was split across three contigs and therefore requires long-read sequencing to confirm such an insertion.

2.3. Distribution of Major Surface Proteins

Most notably, a surface protein required for hypervirulence, HvgA, was present in the genomes of all serotype III ST17 isolates (Figure 1). All except one isolate carried genes for one of the four proteins, alpha-C, alp1, alp2/3, or rib (Table 1), targeted by the protein-based GBS vaccine currently in trial. All serotype Ia isolates (n = 18) and most of the serotype V isolates (n = 8/9) carried a serine-rich repeat glycoprotein determinant, ssr1, while ssr2 was identified only in serotype III ST17 isolates that were also positive for hvgA (Table 1). 66% (28/42) of isolates carried two of the pilus genes, while 29% (12/42) of isolates carried PI-2A only (Table 1).

3. Discussion

This study represents the first in-depth genomic characterisation of GBS isolates from Lithuania. Here, we have identified four known serotypes, with serotypes Ia, III, and V predominating in agreement with studies of GBS carriage and disease in pregnant women [14], and of GBS carriage in pregnant women and neonates in Lithuania [20]. The isolates in our dataset showed high genetic diversity, covering 18 sequence types based on MLST profiles. These sequence types clustered into seven clonal complexes, of which CC23, CC17, and CC1 were most common, corroborating findings from previous studies of GBS carriage [7,14]. Notably, one serotype Ib ST1 isolate carried the macrolide resistance gene, ermB, with similar isolates having been characterised in Portugal as an emergent macrolide resistant lineage [25]. Furthermore, we identified one serotype Ia ST1 isolate and one serotype Ia ST1387 isolate belonging to CC1, indicating some serotype Ia isolates may have participated in a capsular switching event similar to other recombination events reported in the ST1 lineage [26]. However, none of the serotype Ia CC1 isolates in our dataset carried macrolide resistance genes, nor was serotype Ia the most common serotype belonging to CC1, with this cluster consisting primarily of serotype V isolates, as may have been anticipated [7,25].
Following local guidelines, antimicrobial susceptibility testing was not requested for every isolate. Of those tested, all showed susceptibility to benzylpenicillin, with no mutations in pbp1a or pbp2X detected, supporting the general assumption that GBS’s susceptibility to this antibiotic is high, despite studies having reported varying degrees of penicillin resistance in GBS [11,27]. Nine GBS isolates carried genes conveying resistance to erythromycin and clindamycin, in line with reports from other European countries including Serbia [28], Portugal [25], and the UK [29]. In one serotype III ST17 isolate, we detected a multi-drug resistant MGE carrying tetracycline (tetO), macrolide (ermB), and aminoglycoside (aph3-III and ant6-Ia) resistance genes. This MGE had a backbone of the well described ICE, ICESag37 [24], which has been previously identified in multi-drug resistant serotype III ST17 isolates [30]. That we identified a MGE carrying AMR and MDR determinants in our isolates underscores the capacity of GBS populations to exchange genetic material facilitating their adaptation to the selection pressures imposed by antimicrobials. This phenomenon has been described in other streptococcal species, such as Streptococcus pneumoniae [31]. Therefore, further studies are needed to track the genomic changes occurring in GBS populations.
Currently, two strong GBS vaccine candidates are being considered for implementation in vaccination programs for pregnant women to prevent infections in their neonates. One of these is a hexavalent polysaccharide vaccine, covering serotypes Ia, Ib, II, III, IV, and V [16,17], and the other is a multivalent adjuvanted protein vaccine (NCT03807245). Our findings suggest the hexavalent vaccine would offer full coverage in Lithuania, supporting findings from a recent review that reported 93–99% of cases of neonatal and adult infections would be covered by this vaccine [14]. A Minervax protein-based GBS vaccine is also in phase II clinical trials (https://www.minervax.com/product-pipeline/ accessed on 9 September 2022). As this vaccine exploits the N-terminal domains of the most prevalent alpha-like surface proteins (Alpha-C (bca), Rib, Alp1, and Alp2/3), and 97.6% (41/42) of our GBS isolates carried genes for one of these four proteins, it should offer a similarly high level of coverage in Lithuania. However, alpha-like surface proteins were detected here using only in silico genotyping methods that do not necessarily predict the expression of these proteins in isolates.
Unfortunately, too few isolates were investigated in this study to draw any consequential conclusions regarding the overall prevalence of specific AMR profiles and MGEs in the GBS population in Lithuania. Indeed, the serotypes and genotypes presented here may not represent the national picture of disease-causing GBS in Lithuania. A second limitation of this study relates to the ability to fully characterise MGEs using Illumina short reads. In future work, long-read sequencing technologies such as Nanopore sequencing, which can generate ultra-long (>4Mb) reads, may be implemented to better reconstruct and characterise the MGEs present in GBS isolates. Nonetheless, this study contributes a valuable first perspective on the MGEs present in the GBS population in Lithuania.
Taken together, our findings lay the foundations for future GBS surveillance studies in Lithuania, enabling future genomic and epidemiological comparisons to other GBS isolates that may inform national antimicrobial procedures to prevent GBS diseases in the same region. All such work is vital to better understanding the global spread of AMR among GBS and estimating the efficacy of GBS vaccines once available.

4. Materials and Methods

4.1. Bacterial Isolates, Microbiological and Demographic Data

Isolates were collected between July 2020–September 2020 from participating diagnostic laboratories in Kaunas region (Lithuania) as a part of Surveillance and Molecular Epidemiology of Maternal Infections Consortium. In total, 35 isolates were retrieved from vaginal swab samples, six from urine samples and one from a blood sample of an individual with confirmed bacteraemia (Table S1). Vaginal swabs, urine and blood samples were processed following Clinical Laboratory Standards Institute (CLSI) guidelines and local microbiology diagnostic procedures at the participating hospitals. GBS was identified by distinct β-haemolysis on Columbia agar with 5% Sheep Blood (Becton Dickinson GmbH, Heidelberg, Germany) and further confirmed by latex agglutination with a Streptococcal grouping kit (Oxoid, Dublin, Ireland). The following demographic information was collected (where available): age of patient, pregnancy status, sample type from which GBS was isolated, date of GBS isolation and results from routine antibiotic susceptibility testing. Antimicrobial susceptibility testing for penicillin, erythromycin, clindamycin and sulfadiazine-trimethoprim was performed using a Kirby-Bauer disk diffusion assays following the standard protocols [32] and results were interpreted using the EUCAST breakpoint table v.10 (https://www.eucast.org/fileadmin/src/media/PDFs/EUCAST_files/Breakpoint_tables/v_10.0_Breakpoint_Tables.pdf accessed on 1 November 2020). Phenotypic antibiotic susceptibility testing was carried out following local microbiology diagnostic procedures. Only some isolates were tested against all four antibiotics. Following the local Lithuanian screening routines, all pregnant females are routinely screened for Group B Streptococcus colonization between 35th and 37th week of pregnancy. In the absence of symptoms or risk factors, the attending physician, or the patients themselves, have the right to omit any further antimicrobial treatment at the time of testing and therefore the antimicrobial testing can be not requested by the physician. This is, however, subject to change in case new risk factors or symptoms emerge.

4.2. DNA Extraction and Whole Genome Sequencing

GBS isolates were streaked on Columbia Blood agar plates (Oxoid, Basingstoke, UK) and incubated at 37 °C 5% CO2 overnight. Genomic DNA from GBS bacterial isolates was extracted using the GenEluteTM bacterial Genomic DNA kit (Sigma-Aldrich, Burlington, MA, USA) following the manufacturer’s instructions for gram-positive bacteria with modifications as follows: GBS was lysed in 180 μL G+ lysis solution (comes with the GenEluteTM bacterial Genomic DNA kit) with added 20 μL mutanolysin (Sigma-Aldrich, US; prepared at 3000 U/mL) and 20 μL lysozyme (Sigma-Aldrich, USA; prepared at 100 mg/mL) prior to incubation at 37 °C for 1 h. Subsequent steps followed the manufacturer’s instructions. DNA was quantified using a NanoDrop spectrophotometer (Thermo, Waltham, MA, USA).
Multiplexed DNA libraries were prepared using Nextera XT DNA Library Preparation kit (Illumina, San Diego, CA, USA) following manufacturer’s instructions and whole genome sequenced on Illumina HiSeq 2500 System (Illumina) using 2 × 150 bp paired-end mode. Read demultiplexing and adapter trimming were performed by the Imperial BRC Genomics Facility using their production pipeline.

4.3. Genomic Analyses

Quality checks. Raw sequencing reads were checked with FastQC v. 0.11.9 (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) before raw reads were trimmed with Trimmomatic v.0.39 [33] using the following parameters: remove trailing low quality or N bases (below quality three); scan the read with a 5-base wide sliding window, cutting when the average quality per base drops below 30 and drop reads below 50 bases long. Bacterial species from assembled genomes was confirmed using Kraken2 v2.0.8-beta and its full bacterial database [34].
De novo assembly and annotation. Unicycler 0.4.8 [35] was used to generate consensus de novo assemblies. Assembly statistics were checked using QUAST v5.0.2 [36] and assembly graphs were visualised using Bandage [37]. Assemblies were annotated using Prokka v.1.13.4 [38].
Serotyping and multi-locus sequence typing (MLST). In-silico serotyping was done by iPCRess v2.2 (https://www.ebi.ac.uk/about/vertebrate-genomics/software/ipcress-manual) using GBS multiplex primers as reported by Imperi et al., 2010 [39]. MLST was assigned using mlst (Seemann T, mlst Github https://github.com/tseemann/mlst and the PubMLST database (https://pubmlst.org/). Novel allele sequences and allelic combinations were submitted to the PubMLST database (https://pubmlst.org/). Clonal complexes were assigned using the PubMLST database (https://pubmlst.org/).
Characterisation of surface proteins. Blastn [40] was used to screen for surface protein genes: alpha-C, alp1, alp2/3, rib, ssr1, ssr2 and PI-1, PI-1B, PI-2A (including all four allelic variants PI-2A1, PI-2A2, PI-2A3 and PI-2A4), PI-2B, and PI-2B2 using surface protein gene sequences described by McGee et al [7] as references. During Blastn analysis, sequence queries of ≥95% identity were reported as positive results, but we have also recorded surface protein matches with <95% identify, though we did not interpret these as positive results in this study.
Detection of antimicrobial resistance determinants. Acquired antimicrobial resistance genes were detected from assemblies using ABRicate (Seemann T, Abricate, Github https://github.com/tseemann/abricate) with the ARG-ANNOT [41] and ResFinder [42] databases. Chromosomal genes: gyrA, parC, pbp1a, and pbp2x, were investigated using Clustal Omega [43] for the presence of mutations conferring resistance to fluoroquinolones and beta-lactams, as previously reported by McGee et al. [7].
Characterisation of AMR-carrying mobile genetic elements. Annotated genomes were visualised with ARTEMIS [44]. Classification of known mobile genetic elements (MGEs) was completed by manually scanning the regions proximal to acquired AMR genes and subsequently identifying characteristic genes of previously described integrative and conjugative elements (ICEs): insertion sequences (IS), transposases (tnp), integrases (int), excisionases (xis), relaxases (mobC), helicases (recB), and/or recombinases. Characteristic genes involved in conjugation included those of Type IV-secretion systems and the vir family. Regions were also scanned for rumA and rplL genes which are recombination hotspots for ICE integration. Sequences of potential MGEs were compared to a reference found in TnRegistry (https://transposon.lstmed.ac.uk/tn-registry). In genomes where no ICE characterising genes were identified, genes adjacent to acquired AMR genes were investigated manually using ARTEMIS and their annotations checked using BLASTn. Figures were created with EasyFig [45] and the built-in BLASTn alignment tool.
Phylogenetic analysis. Single nucleotide polymorphisms (SNPs) were detected by mapping trimmed raw reads to a reference genome SS1 (CP010867.1) using snippy v4.6.0 (https://github.com/tseemann/snippy). The whole genome alignment was cleaned using the snippy-clean_full_aln function available as part of Snippy. Gubbins v2.4.1 [46] was used to identify recombination regions before processing alignment with SNP-sites [47] to obtain a recombination-free core SNPs file. The maximum-likelihood phylogenetic tree was constructed using IQ-TREE v.2.0.3 [48] with the “GTR+G+ASC” model (which combines the following models: general time reversible model with unequal rates and unequal base frequency (GTR model); discrete Gamma model with default 4 rate categories (G model) and an ascertainment bias correction model (ASC)) and 1000 bootstrap replicates. The maximum-likelihood phylogenetic tree was visualised using iTOL v.5 [49]. SNP pairwise distances were calculated using MEGA X v.10.2.4 [50]. To check for potential capsular switching for ST1 serotype Ia isolate, all available ST1 isolates reported as serotype Ia (by genotype) on the PubMLST database (https://pubmlst.org/) were collected (Table S2) and analysed in combination with serotype Ib ST1, serotype V ST1 and serotype Ia ST1387 isolates, constituting CC1. All isolates underwent the phylogenetic analysis described above. Recombination regions were visualised using Phandango [51] and BRIG [52].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens11091060/s1, Figure S1: A maximum-likelihood phylogenetic tree and comparison of recombination regions of CC1 isolates; Figure S2: A comparison of representative Tn916 and Tn3872 sequences identified in our dataset and reference Tn916 and Tn917 from Enterococcus faecalis; Table S1: Summary of the genomic, demographic and phenotypic characteristics of GBS isolates from Lithuania; Table S2: A list of nine serotype Ia ST1 genomes downloaded from the PubMLST database for comparison analysis of ST1/CC1 genomes from this study.

Author Contributions

Conceptualization, E.J.; methodology, E.J.; formal analysis, J.R., P.M., E.J.; investigation, R.P. (Ruta Prakapaite), R.G., R.P. (Rita Planciuniene), J.R., P.M.; resources, R.P. (Rita Planciuniene), R.G., P.K., E.J.; data curation, R.P. (Ruta Prakapaite), J.R., P.M., P.K., E.J.; writing—original draft preparation, J.R., E.J.; writing—review and editing, all co-authors; visualization, J.R., P.M., E.J.; supervision, P.K., E.J.; project administration, R.P. (Ruta Prakapaite), P.K., E.J.; funding acquisition, E.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Rosetrees Trust and the Stoneygate Trust (M683).

Institutional Review Board Statement

The isolates were collected as a discarded material after routine microbiological diagnostic analysis. The de-identification and double-coding practice was performed to ensure the anonymity of the patients. The study was approved by the Institutional Bioethics Committee (ILPMI-2020-12).

Informed Consent Statement

Not applicable.

Data Availability Statement

The whole genome raw sequence data and assemblies have been submitted to the European Nucleotide Archive (https://www.ebi.ac.uk/ena/), under the project accession PRJEB55405 (Table S1).

Acknowledgments

The authors want to thank to all the participating hospital centres for providing their isolates for this study. This publication made use of the PubMLST website (https://pubmlst.org/) developed by Keith Jolley (Jolley & Maiden 2010, BMC Bioinformatics, 11:595) and sited at the University of Oxford. The development of that website was funded by the Wellcome Trust. The Imperial BRC Genomics Facility has provided resources and support that have contributed to the research results reported within this paper. The Imperial BRC Genomics Facility is supported by NIHR funding to the Imperial Biomedical Research Centre. EJ is a Rosetrees/Stoneygate 2017 Imperial College Research Fellow, funded by Rosetrees Trust and the Stoneygate Trust (Fellowship no. M683). EJ and JR are affiliated with the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Healthcare Associated Infections and Antimicrobial Resistance at Imperial College London in partnership with the UK Health Security Agency (previously PHE), in collaboration with, Imperial Healthcare Partners, University of Cambridge and University of Warwick.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Barcaite, E.; Bartusevicius, A.; Tameliene, R.; Kliucinskas, M.; Maleckiene, L.; Nadisauskiene, R. Prevalence of Maternal Group B Streptococcal Colonisation in European Countries. Acta Obstet. Gynecol. Scand. 2008, 87, 260–271. [Google Scholar] [CrossRef] [PubMed]
  2. Russell, N.J.; Seale, A.C.; O’Driscoll, M.; O’Sullivan, C.; Bianchi-Jassir, F.; Gonzalez-Guarin, J.; Lawn, J.E.; Baker, C.J.; Bartlett, L.; Cutland, C.; et al. Maternal Colonization with Group B Streptococcus and Serotype Distribution Worldwide: Systematic Review and Meta-Analyses. Clin. Infect. Dis. 2017, 65, S100–S111. [Google Scholar] [CrossRef] [PubMed]
  3. Le Doare, K.; Heath, P.T. An Overview of Global GBS Epidemiology. Vaccine 2013, 31, D7–D12. [Google Scholar] [CrossRef]
  4. Puopolo, K.M.; Lynfield, R.; Cummings, J.J. Management of Infants at Risk for Group B Streptococcal Disease. Pediatrics 2019, 144, e20191881. [Google Scholar] [CrossRef] [PubMed]
  5. Hayes, K.; O’Halloran, F.; Cotter, L. A Review of Antibiotic Resistance in Group B Streptococcus: The Story so Far. Crit. Rev. Microbiol. 2020, 46, 253–269. [Google Scholar] [CrossRef]
  6. Metcalf, B.J.; Chochua, S.; Gertz, R.E.; Hawkins, P.A.; Ricaldi, J.; Li, Z.; Walker, H.; Tran, T.; Rivers, J.; Mathis, S.; et al. Short-Read Whole Genome Sequencing for Determination of Antimicrobial Resistance Mechanisms and Capsular Serotypes of Current Invasive Streptococcus Agalactiae Recovered in the USA. Clin. Microbiol. Infect. 2017, 23, 574.e7–574.e14. [Google Scholar] [CrossRef] [PubMed]
  7. McGee, L.; Chochua, S.; Li, Z.; Mathis, S.; Rivers, J.; Metcalf, B.; Ryan, A.; Alden, N.; Farley, M.M.; Harrison, L.H.; et al. Multistate, Population-Based Distributions of Candidate Vaccine Targets, Clonal Complexes, and Resistance Features of Invasive Group B Streptococci within the United States, 2015-2017. Clin. Infect. Dis. 2021, 72, 1004–1013. [Google Scholar] [CrossRef]
  8. Chu, Y.W.; Tse, C.; Tsang, G.K.L.; So, D.K.S.; Fung, J.T.L.; Lo, J.Y.C. Invasive Group B Streptococcus Isolates Showing Reduced Susceptibility to Penicillin in Hong Kong. J. Antimicrob. Chemother. 2007, 60, 1407–1409. [Google Scholar] [CrossRef]
  9. Kimura, K.; Suzuki, S.; Wachino, J.I.; Kurokawa, H.; Yamane, K.; Shibata, N.; Nagano, N.; Kato, H.; Shibayama, K.; Arakawa, Y. First Molecular Characterization of Group B Streptococci with Reduced Penicillin Susceptibility. Antimicrob. Agents Chemother. 2008, 52, 2890–2897. [Google Scholar] [CrossRef]
  10. Nagano, N.; Nagano, Y.; Kimura, K.; Tamai, K.; Yanagisawa, H.; Arakawa, Y. Genetic Heterogeneity in Pbp Genes among Clinically Isolated Group B Streptococci with Reduced Penicillin Susceptibility. Antimicrob. Agents Chemother. 2008, 52, 4258–4267. [Google Scholar] [CrossRef] [Green Version]
  11. Seki, T.; Kimura, K.; Reid, M.E.; Miyazaki, A.; Banno, H.; Jin, W.; Wachino, J.I.; Yamada, K.; Arakawa, Y. High Isolation Rate of MDR Group B Streptococci with Reduced Penicillin Susceptibility in Japan. J. Antimicrob. Chemother. 2015, 70, 2725–2728. [Google Scholar] [CrossRef] [PubMed]
  12. Dahesh, S.; Hensler, M.E.; Van Sorge, N.M.; Gertz, R.E.; Schrag, S.; Nizet, V.; Beall, B.W. Point Mutation in the Group B Streptococcal Pbp2x Gene Conferring Decreased Susceptibility to β-Lactam Antibiotics. Antimicrob. Agents Chemother. 2008, 52, 2915–2918. [Google Scholar] [CrossRef] [PubMed]
  13. Kimura, K.; Nagano, N.; Arakawa, Y. Classification of Group B Streptococci with Reduced β-Lactam Susceptibility (GBS-RBS) Based on the Amino Acid Substitutions in PBPs. J. Antimicrob. Chemother. 2014, 70, 1601–1603. [Google Scholar] [CrossRef] [PubMed]
  14. Bianchi-Jassir, F.; Paul, P.; To, K.N.; Carreras-Abad, C.; Seale, A.C.; Jauneikaite, E.; Madhi, S.A.; Russell, N.J.; Hall, J.; Madrid, L.; et al. Systematic Review of Group B Streptococcal Capsular Types, Sequence Types and Surface Proteins as Potential Vaccine Candidates. Vaccine 2020, 38, 6682–6694. [Google Scholar] [CrossRef] [PubMed]
  15. Madhi, S.A.; Cutland, C.L.; Jose, L.; Koen, A.; Govender, N.; Wittke, F.; Olugbosi, M.; Meulen, A.S.-t.; Baker, S.; Dull, P.M.; et al. Safety and Immunogenicity of an Investigational Maternal Trivalent Group B Streptococcus Vaccine in Healthy Women and Their Infants: A Randomised Phase 1b/2 Trial. Lancet Infect. Dis. 2016, 16, 923–934. [Google Scholar] [CrossRef]
  16. Buurman, E.T.; Timofeyeva, Y.; Gu, J.; Kim, J.H.; Kodali, S.; Liu, Y.; Mininni, T.; Moghazeh, S.; Pavliakova, D.; Singer, C.; et al. A Novel Hexavalent Capsular Polysaccharide Conjugate Vaccine (GBS6) for the Prevention of Neonatal Group B Streptococcal Infections by Maternal Immunization. J. Infect. Dis. 2019, 220, 105–115. [Google Scholar] [CrossRef] [PubMed]
  17. Absalon, J.; Segall, N.; Block, S.L.; Center, K.J.; Scully, I.L.; Giardina, P.C.; Peterson, J.; Watson, W.J.; Gruber, W.C.; Jansen, K.U.; et al. Safety and Immunogenicity of a Novel Hexavalent Group B Streptococcus Conjugate Vaccine in Healthy, Non-Pregnant Adults: A Phase 1/2, Randomised, Placebo-Controlled, Observer-Blinded, Dose-Escalation Trial. Lancet Infect. Dis. 2021, 21, 263–274. [Google Scholar] [CrossRef]
  18. Fischer, P.; Pawlowski, A.; Cao, D.; Bell, D.; Kitson, G.; Darsley, M.; Johansson-Lindbom, B. Safety and Immunogenicity of a Prototype Recombinant Alpha-like Protein Subunit Vaccine (GBS-NN) against Group B Streptococcus in a Randomised Placebo-Controlled Double-Blind Phase 1 Trial in Healthy Adult Women. Vaccine 2021, 39, 4489–4499. [Google Scholar] [CrossRef]
  19. Pawlowski, A.; Lannergård, J.; Gonzalez-Miro, M.; Cao, D.; Larsson, S.; Persson, J.J.; Kitson, G.; Darsley, M.; Rom, A.L.; Hedegaard, M.; et al. A Group B Streptococcus Alpha-like Protein Subunit Vaccine Induces Functionally Active Antibodies in Humans Targeting Homotypic and Heterotypic Strains. Cell Reports Med. 2022, 3, 100511. [Google Scholar] [CrossRef]
  20. Barcaite, E.; Bartusevicius, A.; Tameliene, R.; Maleckiene, L.; Vitkauskiene, A.; Nadisauskiene, R. Group B Streptococcus and Escherichia Coli Colonization in Pregnant Women and Neonates in Lithuania. Int. J. Gynecol. Obstet. 2012, 117, 69–73. [Google Scholar] [CrossRef]
  21. Mockevičiūtė, E. Naujagimių Įgimtos B Grupės Beta Hemolizinio Streptokoko Sukeltos Infekcijos Skirtingų Profilaktikos Metodikų Taikymo Analizė. (Analysis of different prophylaxis measures against early neonate GBS infection). Master’s Thesis, Lithuanian Health Sciences University, Kaunas, Lithuania, 2017; 25p. [Google Scholar]
  22. Abraitis, V.; Arlauskienė, A.; Zakarevičienė, J.; Bagušytė, L.; Valkerienė, G.; Barčaitė, E.; Bartkevičienė, D.; Biržietis, T.; Bumbulienė, Ž.; Drejerienė, E.; et al. Metodika B Grupės β Hemolizinis Streptokokas (BGS) Ir Kita Bakterinė Prenatalinė Infekcija. (Methods for Group B Streptococcus (GBS) and other neonatal bacterial infections). Liet. Ir Šveicarijos Bendradarbiavimo Programa 2014, 58, 1–53. [Google Scholar]
  23. Kawamura, Y.; Fujiwara, H.; Mishima, N.; Tanaka, Y.; Tanimoto, A.; Ikawa, S.; Itoh, Y.; Ezaki, T. First Streptococcus Agalactiae Isolates Highly Resistant to Quinolones, with Point Mutations in GyrA and ParC. Antimicrob. Agents Chemother. 2003, 47, 3605–3609. [Google Scholar] [CrossRef] [PubMed]
  24. Zhou, K.; Xie, L.; Han, L.; Guo, X.; Wang, Y.; Sun, J. ICESag37, a Novel Integrative and Conjugative Element Carrying Antimicrobial Resistance Genes and Potential Virulence Factors in Streptococcus Agalactiae. Front. Microbiol. 2017, 8, 1921. [Google Scholar] [CrossRef] [PubMed]
  25. Lopes, E.; Fernandes, T.; Machado, M.P.; Carriço, J.A.; Melo-Cristino, J.; Ramirez, M.; Martins, E.R.; Oliveira, H.; Vaz, T.; Gião, M.; et al. Increasing Macrolide Resistance among Streptococcus Agalactiae Causing Invasive Disease in Non-Pregnant Adults Was Driven by a Single Capsular-Transformed Lineage, Portugal, 2009 to 2015. Eurosurveillance 2018, 23, pii=1700473. [Google Scholar] [CrossRef] [PubMed]
  26. Neemuchwala, A.; Teatero, S.; Athey, T.B.T.; McGeer, A.; Fittipaldi, N. Capsular Switching and Other Large-Scale Recombination Events in Invasive Sequence Type 1 Group B Streptococcus. Emerg. Infect. Dis. 2016, 22, 1941–1944. [Google Scholar] [CrossRef]
  27. Kasahara, K.; Baltus, A.J.; Lee, S.H.; Edelstein, M.A.; Edelstein, P.H. Prevalence of Non-Penicillin-Susceptible Group B Streptococcus in Philadelphia and Specificity of Penicillin Resistance Screening Methods. J. Clin. Microbiol. 2010, 48, 1468–1469. [Google Scholar] [CrossRef]
  28. Kekic, D.; Gajic, I.; Opavski, N.; Kojic, M.; Vukotic, G.; Smitran, A.; Boskovic, L.; Stojkovic, M.; Ranin, L. Trends in Molecular Characteristics and Antimicrobial Resistance of Group B Streptococci: A Multicenter Study in Serbia, 2015–2020. Sci. Rep. 2021, 11, 540. [Google Scholar] [CrossRef]
  29. Khan, U.B.; Jauneikaite, E.; Andrews, R.; Chalker, V.J.; Spiller, O.B. Identifying Large-Scale Recombination and Capsular Switching Events in Streptococcus Agalactiae Strains Causing Disease in Adults in the UK between 2014 and 2015. Microb. Genomics 2022, 8, 000783. [Google Scholar] [CrossRef]
  30. Hsu, J.F.; Tsai, M.H.; Lin, L.C.; Chu, S.M.; Lai, M.Y.; Huang, H.R.; Chiang, M.C.; Yang, P.H.; Lu, J.J. Genomic Characterization of Serotype III/ST-17 Group b Streptococcus Strains with Antimicrobial Resistance Using Whole Genome Sequencing. Biomedicines 2021, 9, 1477. [Google Scholar] [CrossRef]
  31. Croucher, N.J.; Chewapreecha, C.; Hanage, W.P.; Harris, S.R.; McGee, L.; Van Der Linden, M.; Song, J.H.; Ko, K.S.; De Lencastre, H.; Turner, C.; et al. Evidence for Soft Selective Sweeps in the Evolution of Pneumococcal Multidrug Resistance and Vaccine Escape. Genome Biol. Evol. 2014, 6, 1589–1602. [Google Scholar] [CrossRef]
  32. Hudzicki, J.; Kirby, B. Disk Diffusion Susceptibility Test Protocol. Am. Soc. Microbiol. 2016, 1–23. [Google Scholar]
  33. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A Flexible Trimmer for Illumina Sequence Data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wood, D.E.; Lu, J.; Langmead, B. Improved Metagenomic Analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [PubMed]
  35. Wick, R.R.; Judd, L.M.; Gorrie, C.L.; Holt, K.E. Unicycler: Resolving Bacterial Genome Assemblies from Short and Long Sequencing Reads. PLoS Comput. Biol. 2017, 13, e1005595. [Google Scholar] [CrossRef]
  36. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality Assessment Tool for Genome Assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  37. Wick, R.R.; Schultz, M.B.; Zobel, J.; Holt, K.E. Bandage: Interactive Visualization of de Novo Genome Assemblies. Bioinformatics 2015, 31, 3350–3352. [Google Scholar] [CrossRef]
  38. Seemann, T. Prokka: Rapid Prokaryotic Genome Annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef]
  39. Imperi, M.; Pataracchia, M.; Alfarone, G.; Baldassarri, L.; Orefici, G.; Creti, R. A Multiplex PCR Assay for the Direct Identification of the Capsular Type (Ia to IX) of Streptococcus Agalactiae. J. Microbiol. Methods 2010, 80, 212–214. [Google Scholar] [CrossRef]
  40. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  41. Gupta, S.K.; Padmanabhan, B.R.; Diene, S.M.; Lopez-Rojas, R.; Kempf, M.; Landraud, L.; Rolain, J.M. ARG-Annot, a New Bioinformatic Tool to Discover Antibiotic Resistance Genes in Bacterial Genomes. Antimicrob. Agents Chemother. 2014, 58, 212–220. [Google Scholar] [CrossRef]
  42. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of Acquired Antimicrobial Resistance Genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef] [PubMed]
  43. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, Scalable Generation of High-Quality Protein Multiple Sequence Alignments Using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef] [PubMed]
  44. Carver, T.; Harris, S.R.; Berriman, M.; Parkhill, J.; McQuillan, J.A. Artemis: An Integrated Platform for Visualization and Analysis of High-Throughput Sequence-Based Experimental Data. Bioinformatics 2012, 28, 464–469. [Google Scholar] [CrossRef]
  45. Sullivan, M.J.; Petty, N.K.; Beatson, S.A. Easyfig: A Genome Comparison Visualizer. Bioinformatics 2011, 27, 1009–1010. [Google Scholar] [CrossRef] [PubMed]
  46. Croucher, N.J.; Page, A.J.; Connor, T.R.; Delaney, A.J.; Keane, J.A.; Bentley, S.D.; Parkhill, J.; Harris, S.R. Rapid Phylogenetic Analysis of Large Samples of Recombinant Bacterial Whole Genome Sequences Using Gubbins. Nucleic Acids Res. 2015, 43, e15. [Google Scholar] [CrossRef]
  47. Page, A.J.; Taylor, B.; Delaney, A.J.; Soares, J.; Seemann, T.; Keane, J.A.; Harris, S.R. SNP-Sites: Rapid Efficient Extraction of SNPs from Multi-FASTA Alignments. Microb. genomics 2016, 2, e000056. [Google Scholar] [CrossRef]
  48. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; Von Haeseler, A.; Lanfear, R.; Teeling, E. IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef]
  49. Letunic, I.; Bork, P. Interactive Tree of Life (ITOL) v4: Recent Updates and New Developments. Nucleic Acids Res. 2019, 47, W256–W259. [Google Scholar] [CrossRef]
  50. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  51. Hadfield, J.; Croucher, N.J.; Goater, R.J.; Abudahab, K.; Aanensen, D.M.; Harris, S.R. Phandango: An Interactive Viewer for Bacterial Population Genomics. Bioinformatics 2018, 34, 292–293. [Google Scholar] [CrossRef]
  52. Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. BLAST Ring Image Generator (BRIG): Simple Prokaryote Genome Comparisons. BMC Genomics 2011, 12, 402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Maximum likelihood phylogenetic tree for GBS isolates from women in Lithuania. Maximum likelihood mid-point rooted phylogenetic tree constructed using recombination-free alignment of core-genome SNPs and the reference genome, SS1 (accession number: CP010867), representing serotype V ST1. The first data column shows the ST of each isolate per MLST results; the second data column shows the CC of each isolate. Branch colour shows the serotype of each isolate: Ia (blue), Ib (purple), III (red), and V (yellow). Branch shape shows the sample type of each isolate: vaginal (circle), urine (star), and blood (square). The next six data columns show the presence (filled circles) and absence (unfilled circles) in each isolate of the acquired AMR genes, tetM, tetO, ermA, ermB, lsaC, aph3-III, and ant6-Ia. The next three data columns show the susceptibility (S) or resistance (R) of the isolates to the antimicrobials, benzylpenicillin (PEN), erythromycin (ERY), and clindamycin (CLIN). Empty space in these columns indicates that the antibiotic susceptibility testing result was not required for treatment purposes and therefore not available. The final ten data columns show the presence (filled circles) and absence (unfilled circles) in each isolate of the surface proteins and virulence genes, alpha-C (or bca), rib, alp1, alp2/3, srr1, srr2, hvgA, PI-1, PI-2A, and PI-2B. Tree scale indicates nucleotide substitution rate per site.
Figure 1. Maximum likelihood phylogenetic tree for GBS isolates from women in Lithuania. Maximum likelihood mid-point rooted phylogenetic tree constructed using recombination-free alignment of core-genome SNPs and the reference genome, SS1 (accession number: CP010867), representing serotype V ST1. The first data column shows the ST of each isolate per MLST results; the second data column shows the CC of each isolate. Branch colour shows the serotype of each isolate: Ia (blue), Ib (purple), III (red), and V (yellow). Branch shape shows the sample type of each isolate: vaginal (circle), urine (star), and blood (square). The next six data columns show the presence (filled circles) and absence (unfilled circles) in each isolate of the acquired AMR genes, tetM, tetO, ermA, ermB, lsaC, aph3-III, and ant6-Ia. The next three data columns show the susceptibility (S) or resistance (R) of the isolates to the antimicrobials, benzylpenicillin (PEN), erythromycin (ERY), and clindamycin (CLIN). Empty space in these columns indicates that the antibiotic susceptibility testing result was not required for treatment purposes and therefore not available. The final ten data columns show the presence (filled circles) and absence (unfilled circles) in each isolate of the surface proteins and virulence genes, alpha-C (or bca), rib, alp1, alp2/3, srr1, srr2, hvgA, PI-1, PI-2A, and PI-2B. Tree scale indicates nucleotide substitution rate per site.
Pathogens 11 01060 g001
Table 1. Summary of the surface protein genes identified in the 42 Lithuanian GBS isolates.
Table 1. Summary of the surface protein genes identified in the 42 Lithuanian GBS isolates.
GBS Serotype (n)alpha-Calp1alp2/3ribsrr1srr2PI-1PI-2API-2BhvgA
Ia (n = 18)2475180101710
Ib (n = 1)0010101100
III (n = 14)000134912599
V (n = 9)4041805900
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rodgus, J.; Prakapaite, R.; Mitsidis, P.; Grigaleviciute, R.; Planciuniene, R.; Kavaliauskas, P.; Jauneikaite, E. Molecular Epidemiology of Group B Streptococci in Lithuania Identifies Multi-Drug Resistant Clones and Sporadic ST1 Serotypes Ia and Ib. Pathogens 2022, 11, 1060. https://doi.org/10.3390/pathogens11091060

AMA Style

Rodgus J, Prakapaite R, Mitsidis P, Grigaleviciute R, Planciuniene R, Kavaliauskas P, Jauneikaite E. Molecular Epidemiology of Group B Streptococci in Lithuania Identifies Multi-Drug Resistant Clones and Sporadic ST1 Serotypes Ia and Ib. Pathogens. 2022; 11(9):1060. https://doi.org/10.3390/pathogens11091060

Chicago/Turabian Style

Rodgus, Jonah, Ruta Prakapaite, Panagiotis Mitsidis, Ramune Grigaleviciute, Rita Planciuniene, Povilas Kavaliauskas, and Elita Jauneikaite. 2022. "Molecular Epidemiology of Group B Streptococci in Lithuania Identifies Multi-Drug Resistant Clones and Sporadic ST1 Serotypes Ia and Ib" Pathogens 11, no. 9: 1060. https://doi.org/10.3390/pathogens11091060

APA Style

Rodgus, J., Prakapaite, R., Mitsidis, P., Grigaleviciute, R., Planciuniene, R., Kavaliauskas, P., & Jauneikaite, E. (2022). Molecular Epidemiology of Group B Streptococci in Lithuania Identifies Multi-Drug Resistant Clones and Sporadic ST1 Serotypes Ia and Ib. Pathogens, 11(9), 1060. https://doi.org/10.3390/pathogens11091060

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop