Transcriptional Bursting Research Articles

Cytoplasmic Amplification of Transcriptional Noise Generates Substantial Cell-to-Cell Variability.

Transcription is an episodic process characterized by probabilistic bursts, but how the transcriptional noise from these bursts is modulated by cellular physiology remains unclear. Using simulations and single-molecule RNA counting, we examined how cellular processes influence cell-to-cell variability (noise). The results show that RNA noise is higher in the cytoplasm than the nucleus in ∼85% of genes across diverse promoters, genomic loci, and cell types (human and mouse). Measurements show further amplification of RNA noise in the cytoplasm, fitting a model of biphasic mRNA conversion between translation- and degradation-competent states. This multi-state translation-degradation of mRNA also causes substantial noise amplification in protein levels, ultimately accounting for ∼74% of intrinsic protein variability in cell populations. Overall, the results demonstrate how noise from transcriptional bursts is intrinsically amplified by mRNA processing, leading to a large super-Poissonian variability in protein levels.

How does supercoiling regulation on a battery of RNA polymerases impact on bacterial transcription bursting?

Transcription plays an essential role in gene expression. The transcription bursting in bacteria has been suggested to be regulated by positive supercoiling accumulation in front of a transcribing RNA polymerase (RNAP) together with gyrase binding on DNA to release the supercoiling. In this work, we study the supercoiling regulation in the case of a battery of RNAPs working together on DNA by constructing a multi-state quantitative model, which allows gradual and stepwise supercoiling accumulation and release in the RNAP transcription. We solved for transcription characteristics under the multi-state bursting model for a single RNAP transcription, and then simulated for a battery of RNAPs on DNA with T7 and Escherichia coli RNAP types of traffic, respectively, probing both the average and fluctuation impacts of the supercoiling regulation. Our studies show that due to the supercoiling accumulation and release, the number of RNAP molecules loaded onto the DNA vary significantly along time in the traffic condition. Though multiple RNAPs in transcription promote the mRNA production, they also enhance the supercoiling accumulation to suppress the production. In particular, the fluctuations of the mRNA transcripts become highly pronounced for a battery of RNAPs transcribing together under the supercoiling regulation, especially for a long process of transcription elongation. In such an elongation process, though a single RNAP can work at a high duty ratio, multiple RNAPs are hardly able to do so. Our multi-state model thus provides a systematical characterization of the quantitative features of the bacterial transcription bursting; it also supports improved physical examinations on top of this general modeling framework.

Single molecule approaches for studying gene regulation in metabolic tissues.

Gene expression in metabolic tissues can be regulated at multiple levels, ranging from the control of promoter accessibilities, transcription rates, mRNA degradation rates and mRNA localization. Modulating these processes can differentially affect important performance criteria of cells. These include precision, cellular economy, rapid response and maintenance of DNA integrity. In this review we will describe how distinct strategies of gene regulation impact the trade-offs between the cells' performance criteria. We will highlight tools based on single molecule visualization of transcripts that can be used to measure promoter states, transcription rates and mRNA degradation rates in intact tissues. These approaches revealed surprising recurrent patterns in mammalian tissues, that include transcriptional bursting, nuclear retention of mRNA, and coordination of mRNA lifetimes to facilitate rapid adaptation to changing metabolic inputs. The ability to characterize gene expression at the single molecule level can uncover the design principles of gene regulation in metabolic tissues such as the liver and the pancreas.

Transcriptional burst fraction and size dynamics during lens fiber cell differentiation and detailed insights into the denucleation process

Genes are transcribed in irregular pulses of activity termed transcriptional bursts. Cellular differentiation requires coordinated gene expression; however, it is unknown whether the burst fraction (i.e. the number of active phases of transcription) or size/intensity (the number of RNA molecules produced within a burst) changes during cell differentiation. In the ocular lens, the positions of lens fiber cells correlate precisely with their differentiation status, and the most advanced cells degrade their nuclei. Here, we examined the transcriptional parameters of the β-actin and lens differentiation-specific α-, β-, and γ-crystallin genes by RNA fluorescent in situ hybridization (FISH) in the lenses of embryonic day (E) E12.5, E14.5, and E16.5 mouse embryos and newborns. We found that cellular differentiation dramatically alters the burst fraction in synchronized waves across the lens fiber cell compartment with less dramatic changes in burst intensity. Surprisingly, we observed nascent transcription of multiple genes in nuclei just before nuclear destruction. Nuclear condensation was accompanied by transfer of nuclear proteins, including histone and nonhistone proteins, to the cytoplasm. Although lens-specific deletion of the chromatin remodeler SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A member 5 (Smarca5/Snf2h) interfered with denucleation, persisting nuclei remained transcriptionally competent and exhibited changes in both burst intensity and fraction depending on the gene examined. Our results uncover the mechanisms of nascent transcriptional control during differentiation and chromatin remodeling, confirm the burst fraction as the major factor adjusting gene expression levels, and reveal transcriptional competence of fiber cell nuclei even as they approach disintegration.

Promoter-mediated diversification of transcriptional bursting dynamics following gene duplication

During the evolution of gene families, functional diversification of proteins often follows gene duplication. However, many gene families expand while preserving protein sequence. Why do cells maintain multiple copies of the same gene? Here we have addressed this question for an actin family with 17 genes encoding an identical protein. The genes have divergent flanking regions and are scattered throughout the genome. Surprisingly, almost the entire family showed similar developmental expression profiles, with their expression also strongly coupled in single cells. Using live cell imaging, we show that differences in gene expression were apparent over shorter timescales, with family members displaying different transcriptional bursting dynamics. Strong "bursty" behaviors contrasted steady, more continuous activity, indicating different regulatory inputs to individual actin genes. To determine the sources of these different dynamic behaviors, we reciprocally exchanged the upstream regulatory regions of gene family members. This revealed that dynamic transcriptional behavior is directly instructed by upstream sequence, rather than features specific to genomic context. A residual minor contribution of genomic context modulates the gene OFF rate. Our data suggest promoter diversification following gene duplication could expand the range of stimuli that regulate the expression of essential genes. These observations contextualize the significance of transcriptional bursting.

Gene transcription in bursting: a unified mode for realizing accuracy and stochasticity.

ABSTRACTThere is accumulating evidence that, from bacteria to mammalian cells, messenger RNAs (mRNAs) are produced in intermittent bursts – a much ‘noisier’ process than traditionally thought. Based on quantitative measurements at individual promoters, diverse phenomenological models have been proposed for transcriptional bursting. Nevertheless, the underlying molecular mechanisms and significance for cellular signalling remain elusive. Here, we review recent progress, address the above issues and illuminate our viewpoints with simulation results. Despite being widely used in modelling and in interpreting experimental data, the traditional two‐state model is far from adequate to describe or infer the molecular basis and stochastic principles of transcription. In bacteria, DNA supercoiling contributes to the bursting of those genes that express at high levels and are topologically constrained in short loops; moreover, low‐affinity cis‐regulatory elements and unstable protein complexes can play a key role in transcriptional regulation. Integrating data on the architecture, kinetics, and transcriptional input–output function is a promising approach to uncovering the underlying dynamic mechanism. For eukaryotes, distinct bursting features described by the multi‐scale and continuum models coincide with those predicted by four theoretically derived principles that govern how the transcription apparatus operates dynamically. This consistency suggests a unified framework for comprehending bursting dynamics at the level of the structural and kinetic basis of transcription. Moreover, the existing models can be unified by a generic model. Remarkably, transcriptional bursting enables regulatory information to be transmitted in a digital manner, with the burst frequency representing the strength of regulatory signals. Such a mode guarantees high fidelity for precise transcriptional regulation and also provides sufficient randomness for realizing cellular heterogeneity.

Modulation of transcriptional burst frequency by histone acetylation

Many mammalian genes are transcribed during short bursts of variable frequencies and sizes that substantially contribute to cell-to-cell variability. However, which molecular mechanisms determine bursting properties remains unclear. To probe putative mechanisms, we combined temporal analysis of transcription along the circadian cycle with multiple genomic reporter integrations, using both short-lived luciferase live microscopy and single-molecule RNA-FISH. Using the Bmal1 circadian promoter as our model, we observed that rhythmic transcription resulted predominantly from variations in burst frequency, while the genomic position changed the burst size. Thus, burst frequency and size independently modulated Bmal1 transcription. We then found that promoter histone-acetylation level covaried with burst frequency, being greatest at peak expression and lowest at trough expression, while remaining unaffected by the genomic location. In addition, specific deletions of ROR-responsive elements led to constitutively elevated histone acetylation and burst frequency. We then investigated the suggested link between histone acetylation and burst frequency by dCas9p300-targeted modulation of histone acetylation, revealing that acetylation levels influence burst frequency more than burst size. The correlation between acetylation levels at the promoter and burst frequency was also observed in endogenous circadian genes and in embryonic stem cell fate genes. Thus, our data suggest that histone acetylation-mediated control of transcription burst frequency is a common mechanism to control mammalian gene expression.

Properties of gene expression and chromatin structure with mechanically regulated elongation.

In recent years, physical elements of transcription have emerged as central in our understanding of gene expression. Recent work has been done introducing a simple description of the basic physical elements of transcription where RNA elongation, RNA polymerase (RNAP) rotation and DNA supercoiling are coupled (1). Here we generalize this framework to accommodate the behavior of many RNAPs operating on multiple genes on a shared piece of DNA. The resulting framework is combined with well-established stochastic processes of transcription resulting in a model which characterizes the impact of the mechanical properties of transcription on gene expression and DNA structure. Transcriptional bursting readily emerges as a common phenomenon with origins in the geometric nature of the genetic system and results in the bounding of gene expression statistics. Properties of a multiple gene system are examined with special attention paid to the role that genome composition (gene orientation, size and intergenic distance) plays in the ability of genes to transcribe. The role of transcription in shaping DNA structure is examined and the possibility of transcription driven domain formation is discussed.

Fine-tuning of noise in gene expression with nucleosome remodeling.

Engineering stochastic fluctuations of gene expression (or “noise”) is integral to precisely bias cellular-fate decisions and statistical phenotypes in both single-cell and multi-cellular systems. Epigenetic regulation has been shown to constitute a large source of noise, and thus, engineering stochasticity is deeply intertwined with epigenetics. Here, utilizing chromatin remodeling, we report that Caffeic acid phenethyl ester (CA) and Pyrimethamine (PYR), two inhibitors of BAF250a, a subunit of the Brahma-associated factor (BAF) nucleosome remodeling complex, enable differential and tunable control of noise in transcription and translation from the human immunodeficiency virus long terminal repeat promoter in a dose and time-dependent manner. CA conserves noise levels while increasing mean abundance, resulting in direct tuning of the transcriptional burst size, while PYR strictly increases transcriptional initiation frequency while conserving a constant transcriptional burst size. Time-dependent treatment with CA reveals non-continuous tuning with noise oscillating at a constant mean abundance at early time points and the burst size increasing for treatments after 5 h. Treatments combining CA and Protein Kinase C agonists result in an even larger increase of abundance while conserving noise levels with a highly non-linear increase in variance of up to 63× untreated controls. Finally, drug combinations provide non-antagonistic combinatorial tuning of gene expression noise and map a noise phase space for future applications with viral and synthetic gene vectors. Active remodeling of nucleosomes and BAF-mediated control of gene expression noise expand a toolbox for the future design and engineering of stochasticity in living systems.

Who Needs This Junk, or Genomic Dark Matter.

Centromeres (CEN), pericentromeric regions (periCEN), and subtelomeric regions (subTel) comprise the areas of constitutive heterochromatin (HChr). Tandem repeats (TRs or satellite DNA) are the main components of HChr forming no less than 10% of the mouse and human genome. HChr is assembled within distinct structures in the interphase nuclei of many species - chromocenters. In this review, the main classes of HChr repeat sequences are considered in the order of their number increase in the sequencing reads of the mouse chromocenters (ChrmC). TRs comprise ~70% of ChrmC occupying the first place. Non-LTR (-long terminal repeat) retroposons (mainly LINE, long interspersed nuclear element) are the next (~11%), and endogenous retroviruses (ERV; LTR-containing) are in the third position (~9%). HChr is not enriched with ERV in comparison with the whole genome, but there are differences in distribution of certain elements: while MaLR-like elements (ERV3) are dominant in the whole genome, intracisternal A-particles and corresponding LTR (ERV2) are prevalent in HChr. Most of LINE in ChrmC is represented by the 2-kb fragment at the end of the 2nd open reading frame and its flanking regions. Almost all tandem repeats classified as CEN or periCEN are contained in ChrmC. Our previous classification revealed 60 new mouse TR families with 29 of them being absent in ChrmC, which indicates their location on chromosome arms. TR transcription is necessary for maintenance of heterochromatic status of the HChr genome part. A burst of TR transcription is especially important in embryogenesis and other cases of radical changes in the cell program, including carcinogenesis. The recently discovered mechanism of epigenetic regulation with noncoding sequences transcripts, long noncoding RNA, and its role in embryogenesis and pluripotency maintenance is discussed.

HTLV-1: Regulating the Balance Between Proviral Latency and Reactivation.

HTLV-1 plus-strand transcription begins with the production of doubly-spliced tax/rex transcripts, the levels of which are usually undetectable in freshly isolated peripheral blood mononuclear cells (PBMCs) from HTLV-1-infected individuals. However, the presence of a sustained chronically active cytotoxic T-cell response to HTLV-1 antigens in virtually all HTLV-1-infected individuals, regardless of their proviral load, argues against complete latency of the virus in vivo. There is an immediate burst of plus-strand transcription when blood from infected individuals is cultured ex vivo. How is the HTLV-1 plus strand silenced in PBMCs? Is it silenced in other anatomical compartments within the host? What reactivates the latent provirus in fresh PBMCs? While plus-strand transcription of the provirus appears to be intermittent, the minus-strand hbz transcripts are present in a majority of cells, albeit at low levels. What regulates the difference between the 5′- and 3′-LTR promoter activities and thereby the tax-hbz interplay? Finally, T lymphocytes are a migratory population of cells that encounter variable environments in different compartments of the body. Could these micro-environment changes influence the reactivation kinetics of the provirus? In this review we discuss the questions raised above, focusing on the early events leading to HTLV-1 reactivation from latency, and suggest future research directions.

Flying the RNA Nest: Drosophila Reveals Novel Insights into the Transcriptome Dynamics of Early Development.

Early development is punctuated by a series of pervasive and fast paced transitions. These events reshape a differentiated oocyte into a totipotent embryo and allow it to gradually mount a genetic program of its own, thereby framing a new organism. Specifically, developmental transitions that ensure the maternal to embryonic control of developmental events entail a deep remodeling of transcriptional and transcriptomic landscapes. Drosophila provides an elegant and genetically tractable system to investigate these conserved changes at a dazzling developmental pace. Here, we review recent studies applying emerging technologies such as ribosome profiling, in situ Hi-C chromatin probing and live embryo RNA imaging to investigate the transcriptional dynamics at play during Drosophila embryogenesis. In light of this new literature, we revisit the main models of zygotic genome activation (ZGA). We also review the contributions played by zygotic transcription in shaping embryogenesis and explore emerging concepts of processes such as transcriptional bursting and transcriptional memory.

Post-transcriptional bursting in genes regulated by small RNA molecules.

Gene expression programs in living cells are highly dynamic due to spatiotemporal molecular signaling and inherent biochemical stochasticity. Here we study a mechanism based on molecule-to-molecule variability at the RNA level for the generation of bursts of protein production, which can lead to heterogeneity in a cell population. We develop a mathematical framework to show numerically and analytically that genes regulated post transcriptionally by small RNA molecules can exhibit such bursts due to different states of translation activity (on or off), mostly revealed in a regime of few molecules. We exploit this framework to compare transcriptional and post-transcriptional bursting and also to illustrate how to tune the resulting protein distribution with additional post-transcriptional regulations. Moreover, because RNA-RNA interactions are predictable with an energy model, we define the kinetic constants of on-off switching as functions of the two characteristic free-energy differences of the system, activation and formation, with a nonequilibrium scheme. Overall, post-transcriptional bursting represents a distinctive principle linking gene regulation to gene expression noise, which highlights the importance of the RNA layer beyond the simple information transfer paradigm and significantly contributes to the understanding of the intracellular processes from a first-principles perspective.

The Hidden Genomic and Transcriptomic Plasticity of Giant Marker Chromosomes in Cancer.

Genome amplification in the form of rings or giant rod-shaped marker chromosomes (RGMs) is a common genetic alteration in soft tissue tumors. The mitotic stability of these structures is often rescued by perfectly functioning analphoid neocentromeres, which therefore significantly contribute to cancer progression. Here, we disentangled the genomic architecture of many neocentromeres stabilizing marker chromosomes in well-differentiated liposarcoma and lung sarcomatoid carcinoma samples. In cells carrying heavily rearranged RGMs, these structures were assembled as patchworks of multiple short amplified sequences, disclosing an extremely high level of complexity and definitely ruling out the existence of regions prone to neocentromere seeding. Moreover, by studying two well-differentiated liposarcoma samples derived from the onset and the recurrence of the same tumor, we documented an expansion of the neocentromeric domain that occurred during tumor progression, which reflects a strong selective pressure acting toward the improvement of the neocentromeric functionality in cancer. In lung sarcomatoid carcinoma cells we documented, extensive "centromere sliding" phenomena giving rise to multiple, closely mapping neocentromeric epialleles on separate coexisting markers occur, likely due to the instability of neocentromeres arising in cancer cells. Finally, by investigating the transcriptional activity of neocentromeres, we came across a burst of chimeric transcripts, both by extremely complex genomic rearrangements, and cis/trans-splicing events. Post-transcriptional editing events have been reported to expand and variegate the genetic repertoire of higher eukaryotes, so they might have a determining role in cancer. The increased incidence of fusion transcripts, might act as a driving force for the genomic amplification process, together with the increased transcription of oncogenes.

Clock-dependent chromatin topology modulates circadian transcription and behavior.

The circadian clock in animals orchestrates widespread oscillatory gene expression programs, which underlie 24-h rhythms in behavior and physiology. Several studies have shown the possible roles of transcription factors and chromatin marks in controlling cyclic gene expression. However, how daily active enhancers modulate rhythmic gene transcription in mammalian tissues is not known. Using circular chromosome conformation capture (4C) combined with sequencing (4C-seq), we discovered oscillatory promoter-enhancer interactions along the 24-h cycle in the mouse liver and kidney. Rhythms in chromatin interactions were abolished in arrhythmic Bmal1 knockout mice. Deleting a contacted intronic enhancer element in the Cryptochrome 1 (Cry1) gene was sufficient to compromise the rhythmic chromatin contacts in tissues. Moreover, the deletion reduced the daily dynamics of Cry1 transcriptional burst frequency and, remarkably, shortened the circadian period of locomotor activity rhythms. Our results establish oscillating and clock-controlled promoter-enhancer looping as a regulatory layer underlying circadian transcription and behavior.

Isoform-level gene expression patterns in single-cell RNA-sequencing data.

MotivationRNA sequencing of single cells enables characterization of transcriptional heterogeneity in seemingly homogeneous cell populations. Single-cell sequencing has been applied in a wide range of researches fields. However, few studies have focus on characterization of isoform-level expression patterns at the single-cell level. In this study, we propose and apply a novel method, ISOform-Patterns (ISOP), based on mixture modeling, to characterize the expression patterns of isoform pairs from the same gene in single-cell isoform-level expression data.ResultsWe define six principal patterns of isoform expression relationships and describe a method for differential-pattern analysis. We demonstrate ISOP through analysis of single-cell RNA-sequencing data from a breast cancer cell line, with replication in three independent datasets. We assigned the pattern types to each of 16 562 isoform-pairs from 4929 genes. Among those, 26% of the discovered patterns were significant (P<0.05), while remaining patterns are possibly effects of transcriptional bursting, drop-out and stochastic biological heterogeneity. Furthermore, 32% of genes discovered through differential-pattern analysis were not detected by differential-expression analysis. Finally, the effects of drop-out events and expression levels of isoforms on ISOP's performances were investigated through simulated datasets. To conclude, ISOP provides a novel approach for characterization of isoform-level preference, commitment and heterogeneity in single-cell RNA-sequencing data.Availability and implementationThe ISOP method has been implemented as a R package and is available at https://github.com/nghiavtr/ISOP under a GPL-3 license.Supplementary information Supplementary data are available at Bioinformatics online.

Frequency Modulation of Transcriptional Bursting Enables Sensitive and Rapid Gene Regulation.

Gene regulation is a complex non-equilibrium process. Here, we show that quantitating the temporal regulation of key gene states (transcriptionally inactive, active, and refractory) provides a parsimonious framework for analyzing gene regulation. Our theory makes two non-intuitive predictions. First, for transcription factors (TFs) that regulate transcription burst frequency, as opposed to amplitude or duration, weak TF binding is sufficient to elicit strong transcriptional responses. Second, refractoriness of a gene after a transcription burst enables rapid responses to stimuli. We validate both predictions experimentally by exploiting the natural, optogenetic-like responsiveness of the Neurospora GATA-type TF White Collar Complex (WCC) to blue light. Further, we demonstrate that differential regulation of WCC target genes is caused by different gene activation rates, not different TF occupancy, and that these rates are tuned by both the core promoter and the distance between TF-binding site and core promoter. In total, our work demonstrates the relevance of a kinetic, non-equilibrium framework for understanding transcriptional regulation.

Assessing sufficiency and necessity of enhancer activities for gene expression and the mechanisms of transcription activation.

Enhancers are important genomic regulatory elements directing cell type-specific transcription. They assume a key role during development and disease, and their identification and functional characterization have long been the focus of scientific interest. The advent of next-generation sequencing and clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9-based genome editing has revolutionized the means by which we study enhancer biology. In this review, we cover recent developments in the prediction of enhancers based on chromatin characteristics and their identification by functional reporter assays and endogenous DNA perturbations. We discuss that the two latter approaches provide different and complementary insights, especially in assessing enhancer sufficiency and necessity for transcription activation. Furthermore, we discuss recent insights into mechanistic aspects of enhancer function, including findings about cofactor requirements and the role of post-translational histone modifications such as monomethylation of histone H3 Lys4 (H3K4me1). Finally, we survey how these approaches advance our understanding of transcription regulation with respect to promoter specificity and transcriptional bursting and provide an outlook covering open questions and promising developments.

Estrogen‐dependent control and cell‐to‐cell variability of transcriptional bursting

Cellular decision‐making and environmental adaptation are dependent upon a heterogeneous response of gene expression to external cues. Heterogeneity arises in transcription from random switching between transcriptionally active and inactive states, resulting in bursts of RNA synthesis. Furthermore, the cellular state influences the competency of transcription, thereby globally affecting gene expression in a cell‐specific manner. We determined how external stimuli interplay with cellular state to modulate the kinetics of bursting. To this end, single‐cell dynamics of nascent transcripts were monitored at the endogenous estrogen‐responsive GREB1 locus. Stochastic modeling of gene expression implicated a two‐state promoter model in which the estrogen stimulus modulates the frequency of transcriptional bursting. The cellular state affects transcriptional dynamics by altering initiation and elongation kinetics and acts globally, as GREB1 alleles in the same cell correlate in their transcriptional output. Our results suggest that cellular state strongly affects the first step of the central dogma of gene expression, to promote heterogeneity in the transcriptional output of isogenic cells.

NET works after all? Engineering robustness through diversity

Classical Thermodynamics restricts engineering. These restrictions are independent of mechanism and kinetics, and thereby inescapable. Forgetting these restrictions can lead to over-optimistic designs for making bio-plastics from waste, and to erroneous ideas on early or new Life on this or other planets. This can be rectified by putting the thermodynamics in place. Is every biochemical network design feasible, provided one puts classical thermodynamics in place? Or, are there other, ill-recognized, generic restrictions to bioengineering?For a while a Non-Equilibrium Thermodynamics (NET) has been trying to discover behaviors of dynamical systems away from equilibrium that are completely independent of kinetics and mechanism and thereby not engineerable. The principles discovered were of limited use to bioengineering however.We here show that processes away from equilibrium must indeed depend on kinetics and mechanism, but, importantly, not on all kinetic and mechanistic details: There are limitations to what the engineering of mechanisms and kinetics can achieve. It is of course better to recognize what is impossible before trying to engineer the impossible. Importantly, the new NET methodology also shows that system properties that are possible, can be engineered only in certain ways.The new NET methodology also enables to understand, and perhaps engineer towards, a performance that, by adjusting the network, remains optimal when conditions are changing. Using our in silico discovery tool, we show that this may indeed occur in the Archeon S. solfataricus. What we call ‘variomatic’ gear shifting is a way that some cells may use to self-engineer their ways to maximal growth rates in environments that lack robust resources, such as in environments with fluctuating oxygen levels.Population heterogeneity is another mechanism that can increase the robustness of a cell factory. We discuss a NET principle that suggests ways in which one can engineer the cells’ diversity. Transcription burst size rather than kinetics should be modulated for making a diverse population perform much better than its average.