- Research article
- Open Access
Impact of growth pH and glucose concentrations on the CodY regulatory network in Streptococcus salivarius
BMC Genomics volume 19, Article number: 386 (2018)
Streptococcus salivarius is an abundant isolate of the human oral microbiota. Since both pH and glucose availability fluctuate frequently in the oral cavity, the goal of this study was to investigate regulation by CodY, a conserved pleiotropic regulator of Gram positive bacteria, in response to these two signals. The chemostat culture system was employed to precisely control the growth parameters, and the transcriptomes of wild-type S. salivarius 57.I and its CodY-null derivative (ΔcodY) grown at pH 7 and 5.5, with limited and excessive glucose supply were determined.
The transcriptomic analysis revealed that CodY was most active at pH 7 under conditions of glucose limitation. Based on whether a CodY binding consensus could be located in the 5′ flanking region of the identified target, the transcriptomic analysis also found that CodY shaped the transcriptome via both direct and indirect regulation. Inactivation of codY reduced the glycolytic capacity and the viability of S. salivarius at pH 5.5 or in the presence of H2O2. Studies using the Galleria mellonella larva model showed that CodY was essential for the toxicity generated from S. salivarius infection, suggesting that CodY regulation was critical for immune evasion and systemic infections. Furthermore, the CodY-null mutant strain exhibited a clumping phenotype and reduced attachment in biofilm assays, suggesting that CodY also modulates cell wall metabolism. Finally, the expression of genes belonging to the CovR regulon was affected by codY inactivation, but CodY and CovR regulated these genes in opposite directions.
Metabolic adaptation in response to nutrient availability and growth pH is tightly linked to stress responses and virulence expression in S. salivarius. The regulation of metabolism by CodY allows for the maximal utilization of available nutrients and ATP production. The counteractive regulation of the CovR regulon could fine tune the transcriptomes in response to environmental changes.
Global transcriptional regulators execute the co-regulation of multiple pathways simultaneously, which allows microbes to adapt quickly to environmental changes. CodY is one such conserved regulator that controls the expression of both metabolic and virulence genes in response to nutrient starvation in Gram positive bacteria with low GC content [1, 2]. CodY functions primarily as a repressor during exponential growth, and its repressive action is inhibited when nutrients become scarce . The crystal structures of Bacillus subtilis CodY fragments show that CodY harbors a GAF domain at the N-terminal region which interacts with branched-chain amino acids (BCAAs, e.g. isoleucine, leucine, and valine) . The binding of the effectors leads to a change in conformation and activation of the C-terminal helix-turn-helix DNA binding domain [4, 5]. Functional analysis indicates that both BCAAs and GTP could activate the DNA binding activity of CodY in B. subtilis [6,7,8,9], whereas BCAAs are the only known effectors for CodY in Lactococcus lactis , Streptococcus mutans , and Streptococcus pneumoniae . A recent crystal structure study of full-length CodY demonstrates that GTP binds to a site located in the hinge region connecting the N- and C-terminal domains. In comparison to the CodY proteins that are activated by GTP, the CodY proteins of L. lactis and S. pneumoniae exhibit the substitution of three amino acids in the proposed GTP binding site, which may explain the specificity of the cofactors . The CodY binding consensus sequence, 5’-AATTTTCWGAAAATT, has been defined in B. subtilis and L. lactis [14, 15].
Although CodY is activated by BCAAs and mainly acts as a repressor, Listeria monocytogenes CodY remains active in the absence of BCAAs. Under such conditions CodY can activate the expression of prfA, encoding a master virulence regulator , and subsequently initiate a multifaceted regulatory network. In Bacillus thuringiensis, CodY upregulates the PapR-PlcR quorum sensing system by optimizing the uptake of PapR, the signaling peptide of the transcriptional regulator PlcR. The surplus PapR activates the DNA binding activity of PlcR, and subsequently upregulates the expression of the PlcR regulon and virulence of B. thuringiensis . The above examples demonstrate the diversity of CodY regulation.
Streptococcus salivarius is the primary ureolytic species in the human mouth. Since urea is the most abundant nitrogen source in the saliva and crevicular fluids , the presence of ureolytic S. salivarius is critical for maintaining the pH homeostasis of the oral ecosystem and preventing the development of dental caries . Past studies revealed that the expression of S. salivarius 57.I urease operon was upregulated during growth at acidic pH values and in the presence of excess carbohydrate . Under these conditions, ureolysis is beneficial for the survival of S. salivarius . In vivo and in vitro analyses have confirmed that the repression of the urease operon at pH 7 with limited glucose supply is mainly governed by CodY . It is intriguing that the activity of CodY is modulated by growth pH and carbohydrate availability in S. salivarius, as neither factor is known to affect CodY activation. This observation leads to the hypothesis that growth pH and carbohydrate availability modulate the regulatory activity of CodY by altering the intracellular metabolite profile. This study analyzed the transcriptomes of chemostat-grown S. salivarius 57.I and its CodY-null derivative (strain ΔcodY) to construct the regulatory network of CodY in response to different growth conditions. The impact of CodY regulation on the stress responses and pathogenic capacity of S. salivarius was examined.
The binding consensus sequence of CodY in streptococci
To define the CodY binding consensus in streptococcal species, the CodY regulons of Streptococcus agalactiae, Streptococcus dysgalactiae subsp. equisimilis, Streptococcus gallolyticus UCN34, Streptococcus gordonii CH1, Streptococcus mitis B6, S. mutans UA159, S. pneumoniae TIGR4, Streptococcus pyogenes M1, Streptococcus sanguinis SK36, Streptococcus suis 05ZYH33, Streptococcus thermophilus CNRZ1066, and Streptococcus uberis 0140 J were downloaded from the database and their putative CodY binding consensus sequence was determined. The consensus sequence (Additional file 1: Figure S1A), 5’-AATTTTCAGAAAATT, is close to the known CodY binding consensus sequence (Additional file 1: Figure S1B).
An overview of CodY-regulated genes and the role of CodY in global regulation
CodY regulates the expression of S. salivarius urease operon at both pH 7 and 5.5, with a limited (20 mM) or an excessive (100 mM) amount of glucose , suggesting that different levels of CodY activity would be detected in the growth conditions in this study. To test this hypothesis, the codY gene of S. salivarius 57.I was insertionally inactivated to generate strain ΔcodY. Genes that are differentially expressed between wild-type 57.I and strain ΔcodY under each growth condition were identified by using edgeR  and the regulation pattern of these genes was examined by using hierarchical clustering analysis (Fig. 1a). The expression patterns of wild-type 57.I and strain ΔcodY under different growth conditions were located in two clusters, demonstrating the global impact of CodY regulation in response to growth pH values and glucose concentrations (Fig. 1a). Furthermore, the expression patterns of wild-type 57.I grown in 20 mM and 100 mM glucose were separated, suggesting that the amount of glucose in the culture medium modulates CodY activity more effectively than pH. Finally, the expression pattern of strain ΔcodY under 20 mM glucose at pH 7 was most distant from that of the others, suggesting that CodY is most active at pH 7 with limited glucose supply.
As both activation and repression were observed (Fig. 1a), whether CodY affects the transcription profile of S. salivarius by both direct and indirect regulation was analyzed. Genes whose expression was affected significantly by codY deletion under at least one condition were identified. The direct targets were defined based on the presence of a CodY binding consensus in the 5′ flanking region of the gene or operon. A total of 153 genes was found; 83 and 70 genes are defined as direct and indirect targets of CodY, respectively. Hierarchical clustering analysis also indicated that the expression profile of strain ΔcodY under 20 mM glucose at pH 7 was most distant from the profiles of the other three conditions, in both direct (Fig. 1b) and indirect targets (Fig. 1c). While all direct targets were upregulated in strain ΔcodY, both upregulation and downregulation were observed in the indirect targets. The RPKM (Reads Per Kb per Million mapped reads) values of direct and indirect targets are listed in Additional file 2: Table S1 and S2, respectively. The targets are discussed below.
Genes regulated directly by CodY
In agreement with observations in other bacteria, genes encoding proteins catalyzing the uptake and degradation of peptides, the uptake and biosynthesis of amino acids and urease were repressed directly by CodY in S. salivarius 57.I (Table 1). Although repression by CodY was generally more active under glucose limitation compared to that under excessive glucose, the effect of pH on the regulation of oppA (Ssal_00190), pepC (Ssal_00280) and pepP (Ssal_01946) was unclear. On the other hand, in most of the genes involved in the biosynthesis or uptake of amino acids, the expression level between wild-type 57.I and strain ΔcodY was significantly different under at least three conditions, indicating that the differential expression of these genes in response to growth pH values and carbohydrate concentrations optimizes the fitness of S. salivarius under different growth conditions. Finally, all the genes of the urease operon were derepressed at pH 5.5 and upregulation was further enhanced under excessive glucose in S. salivarius 57.I (Additional file 2: Table S1). Additionally, the differences in the expression level between wild-type 57.I and strain ΔcodY under glucose limitation were generally greater than that under excessive glucose, and the effect of codY deletion at pH 5.5 was less evident under excessive glucose, supporting the hypothesis that CodY is least active at pH 5.5 under excessive glucose.
The expression of genes encoding proteins involved in regulation, metal ions transportation, carbon metabolism, virulence, and other cellular functions was also regulated by CodY (Table 1). Specifically, the expression of Ssal_00908, encoding a putative Mn2+/Fe2+ transporter, and of hmuU-fhuC-Ssal_01464, an operon encoding a putative ferric ATP-binding cassette transporter was repressed by CodY under all conditions, suggesting that CodY modulates metal homeostasis in S. salivarius 57.I. CodY also repressed the expression of gdhA (Ssal_01768), encoding glutamate dehydrogenase (GDH), and the nrgA (Ssal_01668)-glnB (Ssal_01667) operon, encoding an ammonium transporter and a PII-like nitrogen regulator, respectively. The repression of these genes was maximized under glucose limitation, indicating that CodY regulates nitrogen utilization in response to carbon availability. Furthermore, the expression of pgmB (Ssal_01810), encoding β-phosphoglucomutase for converting β-D-glucose-1-phosphate to glucose 6-phosphate, was up regulated at pH 7 under excessive glucose in strain ΔcodY. Finally, the identification of lysM (Ssal_01428, encoding a peptidoglycan binding protein) and the cidAB operon (Ssal_00555 and Ssal_00556, encoding a putative murein hydrolase exporter and a regulator) suggests that CodY regulates cell wall biosynthesis.
Indirect targets of CodY
Forty-eight indirect targets were downregulated in strain ΔcodY under at least 3 growth conditions, while 22 genes were upregulated (Table 2). Among the upregulated genes, Ssal_00630, Ssal_00982 and Ssal_01876 are annotated in GenBank as transcriptional regulators. Particularly, the protein product of Ssal_01876, CovR, is a dominant global regulator controlling virulence expression in streptococci [24,25,26]. Sequence analysis revealed that 33 of the indirect targets harbor a CovR binding consensus sequence (5′-ATTARA)  at the 5′ flanking region, and both upregulation (11 genes) and downregulation (22 genes) were observed in strain ΔcodY. Specifically, a putative cell wall-binding hydrolase gene (Ssal_00236) was downregulated in ΔcodY, suggesting that CodY indirectly modulates cell wall structure.
Inactivation of CodY also led to the down-expression of a gene cluster (fabH-acpP-fabKDGF-accB-fabZ-Ssal_01799-accCDA) encoding proteins involved in membrane fatty acid biosynthesis. Although the effects of growth pH values and glucose concentrations on the regulation of these genes are unclear, the differences between wild-type 57.I and ΔcodY in cells grown at pH 7 were generally greater than in those grown at pH 5.5 under the same concentration of glucose, suggesting that CodY was essential for optimal membrane biosynthesis at pH 7.
Inactivation of CodY promoted clumping, prolonged the doubling time, and reduced the biofilm formation in S. salivarius 57.I
To analyze the impact of CodY regulation, a codY complemented strain, CΔcodY, was generated as detailed in the methods, in which a single copy of codY was established on the chromosome of strain ΔcodY. During the preparation of S. salivarius cultures for growth studies, strain ΔcodY formed aggregates in overnight cultures, while wild-type 57.I did not (Fig. 2a). The colony size of strain ΔcodY was smaller than that of the wild-type 57.I (Fig. 2b). To reduce clumping during the course of the growth kinetic studies, the inocula of all strains were gently sonicated prior to inoculation and the cultures were shaken before each reading. The estimated generation time of wild-type 57.I in BHI was approximately 22 min, whereas strain ΔcodY exhibited a generation time of approximately 44 min (Fig. 3). A wild-type doubling time was observed in the complemented strain (CΔcodY), although the lag phase in strain CΔcodY was longer than that of wild-type 57.I under the experimental condition for unknown reasons.
In both the static biofilm cultures (Fig. 4a) and the flow-cell system (Fig. 4b), cells of strain ΔcodY failed to adhere to the plastic surface effectively. Presumably, self-aggregation impeded colonization, leading to minimal biofilm formation. Both wild-type 57.I and strain CΔcodY formed thick biofilms containing both live and dead cells (Fig. 4b).
Inactivation of CodY reduced the glycolytic rate, acid tolerance and the capacity against oxidative stress in S. salivarius 57.I
Regulation by CodY is likely to play a critical role in energy conservation, which could then influence acid tolerance via the activity of F-ATPase . Since bacterial glycolysis is normally acid-limited, and the terminal pH value reflects the sensitivity of cells to acidification, the acid production from glucose and acid tolerance of all strains were evaluated to test this hypothesis. The acid production of strain ΔcodY was lower than that of wild-type 57.I and strain CΔcodY, and the final pH of strain ΔcodY was 0.5 unit higher than that of wild-type 57.I (pH 4.78 vs. pH 4.26) and strain CΔcodY (Fig. 5), suggesting that the acid tolerance of strain ΔcodY was lower than that of wild-type 57.I. Additionally, the expression levels of glk (Ssal_01383), pfk (Ssal_01270) pyk (Ssal_01268), and pgk (Ssal_00232), encoding glucose kinase, 6-phosphofructokinase, pyruvate kinase and phosphoglycerate kinase in the Embden-Meyerhof-Parnas pathway, in wild-type 57.I were higher than those in strain ΔcodY under at least three growth conditions (Additional file 2: Table S3), indicating that CodY positively modulates the glycolytic activity via an unknown mechanism. Finally, the growth of S. salivarius at pH 5.5 was compromised significantly by codY deletion (Fig. 6).
CodY also participates in the regulation of the oxidative stress response. A comparable level of growth was observed in wild-type 57.I and strain CΔcodY in medium containing up to 3 mM H2O2, while reduction in growth was detected with strain ΔcodY (Fig. 7), indicating that inactivation of CodY sensitized S. salivarius to oxidative stress.
CodY regulation was essential for the toxicity generated from S. salivarius infection in the G. mellonella larva model
To analyze whether codY deletion would affect S. salivarius against host innate immune clearance, the toxicity of S. salivarius was evaluated in the G. mellonella larva model (Fig. 8). Seven out of 12 larvae infected with S. salivarius 57.I died 2 days post infection, and greater than 107 CFUs were recovered from one dead larva. A similar death rate was observed with the group of larva infected with strain CΔcodY. Less than 102 CFUs was recovered from one live larva infected with strain ΔcodY 7 days post infection. Thus, CodY regulation was crucial for the survival of S. salivarius and the death of the larvae.
Microbes in the oral cavity frequently encounter fluctuations in carbohydrate availability and pH. This study was designed to gain insight into the impact of these two environmental factors on CodY regulation. Using the chemostat culture system, we found that CodY is most active under glucose limitation. Reduced CodY repression under excessive glucose, when the level of nitrogen containing nutrients is likely to be limiting, could enhance the overall nutrient utilization and energy generation by promoting ammonium uptake and amino acid biosynthesis. Specifically, S. salivarius is able to utilize ammonium as a sole nitrogen source , therefore the expression of the nrgA-glnB operon under excessive glucose allows for the utilization of an alternative nitrogen source. In S. thermophilus , GDH mediates the conversion of carbon nutrients into amino acids to promote lactose utilization. Conceivably, the reduced repression of gdhA by CodY in S. salivarius under excessive glucose might promote the utilization of carbon nutrients by the upregulation of glutamate biosynthesis. On the other hand, variations in the effects of pH between CodY-regulated genes may result from additional pH-dependent regulation by other regulator(s), as seen in the regulation of the urease operon [22, 30]. Taken together, the current data suggest that the profile of intracellular metabolites, rather than the growth pH plays major roles in modulating CodY activity in S. salivarius.
In S. pyogenes, CovR is a master regulator of the CovRS two-component system, which predominantly acts as a transcriptional repressor . Among the identified targets that potentially belong to the CovR regulon, some are activated by CodY in S. salivarius, suggesting that CodY and CovR regulate common genes in opposite directions. Similar regulation has been reported in S. pyogenes, in which CodY modulates CovRS-regulated transcriptomes in response to nutritional conditions . The opposite regulation in response to specific growth conditions may fine tune the transcriptomes to enhance the fitness and virulence of S. salivarius. Thus, regulation by CovR in the absence of CodY modulation contributes, at least in part, to some of the phenotypes of ΔcodY.
The clumping of strain ΔcodY cells in liquid culture, which has not been reported in other CodY-null bacterial strains, is likely to result from CovR regulation in the absence of CodY modulation. It is established that defects in enzymes catalyzing cell wall metabolism could alter cell morphology, the chain length of streptococci , or the cluster size of staphylococci . Among CodY-regulated genes that potentially participate in cell wall metabolism in S. salivarius 57.I, the expression of Ssal_00236, a potential member of the CovR regulon, is most affected by CodY deletion. Ssal_00236 is highly expressed in S. salivarius 57.I, the expression level of which is approximately 3–8 fold higher than that of dnaA under the same growth conditions (data not shown). Inactivation of CodY almost abolished the expression of Ssal_00236. Although the function of the protein encoded by Ssal_00236 and its homologs has not been characterized experimentally in streptococci, this protein is predicted to harbor a LysM motif and a CHAP domain in a manner similar to the Cse protein of S. thermophilus . Cse interacts with the cell wall via the LysM motif and catalyzes cell wall hydrolysis via the endopeptidase activity of the CHAP domain. Inactivation of this locus results in abnormal daughter cell segregation and long chains of cells. The indirct regulation of Ssal_00236 by CodY might be required for the development of proper chain morphology. While the function of the CidAB system in cell wall hydrolysis has been demonstrated in both S. mutans  and S. aureus , in both cases the hydrolytic activity of the CidAB system is counteracted by the activity of the LrgAB system in response to the growth conditions, and neither system is regulated by CodY. Thus, the contribution of the up-expression of the cidAB operon and lysM (Ssal_01428) in the clumping of ΔcodY cells is unclear.
It is likely that CodY participates in acid tolerance of S. salivarius 57.I by regulating the metabolic status, urease production, and possibly, membrane integrity. Regulation of metabolic pathways by CodY could augment ATP generation for acid tolerance, as seen in non-ureolytic S. mutans [11, 37,38,39]. Additionally, the regulation of the urease operon by CodY not only contributes to pH homeostasis but also provides a surplus of nitrogen-containing nutrients, which could again, promote ATP generation for pumping intracellular H+ via the activity of F-ATPase at acidic growth pH values. Finally, studies in S. mutans have demonstrated that the upregulation of unsaturated membrane fatty acid biosynthesis enhances acid tolerance [40, 41], indicating that membrane integrity is essential for optimal acid resistance. Thus, the down regulation of the fabH-acpP-fabKDGF-accB-fabZ-Ssal_01799-accCDA gene cluster in strain ΔcodY could potentially lead to reduced membrane integrity and acid tolerance in S. salivarius.
CodY participates in the oxidative stress response in S. thermophilus  and S. pneumoniae  by activating the expression of genes encoding glutathione synthetase (gshF) and thiol peroxidase (tpxD), respectively. However, neither gshA (Ssal_00663, an ortholog of gshF) nor tpx (Ssal_01091) is regulated by CodY in S. salivarius 57.I. Instead, CodY regulates directly two iron uptake systems (the products of Ssal_00908 and hmuU-fhuC-Ssal_01464). Intracellular iron homeostasis is critical for controlling the generation of reactive oxygen species resulting from the Fenton reaction. In S. pyogenes, both upregulation of an Fe(II) exporter (PmtA)  or downregulation of iron uptake  enhances the oxidative stress resistance. Thus, it is possible that the elevated iron uptake in strain ΔcodY could lead to the perturbation of iron homeostasis, leading to reduced resistance against H2O2.
A putative fibronectin-binding protein (FBP) (encoded by Ssal_00672) was more actively expressed in strain ΔcodY than in wild-type 57.I under all growth conditions. Although homologs of FBP in Gram positive pathogens are known to participate in adherence, invasion and immune evasion [46, 47], the toxicity of strain ΔcodY in the G. mellonella infection model was significantly lower than that of wild-type 57.I. Thus, either Ssal_00672 is not required for the survival against innate immune responses, or the upregulation of Ssal_00672 enhances recognition by host immune cells, and subsequently enhances clearance.
In conclusion, this study demonstrated that metabolic adaptation in response to nutrient availability and growth pH is tightly linked to stress responses and virulence expression in S. salivarius (Fig. 9). The regulation of metabolism by CodY allows for the maximal utilization of available nutrients and ATP production. The counteractive regulation of the CovR regulon fine tunes the expression of the regulon in response to the nutrient availability and growth pH.
Construction of a codY-deficient and a codY-complemented derivative of S. salivarius 57.I
Primers used in this study are listed in Additional file 2: Table S5. Plasmid pSky2 , in which the region encoding the 25th to 78th amino acids (aa) of CodY was replaced by an erythromycin (Em) resistance gene (erm) , was used to insertionally inactivate codY in S. salivarius 57.I. The genotype of the Em resistant isolates was verified by colony PCR using primers Ssal_00403_S and Ssal_00405_AS. The correct strain was designated ΔcodY. An Ωkan-tagged codY gene was then inserted into the intergenic region between Ssal_00777 and Ssal_00779 of strain ΔcodY using ligation mutagenesis  to generate a codY complemented strain. A 1.1-kbp DNA fragment containing the codY gene, its 5′ flanking region of 250 bp and its 3′ flanking region of 90 bp was generated from S. salivarius 57.I with primers 57.I_CodY_800_XhoI_S and 57.I_CodY_1900_SphI_AS. A DNA fragment containing the Ωkan cassette was generated from plasmid pVT924  with primers kan_BamHI_S and kan_XhoI_AS. Two DNA fragments for integration onto the chromosome were generated from S. salivarius 57.I using primer pairs 57.I_lacZ_1241_S/57.I_lacZ_4970_BamHI_AS and 57.I_lacZ_4970_SphI_S/57.I_ssal_00779_5950_AS, respectively. All four PCR products were digested, mixed in a ligation reaction and then used to transform S. salivarius ΔcodY by natural transformation . The genotype of the kanamycin resistant isolates was verified by colony PCR. The correct strain was designated CΔcodY.
Bacterial strains and growth conditions
Bacterial strains used in this study are listed in Additional file 2: Table S4. S. salivarius 57.I  and ΔcodY were grown in a chemostat culture system (Sartorius) in 3% tryptone-0.5% yeast extract medium (TY) supplemented with 20 mM or 100 mM glucose at a dilution rate of 0.3 h− 1, corresponding to a generation time of 2.3 h. Cultures were maintained at pH 7 or pH 5.5 by the addition of 2 N KOH. Cells were grown for a minimum of 10 generations to achieve the steady state at each condition.
RNA-seq library construction, sequencing and gene expression analysis
Total cellular RNA was isolated from streptococcal strains as described  and further purified using the RNeasy kit (Qiagen). Ribosomal RNA was removed from the RNA preparation using the RiboMinus transcriptome isolation kits (Invitrogen). The RNA-seq library was constructed using SOLiD™ small RNA expression kit (ABI). All procedures were performed per manufacturer’s instruction. The sequence of the sample was determined on an ABI SOLiD 4.0 sequencer. Reads with a quality value greater than 8 were mapped to the genome sequence of S. salivarius 57.I (accession number CP002888.1) by using the SOLiD system analysis pipeline tool (Corona Lite) with mismatches up to 5 bases. RPKM were used to calculate the expression level of each gene. The RPKM values were log10 transformed and the correlation between biological and technical replicates under the same growth condition was examined to verify the quality of the analysis. The significant difference in the expression level of each gene between wild-type 57.I and strain ΔcodY under one specific growth condition was analyzed using a Bioconductor software package edgeR with a p-value cut off at 0.05 .
Identification and characterization of the targets of CodY
Information on the CodY regulons of various streptococcal species was retrieved from RegPrecise version 3.0 database . A position weight matrix (PWM) of the 5′ flanking regions of the CodY regulon was constructed with the MEME software at the default parameter settings. The PWM model was used to search for potential CodY binding motif in the 400-bp region 5′ to the translation start site of all genes in S. salivarius 57.I by using MAST with a p value less than 10− 4 . The transcriptomes of wild-type 57.I and strain ΔcodY under each growth condition were used to seek targets of CodY by using the following criteria. Firstly, a significant difference in the expression level was observed between strains in at least one growth condition, and the regulation pattern (repression or activation) was identical in all conditions. Secondly, the regulation by CodY was consistent in genes of one operon. Thirdly, genes that are differentially expressed between strains and are with a CodY binding consensus were considered direct targets; genes without a CodY binding consensus were considered indirect targets.
The growth kinetics of S. salivarius strains was analyzed using the Bioscreen C Microbiology Reader (Oy Growth Curve AB Ltd.) as described . Overnight cultures were diluted at 1:20 in brain heart infusion (BHI) and grown to an optical density at 600 nm (OD600) of 0.8. The culture was sonicated at 10% amplitude twice, 5 s each, diluted at 1:1000 in BHI and then 300 μl of the diluted suspension was inoculated in a microtiter plate. The OD600 value was monitored at 10-min intervals for 10 h. The plate was shaken for 10 s before each reading.
The static and the flow-cell biofilm systems
The static biofilm formation of S. salivarius strains was examined by the described method  with modifications . Overnight cultures in biofilm medium (BM)  containing 20 mM glucose (BMG) were diluted at 1:20 in BMG and grown to an OD600 of 0.8. The culture was sonicated as described above, followed by dilution at 1:100 in BMG. The diluted suspension was then inoculated in the wells of a 96-well microtiter plate. The plate was incubated at 37 °C in 10% CO2 for 24 h. The biofilm was stained with crystal violet and the absorbance at 562 (OD562) nm was determined with a microplate reader (SoftMax®Pro, Molecular Devices).
The biofilm grown in a flow-cell system was examined as described . Overnight cultures in BMG were diluted at 1:20 and grown to an OD600 of 0.8. The cultures were diluted at 1:800 in BMG and 300 μl of the diluted cell suspension was injected into the flow cell chamber. The chamber was kept at 37 °C in an upside down position without medium flow for 4 h. The medium flow rate was set up at 5 ml h− 1 channel− 1 and the biofilm was grown for 10 h. The biofilm was stained with SYTO 9/propidium iodide (PI) reagent (Invitrogen) and observed by using confocal laser scanning microscope (CLSM).
The pH drop assay
The glycolytic activity of S. salivarius 57.I and its derivatives was evaluated by the acid drop analysis . Cultures at an OD600 of 0.6 in BHI were harvested, washed once with 50 mM KCl-1 mM MgCl2, and concentrated in the same solution with an OD600 of 1.0. The cell suspension was equilibrated at pH 7.2 with 0.01 N KOH at room temperature, and the pH drop was initiated with the addition of 55.5 mM glucose. The pH value was recorded every 30 s for 1 h.
Assays for acid resistance and oxidative stress tolerance
The growth of S. salivarius 57.I and it derivatives at pH 5.5 was evaluated using the described method  with modifications. Overnight cultures were diluted at 1:20 in BHI and grown to an OD600 of 0.8. S. salivarius ΔcodY was subject to sonication three times prior to reading. The cultures were serially diluted and 10 μl of each dilution was spotted on BHI agar (pH 7.2) or BHI agar that has been adjusted to pH 5.5 by the addition of 2 N HCl (BHI-HCl). The plates were incubated at 37 °C in 10% CO2 for 24 h.
A similar approach was use to evaluate oxidative stress tolerance. Cultures were prepared as described above, serially diluted and spotted onto BHI agar plates supplemented with 0.5, 1, 2 or 3 mM H2O2. Plates were incubated at 37 °C in 10% CO2 for 24 h.
Galleria mellonella larva model
The toxicity of S. salivarius strains in the G. mellonella larva model was determined using the described method  with modifications. Bacterial cells were grown to an OD600 of 0.8 in BHI, harvested, washed twice with PBS, and suspended in PBS at 1.0 × 109 CFU ml− 1. 10 μl of the suspension was injected into the proleg of each larva with a 26 s ga needle (Hamilton). PBS and heat inactivated (75 °C for 30 min) bacterial suspensions with the same CFU were used as negative controls. For each bacterial strain, 12 larvae were infected. The infected larvae were kept at 37 °C without feeding and recorded daily until the death of all larvae. The significant difference between infection groups was analyzed by using Kaplan Meier curve with log-rank test.
To determine the CFU in a larva, the larva was surface-cleaned with 70% ethanol, placed in a 2-ml microfuge tube containing 1 ml PBS and 12 zirconia/silica beads (2.3 mm diameter). The larva was homogenized in a Mini-Beadbeater (Biospec Products) for a total of 40 s at room temperature. The bacteria count in the homogenate was determined by serial dilutions and plating.
Branched-chain amino acids
Brain heart infusion medium
Biofilm medium containing 20 mM glucose
Confocal laser scanning microscope
Position weight matrix
Reads per Kb per million mapped reads
Tryptone-yeast extract medium
Brinsmade SR. CodY, a master integrator of metabolism and virulence in gram-positive bacteria. Curr Genet. 2017;63:417–25.
Stenz L, Francois P, Whiteson K, Wolz C, Linder P, Schrenzel J. The CodY pleiotropic repressor controls virulence in gram-positive pathogens. FEMS Immunol Med Microbiol. 2011;62:123–39.
Sonenshein AL. CodY, a global regulator of stationary phase and virulence in gram-positive bacteria. Curr Opin Microbiol. 2005;8:203–7.
Levdikov VM, Blagova E, Joseph P, Sonenshein AL, Wilkinson AJ. The structure of CodY, a GTP- and isoleucine-responsive regulator of stationary phase and virulence in gram-positive bacteria. J Biol Chem. 2006;281:11366–73.
Joseph P, Ratnayake-Lecamwasam M, Sonenshein AL. A region of Bacillus subtilis CodY protein required for interaction with DNA. J Bacteriol. 2005;187:4127–39.
Handke LD, Shivers RP, Sonenshein AL. Interaction of Bacillus subtilis CodY with GTP. J Bacteriol. 2008;190:798–806.
Ratnayake-Lecamwasam M, Serror P, Wong KW, Sonenshein AL. Bacillus subtilis CodY represses early-stationary-phase genes by sensing GTP levels. Genes Dev. 2001;15:1093–103.
Shivers RP, Sonenshein AL. Activation of the Bacillus subtilis global regulator CodY by direct interaction with branched-chain amino acids. Mol Microbiol. 2004;53:599–611.
Villapakkam AC, Handke LD, Belitsky BR, Levdikov VM, Wilkinson AJ, Sonenshein AL. Genetic and biochemical analysis of the interaction of Bacillus subtilis CodY with branched-chain amino acids. J Bacteriol. 2009;191:6865–76.
Guedon E, Serror P, Ehrlich SD, Renault P, Delorme C. Pleiotropic transcriptional repressor CodY senses the intracellular pool of branched-chain amino acids in Lactococcus lactis. Mol Microbiol. 2001;40:1227–39.
Lemos JA, Nascimento MM, Lin VK, Abranches J, Burne RA. Global regulation by (p)ppGpp and CodY in Streptococcus mutans. J Bacteriol. 2008;190:5291–9.
Hendriksen WT, Bootsma HJ, Estevao S, Hoogenboezem T, de Jong A, de Groot R, Kuipers OP, Hermans PW. CodY of Streptococcus pneumoniae: link between nutritional gene regulation and colonization. J Bacteriol. 2008;190:590–601.
Han AR, Kang HR, Son J, Kwon DH, Kim S, Lee WC, Song HK, Song MJ, Hwang KY. The structure of the pleiotropic transcription regulator CodY provides insight into its GTP-sensing mechanism. Nucleic Acids Res. 2016;44:9483–93.
Belitsky BR, Sonenshein AL. Genetic and biochemical analysis of CodY-binding sites in Bacillus subtilis. J Bacteriol. 2008;190:1224–36.
den Hengst CD, van Hijum SA, Geurts JM, Nauta A, Kok J, Kuipers OP. The Lactococcus lactis CodY regulon: identification of a conserved cis-regulatory element. J Biol Chem. 2005;280:34332–42.
Lobel L, Herskovits AA. Systems level analyses reveal multiple regulatory activities of CodY controlling metabolism, motility and virulence in Listeria monocytogenes. PLoS Genet. 2016;12:e1005870.
Slamti L, Lemy C, Henry C, Guillot A, Huillet E, Lereclus D. CodY regulates the activity of the virulence quorum sensor PlcR by controlling the import of the signaling peptide PapR in Bacillus thuringiensis. Front Microbiol. 2015;6:1501.
Golub LM, Borden SM, Kleinberg I. Urea content of gingival crevicular fluid and its relation to periodontal diseases in humans. J Periodontal Res. 1971;6:243–51.
Burne RA, Marquis RE. Alkali production by oral bacteria and protection against dental caries. FEMS Microbiol Lett. 2000;193:1–6.
Chen YY, Burne RA. Analysis of Streptococcus salivarius urease expression using continuous chemostat culture. FEMS Microbiol Lett. 1996;135:223–9.
Chen YY, Weaver CA, Burne RA. Dual functions of Streptococcus salivarius urease. J Bacteriol. 2000;182:4667–9.
Huang SC, Burne RA, Chen YY. The pH-dependent expression of the urease operon in Streptococcus salivarius is mediated by CodY. Appl Environ Microbiol. 2014;80:5386–93.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Alves LA, Nomura R, Mariano FS, Harth-Chu EN, Stipp RN, Nakano K, Mattos-Graner RO. CovR regulates Streptococcus mutans susceptibility to complement immunity and survival in blood. Infect Immun. 2016;84:3206–19.
Churchward G. The two faces of Janus: virulence gene regulation by CovR/S in group a streptococci. Mol Microbiol. 2007;64:34–41.
Landwehr-Kenzel S, Henneke P. Interaction of Streptococcus agalactiae and cellular innateiImmunity in colonization and disease. Front Immunol. 2014;5:519.
Federle MJ, Scott JR. Identification of binding sites for the group a streptococcal global regulator CovR. Mol Microbiol. 2002;43:1161–72.
Kuhnert WL, Zheng G, Faustoferri RC, Quivey RG Jr. The F-ATPase operon promoter of Streptococcus mutans is transcriptionally regulated in response to external pH. J Bacteriol. 2004;186:8524–8.
Lu WW, Wang Y, Wang T, Kong J. The global regulator CodY in Streptococcus thermophilus controls the metabolic network for escalating growth in the milk environment. Appl Environ Microbiol. 2015;81:2349–58.
Huang SC, Chen YY. Role of VicRKX and GlnR in pH-dependent regulation of the Streptococcus salivarius 57.I urease operon. mSphere. 2016;1:e00033–16.
Kreth J, Chen Z, Ferretti J, Malke H. Counteractive balancing of transcriptome expression involving CodY and CovRS in Streptococcus pyogenes. J Bacteriol. 2011;193:4153–65.
Yoshimura G, Komatsuzawa H, Hayashi I, Fujiwara T, Yamada S, Nakano Y, Tomita Y, Kozai K, Sugai M. Identification and molecular characterization of an N-acetylmuraminidase, Aml, involved in Streptococcus mutans cell separation. Microbiol Immunol. 2006;50:729–42.
Bose JL, Lehman MK, Fey PD, Bayles KW. Contribution of the Staphylococcus aureus Atl AM and GL murein hydrolase activities in cell division, autolysis, and biofilm formation. PLoS One. 2012;7:e42244.
Layec S, Gerard J, Legue V, Chapot-Chartier MP, Courtin P, Borges F, Decaris B, Leblond-Bourget N. The CHAP domain of Cse functions as an endopeptidase that acts at mature septa to promote Streptococcus thermophilus cell separation. Mol Microbiol. 2009;71:1205–17.
Ahn SJ, Rice KC, Oleas J, Bayles KW, Burne RA. The Streptococcus mutans Cid and Lrg systems modulate virulence traits in response to multiple environmental signals. Microbiol. 2010;156:3136–47.
Sharma-Kuinkel BK, Mann EE, Ahn JS, Kuechenmeister LJ, Dunman PM, Bayles KW. The Staphylococcus aureus LytSR two-component regulatory system affects biofilm formation. J Bacteriol. 2009;191:4767–75.
Kim JN, Burne RA. CcpA and CodY coordinate acetate metabolism in Streptococcus mutans. Appl Environ Microbiol. 2017;83:e03274–16.
Santiago B, MacGilvray M, Faustoferri RC, Quivey RG Jr. The branched-chain amino acid aminotransferase encoded by ilvE is involved in acid tolerance in Streptococcus mutans. J Bacteriol. 2012;194:2010–9.
Santiago B, Marek M, Faustoferri RC, Quivey RG Jr. The Streptococcus mutans aminotransferase encoded by ilvE is regulated by CodY and CcpA. J Bacteriol. 2013;195:3552–62.
Fozo EM, Quivey RG Jr. The fabM gene product of Streptococcus mutans is responsible for the synthesis of monounsaturated fatty acids and is necessary for survival at low pH. J Bacteriol. 2004;186:4152–8.
Fozo EM, Quivey RG Jr. Shifts in the membrane fatty acid profile of Streptococcus mutans enhance survival in acidic environments. Appl Environ Microbiol. 2004;70:929–36.
Wang Y, He HY, Li HH, Lu WW, Guo TT, Kong J. The global regulator CodY responds to oxidative stress by the regulation of glutathione biosynthesis in Streptococcus thermophilus. J Dairy Sci. 2017;100:8768–75.
Hajaj B, Yesilkaya H, Shafeeq S, Zhi X, Benisty R, Tchalah S, Kuipers OP, Porat N. CodY regulates thiol peroxidase expression as part of the pneumococcal defense mechanism against H2O2 stress. Front Cell Infect Microbiol. 2017;7:210.
VanderWal AR, Makthal N, Pinochet-Barros A, Helmann JD, Olsen RJ, Kumaraswami M. Iron efflux by PmtA is critical for oxidative stress resistance and contributes significantly to group a Streptococcus virulence. Infect Immun. 2017;85:e00091–17.
Ricci S, Janulczyk R, Bjorck L. The regulator PerR is involved in oxidative stress response and iron homeostasis and is necessary for full virulence of Streptococcus pyogenes. Infect Immun. 2002;70:4968–76.
Yamaguchi M, Terao Y, Kawabata S. Pleiotropic virulence factor - Streptococcus pyogenes fibronectin-binding proteins. Cell Microbiol. 2013;15:503–11.
Foster TJ. The remarkably multifunctional fibronectin binding proteins of Staphylococcus aureus. Eur J Clin Microbiol Infect Dis. 2016;35:1923–31.
Rubens CE, Heggen LM. Tn916 delta E: a Tn916 transposon derivative expressing erythromycin resistance. Plasmid. 1988;20:137–42.
Lau PC, Sung CK, Lee JH, Morrison DA, Cvitkovitch DG. PCR ligation mutagenesis in transformable streptococci: application and efficiency. J Microbiol Methods. 2002;49:193–205.
Perez-Casal J, Caparon MG, Scott JR. Mry, a trans-acting positive regulator of the M protein gene of Streptococcus pyogenes with similarity to the receptor proteins of two-component regulatory systems. J Bacteriol. 1991;173:2617–24.
Fontaine L, Boutry C, de Frahan MH, Delplace B, Fremaux C, Horvath P, Boyaval P, Hols P. A novel pheromone quorum-sensing system controls the development of natural competence in Streptococcus thermophilus and Streptococcus salivarius. J Bacteriol. 2010;192:1444–54.
Sissons CH, Hancock EM, Perinpanayagam HE, Cutress TW. The bacteria responsible for ureolysis in artificial dental plaque. Arch Oral Biol. 1988;33:727–33.
Chen YY, Weaver CA, Mendelsohn DR, Burne RA. Transcriptional regulation of the Streptococcus salivarius 57.I urease operon. J Bacteriol. 1998;180:5769–75.
Novichkov PS, Kazakov AE, Ravcheev DA, Leyn SA, Kovaleva GY, Sutormin RA, Kazanov MD, Riehl W, Arkin AP, Dubchak I, Rodionov DA. RegPrecise 3.0--a resource for genome-scale exploration of transcriptional regulation in bacteria. BMC Genomics. 2013;14:745.
Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, Ren J, Li WW, Noble WS. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–8.
Chen YY, Chen YY, Hung JL, Chen PM, Chia JS. The GlnR regulon in Streptococcus mutans is differentially regulated by GlnR and PmrA. PLoS One. 2016;11:e0159599.
O'Toole GA, Kolter R. Initiation of biofilm formation in Pseudomonas fluorescens WCS365 proceeds via multiple, convergent signalling pathways: a genetic analysis. Mol Microbiol. 1998;28:449–61.
Loo CY, Corliss DA, Ganeshkumar N. Streptococcus gordonii biofilm formation: identification of genes that code for biofilm phenotypes. J Bacteriol. 2000;182:1374–82.
Belli WA, Marquis RE. Adaptation of Streptococcus mutans and Enterococcus hirae to acid stress in continuous culture. Appl Environ Microbiol. 1991;57:1134–8.
Abranches J, Miller JH, Martinez AR, Simpson-Haidaris PJ, Burne RA, Lemos JA. The collagen-binding protein Cnm is required for Streptococcus mutans adherence to and intracellular invasion of human coronary artery endothelial cells. Infect Immun. 2011;79:2277–84.
We are grateful to Zong-Xian Ye for providing the G. mellonella system and Chiao-Wen Wu for the assistance with the growth study. We thank Simon Silver and Le Phung for reviewing this manuscript.
This work was supported by the Chang Gung Memorial Hospital (CMRPD1G0151) to YMC and Ministry of Science and Technology of ROC (105–2320-B-182-015-MY3) to YMC. The funding body had no role in study design, collection, analysis, and interpretation of the data, preparation of the manuscript or decision to publish.
Availability of data and materials
The transcriptome data are in the additional files.
Ethics approval and consent to participate
This study has not directly involved humans, animals, plants or cell lines. The source of all bacterial strains is listed in Additional file 2: Table S4.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. The CodY binding consensus sequences in Gram positive bacteria. (A) The binding consensus derived from streptococcal species. A search for the consensus of the CodY binding motifs from 12 streptococcal species was conducted using the MEME suites. (B) An alignment of the CodY binding consensus from Streptococcus spp., Lactococcus lactis, Bacillus spp., C. difficile and S. aureus. N, A or C or G or T; R, A or G; W, A or T; Y, C or T. (PPTX 59 kb)
Table S1. The RPKM values of the direct targets of CodY in S. salivarius 57.I and ΔcodY. Table S2. The RPKM values of the indirect targets of CodY in S. salivarius 57.I and ΔcodY. Table S3. The RPKM values of genes encoding enzymes of the EMP pathway in S. salivarius 57.I and ΔcodY. Table S4. Bacterial strains used in this study. Table S5. Oligonucleotides used in this study. (DOCX 30 kb)
About this article
Cite this article
Geng, J., Huang, S., Chen, Y. et al. Impact of growth pH and glucose concentrations on the CodY regulatory network in Streptococcus salivarius. BMC Genomics 19, 386 (2018). https://doi.org/10.1186/s12864-018-4781-z
- Streptococcus salivarius
- Growth pH
- Carbohydrate availability
- Stress response