Reductive evolution in Streptococcus agalactiaeand the emergence of a host adapted lineage
© Rosinski-Chupin et al.; licensee BioMed Central Ltd. 2013
Received: 31 December 2012
Accepted: 1 April 2013
Published: 15 April 2013
During host specialization, inactivation of genes whose function is no more required is favored by changes in selective constraints and evolutionary bottlenecks. The Gram positive bacteria Streptococcus agalactiae (also called GBS), responsible for septicemia and meningitis in neonates also emerged during the seventies as a cause of severe epidemics in fish farms. To decipher the genetic basis for the emergence of these highly virulent GBS strains and of their adaptation to fish, we have analyzed the genomic sequence of seven strains isolated from fish and other poikilotherms.
Comparative analysis shows that the two groups of GBS strains responsible for fish epidemic diseases are only distantly related. While strains belonging to the clonal complex 7 cannot be distinguished from their human CC7 counterparts according to their gene content, strains belonging to the ST260-261 types probably diverged a long time ago. In this lineage, specialization to the fish host was correlated with a massive gene inactivation and broad changes in gene expression. We took advantage of the low level of sequence divergence between GBS strains and of the emergence of sublineages to reconstruct the different steps involved in this process. Non-homologous recombination was found to have played a major role in the genome erosion.
Our results show that the early phase of genome reduction during host specialization mostly involves accumulation of small and likely reversible indels, followed by a second evolutionary step marked by a higher frequency of large deletions.
KeywordsStreptococcus agalactiae Host-adaptation Non-homologous recombination Gene inactivation Virulence
Comparative genomics of strains belonging to species with a large spectrum of hosts, such as Staphylococcus aureus, have highlighted two main evolutionary trends linked to the adaptation to a new host: acquisition of new functions through lateral gene transfer facilitating colonization of new niches and gene loss associated with the host specialization [1–3]. These trends were also recognized when considering the emergence of highly virulent host-specialized pathogens from bacterial species with a broader host range. For instance massive gene losses were associated to the emergence of human or equine pathogens such as Salmonella enterica sv Typhi and Paratyphi, Bordetella pertussis and parapertussishom and Burkholderia mallei that respectively derive from Salmonella Typhimurium, Bordetella bronchiseptica and Burkholderia pseudomallei[4–9]. During these transitions, it was postulated that gene inactivation and deletions were probably favored by genetic drift and evolutionary bottlenecks. In line with this model the host-specialized pathogens often showed a much higher number of insertions sequences (IS) than their parental strains and IS expansion was proposed to be largely responsible for gene deletions and genome rearrangements observed in these species .
Streptococcus agalactiae also referred to as Group B streptococcus (GBS) is a Gram-positive bacterium that has emerged as a leading cause of neonatal infections during the sixties and represents an increasing cause of infections in the elderly and in adults with underlying diseases [11–13]. As a commensal it colonizes the digestive and genitourinary tracts of up to 30% of the human adult population . However, S. agalactiae was initially described as an animal pathogen causing mastitis in ruminant . Since the 70’s, S. agalactiae was found to be responsible for epidemic events of invasive diseases in fish farms, leading to a mortality of up to 30% [16–18]. Cases of infection were also reported for other aquatic poikilotherms such as frogs  and aquatic mammals such as dolphins . How GBS is able to adapt to its different hosts remains poorly understood. The genetic diversity of GBS populations has been studied using different methods including multilocus sequence typing (MLST) , which led to the recognition of different clonal complexes (CC). Some of these clonal complexes display host preference. For instance, CC67 is essentially associated with the bovine host and the hypervirulent sequence type (ST) 17 strains are mainly isolated from humans. However incidentally strains belonging to human-associated clonal complexes are also isolated from bovines  suggesting that relationships between clonal complexes and host specificity are not so strict. Further analysis of the complete genome sequences of eight isolates of human origin and one of bovine origin has highlighted the composite organization of S. agalactiae genomes with a conserved backbone (representing the core genome of the species) and a dispensable genome composed of genomic islands that are highly variable between the different strains [23–26].
The ST261 strain 2-22 (or ATCC 51487) was initially isolated as responsible for several epidemics in fish farms in Israel [27, 28]. This strain, which showed a restricted metabolic pattern, thermosensitivity and lack of ß hemolytic and CAMP activities was first classified as a different species, Streptococcus difficile, but proved later to be a genuine serotype Ib S. agalactiae strain . S. agalactiae strains were repeatedly isolated from fish infections and found to cluster into two main groups [20, 29]. The first group corresponds to strains belonging to the clonal complex 7, also displaying strains isolated from human and bovine hosts. The other strains share two or more common MLST alleles with strain 2-22 and are classified in ST246, 257, 259, 260, 552 and 553. As ST261 strains, these STs were until now never isolated from humans . In addition, a third group of strains belonging to clonal complex 283 has recently been described in fish and humans . The draft genome sequences of one ST260 strain, strain STIR-CD-17, and two ST7 strains, strain ZQ0910 and GD201008-001, respectively isolated from disease outbreaks affecting farmed tilapia in Honduras and Nile Tilapia in China were recently published [30–32].
To decipher the phylogenomic relationships between S. agalactiae strains isolated from fish or other poikilotherm animals and strains isolated from human or bovine, we analyzed the genomic sequence of seven isolates belonging to ST260-261 and ST6-7, including the strain 2-22. Comparative analysis confirmed that the two groups of GBS strains responsible for fish epidemic diseases are distantly related. We found that adaptation to fish does not involve any specific function compared to human CC7 isolates. Conversely, specialization to the fish host of the ST260-261 strains was associated with massive gene inactivation and deep remodeling of metabolic and regulatory networks that we also characterized at the transcriptome level. This genome reduction likely occurred through RecA independent recombination.
Results and discussion
Fish ST7 strains are closely related to human strains
Characteristics of the different strains used in the study
Genome size (kbp)
tilapia (Oreochromis sp.)
Hybrid striped bass
Gene content was similar between strains CF01173 and A909 except for 13 genes that were disrupted in CF01173 and five in A909 (Additional file 2: Table S3). Seven genes were specifically disrupted in strain SS1014. In addition, we found that the three ST7 strains CF01173, ZQ0910 and GD201008-001 isolated from fish shared one short genomic island that was absent from other GBS and probably resulted from lateral gene transfer. This genomic island encodes proteins 70-95% identical with proteins of Streptococcus anginosus, including a protein with a LPXTG-motif cell-wall anchor domain (GBS1173_1788). Sharing of this island in addition to the low number of SNPs between strains CF01173, ZQ0910 and GD201008-001 suggests that the three strains have a common recent ancestor that propagated worldwide. In contrast, emergence of the ST6 strain SS1014 can be considered as an independent event.
Overall, the low level of polymorphism with human strain indicates that CC7 strains infecting fish recently diverged from strains isolated in humans and bovines. Since at least two independent events of emergence were observed, this suggests that CC7 strains might be more amenable to fish colonization/infection than other GBS clonal complexes isolated from humans or bovines. In agreement with this hypothesis, it was recently shown that a S. agalactiae ST7 strain isolated from human was able to cause disease in Nile tilapia .
ST260-261 strains form an independent lineage which underwent reductive evolution
Whole-genome sequence comparison of the five strains showed that they clustered into two distinct subgroups correlating with the MLST classification (Additional file 1: Table S2). The ST261 strain SS1218 isolated from frog in Louisiana differed from the strain 2-22 by only 30 SNPs. Strains 90-503 isolated in Louisiana in 1990 and 05-108A isolated in 2005 in Honduras, with 49 SNPs can be considered as variants of the same clone. Strain SS1219 isolated from frog in Taiwan diverged from the 2 former ST260 strains by 130 SNPs. On average, the ST260 strains showed 3,100 SNPs with ST261 strains. ST260-261 strains were also related to the ST552 strain Sa20-06 (3,400 and 1,700 SNPs respectively). In contrast, ST260-261 displayed 15,000 SNP with CC7 strains. Analysis of the SNP distribution along the genome sequences of ST260 and ST261 strains revealed a uniform pattern of 3 SNPs per kb (Figure 2c), suggesting that no recombination occurred in this lineage.
The low proportion of nucleotides that did not align with the genomes of human GBS strains suggested that ST260-261 strains encode only few specific functions. Indeed, only four genomic islands were characterized in strain 2-22, that represented 25 kb in total (in violet on Figure 2b). However these four regions essentially contained pseudogenes and one of them was absent from the ST260 strains. We also identified in ST260-261 strains eleven to twenty copies of ISSag1, an insertion sequence previously described in the genomes of human isolates and other streptococci  (Additional file 3: Table S4). In addition, the genomes of ST260-261 strains contained 10 regions categorized as genomic islands in human S. agalactiae genomes , but shared by most GBS characterized so far (shown for strain 2-22 in Figure 2b). They did not contain any Integrative and Conjugative Elements (ICE) or prophages and only two inactivated copies of integrases were identified. They also lacked two mobile genetic elements encoding important virulence loci: the cis-mobilizable element encoding the major surface antigen Alpha like protein/Rib  and the composite transposon coding for the C5A peptidase and the laminin binding protein . In total the difference in genomic island content between ST260-261 strains and GBS strains isolated from humans or bovines accounted for about 70-80% of the genome reduction (210-230 kb). An interesting exception was the genomic island 3.2 that was previously described only in strains A909 (ST7) and H36B (ST6) and therefore is also shared by CC7 strains isolated from fish. This conservation is in favor of a role in fish colonization. This GI encodes two phosphotransferase systems (PTS) for galactitol, sugar ABC transporters and genes for galactose utilization (GBS222_0398-0414 in strain 2-22) (Additional file 4: Figure S1). Sharing of this genomic island explained the higher number of aligned nucleotides between strains 2-22 and A909.
Genome annotation predicted 1547, 1568, 1560 and 1569 protein coding genes for strains 2-22, 90-503, SS1219 and 05-108A respectively, which is significantly less compared to 2096 in NEM316, 1990 in A909 and 2135 in 2603V/R. Furthermore, 190-220 pseudogenes were identified in each strain (Additional file 5: Table S5) compared to 27 to 41 in the human strains. This revealed a massive reduction of the functional genome during the time-course of adaptation to fish in the ST260-261 lineage.
Reductive evolution is an ongoing process in the ST260-261 lineage
In addition, 235 indels and 6363 SNPs relative to strain A909 sequence were common to the five ST260-261 strains but were not observed in A909 vs other GBS strains. Forty per cent of these indels were intragenic and led to 91 gene disruptions (Figure 4c). Finally, 34 genes disrupted by a small indel in one sublineage were deleted in the other, suggesting that gene disruption could be a preliminary step to gene loss and genome reduction by secondary longer deletions. As a consequence some loci were observed under different states of decay in the two sublineages. This is the case for instance for the cyl locus, the CRISPR2 locus and the pil2 pilus locus (Additional file 4: Figure S2). Furthermore, the average size of the indels was found to increase following the divergence of the two sublineages (Figure 4c) compared to the common branch. Therefore, while the process of reductive evolution was already evident in the ancestor of the ST260-261 strains, it was even more pronounced after divergence of the two sublineages.
To better evaluate the specificity of the evolutionary process occurring in ST260-261 strains, we compared the core genome sequence of strains A909 and NEM316, two strains of human origin differing by approximately 11,000 SNP. We detected 314 indels, which represents half the number of indels in the 2-22 vs A909 comparison. These indels were essentially short and as much as 80% of them occurred in intergenic regions while less than 4% led to gene inactivations. Therefore, gain or loss of nucleotides was also frequent in other S. agalactiae strains, but most indels in coding sequences were probably eliminated by purifying selection. Fitting with this hypothesis, we observed one indel per 5-10 SNPs during the recent evolution of ST7 strains where purifying selection is probably not effective. A similar proportion of indels vs SNP was also reported for two strains of S. enterica sv Paratyphi .
From this analysis we propose that accumulation of small and easily reversible indels is common in S. agalactiae, mostly occurring as a consequence of RecA independent recombination. Most of these indels are probably removed from the population due to purifying selection during long-term evolution. Alternatively they may be compensated by a second loss or gain of nucleotides. In the ST260-261 lineages host specialization led to the relaxation of the negative selection on genes that became dispensable, allowing the accumulation of larger deletions.
Functional adaptation to a host restricted way of living
Bacterial-host interactions and virulence factors
Potential virulence factors in the ST260-261 strains
Mechanism of pathogenicity
Capsule synthesis enzymes
Penicillin-binding protein 1
Resistance to host antimicrobial peptides
fibronectin/fibrinogen binding protein
Adherence and invasion
Fibrinogen-binding protein A
Serin-rich repeat protein
Surface immunogenic protein
Surface protein, BibA family
Cell-wall anchored proteins of unknown functions
serine-rich surface protein
surface protein, Amidase family
Among the other proteins that are virulence factors in human GBS isolates, the five strains lack the C5A peptidase, the laminin binding protein and components of the pilus. The cyl locus has also been inactivated in the two sublineages (Additional file 4: Figure S2) in agreement with the absence of detectable hemolytic activity on blood agar plates. Another major test for S. agalactiae identification is the detection of CAMP activity [63, 64]. While both ST260 and 261 strains were reported to be negative for the CAMP factor reaction  we found that the gene encoding the CAMP factor was disrupted only in ST261 strains. We experimentally confirmed that strain SS1219 was negative for the CAMP test detection, probably because of a lower level of gene expression (see below).
Finally, the five ST260-261 strains express a fibrinogen/fibronectin binding protein of the PavA family and a hyaluronidase that were proposed to be virulence factors in GBS or other Streptococcus species [62, 65, 66] and could also have a role in fish infection.
Altogether, approximately 60% of the genes for proteins involved in pathogenesis and considered as important virulence genes in human stains  were affected by genome reduction. In this context the capsule could be a major virulence factor in the fish host, as also observed for Streptococcus iniae, another streptococcus species pathogenic for fish [68, 69].
Analysis of the missing functions revealed a profound remodeling of the metabolism of strains of this lineage. Numerous transport systems for carbon sources (ABC transporters and phosphoenolpyruvate/carbohydrate phospho-transferase systems) and enzymes for degradation of polysaccharides (amylase, extracellular pullulanase, enzyme for degradation of arbutin) were missing or inactivated, reflecting a reduced capacity to utilize diverse carbon sources. In addition ST260-261 strains seem unable to utilize glycerol and glycerol phosphate as genes encoding the glycerol kinase, the glycerol dehydrogenase and the glycerol-phosphate permease are missing or mutated. Fermentative pathways utilizing pyruvate/acetate conversion are also altered, as the phosphotransacetylase gene is missing. In mammals, GBS is primarily considered as a commensal of the digestive tract, an environment rich in diverse C-sources. In contrast, reduction in the catabolic capacities observed in ST260-261 strains is in favor of a transition to an obligate pathogen style-of-life.
Approximately 20% of the genes associated with transport systems were missing or inactivated. The targeted functions were mostly the import of nutrients, in particular of C-sources and the transport of inorganic and metal ions. The Na+/H + antiporter, the K + uptake permease were inactivated in both sublineages. Therefore ST260-261 strains may be affected in ionic exchange and would have a reduced capacity to maintain their homeostasis in the face of a changing external environment.
Globally, gene disruption and deletion events affected 19 out of the 93 transcriptional regulators predicted in human strains (20%). As much as 13 out of the 21 two-component systems (TCS) (60%) found in GBS were inactivated, either in one (six TCS) or in both (seven TCS) sublineages (Additional file 6: Table S6). Interestingly, analysis of the genome sequence of strain 2-22 revealed that the Rgf TCS, which is involved in the control of virulence in the human ST17 strains , was associated with two putative bacteriocins with a double glycine leader peptide and with a bacteriocin export transporter (Additional file 4: Figure S3). This suggests that this TCS originates from a bacteriocin operon that has been partially deleted in human strains. Altogether our observations show that S. agalactiae strains adapted to fish may have a reduced capacity to respond to environmental changes compared to human strains and only eight TCS, including the two major systems, CiaRH and CovRS, may be sufficient to allow GBS adaptation to the different environments encountered in fish. Interestingly, both CiaRH and CovRS were also involved in the regulation of virulence genes in GBS human strains [71, 72]. The higher virulence of ST260-261 strains might also be due to the deregulated expression of some virulence genes, as observed in the transition from local to systemic infections in Group A Streptococci [73, 74].
Adaptation to fish is associated with broad changes in gene expression
Both the number of genes for which expression was modified and the amplitude of these modifications were even larger in ST260-261 isolates (Figure 7a; Additional file 7: Table S7). For instance two genes encoding putative glyoxalases are expressed at 6-45 higher levels in ST260-261 strains than in strain A909, suggesting higher needs in methylglyoxal detoxification. Genes encoding the enzymes for the synthesis of arginine (arginosuccinate synthase and arginosuccinate lyase) were 10-15 fold more expressed. Four functional categories were more affected by variations in gene expression in both ST260-261 sublineages: metabolism, transport, regulation and signal transduction (Figure 7b). Interestingly these functional categories were also particularly affected by gene erosion, suggesting a remodeling of metabolic networks.
In addition to these global tendencies, ST260 and ST261 strains also harbor sublineage-specific variations in gene expression that probably result from differential gene decay and might lead to specific adaptations. For instance, in strain SS1219, genes encoding phosphoglycerate kinase, phosphoglyceromutase and pyruvate kinase involved in glycolysis and neoglucogenesis are expressed at a lower level than in strain A909. In parallel, genes coding for the arginine deiminase, ornithine carbamoylkinase and carbamate kinase involved in the deiminase pathway are up-regulated, suggesting a shift towards arginine fermentation as the main pathway for energy production in ST260 strains . This up-regulation of the deiminase pathway is associated with a decreased expression of the negative regulator ArgR and the increased expression of the arginine/ornithine antiporter. In contrast, in the same growth conditions, upregulation of the genes for Hpr-kinase and CodY involved in catabolic repression and for lactate dehydrogenase suggests that strain 2-22 might generate energy essentially through glucose utilization and lactic fermentation. Iron import might be more efficient in ST260 than in ST261 strains since strain SS1219 expresses genes for iron ABC transporter and ferrichrome ABC transporter at higher levels than strain A909, while these genes are down-regulated or unchanged in strain 2-22. Both sublineages express lower amounts of ATP synthase, the enzyme responsible for ATP-proton motive force interconversion, than strain A909, indicating that they may be impaired in pH homeostasis in acidic conditions .
Since expression of pseudogenes is energetically costly and may generate deleterious products, we took advantage of the massive genome erosion in the ST260-261 lineage to look for evidence for negative selective pressure acting on pseudogene expression. Among the 184 pseudogenes represented by at least one probe on the array, only 22% were down-regulated compared to strain A909, whereas 64% did not present a significant change in mRNA levels and 14% were even up-regulated. Furthermore silencing of pseudogenes in the ST260-261 lineage did not seem to markedly increase with time. Indeed the proportion of down-regulated pseudogenes was similar among pseudogenes arising before (28%) and after the differentiation of sublineages (16% and 25%). Altogether this indicates that no strong evolutionary pressure acts to silence pseudogene expression.
Although modifications in gene expression generally affected different genes in ST260-261 and CC7 strains, 29 genes were found to vary in the same way in the 3 strains isolated from fish. Interestingly, these 29 genes included the operon for the capsule synthesis, the genes for zoocin and hyaluronidase, two response regulators of TCS systems (ciaR and relR) and three targets of the CiaRH TCS that were up-regulated. Sixteen genes were downregulated or inactivated in the three fish strains compared to A909, among which 11 were involved in energy metabolism or in carbohydrate uptake. Whether these regulations reflect a common mechanism of adaptation to fish environment remains to be established. However, the higher expression of the genes for the capsule synthesis might favor resistance to environmental conditions and to fish immune system, as it has also been reported for S. iniae. Similarly the high hyaluronidase expression may help S. agalactiae strains to break through fish tissues and be involved in virulence to fish .
Our results show that S. agalactiae strains leading to epidemic diseases in fish farms and cold blood animals belong to at least two distinct groups that differ by their strategies of host adaptation. CC7 strains have the potency to colonize and infect multiple hosts such as fish, human and cattle. From a genomic point of view, these CC7 fish strains are not distinguished from their human counterparts by any significant genomic island. However, contrasting with this genomic relatedness, large differences in gene expression were observed and could participate to the adaptation to the fish host. Conversely, our genome analysis indicates that strains of the ST260-261 complex diverged anciently from human and cattle strains and subsequently accumulated specific adaptations leading to the emergence of sublineages.
ST260-261 strains exhibit a striking pattern of genome reduction and we took advantage of the emergence of sublineages to reconstruct the different steps involved in this process. We found that accumulation of short indels can be observed all along the evolution of the GBS species, participating to strain-specific gene disruptions. Therefore even in the core genome of human GBS strains some genes are dispensable. Nevertheless the number of inactivated genes greatly increased during specialization of ST260-261 strains to fish. These gene inactivations mainly result from the ongoing accumulation of short indels, but a tendency to eliminate inactivated genes by deleting longer sequences is more noticeable in the sublineage specific branches. In contrast with what has previously been observed in the course of genomic reduction associated to host specialization or to intracellular symbiosis, deletion events are not correlated with an amplification of insertion sequences. Neither could complex genome rearrangements be noticed, suggesting that recombination between IS is not a general mechanism for genome reduction. Indeed our results point to non-homologous recombination as an alternative mechanism of genome reductive evolution.
Bacterial strains and growth conditions
S. agalactiae strains used in this study are described in Table 1. They were grown in Todd Hewitt medium (Difco) at 37°C except strains belonging to ST260-261 that were grown at 30°C.
Genome sequencing and assembly methods
To perform the complete sequence of strain 2-22 (ATCC 51487), a mix of capillary Sanger and 454 pyrosequencing (Roche, http://www.roche.com) was carried out until ~24-fold coverage. A single and a 3kb insert size mate-paired libraries, leading to around 23-fold coverage of 454 GSflx reads (11-fold coverage and 12-fold coverage respectively), were added to Sanger reads, which was derived from a 10 kb insert fragment size library. This library was constructed after mechanical shearing of genomic DNA and cloning of generated inserts into plasmid pCNS (pSU18-derived). Plasmid DNAs were purified and end-sequenced (3633 reads) by dye-terminator chemistry with ABI3730 sequencers (Applied Biosystems, Foster City, USA) allowing an approximately 1-fold supplementary coverage of the genome. The reads were assembled by Newbler (Roche) and validated via the Consed interface (http://www.phrap.org). A first finishing step was performed using primer walking of clones and polymerase chain reactions (PCRs) (92 and 119 additional reads respectively). Around 70-fold coverage of Illumina reads (36 bp) were mapped, using SOAP (http://soap.genomics.org.cn), for the polishing phase as previously described . Remaining gaps were then covered using PCRs with Platinum Hi-fidelity (InVitrogen) and primers specific to gap edges. The PCR fragments were sequenced by primer walking. Order and orientation of the sequences contiguous to ribosomal DNA was determined using long range PCR with PlatinumHI-fi (Invitrogen).
The genomes of the six other strains isolated from fish or frog were sequenced on the Illumina GAIIx with 36-cycle single-end chemistry and coverages of 100-300×. Reads were assembled using Velvet  to produce draft sequences. The number of contigs in each sequence is detailed in Additional file 8: Table S8. Contigs were reordered by Mauve 2.3.1 software  using the strain 2-22 genome as a reference for ST260-261 strains and the strain A909 genome for CC7 strains. When possible, overlapping ends of contigs were further assembled. Draft sequences were used to realign the reads using Bowtie , followed by alignment visualization with Tablet  to detect sequencing errors. Potential sequence ambiguities were only noticed for the ends of the contigs and therefore SNPs determined at ends of contigs were not considered as reliable.
Nucleotide sequence accession number
The S. agalactiae 2-22 complete genome sequence is available from DDBJ/GenBank/EMBL under accession number FO393392. Draft genome sequences of strains 90-503, SS1219, CF01173, SS1014, SS1218 and 05-108A are available under accession numbers CAPZ01000001-CAPZ01000082, CAQA01000001-CAQA01000070, CAQB01000001-CAQB01000131, CAQC01000001-CAQC01000099, CAUB01000001-CAUB01000070 and CATH01000001-CATH01000087 respectively.
Annotation of strain 2-22 was performed using the CAAT-box environment as previously described . Briefly, coding sequences (CDS) were defined by combining Genmark predictions with visual inspection of each open reading frame (ORF) for the presence of a start codon with an upstream ribosome-binding site and blastp similarity searches on Uniref 90, Trembl and Swissprot databases. Function predictions were based on blastp similarity searches and on the analysis of motifs using the PFAM databases . For the six other strains putative CDS were determined using Artemis . Predicted CDS were tested in reciprocal pair-wise blastp comparison against NEM316, A909 and 2-22 sequences. ORFs longer than 100 codons or encoding a protein > 90% identity by blastp with a protein annotated in NEM316, A909 or 2-22 were retained for further analyses. For genes with an orthologous sequence in other S. agalactiae genomes, assignment of the first codon and annotation were transferred from the orthologous genes. The remaining genes were manually annotated. Open reading frames whose lengths differed by more than 20% from orthologous sequences in NEM316 or A909 sequences were considered as putative pseudogenes and individually analyzed. Differences caused by variant start predictions were eliminated, as well as pseudogenes resulting from a mutation localized in the first or in the last 50 nucleotides of a contig, considered as more error-prone.
Whole genome sequence comparisons
Determination of insertion/deletion events in ST260-261 strains compared to strain A909 [GenBank: NC_007432] was performed by combining blastn and nucmer analyses . Each indel was individually inspected to differentiate possible misalignments from true indels and studied for the presence of repeat sequences at its extremities. We also aligned to A909 sequence the genome sequences of three S. agalactiae strains of human and bovine origins (NEM316, GenBank: [NC_004368]; 2603V/R [GenBank: AE009948] and FSL S3-O26/AEXT01000000 [GenBank: AEXT00000000] ). These alignments were used to produce tables of SNPs, indels and gaps relative to A909 genome sequence. Only indels that occurred internally to Velvet contigs in draft genome sequences and were repetitively observed for all strains of the same sequence type were taken into account. No significant difference in indel detection was noted using the whole-genome sequence of strain 2-22 or the Velvet-assembled contigs of the related ST261 strain SS1218. Indels relative to A909 core genome sequence were therefore classified into three categories related to strain phylogeny: i) detected in ST260-261 strains and in one or several other GBS strains, ii) observed for all strains belonging to the ST260-261 lineage but not detected in strains of bovine or human origins and iii) specific to one of the ST260-261 sublineages. For this last category, the allelic sequence identical to that of the human and bovine strains was considered to be the ancestral sequence, allowing to discriminate the mutated sequence and the nature of the mutation (insertion or deletion of nucleotides).
The genome sequence alignment was used to generate a table of SNPs from the S. agalactiae core genome. Regions that were recombined between strains A909 and SS1014 were removed from the analysis, as well as sequence gaps. In total 25,597 polymorphic positions in a 1.34 Mb core genome sequence were used to generate a consensus phylogenetic tree by the Maximum Likelihood method based on the Tamura-Nei model with the MEGA 5.04 software . A bootstrap consensus tree was inferred from 1000 replicates. The same tree was obtained using a Neighbor-Joining method and computation of evolutionary distances by the p-distance method.
The SNP distribution between two genome sequences was inferred from genome alignment using nucmer and SNP were counted on a 1000 nt window. Indels were considered as unique mutation events and therefore counted as one SNP. The goodness of fit to a uniform distribution was tested using a Chi square test on class intervals of 50000 or 100000 nt as indicated in the figure legends.
RNA preparation and transcriptome analyses
S. agalactiae strains were grown in TH medium at 30°C (strains 2-22 and SS1219) or 37°C (strain CF01173). For each transcriptome analysis, the reference strain A909 was grown at the same temperature as the tested strain. All bacterial cultures were harvested for RNA isolation at mid-exponential growth phase (OD 0.35-0.4). Total RNA was extracted as previously described . RNA was prepared from three independent cultures and each RNA sample was hybridized twice to the microarrays (dye swap). RNA was reverse-transcribed with Superscript indirect cDNA kit (Invitrogen) and labeled with Cy5 or Cy3 (Amersham Biosciences) according to the supplier’s instructions. The microarray contains 12889 45-60mer oligonucleotides designed on the predicted gene sequences of strains NEM316, A909 and COH1. The oligonucleotide design was carried out with the OligoArray server  (http://berry.engin.umich.edu/oligoarray/). The microarray was manufactured by Agilent Technologies. Only probes with 100% identity to A909, CF00173, 2-22 and SS1219 sequences were taken into account in the data analysis (ie 4401 probes on 1389 genes). Analysis was performed as described . The transcriptome data are MIAME compliant and have been submitted to the ArrayExpress database maintained at http://www.ebi.ac.uk/microarray-as/ae/ under the Accession Numbers E-MEXP-3828, E-MEXP-3830 and E-MEXP-3829.
Group B Streptococcus
Mobile genetic element
Integrative conjugative element
Multilocus sequence typing
- CAMP factor:
Christie Atkins Munch-Petersen factor
Single nucleotide polymorphism
Open reading frame.
We specially thank Drs J. Hawke, B. Beall and E. Duchaud for providing us the GBS strains used in the study. We thank Drs. E. Souche and P. Lechat for their help in bioinformatics analysis and Pr P. Trieu-Cuot and Dr. P.E. Douarre and R. Guérillot for helpful discussions. This work was supported by the French National Research Agency (grants ANR-08-GENM-027-001) and the LabEx project IBEID. The Institut Pasteur Genopole is a member of France Génomique (ANR10-IBNS-09-08).
- Ben Zakour NL, Guinane CM, Fitzgerald JR: Pathogenomics of the staphylococci: insights into niche adaptation and the emergence of new virulent strains. FEMS Microbiol Lett. 2008, 289 (1): 1-12. 10.1111/j.1574-6968.2008.01384.x.View ArticlePubMedGoogle Scholar
- Lowder BV, Guinane CM, Ben Zakour NL, Weinert LA, Conway-Morris A, Cartwright RA, Simpson AJ, Rambaut A, Nübel U, Fitzgerald JR: Recent human-to-poultry host jump, adaptation, and pandemic spread of Staphylococcus aureus. Proc Natl Acad Sci USA. 2009, 106 (46): 19545-19550. 10.1073/pnas.0909285106.PubMed CentralView ArticlePubMedGoogle Scholar
- Guinane CM, Penades JR, Fitzgerald JR: The role of horizontal gene transfer in Staphylococcus aureus host adaptation. Virulence. 2011, 2 (3): 241-243. 10.4161/viru.2.3.16193.View ArticlePubMedGoogle Scholar
- Holt KE, Thomson NR, Wain J, Langridge GC, Hasan R, Bhutta ZA, Quail MA, Norbertczak H, Walker D, Simmonds M: Pseudogene accumulation in the evolutionary histories of Salmonella enterica serovars Paratyphi A and Typhi. BMC Genomics. 2009, 10: 36-10.1186/1471-2164-10-36.PubMed CentralView ArticlePubMedGoogle Scholar
- Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill FX, Goodhead I, Rance R, Baker S, Maskell DJ, Wain J: High-throughput sequencing provides insights into genome variation and evolution in Salmonella Typhi. Nat Genet. 2008, 40 (8): 987-993. 10.1038/ng.195.PubMed CentralView ArticlePubMedGoogle Scholar
- Losada L, Ronning CM, DeShazer D, Woods D, Fedorova N, Kim HS, Shabalina SA, Pearson TR, Brinkac L, Tan P: Continuing evolution of Burkholderia mallei through genome reduction and large-scale rearrangements. Genome Biol Evol. 2010, 2: 102-116. 10.1093/gbe/evq003.PubMed CentralView ArticlePubMedGoogle Scholar
- Parkhill J, Sebaihia M, Preston A, Murphy LD, Thomson N, Harris DE, Holden MT, Churcher CM, Bentley SD, Mungall KL: Comparative analysis of the genome sequences of Bordetella pertussis, Bordetella parapertussis and Bordetella bronchiseptica. Nat Genet. 2003, 35 (1): 32-40. 10.1038/ng1227.View ArticlePubMedGoogle Scholar
- Moran NA, Plague GR: Genomic changes following host restriction in bacteria. Curr Opin Genet Dev. 2004, 14 (6): 627-633. 10.1016/j.gde.2004.09.003.View ArticlePubMedGoogle Scholar
- Stinear TP, Seemann T, Pidot S, Frigui W, Reysset G, Garnier T, Meurice G, Simon D, Bouchier C, Ma L: Reductive evolution and niche adaptation inferred from the genome of Mycobacterium ulcerans, the causative agent of Buruli ulcer. Genome Res. 2007, 17 (2): 192-200. 10.1101/gr.5942807.PubMed CentralView ArticlePubMedGoogle Scholar
- Song H, Hwang J, Yi H, Ulrich RL, Yu Y, Nierman WC HSK: The early stage of bacterial genome-reductive evolution in the host. PLoS Pathog. 2010, 6: e1000922-10.1371/journal.ppat.1000922.PubMed CentralView ArticlePubMedGoogle Scholar
- Farley M: Group B streptococcal disease in nonpregnant adults. Clin Infect Dis. 2001, 33: 556-561. 10.1086/322696.View ArticlePubMedGoogle Scholar
- Schuchat A: Epidemiology of group B streptococcal disease in the United States: shifting paradigms. Clin Microbiol Rev. 1998, 11: 497-PubMed CentralPubMedGoogle Scholar
- Dermer P, Lee C, Eggert J, Few B: A history of neonatal group B streptococcus with its related morbidity and mortality rates in the United States. J Pediatr Nurs. 2004, 19: 357-363. 10.1016/j.pedn.2004.05.012.View ArticlePubMedGoogle Scholar
- van der Mee-Marquet N, Fourny L, Arnault L, Domelier A, Salloum M, Lartigue M, Quentin R: Molecular characterization of human-colonizing Streptococcus agalactiae strains isolated from throat, skin, anal margin, and genital body sites. J Clin Microbiol. 2008, 46: 2906-2911. 10.1128/JCM.00421-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Keefe G: Streptococcus agalactiae mastitis: a review. Can Vet J. 1997, 38: 429-437.PubMed CentralPubMedGoogle Scholar
- Mian G, Godoy D, Leal C, Yuhara T, Costa G, Figueiredo H: Aspects of the natural history and virulence of S. agalactiae infection in Nile tilapia. Vet Microbiol. 2009, 136: 180-183. 10.1016/j.vetmic.2008.10.016.View ArticlePubMedGoogle Scholar
- Vandamme P, Devriese L, Pot B, Kersters K, Melin P: Streptococcus difficile is a nonhemolytic group B, type Ib Streptococcus. Int J Syst Bacteriol. 1997, 47: 81-85. 10.1099/00207713-47-1-81.View ArticlePubMedGoogle Scholar
- Wilkinson H, Thacker LG, Facklam RR: Nonhemolytic group B streptococci of human, bovine, and ichthyic origin. Infect Immun. 1973, 7: 496-498.PubMed CentralPubMedGoogle Scholar
- Amborski RL, Snider TG, Thune RL, Culley DD: A non-hemolytic, group B streptococcus infection in cultured bullfrogs, Rana catesbeiana, in Brazil. J Wildlife Dis. 1983, 19: 180-184.View ArticleGoogle Scholar
- Evans JJ, Bohnsack JF, Klesius PH, Whiting AA, Garcia JC, Shoemaker CA, Takahashi S: Phylogenetic relationships among Streptococcus agalactiae isolated from piscine, dolphin, bovine and human sources: a dolphin and piscine lineage associated with a fish epidemic in Kuwait is also associated with human neonatal infections in Japan. J Med Microbiol. 2008, 57: 1369-1376. 10.1099/jmm.0.47815-0.View ArticlePubMedGoogle Scholar
- Jones N, Bohnsack JF, Takahashi S, Oliver KA, Chan MS, Kunst F, Glaser P, Rusniok C, Crook DW, Harding RM: Multilocus sequence typing system for group B streptococcus. J Clin Microbiol. 2003, 41: 2530-2536. 10.1128/JCM.41.6.2530-2536.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Sørensen UB, Poulsen K, Ghezzo C, Margarit I, Kilian M: Emergence and global dissemination of host-specific Streptococcus agalactiae clones. MBio. 2010, 1: e00178-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Glaser P, Rusniok C, Buchrieser C, Chevalier F, Frangeul L, Msadek T, Zouine M, Couvé E, Lalioui L, Poyart C: Genome sequence of Streptococcus agalactiae, a pathogen causing invasive neonatal disease. Mol Microbiol. 2002, 45: 1499-1513. 10.1046/j.1365-2958.2002.03126.x.View ArticlePubMedGoogle Scholar
- Richards VP, Lang P, Pavinski Bitar PD, Lefebure T, Schukken YH, Zadoks RN, Stanhope MJ: Comparative genomics and the role of lateral gene transfer in the evolution of bovine adapted Streptococcus agalactiae. Infect Genet Evol. 2011, 11: 1263-1275. 10.1016/j.meegid.2011.04.019.PubMed CentralView ArticlePubMedGoogle Scholar
- Tettelin H, Masignani V, Cieslewicz M, Eisen J, Peterson S, Wessels M, Paulsen I, Nelson K, Margarit I, Read T: Complete genome sequence and comparative genomic analysis of an emerging human pathogen, serotype V Streptococcus agalactiae. Proc Natl Acad Sci U S A. 2002, 99: 12391-12396. 10.1073/pnas.182380799.PubMed CentralView ArticlePubMedGoogle Scholar
- Tettelin H, Masignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, Angiuoli SV, Crabtree J, Jones AL, Durkin AS: Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial “pan-genome”. Proc Natl Acad Sci U S A. 2005, 102 (39): 13950-13955. 10.1073/pnas.0506758102.PubMed CentralView ArticlePubMedGoogle Scholar
- Eldar A, Bejerano Y, Bercovier H: Streptococcus shiloi and Streptococcus difficile: two new streptococcal species causing a meningoencephalitis in fish. Curr Microbiol. 1994, 28: 139-143. 10.1007/BF01571054.View ArticleGoogle Scholar
- Eldar A, Bejerano Y, Livoff A, Horovitcz A, Bercovier H: Experimental streptococcal meningo-encephalitis in cultured fish. Vet Microbiol. 1995, 43: 33-40. 10.1016/0378-1135(94)00052-X.View ArticlePubMedGoogle Scholar
- Delannoy CM, Crumlish M, Fontaine MC, Pollock J, Foster G, Dagleish MP, Turnbull JF, Zadoks RN: Human Streptococcus agalactiae strains in aquatic mammals and fish. BMC Microbiol. 2013, 13: 41-10.1186/1471-2180-13-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Delannoy CM, Zadoks RN, Lainson FA, Ferguson HW, Crumlish M, Turnbull JF, Fontaine MC: Draft genome sequence of a nonhemolytic fish-pathogenic Streptococcus agalactiae strain. J Bacteriol. 2012, 194: 6341-6342. 10.1128/JB.01552-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu G, Zhang W, Lu C: Complete genome sequence of Streptococcus agalactiae GD201008-001, isolated in China from tilapia with meningoencephalitis. J Bacteriol. 2012, 194: 6653-10.1128/JB.01788-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang B, Jian J, Lu Y, Cai S, Huang Y, Tang J, Wu Z: Complete genome sequence of streptococcus agalactiae ZQ0910, a pathogen causing meningoencephalitis in the GIFT strain of nile tilapia (Oreochromis niloticus). J Bacteriol. 2012, 194: 5132-10.1128/JB.01080-12.PubMed CentralView ArticlePubMedGoogle Scholar
- Brochet M, Rusniok C, Couvé E, Dramsi S, Poyart C, Trieu-Cuot P, Kunst F, Glaser P: Shaping a bacterial genome by large chromosomal replacements, the evolutionary history of Streptococcus agalactiae. Proc Natl Acad Sci U S A. 2008, 105: 15961-15966. 10.1073/pnas.0803654105.PubMed CentralView ArticlePubMedGoogle Scholar
- Elliott JA, Facklam RR, Richter CB: Whole-cell protein patterns of nonhemolytic group B, type lb, streptococci isolated from humans, mice, cattle, frogs, and fish. J Clin Microbiol. 1990, 28: 628-630.PubMed CentralPubMedGoogle Scholar
- Olivares-Fuster O, Klesius PH, Evans J, Arias CR: Molecular typing of Streptococcus agalactiae isolates from fish. J Fish Dis. 2008, 31: 277-283. 10.1111/j.1365-2761.2007.00900.x.View ArticlePubMedGoogle Scholar
- Michel C, Pelletier C, Boussaha M, Douet DG, Lautraite A, Tailliez P: Diversity of lactic acid bacteria associated with fish and the fish farm environment, established by amplified rRNA gene restriction analysis. Appl Environ Microbiol. 2007, 73: 2947-2955. 10.1128/AEM.01852-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Evans JJ, Klesius PH, Pasnik DJ, Bohnsack JF: Human Streptococcus agalactiae isolate in Nile tilapia (Oreochromis niloticus). Emerg Infect Dis. 2009, 15: 774-776. 10.3201/eid1505.080222.PubMed CentralView ArticlePubMedGoogle Scholar
- Helm RA, Lee AG, Christman HD, Maloy S: Genomic rearrangements at rrn operons in Salmonella. Genetics. 2003, 165: 951-959.PubMed CentralPubMedGoogle Scholar
- Kothapalli S, Nair S, Alokam S, Pang T, Khakhria R, Woodward D, Johnson W, Stocker BA, Sanderson KE, Liu SL: Diversity of genome structure in Salmonella enterica serovar Typhi populations. J Bacteriol. 2005, 187: 2638-2650. 10.1128/JB.187.8.2638-2650.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, Salzberg SL: Versatile and open software for comparing large genomes. Genome Biol. 2004, 5: R12-10.1186/gb-2004-5-2-r12.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralView ArticlePubMedGoogle Scholar
- Franken C, Brandt C, Bröker G, Spellerberg B: ISSag1 in streptococcal strains of human and animal origin. Int J Med Microbiol. 2004, 294: 247-254. 10.1016/j.ijmm.2004.04.002.View ArticlePubMedGoogle Scholar
- Lachenauer CS, Creti R, Michel JL, Madoff LC: Mosaicism in the alpha-like protein genes of group B streptococci. Proc Natl Acad Sci U S A. 2000, 97: 9630-9635. 10.1073/pnas.97.17.9630.PubMed CentralView ArticlePubMedGoogle Scholar
- Franken C, Haase G, Brandt C, Weber-Heynemann J, Martin S, Lämmler C, Podbielski A, Lütticken R, Spellerberg B: Horizontal gene transfer and host specificity of beta-haemolytic streptococci: the role of a putative composite transposon containing scpB and lmb. Mol Microbiol. 2001, 41: 925-935.View ArticlePubMedGoogle Scholar
- Shen P, Huang HV: Homologous recombination in Escherichia coli: dependence on substrate length and homology. Genetics. 1986, 112: 441-457.PubMed CentralPubMedGoogle Scholar
- Pierce JC, Kong D, Masker W: The effect of the length of direct repeats and the presence of palindromes on deletion between directly repeated DNA sequences in bacteriophage T7. Nucleic Acids Res. 1991, 19: 3901-3905. 10.1093/nar/19.14.3901.PubMed CentralView ArticlePubMedGoogle Scholar
- Bi X, Liu LF: recA-independent and recA-dependent intramolecular plasmid recombination. Differential homology requirement and distance effect. J Mol Biol. 1994, 235: 414-423. 10.1006/jmbi.1994.1002.View ArticlePubMedGoogle Scholar
- Chedin F, Dervyn E, Dervyn R, Ehrlich SD, Noirot P: Frequency of deletion formation decreases exponentially with distance between short direct repeats. Mol Microbiol. 1994, 12: 561-569. 10.1111/j.1365-2958.1994.tb01042.x.View ArticlePubMedGoogle Scholar
- Lovett ST, Gluckman TJ, Simon PJ, Sutera VA, Drapkin PT: Recombination between repeats in Escherichia coli by a recA-independent, proximity-sensitive mechanism. Mol Gen Genet. 1994, 245: 294-300. 10.1007/BF00290109.View ArticlePubMedGoogle Scholar
- Jonsson IM, Pietrocola G, Speziale P, Verdrengh M, Tarkowski A: Role of fibrinogen-binding adhesin expression in septic arthritis and septicemia caused by Streptococcus agalactiae. J Infect Dis. 2005, 192 (8): 1456-1464. 10.1086/491478.View ArticlePubMedGoogle Scholar
- Schubert A, Zakikhany K, Pietrocola G, Meinke A, Speziale P, Eikmanns BJ, Reinscheid DJ: The fibrinogen receptor FbsA promotes adherence of Streptococcus agalactiae to human epithelial cells. Infect Immun. 2004, 72 (11): 6197-6205. 10.1128/IAI.72.11.6197-6205.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Pietrocola G, Schubert A, Visai L, Torti M, Fitzgerald JR, Foster TJ, Reinscheid DJ, Speziale P: FbsA, a fibrinogen-binding protein from Streptococcus agalactiae, mediates platelet aggregation. Blood. 2005, 105 (3): 1052-1059.View ArticlePubMedGoogle Scholar
- Schubert A, Zakikhany K, Schreiner M, Frank R, Spellerberg B, Eikmanns BJ, Reinscheid DJ: A fibrinogen receptor from group B streptococcus interacts with fibrinogen by repetitive units with novel ligand binding sites. Mol Microbiol. 2002, 46 (2): 557-569. 10.1046/j.1365-2958.2002.03177.x.View ArticlePubMedGoogle Scholar
- Lang S, Xue J, Guo Z, Palmer M: Streptococcus agalactiae CAMP factor binds to GPI-anchored proteins. Med Microbiol Immunol. 2007, 196 (1): 1-10. 10.1007/s00430-006-0021-2.View ArticlePubMedGoogle Scholar
- Yim HH, Nittayarin A, Rubens CE: Analysis of the capsule synthesis locus, a virulence factor in group B streptococci. Adv Exp Med Biol. 1997, 418: 995-997.View ArticlePubMedGoogle Scholar
- Poyart C, Quesnes G, Trieu-Cuot P: Sequencing the gene encoding manganese-dependent superoxide dismutase for rapid species identification of enterococci. J Clin Microbiol. 2000, 38 (1): 415-418.PubMed CentralPubMedGoogle Scholar
- Hamilton A, Popham DL, Carl DJ, Lauth X, Nizet V, Jones AL: Penicillin-binding protein 1a promotes resistance of group B streptococcus to antimicrobial peptides. Infect Immun. 2006, 74 (11): 6179-6187. 10.1128/IAI.00895-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Holmes AR, McNab R, Millsap KW, Rohde M, Hammerschmidt S, Mawdsley JL, Jenkinson HF: The pavA gene of Streptococcus pneumoniae encodes a fibronectin-binding protein that is essential for virulence. Mol Microbiol. 2001, 41: 1395-1408. 10.1046/j.1365-2958.2001.02610.x.View ArticlePubMedGoogle Scholar
- van Sorge NM, Quach D, Gurney MA, Sullam PM, Nizet V, Doran KS: The group B streptococcal serine-rich repeat 1 glycoprotein mediates penetration of the blood-brain barrier. J Infect Dis. 2009, 199 (10): 1479-1487. 10.1086/598217.PubMed CentralView ArticlePubMedGoogle Scholar
- Mello LV, De Groot BL, Li S, Jedrzejas MJ: Structure and flexibility of Streptococcus agalactiae hyaluronate lyase complex with its substrate. Insights into the mechanism of processive degradation of hyaluronan. J Biol Chem. 2002, 277 (39): 36678-36688. 10.1074/jbc.M205140200.View ArticlePubMedGoogle Scholar
- Santi I, Scarselli M, Mariani M, Pezzicoli A, Masignani V, Taddei A, Grandi G, Telford JL, Soriani M: BibA: a novel immunogenic bacterial adhesin contributing to group B streptococcus survival in human blood. Mol Microbiol. 2007, 63 (3): 754-767.PubMedGoogle Scholar
- Li S, Jedrzejas MJ: Hyaluronan binding and degradation by Streptococcus agalactiae hyaluronate lyase. J Biol Chem. 2001, 276 (44): 41407-41416. 10.1074/jbc.M106634200.View ArticlePubMedGoogle Scholar
- Lang S, Palmer M: Characterization of Streptococcus agalactiae CAMP factor as a pore-forming toxin. J Biol Chem. 2003, 278 (40): 38167-38173. 10.1074/jbc.M303544200.View ArticlePubMedGoogle Scholar
- Phillips EA, Tapsall JW, Smith DD: Rapid tube CAMP test for identification of Streptococcus agalactiae (Lancefield group B). J Clin Microbiol. 1980, 12 (2): 135-137.PubMed CentralPubMedGoogle Scholar
- Christie J, McNab R, Jenkinson HF: Expression of fibronectin-binding protein FbpA modulates adhesion in Streptococcus gordonii. Microbiology. 2002, 148 (Pt 6): 1615-1625.View ArticlePubMedGoogle Scholar
- Jedrzejas MJ, Mello LV, de Groot BL, Li S: Mechanism of hyaluronan degradation by Streptococcus pneumoniae hyaluronate lyase. Structures of complexes with the substrate. J Biol Chem. 2002, 277 (31): 28287-28297. 10.1074/jbc.M112009200.View ArticlePubMedGoogle Scholar
- Maisey HC, Doran KS, Nizet V: Recent advances in understanding the molecular basis of group B streptococcus virulence. Expert Rev Mol Med. 2008, 10: e27-PubMed CentralView ArticlePubMedGoogle Scholar
- Locke JB, Colvin KM, Datta AK, Patel SK, Naidu NN, Neely MN, Nizet V, Buchanan JT: Streptococcus iniae capsule impairs phagocytic clearance and contributes to virulence in fish. J Bacteriol. 2007, 189 (4): 1279-1287. 10.1128/JB.01175-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller JD, Neely MN: Large-scale screen highlights the importance of capsule for virulence in the zoonotic pathogen Streptococcus iniae. Infect Immun. 2005, 73: 921-934. 10.1128/IAI.73.2.921-934.2005.PubMed CentralView ArticlePubMedGoogle Scholar
- Spellerberg B, Rozdzinski E, Martin S, Weber-Heynemann J, Lutticken R: rgf encodes a novel two-component signal transduction system of Streptococcus agalactiae. Infect Immun. 2002, 70 (5): 2434-2440. 10.1128/IAI.70.5.2434-2440.2002.PubMed CentralView ArticlePubMedGoogle Scholar
- Quach D, van Sorge NM, Kristian SA, Bryan JD, Shelver DW, Doran KS: The CiaR response regulator in group B streptococcus promotes intracellular survival and resistance to innate immune defenses. J Bacteriol. 2009, 191: 2023-2032. 10.1128/JB.01216-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Lamy MC, Zouine M, Fert J, Vergassola M, Couve E, Pellegrini E, Glaser P, Kunst F, Msadek T, Trieu-Cuot P: CovS/CovR of group B streptococcus: a two-component global regulatory system involved in virulence. Mol Microbiol. 2004, 54 (5): 1250-1268. 10.1111/j.1365-2958.2004.04365.x.View ArticlePubMedGoogle Scholar
- Walker MJ, Hollands A, Sanderson-Smith ML, Cole JN, Kirk JK, Henningham A, McArthur JD, Dinkla K, Aziz RK, Kansal RG: DNase Sda1 provides selection pressure for a switch to invasive group A streptococcal infection. Nat Med. 2007, 13 (8): 981-985. 10.1038/nm1612.View ArticlePubMedGoogle Scholar
- Sumby P, Whitney AR, Graviss EA, DeLeo FR, Musser JM: Genome-wide analysis of group a streptococci reveals a mutation that modulates global phenotype and disease specificity. PLoS Pathog. 2006, 2 (1): e5-10.1371/journal.ppat.0020005.PubMed CentralView ArticlePubMedGoogle Scholar
- Cunin R, Glansdorff N, Piérard A, Stalon V: Biosynthesis and metabolism of arginine in bacteria. Microbiol Rev. 1986, 50: 314-352.PubMed CentralPubMedGoogle Scholar
- Kakinuma Y: Inorganic cation transport and energy transduction in Enterococcus hirae and other streptococci. Microbiol Mol Biol Rev. 1998, 62: 1021-1045.PubMed CentralPubMedGoogle Scholar
- Aury JM, Cruaud C, Barbe V, Rogier O, Mangenot S, Samson G, Poulain J, Anthouard V, Scarpelli C, Artiguenave F: High quality draft sequences for prokaryotic genomes using a mix of new sequencing technologies. BMC Genomics. 2008, 9: 603-10.1186/1471-2164-9-603.PubMed CentralView ArticlePubMedGoogle Scholar
- Zerbino DR, Birney E: Velvet: Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008, 18: 821-829. 10.1101/gr.074492.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Darling ACE, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14 (7): 1394-1403. 10.1101/gr.2289704.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Milne I, Bayer M, Cardle L, Shaw P, Stephen G, Wright F, Marshall D: Tablet-next generation sequence assembly visualization. Bioinformatics. 2010, 26: 401-402. 10.1093/bioinformatics/btp666.PubMed CentralView ArticlePubMedGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunesekaran P, Ceric G, Forslund K: The Pfam protein families database. Nucleic Acids Res. 2010, 38: D211-D222. 10.1093/nar/gkp985.PubMed CentralView ArticlePubMedGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.View ArticlePubMedGoogle Scholar
- Rouillard JM, Herbert CJ, Zuker M: OligoArray: genome-scale oligonucleotide design for microarrays. Bioinformatics. 2002, 18: 486-487. 10.1093/bioinformatics/18.3.486.View ArticlePubMedGoogle Scholar
- Dramsi S, Dubrac S, Konto-Ghiorghi Y, Da Cunha V, Couvé E, Glaser P, Caliot E, Débarbouillé M, Bellais S, Trieu-Cuot P: Rga, a RofA-like regulator, is the major transcriptional activator of the PI-2a pilus in Streptococcus agalactiae. Microb Drug Resist. 2012, 18: 286-297. 10.1089/mdr.2012.0005.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.