Analysis of growth-phase regulated genes in Streptococcus agalactiaeby global transcript profiling
BMC Microbiology volume 9, Article number: 32 (2009)
Bacteria employ multiple mechanisms to control gene expression and react to their constantly changing environment. Bacterial growth in rich laboratory medium is a dynamic process in which bacteria utilize nutrients from simple to complex and change physical properties of the medium, as pH, during the process. To determine which genes are differentially expressed throughout growth from mid log to stationary phase, we performed global transcript analysis.
The S. agalactiae transcriptome is dynamic in response to growth conditions. Several genes and regulons involved in virulence factor production and utilization of alternate carbon sources were differentially expressed throughout growth.
These data provide new information about the magnitude of plasticity of the S. agalactiae transcriptome and its adaptive response to changing environmental conditions. The resulting information will greatly assist investigators studying S. agalactiae physiology and pathogenesis.
Bacteria employ multiple mechanisms to control gene expression and react to their constantly changing environment. These processes are especially critical for bacterial pathogens to survive and cause disease in humans and other hosts. Global control of gene expression is achieved using alternative sigma factors, two-component systems (TCSs), small regulatory RNAs, regulators such as RelA and LuxS, or concerted action of regulons (for a review see [1–6] and references therein). Gram positive pathogens such as group A Streptococcus (S. pyogenes, GAS) and group B Streptococcus (S. agalactiae, GBS) lack (or have limited number) of alternative sigma factors of fully confirmed function [7–9]. Analyses of global transcription in GAS under various growth conditions including saliva, blood, and tissue has shown that environmental response regulation is achieved using other mechanisms such RNA stability , "stand alone" regulators such as mga , or TCSs [12–15]. These transcriptome analyses have been especially useful in providing new information about microbial physiology and leads for pathogenesis research. However, the transcriptional response of GBS to changing growth conditions has not been fully analyzed, only single reports were recently published . GBS is an important human and cow pathogen, responsible for thousands of severe invasive infections in man and large economic loss attributable to bovine mastitis (see [17, 18] and references therein).
One of the best examples of sequential gene regulation is bacterial growth in complex medium and activation of stationary phase genes. During growth, bacteria utilize available nutrients, presumably from simple to more complex, and alter their environment (e.g. decrease or increase in pH) as a result of metabolic byproduct release. Therefore, stationary phase can be considered the acid/alkali stress, depending on the type of metabolism and nutrients utilized. GAS grown to stationary phase sequentially expresses genes involved in various aspects of GAS physiology, metabolism and virulence, many genes activated or repressed during the transition to stationary phase have also been shown to play a role in GAS virulence . The purpose of the present study was to identify growth phase-regulated genes in GBS, with a special interest in providing new information about virulence factor gene expression.
Sample collection for microarray analysis
GBS strain NEM316  was grown as three static cultures (3 biological replicates) in liquid Todd Hewitt medium with 0.5% yeast extract in the 5% CO2 atmosphere at 37°C . Samples were collected at OD600 approximately 0.5, 1.0, 2.0 and 2.5, representing mid-logarithmic (ML), late-logarithmic (LL), early stationary (ES) and stationary (S, about 3 h from entering the phase) growth phases, respectively. Growth curve of bacterial cultures used for data collection is presented as Figure 1. Five ml of each sample were immediately mixed after collection with 10 ml of RNAProtect (Qiagen), centrifuged and stored at -80°C until processing.
Glucose content of the medium at the beginning and end of the culture was measured using Optium Xido glucometer (Abbot) and pH was checked using pH test strips (Macherey Nagel).
GBS cells were mechanically opened by shaking with glass beads (Lysing Matrix B, MPBio) and TRIZOL (Invitrogen). RNA was isolated according to Chomczynski and Sacchi , with an additional purification step using RNeasy columns (Qiagen). Targets for hybridization with the array were prepared according to array manufacturer (Affymetrix) as described previously .
Array hybridization and data acquisition
The custom expression array manufactured by Affymetrix  contained redundant sets of probes representing 1,994 open reading frames (ORFs) of previously sequenced GBS strain NEM316 . Arrays were hybridized and scanned according to the manufacturers instructions. The results of hybridization were normalized to mean of total intensity of GBS probes to allow multiple time point comparison. Array hybridization results are presented as Additional file 1 and are deposited in GEO database http://www.ncbi.nlm.nih.gov/projects/geo/ under GSE12238 accession number.
Results and Discussion
General trends in transcription
After determining transcript levels for all probe sets, the 1,994 transcripts were grouped into 15 clusters based on their behavior during growth (Figure 2) (self organizing map algorithm; Array Assist 5.1.0 package, Stratagene). The clusters were grouped into five main categories. The first 3 categories contain genes whose transcription did not correspond to growth phase, and were either expressed at low (cluster 0), medium (clusters 6, 7), or high (clusters 8, 9) levels in all phases of growth. Category 4 genes (clusters 1–4) exhibited increased transcription in ES or S phase, and category 5 genes (clusters 5, 10–14) had transcription levels that peaked in ML phase and decreased into S phase.
Genes exhibiting growth phase-independent transcription
Genes in clusters 6, 7, 8, and 9 did not show growth phase-dependent transcriptional regulation. The genes are clustered instead based on their transcript level and general profile. Clusters 6 and 7 contain genes that are expressed at the same level until ES phase to slightly lower expression in S phase. Clusters 8 and 9 contain genes, which the transcript level is steady or slightly increases over time. Cluster 9 is especially interesting in that it contains a group of highly expressed genes that includes fabF, fabZ, fabH, and accBCD (gbs0331, 0336–0341) encoding subunits of beta-ketoacyl-ACP synthase, subunits of acetyl-CoA carboxylase, 3-hydroxydecanoyl-ACP dehydratase, and biotin carboxylase. Other genes in cluster 9 involved in energy production are ATP synthase subunits (atpABEF, gbs 0875–7 and 9). Interestingly, cluster 9 contains a transcript of putative catabolite control protein A (ccpA), and the amount grows steadily to increase about three-fold in S phase in comparison with ML (Table 1). CcpA is a major mediator of carbon catabolite repression – the control mechanism of nutrient utilization. In GAS, CcpA has recently been shown to be a critical direct link between carbohydrate utilization and virulence . Function of CcpA in GBS has been not experimentally confirmed yet. Based on the consensus CcpA binding site (cre sequence), we detected that genome of NEM316 strain contains multiple putative cre sites in promoter sequences of multiple genes (Table 2), what might be correlated with changes in expression of genes involved in arginine and carbohydrate metabolism (see below). The transcript encoding HPr carrier protein, another element of the CcpA regulatory pathway in Gram-positive bacteria, also belongs to cluster 9. HPr kinase, however, is an S phase-related gene (see below).
Log phase related genes
Almost 50% of all GBS transcripts that were represented on the chip had similar patterns of expression and were classified into clusters 5, 10, 11, 12, 13, and 14. Transcript levels peaked in ML phase and decreased gradually to their lowest levels in S phase. These six clusters differ in their basal level of expression in L phase. The genes assigned to cluster 5 were expressed at low levels in ML phase, whereas genes in cluster 14 had very high transcripts in ML phase. Cluster 5 contains genes involved in multiple cellular and metabolic processes, whereas cluster 14 genes are involved predominantly in synthesis of ribosomal proteins. Clusters 12–14 contain genes encoding RNA polymerase subunits (gbs0084, gbs0105, gbs0156/7, gbs0302) that are down regulated from -2.3 to -10 times, which likely indicates a slowing of gene transcription. RpoD (gbs1496, encoding the major σ70) is also down regulated (~-3×). The RpoE subunit (gbs0105) plays a role in the development of sepsis during GBS infection , and its down regulation during growth might have consequences for GBS virulence.
S phase related genes
We identified a group of known stress response genes present in clusters 1–4 that were significantly up-regulated in S phase, including hrcA, grpE, dnaK (gbs0094–96), clpE, and clpL (gbs0535 and gbs1376). Transcription of genes putatively involved in the stress response such as Gls24 and universal stress response family proteins gbs1202/1204, gbs1721, and gbs1778 were also elevated in S phase compared to ML phase (Table 1).
Two apparent operons responsible for arginine/ornithine transport and metabolism were also among the group of highly transcribed S phase genes. One operon (gbs2083–2085) encodes an arginine/ornithine antiporter, carbamate kinase, and ornithine carbamoyltransferase, respectively, and is up-regulated 350 to >1,000 times. The second operon (gbs2122–2126) encodes arginine deiminase, a second ornithine carbamoyltransferase, a second arginine/ornithine antiporter, and another carbamate kinase and is up-regulated ~55 to 150 times. Enzymes encoded by genes in these apparent operons are involved in arginine fermentation via the arginine deiminase pathway. They allow GBS to use arginine as an energy source after simple carbohydrates are exhausted from the medium, as would occur during stationary phase. On the other hand, activation of arginine deiminase pathway might have protective function against acidic conditions in a way similar to oral Streptococci  as we observed decrease of pH from about 7.9 to 5.5 between ML and S growth phases.
Metabolic changes toward the utilization of increasingly complex nutrient and carbon sources (see below) can be reflected by changes in utilization of simple carbohydrates (drop in the glucose concentration in the medium from ~300 mg/ml in ML to non detectable level in S) and by changes in transcription of the glpKDF (gbs0263–5, +45 to +63 times), a putative operon responsible for glycerol uptake and utilization.
Trends in expression of genes involved in nutrient utilization and energy metabolism
In contrast to genes involved in other aspects of GBS metabolism and physiology, the only genes significantly up-regulated in S phase compared to ML were involved in carbohydrate metabolism (Figure 3). For example, we observed increased levels of certain glycolytic enzymes such as fructose-bisphosphate aldolase (gbs0125), glyceraldehyde 3P-dehydrogenase (gbs1811), phosphoglycerate kinase (gbs1809), enolase (gbs0608), pyruvate dehydrogenase (acoAB), and L-lactate dehydrogenase (gbs0947) (Table 1). This finding is similar to the results reported recently by Chaussee et al  showing that transcripts encoding proteins involved in carbohydrate utilization and transport were more abundant in S phase, presumably to maximize carbohydrate utilization. The authors suggested that increased transcription of genes involved in central metabolism and sequential utilization of more complex carbohydrates might be a particularly useful adaptation during infection of tissues where the concentration of carbohydrates is low . In GAS, transcripts of genes involved in transport and metabolism of lactose, sucrose, mannose, and amylase were also more abundant during the stationary phase of growth , similar to our findings in GBS (Additional file 2). Similar to links between carbohydrate metabolism and virulence in GAS , also carbohydrate metabolism in GBS might be connected to strain invasiveness and strain tissue-disease specificity .
Changes in expression of regulators and signal transduction systems
TCSs are especially important in the control of global gene expression, especially in the absence of alternative sigma factors. Of the multiple TCSs in GBS, only covR/S (gbs 1671/2) has been well characterized. CovR/S in GBS controls expression of multiple virulence factors, such as hemolysin, CAMP factor, and multiple adhesins . The transcript levels of covR/S are down regulated in S phase, which may be responsible for the observed changes in transcription of virulence factors such as cyl genes encoding hemolysin. However, because the putative effect of CovRS on the camp and cyl genes seems to be opposite to those observed in covRS NEM316 mutant  it suggests that these genes are under influence of additional regulators.
Several other GBS genes encoding putative TCSs and regulators had significant changes in transcript levels during the growth phases studied. For example, transcript levels of gbs1908/9 increased 10/14 times between ML and S phases. The GAS homologs (M5005_Spy_0830/1 in strain MGAS5005 and Spy1106/7 in strain SF370) regulate an operon located downstream that encodes NAD-dependent malic enzyme and malate-sodium symport. In a ΔM5005_Spy_0830 deletion strain, the transcript of these downstream genes is decreased 23/40 times , indicating positive regulation. Organization of this chromosomal region in GBS is very similar to GAS, and gbs1906 and gbs1907 encode putative homologues to the GAS NAD-dependent malic enzyme and malate-sodium symport proteins, respectively. Genes gbs1906/7 are 63/81 times up-regulated in S phase; therefore this operon appears to be regulated in a similar manner in both GBS and GAS.
The transcript level of another GAS TCS homolog, gbs1934/5, is also elevated. Gbs1934/5 has close identity (~85%) with GAS M5005_Spy_0785/6 (Spy1061/2 in strain SF370), a TCS that has been implicated in the regulation of the mannose/fructose-specific phosphotransferase (PTS) system . Interestingly, in GBS there is also a homolog of this PTS system located directly downstream of gbs1934/5 that is highly up-regulated (46.5 to 468 times) in S phase. Therefore, based on gene position, homology, and transcription regulation patterns, it is reasonable to speculate that these genes function similarly in GBS and GAS.
The possible functions of other TCSs can be inferred from their position. Two sets of TCSs are located directly upstream (gbs2081/2) and downstream (gbs2086/7) of an operon with arginine catabolism genes that are highly up regulated in S phase (see above). The transcript levels of both TCSs change dynamically during growth (Table 1 and Additional file 2). It is probable that genes encoding arginine catabolism proteins might be under tight control of both or either TCS. However, this needs to be confirmed experimentally. Thus, our transcript profiling results are consistent with the hypothesis that in the absence of global response gene regulation medicated by alternative sigma factors, GBS uses multiple TCSs as key mediators regulating the response to changes in the environment (Table 1).
Among putative regulators of unknown function, the highest changes were observed for gbs0191 encoding a transcriptional antiterminator of the BglG family (+50 times, putative CcpA binding site) and gbs0469 (-34 times). Surprisingly, we observed down regulation of expression of other global regulators that are associated with stress and the stringent response to starvation. These include the gene relA (gbs1928) that encodes a putative GTP pyrophosphokinase (-50), codY (gbs1719; -8), the cell density dependent regulator luxS (-3), and the putative mecA (gbs0135) homolog (-20). This result was unexpected given that relA, codY, and luxS are up-regulated in S phase GAS .
Transcripts of proven or putative virulence genes
We observed changes in the transcript level of multiple genes encoding proteins with a carboxyterminus cell-wall anchoring motif. The putative location off the proteins on the cell surface suggests that they may play a role in GBS virulence or pathogen-host interaction. Four transcripts were significantly up-regulated in S phase gbs1420 (+6.3), encoding choline-binding protein, gbs1539 (+4.7) and gbs1929 (+5.5) encoding a putative nucleotidase, and gbs1143 (+2.6). We also observed down regulation in S phase of transcripts for several cell wall anchored proteins including a paralog of C5A peptidase precursor gbs0451 (-2), gbs1104 (-6.2), putative adhesin gbs1529 (-11) and fbp (gbs0850, -3), and putative laminin binding proteins (gbs1307, gbs1926; -3). Down regulation in S phase of proteins involved in bacterial attachment is consistent with results reported for GAS [14, 15, 19]. It is believed that several cell surface proteins are produced during the initial stages of infection to promote adhesion, and later are down-regulated to avoid immune detection.
Other known virulence factors of GBS that showed decreased transcription in S phase included an operon encoding hemolysin (gbs0644–0654), genes encoded on the putative pathogenicity island IX (gbs1061–1076), the putative group B antigen (gbs1478/9, gbs1481, gbs1484/5, gbs1492–1494), and genes involved in capsule synthesis (gbs1233–1247). The putative kinase cpsX (gbs1250) was upregulated 4.4 times (Table 1). Down regulation of capsule and putative and known surface antigens is known to occur in GAS [14, 15, 19]. For example, capsule, an antiphagocytic factor, is expressed during establishment of GAS infection and is later down-regulated once the infection is established [14, 15]. Our results imply a similar scenario could be occurring in GBS. The only transcript encoding a proven virulence factor that was increased in S phase was CAMP factor (+11.6, cfa, gbs2000).
Our results demonstrate that GBS gene transcript levels are highly dynamic throughout the growth cycle in vitro, likely reflecting exposure to an environment that is altering significantly during growth. The organism activates genes involved in metabolism of nutrients and carbon sources other than glucose such as complex carbohydrates and arginine and protect against changing pH. GBS slows down cell division and decreases transcription and translation. Production of virulence factors involved in establishment of the infection is reduced during growth. The global changes of transcript profiles we identified in GBS grown in rich medium are similar to patterns exhibited by GAS. Our results provide new information useful for the study of pathogen-host interactions and gene regulation in pathogenic bacteria.
Commichau FM, Forchhammer K, Stulke J: Regulatory links between carbon and nitrogen metabolism. Curr Opin Microbiol. 2006, 9: 167-172. 10.1016/j.mib.2006.01.001.
Gruber TM, Gross CA: Multiple sigma subunits and the partitioning of bacterial transcription space. Annu Rev Microbiol. 2003, 57: 441-466. 10.1146/annurev.micro.57.030502.090913.
Laub MT, Goulian M: Specificity in two-component signal transduction pathways. Annu Rev Genet. 2007, 41: 121-145. 10.1146/annurev.genet.41.042007.170548.
Nascimento MM, Lemos JA, Abranches J, Lin VK, Burne RA: Role of RelA of Streptococcus mutans in global control of gene expression. J Bacteriol. 2008, 190: 28-36. 10.1128/JB.01395-07.
Storz G: An expanding universe of noncoding RNAs. Science. 2002, 296: 1260-1263. 10.1126/science.1072249.
Xavier KB, Bassler BL: LuxS quorum sensing: more than just a numbers game. Curr Opin Microbiol. 2003, 6: 191-197. 10.1016/S1369-5274(03)00028-6.
Glaser P, Rusniok C, Buchrieser C, Chevalier F, Frangeul L, Msadek T, Zouine M, Couve E, Lalioui L, Poyart C, Trieu-Cuot P, Kunst F: 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.
Opdyke JA, Scott JR, Moran CP: A secondary RNA polymerase sigma factor from Streptococcus pyogenes. Mol Microbiol. 2001, 42: 495-502. 10.1046/j.1365-2958.2001.02657.x.
Tettelin H, Masignani V, Cieslewicz MJ, Donati C, Medini D, Ward NL, Angiuoli SV, Crabtree J, Jones AL, Durkin AS, Deboy RT, Davidsen TM, Mora M, Scarselli M, Ros I, Peterson JD, Hauser CR, Sundaram JP, Nelson WC, Madupu R, Brinkac LM, Dodson RJ, Rosovitz MJ, Sullivan SA, Daugherty SC, Haft DH, Selengut J, Gwinn ML, Zhou L, Zafar N, Khouri H, Radune D, Dimitrov G, Watkins K, O'Connor KJ, Smith S, Utterback TR, White O, Rubens CE, Grandi G, Madoff LC, Kasper DL, Telford JL, Wessels MR, Rappuoli R, Fraser CM: Genome analysis of multiple pathogenic isolates of Streptococcus agalactiae: implications for the microbial "pan-genome". Proc Natl Acad Sci USA. 2005, 102: 13950-13955. 10.1073/pnas.0506758102.
Barnett TC, Bugrysheva JV, Scott JR: Role of mRNA stability in growth phase regulation of gene expression in the group A streptococcus. J Bacteriol. 2007, 189: 1866-1873. 10.1128/JB.01658-06.
Hondorp ER, McIver KS: The Mga virulence regulon: infection where the grass is greener. Mol Microbiol. 2007, 66: 1056-1065. 10.1111/j.1365-2958.2007.06006.x.
Sitkiewicz I, Musser JM: Expression microarray and mouse virulence analysis of four conserved two-component gene regulatory systems in group a streptococcus. Infect Immun. 2006, 74: 1339-1351. 10.1128/IAI.74.2.1339-1351.2006.
Shelburne SA, Sumby P, Sitkiewicz I, Granville C, DeLeo FR, Musser JM: Central role of a bacterial two-component gene regulatory system of previously unknown function in pathogen persistence in human saliva. Proc Natl Acad Sci USA. 2005, 102: 16037-16042. 10.1073/pnas.0505839102.
Graham MR, Virtaneva K, Porcella SF, Barry WT, Gowen BB, Johnson CR, Wright FA, Musser JM: Group A Streptococcus transcriptome dynamics during growth in human blood reveals bacterial adaptive and survival strategies. Am J Pathol. 2005, 166: 455-465.
Graham MR, Virtaneva K, Porcella SF, Gardner DJ, Long RD, Welty DM, Barry WT, Johnson CA, Parkins LD, Wright FA, Musser JM: Analysis of the transcriptome of group A Streptococcus in mouse soft tissue infection. Am J Pathol. 2006, 169: 927-942. 10.2353/ajpath.2006.060112.
Johri AK, Margarit I, Broenstrup M, Brettoni C, Hua L, Gygi SP, Telford JL, Grandi G, Paoletti LC: Transcriptional and proteomic profiles of group B Streptococcus type V reveal potential adherence proteins associated with high-level invasion. Infect Immun. 2007, 75: 1473-1483. 10.1128/IAI.00638-06.
Keefe GP: Streptococcus agalactiae mastitis: a review. Can Vet J. 1997, 38: 429-437.
Muller AE, Oostvogel PM, Steegers EA, Dorr PJ: Morbidity related to maternal group B streptococcal infections. Acta Obstet Gynecol Scand. 2006, 85: 1027-1037.
Chaussee MA, Dmitriev AV, Callegari EA, Chaussee MS: Growth phase-associated changes in the transcriptome and proteome of Streptococcus pyogenes. Arch Microbiol. 2008, 189: 27-41. 10.1007/s00203-007-0290-1.
Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2.
Shelburne SA, Keith D, Horstmann N, Sumby P, Davenport MT, Graviss EA, Brennan RG, Musser JM: A direct link between carbohydrate utilization and virulence in the major human pathogen group A Streptococcus. Proc Natl Acad Sci USA. 2008, 105: 1698-1703. 10.1073/pnas.0711767105.
Jones AL, Needham RH, Rubens CE: The Delta subunit of RNA polymerase is required for virulence of Streptococcus agalactiae. Infect Immun. 2003, 71: 4011-4017. 10.1128/IAI.71.7.4011-4017.2003.
Quivey RG, Kuhnert WL, Hahn K: Genetics of acid adaptation in oral streptococci. Crit Rev Oral Biol Med. 2001, 12: 301-314. 10.1177/10454411010120040201.
Domelier AS, van dM-M, Grandet A, Mereghetti L, Rosenau A, Quentin R: Loss of catabolic function in Streptococcus agalactiae strains and its association with neonatal meningitis. J Clin Microbiol. 2006, 44: 3245-3250. 10.1128/JCM.02550-05.
Lamy MC, Zouine M, Fert J, Vergassola M, Couve E, Pellegrini E, Glaser P, Kunst F, Msadek T, Trieu-Cuot P, Poyart C: CovS/CovR of group B streptococcus: a two-component global regulatory system involved in virulence. Mol Microbiol. 2004, 54: 1250-1268. 10.1111/j.1365-2958.2004.04365.x.
Jiang SM, Ishmael N, Hotopp JD, Puliti M, Tissi L, Kumar N, Cieslewicz MJ, Tettelin H, Wessels MR: Variation in the group B Streptococcus CsrRS regulon and effects on pathogenicity. J Bacteriol. 2008, 190: 1956-1965. 10.1128/JB.01677-07.
Authors would like to thank Kathryn Stockbauer for critical reading of the manuscript.
IS performed the research, IS and JMM analyzed the data and wrote the paper.
Electronic supplementary material
Additional File 1: Supplemental table 1- Normalized hybridization values. File contains normalized hybridization values for each array used in the study. ML-mid logarithmic, LL-late logarithmic, ES-early stationary, S-stationary. P-"present" signal (detected in sample), M-"marginal" signal, A-"absent" signal (not detected). (XLS 1 MB)
Additional File 2: Supplemental table 2 Changes in transcription of 1,994 transcripts present on array (S/ML ratio). Green- genes down regulated in S phase, Red – genes up regulated in S phase, Gray – P values below 0.05. (XLS 416 KB)
About this article
Cite this article
Sitkiewicz, I., Musser, J.M. Analysis of growth-phase regulated genes in Streptococcus agalactiaeby global transcript profiling. BMC Microbiol 9, 32 (2009). https://doi.org/10.1186/1471-2180-9-32