@article{11713, abstract = {Objective: MazF is a sequence-specific endoribonuclease-toxin of the MazEF toxin–antitoxin system. MazF cleaves single-stranded ribonucleic acid (RNA) regions at adenine–cytosine–adenine (ACA) sequences in the bacterium Escherichia coli. The MazEF system has been used in various biotechnology and synthetic biology applications. In this study, we infer how ectopic mazF overexpression affects production of heterologous proteins. To this end, we quantified the levels of fluorescent proteins expressed in E. coli from reporters translated from the ACA-containing or ACA-less messenger RNAs (mRNAs). Additionally, we addressed the impact of the 5′-untranslated region of these reporter mRNAs under the same conditions by comparing expression from mRNAs that comprise (canonical mRNA) or lack this region (leaderless mRNA). Results: Flow cytometry analysis indicates that during mazF overexpression, fluorescent proteins are translated from the canonical as well as leaderless mRNAs. Our analysis further indicates that longer mazF overexpression generally increases the concentration of fluorescent proteins translated from ACA-less mRNAs, however it also substantially increases bacterial population heterogeneity. Finally, our results suggest that the strength and duration of mazF overexpression should be optimized for each experimental setup, to maximize the heterologous protein production and minimize the amount of phenotypic heterogeneity in bacterial populations, which is unfavorable in biotechnological processes.}, author = {Nikolic, Nela and Sauert, Martina and Albanese, Tanino G. and Moll, Isabella}, issn = {1756-0500}, journal = {BMC Research Notes}, keywords = {General Biochemistry, Genetics and Molecular Biology, General Medicine}, publisher = {Springer Nature}, title = {{Quantifying heterologous gene expression during ectopic MazF production in Escherichia coli}}, doi = {10.1186/s13104-022-06061-9}, volume = {15}, year = {2022}, } @article{11843, abstract = {A key attribute of persistent or recurring bacterial infections is the ability of the pathogen to evade the host’s immune response. Many Enterobacteriaceae express type 1 pili, a pre-adapted virulence trait, to invade host epithelial cells and establish persistent infections. However, the molecular mechanisms and strategies by which bacteria actively circumvent the immune response of the host remain poorly understood. Here, we identified CD14, the major co-receptor for lipopolysaccharide detection, on mouse dendritic cells (DCs) as a binding partner of FimH, the protein located at the tip of the type 1 pilus of Escherichia coli. The FimH amino acids involved in CD14 binding are highly conserved across pathogenic and non-pathogenic strains. Binding of the pathogenic strain CFT073 to CD14 reduced DC migration by overactivation of integrins and blunted expression of co-stimulatory molecules by overactivating the NFAT (nuclear factor of activated T-cells) pathway, both rate-limiting factors of T cell activation. This response was binary at the single-cell level, but averaged in larger populations exposed to both piliated and non-piliated pathogens, presumably via the exchange of immunomodulatory cytokines. While defining an active molecular mechanism of immune evasion by pathogens, the interaction between FimH and CD14 represents a potential target to interfere with persistent and recurrent infections, such as urinary tract infections or Crohn’s disease.}, author = {Tomasek, Kathrin and Leithner, Alexander F and Glatzová, Ivana and Lukesch, Michael S. and Guet, Calin C and Sixt, Michael K}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Type 1 piliated uropathogenic Escherichia coli hijack the host immune response by binding to CD14}}, doi = {10.7554/eLife.78995}, volume = {11}, year = {2022}, } @article{12333, abstract = {Together, copy-number and point mutations form the basis for most evolutionary novelty, through the process of gene duplication and divergence. While a plethora of genomic data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic reporter system that can distinguish between copy-number and point mutations, we study their early and transient adaptive dynamics in real time in Escherichia coli. We find two qualitatively different routes of adaptation, depending on the level of functional improvement needed. In conditions of high gene expression demand, the two mutation types occur as a combination. However, under low gene expression demand, copy-number and point mutations are mutually exclusive; here, owing to their higher frequency, adaptation is dominated by copy-number mutations, in a process we term amplification hindrance. Ultimately, due to high reversal rates and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation, but also constrain sequence divergence over evolutionary time scales.}, author = {Tomanek, Isabella and Guet, Calin C}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Adaptation dynamics between copynumber and point mutations}}, doi = {10.7554/ELIFE.82240}, volume = {11}, year = {2022}, } @misc{12339, abstract = {Copy-number and point mutations form the basis for most evolutionary novelty through the process of gene duplication and divergence. While a plethora of genomic sequence data reveals the long-term fate of diverging coding sequences and their cis-regulatory elements, little is known about the early dynamics around the duplication event itself. In microorganisms, selection for increased gene expression often drives the expansion of gene copy-number mutations, which serves as a crude adaptation, prior to divergence through refining point mutations. Using a simple synthetic genetic system that allows us to distinguish copy-number and point mutations, we study their early and transient adaptive dynamics in real-time in Escherichia coli. We find two qualitatively different routes of adaptation depending on the level of functional improvement selected for: In conditions of high gene expression demand, the two types of mutations occur as a combination. Under low gene expression demand, negative epistasis between the two types of mutations renders them mutually exclusive. Thus, owing to their higher frequency, adaptation is dominated by copy-number mutations. Ultimately, due to high rates of reversal and pleiotropic cost, copy-number mutations may not only serve as a crude and transient adaptation but also constrain sequence divergence over evolutionary time scales.}, author = {Tomanek, Isabella and Guet, Calin C}, publisher = {Dryad}, title = {{Flow cytometry YFP and CFP data and deep sequencing data of populations evolving in galactose}}, doi = {10.5061/dryad.rfj6q57ds}, year = {2022}, } @article{10736, abstract = {Predicting function from sequence is a central problem of biology. Currently, this is possible only locally in a narrow mutational neighborhood around a wildtype sequence rather than globally from any sequence. Using random mutant libraries, we developed a biophysical model that accounts for multiple features of σ70 binding bacterial promoters to predict constitutive gene expression levels from any sequence. We experimentally and theoretically estimated that 10–20% of random sequences lead to expression and ~80% of non-expressing sequences are one mutation away from a functional promoter. The potential for generating expression from random sequences is so pervasive that selection acts against σ70-RNA polymerase binding sites even within inter-genic, promoter-containing regions. This pervasiveness of σ70-binding sites implies that emergence of promoters is not the limiting step in gene regulatory evolution. Ultimately, the inclusion of novel features of promoter function into a mechanistic model enabled not only more accurate predictions of gene expression levels, but also identified that promoters evolve more rapidly than previously thought.}, author = {Lagator, Mato and Sarikas, Srdjan and Steinrück, Magdalena and Toledo-Aparicio, David and Bollback, Jonathan P and Guet, Calin C and Tkačik, Gašper}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Predicting bacterial promoter function and evolution from random sequences}}, doi = {10.7554/eLife.64543}, volume = {11}, year = {2022}, } @article{10812, abstract = {Several promising strategies based on combining or cycling different antibiotics have been proposed to increase efficacy and counteract resistance evolution, but we still lack a deep understanding of the physiological responses and genetic mechanisms that underlie antibiotic interactions and the clinical applicability of these strategies. In antibiotic-exposed bacteria, the combined effects of physiological stress responses and emerging resistance mutations (occurring at different time scales) generate complex and often unpredictable dynamics. In this Review, we present our current understanding of bacterial cell physiology and genetics of responses to antibiotics. We emphasize recently discovered mechanisms of synergistic and antagonistic drug interactions, hysteresis in temporal interactions between antibiotics that arise from microbial physiology and interactions between antibiotics and resistance mutations that can cause collateral sensitivity or cross-resistance. We discuss possible connections between the different phenomena and indicate relevant research directions. A better and more unified understanding of drug and genetic interactions is likely to advance antibiotic therapy.}, author = {Römhild, Roderich and Bollenbach, Mark Tobias and Andersson, Dan I.}, issn = {1740-1534}, journal = {Nature Reviews Microbiology}, keywords = {General Immunology and Microbiology, Microbiology, Infectious Diseases}, pages = {478--490}, publisher = {Springer Nature}, title = {{The physiology and genetics of bacterial responses to antibiotic combinations}}, doi = {10.1038/s41579-022-00700-5}, volume = {20}, year = {2022}, } @article{10939, abstract = {Understanding and characterising biochemical processes inside single cells requires experimental platforms that allow one to perturb and observe the dynamics of such processes as well as computational methods to build and parameterise models from the collected data. Recent progress with experimental platforms and optogenetics has made it possible to expose each cell in an experiment to an individualised input and automatically record cellular responses over days with fine time resolution. However, methods to infer parameters of stochastic kinetic models from single-cell longitudinal data have generally been developed under the assumption that experimental data is sparse and that responses of cells to at most a few different input perturbations can be observed. Here, we investigate and compare different approaches for calculating parameter likelihoods of single-cell longitudinal data based on approximations of the chemical master equation (CME) with a particular focus on coupling the linear noise approximation (LNA) or moment closure methods to a Kalman filter. We show that, as long as cells are measured sufficiently frequently, coupling the LNA to a Kalman filter allows one to accurately approximate likelihoods and to infer model parameters from data even in cases where the LNA provides poor approximations of the CME. Furthermore, the computational cost of filtering-based iterative likelihood evaluation scales advantageously in the number of measurement times and different input perturbations and is thus ideally suited for data obtained from modern experimental platforms. To demonstrate the practical usefulness of these results, we perform an experiment in which single cells, equipped with an optogenetic gene expression system, are exposed to various different light-input sequences and measured at several hundred time points and use parameter inference based on iterative likelihood evaluation to parameterise a stochastic model of the system.}, author = {Davidović, Anđela and Chait, Remy P and Batt, Gregory and Ruess, Jakob}, issn = {1553-7358}, journal = {PLoS Computational Biology}, number = {3}, publisher = {Public Library of Science}, title = {{Parameter inference for stochastic biochemical models from perturbation experiments parallelised at the single cell level}}, doi = {10.1371/journal.pcbi.1009950}, volume = {18}, year = {2022}, } @article{11339, abstract = {The interaction between a cell and its environment shapes fundamental intracellular processes such as cellular metabolism. In most cases growth rate is treated as a proximal metric for understanding the cellular metabolic status. However, changes in growth rate might not reflect metabolic variations in individuals responding to environmental fluctuations. Here we use single-cell microfluidics-microscopy combined with transcriptomics, proteomics and mathematical modelling to quantify the accumulation of glucose within Escherichia coli cells. In contrast to the current consensus, we reveal that environmental conditions which are comparatively unfavourable for growth, where both nutrients and salinity are depleted, increase glucose accumulation rates in individual bacteria and population subsets. We find that these changes in metabolic function are underpinned by variations at the translational and posttranslational level but not at the transcriptional level and are not dictated by changes in cell size. The metabolic response-characteristics identified greatly advance our fundamental understanding of the interactions between bacteria and their environment and have important ramifications when investigating cellular processes where salinity plays an important role.}, author = {Glover, Georgina and Voliotis, Margaritis and Łapińska, Urszula and Invergo, Brandon M. and Soanes, Darren and O’Neill, Paul and Moore, Karen and Nikolic, Nela and Petrov, Peter and Milner, David S. and Roy, Sumita and Heesom, Kate and Richards, Thomas A. and Tsaneva-Atanasova, Krasimira and Pagliara, Stefano}, issn = {2399-3642}, journal = {Communications Biology}, publisher = {Springer Nature}, title = {{Nutrient and salt depletion synergistically boosts glucose metabolism in individual Escherichia coli cells}}, doi = {10.1038/s42003-022-03336-6}, volume = {5}, year = {2022}, } @article{15270, abstract = {Various toxic compounds disrupt bacterial physiology. While bacteria harbor defense mechanisms to mitigate the toxicity, these mechanisms are often coupled to the physiological state of the cells and become ineffective when the physiology is severely disrupted.}, author = {Le, Dai and Krasnopeeva, Ekaterina and Sinjab, Faris and Pilizota, Teuta and Kim, Minsu}, issn = {2150-7511}, journal = {mBio}, keywords = {Virology, Microbiology}, number = {4}, publisher = {American Society for Microbiology}, title = {{Active efflux leads to heterogeneous dissipation of proton motive force by protonophores in bacteria}}, doi = {10.1128/mbio.00676-21}, volume = {12}, year = {2021}, } @article{10363, abstract = {Erythropoietin enhances oxygen delivery and reduces hypoxia-induced cell death, but its pro-thrombotic activity is problematic for use of erythropoietin in treating hypoxia. We constructed a fusion protein that stimulates red blood cell production and neuroprotection without triggering platelet production, a marker for thrombosis. The protein consists of an anti-glycophorin A nanobody and an erythropoietin mutant (L108A). The mutation reduces activation of erythropoietin receptor homodimers that induce erythropoiesis and thrombosis, but maintains the tissue-protective signaling. The binding of the nanobody element to glycophorin A rescues homodimeric erythropoietin receptor activation on red blood cell precursors. In a cell proliferation assay, the fusion protein is active at 10−14 M, allowing an estimate of the number of receptor–ligand complexes needed for signaling. This fusion protein stimulates erythroid cell proliferation in vitro and in mice, and shows neuroprotective activity in vitro. Our erythropoietin fusion protein presents a novel molecule for treating hypoxia.}, author = {Lee, Jungmin and Vernet, Andyna and Gruber, Nathalie and Kready, Kasia M. and Burrill, Devin R. and Way, Jeffrey C. and Silver, Pamela A.}, issn = {1741-0134}, journal = {Protein Engineering, Design and Selection}, publisher = {Oxford University Press}, title = {{Rational engineering of an erythropoietin fusion protein to treat hypoxia}}, doi = {10.1093/protein/gzab025}, volume = {34}, year = {2021}, } @article{9046, author = {Römhild, Roderich and Andersson, Dan I.}, issn = {1553-7374}, journal = {PLoS Pathogens}, number = {1}, publisher = {Public Library of Science}, title = {{Mechanisms and therapeutic potential of collateral sensitivity to antibiotics}}, doi = {10.1371/journal.ppat.1009172}, volume = {17}, year = {2021}, } @article{9283, abstract = {Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.}, author = {Nagy-Staron, Anna A and Tomasek, Kathrin and Caruso Carter, Caroline and Sonnleitner, Elisabeth and Kavcic, Bor and Paixão, Tiago and Guet, Calin C}, issn = {2050-084X}, journal = {eLife}, keywords = {Genetics and Molecular Biology}, publisher = {eLife Sciences Publications}, title = {{Local genetic context shapes the function of a gene regulatory network}}, doi = {10.7554/elife.65993}, volume = {10}, year = {2021}, } @article{9647, abstract = {Gene expression is regulated by the set of transcription factors (TFs) that bind to the promoter. The ensuing regulating function is often represented as a combinational logic circuit, where output (gene expression) is determined by current input values (promoter bound TFs) only. However, the simultaneous arrival of TFs is a strong assumption, since transcription and translation of genes introduce intrinsic time delays and there is no global synchronisation among the arrival times of different molecular species at their targets. We present an experimentally implementable genetic circuit with two inputs and one output, which in the presence of small delays in input arrival, exhibits qualitatively distinct population-level phenotypes, over timescales that are longer than typical cell doubling times. From a dynamical systems point of view, these phenotypes represent long-lived transients: although they converge to the same value eventually, they do so after a very long time span. The key feature of this toy model genetic circuit is that, despite having only two inputs and one output, it is regulated by twenty-three distinct DNA-TF configurations, two of which are more stable than others (DNA looped states), one promoting and another blocking the expression of the output gene. Small delays in input arrival time result in a majority of cells in the population quickly reaching the stable state associated with the first input, while exiting of this stable state occurs at a slow timescale. In order to mechanistically model the behaviour of this genetic circuit, we used a rule-based modelling language, and implemented a grid-search to find parameter combinations giving rise to long-lived transients. Our analysis shows that in the absence of feedback, there exist path-dependent gene regulatory mechanisms based on the long timescale of transients. The behaviour of this toy model circuit suggests that gene regulatory networks can exploit event timing to create phenotypes, and it opens the possibility that they could use event timing to memorise events, without regulatory feedback. The model reveals the importance of (i) mechanistically modelling the transitions between the different DNA-TF states, and (ii) employing transient analysis thereof.}, author = {Petrov, Tatjana and Igler, Claudia and Sezgin, Ali and Henzinger, Thomas A and Guet, Calin C}, issn = {0304-3975}, journal = {Theoretical Computer Science}, pages = {1--16}, publisher = {Elsevier}, title = {{Long lived transients in gene regulation}}, doi = {10.1016/j.tcs.2021.05.023}, volume = {893}, year = {2021}, } @article{9746, abstract = {Evolutionary adaptation is a major source of antibiotic resistance in bacterial pathogens. Evolution-informed therapy aims to constrain resistance by accounting for bacterial evolvability. Sequential treatments with antibiotics that target different bacterial processes were previously shown to limit adaptation through genetic resistance trade-offs and negative hysteresis. Treatment with homogeneous sets of antibiotics is generally viewed to be disadvantageous, as it should rapidly lead to cross-resistance. We here challenged this assumption by determining the evolutionary response of Pseudomonas aeruginosa to experimental sequential treatments involving both heterogenous and homogeneous antibiotic sets. To our surprise, we found that fast switching between only β-lactam antibiotics resulted in increased extinction of bacterial populations. We demonstrate that extinction is favored by low rates of spontaneous resistance emergence and low levels of spontaneous cross-resistance among the antibiotics in sequence. The uncovered principles may help to guide the optimized use of available antibiotics in highly potent, evolution-informed treatment designs.}, author = {Batra, Aditi and Römhild, Roderich and Rousseau, Emilie and Franzenburg, Sören and Niemann, Stefan and Schulenburg, Hinrich}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{High potency of sequential therapy with only beta-lactam antibiotics}}, doi = {10.7554/elife.68876}, volume = {10}, year = {2021}, } @article{9822, abstract = {Attachment of adhesive molecules on cell culture surfaces to restrict cell adhesion to defined areas and shapes has been vital for the progress of in vitro research. In currently existing patterning methods, a combination of pattern properties such as stability, precision, specificity, high-throughput outcome, and spatiotemporal control is highly desirable but challenging to achieve. Here, we introduce a versatile and high-throughput covalent photoimmobilization technique, comprising a light-dose-dependent patterning step and a subsequent functionalization of the pattern via click chemistry. This two-step process is feasible on arbitrary surfaces and allows for generation of sustainable patterns and gradients. The method is validated in different biological systems by patterning adhesive ligands on cell-repellent surfaces, thereby constraining the growth and migration of cells to the designated areas. We then implement a sequential photopatterning approach by adding a second switchable patterning step, allowing for spatiotemporal control over two distinct surface patterns. As a proof of concept, we reconstruct the dynamics of the tip/stalk cell switch during angiogenesis. Our results show that the spatiotemporal control provided by our “sequential photopatterning” system is essential for mimicking dynamic biological processes and that our innovative approach has great potential for further applications in cell science.}, author = {Zisis, Themistoklis and Schwarz, Jan and Balles, Miriam and Kretschmer, Maibritt and Nemethova, Maria and Chait, Remy P and Hauschild, Robert and Lange, Janina and Guet, Calin C and Sixt, Michael K and Zahler, Stefan}, issn = {1944-8252}, journal = {ACS Applied Materials and Interfaces}, number = {30}, pages = {35545–35560}, publisher = {American Chemical Society}, title = {{Sequential and switchable patterning for studying cellular processes under spatiotemporal control}}, doi = {10.1021/acsami.1c09850}, volume = {13}, year = {2021}, } @phdthesis{10307, abstract = {Bacteria-host interactions represent a continuous trade-off between benefit and risk. Thus, the host immune response is faced with a non-trivial problem – accommodate beneficial commensals and remove harmful pathogens. This is especially difficult as molecular patterns, such as lipopolysaccharide or specific surface organelles such as pili, are conserved in both, commensal and pathogenic bacteria. Type 1 pili, tightly regulated by phase variation, are considered an important virulence factor of pathogenic bacteria as they facilitate invasion into host cells. While invasion represents a de facto passive mechanism for pathogens to escape the host immune response, we demonstrate a fundamental role of type 1 pili as active modulators of the innate and adaptive immune response.}, author = {Tomasek, Kathrin}, issn = {2663-337X}, pages = {73}, publisher = {Institute of Science and Technology Austria}, title = {{Pathogenic Escherichia coli hijack the host immune response}}, doi = {10.15479/at:ista:10307}, year = {2021}, } @unpublished{10316, abstract = {A key attribute of persistent or recurring bacterial infections is the ability of the pathogen to evade the host’s immune response. Many Enterobacteriaceae express type 1 pili, a pre-adapted virulence trait, to invade host epithelial cells and establish persistent infections. However, the molecular mechanisms and strategies by which bacteria actively circumvent the immune response of the host remain poorly understood. Here, we identified CD14, the major co-receptor for lipopolysaccharide detection, on dendritic cells as a previously undescribed binding partner of FimH, the protein located at the tip of the type 1 pilus of Escherichia coli. The FimH amino acids involved in CD14 binding are highly conserved across pathogenic and non-pathogenic strains. Binding of pathogenic bacteria to CD14 lead to reduced dendritic cell migration and blunted expression of co-stimulatory molecules, both rate-limiting factors of T cell activation. While defining an active molecular mechanism of immune evasion by pathogens, the interaction between FimH and CD14 represents a potential target to interfere with persistent and recurrent infections, such as urinary tract infections or Crohn’s disease.}, author = {Tomasek, Kathrin and Leithner, Alexander F and Glatzová, Ivana and Lukesch, Michael S. and Guet, Calin C and Sixt, Michael K}, booktitle = {bioRxiv}, publisher = {Cold Spring Harbor Laboratory}, title = {{Type 1 piliated uropathogenic Escherichia coli hijack the host immune response by binding to CD14}}, doi = {10.1101/2021.10.18.464770}, year = {2021}, } @misc{8951, abstract = {Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions, such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks remains a major challenge. Here, we use a well-defined synthetic gene regulatory network to study how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one gene regulatory network with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Our results demonstrate that changes in local genetic context can place a single transcriptional unit within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual transcriptional units, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of gene regulatory networks.}, author = {Nagy-Staron, Anna A}, keywords = {Gene regulatory networks, Gene expression, Escherichia coli, Synthetic Biology}, publisher = {Institute of Science and Technology Austria}, title = {{Sequences of gene regulatory network permutations for the article "Local genetic context shapes the function of a gene regulatory network"}}, doi = {10.15479/AT:ISTA:8951}, year = {2020}, } @misc{7383, abstract = {Organisms cope with change by employing transcriptional regulators. However, when faced with rare environments, the evolution of transcriptional regulators and their promoters may be too slow. We ask whether the intrinsic instability of gene duplication and amplification provides a generic alternative to canonical gene regulation. By real-time monitoring of gene copy number mutations in E. coli, we show that gene duplications and amplifications enable adaptation to fluctuating environments by rapidly generating copy number, and hence expression level, polymorphism. This ‘amplification-mediated gene expression tuning’ occurs on timescales similar to canonical gene regulation and can deal with rapid environmental changes. Mathematical modeling shows that amplifications also tune gene expression in stochastic environments where transcription factor-based schemes are hard to evolve or maintain. The fleeting nature of gene amplifications gives rise to a generic population-level mechanism that relies on genetic heterogeneity to rapidly tune expression of any gene, without leaving any genomic signature.}, author = {Grah, Rok}, keywords = {Matlab scripts, analysis of microfluidics, mathematical model}, publisher = {Institute of Science and Technology Austria}, title = {{Matlab scripts for the Paper: Gene Amplification as a Form of Population-Level Gene Expression regulation}}, doi = {10.15479/AT:ISTA:7383}, year = {2020}, } @article{7569, abstract = {Genes differ in the frequency at which they are expressed and in the form of regulation used to control their activity. In particular, positive or negative regulation can lead to activation of a gene in response to an external signal. Previous works proposed that the form of regulation of a gene correlates with its frequency of usage: positive regulation when the gene is frequently expressed and negative regulation when infrequently expressed. Such network design means that, in the absence of their regulators, the genes are found in their least required activity state, hence regulatory intervention is often necessary. Due to the multitude of genes and regulators, spurious binding and unbinding events, called “crosstalk”, could occur. To determine how the form of regulation affects the global crosstalk in the network, we used a mathematical model that includes multiple regulators and multiple target genes. We found that crosstalk depends non-monotonically on the availability of regulators. Our analysis showed that excess use of regulation entailed by the formerly suggested network design caused high crosstalk levels in a large part of the parameter space. We therefore considered the opposite ‘idle’ design, where the default unregulated state of genes is their frequently required activity state. We found, that ‘idle’ design minimized the use of regulation and thus minimized crosstalk. In addition, we estimated global crosstalk of S. cerevisiae using transcription factors binding data. We demonstrated that even partial network data could suffice to estimate its global crosstalk, suggesting its applicability to additional organisms. We found that S. cerevisiae estimated crosstalk is lower than that of a random network, suggesting that natural selection reduces crosstalk. In summary, our study highlights a new type of protein production cost which is typically overlooked: that of regulatory interference caused by the presence of excess regulators in the cell. It demonstrates the importance of whole-network descriptions, which could show effects missed by single-gene models.}, author = {Grah, Rok and Friedlander, Tamar}, issn = {1553-7358}, journal = {PLOS Computational Biology}, number = {2}, publisher = {Public Library of Science}, title = {{The relation between crosstalk and gene regulation form revisited}}, doi = {10.1371/journal.pcbi.1007642}, volume = {16}, year = {2020}, }