Long-Distance Cooperative and Antagonistic RNA Polymerase Dynamics via DNA Supercoiling.

Genes are often transcribed by multiple RNA polymerases (RNAPs) at densities that can vary widely across genes and environmental conditions. Here, we provide in vitro and in vivo evidence for a built-in mechanism by which co-transcribing RNAPs display either collaborative or antagonistic dynamics over long distances (>2 kb) through transcription-induced DNA supercoiling. In Escherichia coli, when the promoter is active, co-transcribing RNAPs translocate faster than a single RNAP, but their average speed is not altered by large variations in promoter strength and thus RNAP density. Environmentally induced promoter repression reduces the elongation efficiency of already-loaded RNAPs, causing premature termination and quick synthesis arrest of no-longer-needed proteins. This negative effect appears independent of RNAP convoy formation and is abrogated by topoisomerase I activity. Antagonistic dynamics can also occur between RNAPs from divergently transcribed gene pairs. Our findings may be broadly applicable given that transcription on topologically constrained DNA is the norm across organisms.

Genetic switching by the Lac repressor is based on two-state Monod-Wyman-Changeux allostery.

High-resolution NMR spectroscopy enabled us to characterize allosteric transitions between various functional states of the dimeric Escherichia coli Lac repressor. In the absence of ligands, the dimer exists in a dynamic equilibrium between DNA-bound and inducer-bound conformations. Binding of either effector shifts this equilibrium toward either bound state. Analysis of the ternary complex between repressor, operator DNA, and inducer shows how adding the inducer results in allosteric changes that disrupt the interdomain contacts between the inducer binding and DNA binding domains and how this in turn leads to destabilization of the hinge helices and release of the Lac repressor from the operator. Based on our data, the allosteric mechanism of the induction process is in full agreement with the well-known Monod-Wyman-Changeux model.

Molecular mechanisms controlling fructose-specific memory and catabolite repression in lactose metabolism by Streptococcus mutans.

Lactose is an abundant dietary carbohydrate metabolized by the dental pathogen Streptococcus mutans. Lactose metabolism presents both classic diauxic behaviors and long-term memory, where the bacteria can pause for >11 h before initiating growth on lactose. Here, we explored mechanisms contributing to unusual aspects of regulation of the lac operon. The fructose-phosphate metabolites, F-1-P and F-6-P, could modulate the DNA-binding activities of the lactose repressor. Recombinant LacR proteins bound upstream of lacA and Gal-6-P induced the formation of different LacR-DNA complexes. Deletion of lacR resulted in strain-specific growth phenotypes on lactose, but also on a number of mono- and di-saccharides that involve the glucose-PTS or glucokinase in their catabolism. The phenotypes were consistent with the novel findings that loss of LacR altered glucose-PTS activity and expression of the gene for glucokinase. CcpA was also shown to affect lactose metabolism in vivo and to bind to the lacA promoter region in vitro. Collectively, our study reveals complex molecular circuits controlling lactose metabolism in S. mutans, where LacR and CcpA integrate cellular and environmental cues to regulate metabolism of a variety of carbohydrates that are critical to persistence and pathogenicity of S. mutans.

Subpopulation behaviors in lactose metabolism by Streptococcus mutans.

When Streptococcus mutans is transferred from a preferred carbohydrate (glucose or fructose) to lactose, initiation of growth can take several hours, and substantial amounts of glucose are released during growth. Here, S. mutans strains UA159 and GS-5 were examined for stochastic behaviors in transcription of the lac operon. Using a gfp reporter fusion, we demonstrated that induction of the lac operon occurs in only a fraction of the population, with prior exposure to carbohydrate source and strain influencing the magniture of the sub-population response. Lower glucokinase activity in GS-5 was associated with release of substantially more glucose than UA159 and significantly lower lac expression. Mutants unable to use lactose grew on lactose as the sole carbohydrate when strains with an intact lac operon were also present in the cultures, indicative of the potential for population cheating. Utilizing a set of engineered obligate cheating and non-cheating strains, we confirmed that cheating can sustain a heterogeneous population. Futher, obligate cheaters of GS-5 competed well with the non-cheaters and showed a high degree of competitive fitness in a human-derived consortium biofilm model. The results show that bet-hedging behaviors in carbohydrate metabolism may substantially influence the composition and pathogenic potential of oral biofilms.

A lacZ reporter with high activity in the human fungal pathogen Candida albicans.

Reporter genes are useful tools to study gene transcription in various organisms. For example, the lacZ gene encoding ß-galactosidase has been extensively used as a reporter in bacteria, budding yeast, fruit fly, mouse etc. However, use of this gene in the human fungal pathogen Candida albicans has been limited, probably due to low ß-galactosidase activity. Here, we describe a reporter derived from the Vibrio cholerae lacZ gene in which codons have been optimized for expression in C. albicans. The constitutively active ACT1 promoter was fused to this synthetic lacZ reporter and integrated in the C. albicans genome. High ß-galactosidase activity in liquid assays was observed for this reporter as well as coloration on X-gal plates. When the lacZ reporter expression was driven by the MET3 promoter, ß-galactosidase activity in liquid assays and coloration on X-gal plates was higher in the absence of methionine, thus recapitulating the regulation of the native MET3 gene. This synthetic lacZ gene extends the toolbox of C. albicans reagents by providing a useful reporter for analysis of promoter activity in this organism of medical importance.

CRISPR-Cas experiments for schools and the public.

The "gene scissors" CRISPR-Cas currently revolutionize the field of molecular biology with an enormous impact on society due to the broad application potentials in biomedicine, biotechnology and agriculture. We have developed simple CRISPR-Cas experiments that can serve to introduce pupils, students and non-scientists alike to the fascinating power of targeted gene editing. The experimental course is divided into two parts. In part 1, we target plasmid borne lacZ to convert blue E. coli to white E. coli. In part 2, we analyse the CRISPR-Cas9 mediated double strand breaks in the lacZ gene by a) colony PCR, b) colony cracking gel or c) restriction digest of the plasmids. Experimental work is embedded in short theoretical lecture parts that provide background of CRISPR-Cas and a step-by-step tutorial for the practical work. Though the experiment is robust, inexpensive and simple it should be noted that guidance by an expert instructor is required. Based on our experience, a full day lab course has a positive influence on the participants' attitude towards research in general. This is true for high school students as well as non-scientists (age groups 16-70 years).

Characterization of Gene Repression by Designed Transcription Activator-like Effector Dimer Proteins.

Gene regulation by control of transcription initiation is a fundamental property of living cells. Much of our understanding of gene repression originated from studies of the Escherichia coli lac operon switch, in which DNA looping plays an essential role. To validate and generalize principles from lac for practical applications, we previously described artificial DNA looping driven by designed transcription activator-like effector dimer (TALED) proteins. Because TALE monomers bind the idealized symmetrical lac operator sequence in two orientations, our prior studies detected repression due to multiple DNA loops. We now quantitatively characterize gene repression in living E. coli by a collection of individual TALED loops with systematic loop length variation. Fitting of a thermodynamic model allows unequivocal demonstration of looping and comparison of the engineered TALED repression system with the natural lac repressor system.

Genome re-sequencing and reannotation of the Escherichia coli ER2566 strain and transcriptome sequencing under overexpression conditions.

BACKGROUND: The Escherichia coli ER2566 strain (NC_CP014268.2) was developed as a BL21 (DE3) derivative strain and had been widely used in recombinant protein expression. However, like many other current RefSeq annotations, the annotation of the ER2566 strain was incomplete, with missing gene names and miscellaneous RNAs, as well as uncorrected annotations of some pseudogenes. Here, we performed a systematic reannotation of the ER2566 genome by combining multiple annotation tools with manual revision to provide a comprehensive understanding of the E. coli ER2566 strain, and used high-throughput sequencing to explore how the strain adapted under external pressure. RESULTS: The reannotation included noteworthy corrections to all protein-coding genes, led to the exclusion of 190 hypothetical genes or pseudogenes, and resulted in the addition of 237 coding sequences and 230 miscellaneous noncoding RNAs and 2 tRNAs. In addition, we further manually examined all 194 pseudogenes in the Ref-seq annotation and directly identified 123 (63%) as coding genes. We then used whole-genome sequencing and high-throughput RNA sequencing to assess mutational adaptations under consecutive subculture or overexpression burden. Whereas no mutations were detected in response to consecutive subculture, overexpression of the human papillomavirus 16 type capsid led to the identification of a mutation (position 1,094,824 within the 3' non-coding region) positioned 19-bp away from the lacI gene in the transcribed RNA, which was not detected at the genomic level by Sanger sequencing. CONCLUSION: The ER2566 strain was used by both the general scientific community and the biotechnology industry. Reannotation of the E. coli ER2566 strain not only improved the RefSeq data but uncovered a key site that might be involved in the transcription and translation of genes encoding the lactose operon repressor. We proposed that our pipeline might offer a universal method for the reannotation of other bacterial genomes with high speed and accuracy. This study might facilitate a better understanding of gene function for the ER2566 strain under external burden and provided more clues to engineer bacteria for biotechnological applications.

H2Av facilitates H3S10 phosphorylation but is not required for heat shock-induced chromatin decondensation or transcriptional elongation.

A model has been proposed in which JIL-1 kinase-mediated H3S10 and H2Av phosphorylation is required for transcriptional elongation and heat shock-induced chromatin decondensation. However, here we show that although H3S10 phosphorylation is indeed compromised in the H2Av null mutant, chromatin decondensation at heat shock loci is unaffected in the absence of JIL-1 as well as of H2Av and that there is no discernable decrease in the elongating form of RNA polymerase II in either mutant. Furthermore, mRNA for the major heat shock protein Hsp70 is transcribed at robust levels in both H2Av and JIL-1 null mutants. Using a different chromatin remodeling paradigm that is JIL-1 dependent, we provide evidence that ectopic tethering of JIL-1 and subsequent H3S10 phosphorylation recruits PARP-1 to the remodeling site independently of H2Av phosphorylation. These data strongly suggest that H2Av or H3S10 phosphorylation by JIL-1 is not required for chromatin decondensation or transcriptional elongation in Drosophila.

Escherichia coli σ70 promoters allow expression rate control at the cellular level in genome-integrated expression systems.

BACKGROUND: The genome-integrated T7 expression system offers significant advantages, in terms of productivity and product quality, even when expressing the gene of interest (GOI) from a single copy. Compared to plasmid-based expression systems, this system does not incur a plasmid-mediated metabolic load, and it does not vary the dosage of the GOI during the production process. However, long-term production with T7 expression system leads to a rapidly growing non-producing population, because the T7 RNA polymerase (RNAP) is prone to mutations. The present study aimed to investigate whether two σ70 promoters, which were recognized by the Escherichia coli host RNAP, might be suitable in genome-integrated expression systems. We applied a promoter engineering strategy that allowed control of expressing the model protein, GFP, by introducing lac operators (lacO) into the constitutive T5 and A1 promoter sequences. RESULTS: We showed that, in genome-integrated E. coli expression systems that used σ70 promoters, the number of lacO sites must be well balanced. Promoters containing three and two lacO sites exhibited low basal expression, but resulted in a complete stop in recombinant protein production in partially induced cultures. In contrast, expression systems regulated by a single lacO site and the lac repressor element, lacIQ, on the same chromosome caused very low basal expression, were highly efficient in recombinant protein production, and enables fine-tuning of gene expression levels on a cellular level. CONCLUSIONS: Based on our results, we hypothesized that this phenomenon was associated with the autoregulation of the lac repressor protein, LacI. We reasoned that the affinity of LacI for the lacO sites of the GOI must be lower than the affinity of LacI to the lacO sites of the endogenous lac operon; otherwise, LacI autoregulation could not take place, and the lack of LacI autoregulation would lead to a disturbance in lac repressor-mediated regulation of transcription. By exploiting the mechanism of LacI autoregulation, we created a novel E. coli expression system for use in recombinant protein production, synthetic biology, and metabolic engineering applications.

Stochasticity of metabolism and growth at the single-cell level.

Elucidating the role of molecular stochasticity in cellular growth is central to understanding phenotypic heterogeneity and the stability of cellular proliferation. The inherent stochasticity of metabolic reaction events should have negligible effect, because of averaging over the many reaction events contributing to growth. Indeed, metabolism and growth are often considered to be constant for fixed conditions. Stochastic fluctuations in the expression level of metabolic enzymes could produce variations in the reactions they catalyse. However, whether such molecular fluctuations can affect growth is unclear, given the various stabilizing regulatory mechanisms, the slow adjustment of key cellular components such as ribosomes, and the secretion and buffering of excess metabolites. Here we use time-lapse microscopy to measure fluctuations in the instantaneous growth rate of single cells of Escherichia coli, and quantify time-resolved cross-correlations with the expression of lac genes and enzymes in central metabolism. We show that expression fluctuations of catabolically active enzymes can propagate and cause growth fluctuations, with transmission depending on the limitation of the enzyme to growth. Conversely, growth fluctuations propagate back to perturb expression. Accordingly, enzymes were found to transmit noise to other unrelated genes via growth. Homeostasis is promoted by a noise-cancelling mechanism that exploits fluctuations in the dilution of proteins by cell-volume expansion. The results indicate that molecular noise is propagated not only by regulatory proteins but also by metabolic reactions. They also suggest that cellular metabolism is inherently stochastic, and a generic source of phenotypic heterogeneity.

Analytical Expressions and Physics for Single-Cell mRNA Distributions of the lac Operon of E. coli.

Mechanistic models of stochastic gene expression are of considerable interest, but their complexity often precludes tractable analytical expressions for messenger RNA (mRNA) and protein distributions. The lac operon of Escherichia coli is a model system with regulatory elements such as multiple operators and DNA looping that are shared by many operons. Although this system is complex, intuition suggests that fast DNA looping may simplify it by causing the repressor-bound states of the operon to equilibrate rapidly, thus ensuring that the subsequent dynamics are governed by slow transitions between the repressor-free and the equilibrated repressor-bound states. Here, we show that this intuition is correct by applying singular perturbation theory to a mechanistic model of lac transcription with the scaled time constant of DNA looping as the perturbation parameter. We find that at steady state, the repressor-bound states satisfy detailed balance and are dominated by the looped states; moreover, the interaction between the repressor-free and the equilibrated repressor-bound states is described by an extension of the Peccoud-Ycart two-state model in which both (repressor-free and repressor-bound) states support transcription. The solution of this extended two-state model reveals that the steady-state mRNA distribution is a mixture of the Poisson and negative hypergeometric distributions, which reflects mRNAs obtained by transcription from the repressor-bound and repressor-free states. Finally, we show that the physics revealed by perturbation theory makes it easy to derive the extended two-state model equations for complex regulatory architectures.

Gene expression in E. coli influences the position and motion of the lac operon and vicinal loci.

Transcription and translation of active genes play an important role in determining the global organization of the chromosome. To further elucidate this phenomenon, we examined how the expression of either the lacY or the cfp gene in the native lac operon influences adjacent chromosomal segments by fluorescently labeling loci upstream and downstream of the expressed gene. Based on the positions and motile behaviors of these loci, our results reveal that the local organization of the vicinal chromosomal segments and its position in the nucleoid are both influenced by gene expression. Furthermore, we found that the effects on local organization depend on whether the expressed gene encodes a membrane protein or a cytoplasmic protein. Our measurements showing the movement of loci toward the membrane and the correlation between the motions of the upstream and downstream loci support the conclusion that the expression of genes encoding membrane proteins greatly influences chromosome dynamics.

Progressive developmental restriction, acquisition of left-right identity and cell growth behavior during lobe formation in mouse liver development.

To identify cell-based decisions implicated in morphogenesis of the mammalian liver, we performed clonal analysis of hepatocytes/hepatoblasts in mouse liver development, using a knock-in allele of Hnf4a/laacZ This transgene randomly undergoes a low frequency of recombination that generates a functional lacZ gene that produces ß-galactosidase in tissues in which Hnf4a is expressed. Two types of ß-galactosidase-positive clones were found. Most have undergone three to eight cell divisions and result from independent events (Luria-Delbrück fluctuation test); we calculate that they arose between E8.5 and E13.5. A second class was mega-clones derived from early endoderm progenitors, generating many descendants. Some originated from multi-potential founder cells, with labeled cells in the liver, pancreas and/or intestine. A few mega-clones populate only one side of the liver, indicating hepatic cell chirality. The patterns of labeled cells indicate cohesive and often oriented growth, notably in broad radial stripes, potentially implicated in the formation of liver lobes. This retrospective clonal analysis gives novel insights into clonal origins, cell behavior of progenitors and distinct properties of endoderm cells that underlie the formation and morphogenesis of the liver.

An Oxygen-Sensing Two-Component System in the Burkholderia cepacia Complex Regulates Biofilm, Intracellular Invasion, and Pathogenicity.

Burkholderia dolosa is a member of the Burkholderia cepacia complex (BCC), which is a group of bacteria that cause chronic lung infection in patients with cystic fibrosis (CF) and can be associated with outbreaks carrying high morbidity and mortality. While investigating the genomic diversity of B. dolosa strains collected from an outbreak among CF patients, we previously identified fixL as a gene showing signs of strong positive selection. This gene has homology to fixL of the rhizobial FixL/FixJ two-component system. The goals of this study were to determine the functions of FixLJ and their role in virulence in B. dolosa. We generated a fixLJ deletion mutant and complemented controls in B. dolosa strain AU0158. Using a fixK-lacZ reporter we found that FixLJ was activated in low oxygen in multiple BCC species. In a murine pneumonia model, the B. dolosa fixLJ deletion mutant was cleared faster from the lungs and spleen than wild-type B. dolosa strain AU0158 at 7 days post infection. Interestingly, the fixLJ deletion mutant made more biofilm, albeit with altered structure, but was less motile than strain AU0158. Using RNA-seq with in vitro grown bacteria, we found ~11% of the genome was differentially expressed in the fixLJ deletion mutant relative to strain AU0158. Multiple flagella-associated genes were down-regulated in the fixLJ deletion mutant, so we also evaluated virulence of a fliC deletion mutant, which lacks a flagellum. We saw no difference in the ability of the fliC deletion mutant to persist in the murine model relative to strain AU0158, suggesting factors other than flagella caused the phenotype of decreased persistence. We found the fixLJ deletion mutant to be less invasive in human lung epithelial and macrophage-like cells. In conclusion, B. dolosa fixLJ is a global regulator that controls biofilm formation, motility, intracellular invasion/persistence, and virulence.

Exploitation of the Escherichia coli lac operon promoter for controlled recombinant protein production.

The Escherichia coli lac operon promoter is widely used as a tool to control recombinant protein production in bacteria. Here, we give a brief review of how it functions, how it is regulated, and how, based on this knowledge, a suite of lac promoter derivatives has been developed to give a controlled expression that is suitable for diverse biotechnology applications.

The Escherichia coli multiple antibiotic resistance activator protein represses transcription of the lac operon.

In Escherichia coli, the marRAB operon is a determinant for antibiotic resistance. Such phenotypes require the encoded transcription factor MarA that activates efflux pump expression. To better understand all genes controlled by MarA, we recently mapped binding of the regulator across the E. coli genome. As expected, many MarA targets were adjacent to genes encoding stress response systems. Surprisingly, one MarA-binding site overlapped the lac operon regulatory region. Here, we show that MarA specifically targets this locus and can block transcription of the lac genes. Repression is mediated by binding of MarA to a site overlapping the lacP1 promoter -35 element. Control of the lac operon by MarA does not impact antibiotic resistance.

Long-term evolution of the natural isolate of Escherichia coli 536 in the mouse gut colonized after maternal transmission reveals convergence in the constitutive expression of the lactose operon.

In vitro experimental evolution has taught us many lessons on the molecular bases of adaptation. To move towards more natural settings, evolution in the mice gut has been successfully performed. Yet, these experiments suffered from the use of laboratory strains as well as the use of axenic or streptomycin-treated mice to maintain the inoculated strains. To circumvent these limitations, we conducted a one-year experimental evolution in vivo using a natural isolate of E. coli, strain 536, in conditions mimicking as much as possible natural environment with mother-to-offspring microbiota transmission. Mice were then distributed in 24 independent cages and separated into two different diets: a regular one (chow diet, CD) and high-fat and high-sugar one (Western Diet, WD). Genome sequences revealed an early and rapid selection during the breastfeeding period that selected the constitutive expression of the well-characterized lactose operon. E. coli was lost significantly more in CD than WD; however, we could not detect any genomic signature of selection, nor any diet specificities during the later part of the experiments. The apparently neutral evolution presumably due to low population size maintained nevertheless at high frequency the early selected mutations affecting lactose regulation. The rapid loss of lactose operon regulation challenges the idea that plastic gene expression is both optimal and stable in the wild.

Price of disorder in the lac repressor hinge helix.

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.

Quantifying heterogeneity of stochastic gene expression.

The heterogeneity of stochastic gene expression, which refers to the temporal fluctuation in a gene product and its cell-to-cell variation, has attracted considerable interest from biologists, physicists, and mathematicians. The dynamics of protein production and degradation have been modeled as random processes with transition probabilities. However, there is a gap between theory and phenomena, particularly in terms of analytical formulation and parameter estimation. In this study, we propose a theoretical framework in which we present a basic model of a gene regulatory system, derive a steady-state solution, and provide a Bayesian approach for estimating the model parameters from single-cell experimental data. The proposed framework is demonstrated to be applicable for various scales of single-cell experiments at both the mRNA and protein levels and is useful for comparing kinetic parameters across species, genomes, and cell strains.