Permalink
Switch branches/tags
Nothing to show
Find file
Fetching contributors…
Cannot retrieve contributors at this time
1 lines (1 sloc) 139 KB
<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">01501</article-id><article-id pub-id-type="doi">10.7554/eLife.01501</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Microbiology and infectious disease</subject></subj-group></article-categories><title-group><article-title>A serine sensor for multicellularity in a bacterium</article-title></title-group><contrib-group><contrib contrib-type="author" id="author-7695"><name><surname>Subramaniam</surname><given-names>Arvind R</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="other" rid="par-4"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf3"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-7696"><name><surname>DeLoughery</surname><given-names>Aaron</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf3"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-7697"><name><surname>Bradshaw</surname><given-names>Niels</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf3"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-7702"><name><surname>Chen</surname><given-names>Yun</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf3"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-7699"><name><surname>O’Shea</surname><given-names>Erin</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="aff" rid="aff3"/><xref ref-type="aff" rid="aff4"/><xref ref-type="other" rid="par-2"/><xref ref-type="fn" rid="con7"/><xref ref-type="fn" rid="conf1"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-1019"><name><surname>Losick</surname><given-names>Richard</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-7701"><name><surname>Chai</surname><given-names>Yunrong</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff5"/><xref ref-type="corresp" rid="cor2">*</xref><xref ref-type="other" rid="par-3"/><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf3"/><xref ref-type="other" rid="dataro1"/></contrib><aff id="aff1"><institution content-type="dept">Department of Molecular and Cellular Biology</institution>, <institution>Harvard University</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff2"><institution content-type="dept">Faculty of Arts and Sciences Center for Systems Biology</institution>, <institution>Harvard University</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff3"><institution content-type="dept">Department of Chemistry and Chemical Biology</institution>, <institution>Harvard University</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff4"><institution>Howard Hughes Medical Institute, Harvard University</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff5"><institution content-type="dept">Department of Biology</institution>, <institution>Northeastern University</institution>, <addr-line><named-content content-type="city">Boston</named-content></addr-line>, <country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Storz</surname><given-names>Gisela</given-names></name><role>Reviewing editor</role><aff><institution>National Institute of Child Health and Human Development</institution>, <country>United States</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>losick@mcb.harvard.edu</email> (RL);</corresp><corresp id="cor2"><label>*</label>For correspondence: <email>y.chai@neu.edu</email> (YC)</corresp></author-notes><pub-date date-type="pub" publication-format="electronic"><day>17</day><month>12</month><year>2013</year></pub-date><pub-date pub-type="collection"><year>2013</year></pub-date><volume>2</volume><elocation-id>e01501</elocation-id><history><date date-type="received"><day>09</day><month>09</month><year>2013</year></date><date date-type="accepted"><day>05</day><month>11</month><year>2013</year></date></history><permissions><copyright-statement>© 2013, Subramaniam et al</copyright-statement><copyright-year>2013</copyright-year><copyright-holder>Subramaniam et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife01501.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.01501.001</object-id><p>We report the discovery of a simple environmental sensing mechanism for biofilm formation in the bacterium <italic>Bacillus subtilis</italic> that operates without the involvement of a dedicated RNA or protein. Certain serine codons, the four TCN codons, in the gene for the biofilm repressor SinR caused a lowering of SinR levels under biofilm-inducing conditions. Synonymous substitutions of these TCN codons with AGC or AGT impaired biofilm formation and gene expression. Conversely, switching AGC or AGT to TCN codons upregulated biofilm formation. Genome-wide ribosome profiling showed that ribosome density was higher at UCN codons than at AGC or AGU during biofilm formation. Serine starvation recapitulated the effect of biofilm-inducing conditions on ribosome occupancy and SinR production. As serine is one of the first amino acids to be exhausted at the end of exponential phase growth, reduced translation speed at serine codons may be exploited by other microbes in adapting to stationary phase.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.001">http://dx.doi.org/10.7554/eLife.01501.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.01501.002</object-id><title>eLife digest</title><p>Bacteria use several different mechanisms to recognize and respond to changes in their environment. Protein sensors, for example, relay signals from the cell surface to target molecules within the cell. Additionally, RNA sensors can respond to internal levels of chemicals by regulating the expression of genes. Now, Subramaniam et al. have discovered a sensing mechanism in bacteria that does not rely on either protein sensors or RNA sensors.</p><p>Bacteria normally live as free swimming cells in water. But sometimes, in response to certain environmental conditions, they form a multicellular community called a biofilm. In this biofilm, bacteria encase themselves with layers of carbohydrates and proteins, which then protect the bacteria from adverse chemicals such as antibiotics.</p><p>A protein known as SinR plays a key role during biofilm formation in the model bacterium <italic>Bacillus subtilis</italic>. SinR normally prevents the formation of the carbohydrates and proteins that make up the biofilm. Upon the decision to form a biofilm, <italic>B. subtilis</italic> counters the effect of SinR by producing an anti-SinR protein called SinI.</p><p>Now Subramaniam et al. have found that as well as producing the SinI protein, <italic>B. subtilis</italic> use an additional mechanism to promote biofilm formation. This mechanism relies on codons, the elements within genes that correspond to specific amino acids. Six different codons correspond to the amino acid serine, and the gene for the SinR protein contains above average numbers of four of them. These four codons are highly sensitive to serine levels, and they decrease the levels of the SinR protein when there is less serine in the environment, as happens to be the case in biofilms. And since SinR prevents the production of the carbohydrates and proteins that make up the biofilms, a decrease in the levels of SinR leads to an increase in the production of biofilms.</p><p>Since the serine codons at the heart of the sensing mechanism discovered by Subramaniam et al. are present in all forms of life, from viruses to humans, it is possible that similar sensing mechanisms might be found in contexts other than bacterial biofilms, such as in viral infection and cancer.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.002">http://dx.doi.org/10.7554/eLife.01501.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>codon usage</kwd><kwd>ribosome pausing</kwd><kwd>starvation</kwd><kwd>biofilm</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd><italic>B. subtilis</italic></kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>GM18568</award-id><principal-award-recipient><name><surname>Losick</surname><given-names>Richard</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution>Howard Hughes Medical Institute</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>O’Shea</surname><given-names>Erin</given-names></name></principal-award-recipient></award-group><award-group id="par-3"><funding-source><institution-wrap><institution>Northeastern University</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Chai</surname><given-names>Yunrong</given-names></name></principal-award-recipient></award-group><award-group id="par-4"><funding-source><institution-wrap><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>1K99 GM107113</award-id><principal-award-recipient><name><surname>Subramaniam</surname><given-names>Arvind R</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta><meta-name>elife-xml-version</meta-name><meta-value>2</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>Biofilm formation by Bacillus subtilis relies on a sensing mechanism for the amino acid serine that does not depend on a dedicated protein or RNA.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Bacteria constantly monitor their environment and internal physiological state so that they can adapt to changing conditions. A wide variety of sensing mechanisms are deployed for this purpose, including dedicated protein sensors, such as histidine kinases, which mediate changes in gene expression by controlling the phosphorylation of cognate response regulators in response to environmental cues (<xref ref-type="bibr" rid="bib49">West and Stock, 2001</xref>). Bacteria also sense changes in their environment and physiology by means of dedicated RNAs, such as the highly structured, leader RNA for the tryptophan operon, which controls the transcription of downstream genes in the operon by a mechanism involving ribosome stalling at tryptophan codons (<xref ref-type="bibr" rid="bib18">Henkin and Yanofsky, 2002</xref>). Here we report the discovery of an unusually simple mechanism of environmental sensing involved in the process of biofilm formation by the bacterium <italic>B. subtilis</italic> that does not require a dedicated RNA or protein.</p><p>Biofilm formation involves a switch from planktonic growth as individual cells to the formation of complex, multicellular communities in response to environmental cues (<xref ref-type="bibr" rid="bib26">Kolter and Greenberg, 2006</xref>). In <italic>B. subtilis</italic>, these communities are embedded in a self-produced matrix consisting of polysaccharide and an amyloid-like protein, which are specified by the <italic>epsA-O</italic> and the <italic>tapA-sipW-tasA</italic> operons, respectively (<xref ref-type="bibr" rid="bib4">Branda et al., 2001</xref>; <xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>). The transition to multicellularity is governed in part by four histidine kinases (KinA, KinB, KinC and KinD) that control the phosphorylation of the response regulator, Spo0A, a master regulator of post-exponential phase gene expression (<xref ref-type="fig" rid="fig1">Figure 1A</xref>) (<xref ref-type="bibr" rid="bib22">Jiang et al., 2000</xref>; <xref ref-type="bibr" rid="bib46">Vlamakis et al., 2013</xref>). Recent studies suggest that KinA and KinB respond to impaired respiration (<xref ref-type="bibr" rid="bib25">Kolodkin-Gal et al., 2013</xref>), whereas KinC responds to membrane perturbations and KinD to unknown chemical signals (<xref ref-type="bibr" rid="bib31">López et al., 2009</xref>; <xref ref-type="bibr" rid="bib43">Shemesh et al., 2010</xref>; <xref ref-type="bibr" rid="bib10">Chen et al., 2012</xref>; <xref ref-type="bibr" rid="bib2">Beauregard et al., 2013</xref>). Once phosphorylated, Spo0A turns on <italic>sinI</italic>, a gene encoding a small protein antagonist of the biofilm-specific regulatory protein SinR (<xref ref-type="bibr" rid="bib34">Molle et al., 2003</xref>; <xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>). SinR, which is produced constitutively, is a repressor of the matrix operons, <italic>epsA-O</italic> and the <italic>tapA-sipW-tasA</italic>, as well as other biofilm-related genes (<xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>; <xref ref-type="bibr" rid="bib11">Chu et al., 2006</xref>; <xref ref-type="bibr" rid="bib6">Chai et al., 2009</xref>). SinR is also a repressor of the gene for SlrR (<xref ref-type="bibr" rid="bib12">Chu et al., 2008</xref>), which together with SinR sets up a self-reinforcing, double-negative feedback loop for matrix gene expression (<xref ref-type="fig" rid="fig1">Figure 1A</xref>) (<xref ref-type="bibr" rid="bib7">Chai et al., 2010</xref>; <xref ref-type="bibr" rid="bib35">Norman et al., 2013</xref>). A special feature of SinR of relevance to this investigation is that the expression of matrix genes is hypersensitive to small perturbations in the level of the protein (<xref ref-type="bibr" rid="bib8">Chai et al., 2011</xref>). This hypersensitivity is attributed to molecular titration of SinR by SinI and cooperativity among SinR molecules bound to tandem target sequences at regulatory sites for the matrix operons (<xref ref-type="bibr" rid="bib8">Chai et al., 2011</xref>; <xref ref-type="bibr" rid="bib5">Chai et al., 2008</xref>).<fig-group><fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.01501.003</object-id><label>Figure 1.</label><caption><title>Switching synonymous serine codons in <italic>sinR</italic> affects biofilm formation.</title><p>(<bold>A</bold>) Regulatory circuit controlling biofilm formation in <italic>B. subtilis</italic>. (<bold>B</bold>) Top: Serine codon usage in the <italic>sinR</italic> coding sequence. Number within parenthesis indicates the frequency of the corresponding codon in <italic>sinR</italic>. Bottom: Average serine codon usage across 4153 protein-coding sequences in the <italic>B. subtilis</italic> genome. Number within parenthesis indicates the relative frequency of each codon in the genome. (<bold>C</bold>) Colony morphology for the wild-type strain and the indicated <italic>sinR</italic> synonymous variants grown on solid biofilm-inducing medium. Three TCA codons in the wild-type sequence of <italic>sinR</italic> were switched to each of the other five serine codons. The wild-type (WT) <italic>sinR</italic> sequence was replaced by the <italic>sinR</italic> synonymous mutant at the native <italic>sinR</italic> locus of the strain 3610. (<bold>D</bold>) SinR protein level during entry into biofilm formation (OD<sub>600</sub> = 2) measured using an anti-SinR antibody that also cross-reacts with SlrR, a protein that is 85% identical to SinR. Western blot against the RNA polymerase subunit SigA was used as the loading control. Whole cell lysates were loaded at different dilutions (indicated as X, X/2, and X/3). (<bold>E</bold>) Densitometry of SinR bands in (<bold>D</bold>) after normalization by SigA. (<bold>F</bold>) Top panel: Western blot against SinR and SlrR using anti-SinR antibody. Bottom panel: Densitometry ratio of the SlrR and SinR bands in the top panel. Error bars represent standard error over three replicate Western blots. The SlrR/SinR ratio for each blot was normalized such that the wild-type strain had a ratio of 1. (<bold>G</bold>) Matrix gene expression monitored using a P<sub><italic>epsA</italic></sub>–<italic>lacZ</italic> transcriptional reporter inserted at the chromosomal <italic>amyE</italic> locus. <italic>β</italic>-galactosidase activity was measured at OD<sub>600</sub> = 2 in liquid biofilm-inducing medium. Error bars represent standard error of three measurements.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.003">http://dx.doi.org/10.7554/eLife.01501.003</ext-link></p></caption><graphic xlink:href="elife01501f001"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01501.004</object-id><label>Figure 1—figure supplement 1.</label><caption><title><italic>sinR</italic> coding sequence.</title><p>The three TCA codons (switched in <xref ref-type="fig" rid="fig1">Figure 1</xref>) are highlighted in red. The three TCC codons and the two AGC/AGT codons (switched in <xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>) are highlighted in green and blue respectively. The remaining serine codons are shown in yellow.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.004">http://dx.doi.org/10.7554/eLife.01501.004</ext-link></p></caption><graphic xlink:href="elife01501fs001"/></fig><fig id="fig1s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01501.005</object-id><label>Figure 1—figure supplement 2.</label><caption><title>Effect of TCC and AGC/AGT synonymous substitutions in the <italic>sinR</italic> gene on colony morphology and biofilm reporter activity.</title><p>(<bold>A</bold>) Colony morphology for the indicated <italic>sinR</italic> synonymous variants grown on solid biofilm-inducing medium. Either three TCC codons or two AGC/AGT (AGY) codons in the wild-type sequence of <italic>sinR</italic> were switched to remaining serine synonymous codons. The wild-type (WT) <italic>sinR</italic> sequence was replaced by the <italic>sinR</italic> synonymous variant at the native <italic>sinR</italic> locus of the strain 3610. Colony morphology of the wild-type strain is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>. (<bold>B</bold> and <bold>C</bold>) Matrix gene expression monitored using a P<sub><italic>epsA</italic></sub>–<italic>lacZ</italic> transcriptional reporter inserted at the chromosomal <italic>amyE</italic> locus. Strains were grown in liquid biofilm-inducing medium and <italic>β</italic>-galactosidase activity was measured at an OD<sub>600</sub> = 2. Error bars represent standard error of three measurements. The synonymous variants highlighted in red do not follow the hierarchy between TCN and AGC/AGT codons seen for the six TCA synonymous variants in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.005">http://dx.doi.org/10.7554/eLife.01501.005</ext-link></p></caption><graphic xlink:href="elife01501fs002"/></fig></fig-group></p><p>In the current study, we present evidence for the existence of a novel cellular sensing mechanism controlling biofilm formation. Rather than relying on regulation by a dedicated RNA or protein, the translation speed of ribosomes decreases at certain serine codons, resulting in lower SinR levels, which as a consequence, contributes to derepression of SinR-controlled genes. We propose that specific serine codons in the <italic>sinR</italic> mRNA act as a simple sensor for monitoring, and triggering a response to, serine depletion under biofilm-inducing conditions.</p></sec><sec id="s2" sec-type="results"><title>Results</title><sec id="s2-1"><title>Switching synonymous serine codons in <italic>sinR</italic> affects biofilm formation</title><p>In a genetic screen to identify suppressor mutations that rescued the biofilm-defective phenotype of a <italic>B. subtilis</italic> mutant (‘Materials and methods’), we unexpectedly recovered a variant that contained a ‘silent’ mutation that resulted in a switch from one serine codon to a synonymous codon in <italic>sinR</italic>. This observation prompted us to ask whether switching other serine codons might also influence biofilm formation. Serine is specified by six codons: AGC, AGT, TCA, TCC, TCG and TCT (where T is U in the mRNA). We noticed that the <italic>sinR</italic> coding sequence has a slightly higher frequency of the four TCN serine codons as compared to the average frequency of TCN codons in the <italic>B. subtilis</italic> genome (<xref ref-type="fig" rid="fig1">Figure 1B</xref>, p=0.22, N = 12). To test whether this bias towards TCN codons has an effect on biofilm formation, we systematically replaced the three TCA codons in <italic>sinR</italic> (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>) with each of the other five serine codons. Replacing the TCA codons with AGT or AGC resulted in flat, featureless colonies on solid biofilm-inducing medium, indicating severely impaired biofilm formation (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). In contrast, replacing the TCA codons with either TCC or TCT had little or no effect on colony morphology whereas switching to TCG increased the wrinkled appearance of the colonies (<xref ref-type="fig" rid="fig1">Figure 1C</xref>), which is indicative of robust biofilm formation.</p><p>Next, we asked whether switching serine codons was altering the level of SinR protein in liquid biofilm-inducing medium. Immunoblot analysis with anti-SinR antibodies revealed slightly yet consistently higher SinR levels in the strain with the AGT variant of <italic>sinR</italic> when compared to either the wild-type strain with three TCA codons or the TCG variant of <italic>sinR</italic> (<xref ref-type="fig" rid="fig1">Figure 1D,E</xref>).</p><p>SinR is highly similar (85% identity) to SlrR, which also plays a critical role in biofilm formation and whose gene (<italic>slrR</italic>) is under the direct negative control of SinR (<xref ref-type="bibr" rid="bib12">Chu et al., 2008</xref>). Because SlrR cross reacts with the anti-SinR antibodies, we were also able to detect SlrR in our immunoblot analysis. Strikingly, the levels of SlrR were almost perfectly anti-correlated with those of SinR, with the differences in the SlrR protein levels among the <italic>sinR</italic> synonymous variants being much higher than the corresponding differences in SinR protein levels (<xref ref-type="fig" rid="fig1">Figure 1F</xref>). Because repression by SinR is ultrasensitive to SinR levels (<xref ref-type="bibr" rid="bib8">Chai et al., 2011</xref>), small differences in SinR protein levels among <italic>sinR</italic> synonymous variants might be sufficient to cause large differences in the levels of expression of SinR-repressed genes such as <italic>slrR</italic>. Consistent with this idea, an <italic>eps-lacZ</italic> transcriptional fusion reporter for the SinR-repressed <italic>epsA-O</italic> matrix operon showed that the four TCN <italic>sinR</italic> variants had 3- to 19-fold higher <italic>β</italic>-galactosidase activity than the AGT and AGC variants (<xref ref-type="fig" rid="fig1">Figure 1G</xref>).</p><p>To test the generality of the observed hierarchy between the synonymous variants of <italic>sinR</italic>, we generated an additional set of eleven <italic>sinR</italic> synonymous variants in which we replaced either three TCC codons or two AGC/AGT codons (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1</xref>) with their synonymous counterparts. Eight of these variants conformed to the hierarchy described above, namely, the four TCN variants behaved oppositely to the two AGC/AGT variants in colony morphology and in <italic>eps-lacZ</italic> reporter expression (<xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>). The three variants that did not conform to the hierarchy could potentially reflect alterations to the mRNA sequence context near the mutation rather than the effect of a synonymous substitution per se. Taken together, the above results suggest that serine synonymous codons in the <italic>sinR</italic> coding sequence have a stereotypical effect on biofilm formation that is primarily determined by the differential usage of the four TCN and the two AGC/AGT codons.</p></sec><sec id="s2-2"><title>Entry into biofilm formation is accompanied by codon-specific increase in ribosome density</title><p>What is the mechanism by which serine codon usage affects SinR protein levels and biofilm formation? Synonymous codon changes can alter the synthesis of the encoded protein through changes in the translation initiation rate, mRNA levels or the ribosome elongation rate (<xref ref-type="bibr" rid="bib37">Plotkin and Kudla, 2010</xref>). However, the effect of synonymous codon usage on the initiation rate and mRNA levels is context-specific; only codons near the AUG start site affect translation initiation (<xref ref-type="bibr" rid="bib27">Kudla et al., 2009</xref>), whereas only codons that are located in certain regions of secondary structure or at ribonuclease cleavage sites affect mRNA levels (<xref ref-type="bibr" rid="bib3">Bernstein et al., 2002</xref>). Our observation that synonymous codon replacements at multiple locations along <italic>sinR</italic> have a stereotypical effect on biofilm formation argues against such context-specific mechanisms (except for the three exceptional cases noted above).</p><p>To test the alternative hypothesis that serine codon usage might alter the ribosome elongation rate, and given that the ribosome elongation rate at a codon varies inversely with the average ribosome density at that codon, we measured ribosome density on mRNAs at single-codon resolution using the ribosome profiling method (<xref ref-type="bibr" rid="bib19">Ingolia et al., 2009</xref>; <xref ref-type="bibr" rid="bib36">Oh et al., 2011</xref>; <xref ref-type="bibr" rid="bib21">Ingolia et al., 2012</xref>). We grew <italic>B. subtilis</italic> in liquid biofilm-inducing medium, harvested cells either during exponential phase growth (OD<sub>600</sub> = 0<italic>.</italic>6) or during stationary phase when biofilm formation is induced (OD<sub>600</sub> = 1<italic>.</italic>4), and performed deep-sequencing of ribosome protected mRNA fragments and size-matched total mRNA fragments. Ribosome profiling yielded 3.75 and 2.55 million sequencing reads aligning to annotated protein-coding sequences for the exponential phase sample and the biofilm entry sample, respectively. The number of reads aligning to a single codon on individual mRNAs was too low for accurate quantification of ribosome density, and was not sufficient to directly detect increased ribosome density on the <italic>sinR</italic> transcript. However, we reasoned that the global pattern of ribosome density at codons across all mRNAs should reflect the ribosome density on individual transcripts such as <italic>sinR</italic>. Therefore, we calculated the median ribosome density at each of the 61 sense codons across the 1556 protein-coding sequences in the exponential phase sample and the 1148 sequences in the biofilm entry sample that had an average coverage greater than one sequencing read per codon (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>, ‘Materials and methods’). This analysis reproduced the previous observation (<xref ref-type="bibr" rid="bib29">Li et al., 2012</xref>) of increased ribosome density 8 to 11 nt downstream of Shine-Dalgarno-like trinucleotide sequences both during exponential phase and during biofilm entry (<xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2</xref>).</p><p>During exponential phase, median ribosome density varied over a twofold range, with no systematic difference between serine codons and the remaining codons (<xref ref-type="fig" rid="fig2">Figure 2A</xref>). By contrast, during biofilm entry, median ribosome density was significantly higher at serine and cysteine codons as compared to the remaining codons (<xref ref-type="fig" rid="fig2">Figure 2B</xref>), suggesting that the translation speed of ribosomes is selectively reduced at these codons during biofilm entry. Notably, the ribosome density at serine codons was not uniform: the four UCN codons had 1.9 to 2.1-fold higher ribosome density whereas the AGC and AGU codons had only 1.1 and 1.3-fold higher ribosome density respectively, relative to the median value across 61 sense codons. Further, this difference in ribosome density between UCN codons and the AGC/AGU codons was essentially identical when computed separately for codons located in the first half or in the second half of each gene (<xref ref-type="fig" rid="fig2s3">Figure 2—figure supplement 3</xref>), a finding that underscores the statistical robustness and the context independence of the observed difference.<fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.01501.006</object-id><label>Figure 2.</label><caption><title>Entry into biofilm formation is accompanied by codon-specific increase in ribosome density.</title><p>Genome-wide median ribosome density and total mRNA density at 61 sense codons (excluding start and stop codons) (<bold>A</bold>) during exponential phase growth (OD<sub>600</sub> = 0.6), and (<bold>B</bold>) during stationary phase when biofilm formation is induced (OD<sub>600</sub> = 1.4). The six serine (red) and two cysteine (green) codons are highlighted. Genome-wide ribosome density and total mRNA density were measured by deep-sequencing of ribosome-protected mRNA fragments and total mRNA fragments respectively, of a <italic>B. subtilis</italic> 3610 derivative (Δ<italic>epsH</italic>) grown in liquid biofilm-inducing medium.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.006">http://dx.doi.org/10.7554/eLife.01501.006</ext-link></p></caption><graphic xlink:href="elife01501f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01501.007</object-id><label>Figure 2—figure supplement 1.</label><caption><title>Computational workflow for deep-sequencing data analysis.</title><p>All steps outlined here were performed in Bash and Python programming languages. For further details on individual steps, see ‘Materials and methods’.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.007">http://dx.doi.org/10.7554/eLife.01501.007</ext-link></p></caption><graphic xlink:href="elife01501fs003"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01501.008</object-id><label>Figure 2—figure supplement 2.</label><caption><title>Increase in ribosome density downstream of Shine-Dalgarno-like trinucleotide sequences.</title><p>Median ribosome density across all protein coding sequences was computed for the 60 nt region around each of six Shine-Dalgarno-like trinucleotide sequences (<xref ref-type="bibr" rid="bib29">Li et al., 2012</xref>) for the exponential phase sample (left-hand panel) and the biofilm entry sample (right-hand panel).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.008">http://dx.doi.org/10.7554/eLife.01501.008</ext-link></p></caption><graphic xlink:href="elife01501fs004"/></fig><fig id="fig2s3" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01501.009</object-id><label>Figure 2—figure supplement 3.</label><caption><title>Context independence of ribosome and mRNA densities during biofilm formation.</title><p>Each gene was conceptually divided into two equal halves and the ribosome density and mRNA density was computed separately for codons located either in the first half (left-hand panel) or in the second half (right-hand panel) of each gene. All other analysis steps were identical to those in <xref ref-type="fig" rid="fig2">Figure 2</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.009">http://dx.doi.org/10.7554/eLife.01501.009</ext-link></p></caption><graphic xlink:href="elife01501fs005"/></fig></fig-group></p></sec><sec id="s2-3"><title>Serine starvation reduces translation speed and inhibits SinR synthesis in a codon-specific manner</title><p>Because ribosomes density increased at serine and cysteine codons only during biofilm formation and not during exponential phase growth, we hypothesized that the increased ribosome density was caused by the depletion of intracellular pools of these two amino acids during biofilm entry rather than by any intrinsic feature of the mRNA (<xref ref-type="bibr" rid="bib29">Li et al., 2012</xref>) or the nascent polypeptide (<xref ref-type="bibr" rid="bib9">Charneski and Hurst, 2013</xref>). Our hypothesis is also supported by the previous observation that synonymous codon usage can have a starvation-specific effect on protein levels (<xref ref-type="bibr" rid="bib44">Subramaniam et al., 2013</xref>). Further, serine is a precursor metabolite for the biosynthesis of cysteine (<xref ref-type="bibr" rid="bib16">Gagnon et al., 1994</xref>); hence cysteine depletion was likely the result of a decrease in intracellular serine concentration. Consistent with this hypothesis, we also observed an increase in ribosome density at both serine and cysteine codons during serine starvation of a <italic>B. subtilis</italic> serine-auxotrophic mutant (<xref ref-type="fig" rid="fig3">Figure 3C</xref>). Importantly, serine starvation resulted in an increase in ribosome density at only the four UCN serine codons but not at the AGC and AGU codons, matching the hierarchy seen during biofilm entry (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). Serine starvation also resulted in differential levels of production of SinR-YFP protein fusions bearing different synonymous serine codons, whereas serine-rich growth resulted in identical levels of fusion protein production from these variants (<xref ref-type="fig" rid="fig3">Figure 3A,B</xref>). Finally, the addition of excess serine or cysteine (but not any of the other 18 amino acids) blocked biofilm formation in wild type cells as judged after 48 hr of growth in biofilm-inducing medium (<xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref>, and data not shown).<fig-group><fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.01501.010</object-id><label>Figure 3.</label><caption><title>Serine starvation reduces translation speed and inhibits SinR synthesis in a codon-specific manner.</title><p>(<bold>A</bold> and <bold>B</bold>) Three <italic>sinR</italic> synonymous variants were synthesized with 10 serine codons switched to AGC, TCA or TCG. The variants were expressed as SinR-YFP fusions from the <italic>amyE</italic> locus under the control of a <italic>lac</italic> promoter in a 3610-Δ<italic>serA</italic> serine auxotroph strain growing in serine-limited medium. Black arrow around 300 min indicates the onset of serine starvation caused by depletion of exogenously-added serine in the growth medium. Cell density (<bold>A</bold>) and the corresponding SinR-YFP protein level (<bold>B</bold>) were monitored using a 96-well spectrophotometer. (<bold>C</bold>) Genome-wide median ribosome density for 61 sense codons (excluding start and stop codons) during serine starvation (vertical axis) and serine-rich growth (horizontal axis) of a serine auxotrophic strain. (<bold>D</bold>) Fold-change in average ribosome density for individual genes upon biofilm entry (vertical axis) or serine starvation (horizontal axis). Genes that were preferentially up-regulated at least 10-fold upon biofilm entry in comparison to serine starvation are highlighted in red (68 genes, <xref ref-type="table" rid="tbl1">Table 1</xref>). Only genes with a minimum of 100 ribosome profiling reads in at least one of the samples were included in this analysis (1724 genes) and the reported log2 fold-changes are median-subtracted values across this gene set.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.010">http://dx.doi.org/10.7554/eLife.01501.010</ext-link></p></caption><graphic xlink:href="elife01501f003"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01501.011</object-id><label>Figure 3—figure supplement 1.</label><caption><title>Addition of excess serine or cysteine blocks pellicle formation by <italic>B. subtilis</italic>.</title><p>Single amino acids were added at 300 µg ml<sup>−1</sup> to liquid MSgg medium. Biofilm formation of 3610 was assayed visually by pellicle formation at the air-liquid interface 48 hr after inoculation. Serine and cysteine were found to block pellicle formation out of all 20 amino acids tested.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.011">http://dx.doi.org/10.7554/eLife.01501.011</ext-link></p></caption><graphic xlink:href="elife01501fs006"/></fig></fig-group></p><p>Given that both biofilm entry and serine starvation resulted in increased ribosome density at serine and cysteine codons, we asked whether these two apparently different conditions invoke the same gene expression program in <italic>B. subtilis</italic>. The fold-change in average ribosome density of <italic>B. subtilis</italic> genes was positively correlated between biofilm entry and serine starvation (<xref ref-type="fig" rid="fig3">Figure 3D</xref>, <italic>R</italic><sup><italic>2</italic></sup> = 0<italic>.</italic>27<italic>,</italic> p=10<sup>−5</sup>, 1724 genes). However, a subset of 68 genes was induced at least 10-fold higher upon biofilm entry than during serine starvation (indicated by red markers in <xref ref-type="fig" rid="fig3">Figure 3D</xref>, <xref ref-type="table" rid="tbl1">Table 1</xref>). This subset included genes for anaerobic metabolism such as <italic>lctEP</italic>, <italic>nasDE</italic>, and <italic>cydAB</italic>, and is consistent with the recently proposed role of impaired respiration in biofilm formation (<xref ref-type="bibr" rid="bib25">Kolodkin-Gal et al., 2013</xref>). We observed that sulfur metabolism genes were also enriched in this subset, possibly indicating a stronger response to cysteine depletion during biofilm entry than during serine starvation.<table-wrap id="tbl1" position="float"><object-id pub-id-type="doi">10.7554/eLife.01501.012</object-id><label>Table 1.</label><caption><p><italic>B. subtilis</italic> genes that have greater than 10-fold difference in expression ratio between biofilm formation and serine starvation</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.012">http://dx.doi.org/10.7554/eLife.01501.012</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><th>Gene</th><th>A</th><th>B</th><th>C</th><th>D</th><th>E</th><th>F</th><th>Function</th></tr></thead><tbody><tr><td><italic>albA</italic></td><td align="char" char=".">9.35</td><td>−0.69</td><td align="char" char=".">32</td><td align="char" char=".">10129</td><td align="char" char=".">175</td><td align="char" char=".">144</td><td>antilisterial bacteriocin (subtilosin) production protein</td></tr><tr><td><italic>alsD</italic></td><td align="char" char=".">6.38</td><td>0.95</td><td align="char" char=".">34</td><td align="char" char=".">1384</td><td align="char" char=".">53</td><td align="char" char=".">135</td><td>alpha-acetolactate decarboxylase</td></tr><tr><td><italic>alsS</italic></td><td align="char" char=".">6.96</td><td>1.79</td><td align="char" char=".">76</td><td align="char" char=".">4613</td><td align="char" char=".">86</td><td align="char" char=".">396</td><td>acetolactate synthase</td></tr><tr><td><italic>cah</italic></td><td align="char" char=".">2.08</td><td>−2.74</td><td align="char" char=".">2531</td><td align="char" char=".">5254</td><td align="char" char=".">1486</td><td align="char" char=".">295</td><td>S-deacylase</td></tr><tr><td><italic>ctc</italic></td><td align="char" char=".">3.78</td><td>−0.73</td><td align="char" char=".">327</td><td align="char" char=".">2206</td><td align="char" char=".">322</td><td align="char" char=".">257</td><td>50S ribosomal protein L25</td></tr><tr><td><italic>cydA</italic></td><td align="char" char=".">5.36</td><td>−3.18</td><td align="char" char=".">234</td><td align="char" char=".">4720</td><td align="char" char=".">408</td><td align="char" char=".">60</td><td>cytochrome bd ubiquinol oxidase subunit I</td></tr><tr><td><italic>cydB</italic></td><td align="char" char=".">9.6</td><td>−1.88</td><td align="char" char=".">9</td><td align="char" char=".">3374</td><td align="char" char=".">129</td><td align="char" char=".">46</td><td>cytochrome bd ubiquinol oxidase subunit II</td></tr><tr><td><italic>cysK</italic></td><td align="char" char=".">1.52</td><td>−2.21</td><td align="char" char=".">15760</td><td align="char" char=".">22120</td><td align="char" char=".">8991</td><td align="char" char=".">2580</td><td>cysteine synthase</td></tr><tr><td><italic>gcvPA</italic></td><td align="char" char=".">0.97</td><td>−2.38</td><td align="char" char=".">1500</td><td align="char" char=".">1440</td><td align="char" char=".">1001</td><td align="char" char=".">255</td><td>glycine dehydrogenase subunit 1</td></tr><tr><td><italic>gcvPB</italic></td><td align="char" char=".">1.05</td><td>−2.28</td><td align="char" char=".">2029</td><td align="char" char=".">2057</td><td align="char" char=".">1224</td><td align="char" char=".">335</td><td>glycine dehydrogenase subunit 2</td></tr><tr><td><italic>gspA</italic></td><td align="char" char=".">4.54</td><td>−0.75</td><td align="char" char=".">111</td><td align="char" char=".">1261</td><td align="char" char=".">135</td><td align="char" char=".">106</td><td>glycosyl transferase (general stress protein)</td></tr><tr><td><italic>iseA</italic></td><td align="char" char=".">2.24</td><td>−1.14</td><td align="char" char=".">1117</td><td align="char" char=".">2579</td><td align="char" char=".">2521</td><td align="char" char=".">1516</td><td>inhibitor of cell-separation enzymes</td></tr><tr><td><italic>lctP</italic></td><td align="char" char=".">6.62</td><td>−2.9</td><td align="char" char=".">62</td><td align="char" char=".">3009</td><td align="char" char=".">247</td><td align="char" char=".">44</td><td>L-lactate permease</td></tr><tr><td><italic>ldh</italic></td><td align="char" char=".">3.47</td><td>−3.88</td><td align="char" char=".">2810</td><td align="char" char=".">15,254</td><td align="char" char=".">1732</td><td align="char" char=".">157</td><td>L-lactate dehydrogenase</td></tr><tr><td><italic>maeN</italic></td><td align="char" char=".">4.08</td><td>−1.61</td><td align="char" char=".">148</td><td align="char" char=".">1226</td><td align="char" char=".">347</td><td align="char" char=".">151</td><td>Na+/malate symporter</td></tr><tr><td><italic>mccA</italic></td><td align="char" char=".">3.23</td><td>−2.72</td><td align="char" char=".">759</td><td align="char" char=".">3494</td><td align="char" char=".">386</td><td align="char" char=".">78</td><td>cystathionine beta-synthase</td></tr><tr><td><italic>metE</italic></td><td align="char" char=".">2.72</td><td>−2.13</td><td align="char" char=".">33,027</td><td align="char" char=".">106473</td><td align="char" char=".">32,231</td><td align="char" char=".">9760</td><td>5-methyltetrahydropteroyltriglutamate/homocysteine S-methyltransferase</td></tr><tr><td><italic>mgsR</italic></td><td align="char" char=".">4.72</td><td>−1.53</td><td align="char" char=".">163</td><td align="char" char=".">2099</td><td align="char" char=".">202</td><td align="char" char=".">93</td><td>stress transcriptional regulator</td></tr><tr><td><italic>mtlA</italic></td><td align="char" char=".">3.92</td><td>−0.22</td><td align="char" char=".">126</td><td align="char" char=".">934</td><td align="char" char=".">104</td><td align="char" char=".">119</td><td>PTS system mannitol-specific transporter subunit IICB</td></tr><tr><td><italic>mtnA</italic></td><td align="char" char=".">1.84</td><td>−2.79</td><td align="char" char=".">2522</td><td align="char" char=".">4425</td><td align="char" char=".">3501</td><td align="char" char=".">671</td><td>methylthioribose-1-phosphate isomerase</td></tr><tr><td><italic>mtnD</italic></td><td align="char" char=".">1.72</td><td>−1.63</td><td align="char" char=".">4095</td><td align="char" char=".">6600</td><td align="char" char=".">3182</td><td align="char" char=".">1361</td><td>acireductone dioxygenase</td></tr><tr><td><italic>mtnK</italic></td><td align="char" char=".">2.54</td><td>−2.64</td><td align="char" char=".">4290</td><td align="char" char=".">12,196</td><td align="char" char=".">5866</td><td align="char" char=".">1250</td><td>methylthioribose kinase</td></tr><tr><td><italic>nasD</italic></td><td align="char" char=".">6.21</td><td>−0.65</td><td align="char" char=".">217</td><td align="char" char=".">7872</td><td align="char" char=".">635</td><td align="char" char=".">536</td><td>assimilatory nitrite reductase subunit</td></tr><tr><td><italic>nasE</italic></td><td align="char" char=".">6.05</td><td>0.38</td><td align="char" char=".">34</td><td align="char" char=".">1106</td><td align="char" char=".">63</td><td align="char" char=".">109</td><td>assimilatory nitrite reductase subunit</td></tr><tr><td><italic>rbfK</italic></td><td align="char" char=".">3.74</td><td>−2.57</td><td align="char" char=".">1329</td><td align="char" char=".">8681</td><td align="char" char=".">434</td><td align="char" char=".">97</td><td>RNA-binding riboflavin kinase</td></tr><tr><td><italic>sboA</italic></td><td align="char" char=".">7.77</td><td>0.64</td><td align="char" char=".">121</td><td align="char" char=".">12,895</td><td align="char" char=".">256</td><td align="char" char=".">530</td><td>subtilosin A</td></tr><tr><td><italic>sboX</italic></td><td align="char" char=".">8.26</td><td>0.19</td><td align="char" char=".">21</td><td align="char" char=".">3128</td><td align="char" char=".">76</td><td align="char" char=".">115</td><td>bacteriocin-like product</td></tr><tr><td><italic>ssuA</italic></td><td align="char" char=".">2.09</td><td>−2.3</td><td align="char" char=".">3691</td><td align="char" char=".">7682</td><td align="char" char=".">877</td><td align="char" char=".">236</td><td>aliphatic sulfonate ABC transporter binding lipoprotein</td></tr><tr><td><italic>ssuB</italic></td><td align="char" char=".">2.57</td><td>−2.09</td><td align="char" char=".">2416</td><td align="char" char=".">7004</td><td align="char" char=".">779</td><td align="char" char=".">242</td><td>aliphatic sulfonate ABC transporter ATP-binding protein</td></tr><tr><td><italic>ssuC</italic></td><td align="char" char=".">2.17</td><td>−2.16</td><td align="char" char=".">3784</td><td align="char" char=".">8362</td><td align="char" char=".">692</td><td align="char" char=".">206</td><td>aliphatic sulfonate ABC transporter permease</td></tr><tr><td><italic>ssuD</italic></td><td align="char" char=".">2.2</td><td>−2.19</td><td align="char" char=".">12,961</td><td align="char" char=".">29,135</td><td align="char" char=".">2812</td><td align="char" char=".">817</td><td>alkanesulfonate monooxygenase</td></tr><tr><td><italic>tcyJ</italic></td><td align="char" char=".">3.25</td><td>−3.26</td><td align="char" char=".">1046</td><td align="char" char=".">4856</td><td align="char" char=".">427</td><td align="char" char=".">59</td><td>sulfur-containing amino acid ABC transporter binding lipoprotein</td></tr><tr><td><italic>tcyK</italic></td><td align="char" char=".">3.79</td><td>−3.55</td><td align="char" char=".">2815</td><td align="char" char=".">19,087</td><td align="char" char=".">1095</td><td align="char" char=".">124</td><td>sulfur-containing amino acid ABC transporter binding lipoprotein</td></tr><tr><td><italic>tcyL</italic></td><td align="char" char=".">3.33</td><td>−3.05</td><td align="char" char=".">855</td><td align="char" char=".">4223</td><td align="char" char=".">387</td><td align="char" char=".">62</td><td>sulfur-containing amino acid ABC transporter permease</td></tr><tr><td><italic>tcyM</italic></td><td align="char" char=".">3.54</td><td>−2.97</td><td align="char" char=".">1859</td><td align="char" char=".">10,556</td><td align="char" char=".">581</td><td align="char" char=".">99</td><td>sulfur-containing amino acid ABC transporter permease</td></tr><tr><td><italic>tcyN</italic></td><td align="char" char=".">3.38</td><td>−2.74</td><td align="char" char=".">3363</td><td align="char" char=".">17,149</td><td align="char" char=".">1119</td><td align="char" char=".">223</td><td>sulfur-containing amino acid ABC transporter ATP-binding protein</td></tr><tr><td><italic>ureA</italic></td><td align="char" char=".">4.36</td><td>0.78</td><td align="char" char=".">124</td><td align="char" char=".">1252</td><td align="char" char=".">77</td><td align="char" char=".">176</td><td>urease subunit gamma</td></tr><tr><td><italic>ycgL</italic></td><td align="char" char=".">0.99</td><td>−3.01</td><td align="char" char=".">512</td><td align="char" char=".">500</td><td align="char" char=".">278</td><td align="char" char=".">46</td><td>hypothetical protein</td></tr><tr><td><italic>ycgM</italic></td><td align="char" char=".">2.46</td><td>−1.35</td><td align="char" char=".">39</td><td align="char" char=".">105</td><td align="char" char=".">178</td><td align="char" char=".">93</td><td>proline oxidase</td></tr><tr><td><italic>ycgN</italic></td><td align="char" char=".">2.45</td><td>−1.6</td><td align="char" char=".">928</td><td align="char" char=".">2475</td><td align="char" char=".">1267</td><td align="char" char=".">556</td><td>1-pyrroline-5-carboxylate dehydrogenase</td></tr><tr><td><italic>ycnJ</italic></td><td align="char" char=".">0.75</td><td>−2.63</td><td align="char" char=".">168</td><td align="char" char=".">138</td><td align="char" char=".">121</td><td align="char" char=".">26</td><td>copper import protein</td></tr><tr><td><italic>ydaG</italic></td><td align="char" char=".">4.14</td><td>0.49</td><td align="char" char=".">70</td><td align="char" char=".">604</td><td align="char" char=".">80</td><td align="char" char=".">150</td><td>general stress protein</td></tr><tr><td><italic>ydbL</italic></td><td align="char" char=".">3.07</td><td>−0.31</td><td align="char" char=".">294</td><td align="char" char=".">1210</td><td align="char" char=".">130</td><td align="char" char=".">139</td><td>hypothetical protein</td></tr><tr><td><italic>yeaA</italic></td><td align="char" char=".">1.34</td><td>−2.59</td><td align="char" char=".">112</td><td align="char" char=".">138</td><td align="char" char=".">139</td><td align="char" char=".">31</td><td>hypothetical protein</td></tr><tr><td><italic>yezD</italic></td><td align="char" char=".">2.52</td><td>−4.17</td><td align="char" char=".">145</td><td align="char" char=".">406</td><td align="char" char=".">340</td><td align="char" char=".">25</td><td>hypothetical protein</td></tr><tr><td><italic>yitJ</italic></td><td align="char" char=".">3.01</td><td>−2.43</td><td align="char" char=".">2785</td><td align="char" char=".">10,990</td><td align="char" char=".">4687</td><td align="char" char=".">1153</td><td>bifunctional homocysteine S-methyltransferase/5,10-methylenetetrahydrofolate reductase</td></tr><tr><td><italic>yjbC</italic></td><td align="char" char=".">3.67</td><td>−0.57</td><td align="char" char=".">104</td><td align="char" char=".">647</td><td align="char" char=".">248</td><td align="char" char=".">222</td><td>thiol oxidation management factor; acetyltransferase</td></tr><tr><td><italic>yjnA</italic></td><td align="char" char=".">1.5</td><td>−2.81</td><td align="char" char=".">972</td><td align="char" char=".">1342</td><td align="char" char=".">790</td><td align="char" char=".">149</td><td>hypothetical protein</td></tr><tr><td><italic>yoaB</italic></td><td align="char" char=".">2.28</td><td>−2.26</td><td align="char" char=".">2206</td><td align="char" char=".">5244</td><td align="char" char=".">2862</td><td align="char" char=".">792</td><td>negatively charged metabolite transporter</td></tr><tr><td><italic>yoaC</italic></td><td align="char" char=".">2.92</td><td>−1.89</td><td align="char" char=".">1200</td><td align="char" char=".">4459</td><td align="char" char=".">1436</td><td align="char" char=".">513</td><td>hydroxylated metabolite kinase</td></tr><tr><td><italic>yrhB</italic></td><td align="char" char=".">2.92</td><td>−2.93</td><td align="char" char=".">3806</td><td align="char" char=".">14,096</td><td align="char" char=".">1904</td><td align="char" char=".">333</td><td>cystathionine beta-lyase</td></tr><tr><td><italic>yrrT</italic></td><td align="char" char=".">2.97</td><td>−3.12</td><td align="char" char=".">546</td><td align="char" char=".">2089</td><td align="char" char=".">443</td><td align="char" char=".">68</td><td>AdoMet-dependent methyltransferase</td></tr><tr><td><italic>ytlI</italic></td><td align="char" char=".">1.66</td><td>−3.2</td><td align="char" char=".">206</td><td align="char" char=".">318</td><td align="char" char=".">141</td><td align="char" char=".">20</td><td>LysR family transcriptional regulator</td></tr><tr><td><italic>ytmI</italic></td><td align="char" char=".">3.34</td><td>−3.27</td><td align="char" char=".">3452</td><td align="char" char=".">17,173</td><td align="char" char=".">1640</td><td align="char" char=".">226</td><td>N-acetyltransferase</td></tr><tr><td><italic>ytmO</italic></td><td align="char" char=".">3.4</td><td>−2.88</td><td align="char" char=".">3866</td><td align="char" char=".">20,007</td><td align="char" char=".">1179</td><td align="char" char=".">213</td><td>monooxygenase</td></tr><tr><td><italic>ytnI</italic></td><td align="char" char=".">3</td><td>−2.6</td><td align="char" char=".">3522</td><td align="char" char=".">13,770</td><td align="char" char=".">867</td><td align="char" char=".">189</td><td>redoxin</td></tr><tr><td><italic>ytnJ</italic></td><td align="char" char=".">3.14</td><td>−2.86</td><td align="char" char=".">10,645</td><td align="char" char=".">45,997</td><td align="char" char=".">2495</td><td align="char" char=".">456</td><td>monooxygenase</td></tr><tr><td><italic>ytnL</italic></td><td align="char" char=".">3.56</td><td>−2.45</td><td align="char" char=".">1281</td><td align="char" char=".">7371</td><td align="char" char=".">354</td><td align="char" char=".">86</td><td>aminohydrolase</td></tr><tr><td><italic>ytnM</italic></td><td align="char" char=".">3.5</td><td>−2.6</td><td align="char" char=".">4542</td><td align="char" char=".">25,202</td><td align="char" char=".">1264</td><td align="char" char=".">277</td><td>transporter</td></tr><tr><td><italic>yuaF</italic></td><td align="char" char=".">1.89</td><td>−1.74</td><td align="char" char=".">87</td><td align="char" char=".">157</td><td align="char" char=".">158</td><td align="char" char=".">63</td><td>membrane integrity integral inner membrane protein</td></tr><tr><td><italic>yvzB</italic></td><td align="char" char=".">0.95</td><td>−2.48</td><td align="char" char=".">125</td><td align="char" char=".">118</td><td align="char" char=".">168</td><td align="char" char=".">40</td><td>flagellin</td></tr><tr><td><italic>yxaL</italic></td><td align="char" char=".">3.54</td><td>−0.41</td><td align="char" char=".">1069</td><td align="char" char=".">6077</td><td align="char" char=".">442</td><td align="char" char=".">442</td><td>membrane associated protein kinase</td></tr><tr><td><italic>yxbB</italic></td><td align="char" char=".">3.7</td><td>−0.01</td><td align="char" char=".">108</td><td align="char" char=".">685</td><td align="char" char=".">118</td><td align="char" char=".">155</td><td>S-adenosylmethionine-dependent methyltransferase</td></tr><tr><td><italic>yxeK</italic></td><td align="char" char=".">0.86</td><td>−2.72</td><td align="char" char=".">2702</td><td align="char" char=".">2406</td><td align="char" char=".">1073</td><td align="char" char=".">216</td><td>monooxygenase</td></tr><tr><td><italic>yxeL</italic></td><td align="char" char=".">1.29</td><td>−2.94</td><td align="char" char=".">437</td><td align="char" char=".">525</td><td align="char" char=".">202</td><td align="char" char=".">35</td><td>acetyltransferase</td></tr><tr><td><italic>yxeM</italic></td><td align="char" char=".">0.87</td><td>−2.57</td><td align="char" char=".">3003</td><td align="char" char=".">2692</td><td align="char" char=".">1047</td><td align="char" char=".">233</td><td>ABC transporter binding lipoprotein</td></tr><tr><td><italic>yxeP</italic></td><td align="char" char=".">1.75</td><td>−2.24</td><td align="char" char=".">2577</td><td align="char" char=".">4246</td><td align="char" char=".">736</td><td align="char" char=".">207</td><td>amidohydrolase</td></tr><tr><td><italic>yxjH</italic></td><td align="char" char=".">2.02</td><td>−1.82</td><td align="char" char=".">4162</td><td align="char" char=".">8243</td><td align="char" char=".">3811</td><td align="char" char=".">1432</td><td>methyl-tetrahydrofolate methyltransferase</td></tr></tbody></table><table-wrap-foot><fn><p>A—median-subtracted log2 fold-change: biofilm/exponential-phase, B—median-subtracted log2 fold-change: serine starvation/serine rich, C—raw counts: biofilm entry, D—raw counts: exponential phase, E—raw counts: serine rich, F—raw counts: serine starvation.</p></fn></table-wrap-foot></table-wrap></p><p><italic>In toto</italic>, the results with the serine auxotroph support the inference that ribosome stalling observed during biofilm formation is due to a drop in intracellular serine levels. Our efforts to measure intracellular serine levels directly during growth in minimal, biofilm-inducing medium (MSgg) have been unsuccessful. Hence, we cannot rule out the less likely possibilities that biofilm entry causes serine and cysteine to be sequestered away from protein synthesis or that aminoacylation rate of the corresponding tRNAs decreases without changes in the intracellular serine and cysteine pools.</p></sec><sec id="s2-4"><title>Serine codon bias in biofilm-regulated genes reflects their expression during serine starvation</title><p>Here we found that the translation speed of ribosomes decreases at UCN serine codons and thereby modulates production of SinR and entry into biofilm formation. Based on our genome-wide measurements of ribosome density, we expect that the translation speed of ribosomes should decrease during biofilm entry not only on the <italic>sinR</italic> mRNA, but also on other mRNAs that are enriched for any of the four UCN codons. For example, we found that the four TCN codons are over-represented in two nucleotide biosynthesis genes, <italic>pyrAA</italic> and <italic>purB</italic> (<xref ref-type="fig" rid="fig4">Figure 4A</xref>). These two genes are also transcriptionally down regulated during biofilm entry (<xref ref-type="fig" rid="fig4">Figure 4B</xref>). A high frequency of TCN codons in these genes might serve to reinforce their transcriptional down regulation by reducing translation speed. Conversely, TCN codons are under-represented in the genes encoding lactate dehydrogenase, <italic>ldh</italic> (<xref ref-type="bibr" rid="bib39">Cruz Ramos et al., 2000</xref>) and a master regulator of post-exponential phase gene expression, <italic>spo0A</italic> (<xref ref-type="bibr" rid="bib34">Molle et al., 2003</xref>) (<xref ref-type="fig" rid="fig4">Figure 4A</xref>). The two genes are transcriptionally up regulated upon biofilm entry (<xref ref-type="fig" rid="fig4">Figure 4B</xref>). The low frequency of TCN codons in these genes might represent a mechanism for optimizing the production of their protein products by minimizing the slowing down of translation during biofilm entry. Consistent with this idea, replacement of AGC/AGT codons by TCN codons in <italic>spo0A</italic>, whose protein product positively regulates biofilm entry (by turning on the synthesis of the SinR antagonist SinI), resulted in defective biofilm formation (<xref ref-type="fig" rid="fig4">Figure 4C</xref>) in contrast to the stimulatory effect of replacing AGC/AGT codons with TCN codons in <italic>sinR</italic> (<xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>).<fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.01501.013</object-id><label>Figure 4.</label><caption><title>Serine codon bias of biofilm-regulated genes reflects their expression under serine starvation.</title><p>(<bold>A</bold>) Relative serine codon fraction in genes for nucleotide biosynthesis (<italic>pyrAA</italic>, <italic>purB</italic>), lactate dehydrogenase (<italic>ldh</italic>) and a sporulation regulator (<italic>spo0A</italic>). Numbers in parentheses indicate the number of serine codons in each gene. Relative fraction of serine codons across the <italic>B. subtilis</italic> genome is shown for comparison. (<bold>B</bold>) Fold-change (expressed in log2 units) in average ribosome density upon biofilm entry for the four genes shown in <bold>A</bold>. (<bold>C</bold>) Colony morphology of a wild-type strain and two <italic>spo0A</italic> synonymous variants grown on solid biofilm-inducing medium. Seven AGC/AGT codons in wild-type <italic>spo0A</italic> were replaced by either 7 TCC codons or 3 TCC and 4 TCG codons and inserted at the chromosomal <italic>spo0A</italic> locus. Both the wild-type <italic>spo0A</italic> and the synonymous <italic>spo0A</italic> variants were inserted with a downstream selection marker. (<bold>D</bold>) Left: Codon Adaptation Index (CAI) for the four genes shown in <bold>A</bold>. Right: Distribution of CAI values for 4153 protein-coding sequences of <italic>B. subtilis</italic>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.013">http://dx.doi.org/10.7554/eLife.01501.013</ext-link></p></caption><graphic xlink:href="elife01501f004"/></fig></p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>Together, our results implicate serine depletion as an environmental cue that contributes to promoting biofilm formation in <italic>B. subtilis</italic> together with other cues that are sensed by the histidine kinases KinA–D. Serine depletion is sensed through a remarkably simple mechanism based on reduced translation speed at UCN serine codons in the mRNA for a regulatory protein, SinR, whose repressive effects are highly sensitive to small changes in the level of the protein. We presume that UCN codons lower SinR levels simply by slowing the rate of ribosome movement along the mRNA (elongation). However, it is possible that the reduced translation speed at UCN codons during biofilm entry could be followed by downstream events such as ribosome rescue (<xref ref-type="bibr" rid="bib24">Keiler et al., 1996</xref>) or mRNA decay (<xref ref-type="bibr" rid="bib17">Hayes and Sauer, 2003</xref>) that might also contribute to lowering the levels of the SinR protein.</p><p>The serine sensing mechanism uncovered here operates through over-representation of the TCN serine codons in the <italic>sinR</italic> gene without the necessity for any other dedicated protein or RNA for sensing serine depletion. By contrast, transcriptional attenuation, a widespread mechanism among bacteria for sensing amino acids that also relies on changes in translation speed, involves a translation-transcription coupling mechanism that is mediated by highly structured mRNAs and leader peptides (<xref ref-type="bibr" rid="bib18">Henkin and Yanofsky, 2002</xref>).</p><p>We note that the biased usage of the four TCN serine codons, which act as starvation sensors during biofilm formation, is not evident from widely-used phenomenological measures of codon bias such as the codon adaptation index (<xref ref-type="fig" rid="fig4">Figure 4D</xref>), which primarily reflects codon preferences during exponential growth (<xref ref-type="bibr" rid="bib42">Sharp and Li, 1987</xref>; <xref ref-type="bibr" rid="bib1">Andersson and Kurland, 1990</xref>). The difference in translation speed between the four UCN codons and the two AGC/AGU codons under biofilm-inducing conditions is likely mediated by differences in concentration of the corresponding aminoacylated tRNAs (<xref ref-type="bibr" rid="bib15">Elf et al., 2003</xref>; <xref ref-type="bibr" rid="bib13">Dittmar et al., 2005</xref>), as was recently observed in serine-starved <italic>E. coli</italic> (<xref ref-type="bibr" rid="bib44">Subramaniam et al., 2013</xref>). Interestingly, the hierarchy between UCN codons and AGC/AGU codons in <italic>B. subtilis</italic> during serine starvation is similar to the one in <italic>E. coli</italic> even though copy numbers of the corresponding tRNA genes have diverged significantly between these two organisms (<xref ref-type="bibr" rid="bib32">Lowe and Eddy, 1997</xref>). Despite different tRNA gene copy numbers, it is possible that the relative abundances of the serine tRNA isoacceptors are similar between the two organisms or that their relative abundances might be regulated in the same manner in response to nutrient deprivation (<xref ref-type="bibr" rid="bib14">Doi et al., 1966</xref>).</p><p>Serine is one of the first amino acids to be completely consumed from the culture medium when either <italic>B. subtilis</italic> or <italic>E. coli</italic> cells are grown in complex rich medium (<xref ref-type="bibr" rid="bib30">Liebs et al., 1988</xref>; <xref ref-type="bibr" rid="bib38">Prüss et al., 1994</xref>; <xref ref-type="bibr" rid="bib41">Sezanov et al., 2007</xref>). Indeed, increased ribosome density has been observed at serine codons during growth of <italic>E. coli</italic> in Luria-Bertani broth (<xref ref-type="bibr" rid="bib29">Li et al., 2012</xref>). Thus the role of synonymous serine codons as starvation sensors discovered here in the specific context of biofilm formation might be a general regulatory strategy in microbes for adapting to nutrient depletion at the end of exponential phase growth. It is noteworthy that depletion of specific amino acids affects developmental transitions in several eukaryotic cells (<xref ref-type="bibr" rid="bib33">Marin, 1976</xref>; <xref ref-type="bibr" rid="bib45">Sundrud et al., 2009</xref>; <xref ref-type="bibr" rid="bib48">Wang et al., 2009</xref>). It will be interesting to test whether a codon-based sensing mechanism, similar to the one found here in bacterial biofilm development, also plays a role in eukaryotic cells during amino acid depletion.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Bacterial strains and media</title><p>For ribosome profiling during biofilm formation, a 3610-Δ<italic>epsH</italic> strain (RL3852) was used to ensure dispersed growth in liquid media (<xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>). For serine starvation experiments, a serine-auxotrophic 3610-Δ<italic>serA</italic> strain (YC865) was used. A list of strains, plasmids, and oligonucleotides used in this work are summarized in <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>.</p><p>For general purposes, <italic>B. subtilis</italic> strains PY79, 3610, and their derivatives were grown in Luria-Bertani (LB) medium (10 g tryptone, 5 g yeast extract, and 5 g NaCl per liter broth) at 37°C. Antibiotics were added to the media at the following concentrations for <italic>B. subtilis</italic> strains: 10 µg ml<sup>−1</sup> of tetracycline, 100 µg ml<sup>−1</sup> of spectinomycin, 10 µg ml<sup>−1</sup> of kanamycin, 5 µg ml<sup>−1</sup> of chloramphenicol, and 1 µg ml<sup>−1</sup> of erythromycin. Minimal MSgg medium was used as the liquid biofilm-inducing medium. The same medium with 1.5% Bacto-agar (Difco, Franklin Lakes, NJ) was used as the solid biofilm-inducing medium. MSgg medium composition: 5 mM potassium phosphate (pH 7), 100 mM MOPS (pH 7), 2 mM MgCl<sub>2</sub>, 700 µM CaCl<sub>2</sub>, 50 µM MnCl<sub>2</sub>, 50 µM FeCl<sub>3</sub>, 1 µM ZnCl<sub>2</sub>, 2 µM thiamine, 0.5% glycerol, 0.5% glutamate, 50 µg ml<sup>−1</sup> tryptophan, 50 µg ml<sup>−1</sup> phenylalanine and 50 µg ml<sup>−1</sup> threonine. For overnight growth of serine auxotrophic 3610 strains, MSgg medium was supplemented with serine to a final concentration of 5 mM. For serine starvation experiments in which YFP fluorescence was measured (<xref ref-type="fig" rid="fig3">Figure 3A,B</xref>), MSgg medium was supplemented with 800 µM serine and 400 µM serine methyl-ester (Sigma, St. Louis, MO). Serine methyl-ester was added to ensure slow growth under serine starvation conditions.</p></sec><sec id="s4-2"><title>Strain construction</title><p>General methods for molecular cloning followed published protocols (<xref ref-type="bibr" rid="bib40">Sambrook 2001</xref>). SPP1 phage-mediated transduction was used to transfer antibiotic-marked DNA fragments between different strains (<xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>). Long-flanking PCR mutagenesis was applied to generate insertional deletion mutations (<xref ref-type="bibr" rid="bib47">Wach 1996</xref>). Synonymous switches in <italic>sinR</italic> and <italic>spo0A</italic> were generated by using synthetic DNA fragments (Genewiz, South Plainfield, NJ) or by applying site-directed mutagenesis (Roche, Switzerland). Sequences of the primers used in constructing mutant <italic>sinR</italic> alleles are described in <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>. Incorporation of synonymous substitutions into the <italic>sinR</italic> or <italic>spo0A</italic> gene at the native locus was done by allele exchange and followed a method described previously (<xref ref-type="bibr" rid="bib8">Chai et al., 2011</xref>).</p></sec><sec id="s4-3"><title>Biofilm assays</title><p><italic>B. subtilis</italic> cells were first grown in LB broth at 37°C to mid-exponential phase. For formation of biofilm colonies, 2 µl of the cells was spotted onto MSgg medium solidified with 1.5% agar. Plates were incubated at 23°C for 3–4 days before analysis. All images were taken using either a Nikon CoolPix 950 digital camera or using a SPOT camera (Diagnostic Instruments, Sterling Heights, MI). Assays for the <italic>β</italic>-galactosidase activities were described previously (<xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>).</p></sec><sec id="s4-4"><title>Genetic screen for suppressor mutants</title><p>Following a previously-published protocol (<xref ref-type="bibr" rid="bib7">Chai et al., 2010</xref>), we set up a genetic screen to search for spontaneous mutations that suppressed the defective biofilm phenotype of a <italic>B. subtilis ΔslrR</italic> mutant (YC131). The defective biofilm phenotype of the <italic>ΔslrR</italic> mutant is manifested as an inability to form robust floating pellicles (<xref ref-type="bibr" rid="bib12">Chu et al., 2008</xref>; <xref ref-type="bibr" rid="bib7">Chai et al., 2010</xref>). Briefly, the <italic>ΔslrR</italic> strain was inoculated into liquid MSgg medium in 6-well plates and incubated at 30°C. After 48 hr, pellicle formation was examined visually. In some wells, robust pellicles appeared possibly due to a second, suppressor mutation elsewhere in the genome. Cells from those wells were picked and streaked out on fresh LB agar plates to isolate single colonies. Cells from the single colonies were then tested for altered colony morphology on solid MSgg medium to confirm the suppressor phenotype. Similar genetic screens in previous studies (<xref ref-type="bibr" rid="bib23">Kearns et al., 2005</xref>; <xref ref-type="bibr" rid="bib7">Chai et al., 2010</xref>) had established that mutations in the <italic>sinR</italic> gene could suppress the defective biofilm phenotype of the <italic>ΔslrR</italic> mutant. Therefore, we isolated genomic DNA from the putative suppressor mutants, amplified the <italic>sinR</italic> gene by PCR, and then sequenced the <italic>sinR</italic> locus. Once a mutation in the <italic>sinR</italic> gene was confirmed by sequencing, the same mutation was reconstituted in the wild type background (3610) following a previous protocol (<xref ref-type="bibr" rid="bib7">Chai et al., 2010</xref>), and assayed for alteration in colony morphology on solid MSgg medium.</p></sec><sec id="s4-5"><title>Bacterial growth for ribosome profiling</title><p>For ribosome profiling during biofilm formation, fresh colonies were inoculated into 8 ml of MSgg liquid medium and grown for 12 hr at 30°C, 200 rpm. Saturated cultures were diluted 1:1000 into 200 ml aliquots of fresh MSgg medium and shaken in a 1L flask at 30°C, 200 rpm. For exponential-phase ribosome profiling (<xref ref-type="fig" rid="fig2">Figure 2A</xref>), cultures were harvested at OD<sub>600</sub> = 0.6. For ribosome profiling during biofilm entry (<xref ref-type="fig" rid="fig2">Figure 2B</xref>), cultures were harvested at OD<sub>600</sub> = 1.4. For the serine repletion experiment (<xref ref-type="fig" rid="fig4">Figure 4</xref>), serine was added to a final concentration of 2.5 mM at OD<sub>600</sub> = 1.4 and harvested after 30 min at 30°C, 200 rpm. For the serine starvation experiment (<xref ref-type="fig" rid="fig3">Figure 3B,C</xref>), pre-cultures were grown in MSgg medium supplemented with 5 mM serine and then diluted into 200 ml of the same medium. At an OD<sub>600</sub> = 0.6, the cultures were filtered and re-suspended either in MSgg medium (starvation) or in MSgg medium with 5 mM serine (control), and harvested after 60 min at 30°C, 200 rpm (<xref ref-type="fig" rid="fig3">Figure 3D</xref>).</p></sec><sec id="s4-6"><title>Western blotting</title><p>Cultures were grown with shaking at 37°C, and 14 ml culture aliquots were harvested at an OD<sub>600</sub> between 2.0 and 2.5. Cell pellets were collected by centrifugation and washed once with 10 ml of lysis buffer (20 mM Tris-HCl, 200 mM NaCl, 1 mM EDTA pH 7.4). Pellets were resuspended in 1.2 ml lysis buffer and incubated with 20 µg/ml of lysozyme (Sigma) for 1 hr on ice. Cells were further lysed by sonication. Cell debris was removed by centrifugation (14,000 rpm, 30 min, 4°C). The concentration of total protein in the lysates was determined by a Bradford assay (Bio-Rad, Hercules, CA). Samples for SDS-PAGE were prepared in Laemmli buffer normalized to equal protein concentration. Samples were ran on an NuPAGE 12% gel (Invitrogen, Carlsbad, CA, 1.0 mm, Bis/Tris, 200 V, ∼50 min) and transferred to a PVDF membrane (Millipore, Billerica, MA) at 100 V for 1 hr. The membrane was blocked in 5% milk-TPBS for 1 hr, and then incubated with anti-SinR antibody (1:2500, polyclonal) and anti-SigA (1:100,000, polyclonal) overnight. The membranes were washed 3 times in TPBS for 5 min each. Blots were incubated with goat anti-rabbit secondary antibody conjugated to Horseradish peroxidase (1:10,000, Bio-Rad). Blots were washed three times in TPBS, and developed with SuperSignal West Dura chemiluminescent substrate (Thermo, Waltham, MA) and imaged on a gel-doc (Bio-Rad). Densitometry analysis of Western blot images was performed using ImageJ software (NIH, <ext-link ext-link-type="uri" xlink:href="http://rsbweb.nih.gov/ij/">http://rsbweb.nih.gov/ij/</ext-link>). Rectangles were drawn around each distinct band and the average pixel intensity in this rectangle was calculated, followed by subtraction of background pixel intensity from an identical rectangle drawn in a region without any band. These background subtracted values were used for calculating SinR/SigA, SlrR/SigA and SinR / SlrR ratios in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p></sec><sec id="s4-7"><title>YFP fluorescence measurements</title><p>A previously published measurement protocol was used (<xref ref-type="bibr" rid="bib44">Subramaniam et al., 2013</xref>). Fresh colonies were inoculated into 1 ml of MSgg liquid medium with 5 mM serine and grown overnight in deep 96-well plates at 30°C, 1400 rpm. In the morning, saturated cultures were diluted 1:100 into 1 ml of MSgg medium with 800 µM serine and 400 µM serine methyl-ester and shaken at 30°C, 1350 rpm for 3 hr. Three aliquots of 150 µl from each culture was pipetted into three wells of three 96-well plates (3799, Costar, Corning, NY). Wallac Victor2 plate reader (PerkinElmer, Waltham, MA) was used to monitor cell density (absorbance at 600 nm) and YFP expression (fluorescence, excitation 504 nm and emission 540 nm). Each plate was read every 15 min using a robotic system (Caliper Life Sciences, Hopkinton, MA) and shaken at 1000 rpm in between readings (Variomag Teleshake shaker, Daytona Beach, FL). 30°C and 60% relative humidity was maintained throughout the experiment.</p></sec><sec id="s4-8"><title>Ribosome profiling and mRNA quantification</title><p>Ribosome profiling protocol was adapted from published literature (<xref ref-type="bibr" rid="bib19">Ingolia et al., 2009</xref>; <xref ref-type="bibr" rid="bib36">Oh et al., 2011</xref>; <xref ref-type="bibr" rid="bib21">Ingolia et al., 2012</xref>) with minor modifications as described below. Briefly, 200 ml of bacterial culture was harvested by filtration. The filter was immediately inserted into a 50 ml conical tube, flash frozen in liquid nitrogen, and stored at −80°C until further processing. Frozen cells were re-suspended in 8 ml of polysome resuspension buffer (20 mM Tris pH 8.0, 10 mM MgCl<sub>2</sub>, 100 mM NH<sub>4</sub>Cl, and 100 µg ml<sup>−1</sup> Chloramphenicol). Re-suspended cells were pelleted by centrifugation (3000 g, 4°C, 5 min) and the supernatant was discarded. The cell pellet was re-suspended in 500 µl of polysome lysis buffer (1X polysome resuspension buffer, 5 mM CaCl<sub>2</sub>, 0.4% TritonX-100, 0.1% NP-40, and 100 U/ml RNase-free DNase [04716728001; Roche]), and transferred to an ice-cold 1.5 ml tube containing 500 µl of 0.2–0.3 µm acid-washed glass beads (G1277; Sigma). Cells were lysed by vortexing at maximum speed on a vortexer in a 4°C room (Vortex Genie 2, 10 × 30 s with 1 min cooling on ice in between). The lysate was clarified by centrifugation (20,000 g, 4°C, 10 min) and the supernatant was transferred to a fresh 1.5 ml tube. 500 µg of total RNA (A<sub>260</sub> units) was digested (25°C, 1400 rpm, 60 min, 150 µl vol) with 2 U/µg of Micrococcal nuclease (LS004797; Worthington, Lakewood, NJ). The digestion was quenched with 1.5 µl of 0.5 M EGTA, loaded on top of a 10–50% sucrose gradient and ultra-centrifuged in a SW41 rotor (35000 rpm, 4°C, 150 min). Monosomes collected by gradient fractionation (Biocomp Instruments, Canada).</p><p>For total RNA extraction, 100 µl of polysome lysate was mixed with 400 µl of RNA extraction buffer (0.3 M sodium acetate, 10 mM EDTA, pH 4.5) and the aqueous phase was extracted twice with phenol-chloroform and once with chloroform. RNA was precipitated with an equal volume of isopropanol. The pellet was washed with 70% ice-cold ethanol and re-suspended in 100 µl of 10 mM Tris pH 7.0. 10 µg of total RNA was DNase-treated and mRNA enriched using the Microbe Express kit (Invitrogen). mRNA was fragmented by heating at 95°C with a bicarbonate buffer (<xref ref-type="bibr" rid="bib19">Ingolia et al., 2009</xref>) for 20 min.</p><p>Collected monosome fractions were purified using the same phenol-chloroform method as used for total RNA extraction above. Monosomes and fragmented mRNA were then used for small RNA sequencing library preparation. Size selection, dephosphorylation, polyadenylation, reverse-transcription, circularization and PCR amplification were performed using the same protocol as in (<xref ref-type="bibr" rid="bib19">Ingolia et al., 2009</xref>). An rRNA subtraction step was carried out between the circularization and PCR amplification steps using the same protocol as in (<xref ref-type="bibr" rid="bib36">Oh et al., 2011</xref>). Typically, several samples were multiplexed for sequencing in an Illumina HiSeq sequencer such that at least 1 million reads were obtained for each sample.</p></sec><sec id="s4-9"><title>Deep-sequencing data analysis</title><p>Deep-sequencing data analysis was carried out in Bash and Python programming languages, and performed on the Harvard research computing Odyssey cluster. Main steps are summarized below and shown schematically in <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>.<list list-type="order"><list-item><p>Each 50 nt single-end read was polyA-trimmed by identifying 10 or greater number of adjacent adenines and discarding all nucleotides starting from −1nt of the polyA run. The first 5<sup>′</sup> nt of the read was also discarded since its identity was ambiguous in several reads. PolyA-trimmed reads were first aligned against all non-coding RNAs in the <italic>B. subtilis</italic> genome using bowtie aligner (ver. 0.12.7, <xref ref-type="bibr" rid="bib28">Langmead et al., 2009</xref>).</p></list-item><list-item><p>The non-coding RNA sequences were downloaded from NCBI (NC_000964.frn).</p></list-item><list-item><p>Reads that did not align to non-coding RNAs were aligned against the whole <italic>B. subtilis</italic> genome using bowtie aligner. The <italic>B. subtilis</italic> genome was downloaded from NCBI (NC_000964.fna). Only reads that had less than three mismatches with the reference genome were considered for further analysis.</p></list-item><list-item><p>Reads that aligned to the <italic>B. subtilis</italic> genome were further trimmed by 8 nt from each end to approximate the ribosome A-site coordinate. The remaining sequence was normalized by its length and assigned to the corresponding genomic coordinate, and this value was designated as the read density at this genomic coordinate during further downstream analyses.</p></list-item><list-item><p>Average ribosome and mRNA density for a single gene was calculated by summing the read density between the start and stop codon, and then normalizing by the length of the gene.</p></list-item><list-item><p>Fold-change in ribosome density for a single gene between two samples was calculated by taking the log2 of the ratio of average ribosome density between the two samples for that gene. The median value of this log2 fold-change across all genes that received a minimum of 100 reads in at least one of the two samples was then subtracted from the fold-change value for each gene. This median-subtracted log2 fold change in reported throughout this work. We note that the average ribosome density on any gene is directly proportional to the corresponding mRNA level in the absence of specific translational regulation. Hence fold-changes in ribosome density (such as the one shown in <xref ref-type="fig" rid="fig3">Figure 3D</xref>) primarily reflect fold-changes in mRNA level.</p></list-item><list-item><p>To calculate the ribosome and mRNA density at individual codons,</p><list list-type="alpha-lower"><list-item><p>The start codon was treated as a single separate codon irrespective of its identity.</p></list-item><list-item><p>Only genes with average ribosome density of at least 1 read per codon were considered.</p></list-item><list-item><p>The ribosome density at the first nucleotide of each codon was assigned as the ribosome density at that codon. For each of the 64 codons and the start codon, read density at the first nucleotide of the codon was averaged across all occurrences of that codon in a single gene and then normalized by the average ribosome density for that gene. Hence a codon without over- or under- representation of ribosome or mRNA density will have a density value equal to 1.</p></list-item><list-item><p>Genome-wide ribosome and mRNA density was calculated as the median of the individual gene read density from Step (3) across all genes that pass the threshold of Step (2).</p></list-item><list-item><p>Ribosome profiling measurement resulted in a high ribosome density at start and stop codons. It is unclear whether this increase is a true biological signal or caused by the measurement protocol (<xref ref-type="bibr" rid="bib20">Ingolia et al., 2011</xref>). Hence these codons were excluded from the ribosome density plots shown in <xref ref-type="fig" rid="fig2 fig3">Figure 2A,B and 3C</xref>. However, the increased ribosome density at these codons resulted in a concomitant decrease in ribosome density at the remaining codons due to normalization by the average ribosome density for each gene (which included the start and stop codons).</p></list-item></list></list-item></list></p></sec><sec id="s4-10"><title>Codon usage analysis</title><p>Codon Adpatiaton Index (CAI) was calculated according to the original prescription of Sharp and Li (<xref ref-type="bibr" rid="bib42">Sharp and Li, 1987</xref>). For this calculation, 68 genes that had an annotation as ‘ribosomal’ were used as the reference set of highly expressed genes. The annotations file used for this analysis was downloaded from NCBI (NC_000964.ptt). The p value for the higher frequency of TCN codons in <italic>sinR</italic> was calculated by assuming a binomial distribution of TCN and AGC/AGT codons. The p value represents the binomial probability that there are 10 or more TCN codons in the <italic>sinR</italic> gene (12 serine codons total) given the genome-wide frequency of serine codons (0.66 for TCN codons and 0.34 for AGC/AGT codons).</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>We thank P Cluzel, AL Sonenshein, T Norman, and J Weissman for valuable comments during manuscript preparation, T Gao for performing some experiments, and B Burton for use of sucrose-gradient fractionator. Preliminary work by AS in the Cluzel lab was supported by a start-up grant from Harvard University. The computations in this paper were run on the Odyssey cluster supported by the FAS Science Division Research Computing Group at Harvard University.</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>EO: Chief Scientific Officer and a Vice President at the Howard Hughes Medical Institute, one of the three founding funders of <italic>eLife</italic>.</p></fn><fn fn-type="conflict" id="conf2"><p>RL: Senior editor, <italic>eLife</italic>.</p></fn><fn fn-type="conflict" id="conf3"><p>The other authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>ARS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con2"><p>ADL, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con3"><p>YC, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con4"><p>NB, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con5"><p>RL, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con6"><p>YC, Conception and design, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con7"><p>EO, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn></fn-group></sec><sec sec-type="supplementary-material"><title>Additional files</title><supplementary-material id="SD1-data"><object-id pub-id-type="doi">10.7554/eLife.01501.014</object-id><label>Supplementary file 1.</label><caption><title>Lists of strains, plasmids, and primers.</title><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01501.014">http://dx.doi.org/10.7554/eLife.01501.014</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife01501s001.xlsx"/></supplementary-material><sec sec-type="datasets"><title>Major dataset</title><p>The following dataset was generated:</p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro1"><name><surname>Subramaniam</surname><given-names>AR</given-names></name>, <name><surname>DeLoughery</surname><given-names>A</given-names></name>, <name><surname>Bradshaw</surname><given-names>N</given-names></name>, <name><surname>Chen</surname><given-names>Y</given-names></name>, <name><surname>O’Shea</surname><given-names>E</given-names></name>, <name><surname>Losick</surname><given-names>R</given-names></name>, <name><surname>Chai</surname><given-names>Y</given-names></name>, <year>2013</year><x>, </x><source>A Serine Sensor for Multicellularity in a Bacterium</source><x>, </x><object-id pub-id-type="art-access-id">GSE50870</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE50870">http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE50870</ext-link><x>, </x><comment>Publicly available at NCBI GEO (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/">http://www.ncbi.nlm.nih.gov/geo/</ext-link>).</comment></related-object></p></sec></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Andersson</surname><given-names>SG</given-names></name><name><surname>Kurland</surname><given-names>CG</given-names></name></person-group><year>1990</year><article-title>Codon preferences in free-living microorganisms</article-title><source>Microbiological Reviews</source><volume>54</volume><fpage>198</fpage><lpage>210</lpage></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Beauregard</surname><given-names>PB</given-names></name><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Vlamakis</surname><given-names>H</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name></person-group><year>2013</year><article-title>Bacillus subtilis biofilm induction by plant polysaccharides</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>110</volume><fpage>E1621</fpage><lpage>E1630</lpage><pub-id pub-id-type="doi">10.1073/pnas.1218984110</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bernstein</surname><given-names>JA</given-names></name><name><surname>Khodursky</surname><given-names>AB</given-names></name><name><surname>Lin</surname><given-names>P-H</given-names></name><name><surname>Lin-Chao</surname><given-names>S</given-names></name><name><surname>Cohen</surname><given-names>SN</given-names></name></person-group><year>2002</year><article-title>Global analysis of mRNA decay and abundance in <italic>Escherichia coli</italic> at single-gene resolution using two-color fluorescent DNA microarrays</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>99</volume><fpage>9697</fpage><lpage>9702</lpage><pub-id pub-id-type="doi">10.1073/pnas.112318199</pub-id></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Branda</surname><given-names>SS</given-names></name><name><surname>González-Pastor</surname><given-names>JE</given-names></name><name><surname>Ben-Yehuda</surname><given-names>S</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name></person-group><year>2001</year><article-title>Fruiting body formation by Bacillus subtilis</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>98</volume><fpage>11621</fpage><lpage>11626</lpage><pub-id pub-id-type="doi">10.1073/pnas.191384198</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Chu</surname><given-names>F</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>Bistability and biofilm formation in <italic>Bacillus subtilis</italic></article-title><source>Molecular Microbiology</source><volume>67</volume><fpage>254</fpage><lpage>263</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2007.06040.x</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>A widely conserved gene cluster required for lactate utilization in <italic>Bacillus subtilis</italic> and its involvement in biofilm formation</article-title><source>Journal of Bacteriology</source><volume>191</volume><fpage>2423</fpage><lpage>2430</lpage><pub-id pub-id-type="doi">10.1128/JB.01464-08</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Norman</surname><given-names>T</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2010</year><article-title>An epigenetic switch governing daughter cell separation in <italic>Bacillus subtilis</italic></article-title><source>Genes and Development</source><volume>24</volume><fpage>754</fpage><lpage>765</lpage><pub-id pub-id-type="doi">10.1101/gad.1915010</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Norman</surname><given-names>T</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2011</year><article-title>Evidence that metabolism and chromosome copy number control mutually exclusive cell fates in <italic>Bacillus subtilis</italic></article-title><source>The EMBO Journal</source><volume>30</volume><fpage>1402</fpage><lpage>1413</lpage><pub-id pub-id-type="doi">10.1038/emboj.2011.36</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Charneski</surname><given-names>CA</given-names></name><name><surname>Hurst</surname><given-names>LD</given-names></name></person-group><year>2013</year><article-title>Positively charged residues are the major determinants of ribosomal velocity</article-title><source>PLOS Biol</source><volume>11</volume><fpage>e1001508</fpage><pub-id pub-id-type="doi">10.1371/journal.pbio.1001508</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>Y</given-names></name><name><surname>Cao</surname><given-names>S</given-names></name><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Clardy</surname><given-names>J</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Guo</surname><given-names>JH</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2012</year><article-title>A <italic>Bacillus subtilis</italic> sensor kinase involved in triggering biofilm formation on the roots of tomato plants</article-title><source>Molecular Microbiology</source><volume>85</volume><fpage>418</fpage><lpage>430</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2012.08109.x</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chu</surname><given-names>F</given-names></name><name><surname>Kearns</surname><given-names>DB</given-names></name><name><surname>Branda</surname><given-names>SS</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2006</year><article-title>Targets of the master regulator of biofilm formation in <italic>Bacillus subtilis</italic></article-title><source>Molecular Microbiology</source><volume>59</volume><fpage>1216</fpage><lpage>1228</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2005.05019.x</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chu</surname><given-names>F</given-names></name><name><surname>Kearns</surname><given-names>DB</given-names></name><name><surname>McLoon</surname><given-names>A</given-names></name><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>A novel regulatory protein governing biofilm formation in <italic>Bacillus subtilis</italic></article-title><source>Molecular Microbiology</source><volume>68</volume><fpage>1117</fpage><lpage>1127</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2008.06201.x</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cruz Ramos</surname><given-names>H</given-names></name><name><surname>Hoffmann</surname><given-names>T</given-names></name><name><surname>Marino</surname><given-names>M</given-names></name><name><surname>Nedjari</surname><given-names>H</given-names></name><name><surname>Presecan-Siedel</surname><given-names>E</given-names></name><name><surname>Dreesen</surname><given-names>O</given-names></name><name><surname>Glaser</surname><given-names>P</given-names></name><name><surname>Jahn</surname><given-names>D</given-names></name></person-group><year>2000</year><article-title>Fermentative metabolism of <italic>Bacillus subtilis</italic>: physiology and regulation of gene expression</article-title><source>Journal of Bacteriology</source><volume>182</volume><fpage>3072</fpage><lpage>3080</lpage><pub-id pub-id-type="doi">10.1128/JB.182.11.3072-3080.2000</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dittmar</surname><given-names>KA</given-names></name><name><surname>Sorensen</surname><given-names>MA</given-names></name><name><surname>Elf</surname><given-names>J</given-names></name><name><surname>Ehrenberg</surname><given-names>M</given-names></name><name><surname>Pan</surname><given-names>T</given-names></name></person-group><year>2005</year><article-title>Selective charging of tRNA isoacceptors induced by amino-acid starvation</article-title><source>EMBO Reports</source><volume>6</volume><fpage>151</fpage><lpage>157</lpage><pub-id pub-id-type="doi">10.1038/sj.embor.7400341</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Doi</surname><given-names>RH</given-names></name><name><surname>Kaneko</surname><given-names>I</given-names></name><name><surname>Goehler</surname><given-names>B</given-names></name></person-group><year>1966</year><article-title>Regulation of a serine transfer RNA of <italic>Bacillus subtilis</italic> under two growth conditions</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>56</volume><fpage>1548</fpage><lpage>1551</lpage><pub-id pub-id-type="doi">10.1073/pnas.56.5.1548</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Elf</surname><given-names>J</given-names></name><name><surname>Nilsson</surname><given-names>D</given-names></name><name><surname>Tenson</surname><given-names>T</given-names></name><name><surname>Ehrenberg</surname><given-names>M</given-names></name></person-group><year>2003</year><article-title>Selective charging of tRNA isoacceptors explains patterns of codon usage</article-title><source>Science</source><volume>300</volume><fpage>1718</fpage><lpage>1722</lpage><pub-id pub-id-type="doi">10.1126/science.1083811</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gagnon</surname><given-names>Y</given-names></name><name><surname>Breton</surname><given-names>R</given-names></name><name><surname>Putzer</surname><given-names>H</given-names></name><name><surname>Pelchat</surname><given-names>M</given-names></name><name><surname>Grunberg-Manago</surname><given-names>M</given-names></name><name><surname>Lapointe</surname><given-names>J</given-names></name></person-group><year>1994</year><article-title>Clustering and co-transcription of the <italic>Bacillus subtilis</italic> genes encoding the aminoacyl-tRNA synthetases specific for glutamate and for cysteine and the first enzyme for cysteine biosynthesis</article-title><source>Journal of Biological Chemistry</source><volume>269</volume><fpage>7473</fpage><lpage>7482</lpage></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hayes</surname><given-names>CS</given-names></name><name><surname>Sauer</surname><given-names>RT</given-names></name></person-group><year>2003</year><article-title>Cleavage of the A site mRNA codon during ribosome pausing provides a mechanism for translational quality control</article-title><source>Molecular Cell</source><volume>12</volume><fpage>903</fpage><lpage>911</lpage><pub-id pub-id-type="doi">10.1016/S1097-2765(03)00385-X</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Henkin</surname><given-names>TM</given-names></name><name><surname>Yanofsky</surname><given-names>C</given-names></name></person-group><year>2002</year><article-title>Regulation by transcription attenuation in bacteria: how RNA provides instructions for transcription termination/antitermination decisions</article-title><source>Bioessays</source><volume>24</volume><fpage>700</fpage><lpage>707</lpage><pub-id pub-id-type="doi">10.1002/bies.10125</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ingolia</surname><given-names>NT</given-names></name><name><surname>Ghaemmaghami</surname><given-names>S</given-names></name><name><surname>Newman</surname><given-names>JRS</given-names></name><name><surname>Weissman</surname><given-names>JS</given-names></name></person-group><year>2009</year><article-title>Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling</article-title><source>Science</source><volume>324</volume><fpage>218</fpage><lpage>223</lpage><pub-id pub-id-type="doi">10.1126/science.1168978</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ingolia</surname><given-names>NT</given-names></name><name><surname>Lareau</surname><given-names>LF</given-names></name><name><surname>Weissman</surname><given-names>JS</given-names></name></person-group><year>2011</year><article-title>Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes</article-title><source>Cell</source><volume>147</volume><fpage>789</fpage><lpage>802</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2011.10.002</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ingolia</surname><given-names>NT</given-names></name><name><surname>Brar</surname><given-names>GA</given-names></name><name><surname>Rouskin</surname><given-names>S</given-names></name><name><surname>McGeachy</surname><given-names>AM</given-names></name><name><surname>Weissman</surname><given-names>JS</given-names></name></person-group><year>2012</year><article-title>The ribosome profiling strategy for monitoring translation in vivo by deep sequencing of ribosome-protected mRNA fragments</article-title><source>Nature Protocols</source><volume>7</volume><fpage>1534</fpage><lpage>1550</lpage><pub-id pub-id-type="doi">10.1038/nprot.2012.086</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jiang</surname><given-names>M</given-names></name><name><surname>Shao</surname><given-names>W</given-names></name><name><surname>Perego</surname><given-names>M</given-names></name><name><surname>Hoch</surname><given-names>JA</given-names></name></person-group><year>2000</year><article-title>Multiple histidine kinases regulate entry into stationary phase and sporulation in <italic>Bacillus subtilis</italic></article-title><source>Molecular Microbiology</source><volume>38</volume><fpage>535</fpage><lpage>542</lpage><pub-id pub-id-type="doi">10.1046/j.1365-2958.2000.02148.x</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kearns</surname><given-names>DB</given-names></name><name><surname>Chu</surname><given-names>F</given-names></name><name><surname>Branda</surname><given-names>SS</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2005</year><article-title>A master regulator for biofilm formation by <italic>Bacillus subtilis</italic></article-title><source>Molecular Microbiology</source><volume>55</volume><fpage>739</fpage><lpage>749</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2958.2004.04440.x</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Keiler</surname><given-names>KC</given-names></name><name><surname>Waller</surname><given-names>PRH</given-names></name><name><surname>Sauer</surname><given-names>RT</given-names></name></person-group><year>1996</year><article-title>Role of a peptide tagging system in degradation of proteins synthesized from damaged messenger RNA</article-title><source>Science</source><volume>271</volume><fpage>990</fpage><lpage>993</lpage><pub-id pub-id-type="doi">10.1126/science.271.5251.990</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kolodkin-Gal</surname><given-names>I</given-names></name><name><surname>Elsholz</surname><given-names>AK</given-names></name><name><surname>Muth</surname><given-names>C</given-names></name><name><surname>Girguis</surname><given-names>PR</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2013</year><article-title>Respiration control of multicellularity in <italic>Bacillus subtilis</italic> by a complex of the cytochrome chain with a membrane-embedded histidine kinase</article-title><source>Genes and Development</source><volume>27</volume><fpage>887</fpage><lpage>899</lpage><pub-id pub-id-type="doi">10.1101/gad.215244.113</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Greenberg</surname><given-names>EP</given-names></name></person-group><year>2006</year><article-title>Microbial sciences: the superficial life of microbes</article-title><source>Nature</source><volume>441</volume><fpage>300</fpage><lpage>302</lpage><pub-id pub-id-type="doi">10.1038/441300a</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kudla</surname><given-names>G</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name><name><surname>Tollervey</surname><given-names>D</given-names></name><name><surname>Plotkin</surname><given-names>JB</given-names></name></person-group><year>2009</year><article-title>Coding-sequence determinants of gene expression in <italic>Escherichia coli</italic></article-title><source>Science</source><volume>324</volume><fpage>255</fpage><lpage>258</lpage><pub-id pub-id-type="doi">10.1126/science.1170160</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Langmead</surname><given-names>B</given-names></name><name><surname>Trapnell</surname><given-names>C</given-names></name><name><surname>Pop</surname><given-names>M</given-names></name><name><surname>Salzberg</surname><given-names>SL</given-names></name></person-group><year>2009</year><article-title>Ultrafast and memory-efficient alignment of short DNA sequences to the human genome</article-title><source>Genome Biol</source><volume>10</volume><fpage>R25</fpage><pub-id pub-id-type="doi">10.1186/gb-2009-10-3-r25</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>G-W</given-names></name><name><surname>Oh</surname><given-names>E</given-names></name><name><surname>Weissman</surname><given-names>JS</given-names></name></person-group><year>2012</year><article-title>The anti-Shine-Dalgarno sequence drives translational pausing and codon choice in bacteria</article-title><source>Nature</source><volume>484</volume><fpage>538</fpage><lpage>541</lpage><pub-id pub-id-type="doi">10.1038/nature10965</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Liebs</surname><given-names>P</given-names></name><name><surname>Riedel</surname><given-names>K</given-names></name><name><surname>Graba</surname><given-names>JP</given-names></name><name><surname>Schrapel</surname><given-names>D</given-names></name><name><surname>Tischler</surname><given-names>U</given-names></name></person-group><year>1988</year><article-title>Formation of some extracellular enzymes during the exponential growth of <italic>Bacillus subtilis</italic></article-title><source>Folia Microbiologica</source><volume>33</volume><fpage>88</fpage><lpage>95</lpage><pub-id pub-id-type="doi">10.1007/BF02928073</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>López</surname><given-names>D</given-names></name><name><surname>Fischbach</surname><given-names>MA</given-names></name><name><surname>Chu</surname><given-names>F</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Structurally diverse natural products that cause potassium leakage trigger multicellularity in <italic>Bacillus subtilis</italic></article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>106</volume><fpage>280</fpage><lpage>285</lpage><pub-id pub-id-type="doi">10.1073/pnas.0810940106</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lowe</surname><given-names>TM</given-names></name><name><surname>Eddy</surname><given-names>SR</given-names></name></person-group><year>1997</year><article-title>tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence</article-title><source>Nucleic Acids Research</source><volume>25</volume><fpage>955</fpage><lpage>964</lpage><pub-id pub-id-type="doi">10.1093/nar/25.5.0955</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marin</surname><given-names>FT</given-names></name></person-group><year>1976</year><article-title>Regulation of development in <italic>Dictyostelium discoideum</italic>: I. Initiation of the growth to development transition by amino acid starvation</article-title><source>Developmental Biology</source><volume>48</volume><fpage>110</fpage><lpage>117</lpage><pub-id pub-id-type="doi">10.1016/0012-1606(76)90050-6</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Molle</surname><given-names>V</given-names></name><name><surname>Fujita</surname><given-names>M</given-names></name><name><surname>Jensen</surname><given-names>ST</given-names></name><name><surname>Eichenberger</surname><given-names>P</given-names></name><name><surname>González-Pastor</surname><given-names>JE</given-names></name><name><surname>Liu</surname><given-names>JS</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2003</year><article-title>The Spo0A regulon of <italic>Bacillus subtilis</italic></article-title><source>Molecular Microbiology</source><volume>50</volume><fpage>1683</fpage><lpage>1701</lpage><pub-id pub-id-type="doi">10.1046/j.1365-2958.2003.03818.x</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Norman</surname><given-names>T</given-names></name><name><surname>Lord</surname><given-names>N</given-names></name><name><surname>Paulsson</surname><given-names>J</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2013</year><article-title>Memory and modularity in cell-fate decision making</article-title><source>Nature</source><volume>503</volume><fpage>481</fpage><lpage>486</lpage><pub-id pub-id-type="doi">10.1038/nature12804</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Oh</surname><given-names>E</given-names></name><name><surname>Becker</surname><given-names>AH</given-names></name><name><surname>Sandikci</surname><given-names>A</given-names></name><name><surname>Huber</surname><given-names>D</given-names></name><name><surname>Chaba</surname><given-names>R</given-names></name><name><surname>Gloge</surname><given-names>F</given-names></name><name><surname>Nichols</surname><given-names>RJ</given-names></name><name><surname>Typas</surname><given-names>A</given-names></name><name><surname>Gross</surname><given-names>CA</given-names></name><name><surname>Kramer</surname><given-names>G</given-names></name><name><surname>Weissman</surname><given-names>JS</given-names></name><name><surname>Bukau</surname><given-names>B</given-names></name></person-group><year>2011</year><article-title>Selective ribosome profiling reveals the cotranslational chaperone action of trigger factor in vivo</article-title><source>Cell</source><volume>147</volume><fpage>1295</fpage><lpage>1308</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2011.10.044</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Plotkin</surname><given-names>JB</given-names></name><name><surname>Kudla</surname><given-names>G</given-names></name></person-group><year>2010</year><article-title>Synonymous but not the same: the causes and consequences of codon bias</article-title><source>Nature Reviews Genetics</source><volume>12</volume><fpage>32</fpage><lpage>42</lpage><pub-id pub-id-type="doi">10.1038/nrg2899</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Prüss</surname><given-names>BM</given-names></name><name><surname>Nelms</surname><given-names>JM</given-names></name><name><surname>Park</surname><given-names>C</given-names></name><name><surname>Wolfe</surname><given-names>AJ</given-names></name></person-group><year>1994</year><article-title>Mutations in NADH:ubiquinone oxidoreductase of <italic>Escherichia coli</italic> affect growth on mixed amino acids</article-title><source>Journal of Bacteriology</source><volume>176</volume><fpage>2143</fpage><lpage>2150</lpage></element-citation></ref><ref id="bib40"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Sambrook</surname><given-names>J</given-names></name></person-group><year>2001</year><source>Molecular cloning: a laboratory Manual</source><edition>3rd ed</edition><publisher-loc>Long Island, NY</publisher-loc><publisher-name>Cold Spring Harbor Laboratory Press</publisher-name></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sezonov</surname><given-names>G</given-names></name><name><surname>Joseleau-Petit</surname><given-names>D</given-names></name><name><surname>D’Ari</surname><given-names>R</given-names></name></person-group><year>2007</year><article-title><italic>Escherichia coli</italic> physiology in Luria-Bertani broth</article-title><source>Journal of Bacteriology</source><volume>189</volume><fpage>8746</fpage><lpage>8749</lpage><pub-id pub-id-type="doi">10.1128/JB.01368-07</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname><given-names>PM</given-names></name><name><surname>Li</surname><given-names>WH</given-names></name></person-group><year>1987</year><article-title>The Codon Adaptation Index–a measure of directional synonymous codon usage bias, and its potential applications</article-title><source>Nucleic Acids Research</source><volume>15</volume><fpage>1281</fpage><lpage>1295</lpage><pub-id pub-id-type="doi">10.1093/nar/15.3.1281</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shemesh</surname><given-names>M</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name></person-group><year>2010</year><article-title>The biocide chlorine dioxide stimulates biofilm formation in <italic>Bacillus subtilis</italic> by activation of the histidine kinase KinC</article-title><source>Journal of Bacteriology</source><volume>192</volume><fpage>6352</fpage><lpage>6356</lpage><pub-id pub-id-type="doi">10.1128/JB.01025-10</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Subramaniam</surname><given-names>AR</given-names></name><name><surname>Pan</surname><given-names>T</given-names></name><name><surname>Cluzel</surname><given-names>P</given-names></name></person-group><year>2013</year><article-title>Environmental perturbations lift the degeneracy of the genetic code to regulate protein levels in bacteria</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>110</volume><fpage>2419</fpage><lpage>2424</lpage><pub-id pub-id-type="doi">10.1073/pnas.1211077110</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sundrud</surname><given-names>MS</given-names></name><name><surname>Koralov</surname><given-names>SB</given-names></name><name><surname>Feuerer</surname><given-names>M</given-names></name><name><surname>Calado</surname><given-names>DP</given-names></name><name><surname>Kozhaya</surname><given-names>AE</given-names></name><name><surname>Rhule-Smith</surname><given-names>A</given-names></name><name><surname>Lefebvre</surname><given-names>RE</given-names></name><name><surname>Unutmaz</surname><given-names>D</given-names></name><name><surname>Mazitschek</surname><given-names>R</given-names></name><name><surname>Waldner</surname><given-names>H</given-names></name><name><surname>Whitman</surname><given-names>M</given-names></name><name><surname>Keller</surname><given-names>T</given-names></name><name><surname>Rao</surname><given-names>A</given-names></name></person-group><year>2009</year><article-title>Halofuginone inhibits TH17 cell differentiation by activating the amino acid starvation response</article-title><source>Science</source><volume>324</volume><fpage>1334</fpage><lpage>1338</lpage><pub-id pub-id-type="doi">10.1126/science.1172638</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vlamakis</surname><given-names>H</given-names></name><name><surname>Chai</surname><given-names>Y</given-names></name><name><surname>Beauregard</surname><given-names>P</given-names></name><name><surname>Losick</surname><given-names>R</given-names></name><name><surname>Kolter</surname><given-names>R</given-names></name></person-group><year>2013</year><article-title>Sticking together: building a biofilm the <italic>Bacillus subtilis</italic> way</article-title><source>Nature Reviews Microbiology</source><volume>11</volume><fpage>157</fpage><lpage>168</lpage><pub-id pub-id-type="doi">10.1038/nrmicro2960</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wach</surname><given-names>A</given-names></name></person-group><year>1996</year><article-title>PCR-synthesis of marker cassettes with long flanking homology regions for gene disruptions in <italic>S. cerevisiae</italic></article-title><source>Yeast</source><volume>12</volume><fpage>259</fpage><lpage>265</lpage><pub-id pub-id-type="doi">10.1002/(SICI)1097-0061(19960315)12:3&lt;259::AID-YEA901&gt;3.0.CO;2-C</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>J</given-names></name><name><surname>Alexander</surname><given-names>P</given-names></name><name><surname>Wu</surname><given-names>L</given-names></name><name><surname>Hammer</surname><given-names>R</given-names></name><name><surname>Cleaver</surname><given-names>O</given-names></name><name><surname>McKnight</surname><given-names>SL</given-names></name></person-group><year>2009</year><article-title>Dependence of mouse embryonic stem cells on threonine catabolism</article-title><source>Science</source><volume>325</volume><fpage>435</fpage><lpage>439</lpage><pub-id pub-id-type="doi">10.1126/science.1173288</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>West</surname><given-names>AH</given-names></name><name><surname>Stock</surname><given-names>AM</given-names></name></person-group><year>2001</year><article-title>Histidine kinases and response regulator proteins in two-component signaling systems</article-title><source>Trends in Biochemical Sciences</source><volume>26</volume><fpage>369</fpage><lpage>376</lpage><pub-id pub-id-type="doi">10.1016/S0968-0004(01)01852-7</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.01501.015</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Storz</surname><given-names>Gisela</given-names></name><role>Reviewing editor</role><aff><institution>National Institute of Child Health and Human Development</institution>, <country>United States</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://elife.elifesciences.org/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled “A Serine Sensor for Multicellularity in a Bacterium” for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by two reviewers and one member of the Board of Reviewing Editors.</p><p>The Reviewing editor and the reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>All agree that the study presents an exciting novel mechanism for the regulation of the transition to biofilm formation by <italic>Bacillus subtilis</italic>. In this mechanism, the cells sense the decrease in serine and cysteine abundance through changes in the speed of translation of genes that contain certain serine codons. It is noteworthy that these codons are not correlated with the usual codon usage, indicating a possible change in the tRNA levels or in the abundance of codons using them. The manuscript raises intriguing questions about the use of similar mechanisms by other organisms for developmental transitions or when cells are confronted with reductions in the availability of specific amino acids. However, the study would be strengthened if the authors addressed the following comments:</p><p>1) The conclusion that slow ribosome elongation at UCN codons reduces production of SinR protein to elicit biofilm formation is strongly implicated but not demonstrated. The authors should show ribosome profiles of genes with both types of codons, before and after Ser starvation or biofilm formation. The reviewers appreciate that expression of some genes might be low, but pausing in specific transcripts can be detected, and pausing at low abundance transcripts such as <italic>sinR</italic> should be detected by increasing sequencing coverage.</p><p>2) The data are largely compelling regarding the effect of different serine codons on SinR level. However the authors should discuss why certain codon substitutions behave in a non-stereotypical manner (<xref ref-type="fig" rid="fig1s2">Figure 1—figure supplement 2</xref>). While the hypothesis that ribosome pausing at these codons decreases translation elongation rate and thus SinR protein levels is the more intriguing possibility, an alternative hypothesis that pausing decreases <italic>sinR</italic> transcription via an attenuation-like mechanism has not been ruled out. The main analysis in <xref ref-type="fig" rid="fig1">Figure 1</xref> involves replacing TCA codons with synonymous mutations. These three codons are toward the 5' end of the transcript where translation may affect transcription. The TCC codons are closer to the 3' end, and replacing them with AGC yields a result opposite to the one expected. Does surrounding mRNA sequence or structure have an impact?</p><p>3) The authors note that changes in expression of genes do not correlate with the codon usage of the specific Ser codons in ribosomal proteins. To understand how this is possible, the authors should determine if the codon usage in the total mRNA sample (based on the abundance of each codon in the total mRNA) changes during biofilm formation or during Ser starvation.</p><p>4) A major conclusion is that the UCN codons act as a sensor for serine starvation that is essential for biofilm formation. How is intracellular serine concentration regulated during biofilm formation? Do cells down-regulate serine biosynthesis during the transition to stationary phase? Do the authors see any transcriptional down-regulation of serine biosynthesis genes? This should be clarified. In addition, the authors mention that serine is among the first amino acids exhausted from the laboratory medium during exponential growth in complex medium. How does that apply to bacteria in non-laboratory conditions, where the presumed mechanism evolved?</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.01501.016</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) The conclusion that slow ribosome elongation at UCN codons reduces production of SinR protein to elicit biofilm formation is strongly implicated but not demonstrated. The authors should show ribosome profiles of genes with both types of codons, before and after Ser starvation or biofilm formation. The reviewers appreciate that expression of some genes might be low, but pausing in specific transcripts can be detected, and pausing at low abundance transcripts such as</italic> sinR <italic>should be detected by increasing sequencing coverage</italic>.</p><p>We thank the reviewers for pointing that we have not explicitly demonstrated ribosome pausing on the <italic>sinR</italic> transcript. To avoid misinterpretation by the reader, we now clearly state this limitation in our revised text. Our ribosome profiling measurement did not allow us to robustly infer ribosome pausing at UCN codons on any individual transcript. At best, the selective ribosome pausing at UCN codons was evident only after calculating the median over the 10 transcripts with the highest ribosome density, while the median over the 3 transcripts with the highest ribosome density was significantly affected by noise (see Author response image 1 below). In our ribosome profiling measurement, averaging over the top 10 transcripts resulted in an average read density of ∼1500 reads per codon. By comparison, the average read density for the <italic>sinR</italic> transcript was ∼3 reads per codon. Therefore, we estimate that a ∼500X higher coverage (implying ∼2 billion total sequencing reads) is required to robustly detect ribosome pausing on the <italic>sinR</italic> transcript, which will be a prohibitively expensive experiment. Nevertheless, we are not aware of any experimental evidence or plausible mechanism by which the ribosome pausing at UCN codons that is inferred from a genome-wide median, is likely to be invalid at the level of individual mRNAs such as the <italic>sinR</italic> transcript. The magnitude of the specific ribosome pauses identified in Ref. 26 relative to the serine-starvation induced pauses observed here is uncertain, making a direct quantitative comparison difficult.<fig id="fig5" position="float"><label>Author response image 1.</label><caption><p>Sequencing read density around six serine codons, calculated as a median over the 3 genes with the highest ribosome density (left) or the 10 genes with the highest ribosome density (right). The sequencing read density for each gene was normalized by the total number of sequencing reads for that gene (same procedure as in <xref ref-type="fig" rid="fig2">Figure 2</xref>).</p></caption><graphic xlink:href="elife01501f005"/></fig></p><p><italic>2) The data are largely compelling regarding the effect of different serine codons on SinR level. However the authors should discuss why certain codon substitutions behave in a non-stereotypical manner (</italic><xref ref-type="fig" rid="fig1s2"><italic>Figure 1–figure supplement 2</italic></xref><italic>). While the hypothesis that ribosome pausing at these codons decreases translation elongation rate and thus SinR protein levels is the more intriguing possibility, an alternative hypothesis that pausing decreases</italic> sinR <italic>transcription via an attenuation-like mechanism has not been ruled out. The main analysis in</italic> <xref ref-type="fig" rid="fig1"><italic>Figure 1</italic></xref> <italic>involves replacing TCA codons with synonymous mutations. These three codons are toward the 5' end of the transcript where translation may affect transcription. The TCC codons are closer to the 3' end, and replacing them with AGC yields a result opposite to the one expected. Does surrounding mRNA sequence or structure have an impact</italic>?</p><p>If an attenuation-like mechanism underlies the observed decrease in SinR protein levels due to UCN codons, such a mechanism is expected to also decrease levels of the <italic>sinR</italic> mRNA. On the contrary, in our genome-wide mRNA measurements, we observed that <italic>sinR</italic> mRNA levels increased slightly upon biofilm formation (median-normalized fold-change: 1.33). This result is also consistent with earlier observation that expression of <italic>sinR</italic> mRNA is almost constant [Gaur et al. J Bacteriol 170, 1046 (1988)]. Further, in addition to a ribosome pause site, transcriptional attenuation typically requires a leader peptide, and an anti-terminator RNA structure, neither of which are known to be present in the <italic>sinR</italic> transcript.</p><p>Introducing synonymous substitutions can affect expression levels of the encoded protein by several mechanisms including changes in mRNA stability or secondary structure. We agree with the reviewers that the local mRNA sequence is a possible determinant of the anomalous result observed with the TCC&gt;AGC <italic>sinR</italic> variant as well as the two other <italic>sinR</italic> variants. We have now indicated this possibility in the main text. Based on our experience, such context-specific effects in a few variants are unavoidable when a large number of synonymous variants of a gene are systematically interrogated.</p><p><italic>3) The authors note that changes in expression of genes do not correlate with the codon usage of the specific Ser codons in ribosomal proteins. To understand how this is possible, the authors should determine if the codon usage in the total mRNA sample (based on the abundance of each codon in the total mRNA) changes during biofilm formation or during Ser starvation</italic>.</p><p>Usage of the six serine codons in the total mRNA pool differed less than 20% between exponential phase and biofilm entry. Further, this change was not specific to the four serine-sensitive TCN codons. More generally, protein expression level does not correlate with codon usage even under nutrient-rich conditions. We expect that protein expression during amino acid limitation is determined by the aminoacylated fraction of tRNA isoacceptors cognate to the limiting amino acid. It is worth noting that the aminoacylated fraction of tRNA isoacceptors during amino acid limitation might even be inversely correlated with codon usage. The exact value of this aminoacylated fraction has been hypothesized to depend on the balance between supply and demand for aminoacylated tRNA isoacceptors during protein synthesis.</p><p><italic>4) A major conclusion is that the UCN codons act as a sensor for serine starvation that is essential for biofilm formation. How is intracellular serine concentration regulated during biofilm formation? Do cells down-regulate serine biosynthesis during the transition to stationary phase? Do the authors see any transcriptional down-regulation of serine biosynthesis genes? This should be clarified. In addition, the authors mention that serine is among the first amino acids exhausted from the laboratory medium during exponential growth in complex medium. How does that apply to bacteria in non-laboratory conditions, where the presumed mechanism evolved</italic>?</p><p>Our RNA-Seq measurements did not reveal an obvious smoking gun for how serine levels might be regulated during biofilm formation. Specifically, we did not observe decreased transcription of serine biosynthesis genes. Serine also serves as a substrate for a wide range of cellular processes such as for purine biosynthesis, biosynthesis of several amino acids and phospholipid biosynthesis. However, none of these processes were transcriptionally up-regulated upon biofilm entry. In addition, conversion of serine to pyruvate is known to be catalyzed by the SdaA/B enzymes in <italic>E. coli</italic> upon oxygen limitation. Nevertheless, the homologous <italic>B. subtilis</italic> enzymes SdaAA/AB were not transcriptionally induced upon biofilm entry. Serine levels are likely regulated by an unidentified metabolic switch that accompanies biofilm entry. Future studies using metabolic labeling with stable isotope tracers could potentially shed light on this metabolic switch. For the sake of clarity, we have added a sentence to our main text that our inability to directly measure intracellular serine levels prevents us from ruling out unlikely possibilities, such as serine sequestration or changes in serine tRNA aminoacylation, which could theoretically occur in the absence of serine starvation.</p><p>We speculate on two plausible natural environments that could underlie the evolution of a codon-based serine-sensing mechanism in bacteria. In the case of enteric bacteria such as <italic>E. coli</italic>, alternating conditions of feast and famine could result in periodic depletion of extracellular serine, similar to that observed in laboratory complex media. In the case of free-living bacteria such as <italic>B. subtilis</italic>, intracellular serine levels could fluctuate as a result of alternating planktonic and biofilm lifestyles without changes in the extracellular serine levels.</p></body></sub-article></article>