Skip to main content

Phylogenetic, comparative genomic and structural analyses of human Streptococcus agalactiae ST485 in China



Streptococcus agalactiae (Group B Streptococcus, GBS) is a common bacteria species infecting both human and bovine. Previous studies have shown that the GBS isolated from human and bovine are mostly unrelated and belong to separate populations. However, recently, the bovine GBS CC103 has become the dominant epidemic strain and frequently isolated from human patients. In particular, the ST485 GBS, a member of CC103, has become the new dominant ST in China and exhibited very high pathogenicity. This phenomenon is not consistent with the established understanding about the relationship between bovine and human GBS, which needs to be re-investigated.


The genome-based phylogenetic analysis showed that the human and bovine GBS CC103 strains had very close genetic relationship and they were alternately distributed on the evolutionary tree. CC103 strains evolved into several branches, including the ST485, which exhibited high pathogenicity and specifically infected human. Compared to other CC103 strains, the ST485 lacked Lac.2 gene structure and acquired the CadDX gene structure in their genomes.


Our results indicate that GBS CC103 could propagate across human and bovine, and GBS ST485 might evolve from the ST103 that could infect both human and bovine. Moreover, the recombination of Lac.2 and CadDX gene structures might play an important role in the formation of highly pathogenic ST485 in China.


Streptococcus agalactiae, also called Group B streptococcus (GBS), is associated with early and late onset diseases in infants. It colonizes the urogenital and gastro-intestinal tracts without causing symptoms, and leads to diseases like septicemia in non-pregnant adults [1, 2]. The multi-locus sequence typing (MLST) of GBS strains isolated from different countries showed that most human-derived strains and clinical isolates were clustered into five major clone complexes (CC) (CC1, CC10, CC17, CC19 and CC23) [3], and the majority of bovine isolates belonged to another clone complex, CC67 [4]. However, recently, in some regions, the human dominant strains (such as CC1 and CC23) began to spread in bovine, and replaced CC67 to become the most widespread CC causing cow mastitis. Moreover, the newly appeared GBS CC103 strains have become the major strains causing cow mastitis in some countries of Europe and Asia [5, 6]. Recently (2015–2017), the isolation frequency of human GBS CC103 significantly increased in China (from 1.25 to 21.74%), especially the sequence type (ST) 485 (from 1.25 to 14.13%), which has become the new dominant ST and exhibits 10–20 times higher pathogenicity than the strains outside of China [7,8,9,10,11]. In this study, we examined the evolution of human and bovine CC103 strains by performing phylogenetic, comparative genomic and structural analyses, and investigated the origin and causes of the highly pathogenic GBS ST485 in China.


GBS isolates, genome sequencing and annotation

The draft genomes of 18 isolated GBS CC103 strains were determined using Illumina HiSeq2000 sequencing platforms, and then they were assembled by the ABySS program [12]. The minimal coverage was 500-fold. Since all the current public databases do not have the full-length genome data of ST485 strains, we used primer walking to close the gaps in the draft genome of ST485 strain BSE009 (one of the above 18 CC103 strains), and the resulting PCR products were sequenced to generate the whole genome. The assembled sequences were uploaded to the RAST website for gene function annotation and metabolic pathway construction ( In addition, the genomes of another 52 GBS strains in major CC groups, which came from different hosts, were selected for evolution analysis. The strains involved in this study are listed in Additional file 1, which includes the information about their CC and ST, the capsular serotype, host, source, associated diseases, year of isolation, geographical origin, and genome accession number. The 18 CC103 strains were obtained directly from patients, and the information about the isolation procedure was reported in our previous article [11]. All subjects provided written informed consent before their inclusion in the study.

Phylogenetic analysis and CRISPRs analysis

OrthoMCL was used to record the orthologous protein sequences among the isolates [13]. MAFFT was used to perform the multi-sequence alignment for single copy homologous proteins [14], and the poorly aligned positions and divergent regions were removed. ProtTest was used to do the Maximum Likelihood (ML) Estimation for phylogenetic tree, and the model parameters were obtained from PhyML [15]. AIC and BIC scores were evaluated to obtain the optimal amino acid substitution model. ML method was used to construct a phylogenetic tree with 1,000 bootstrap replications using RaxML software [16]. CRISPRs finder and CRISPR recognition tool (CRT) were used to identify the CRISPR sequences (clustered regularly interspaced short palindromic repeats) [17, 18]. Each unique spacer was numbered manually, and the result was analyzed and modified according to Lier et al. [19].

Comparison of whole genome sequences and functional genes

The genome sequences and functional coding genes of GBS strains were compared by the sequence alignment and functional gene alignment functions in RAST Server [20]. The results were confirmed by comparing with public databases (nr/nt). The alignment of GBS Lac.2 and CadDX gene structures, and the construction of phylogenetic trees were performed with MEGA7 using the ML method [21].


Genome-based phylogenic analysis found human-bovine co-infecting GBS strains

Considering the strain hosts and their serotypes, we chose at least 3 representative strains from each CC for phylogenetic analysis. Based on the MLST studies, a total of 70 GBS strains were selected to represent the known diversity of GBS population, which included 46 human strains, 12 bovine strains and 9 fish strains. Three strains from the genus Streptococcus were chosen as outgroup strains, and their genomic information was available in GenBank. The strains used in the study are listed in Additional file 1. The Maximum Likelihood (ML) method was used to construct the phylogenetic tree (Fig. 1), with the information from 96,665 amino acids of 391 single copy orthology clusters. We found that the 70 isolates were clustered into nine well-resolved lineages, which corresponded to the MLST defined CCs (Fig. 1). CC103 strains were clustered together into one branch. To more accurately reflect the evolutionary relationship within CC103 strains, we constructed the phylogenetic tree for CC103 (Fig. 2) using the information from 328,013 amino acids of 1,104 single copy orthology clusters of the 26 CC103 strains. The evolution result of CC103 strains was consistent with the CC103 ST clustering results. The bovine and human ST103 strains were alternately distributed on the evolutionary tree. Interestingly, the human dominant strain ST485 that emerged in China had a common ancestor with the bovine ST103 strain MRI-Z1–023. In addition, the ST485 strains were clustered together and had very close genetic distance to each other. Among all the isolated human CC103 clinical strains, ST485 accounted for 67% (12/18), indicating an evolutionary bottleneck for human CC103 strains.

Fig. 1

The phylogenetic relationships among GBS strains from different CC groups. On the left is the phylogenetic tree of 70 GBS strains, which was constructed with 1,000 bootstrap replications and rooted by outgroup. On the right is the information about strain host, serotype and ST. The major clonal complexes (CC) 1, 10, 17, 19, 23, 67, 102 and 552, which were defined by GBS MLST website (, were all separated by distinct branches. H, Homo sapiens. B, Bos Taurus. F, fish. C, Canis lupus familiaris. T, Tursiops. adhP, pheS, atr, glnA, sdhA, glcK and tkt are the names of 7 allelic genes in MLST analysis, and below are the numbers of each allelic gene in MLST website. The colors (or absence of color) are used to emphasize certain regions of the chart

Fig. 2

The phylogenetic analysis and gene structure comparison of 26 GBS CC103 strains. a The phylogenetic analysis of 26 GBS CC103 strains based on genome sequences, and the phylogeny was constructed with 1,000 bootstrap replications and rooted by the result of Fig. 1. A few strains with non-frequent ST and serotypes are highlighted in red. b Comparison of CRISPR1 loci. Internal repeats are not included; only terminal repeats (RT) and spacers are represented. The spacers are numbered, and same number represents the same sequence. The colors are used to emphasize the major spacers for easier viewing, and the same spacers are represented by same colors. c Gene structure comparisons of Lac.2 and CadDX. Lac.2–1 and Lac.2–2 are two types of Lac.2. +, the strain has this gene; −, the strain does not have this gene structure. *, there’s a gap at the CRISPR1 sequence

CRISPRs analysis further revealed the characteristics of GBS CC103 strains

CRISPRs (clustered regularly interspaced short palindromic repeats) are a bacterial adaptive immune defense mechanism to protect against the deadly outcomes from mobile genetic elements (MGEs) [22]. CRISPRs are a family of noncontiguous DNA repeats interspaced by unique spacer sequences with constant length [23]. The repeats are highly conserved within a CRISPR array, but the 3′ terminal repeat (RT) is specific for different GBS CCs. Spacer sequence corresponds to the previous exposure of MGEs, and it is inserted at the CRISPR leading end. The closer the spacer to the 3′ terminal, the more conserved the spacer sequence is [24]. Therefore, CRISPRs analysis can be used to investigate the phylogenetic relationships between subspecies [22]. There are two CRISPRs in GBS genome: Type 1-C CRISPR2, only present in a few strains, and type 2-A CRISPR1, ubiquitously present in all strains [25]. As shown in Fig. 2, except for one strain (IP43) that has sequence gap at CRISPR1, all the other 25 strains had intact CRISPR1 sequence structures. All the CC103 strains had consistent terminal repeats (RT) and terminal spacers (ST), which were also similar to the RT and ST of ST22 strains (Additional file 2). The three strains that had serotype variation exhibited the most distinct CRISPR1 sequences, with only two identical spacer sequences as other CC103 strains. There was no significant host specificity in ST103 strains, and their CRISPR sequences mostly contained No. 196–198 spacers, as compared to ST485 strains. The CRISPR1 sequences of ST485 strains were highly conserved, which again indicates an evolutionary bottleneck.

Genomic coding sequence alignment showed mutation characteristics of GBS CC103 strains in evolution

We compared the genomic coding sequences of GBS major CCs, CC103 strains and ST485 strains, and used BSE009 (ST485), COH1 (ST17) and C001 (ST103) as the reference strains, respectively. As shown in Fig. 3, there were significant differences in gene sequences between different CCs (Fig. 3a and b). According to the genome sequence of C001 strain, there were mainly four mutant gene islands in CC103 strains (Fig. 3c). In addition to these four mutant gene islands, ST485 strains also carried two mutation enriched regions, Y1 and Y2 (Fig. 3d). Based on the genome sequence of BSE009 strain, there were four mutant gene islands in ST485 strains (Fig. 3e). In addition to these four mutant gene islands, CC103 strains also carried four mutation enriched regions (Fig. 3f). The R1-R6, S1-S4, Q2-Q4, P1-P4, X1 and X2 regions mainly contained gene deletions (Fig. 3). The sequences in these mutant gene islands were mainly coded for Phage-related proteins, hypothetical proteins with unknown function, and gene recombination related enzymes. Y1, Y2, X3 and X4 regions mainly contained point mutations and the sequence identity was generally above 97%. Most of the genes in these regions were metabolism and transport related genes. Q1 region contained both gene deletions and point mutations, and it mainly included Lac.2 (lactose operon) related genes and glucose metabolism related genes.

Fig. 3

Genome coding sequence alignment of GBS strains. Each circle represents the genome of one GBS strain, and the strain name and ST/CC are labeled. The circle correspond to reference genomes are not shown. The major variation regions are numbered, such as R1, R2 ... The color bar on the bottom indicates the percentage of protein sequence identity compared to the reference genome. The positions of Lac.2 and CadDX are marked

Fig. 4

Characteristics of Lac.2 from GBS strains. a The phylogenetic tree of Lac.2 was constructed with MEGA7 using ML method. GBS strains are labeled with strain names, other Streptococcus strains are labeled with species and strain names. b Genetic organization of the Lac.2. The coding sequences of all genes in Lac.2–1/2 are depicted as arrows in different colors. The strains in A and B regions have Lac.2–2 and Lac.2–1, respectively. The branch lengths are labeled

The recombination of lac.2 and CadDX gene structures may be the key factors contributing to the emergence of highly pathogenic ST485 strains

We compared the functional genes of metabolically reprogrammed CC103 strains, and screened for different gene structures. The results showed that all ST485 strains did not have the Lac.2 gene structure (a cluster of genes that are responsible for transport and lactose metabolism), but acquired the CadDX gene structure (cadD and cadX that are related to cadmium resistance) (Fig. 2). On the contrary, seven ST103 strains had Lac.2 structure and two had CadDX. We then did sequential evolution analysis on the Lac.2 gene structure of 21 GBS strains and five Streptococcus strains, and confirmed the sequences by comparing to public databases (nr/nt). The results showed that the Lac.2 in CC103 strains had two types, Lac.2–1 and Lac.2–2. The main difference between these two types was that Lac.2–1 had an extra lacT gene between lacD and lacF loci and lacked hpB gene between lacR and lacA loci; additionally, the directions of lacR genes are different between these two types (Fig. 4). As expected, Lac.2–1 and Lac.2–2 formed different evolutionary branches. From all the GBS strains we tested, only four human strains had Lac.2; all the other strains carrying Lac.2 were from bovine (11 strains) or camel (1 strain).

The ST103 and ST485 strains in later evolutionary branches (Fig. 2), including the MRI Z1–023, all had cadmium resistance gene structure CadDX, which encoded for cadmium resistance proteins and cadmium efflux system accessory proteins. We did phylogenetic analysis on the CadDX sequences from 12 GBS strains with different ST and seven strains from the genus Streptococcus. The results showed that eight GBS strains were found to be highly homologous in CadDX gene and clustered together, including ST23, ST1, ST17, ST22, ST6, etc.; but the CadDX sequences of ST485 were more similar to the bovine ST67 strains (Fig. 5).

Fig. 5

The phylogenetic tree of CadDX constructed with MEGA7 using ML method. GBS strains are labeled with strain names, other Streptococcus strains are labeled with species and strain names. The branch lengths are labeled


GBS can infect many species, but human and bovine are the major hosts. Whether there is cross-infection of GBS between human and bovine is important for studying GBS related diseases. Previous studies showed that the GBS strains isolated from human were highly distinct from the bovine strains [4, 26,27,28,29,30,31,32,33,34,35]. However, the GBS CC103 strains were found to be able to infect many animals, including human, dog, cat, cattle, fish, etc. [5, 36,37,38,39]. Moreover, in recent years, GBS CC103 has gradually replaced the CC67 and become the dominant strain causing bovine mastitis [6]. At the same time, GBS CC103 and ST485 showed increased propagation in humans, especially in China [8, 11]. The phylogenetic analysis on GBS strains showed that GBS CC103 strains belonged to different evolutionary branches as compared to other GBS CCs, and they were clustered by themselves. Also, the human and bovine CC103 strains showed very close evolutionary relationships. These results indicate that the CC103 strains have unique evolutionary characteristics. CC103 generated multiple ST related evolutionary branches when propagated in human and bovine, including ST103, ST930, ST651, ST862 and ST485. ST103 strains can infect both humans and bovine, and they did not form independent branches in evolution, suggesting that ST103 strains did not evolve on a single host, and cross infection might happen between human and bovine. Among ST930, ST651, ST862 and ST485 strains, ST485 had the strongest infection ability [11]. We also found that ST485 could specifically infect human, which is consistent with the previous reports [8, 9, 11, 40]. Our results demonstrate a new relationship between human and bovine GBS stains, and indicate that GBS might form new, specific and highly pathogenic strain through propagating in human and bovine. However, the detailed mechanisms about how GBS propagates across human and bovine needs further investigation.

In order to study how the highly invasive ST485 appeared, we did genomic sequence alignment between GBS strains. Previous studies have found that genome recombination was a major driver for GBS genetic diversity [41, 42]. However, from the 202 invasive ST1 strains, only eight showed significant genome recombination; the remaining 194 strains only differed by 97 SNPs on average [43]. Our results also confirmed that there were a lot more mutations between different CCs, but less mutations within the same CC or ST. Compared to other CC103 strains, the genomic mutations of ST485 were enriched in several regions, and the affected genes were mainly coded for phage related proteins, hypothetical proteins with unknown function, and gene recombination related enzymes. The coding sequence alignment within the same STs showed that some of these mutations were random, and some were conserved, indicating that these mutations might affect GBS evolution, but they are not the key factors for the formation of genetic bottlenecks or the emergence and prevalence of new strains (such as ST485). For example, the key to the formation and prevalence of GBS CC17 is the tetracycline resistance gene, and the evolutionary bottleneck is the use of tetracycline, which led to a global replacement of GBS population [42]. According to the genomic coding sequence alignment and functional gene analysis, we found that Lac.2 and the cadmium resistance related genes might play a key role in the emergence and prevalence of ST485.

Carbohydrates are the most common energy source for cell growth, and thereby the carbohydrate metabolism plays an important role in the survival of prokaryotes [44, 45]. In nature, lactose is only found in mammalian milk, and is almost the only carbohydrate source in milk [46]. Therefore, Lac.2, an important gene structure for lactose metabolism, is essential for the bacteria that survive in milk. It has been shown that all the bovine GBS strains carry Lac.2 gene structure, which is necessary for lactose fermentation, while only a few human GBS strains have Lac.2 [35, 47]. The genomes of bovine and human ST1 strains have > 99% similarity, and the most significant difference was the presence of Lac.2 gene, indicating that a single genetic event can cause large phenotypic changes, such as the shift in host adaptability [43]. Our results showed that although both ST103 and ST485 strains belonged to CC103, most of the ST103 strains (both human and bovine) had Lac.2, but all ST485 strains did not have Lac.2, suggesting that the recombination of Lac.2 gene structure played a key role in GBS host adaptation.

We also noticed that the Lac.2 of GBS had two types, Lac.2–1 and Lac.2–2, and Lac.2–2 was mainly present in bovine GBS. However, the human ST103 strains GBS85147 and BSU451 had Lac.2–2, which were closely related to the Lac.2 in bovine strains. In addition, GBS85147 and BSU451 were isolated from human pharynx and respiratory tract, respectively (Additional file 1). Thus, we hypothesized that these two cases of GBS infection are possibly due to milk drinking, which leads to the propagation of bovine ST103 strains in human. Since human GBS were initially isolated from nose or throat, some studies also speculated that drinking milk might cause the infection of bovine GBS strains on human [48,49,50]. After pasteurized milk became commonly available, the chance of human getting infected by bovine GBS greatly reduced. However, the GBS detection rate in bulk tank milk (BTM) from dairy herds was still 90% or more [51]. Currently, there is no evidence confirming that cows are the storage pools for human epidemic GBS strains. Our results suggested that the GBS ST485 strains, which were suitable for propagation in human, might have evolved from GBS ST103, and the increasing prevalence of CC103 strains in human and bovine could be a potential threat for public health.

Whole-genome analysis of 150 bovine GBS isolates revealed that GBS CC61 had replaced the previous GBS population in Portugal, and all of the CC61 isolates acquired the mutations within an iron/manganese transporter, underlining a key adaptive strategy for colonizing in bovine host [52]. Compared to other CC103 strains, all the human ST485 strains had cadmium resistance related gene structure CadDX, which was closely related to bovine CC67 strain, but far from other human GBS strains, indicating that ST485 acquired CadDX gene structure through gene transfer. The transfer of cadmium resistance genes is often associated with the transfer of antibiotic resistance genes, which raises the possibility that antibiotic resistance and heavy-metal resistance were co-selected via the mobile genetic elements [53,54,55]. Cadmium pollution occurs in many countries, such as Japan, Thailand and China [56, 57]. The cadmium in soil and water are absorbed by crops, and accumulates in human body via food or drink [58, 59]. The estimated Cd in-take of each person is 30 mg per day [60, 61]. Cadmium enters bacterial cell via the transport systems that are normally used for essential divalent cations. By binding to sulfhydryl groups of essential proteins, cadmium can inhibit cell respiration and cause bacteria death [62]. With the accumulation of cadmium in environment and human bodies, GBS strains with cadmium resistance are easier to survive and spread. Therefore, we reasoned that acquiring the CadDX gene structure through horizontal gene transfer was another key factor for the emergence and propagation of ST485 strains in human.


Our results suggested that: 1. GBS CC103 strains could spread across human and bovine; 2. GBS ST485 evolved from ST103; 3. The recombination of Lac.2 and CadDX gene structures played an important role in the formation of highly pathogenic GBS ST485 in China. Based on these results, the cow farm workers should be cautious about GBS infection. Moreover, further studies need to be focused on the associations between GBS strains and hosts, which will be important for the prevention and control of GBS disease.



Bulk tank milk


Clonal complex


Clustered regularly interspaced short palindromic repeats


Group B Streptococcus


the Maximum likelihood


Multi-locus sequence typing


Sequence type


  1. 1.

    Jones N, Oliver KA, Barry J, Harding RM, Bisharat N, Spratt BG, et al. Enhanced invasiveness of bovine-derived neonatal sequence type 17 group B streptococcus is independent of capsular serotype. Clin Infect Dis. 2006;42:915–24.

    CAS  Article  Google Scholar 

  2. 2.

    Larppanichpoonphol P, Watanakunakorn C. Group B streptococcal bacteremia in nonpregnant adults at a community teaching hospital. South Med J. 2001;94:1206–11.

    CAS  Article  Google Scholar 

  3. 3.

    Jones N, Bohnsack JF, Takahashi S, Oliver KA, Chan MS, Kunst F, et al. Multilocus sequence typing system for group B streptococcus. J Clin Microbiol. 2003;41:2530–6.

    CAS  Article  Google Scholar 

  4. 4.

    Sorensen UB, Poulsen K, Ghezzo C, Margarit I, Kilian M. Emergence and global dissemination of host-specific Streptococcus agalactiae clones. MBio. 2010;1:e00178–10.

    Article  Google Scholar 

  5. 5.

    Zadoks RN, Middleton JR, McDougall S, Katholm J, Schukken YH. Molecular epidemiology of mastitis pathogens of dairy cattle and comparative relevance to humans. J Mammary Gland Biol Neoplasia. 2011;16:357–72.

    Article  Google Scholar 

  6. 6.

    Yang Y, Liu Y, Ding Y, Yi L, Ma Z, Fan H, et al. Molecular characterization of Streptococcus agalactiae isolated from bovine mastitis in eastern China. PLoS One. 2013;8:e67755.

    CAS  Article  Google Scholar 

  7. 7.

    Lu B, Wang D, Zhou H, Zhu F, Li D, Zhang S, et al. Distribution of pilus islands and alpha-like protein genes of group B Streptococcus colonized in pregnant women in Beijing, China. Eur J Clin Microbiol Infect Dis. 2015;34:1173–9.

    CAS  Article  Google Scholar 

  8. 8.

    Wang P, Tong JJ, Ma XH, Song FL, Fan L, Guo CM, et al. Serotypes, antibiotic susceptibilities, and multi-locus sequence type profiles of Streptococcus agalactiae isolates circulating in Beijing, China. PLoS One. 2015;e0120035:10.

    Google Scholar 

  9. 9.

    Jiang H, Chen M, Li T, Liu H, Gong Y, Li M. Molecular characterization of Streptococcus agalactiae causing community- and hospital-acquired infections in Shanghai China. Front Microbiol. 2016;7:1308.

    PubMed  PubMed Central  Google Scholar 

  10. 10.

    Guo D, Cao X, Li S, Ou Q, Lin D, Yao Z, et al. Neonatal colonization of group B Streptococcus in China: prevalence, antimicrobial resistance, serotypes, and molecular characterization. Am J Infect Control. 2018.

  11. 11.

    Li L, Wang R, Huang Y, Huang T, Luo F, Huang W, et al. High incidence of pathogenic Streptococcus agalactiae ST485 strain in pregnant/puerperal women and isolation of hyper-virulent human CC67 strain. Front Microbiol. 2018;9:1–14.

    Article  Google Scholar 

  12. 12.

    Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJ, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.

    CAS  Article  Google Scholar 

  13. 13.

    Li L, Stoeckert CJ Jr, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.

    CAS  Article  Google Scholar 

  14. 14.

    Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.

    CAS  Article  Google Scholar 

  15. 15.

    Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27:1164–5.

    CAS  Article  Google Scholar 

  16. 16.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.

    CAS  Article  Google Scholar 

  17. 17.

    Grissa I, Vergnaud G, Pourcel C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007;35:W52–7.

    Article  Google Scholar 

  18. 18.

    Bland C, Ramsey TL, Sabree F, Lowe M, Brown K, Kyrpides NC, et al. CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinform. 2007;8:209.

    Article  Google Scholar 

  19. 19.

    Lier C, Baticle E, Horvath P, Haguenoer E, Valentin AS, Glaser P, et al. Analysis of the type II-A CRISPR-Cas system of Streptococcus agalactiae reveals distinctive features according to genetic lineages. Front Genet. 2015;6:214.

    Article  Google Scholar 

  20. 20.

    Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75.

    Article  Google Scholar 

  21. 21.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    CAS  Article  Google Scholar 

  22. 22.

    Horvath P, Barrangou R. CRISPR/Cas, the immune system of bacteria and archaea. Science. 2010;327:167–70.

    CAS  Article  Google Scholar 

  23. 23.

    Bhaya D, Davison M, Barrangou R. CRISPR-Cas systems in bacteria and archaea: versatile small RNAs for adaptive defense and regulation. Annu Rev Genet. 2011;45:273–97.

    CAS  Article  Google Scholar 

  24. 24.

    Horvath P, Romero DA, Coute-Monvoisin AC, Richards M, Deveau H, Moineau S, et al. Diversity, activity, and evolution of CRISPR loci in Streptococcus thermophilus. J Bacteriol. 2008;190:1401–12.

    CAS  Article  Google Scholar 

  25. 25.

    Lopez-Sanchez MJ, Sauvage E, Da Cunha V, Clermont D, Ratsima Hariniaina E, Gonzalez-Zorn B, et al. The highly dynamic CRISPR1 system of Streptococcus agalactiae controls the diversity of its mobilome. Mol Microbiol. 2012;85:1057–71.

    CAS  Article  Google Scholar 

  26. 26.

    Finch LA, Martin DR. Human and bovine group B streptococci: two distinct populations. J Appl Bacteriol. 1984;57:273–8.

    CAS  Article  Google Scholar 

  27. 27.

    Pattison IH, Matthews PR, Howell DG. The type classification of group-B streptococci, with special reference to bovine strains apparently lacking in type polysaccharide. J Pathol Bacteriol. 1955;69:51–60.

    CAS  Article  Google Scholar 

  28. 28.

    Wanger AR, Dunny GM. Identification of a Streptococcus agalactiae protein antigen associated with bovine mastitis isolates. Infect Immun. 1987;55:1170–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Wibawan IW, Lammler C. Properties of group B streptococci with protein surface antigens X and R. J Clin Microbiol. 1990;28:2834–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Bisharat N, Crook DW, Leigh J, Harding RM, Ward PN, Coffey TJ, et al. Hyperinvasive neonatal group B streptococcus has arisen from a bovine ancestor. J Clin Microbiol. 2004;42:2161–7.

    Article  Google Scholar 

  31. 31.

    Bohnsack JF, Whiting AA, Martinez G, Jones N, Adderson EE, Detrick S, et al. Serotype III Streptococcus agalactiae from bovine milk and human neonatal infections. Emerg Infect Dis. 2004;10:1412–9.

    Article  Google Scholar 

  32. 32.

    Dogan B, Schukken YH, Santisteban C, Boor KJ. Distribution of serotypes and antimicrobial resistance genes among Streptococcus agalactiae isolates from bovine and human hosts. J Clin Microbiol. 2005;43:5899–906.

    CAS  Article  Google Scholar 

  33. 33.

    Sukhnanand S, Dogan B, Ayodele MO, Zadoks RN, Craver MP, Dumas NB, et al. Molecular subtyping and characterization of bovine and human Streptococcus agalactiae isolates. J Clin Microbiol. 2005;43:1177–86.

    CAS  Article  Google Scholar 

  34. 34.

    Emaneini M, Khoramian B, Jabalameli F, Abani S, Dabiri H, Beigverdi R. Comparison of virulence factors and capsular types of Streptococcus agalactiae isolated from human and bovine infections. Microb Pathog. 2016;91:1–4.

    CAS  Article  Google Scholar 

  35. 35.

    Richards VP, Lang P, Bitar PD, Lefebure T, Schukken YH, Zadoks RN, et al. Comparative genomics and the role of lateral gene transfer in the evolution of bovine adapted Streptococcus agalactiae. Infect Genet Evol. 2011;11:1263–75.

    CAS  Article  Google Scholar 

  36. 36.

    Springman AC, Lacher DW, Wu G, Milton N, Whittam TS, Davies HD, et al. Selection, recombination, and virulence gene diversity among group B streptococcal genotypes. J Bacteriol. 2009;191:5419–27.

    CAS  Article  Google Scholar 

  37. 37.

    Yildirim AO, Lammler C, Weiss R, Kopp P. Pheno- and genotypic properties of streptococci of serological group B of canine and feline origin. FEMS Microbiol Lett. 2002;212:187–92.

    CAS  Article  Google Scholar 

  38. 38.

    Godoy DT, Carvalho-Castro GA, Leal CA, Pereira UP, Leite RC, Figueiredo HC. Genetic diversity and new genotyping scheme for fish pathogenic Streptococcus agalactiae. Lett Appl Microbiol. 2013;57:476–83.

    CAS  Article  Google Scholar 

  39. 39.

    Brochet M, Couve E, Zouine M, Vallaeys T, Rusniok C, Lamy MC, et al. Genomic diversity and evolution within the species Streptococcus agalactiae. Microbes Infect. 2006;8:1227–43.

    CAS  Article  Google Scholar 

  40. 40.

    Lu B, Chen X, Wang J, Wang D, Zeng J, Li Y, et al. Molecular characteristics and antimicrobial resistance in invasive and noninvasive group B Streptococcus between 2008 and 2015 in China. Diagn Microbiol Infect Dis. 2016;86:351–7.

    CAS  Article  Google Scholar 

  41. 41.

    Brochet M, Rusniok C, Couve E, Dramsi S, Poyart C, Trieu-Cuot P, et al. Shaping a bacterial genome by large chromosomal replacements, the evolutionary history of Streptococcus agalactiae. Proc Natl Acad Sci U S A. 2008;105:15961–6.

    CAS  Article  Google Scholar 

  42. 42.

    Da Cunha V, Davies MR, Douarre PE, Rosinski-Chupin I, Margarit I, Spinali S, et al. Streptococcus agalactiae clones infecting humans were selected and fixed through the extensive use of tetracycline. Nat Commun. 2014;5:4544.

    CAS  Article  Google Scholar 

  43. 43.

    Flores AR, Galloway-Pena J, Sahasrabhojane P, Saldana M, Yao H, Su X, et al. Sequence type 1 group B Streptococcus, an emerging cause of invasive disease in adults, evolves by small genetic changes. Proc Natl Acad Sci U S A. 2015;112:6431–6.

    CAS  Article  Google Scholar 

  44. 44.

    Stulke J, Hillen W. Coupling physiology and gene regulation in bacteria: the phosphotransferase sugar uptake system delivers the signals. Naturwissenschaften. 1998;85:583–92.

    CAS  Article  Google Scholar 

  45. 45.

    Titgemeyer F, Hillen W. Global control of sugar metabolism: a gram-positive solution. Antonie Van Leeuwenhoek. 2002;82:59–71.

    CAS  Article  Google Scholar 

  46. 46.

    Lomer MC, Parkes GC, Sanderson JD. Review article: lactose intolerance in clinical practice--myths and realities. Aliment Pharmacol Ther. 2008;27:93–103.

    CAS  Article  Google Scholar 

  47. 47.

    Richards VP, Choi SC, Pavinski Bitar PD, Gurjar AA, Stanhope MJ. Transcriptomic and genomic evidence for Streptococcus agalactiae adaptation to the bovine environment. BMC Genomics. 2013;14:920.

    Article  Google Scholar 

  48. 48.

    Tongs MS. Hemolytic streptococci in the nose and throat with special reference to their occurrence after tonsillectomy. J Am Med Assoc. 1919;73:1050–3.

    Article  Google Scholar 

  49. 49.

    Meleney FL. Seasonal incidence of hemolytic streptococcus in the nose and throat. J Am Med Assoc. 1927;88:1392–4.

    Article  Google Scholar 

  50. 50.

    Fry RM, Eng MRCS. Fatal infections by hemolytic Streptococcus group B. Lancet. 1938;231:199–201.

    Article  Google Scholar 

  51. 51.

    Bi Y, Wang YJ, Qin Y, Guix Vallverdu R, Maldonado Garcia J, Sun W, et al. Prevalence of bovine mastitis pathogens in bulk tank Milk in China. PLoS One. 2016;11:e0155621.

    Article  Google Scholar 

  52. 52.

    Almeida A, Alves-Barroco C, Sauvage E, Bexiga R, Albuquerque P, Tavares F, et al. Persistence of a dominant bovine lineage of group B Streptococcus reveals genomic signatures of host adaptation. Environ Microbiol. 2016;18:4216–29.

    CAS  Article  Google Scholar 

  53. 53.

    Rojo-Bezares B, Azcona-Gutierrez JM, Martin C, Jareno MS, Torres C, Saenz Y. Streptococcus agalactiae from pregnant women: antibiotic and heavy-metal resistance mechanisms and molecular typing. Epidemiol Infect. 2016;144:3205–14.

    CAS  Article  Google Scholar 

  54. 54.

    Gomez-Sanz E, Kadlec K, Fessler AT, Zarazaga M, Torres C, Schwarz S. Novel erm(T)-carrying multiresistance plasmids from porcine and human isolates of methicillin-resistant Staphylococcus aureus ST398 that also harbor cadmium and copper resistance determinants. Antimicrob Agents Chemother. 2013;57:3275–82.

    CAS  Article  Google Scholar 

  55. 55.

    Silveira E, Freitas AR, Antunes P, Barros M, Campos J, Coque TM, et al. Co-transfer of resistance to high concentrations of copper and first-line antibiotics among enterococcus from different origins (humans, animals, the environment and foods) and clonal lineages. J Antimicrob Chemother. 2014;69:899–906.

    CAS  Article  Google Scholar 

  56. 56.

    Suwazono Y, Nogawa K, Uetani M, Miura K, Sakata K, Okayama A, et al. Application of hybrid approach for estimating the benchmark dose of urinary cadmium for adverse renal effects in the general population of Japan. J Appl Toxicol. 2011;31:89–93.

    CAS  Article  Google Scholar 

  57. 57.

    Nishijo M, Suwazono Y, Ruangyuttikarn W, Nambunmee K, Swaddiwudhipong W, Nogawa K, et al. Risk assessment for Thai population: benchmark dose of urinary and blood cadmium levels for renal effects by hybrid approach of inhabitants living in polluted and non-polluted areas in Thailand. BMC Public Health. 2014;14:702.

    Article  Google Scholar 

  58. 58.

    Seebaugh DR, Goto D, Wallace WG. Bioenhancement of cadmium transfer along a multi-level food chain. Mar Environ Res. 2005;59:473–91.

    CAS  Article  Google Scholar 

  59. 59.

    Abdu N, Agbenin JO, Buerkert A. Phytoavailability, human risk assessment and transfer characteristics of cadmium and zinc contamination from urban gardens in Kano, Nigeria. J Sci Food Agric. 2011;91:2722–30.

    CAS  Article  Google Scholar 

  60. 60.

    Nriagu JO, Pacyna JM. Quantitative assessment of worldwide contamination of air, water and soils by trace metals. Nature. 1988;333:134–9.

    CAS  Article  Google Scholar 

  61. 61.

    Joseph P. Mechanisms of cadmium carcinogenesis. Toxicol Appl Pharmacol. 2009;238:272–9.

    CAS  Article  Google Scholar 

  62. 62.

    Vallee BL, Ulmer DD. Biochemical effects of mercury, cadmium, and lead. Annu Rev Biochem. 1972;41:91–128.

    CAS  Article  Google Scholar 

Download references


This work was supported by National Natural Science Foundation of China (31460695) and Guangxi innovation-driven development special funds (Grant no. AA17204081–3) to MC, Guangxi Natural Science Foundation (2016GXNSFDA380020) to LL, the funds of Guangxi Key Laboratory for Aquatic Genetic Breeding and Healthy Aquaculture (2016–2018). The funding body did not exert influence on the design of the study, and collection, analysis, and interpretation of data or in writing of the manuscript.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the GenBank database,, and all genome accession numbers are listed in Additional file 1.

Author information




RW and MC designed the study. RW and WH performed genome sequencing and annotation; RW, LL, and TH performed phylogenetic analysis and CRISPRs analysis; YH, XY, and AL performed whole genome sequences and functional genes comparison; RW, LL, TH, and MC wrote the manuscript. All authors read, revised extensively, and gave final approval of the manuscript.

Corresponding author

Correspondence to Ming Chen.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Sequenced strains and available genomes used in this study. (XLSX 16 kb)

Additional file 2:

Inventory and distribution of CRISPR1 the terminal repeat (RT) and the terminal spacer (ST) sequences among S. agalactiae CC or ST. (XLSX 9 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, R., Li, L., Huang, T. et al. Phylogenetic, comparative genomic and structural analyses of human Streptococcus agalactiae ST485 in China. BMC Genomics 19, 716 (2018).

Download citation


  • Group B Streptococcus
  • CC103
  • ST485
  • Lac.2
  • CadDX
  • Cross infection