1
|
Stirling F, Bitzan L, O'Keefe S, Redfield E, Oliver JWK, Way J, Silver PA. Rational Design of Evolutionarily Stable Microbial Kill Switches. Mol Cell 2017; 68:686-697.e3. [PMID: 29149596 DOI: 10.1016/j.molcel.2017.10.033] [Citation(s) in RCA: 74] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Revised: 08/11/2017] [Accepted: 10/24/2017] [Indexed: 12/12/2022]
Abstract
The evolutionary stability of synthetic genetic circuits is key to both the understanding and application of genetic control elements. One useful but challenging situation is a switch between life and death depending on environment. Here are presented "essentializer" and "cryodeath" circuits, which act as kill switches in Escherichia coli. The essentializer element induces cell death upon the loss of a bi-stable cI/Cro memory switch. Cryodeath makes use of a cold-inducible promoter to express a toxin. We employ rational design and a toxin/antitoxin titering approach to produce and screen a small library of potential constructs, in order to select for constructs that are evolutionarily stable. Both kill switches were shown to maintain functionality in vitro for at least 140 generations. Additionally, cryodeath was shown to control the growth environment of a population, with an escape frequency of less than 1 in 105 after 10 days of growth in the mammalian gut.
Collapse
Affiliation(s)
- Finn Stirling
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA; Wyss Institute for Biologically Inspired Engineering, Harvard University, 3 Blackfan Circle, 5th Floor, Boston, MA 02115, USA
| | - Lisa Bitzan
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA
| | - Samuel O'Keefe
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA
| | - Elizabeth Redfield
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA
| | - John W K Oliver
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA; Wyss Institute for Biologically Inspired Engineering, Harvard University, 3 Blackfan Circle, 5th Floor, Boston, MA 02115, USA
| | - Jeffrey Way
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA; Wyss Institute for Biologically Inspired Engineering, Harvard University, 3 Blackfan Circle, 5th Floor, Boston, MA 02115, USA
| | - Pamela A Silver
- Department of Systems Biology, Harvard Medical School, 200 Longwood Avenue, Warren Alpert 536, Boston, MA 02115, USA; Wyss Institute for Biologically Inspired Engineering, Harvard University, 3 Blackfan Circle, 5th Floor, Boston, MA 02115, USA.
| |
Collapse
|
2
|
Haustenne L, Bastin G, Hols P, Fontaine L. Modeling of the ComRS Signaling Pathway Reveals the Limiting Factors Controlling Competence in Streptococcus thermophilus. Front Microbiol 2015; 6:1413. [PMID: 26733960 PMCID: PMC4686606 DOI: 10.3389/fmicb.2015.01413] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2015] [Accepted: 11/27/2015] [Indexed: 12/25/2022] Open
Abstract
In streptococci, entry in competence is dictated by ComX abundance. In Streptococcus thermophilus, production of ComX is transient and tightly regulated during growth: it is positively regulated by the cell-cell communication system ComRS during the activation phase and negatively regulated during the shut-off phase by unidentified late competence gene(s). Interestingly, most S. thermophilus strains are not or weakly transformable in permissive growth conditions (i.e., chemically defined medium, CDM), suggesting that some players of the ComRS regulatory pathway are limiting. Here, we combined mathematical modeling and experimental approaches to identify the components of the ComRS system which are critical for both dynamics and amplitude of ComX production in S. thermophilus. We built a deterministic, population-scaled model of the time-course regulation of specific ComX production in CDM growth conditions. Strains LMD-9 and LMG18311 were respectively selected as representative of highly and weakly transformable strains. Results from in silico simulations and in vivo luciferase activities show that ComR concentration is the main limiting factor for the level of comX expression and controls the kinetics of spontaneous competence induction in strain LMD-9. In addition, the model predicts that the poor transformability of strain LMG18311 results from a 10-fold lower comR expression level compared to strain LMD-9. In agreement, comR overexpression in both strains was shown to induce higher competence levels with deregulated kinetics patterns during growth. In conclusion, we propose that the level of ComR production is one important factor that could explain competence heterogeneity among S. thermophilus strains.
Collapse
Affiliation(s)
- Laurie Haustenne
- Biochimie, Biophysique et Génétique des Microorganismes, Institut des Sciences de la Vie, Université catholique de Louvain Louvain-la-Neuve, Belgium
| | - Georges Bastin
- Center for Systems Engineering and Applied Mechanics, ICTEAM, Université catholique de Louvain Louvain-la-Neuve, Belgium
| | - Pascal Hols
- Biochimie, Biophysique et Génétique des Microorganismes, Institut des Sciences de la Vie, Université catholique de Louvain Louvain-la-Neuve, Belgium
| | - Laetitia Fontaine
- Biochimie, Biophysique et Génétique des Microorganismes, Institut des Sciences de la Vie, Université catholique de Louvain Louvain-la-Neuve, Belgium
| |
Collapse
|
3
|
Mina P, di Bernardo M, Savery NJ, Tsaneva-Atanasova K. Modelling emergence of oscillations in communicating bacteria: a structured approach from one to many cells. J R Soc Interface 2013; 10:20120612. [PMID: 23135248 DOI: 10.1098/rsif.2012.0612] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Population-level measurements of phenotypic behaviour in biological systems may not necessarily reflect individual cell behaviour. To assess qualitative changes in the behaviour of a single cell, when alone and when part of a community, we developed an agent-based model describing the metabolic states of a population of quorum-coupled cells. The modelling is motivated by published experimental work of a synthetic genetic regulatory network (GRN) used in Escherichia coli cells that exhibit oscillatory behaviour across the population. To decipher the mechanisms underlying oscillations in the system, we investigate the behaviour of the model via numerical simulation and bifurcation analysis. In particular, we study the effect of an increase in population size as well as the spatio-temporal behaviour of the model. Our results demonstrate that oscillations are possible only in the presence of a high concentration of the coupling chemical and are due to a time scale separation in key regulatory components of the system. The model suggests that the population establishes oscillatory behaviour as the system's preferred stable state. This is achieved via an effective increase in coupling across the population. We conclude that population effects in GRN design need to be taken into consideration and be part of the design process. This is important in planning intervention strategies or designing specific cell behaviours.
Collapse
Affiliation(s)
- Petros Mina
- Bristol Center for Complexity Sciences, University of Bristol, Bristol, UK.
| | | | | | | |
Collapse
|
4
|
Ta LH, Hansen LM, Sause WE, Shiva O, Millstein A, Ottemann KM, Castillo AR, Solnick JV. Conserved transcriptional unit organization of the cag pathogenicity island among Helicobacter pylori strains. Front Cell Infect Microbiol 2012; 2:46. [PMID: 22919637 PMCID: PMC3417554 DOI: 10.3389/fcimb.2012.00046] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2011] [Accepted: 03/17/2012] [Indexed: 12/14/2022] Open
Abstract
The Helicobacter pyloricag pathogenicity island (cag PAI) encodes a type IV secretion system that is more commonly found in strains isolated from patients with gastroduodenal disease than from those with asymptomatic gastritis. Genome-wide organization of the transcriptional units in H. pylori strain 26695 was recently established using RNA sequence analysis (Sharma et al., 2010). Here we used quantitative reverse-transcription polymerase chain reaction of open reading frames and intergenic regions to identify putative cag PAI operons in H. pylori; these operons were analyzed further by transcript profiling after deletion of selected promoter regions. Additionally, we used a promoter-trap system to identify functional cag PAI promoters. The results demonstrated that expression of genes on the H. pyloricag PAI varies by nearly five orders of magnitude and that the organization of cag PAI genes into transcriptional units is conserved among several H. pylori strains, including, 26695, J99, G27, and J166. We found evidence for 20 transcripts within the cag PAI, many of which likely overlap. Our data suggests that there are at least 11 operons: cag1-4, cag3-4, cag10-9, cag8-7, cag6-5, cag11-12, cag16-17, cag19-18, cag21-20, cag23-22, and cag25-24, as well as five monocistronic genes (cag4, cag13, cag14, cag15, and cag26). Additionally, the location of four of our functionally identified promoters suggests they are directing expression of, in one case, a truncated version of cag26 and in the other three, transcripts that are antisense to cag7, cag17, and cag23. We verified expression of two of these antisense transcripts, those antisense to cag17 and cag23, by reverse-transcription polymerase chain reaction. Taken together, our results suggest that the cag PAI transcriptional profile is generally conserved among H. pylori strains, 26695, J99, G27, and J166, and is likely complex.
Collapse
Affiliation(s)
- Linda H Ta
- Departments of Medicine and Microbiology & Immunology, Center for Comparative Medicine, University of California Davis, Davis, CA, USA
| | | | | | | | | | | | | | | |
Collapse
|
5
|
Garcia HG, Lee HJ, Boedicker JQ, Phillips R. Comparison and calibration of different reporters for quantitative analysis of gene expression. Biophys J 2011; 101:535-44. [PMID: 21806921 DOI: 10.1016/j.bpj.2011.06.026] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2010] [Revised: 05/11/2011] [Accepted: 06/09/2011] [Indexed: 01/05/2023] Open
Abstract
Absolute levels of gene expression in bacteria are observed to vary over as much as six orders of magnitude. Thermodynamic models have been proposed as a tool to describe the expression levels of a given transcriptional circuit. In this context, it is essential to understand both the limitations and linear range of the different methods for measuring gene expression and to determine to what extent measurements from different reporters can be directly compared with one aim being the stringent testing of theoretical descriptions of gene expression. In this article, we compare two protein reporters by measuring both the absolute level of expression and fold-change in expression using the fluorescent protein EYFP and the enzymatic reporter β-galactosidase. We determine their dynamic and linear range and show that they are interchangeable for measuring mean levels of expression over four orders of magnitude. By calibrating these reporters such that they can be interpreted in terms of absolute molecular counts, we establish limits for their applicability: autofluorescence on the lower end of expression for EYFP (at ∼10 molecules per cell) and interference with cellular growth on the high end for β-galactosidase (at ∼20,000 molecules per cell). These qualities make the reporters complementary and necessary when trying to experimentally verify the predictions from the theoretical models.
Collapse
Affiliation(s)
- Hernan G Garcia
- Department of Physics, California Institute of Technology, Pasadena, California, USA
| | | | | | | |
Collapse
|
6
|
In vivo expression of Helicobacter pylori virulence genes in patients with gastritis, ulcer, and gastric cancer. Infect Immun 2011; 80:594-601. [PMID: 22124657 DOI: 10.1128/iai.05845-11] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Abstract
The best-studied Helicobacter pylori virulence factor associated with development of peptic ulcer disease or gastric cancer (GC) rather than asymptomatic nonatrophic gastritis (NAG) is the cag pathogenicity island (cagPAI), which encodes a type IV secretion system (T4SS) that injects the CagA oncoprotein into host epithelial cells. Here we used real-time reverse transcription-PCR (RT-PCR) to measure the in vivo expression of genes on the cagPAI and of other virulence genes in patients with NAG, duodenal ulcer (DU), or GC. In vivo expression of H. pylori virulence genes was greater overall in gastric biopsy specimens of patients with GC than in those of patients with NAG or DU. However, since in vitro expression of cagA was not greater in H. pylori strains from patients with GC than in those from patients with NAG or DU, increased expression in GC in vivo is likely a result of environmental conditions in the gastric mucosa, though it may in turn cause more severe pathology. Increased expression of virulence genes in GC may represent a stress response to elevated pH or other environmental conditions in the stomach of patients with GC, which may be less hospitable to H. pylori colonization than the acidic environment in patients with NAG or DU.
Collapse
|
7
|
Abstract
To examine the previously proposed retroregulation model of spc mRNA degradation, two strains of Escherichia coli B/r were used; one has wild-type spc and lac operons and the other has a lac operon deletion, a wild-type spc operon, and a Pspc-rplN-lacZ fusion operon lacking the normal control sites of the spc operon (rplN is the first gene in the spc operon of ribosomal proteins). The decay of rplN mRNA and of lacZ mRNA in these strains was determined during exponential growth at different rates and after transcript initiation was inhibited by the antibiotic rifampicin. Functional decay of lacZ mRNA was monitored by measurements of beta-galactosidase activity and chemical decay was monitored using probes complementary to rplN, rplX, and to the 5' and 3'-terminal sections of lacZ. Analysis of the data was based on the assumption that the decay involves an endonucleolytic cleavage that functionally inactivates the mRNA and that this is followed by exonucleolytic degradation of the cleavage products. The major conclusions were: (1) During exponential growth, lacZ mRNA of the lac operon was translated about twice as frequently as lacZ mRNA of the spc-lac fusion, and both kinds of lacZ mRNA were translated at an elevated rate in the presence of rifampicin. (2) For lacZ mRNA from the lac operon, the endonuclease inactivation reaction was not affected by rifampicin, but the exonuclease reaction was inhibited. (3) The decay of rplN mRNA from the spc operon was accelerated in the presence of rifampicin; the average life was estimated to be six minutes during exponential growth in LB medium, and 2.8 minutes in the presence of rifampicin. (4) The decay of the rplN section of mRNA from the spc-lac operon fusion was coupled to the decay of the downstream lacZ mRNA section and was strongly inhibited (i.e. partially blocked) in the presence of rifampicin. These results show that the decay of spc mRNA differs in some important aspects from the decay of lac mRNA and support the retroregulation model. Moreover, the results indicate that rifampicin can have a significant and selective impact on the kinetics of both mRNA translation and decay.
Collapse
Affiliation(s)
- S T Liang
- Program in Molecular and Cell Biology, University of Texas at Dallas, Richardson, TX, 75083-0688, USA
| | | | | | | |
Collapse
|
8
|
Abstract
The S gene of bacteriophage lambda encodes the holin required for release of the R endolysin at the onset of phage-induced host lysis. S is the promoter-proximal gene on the single lambda late transcript and spans 107 codons. S has a novel translational initiation region with dual start codons, resulting in the production of two protein products, S105 and S107. Although differing only by the Met-1-Lys-2... N-terminal extension present on S107, the two proteins are thought to have opposing functions, with the shorter polypeptide acting as the lysis effector and the longer one acting as an inhibitor. The expression of wild-type and mutant alleles of the holin gene has been assessed quantitatively with respect to the scheduling of lysis. S mRNA accumulates during the late gene expression period to a final level of about 170 molecules per cell and is maintained at that level for at least the last 15 min before lysis. Total S protein synthesis, partitioned at about 2:1 in favor of the S105 protein compared with the other product, S107, accumulates to a final level of approximately 4,600 molecules per cell. The kinetics of accumulation of S is consistent with a constant translational rate of less than one S protein per mRNA per minute. Mutant alleles with alterations in the translational initiation region were studied to determine how the translational initiation region of S achieves the proper partition of initiation events at the two S start codons and how the synthesis of S105 and S107 relates to lysis timing. The results are discussed in terms of a model for the pathway by which the 30S ribosome-fMet-tRNA complex binds to the translational initiation region of S. In addition, analysis of the relationship between lysis timing and the levels of the two S gene products suggests that S107 inhibits S105, the lethal lysis effector, by a stoichiometric titration.
Collapse
Affiliation(s)
- C Y Chang
- Department of Biology, Texas A&M University, College Station 77843-2128, USA
| | | | | |
Collapse
|
9
|
Harayama S. Induction kinetics of RNA and proteins in exponentially growing organisms. Biochem Biophys Res Commun 1991; 180:913-9. [PMID: 1719975 DOI: 10.1016/s0006-291x(05)81152-0] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
A mathematical model of the induction kinetics of RNAs and proteins in exponentially growing organisms is derived, and the cellular concentrations of the induced macromolecules at a given time after induction are related to three parameters: the fraction of the synthesis of these macromolecules in total synthesis, the half life of the inducible macromolecules, and the generation time of the organisms. The model predicts that the concentrations of the inducible macromolecules reach one half of the maximum induction level within one generation time after the onset of the induction. The model also predicts that induction curves of proteins are parabolic when their mRNAs are short-lived, but sigmoid when they are stable. Observed induction curves of beta-galactosidase in Escherichia coli cells fit in the theoretical induction curves.
Collapse
Affiliation(s)
- S Harayama
- Department of Medical Biochemistry, University Medical Center, Geneva, Switzerland
| |
Collapse
|
10
|
Abstract
A method has been developed to determine the synthesis rate and lifetime of bacterial mRNAs, either bulk mRNA or specific mRNAs, with a minimum of physiological disturbance. The method uses hybridization of pulse-labeled RNA to specific probes followed by an evaluation based on a computer simulation of the labeling kinetics of different classes of RNA. The method was applied to the determination of bulk mRNA in Escherichia coli growing in glucose minimal medium: 60% of the instantaneous rate of RNA synthesis, or 2.3% of the total amount of RNA, was found to be mRNA with an average lifetime of 1.0 +/- 0.2 min (= 0.7 min half-life).
Collapse
Affiliation(s)
- E Baracchini
- Biology Program, University of Texas, Dallas 75083-0688
| | | |
Collapse
|
11
|
Little R, Ryals J, Bremer H. Physiological characterization of Escherichia coli rpoB mutants with abnormal control of ribosome synthesis. J Bacteriol 1983; 155:1162-70. [PMID: 6193095 PMCID: PMC217812 DOI: 10.1128/jb.155.3.1162-1170.1983] [Citation(s) in RCA: 19] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
We have previously reported the isolation of Escherichia coli rpoB mutants in which the control of ribosome synthesis by the nucleotide effector guanosine tetraphosphate (ppGpp) is altered, owing to a 20-fold increased sensitivity of the mutant RNA polymerases to ppGpp. In these mutants, the level of ppGpp during exponential growth is decreased about 10-fold, relative to that of rpoB+ wild-type strains, such that a near normal partitioning of RNA polymerase occurs with respect to stable RNA (rRNA and tRNA) gene activity. Here, the physiological effects of two different rpoB alleles in a relA+ and relA background were analyzed in greater detail by comparison with their isogenic rpoB+ wild-type parents. For a given growth medium, the rpoB mutations were found to affect four parameters which resulted in a reduction of growth rate. The results reinforce a previous conclusion that a key element in control of the bacterial growth rate is a mutual relationship between control of ribosome synthesis by ppGpp and control of relA-independent ppGpp metabolism by the concentration and function of ribosomes.
Collapse
|
12
|
|