Skip to main content

Balancing noise and plasticity in eukaryotic gene expression

Abstract

Background

Coupling the control of expression stochasticity (noise) to the ability of expression change (plasticity) can alter gene function and influence adaptation. A number of factors, such as transcription re-initiation, strong chromatin regulation or genome neighboring organization, underlie this coupling. However, these factors do not necessarily combine in equivalent ways and strengths in all genes. Can we identify then alternative architectures that modulate in distinct ways the linkage of noise and plasticity?

Results

Here we first show that strong chromatin regulation, commonly viewed as a source of coupling, can lead to plasticity without noise. The nature of this regulation is relevant too, with plastic but noiseless genes being subjected to general activators whereas plastic and noisy genes experience more specific repression. Contrarily, in genes exhibiting poor transcriptional control, it is translational efficiency what separates noise from plasticity, a pattern related to transcript length. This additionally implies that genome neighboring organization –as modifier– appears only effective in highly plastic genes. In this class, we confirm bidirectional promoters (bipromoters) as a configuration capable to reduce coupling by abating noise but also reveal an important trade-off, since bipromoters also decrease plasticity. This presents ultimately a paradox between intergenic distances and modulation, with short intergenic distances both associated and disassociated to noise at different plasticity levels.

Conclusions

Balancing the coupling among different types of expression variability appears as a potential shaping force of genome regulation and organization. This is reflected in the use of different control strategies at genes with different sets of functional constraints.

Background

Variation in gene expression is observed between closely related species, even when the specific gene coding sequence is largely conserved, e.g.,[1]. Within a species, expression can fluctuate following a perturbation (environmental or genetic) and even in the absence of perturbations variation among individuals is found – this being often interpreted as disturbing noise[2]. What molecular factors determine these fluctuations? Are these factors subjected to selection pressures? And which general trends on expression variability can one identify at the genomic level?

Partial answers to these questions were recently reached by using high-throughput experiments on the budding yeast Saccharomyces cerevisiae. Noise was measured in >2,500 proteins using GFP-tagged yeast strains[3] and this validated the contribution of mRNA dynamics to protein noise. Both protein function (e.g., housekeeping proteins exhibiting low noise, while stress-response proteins being noisy) and chromatin dynamics (transitions between active/inactive states) were also shown to correlate with noise. Moreover, expression plasticity (responsiveness of S. cerevisiae genes to change in external conditions) and divergence (among closely related species) were also quantified with the use of a compendium of genome-wide expression profiles in four yeasts[4]. Genes presenting a TATA box in their promoter showed higher interspecies variability, controlling for function, which suggested the influence of transcription re-initiation mechanisms and bursting expression[5]. Similarly, response to mutations (using mutation-accumulation experiments[6]) identified TATA boxes and trans-mutational target sizes (number of proteins influencing the expression of a focal gene) as determinants of neutral variability. Finally, the production of a complete nucleosome occupancy map covering 81% of the genome[7] helped identify how different (absolute and relative) occupancy levels further controls variability.

These initial findings are leading to new questions. For instance, are all these aspects of variation (short-term –noise/plasticity– vs. long-term –divergence) linked to a unifying promoter structure? This is clearly suggested in recent studies, with an emphasis on the role of chromatin regulation[810]. This strategy could be positive in terms of the economics of regulation[11], but negative in terms of functional conflicts, e.g., need of bipolarity in genome-wide transcription[12, 13], presence of gene classes requiring precise but plastic expression[14], etc. An additional question is to what extent a demand for variation acts as a central force for the organization of genomes and vice versa, i.e., whether structural genomic features constrain variation[1517].

Here, we first revisited the influence of chromatin regulation in the linkage of noise and plasticity. We observe that both regulatory strength and character modulates this linkage. Plastic genes exhibiting relatively strong chromatin regulation can appear independent of noise, but an extra increase in plasticity associates plasticity with noise. This association –or the lack of it–is revealed in the type of chromatin control, with a contrast between global and specific regulation. While these patterns indicate transcriptional initiation as fundamental mechanism of modulation (as previously suggested, e.g.,[5]), we alternatively find that noise uncouples from plasticity in low-plastic genes due to changes in translational efficiency. These distinct modes are confirmed by the differential influence of genomic neighborhood on coupling depending on plasticity. Interestingly, short intergenic distance and bidirectional promoter architecture can both be related to high and low noise.

Results

Chromatin regulation does not always link plasticity to noise

TATA boxes and high nucleosomal occupancy at the proximal regions of transcriptional starting sites (TSSs) have been recognized as fundamental promoter features leading to gene expression variability[4, 6, 12]. Both features were shown to couple two specific forms of variability, i.e., expression noise and plasticity. Linkage between noise and plasticity was additionally associated to a highly dynamic chromatin, as quantified by histone exchange rates[14]. However, histone exchange rates do not fully describe the many trans factors influencing nucleosome dynamics.

To better understand how such factors determine the noise-plasticity coupling, we used a score that assesses chromatin regulation effects (CRE), i.e., how much the expression of a given gene varies when deleting its trans-acting chromatin regulators[810, 18] (Methods). CRE correlated with plasticity as expected (Spearman’s correlation coefficient ρ = 0.57,p < 10−20, n = 2045). We then grouped genes in terms of proximal nucleosome occupancy and computed mean plasticity for those genes exhibiting high or low CRE within each group. Notably, chromatin regulation can induce a relatively high level of plasticity independent of nucleosomal occupancy and presence of TATA promoters (Figure1). Is it possible to identify genes with particularly strong plasticity? This distinguishes genes whose promoters present both high proximal nucleosomal occupancy –increasing sensitivity to regulation (Additional file1: Figure S1)[10]– and presence of TATA box –that results in bursting transcription and increased transcriptional efficiency[5]– which in turn involves coupling to noise (confirming earlier reports[3, 8], Additional file1: Figure S2).

Figure 1
figure 1

Proximal nucleosomal occupancy, chromatin regulation, and the noise-plasticity coupling. We divided the dataset (n = 2045) in ten equally sized bins of increasing proximal nucleosome occupancy. In each bin, we computed the median chromatin regulation effect (CRE). We plotted the difference in plasticity (A, blue curve) or noise (B, red curve) of genes above/below this median and contrasted the observed values with those expected randomly (permutation test in each bin to depict significance, shown as the mean –gray curve– and mean plus two standard deviations –dashed gray curves– obtained with 10000 randomizations). Plasticity is always enhanced by strong chromatin regulation; however, regulation enhances noise only in promoters with high proximal nucleosome occupancy and TATA box (shaded area in A,B). An identical analysis is shown in (C) and (D), but excluding genes with TATA-containing promoters. High occupancy does not lead to increased plasticity/noise in this case (shaded area in C,D).

The crucial effect of the high (proximal) nucleosomal occupancy to enhance coupling is emphasized by the pronounced nucleosome depleted region (NDR) exhibited by a subset of TATA-containing genes with low noise (NDRs are similarly observed in TATAless genes, Additional file1: Figure S3). Moreover, if the noise-plasticity coupling had its origin in the stability of the transcriptional apparatus at the promoter, this would predict the presence of coupling in TATAless genes with a SAGA-dominated initiation (that also produces transcriptional bursting[19]). This is indeed what we observed (noise-plasticity correlation in TATAless and SAGA dominated genes, ρ = 0.51,p = 1.2 × 10−5, n = 66, see also Additional file1: Figure S4).

High plasticity implies different trans- regulation strategies when coupled/uncoupled to noise

To further appreciate what determines the coupling (or uncoupling) of noise with plasticity, we inspected potential differences in the type of chromatin regulation. We computed the mean effect in expression of a compendium of mutations in regulators[18] (CRE score before represents a subset, see Methods) on plastic genes. This analysis highlighted a strong anti-correlation between the effect of perturbations in low-noise high-plasticity genes (LNHP, see Methods for definition of these classes) and high-noise high-plasticity genes (HNHP, ρ=−0.83,p<2.2×10−16n=170, Figure2A; this correlation is much stronger than the expected baseline correlation, Additional file1: Figure S5). This confirms mechanistically a complementary program of regulation between these two groups of genes (enriched by growth –ribosomal– and stress genes, respectively[12, 20]), to be added to the previously observed distinctions in promoter nucleosome occupancy (Additional file1: Figure S6), and histone modification enrichment[13].

More specifically, perturbations affecting LNHP genes commonly cause a decrease in expression (Figure2A, red dots), i.e., regulators activate expression, and dominantly correspond to general transcription factors (20 out of 41). Even if most of these perturbations involved mutations in TAF1 –which is a general transcription factor (also associated to chromatin modulation activities, e.g.,[22, 23]) needed for the expression of nearly 90% of yeast genes–, its perturbation affected LNHP genes more significantly than low plasticity (LP) ones (Figure2B; this suggests these genes as preferred targets of TAF1 broad regulatory action, see Additional file1: Supplement). In contrast, a wider range of different regulators affected HNHP genes, and the majority of these were independent trans-acting chromatin regulators (92 out of 135) whose deletion activated gene expression (Figure2A, blue dots), i.e., regulators repress expression (Additional file2: Table S1).

More specifically, perturbations affecting LNHP genes commonly cause a decrease in expression (Figure2A, red dots), i.e., regulators activate expression, and dominantly correspond to general transcription factors (20 out of 41). Even if most of these perturbations involved mutations in TAF1 –which is a general transcription factor (also associated to chromatin modulation activities, e.g.,[22, 23]) needed for the expression of nearly 90% of yeast genes–, its perturbation affected LNHP genes more significantly than low plasticity (LP) ones (Figure2B; this suggests these genes as preferred targets of TAF1 broad regulatory action, see Additional file1: Supplement). In contrast, a wider range of different regulators affected HNHP genes, and the majority of these were independent trans-acting chromatin regulators (92 out of 135) whose deletion activated gene expression (Figure2A, blue dots), i.e., regulators repress expression (Additional file2: Table S1).

We also found that some regulators whose deletion results in decreasing expression level in the majority of LNHP genes tends in comparison to increase it in the majority of HNHP (Figure3). As expected, this is not related to general transcription factors, whose deletion reduce expression level, but to many specific chromatin regulators and, notably, to histones. As much as 81% of these histone deletions (see also Additional file1: Figure S7) caused an increase in the expression of most HNHP genes, while they decreased the expression level of the majority of LNHP. A recent result can help us understand this[13]. Namely, LNHP genes are greatly enriched in activating marks (mostly acetylations), a strong change in acetylation level being observed when repressed. Therefore, acetylated histones are probably essential for the expression of these genes. On the contrary, HNHP genes do not show such changes in histone acetylation status, but they reduce their occupancy level when activated. Histone deletion is in this case more likely to impede the formation of repressive nucleosomes, resulting in a more frequently open promoter and increasing expression level (see Additional file1: Supplement for more discussion, and Additional file1: Figure S8 for a detailed profile of histone modifications).

Figure 2
figure 2

Distinct chromatin regulation strategies to achieve noisy or quiet plasticity. A) Each dot represents the mean effect in the expression of a set of genes in a subclass (HNHP, x coordinate; LNHP, y coordinate; normalized by effect in HP class) when a particular regulator is mutated[18]. A ratio >1 thus implies that the corresponding subclass is more strongly influenced by certain regulator than the full HP group. A strong negative correlation is found indicating that many regulators are highly specific to either HNHP or LNHP genes. This confirms that these groups are enriched by complementary functional classes (stress and growth related genes, respectively) which are generally regulated in opposite sense[12, 13, 20]. Dot colors denote the dominant effect of the regulator on the HP class (blue; regulator is mostly repressing expression, red; regulator is commonly activating) while sizes describe the strength of the dominant effect; e.g., LNHP genes are frequently affected by strong chromatin activators. B) We examined in detail the effects on LNHP genes (box in A). Except rsc30 (a regulator of ribosomal proteins[21]) all these mutations involved TAF1, which is part of the general transcription factor TFIID[22, 23]. This essential factor regulates 90% of the genes in the genome, not including most of HNHP (which are regulated by SAGA) but including almost all LP genes (see main text). Nevertheless, we observed that all these mutations affected significantly more strongly LNHP than LP genes [K-S tests with FDR-corrected -log(p-value)’s shown at the right].

Figure 3
figure 3

Dual action of histones. Both trans-acting chromatin regulators and histones tend to have opposite effects in HNHP and LNHP genes. For each mutation from[18], we plot the fraction of genes in the HNHP group and in the LNHP group that decreased expression. Consistently with that observed in Figure2, we find that mutating as much as 50% of chromatin regulators results in the de-activation of the majority of the LNHP genes, but de-repression of most of the HNHP genes. In addition, and perhaps more importantly, we observe that 81% of mutations in histones also exhibit this behavior. This is probably crucial, and indicates that histones by themselves are needed for repression of HNHP genes and, at the same time, activation of LNHP.

Noise in low-plasticity genes arises from enhanced translational efficiency

We noted that LP genes also present differential coupling to noise. In contrast to HP genes, this disparity does not seem to respond to transcriptional-based determinants. LP genes hardly present TATA promoters (5.1% –26/513– in LP, 22.4% –343/1529– in the rest, p = 1.6 × 10−18, χ2-test), display pronounced NDRs (mean proximal nucleosomal occupancy LP: - 1.70 –n = 513–, rest: -1.46 – n = 1532– p = 7.6 × 10−7, Kolmogorov-Smirnov KS-test, see also Methods) and are poorly regulated by chromatin (mean CRE LP: 0.58 –n = 513–, rest: 0.70 –n = 1532– p < 2.2 × 10−16, KS-test). A notable feature of these genes is their enrichment in histone H2A.Z at promoters, which has been already noted[13] and is thought to help stabilizing the NDR (with our dataset, mean LP: 0.41, n = 365, rest: 0.12, n = 1144, p = 1.2 × 10−9, KS-test, see also Additional file1: Figure S8). Indeed, there are not observable differences in all these factors when considering low and high noise subgroups within the LP set (data not shown).

We thus inspected if uncoupling could be associated in these genes to translation as this is known to control noise[24, 25]. Our analysis shows that noise in LP genes is correlated with translational efficiency[26] and ribosomal density[27] (ρ = 0.22,p = 7.9 × 10−5, and ρ = 0.21,p = 1.6 × 10−4, respectively; n = 327, and Figure4A,B) while we did not observe this in highly plastic genes (ρ = 0.06,p = 0.32 and ρ = 0.08,p = 0.14, respectively; n = 312). If translation controls noise in LP genes, then noise should also covariate with factors influencing translation efficiency such as ORF length or codon bias (see, for instance,[27, 28]). In LP genes, translational efficiency correlated more strongly with ORF length (ρ = −0.58,p = 6.7 × 10−31,n = 327) than with frequency of optimal codons (FOP, ρ = 0.31,p = 9.9 × 10−9,n = 327). Consistently, noise correlated with ORF length (ρ = −0.18,p = 5.5 × 10−5,n = 513 and Figure4C) but not so with FOP (ρ = 0.07,p = 0.12,n = 513). On the other hand, noise correlated with ORF length in an opposite way in HP genes (ρ = 0.20,p = 5.6 × 10−6,n = 309) which probably reflects complementary constraints on gene length (e.g., low noise genes in the HP class are mostly ribosomal genes, see Additional file1: Figure S9. These genes may exhibit short length due to minimization of biosynthetic costs given their high expression[29]).

Genomic neighborhood modulates the noise-plasticity coupling

What other mechanisms could modulate the noise-plasticity coupling? The specific architecture of the genomic neighborhood of a gene appears as a possible candidate. Indeed, it was recently found evidence for how bidirectional promoters (bipromoters) could reduce noise by favoring nucleosome depletion[16, 17]. We confirm this result for the full dataset (mean noise bipromoter = 0.0796, not bipromoter = 0.0965, Wilcoxon p = 1.2 × 10−8) and validate as well that potentially noise-sensitive gene classes are enriched in bipromoters: essential (465/1062, 43.8%) compared to nonessential (1651/4577, 36.1% p = 3.4 × 10−6, Fisher’s Exact test) and genes coding for protein complex subunits (666/1565, 42.6%) compared to the rest (1486/4195, 35.4% p = 8.1 × 10−7, Fisher’s Exact test).

However, as we have shown, noise in low-plasticity genes is modulated mostly at the translational level, and consequently should not be affected by the presence of bipromoters. Although (nucleosome) depletion is observed independently of plasticity, it does not seem to affect noise in the LP class. In contrast, the noise-reduction effect of bipromoters is observed in HP genes, as expected[16] (Additional file3: Table S2). This is further corroborated by the enrichment of bipromoters observed in LNHP genes (27/66, 41%) compared to HNHP (47/236, 20% p = 1.0 × 10−3, Fisher’s Exact test). Consistently, this enrichment is not significant in LNLP genes (75/145, 52%) compared to HNLP (38/67, 57%, Fisher’s Exact test, p = 0.6). This pattern is confirmed in noise-sensitive genes: while they are enriched in bipromoters in the HP group (79/237, 33.3%, compared to 67/274, 24.5% in non noise-sensitive, i.e., noise-tolerant, p = 0.03, Fisher’s Exact test), in the rest (low and medium plasticity) there are no observable differences in bipromoter frequency (noise-sensitive 348/745, 46.7%; and noise-tolerant 348/789, 44.1% p = 0.33, Fisher’s Exact test). Notably, low and medium plasticity (noise-sensitive) genes do exhibit a difference in noise (noise-sensitive genes conform 257/442, 58.1% of the LN, but only 109/252, 41.6% of the HN, p = 2.4 × 10−5, Fisher’s Exact test). If not to bipromoters, this difference could be attributed, as we discussed, to differences in ORF length (mean length noise-sensitive: 1732.0, noise-tolerant: 1517.2, p = 1.9 × 10−4, Wilcoxon test).

Reducing noise by bipromoters could additionally decrease expression plasticity due to their association to short intergenic distance and nucleosome depletion, and this we actually distinguished (Additional file3: Table S2). This suggests then a limitation on the adequacy of bipromoters for reducing noise. Interestingly, we detect a strong bipromoter-independent effect in noise-sensitive HP genes (Additional file4: Table S3). We thus hypothesized that there could be a tendency to evolve noise-abating mechanisms that affect plasticity more weakly, thereby uncoupling it from noise. Indeed, if we consider only bipromoter HP genes, we find that noise-sensitive ones tend to be TATAless (80%,n = 79) compared to noise-tolerant (61%,n = 67, p = 0.017, Fisher’s Exact test). In agreement with this, the noise-sensitive group has significantly lower noise (mean noise in bipromoter noise-sensitive = 0.083, mean noise in bipromoter noise-tolerant = 0.124, p = 1.0 × 10−4, Wilcoxon test) while the difference in plasticity is not significant (mean plasticity in bipromoter noise-sensitive = 0.105, mean plasticity in bipromoter noise-tolerant = 0.12, p = 0.62, Wilcoxon test) indicating an effective uncoupling of noise from plasticity. Can we identify further features illustrating that noise-sensitive bipromoter genes tend to maintain plasticity levels? Intergenic distances suggest that this could be the case (234 bp for noise-sensitive bipromoters and 190 bp noise-tolerant bipromoters, p = 0.019, Wilcoxon test).

Noncoding transcripts and modulation

The above can be complementary analyzed if we consider all possible local genomic architectures around a focal gene (Figure5A), i.e., parallel, divergent and bipromoters with a coding or non-coding transcript as upstream partner (noncoding partners include “cryptic unstable transcripts”, CUTs, and “stable untranslated transcripts”, SUTs, see[30] and Methods; bipromoters CUTS were recently associated with low noise[16]). We computed the coupling between noise and plasticity for each architecture. Coupling is strong for genes with divergent transcripts (independent of the type of upstream partner) and weak for those with a bipromoter with a coding partner (Figure5A). This further validates the observed absence of bipromoters in HNHP genes and their enrichment in the other three classes (bipromoters are the most commonly found architecture in LNLP, HNLP and LNHP) where they are associated, of course, to short intergenic distances (Figure5A, see also Additional file1: Figure S10). Interestingly, bipromoters of plastic genes with low noise are the ones with the biggest (relative) intergenic distance (with respect to LNLP and HNLP), which suggests again the requirement of a minimal distance to locate the regulatory demands associated to enhance plasticity (mean distance bipromoters of LNHP: 252 bp, in the LNLP and HNLP groups: 178bp, p = 1.1 × 10−3, Wilcoxon test). Overall, this emphasizes bipromoters as noise-abating architecture only when noise and plasticity are transcriptionally modulated.

Discussion

We analyzed the molecular determinants that adjust the linkage between gene expression plasticity and noise in S. cerevisiae. Noise was confirmed to be connected to plasticity when genes exhibit particular modes of transcription initiation (and re-initiation) related to the presence of TATA boxes at the promoter and strong chromatin regulation[5, 710, 12]. This could suggest a model in which intrinsic noise is a byproduct of the need for plasticity[14]. We show, however, that noisy expression can be observed in genes with low plasticity. These genes are generally simple (poor in transcription factor binding sites and generally TATAless), small and poorly regulated by chromatin effectors, and they show a prominent nucleosome depleted region (Additional file1: Figure S6, see also Figure5B).

In this loose regulatory scenario, transcription is likely to be produced by single, isolated in time, initiation events[22]. For such bacteria-like transcription, noise is anticipated to depend on translational efficiency[24, 25] (but see[31]) and this is indeed what we notice. In addition, ORF length appears as a strong determinant of ribosome occupancy in this class and thus of translation efficiency (potentially due to the lack of post-transcriptional regulation, see Additional file1: Figure S11)[32]. Following this model we expect essential genes –usually of low plasticity– to be large as we observe (size essential genes: 1646 bp, size nonessential: 1468 bp, p < 0.001, Wilcoxon test; see also Figure4).

Figure 4
figure 4

Noise determinants in low plasticity genes. Noise in LP genes is related to translational efficiency, which in turn is related to ORF length. We ordered LP genes by increasing noise. We performed a sliding window analysis of translational efficiency (A), ribosomal density (B) and ORF length (C). Shaded regions represent the mean and two standard deviations at each point obtained with the same sliding window analysis over randomized data; the process was repeated 10000 times. See also main text.

On the other hand, we also established that relatively high expression plasticity can be reached by chromatin regulation without necessarily coupling it to noise. These genes that are enriched in growth-related functions (such as ribosomal protein genes, RPs) depend strongly on TAF1, an essential subunit of TFIID, but less strongly on most of the more specific chromatin regulators (as compared to noisy and plastic genes in Figure2, see further discussion in Additional file1: Supplement). This indicates that they respond to general, rather than gene specific, regulatory strategies which partly explain the high degree of co-regulation previously observed[13]. Beyond this, the group exhibits a characteristic pattern of low nucleosome occupancy in both proximal and distal promoter regions[8] possibly caused by the strong enrichment in activating histone modifications[13], and particularly acetylations. We hypothesize that rather than a promoter-localized open-chromatin state these genes could be located at broader open chromatin domains. Indeed, we distinguished that RP genes tend to be localized on broad open-chromatin domains that extend up to at least 40Kb (Additional file1: Figure S12), supporting this view.

The low noise but highly plastic expression is therefore consistent with two (not mutually exclusive) models previously proposed. Firstly, it is consistent with a detailed model[33] in which low nucleosome occupancy at the promoter indicates a stable open state, allowing the high expression levels exhibited by these genes (Additional file3: Figure S13). The concomitant noise reduction would not be possible if the high expression level would be reached by an increase in transcription or translation efficiency[33]. Secondly, at a broader level, the localization in open chromatin, and consequently low noise, genomic regions could also contribute to the small level of noise detected[15].

Finally, we determined how additional enhancement in plasticity fundamentally associates to strong intrinsic noise. This emphasizes two additional promoter features. Firstly, a particular nucleosome located at the promoter[79] allows a fine-tuned –and gene-specific– control of open and close promoter states by external regulators. The repressive effect of this nucleosome is evidenced by the increase in expression level upon histone depletion, an increase only observed in these genes (Figure3). We demonstrate that many specific chromatin regulators act on these genes, in contrast to their low-noise counterpart (Figure2) that could reflect also in the lack of co-regulation reported[13]. Additionally, during the time lapses that the promoter is in open state, the presence of a TATA box allows the pre-initiation complex to stay assembled firing continuous initiation events. This increases the sensitivity to changes –in the time the promoter stays in open state produced by chromatin remodeling– thereby allowing an increase in plasticity. Moreover, and even when repressed, the nucleosome can occasionally be destabilized allowing strong bursts of transcription which result in the observed noise. That coupling is related to an efficient transcription initiation is confirmed by the strong coupling found in SAGA but TATAless genes (Additional file1: Figure S3), what confirms the model in[33] to a genome-wide scale.

As the critical promoter-covering nucleosome is probably stabilized by particular DNA properties, such as high bendability[9, 34], this could potentially increase the number of phenotype-affecting mutational targets, which could be in turn the cause of the increased expression divergence in these genes (Additional file1: Figure S14). Therefore, linkage between different types of variability is mechanistically a consequence of the sophisticated regulatory strategy involving promoter nucleosomes and TATA boxes. This regulation also brings higher sensibility to chromatin regulation (leading to plasticity), to stochastic nucleosome fluctuations (leading to noise) and to mutational effects (leading to expression divergence).

The action of these two distinct strategies to modulate noise and plasticity coupling is further emphasized by the structure of the genomic neighborhood of the focal gene under study. While a seemingly general architecture in which (relative) reduction of intergenic distance and enrichment of bipromoters should indicate noise-plasticity uncoupling, this only applies to transcriptional-based modulation (Figure5A). Indeed, both small integenic distance and bipromoter can be broadly related to high or low noise in poorly plastic genes, where uncoupling is rather associated to other mechanisms of modulating noise (HNLP genes exhibit higher translational efficiency). Moreover, we observe that genes particularly sensitive to noise (e.g., genes specifying proteins in complexes) can separate noise from the requirement of high plasticity (these genes showed larger intergenic distances as bipromoters).

In view of the exposed implications, can we speculate about the evolution of bipromoters? Firstly, the S. cerevisiae genome is highly gene-dense, averaging one gene each 2 Kb with a median intergene distance in our dataset (considering noncoding transcripts) of 204 bp. A high bipromoter frequency seems then plausible in absence of selection. In addition, we observed a correlation between intergene distance and plasticity (ρ = 0.19,p < 10−43,n = 5102), stronger when only transcripts with divergently oriented upstream partners are considered (ρ = 0.27,p < 10−57,n = 3271; Figure5A). This relationship probably responds to the need of a greater genomic space to accommodate a more complex regulatory landscape, which is in turn needed to achieve controlled expression variability. We observed a strong bias in bipromoter frequency towards genes with low plasticity (Additional file1: Figure S15), indicating that these regulatory needs are a major force determining the absence of bipromoters. In contrast, we suggest that their presence could have an almost neutral origin, as it displays an expectable distribution in view of intergenic distances (Additional file1: Figure S16).

Conclusions

In sum, the results reported here reveal that transcriptional- and translational-based regulatory strategies are alternatively used to modulate the balance between noise and plasticity in eukaryotic gene expression (Figure5). These strategies appear clearly associated to distinctive functional (e.g., growth/stress programs in S. cerevisiae[12, 13, 20]), and genomic constrains (e.g., presence of bipromoters, non-coding transcription or length of coding sequence). Future analysis of additional questions, e.g., role of post-transcriptional regulation (Additional file1: Figure S11), potential presence of condition-dependent variation[32], or level and relevance of coupling in higher eukaryotes, should ultimately expose the many aspects of gene expression variation and its evolution.

Figure 5
figure 5

Noise plasticity coupling is modulated by genomic neighborhood and distinguishes four control strategies overall. A). A cartoon depicting the different genomic structures (bipromoter, parallel, divergent) upstream of coding genes is shown in ascending order of proximal nucleosome occupancy, plasticity and noise (which coincide). For each structure, we show the average intergenic distance in blue. In red is shown the Spearman ρ coefficient for the observed noise-plasticity correlation. We also show the percent within each class of a given upstream structure , e.g., HNHP mostly exhibit parallel/divergent coding (C) and divergent non-coding (NC) transcripts. B) Four regulatory strategies broadly adjust the noise-plasticity coupling. These strategies emphasize the alternative transcriptional- or translational-based modes of balancing noise and plasticity in yeast.

Methods

Gene expression variability

To quantify gene expression plasticity, i.e. responsiveness to environmental change, the variability in mRNA levels among >1500 different growth conditions was measured. as the sum of squares of the log2-ratios over all these conditions[4]. Noise, or stochastic variability, was measured by proteomic analysis[3]. We used the “distance to the median”(DM) score, which rules out confounding effects of protein abundance, allowing protein-specific noise levels to be compared. Evolutionary divergence in gene expression was measured as the variation of gene expression between orthologs in in four related yeast species and 32 different conditions[4]. These three gene expression variability measures we scaled between 0 and 1. After scaling, mean values are 0.062, 0.089, and 0.186 for plasticity, noise and expression divergence, respectively. Moreover, for the noise and plasticity measures, we define three categories “high”, “medium”and“low”, using percentiles 25 and 75 in each case as boundaries. Thus we obtain groups of genes with high noise and high plasticity (HNHP), high noise and low plasticity (HNLP) and so on (see also Additional file5: Table S4). Finally, mRNA level in rich media was obtained from[35] (mean=3.915 mRNA copies/cell).

Genomic localization and neighborhood

The coordinates of each transcript (coding ORFs–ORF-T; and noncoding, which can be in turn “cryptic unstable transcripts”–CUTs, and “stable untranslated transcripts”–SUTs) were obtained from a high resolution transcriptomic analysis[30]. For each of these transcripts, we used chromosomic coordinates of transcription start sites (TSS) and transcription end sites (TES), and orientation (strand). This data allowed us to characterize the genomic neighborhood of each transcript, in terms of distance to its upstream partner in bp, orientation of this upstream partner (which can be divergent or parallel). As well, from[30] we obtained data describing for each gene whether it is transcribed from a bi-directional promoter based on the existence of a shared nucleosome depleted region (NDR). For the genomic neighborhood analysis, in order to maximize its reliability, we removed from the dataset gene whose upstream partner was a “pseudogene” or a “dubious ORF”, as well as a few confounding cases where adjacent transcripts were overlapping. For some genes upstream distance could not be calculated as TSS and/or TES coordinates could not be accurately determined in the original source (see[30]).

Promoter characterization and regulation

The presence/absence of TATA boxes at the promoters was obtained from[12]. Nucleosome occupancy data for the whole genome was obtained by DNA digestion with micrococcal nuclease and identification of nucleosome-protected fragments by high resolution microarray analysis[7]. We use the log2-ratios provided in the reference. As suggested in[8], we obtained two different nucleosome occupancy values for each promoter. Taking as reference the TSS, proximal nucleosome occupancy was the average in the -100 to 0bp region, while distal nucleosome occupancy corresponds to the -400 to -150 bp region. For an idea, the highest occupancy for a proximal region in our dataset was 0.27, and the lowest -3.64.

Transcription regulation data

To explore chromatin regulation, we used a compendium, assembled in[18], consisting of 170 expression profiles for chromatin regulation related mutations (expressed in log2-ratios). We classified these mutations in three classes (see Additional file2: Table S1).“Chromatin” tag was assigned to mutations in histone acetyltransferases, deacetyltransferases, methylases, demethylases, ubiquitinating and deubiquitinating enzymes, chromatin remodellers and silencing factors. “General” tag was assigned to genotypes involving at least one mutation in essential, general transcription factors (TAF1). “Histone” tag was assigned to mutations in the very histones. As suggested in[9], we normalized each dataset from the compendium to unit variance. The absolute value of the normalized log2 ratios represented responsiveness measures; the mean responsiveness of each gene represented its “chromatin regulation effect” (CRE) or “histone regulation effect” (HRE). We used also this normalized dataset without taking the absolute value to analyze the sense of the observed regulation. Data for nucleosome-normalized, chromatin modification states at promoter were obtained from ChromatinDB (http://www.bioinformatics2.wsu.edu/chromatindb) which unifies several experimental genome-wide datasets measuring levels of different histone modifications. For dependence of each gene on general TFs, we used categorical data from[22] defining for the expression of each gene if it is dominated by TFIID or SAGA complex.

Translation related measures

We used a measure of translation efficiency obtained in[26] and based for each gene in percent of each transcript in polysomes, its ribosome density, and the relative transcript level (mean=4.35, sd=1.72). We used an additional dataset of ribosome density obtained from[27] (mean=0.53, sd=0.31).

Noise-sensitive genes

We considered essential genes from the Saccharomyces Genome Deletion Project and genes specifying proteins in complexes[36]. Due to a big reduction in sample size, we excluded from this group haploinsufficient genes[37]. However, note that they are virtually not excluded, since 44/46 (95%) of the identified haploinsufficient genes are labeled as either essential or as complex-forming; indeed, 100% of the haploinsufficient genes located in the HP group are so.

Author’s contributions

JFP conceived the project, DB conducted computational experiments, DB and JFP analyzed experiments and wrote the manuscript. Both authors read and approved the final manuscript.

References

  1. Enard W, Khaitovich P, Klose J, Zöllner S, Heissig F, Giavalisco P, Nieselt-Struwe K, Muchmore E, Varki A, Ravid R, Doxiadis GM, Bontrop RE, Pääbo S: Intra- and interspecific variation in primate gene expression patterns. Sci (New York, N.Y.). 2002, 296 (5566): 340-343. 10.1126/science.1068996. [http://www.sciencemag.org/content/296/5566/340.abstract],

    Article  CAS  Google Scholar 

  2. Fraser HB, Hirsh AE, Giaever G, Kumm J, Eisen MB: Noise minimization in eukaryotic gene expression. PLoS Biol. 2004, 2 (6): e137-10.1371/journal.pbio.0020137. [http://www.ncbi.nlm.nih.gov/pubmed/15124029],

    Article  PubMed Central  PubMed  Google Scholar 

  3. Newman JRS, Ghaemmaghami S, Ihmels J, Breslow DK, Noble M, DeRisi JL, Weissman JS: Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise. Nat. 2006, 441 (7095): 840-846. 10.1038/nature04785. [http://www.ncbi.nlm.nih.gov/pubmed/16699522],

    Article  CAS  Google Scholar 

  4. Tirosh I, Weinberger A, Carmi M, Barkai N: SUPL-A genetic signature of interspecies variations in gene expression. Nat Genet. 2006, 38 (7): 830-834. 10.1038/ng1819. [http://www.ncbi.nlm.nih.gov/pubmed/16783381],

    Article  CAS  PubMed  Google Scholar 

  5. Blake WJ, Kaern M, Cantor CR, Collins JJ: Noise in eukaryotic gene expression. Nat. 2003, 422 (6932): 633-637. 10.1038/nature01546. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=400249&tool=pmcentrez&rendertype=abstract],

    Article  CAS  Google Scholar 

  6. Landry CR, Lemos B, Rifkin Sa, Dickinson WJ, Hartl DL: Genetic properties influencing the evolvability of gene expression. Sci(New York, N.Y.). 2007, 317 (5834): 118-121. 10.1126/science.1140247. [http://www.ncbi.nlm.nih.gov/pubmed/17525304],

    Article  CAS  Google Scholar 

  7. Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow C: high-resolution atlas of nucleosome occupancy in yeast. Nat Genet. 2007, 39 (10): 1235-1244. 10.1038/ng2117. [http://www.ncbi.nlm.nih.gov/pubmed/17873876],

    Article  CAS  PubMed  Google Scholar 

  8. Tirosh I, Barkai N: Two strategies for gene regulation by promoter nucleosomes. Genome Res. 2008, 18 (7): 1084-1091. 10.1101/gr.076059.108. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2493397&tool=pmcentrez&rendertype=abstract],

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Choi JK, Kim YJ: Epigenetic regulation and the variability of gene expression. Nat Genet. 2008, 40 (2): 141-147. 10.1038/ng.2007.58. [http://www.ncbi.nlm.nih.gov/pubmed/18227874],

    Article  CAS  PubMed  Google Scholar 

  10. Choi JK, Kim YJ: Intrinsic variability of gene expression encoded in nucleosome positioning sequences. Nat Genet. 2009, 41 (4): 498-503. 10.1038/ng.319. [http://www.ncbi.nlm.nih.gov/pubmed/19252489],

    Article  CAS  PubMed  Google Scholar 

  11. Warner JR: The economics of ribosome biosynthesis in yeast. Trends Biochem Sci. 1999, 24 (11): 437-440. 10.1016/S0968-0004(99)01460-7. [http://www.ncbi.nlm.nih.gov/pubmed/10542411],

    Article  CAS  PubMed  Google Scholar 

  12. Basehoar AD, Zanton SJ, Pugh BF: Identification and distinct regulation of yeast TATA box-containing genes. Cell. 2004, 116 (5): 699-709. 10.1016/S0092-8674(04)00205-3. [http://www.ncbi.nlm.nih.gov/pubmed/15006352],

    Article  CAS  PubMed  Google Scholar 

  13. Kim SC, Choi JK: Controlling transcriptional programs for cellular adaptation by chromatin regulation. Mol BioSyst. 2011, 7 (5): 1713-9. 10.1039/c0mb00319k. [http://www.ncbi.nlm.nih.gov/pubmed/21412562],

    Article  CAS  PubMed  Google Scholar 

  14. Lehner B: Conflict between noise and plasticity in yeast. PLoS Genet. 2010, 6 (11): e1001185-10.1371/journal.pgen.1001185. [http://dx.plos.org/10.1371/journal.pgen.1001185],

    Article  PubMed Central  PubMed  Google Scholar 

  15. Batada NN, Hurst LD: Evolution of chromosome organization driven by selection for reduced gene expression noise. Nat Genet. 2007, 39 (8): 945-949. 10.1038/ng2071. [http://www.ncbi.nlm.nih.gov/pubmed/17660811],

    Article  CAS  PubMed  Google Scholar 

  16. Wang GZ, Lercher MJ, Hurst LD: Transcriptional coupling of neighboring genes and gene expression noise: evidence that gene orientation and noncoding transcripts are modulators of noise. Genome Biol Evol. 2011, 3: 320-331. 10.1093/gbe/evr025. [http://www.ncbi.nlm.nih.gov/pubmed/21402863],

    Article  CAS  PubMed  Google Scholar 

  17. Woo YH, Li WH: Gene clustering pattern, promoter architecture, and gene expression stability in eukaryotic genomes. Proc National Acad Sci U S A. 2011, 108 (8): 3306-3311. 10.1073/pnas.1100210108. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3044413&tool=pmcentrez&rendertype=abstract],

    Article  CAS  Google Scholar 

  18. Steinfeld I, Shamir R, Kupiec M: A genome-wide analysis in Saccharomyces cerevisiae demonstrates the influence of chromatin modifiers on transcription. Nat Genet. 2007, 39 (3): 303-309. 10.1038/ng1965.

    Article  CAS  PubMed  Google Scholar 

  19. Zenklusen D, Larson DR, Singer RH: Single-RNA counting reveals alternative modes of gene expression in yeast. Nat Struct Mol Biol. 2008, 15 (12): 1263-1271. 10.1038/nsmb.1514. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3154325&tool=pmcentrez&rendertype=abstract],

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. López-Maury L, Marguerat S, Bähler J: Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nat Rev Genet. 2008, 9 (8): 583-593. [http://www.ncbi.nlm.nih.gov/pubmed/18591982],

    Article  PubMed  Google Scholar 

  21. Angus-Hill ML, Schlichter A, Roberts D, Erdjument-Bromage H, Tempst P, Cairns BR: A Rsc3/Rsc30 zinc cluster dimer reveals novel roles for the chromatin remodeler RSC in gene expression and cell cycle control. Mol cell. 2001, 7 (4): 741-751. 10.1016/S1097-2765(01)00219-2. [http://www.ncbi.nlm.nih.gov/pubmed/11336698],

    Article  CAS  PubMed  Google Scholar 

  22. Huisinga KL, Pugh BF: A genome-wide housekeeping role for TFIID and a highly regulated stress-related role for SAGA in Saccharomyces cerevisiae. Mol Cell. 2004, 13 (4): 573-585. 10.1016/S1097-2765(04)00087-5. [http://www.ncbi.nlm.nih.gov/pubmed/14992726],

    Article  CAS  PubMed  Google Scholar 

  23. Durant M, Pugh BF: Genome-Wide Relationships between TAF1 and Histone Acetyltransferases in Saccharomyces cerevisiae. Soc. 2006, 26 (7): 2791-2802.

    CAS  Google Scholar 

  24. Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A: Regulation of noise in the expression of a single gene. Nat Genet. 2002, 31: 69-73. 10.1038/ng869. [http://www.ncbi.nlm.nih.gov/pubmed/11967532],

    Article  CAS  PubMed  Google Scholar 

  25. Raser JM, O’Shea EK: Noise in gene expression: origins, consequences, and control. Sci (New York, N.Y.). 2005, 309 (5743): 2010-2013. 10.1126/science.1105891. [http://www.ncbi.nlm.nih.gov/pubmed/16179466],

    Article  CAS  Google Scholar 

  26. MacKay VL, Li X, Flory MR, Turcott E, Law GL, Serikawa KA, Xu XL, Lee H, Goodlett DR, Aebersold R, Zhao LP, Morris DR: Gene expression analyzed by high-resolution state array analysis and quantitative proteomics: response of yeast to mating pheromone. Mol Cell Proteomics : MCP. 2004, 3 (5): 478-489. 10.1074/mcp.M300129-MCP200. [http://www.ncbi.nlm.nih.gov/pubmed/14766929],

    Article  CAS  PubMed  Google Scholar 

  27. Arava Y, Wang Y, Storey JD, Liu CL, Brown PO, Herschlag D: Genome-wide analysis of mRNA translation profiles in Saccharomyces cerevisiae. Proc National Acad Sci U S A. 2003, 100 (7): 3889-3894. 10.1073/pnas.0635171100. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=153018&tool=pmcentrez&rendertype=abstract],

    Article  CAS  Google Scholar 

  28. Lackner DH, Beilharz TH, Marguerat S, Mata J, Watt S, Schubert F, Preiss T, Bähler J: A network of multiple regulatory layers shapes gene expression in fission yeast. Mol Cell. 2007, 26: 145-155. 10.1016/j.molcel.2007.03.002. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1885965&tool=pmcentrez&rendertype=abstract],

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Warringer J, Blomberg A: Evolutionary constraints on yeast protein size. BMC Evolutionary Biol. 2006, 6: 61-10.1186/1471-2148-6-61. [http://www.ncbi.nlm.nih.gov/pubmed/16911784],

    Article  Google Scholar 

  30. Xu Z, Wei W, Gagneur J, Perocchi F, Clauder-Münster S, Camblong J, Guffanti E, Stutz F, Huber W, Steinmetz LM: Bidirectional promoters generate pervasive transcription in yeast. Nat. 2009, 457 (7232): 1033-1037. 10.1038/nature07728. [http://www.ncbi.nlm.nih.gov/pubmed/19169243],

    Article  CAS  Google Scholar 

  31. Golding I, Paulsson J, Zawilski SM, Cox EC: Real-time kinetics of gene activity in individual bacteria. Cell. 2005, 123 (6): 1025-1036. 10.1016/j.cell.2005.09.031. [http://www.ncbi.nlm.nih.gov/pubmed/16360033],

    Article  CAS  PubMed  Google Scholar 

  32. Lee MV, Topper SE, Hubler SL, Hose J, Wenger CD, Coon JJ, Gasch AP: A dynamic model of proteome changes reveals new roles for transcript alteration in yeast. Mol Syst Biol. 2011, 7: 514-[http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3159980&tool=pmcentrez&rendertype=abstract],

    Article  PubMed Central  PubMed  Google Scholar 

  33. Raser JM, O’Shea EK: Control of stochasticity in eukaryotic gene expression. Sci (New York, N.Y.). 2004, 304 (5678): 1811-1814. 10.1126/science.1098641. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1410811&tool=pmcentrez&rendertype=abstract],

    Article  CAS  Google Scholar 

  34. Tirosh I, Berman J, Barkai N: The pattern and evolution of yeast promoter bendability. Trends Genet: TIG. 2007, 23 (7): 318-321. 10.1016/j.tig.2007.03.015. [http://www.ncbi.nlm.nih.gov/pubmed/17418911],

    Article  CAS  PubMed  Google Scholar 

  35. Holstege FC, Jennings EG, Wyrick JJ, Lee TI, Hengartner CJ, Green MR, Golub TR, Lander ES, Young Ra: Dissecting the regulatory circuitry of a eukaryotic genome. Cell. 1998, 95 (5): 717-28. 10.1016/S0092-8674(00)81641-4. [http://www.ncbi.nlm.nih.gov/pubmed/9845373],

    Article  CAS  PubMed  Google Scholar 

  36. Wang H, Kakaradov B, Collins SR, Karotki L, Fiedler D, Shales M, Shokat KM, Walther TC, Krogan NJ, Koller D: A complex-based reconstruction of the Saccharomyces cerevisiae interactome. Mol cell Proteomics : MCP. 2009, 8 (6): 1361-1381. 10.1074/mcp.M800490-MCP200. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2690481&tool=pmcentrez&rendertype=abstract],

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Deutschbauer AM, Jaramillo DF, Proctor M, Kumm J, Hillenmeyer ME, Davis RW, Nislow C, Giaever G: Mechanisms of haploinsufficiency revealed by genome-wide profiling in yeast. Genet. 2005, 169 (4): 1915-1925. 10.1534/genetics.104.036871. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1449596&tool=pmcentrez&rendertype=abstract],

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank Laurence D. Hurst for comments on an earlier version of the manuscript. This research was partially supported by the Spanish Ministerio de Economía y Competitividad BFU2011-24691 grant (JFP).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Juan F Poyatos.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

12864_2012_4554_MOESM1_ESM.pdf

Additional file 1: Supplement. Supplementary material including additional discussions on Taf1 and Tup1 effects and figures S1 to S16. (PDF 284 KB)

12864_2012_4554_MOESM2_ESM.xls

Additional file 2: Table S1. Additional table with detailed analysis of chromatin regulators associated to Figure2. (XLS 52 KB)

12864_2012_4554_MOESM3_ESM.xls

Additional file 3: Table S2. Additional table with analysis of the effect of bipromoters on proximal nucleosomal occupancy, noise and plasticity. (XLS 10 KB)

12864_2012_4554_MOESM4_ESM.xls

Additional file 4: Table S3. Additional table with analysis of bipromoter enrichment on noise sensitive genes. (XLS 20 KB)

12864_2012_4554_MOESM5_ESM.xls

Additional file 5: Table S4. Additional table with functional enrichment analysis on the different genes according to noise-plasticity coupling. (XLS 163 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Bajić, D., Poyatos, J.F. Balancing noise and plasticity in eukaryotic gene expression. BMC Genomics 13, 343 (2012). https://doi.org/10.1186/1471-2164-13-343

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-343

Keywords