The Bacillus subtilis sin operon: an evolvable network motif

The strategy of combining genes from a regulatory protein and its antagonist within the same operon, but controlling their activities differentially, can lead to diverse regulatory functions. This protein-antagonist motif is ubiquitous and present in evolutionarily unrelated regulatory pathways. Using the sin operon from the Bacillus subtilis sporulation pathway as a model system, we built a theoretical model, parameterized it using data from the literature, and used bifurcation analyses to determine the circuit functions it could encode. The model demonstrated that this motif can generate a bistable switch with tunable control over the switching threshold and the degree of population heterogeneity. Further, the model predicted that a small perturbation of a single critical parameter can bias this architecture into functioning like a graded response, a bistable switch, an oscillator, or a pulse generator. By mapping the parameters of the model to specific DNA regions and comparing the genomic sequences of Bacillus species, we showed that phylogenetic variation tends to occur in those regions that tune the switch threshold without disturbing the circuit function. The dynamical plasticity of the protein-antagonist operon motif suggests that it is an evolutionarily convergent design selected not only for particular immediate function but also for its evolvability.     

A simple framework to describe the regulation of gene expression in prokaryotes

Based on the bimolecular mass action law and the derived mass conservation laws, we propose a mathematical framework in order to describe the regulation of gene expression in prokaryotes. It is shown that the derived models have all the qualitative properties of the activation and inhibition regulatory mechanisms observed in experiments. The basic construction considers genes as templates for protein production, where regulation processes result from activators or repressors connecting to DNA binding sites. All the parameters in the models have a straightforward biological meaning. After describing the general properties of the basic mechanisms of positive and negative gene regulation, we apply this framework to the self-regulation of the trp operon and to the genetic switch involved in the regulation of the lac operon. One of the consequences of this approach is the existence of conserved quantities depending on the initial conditions that tune bifurcations of fixed points. This leads naturally to a simple explanation of threshold effects as observed in some experiments.     

Gene Expression Switching of Receptor Subunits in Human Brain Development

Synaptic receptors in the human brain consist of multiple protein subunits, many of which have multiple variants, coded by different genes, and are differentially expressed across brain regions and developmental stages. The brain can tune the electrophysiological properties of synapses to regulate plasticity and information processing by switching from one protein variant to another. Such condition-dependent variant switch during development has been demonstrated in several neurotransmitter systems including NMDA and GABA. Here we systematically detect pairs of receptor-subunit variants that switch during the lifetime of the human brain by analyzing postmortem expression data collected in a population of donors at various ages and brain regions measured using microarray and RNA-seq. To further detect variant pairs that co-vary across subjects, we present a method to quantify age-corrected expression correlation in face of strong temporal trends. This is achieved by computing the correlations in the residual expression beyond a cubic-spline model of the population temporal trend, and can be seen as a nonlinear version of partial correlations. Using these methods, we detect multiple new pairs of context dependent variants. For instance, we find a switch from GLRA2 to GLRA3 that differs from the known switch in the rat. We also detect an early switch from HTR1A to HTR5A whose trends are negatively correlated and find that their age-corrected expression is strongly positively correlated. Finally, we observe that GRIN2B switch to GRIN2A occurs mostly during embryonic development, presumably earlier than observed in rodents. These results provide a systematic map of developmental switching in the neurotransmitter systems of the human brain.     

Variability of the distribution of differentiation pathway choices regulated by a multipotent delayed stochastic switch

We study the plasticity of a delayed stochastic model of a genetic toggle switch as a multipotent differentiation pathway switch, at the single cell and cell population levels, by observing distributions of differentiation pathways choices of genetically homogeneous cell populations. Assuming a model of stochastic pathway determination of cell differentiation that is regulated by the proteins of the switch, we vary the proteins’ expression level and degradation rates, which cells are known to be able to regulate, to vary mean level, noise, and bias of the proteins’ expression levels. It is shown that small changes in each of these dynamical features significantly and distinctively affects the dynamics of the switch at the single cell level and thus, the cell differentiation patterns. The regulation of these features allows cells to regulate their pluripotency and cell populations’ distribution of lineage choice, suggesting that the stochastic switch has high plasticity regarding differentiation pathway choice regulation, thus providing adaptability to environmental stresses and changes.     

Coordinated induction of multi-gene pathways in Saccharomyces cerevisiae

Bacterial operons are nature’s tool for regulating and coordinating multi-gene expression in prokaryotes. They are also a gene architecture commonly used in the biosynthesis of many pharmaceutically important compounds and industrially useful chemicals. Despite being an important eukaryotic production host, Saccharomyces cerevisiae has never had such gene architecture. Here, we report the development of a system to assemble and regulate a multi-gene pathway in S. cerevisiae. Full pathways can be constructed using pre-made parts from a plasmid toolbox. Subsequently, through the use of a yeast strain containing a stably integrated gene switch, the assembled pathway can be regulated using a readily available and inexpensive compound—estradiol—with extremely high sensitivity (10 nM). To demonstrate the use of the system, we assembled the five-gene zeaxanthin biosynthetic pathway in a single step and showed the ligand-dependent coordinated expression of all five genes as well as the tightly regulated production of zeaxanthin. Compared with a previously reported constitutive zeaxanthin pathway, our inducible pathway was shown to have 50-fold higher production level.     

Plasticity Regulators Modulate Specific Root Traits in Discrete Nitrogen Environments

Plant development is remarkably plastic but how precisely can the plant customize its form to specific environments? When the plant adjusts its development to different environments, related traits can change in a coordinated fashion, such that two traits co-vary across many genotypes. Alternatively, traits can vary independently, such that a change in one trait has little predictive value for the change in a second trait. To characterize such “tunability” in developmental plasticity, we carried out a detailed phenotypic characterization of complex root traits among 96 accessions of the model Arabidopsis thaliana in two nitrogen environments. The results revealed a surprising level of independence in the control of traits to environment – a highly tunable form of plasticity. We mapped genetic architecture of plasticity using genome-wide association studies and further used gene expression analysis to narrow down gene candidates in mapped regions. Mutants in genes implicated by association and expression analysis showed precise defects in the predicted traits in the predicted environment, corroborating the independent control of plasticity traits. The overall results suggest that there is a pool of genetic variability in plants that controls traits in specific environments, with opportunity to tune crop plants to a given environment.     

The regulatory roles of the galactose permease and kinase in the induction response of the GAL network in Saccharomyces cerevisiae

The GAL genetic switch of Saccharomyces cerevisiae exhibits an ultrasensitive response to the inducer galactose as well as the "all-or-none" behavior characteristic of many eukaryotic regulatory networks. We have constructed a strain that allows intermediate levels of gene expression from a tunable GAL1 promoter at both the population and the single cell level by altering the regulation of the galactose permease Gal2p. Similar modifications to other feedback loops regulating the Gal80p repressor and the Gal3p signaling protein did not result in similarly tuned responses, indicating that the level of inducer transport is unique in its ability to control the switch response of the network. In addition, removal of the Gal1p galactokinase from the network resulted in a regimed response due to the dual role of this enzyme in galactose catabolism and transport. These two activities have competing effects on the response of the network to galactose such that the transport effects of Gal1p are dominant at low galactose concentrations, whereas its catabolic effects are dominant at high galactose concentrations. In addition, flow cytometry analysis revealed the unexpected phenomenon of multiple populations in the gal1delta strains, which were not present in the isogenic GAL1 background. This result indicates that Gal1p may play a previously undescribed role in the stability of the GAL network response.     

Characterization of the yeast copper-inducible promoter system in Arabidopsis thaliana

Inducible promoters or gene-switches are used to both spatially and temporally regulate gene expression. Such regulation can provide information concerning the function of a gene in a developmental context as well as avoid potential harmful effects due to overexpression. A gfp construct under the control of a copper-inducible promoter was introduced into Arabidopsis thaliana (L.) Heynh. and the regulatory parameters of this inducible promoter were determined. Here, we describe the time-course of up- and down-regulation of GFP expression in response to copper level, the optimal regulatory levels of copper, and the tissue specificity of expression in three transgenic lines. We conclude that the copper-inducible promoter system may be useful in regulating the time and location of gene expression in A. thaliana.     

Simultaneous monitoring of pituitary hormone secretion and gene expression within individual cells.

We have recently developed a method to simultaneously quantitate the level of gene expression and the level of secretion of a peptide from individual cells. Our approach has been to combine the reverse hemolytic plaque assay sequentially with in situ hybridization. We present data to show how we have used the pituitary lactotroph as a model to demonstrate the power of this technique. However, we are particularly excited about the potential application of this strategy to approach a broad spectrum of questions regarding the cellular and molecular mechanisms that regulate the coupling of peptide secretion and gene expression at the single cell level. The method can be used in any system in which an appropriate antibody for the reverse hemolytic plaque assay and probes complementary to the mRNA of interest are available.     

Improvement of a monopartite ecdysone receptor gene switch and demonstration of its utility in regulation of transgene expression in plants

In plants, regulation of transgene expression is typically accomplished through the use of inducible promoter systems. The ecdysone receptor (EcR) gene switch is one of the best inducible systems available to regulate transgene expression in plants. However, the monopartite EcR gene switches developed to date require micromolar concentrations of ligand for activation. We tested several EcR mutants that were generated by changing one or two amino acid residues in the highly flexible ligand-binding domain of Choristoneura fumiferana EcR (CfEcR). Based on the transient expression assays, we selected a double mutant, V395I + Y415E (VY), of CfEcR (CfEcR(VY)) for further testing in stable transformation experiments. The CfEcR(VY) mutant only slightly improved the induction characteristics of the two-hybrid gene switch, whereas the CfEcR(VY) mutant significantly improved the induction characteristics of the monopartite gene switch (VGCfE(VY)). The ligand sensitivity of the VGCfE(VY) switch was improved by 125-15 625-fold in different transgenic lines analyzed, compared to the VGCfE(Wt) switch. The utility of the VGCfE(VY) switch was tested by regulating the expression of an Arabidopsis zinc finger protein gene (AtZFP11) in both tobacco and Arabidopsis plants. These data showed that the VGCfE(VY) switch efficiently regulated the expression of AtZFP11 and that the phenotype of AtZFP11 could be induced by the application of ligand. In addition, the affected plants recovered after withdrawal of the ligand, demonstrating the utility of this gene switch in regulating the expression of critical transgenes in plants.     

A directional strategy for monitoring Cre-mediated recombination at the cellular level in the mouse

Functional redundancies, compensatory mechanisms, and lethal phenotypes often prevent the full analysis of gene functions through generation of germline null mutations in the mouse. The use of site-specific recombinases, such as Cre, which catalyzes recombination between loxP sites, has allowed the engineering of mice harboring targeted somatic mutations, which are both temporally controlled and cell-type restricted. Many Cre-expressing mouse lines exist, but only a few transgenic lines are available that harbor a reporter gene whose expression is dependent on a Cre-mediated event. Moreover, their use to monitor gene ablation at the level of individual cells is often limited, as in some tissues the reporter gene may be silenced, be affected by position-effect variegation, or reside in a chromatin configuration inaccessible for recombination. Thus, one cannot validly extrapolate from the expression of a reporter transgene to an identical ablation pattern for the conditional allele of a given gene. By combining the ability of Cre recombinase to invert or excise a DNA fragment, depending on the orientation of the flanking loxP sites, and the availability of both wild-type (WT) and mutant loxP sites, we designed a Cre-dependent genetic switch (FLEx switch) through which the expression of a given gene is turned off, while the expression of another one is concomitantly turned on. We demonstrate the efficiency and reliability of this switch to readily detect, in the mouse, at the single cell level, Cre-mediated gene ablation. We discuss how this strategy can be used to generate genetic modifications in a conditional manner.