Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Research article

Protein stickiness, rather than number of functional protein-protein interactions, predicts expression noise and plasticity in yeast

Leandra M Brettner12 and Joanna Masel1*

Author Affiliations

1 Present address: Ecology & Evolutionary Biology, University of Arizona, 1041 E Lowell St, Tucson, AZ, 85721, USA

2 Present address: Bioengineering, University of Washington, 3720 15th Ave NE, Seattle, WA, 98195, USA

For all author emails, please log on.

BMC Systems Biology 2012, 6:128  doi:10.1186/1752-0509-6-128


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/6/128


Received:16 July 2012
Accepted:21 September 2012
Published:27 September 2012

© 2012 Brettner and Masel; licensee 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.

Abstract

Background

A hub protein is one that interacts with many functional partners. The annotation of hub proteins, or more generally the protein-protein interaction “degree” of each gene, requires quality genome-wide data. Data obtained using yeast two-hybrid methods contain many false positive interactions between proteins that rarely encounter each other in living cells, and such data have fallen out of favor.

Results

We find that protein “stickiness”, measured as network degree in ostensibly low quality yeast two-hybrid data, is a more predictive genomic metric than the number of functional protein-protein interactions, as assessed by supposedly higher quality high throughput affinity capture mass spectrometry data. In the yeast Saccharomyces cerevisiae, a protein’s high stickiness, but not its high number of functional interactions, predicts low stochastic noise in gene expression, low plasticity of gene expression across different environments, and high probability of forming a homo-oligomer. Our results are robust to a multiple regression analysis correcting for other known predictors including protein abundance, presence of a TATA box and whether a gene is essential. Once the higher stickiness of homo-oligomers is controlled for, we find that homo-oligomers have noisier and more plastic gene expression than other proteins, consistent with a role for homo-oligomerization in mediating robustness.

Conclusions

Our work validates use of the number of yeast two-hybrid interactions as a metric for protein stickiness. Sticky proteins exhibit low stochastic noise in gene expression, and low plasticity in expression across different environments.

Keywords:
Protein-protein interaction networks; Stochastic gene expression; Evolutionary constraint; Correlomics; Cooperativity; Phenotypic plasticity

Background

A protein that functionally interacts with many other proteins may be more sensitive to noise in gene expression [1]. In agreement with this prediction, a negative correlation between noise and protein-protein interaction (PPI) degree has been found [2,3]. However, PPI datasets are notorious for high rates of false positive and false negative interactions [4-7]. Older high throughput datasets rely on yeast two-hybrid (Y2H) studies, which can measure interactions between two proteins that would never even encounter each other in nature. More recently, high throughput affinity capture mass spectrometry (ACMS) data have become available [8], which do not suffer from this drawback and subsequent high false positive rate.

Y2H data may indicate the non-specific “stickiness” of a protein towards a random polypeptide better than it indicates the number of functional protein-protein interactions that the protein is involved in [9,10]. The numbers of PPIs per protein (node degrees) for Y2H vs. ACMS data are only weakly correlated in yeast (Figure 1, R2 = 0.008, p = 2e-05). Y2H and ACMS data clearly measure different things.

thumbnailFigure 1. A protein’s number of PPIs correlates poorly across two high throughput data types. Model 1 regression line is shown for illustrative purposes only, to show the weakness of the correlation.

Here we find that Y2H degree/stickiness predicts gene noise better than ACMS/“real PPI” degree does. It has previously been argued that a protein that functionally interacts with many other proteins will be more sensitive to noise [1]. However, if a protein binds promiscuously to many non-functional partners, variable expression of that protein may be extremely disruptive to diverse processes, also leading to a negative correlation, in this case between the number of false positive PPIs and noise. By itself, ACMS degree is correlated with noise, perhaps due to residual sticky but non-functional false positives in the ACMS data. But ACMS drops out as a statistically significant predictor in a multiple regression when Y2H degree, a better estimator of non-specific protein-protein binding, is included.

Here we also introduce a new metric of plasticity, i.e. the variation in gene expression across different experiments. Again, we find that Y2H degree/stickiness predicts plasticity better than higher quality ACMS data on the number of functional PPIs. For both noise and plasticity, our results are robust to a multiple regression analysis that controls for protein abundance and that controls noise for plasticity and vice versa. Other significant predictors include the presence of a TATA box, whether a protein forms a homo-oligomer, and gene essentiality.

Results

Table 1 summarizes the results of regression analyses predicting noise. After extensive bottom-up and top-down multiple regression model building, the best model (first numeric column: total coefficient of determination R2 = 0.1083) included Y2H, but not ACMS, as a predictor of noise. The second numeric column shows the R2 when a single predictor is considered in isolation. These values are given as a contrast to our primary results in the first numeric column, which show the extent to which the coefficient of determination is reduced if a predictor is removed from the best model. For the purposes of more graphical illustration, Figure 2A also shows the effect of each binary predictor in isolation.

Table 1. Multiple regression results predicting noise

thumbnailFigure 2. Illustration of binary predictors of noise and plasticity, taken in isolation. The presence of a TATA box strongly predicts noise and plasticity. Homo-oligomerization does not, in isolation, predict plasticity, and its effect on noise is only marginally statistically significant (p = 0.0496). However, these effects become significant when confounding factors are accounted for (Tables 1 and 2). Essentiality predicts noise but not plasticity. To better assess effect sizes using more intuitive noise and plasticity measures, back transformations were performed to restore original units. The mean plasticity residual was added to the mean Box-Cox transformed plasticity score, and then the Box-Cox transform was reversed, so that plasticity corresponds simply to the estimated number of experiments for which expression varies. The noise axis corresponds to the DM metric of Newman et al. [13]. Error bars correspond to 95% confidence intervals.

The strongest predictor for noise is the presence of a TATA box, consistent with earlier findings: TATA boxes are associated with higher noise [11-13]. Gene essentiality is also an important predictor of gene noise, again consistent with earlier findings that essential genes have lower noise [1,2,13,14]. We also found a statistically significant interaction term, with genes that are both non-essential and possess a TATA box having higher noise than would be expected from the two factors in isolation.

Genes that interact with themselves (form homo-oligomers) have higher noise than genes that do not self-associate. Explanations for this novel finding are explored further at the end of the Results section and in the Discussion.

Genes with high noise in a single environment tend also to have high variation across different environmental conditions (plasticity), due at least in part to mechanistic coupling at the promoter level [2,11,15-18] (Figure 3). In order to infer variables that affect noise reliably, it is therefore important to correct for plasticity. Previous metrics of plasticity have been based on the average pair-wise ratio between microarray spot densities across a variety of environmental conditions [11,12,19,20]. However, the dynamic range of microarray signals depends on transcript abundance [21], making this plasticity metric dependent on abundance. Here we construct a plasticity metric that is less abundance-dependent by design, and which we then correct for residual effects of protein abundance (see Methods). Note that our estimate of protein noise has already been corrected for protein abundance [13]. It is important to correct noise and plasticity for abundance before testing their correlation with PPI degree, since PPI degree can be confounded with abundance [22]. Indeed, protein abundance is an important constraint on evolution, and so may affect a wide range of properties [23].

thumbnailFigure 3. A gene’s noise and plasticity are correlated. Multiple regression analyses in Tables 1 and 2 use Model 1 regression, but with reversed dependent and independent variables. For such a weak correlation, plasticity as a function of noise is quite different from the inverse function of noise as a function of plasticity: both lines are shown here. In the absence of a correlation, the functions describing these two lines would be horizontal and vertical, respectively. For comparison, the Model 2 Standard Major Axis regression line is also shown. The correlation between noise and plasticity is tighter in the top right corner, where values of both are high [11].

When we correct for plasticity, our main results on predictors of noise all still hold (Table 1, last two columns, total coefficient of determination R2 rises to 0.1609), although many of the R2 values attributable to specific predictors are modestly reduced. Correlations between noise and plasticity are known to be stronger in genes containing a TATA box [11]. In agreement with this, the TATA × Plasticity interaction term is statistically significant. Plasticity predicts noise both for TATA(+) genes and for TATA(−) genes, but the effect size (i.e. regression coefficient or slope) is 2.5 times as large for TATA(+) genes (Table 1).

Next, we considered how PPI metrics and other factors predict plasticity (Table 2, total coefficient of determination R2 =0.1267) for that subset of genes for which noise data were also available. Note that this requirement for the availability of noise data biases analyses towards the properties of higher-abundance proteins. Fortuitously, this makes ACMS a more reliable metric of “true” PPIs [5,24], strengthening our interpretation of the results.

Table 2. Multiple regression results predicting plasticity

As with low noise, high Y2H degree/stickiness predicts low plasticity but many ACMS/“functional” PPIs do not. This plasticity correlation holds true even after correcting for the effects of noise. This may be because promiscuous binding poses a particular challenge when it occurs at different extents in different environments. Or it may be because the successful fulfilment of the function of a plastic gene, which requires different levels of expression in different environments, is more sensitive to the effects of a given quantity of noise.

Gene essentiality predicts noise, but it does not predict plasticity. This could be because many genes are only essential in some environments or cell cycle stages, rather than constantly. In agreement with previous findings [11], the presence of a TATA box predicts plasticity as well as noise.

A novel finding of this paper is that the ability to homo-oligomerize predicts both noise and plasticity. Sticky proteins that bind promiscuously are also more likely to stick to themselves [25]. This means that self-interaction and Y2H are both surrogate metrics for intrinsic protein stickiness. In agreement with this interpretation of Y2H, but not ACMS, as a measure of “stickiness”, self-interaction correlates with Y2H but not ACMS (Figure 4).

thumbnailFigure 4. Proteins that homo-oligomerize are stickier, but do not have more functional PPIs. Analyses were performed on log(PPI) and back-transformed to yield more intuitive PPI metrics. 95% confidence intervals are shown.

However, in our multiple regressions, supposedly sticky self-interacting proteins had higher rather than lower noise and plasticity. In the Discussion, we explore possible causes of this relationship. The relationship can only be seen when stickiness is first controlled for, via Y2H data, in a multiple regression analysis. Deceptively, correlations between homo-oligomerization and noise or plasticity were weak to non-existent in single factor analyses (Table 1, Table 2, Figure 2), where stickiness is a confounding factor.

Unsurprisingly given that both correlate with stickiness, the Self-Interaction × Y2H interaction term is statistically significant in our predictive model of plasticity. Y2H predicts plasticity more strongly (i.e. with a larger slope/coefficient) for the already-sticky self-interacting proteins than for non-self-interacting proteins.

Discussion

Avoiding non-functional PPIs is an important constraint in protein evolution [9,26-28]. Use of the number of Y2H interactions as a validated metric of non-functional PPIs, or “stickiness”, has the promise to reveal more about the nature and consequences of this constraint. Here we have contributed to this validation by showing that Y2H degree is a better predictor of gene expression noise, plasticity, and likelihood of homo-oligomerization than the supposedly superior ACMS data on “true” protein-protein interactions. Given that Y2H data are known to be poor indicators of functional PPIs, our results imply that Y2H data can nevertheless yield a metric with real biological meaning.

We also found that proteins that homo-oligomerize had higher noise and higher plasticity, after confounding factors (including PPI [25]) were controlled for. High variation in protein abundance (noise) does not necessarily correspond linearly with high variation in protein activity. To explain our results, we hypothesize that homo-oligomerization decreases the sensitivity of protein activity to stochastic noise in protein abundance. Plastic genes, which require different levels of activity in different environments, may be more sensitive to the effects of a given quantity of noise, explaining why plasticity is also predicted by homo-oligomerization.

Two very different mechanisms may explain how homo-oligomerization decreases the sensitivity of protein activity to stochastic noise in protein abundance, depending on whether the active form of the protein in question is a monomer or a homo-oligomer. First, consider the case where the monomer is the active form. Homo-oligomerization may act as a sequestration sink that depends in a stronger than linear fashion on concentration. This creates robustness to noise by making the active monomer concentration less dependent on the total level of expression of that protein [29,30]. Sequestration via homo-oligomers rather than hetero-oligomers could help prevent concentration changes from cascading through the PPI network [31].

If the homo-oligomer is the active form, noise in protein abundance can be mitigated by switch-like kinetics (i.e. a sigmoidal dose–response curve) [32]. With a switch, increasing gene expression has a negligible effect until a critical threshold concentration is reached (Figure 5A). The response is then rapidly amplified until near saturation. Sigmoidal kinetics attenuate the effects of noise by allowing the cell to react only to stimuli of an adequate magnitude [32]. Sigmoidal kinetics control noise by controlling the level of activity, rather than by closely regulating the concentration of a signal molecule.

thumbnailFigure 5. Sigmoidal dose–response curves of cooperative proteins.A) In the shaded area, cooperativity suppresses the effects of gene expression noise, preventing inappropriate pathways from being switched on. B) Dose–response curves shown for Hill coefficients of 1, 2, 3, and 4.

Protein cooperativity is a common example of a switch-like system [33,34]. For example, the active form of a protein may consist of several subunits, each with a binding site. In a cooperative protein, the binding of one subunit to one target will lead to a conformational change that has allosteric effects on the other subunits. This will in turn affect the target binding affinity of the other sites. For example, when oxygen binds to one of the four subunits of hemoglobin, the entire complex relaxes, allowing oxygen to bind more easily to the other three subunits. With each successive binding, oxygen is taken up more readily. The binding activity of hemoglobin increases with the amount of available substrate. This leads to the “all or nothing” switching response. Cooperative systems follow the Hill equation θ = [P]n/(Kd + [P]n) where θ is the proportion of active protein and n is the Hill coefficient, which is related to, but often lower than, the number of subunits in the active complex (Figure 5B) [33].

Cooperative binding means that sigmoidal all-or-nothing response curves are likely to be more common for homo-oligomers. This could mean that the evolution of homo-oligomerization is favored for gene products whose activity needs to be robust to accidental activation, for example if gene expression is particularly noisy, or if plasticity is critical, with strong selection against accidently turning on an inappropriate pathway [35]. The de novo evolution of active homo-oligomers is likely to be rare, making causality more plausible in the opposite direction, at least under the cooperativity rather than the sequestration explanation of our results. In other words, given a homo-oligomer that reduces the consequences of any expression noise, relaxed selection allows greater noise and plasticity to evolve, either against a background of constant gene function, or in terms of which new functions are likely to evolve within a gene family. Whichever direction causality operates, we expect homo-oligomerization to be correlated with noisy and plastic gene expression, in agreement with our regression analyses. This effect became statistically detectable only after we accounted for the important confounding factor of protein stickiness, which our results suggest can be measured using Y2H data.

Conclusions

Our work validates use of the number of Y2H interactions as a metric for protein stickiness. Sticky proteins, but not proteins with more functional partners, exhibit low stochastic noise and low plasticity across environments. This presumably indicates the greater evolutionary constraints acting on intrinsically sticky proteins. Homo-oligomers also exhibit low noise and low plasticity, once their high level of stickiness is controlled for. This suggests that homo-oligomers might help mediate robustness to the consequences of noisy expression.

Methods

Protein abundance, TATA status and essentiality

Protein abundance measures were taken from Ghaemmaghami et al. [36] and subjected to a log transform. Classification of a gene’s promoter type as TATA(+) or TATA(−) was taken from Basehoar et al. [37]. The dispensability of each gene was identified using the essentiality classification of Mewes et al. [38].

Noise

Noise values for 2168 genes were taken from Newman et al. [13], who used flow cytometry to measure the fluorescence of individual cells expressing GFP-fusion proteins from their endogenous promoters. The total coefficient of variance includes substantial contributions from variation in cell size and cell cycle state. We used the gated measurements of Newman et al. [13], which minimize the effects of these confounding factors. Newman et al. [13] reported their findings both as coefficients of variance (CV) and as a distance of each CV to a running median of CVs (referred to as DM). The DM values remove the strong and intrinsic effects of protein abundance on noise, and are the most appropriate for the study of evolutionary constraints. We therefore used the DM values, taken from cells grown in rich media. We performed an optimized Box-Cox transform (λ = −1.879) to make the data normal, as assessed by a Shapiro-Wilk test for normality. Note that noise data tended to be unavailable for genes expressed at low levels.

Plasticity

mRNA expression data were downloaded from the Saccharomyces Genome Database [39,40]. We excluded 19 of the listed microarray papers on the grounds that they looked at conditions that wild yeast populations would not be expected to encounter, leaving 11 papers suitable for analysis [41-51].

Most papers included results from several trials. For example, Gasch et al. [44] contains expression data from cells using fructose, galactose, glucose, etc. as a carbon source. In contrast, Roberts et al. [41] contains expression data only from cells exposed to different concentrations and time durations of alpha factor. Our aim was to count one data point per biologically relevant environmental condition. In the first case, each microarray dataset was classified as its own experiment. In the second, since each microarray dataset involved alpha factor exposure, we grouped these as a single experiment. After we classified the microarrays within each paper, 27 independent experiments were obtained from the 11 papers.

If at least one measurement within an experiment showed a change in a gene’s expression by a factor of at least two relative to the experimental control, we counted that as an experiment for which that gene changed. For each gene, we counted the number of experiments in which a gene showed a change in expression at least as large as this arbitrary cut-off. We then performed an optimized Box-Cox transform (λ = 0.303) on this number, followed by a loess regression against protein abundance (Figure 6). Subsequent analyses were performed on the residuals from this loess regression, which are normally distributed as assessed by a Shapiro-Wilk test.

thumbnailFigure 6. Loess regression correcting plasticity for protein abundance. Statistical analyses were performed on transformed plasticity numbers (left vertical axis), untransformed plasticity is shown right for illustration. Further analysis was performed on the deviate of each data point from the red loess regression line. The R loess regression function was used rather than the lowess function because loess returns residuals and better handles larger datasets.

Protein-protein interaction data, including self-interaction

The Y2H data were isolated from downloads of the BioGRID Interaction Database, Database for Interacting Proteins (DIP), IntAct database, and Molecular INTeraction Database (MINT) [52-59], yielding 29096 unique interactions from 1680 publications. No quality filter was applied to the Y2H data. The ACMS data were taken solely from the BioGRID Interaction Database [52,53], as BioGRID provides a comprehensive listing of data for that experiment type. To help reduce the influence of false positives in the ACMS data, we only kept interactions that appeared across 2 or more independent experiments: this is similar to the ACMS filtering procedure performed by Heo et al. [24]. Note that this further biased the ACMS data towards greater accuracy for more abundant proteins. After filtering, the “True” PPI consisted of 16786 unique interactions from 436 ACMS publications. Both Y2H and ACMS data were subjected to log transforms, yielding truncated normal data distributions, as assessed visually.

Self interaction status was assessed using Y2H data. Y2H data frequently contain false positive PPIs between proteins that would never be expressed in the same place and time. This drawback is clearly not a problem for the assessment of self-interactions, and Y2H rather than ACMS data were used to minimize false negatives.

Regression models

Multiple regression models were calculated using a linear regression function (lm) in the R statistical computing environment. Continuous variables (i.e. noise, plasticity, Y2H degree, and ACMS degree) were transformed, as described above, to make the data normal or nearly so and, in the case of noise and plasticity, to control for the confounding factor of protein abundance. Note, however, that both of the log-transformed PPI variables have truncated normal distributions, due to a floor at zero. The coefficients of determination of nested models were compared using an ANOVA. See Figure 7 for a flowchart demonstrating the research procedure.

thumbnailFigure 7. Methods flowchart. Simple illustrative flowchart showing progression of research methods including datasets analysed, data transforms, statistical tests, and regression models.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JM conceived the study. LB and JM contributed to the design of the study, to the analysis and interpretation of the results, and to writing the manuscript. LB carried out the study. All authors read and approved the final manuscript.

Acknowledgements

Work was supported by the National Institutes of Health (GM072733) and by the Howard Hughes Medical Institute (52005889) via the Undergraduate Biology Research Program at the University of Arizona. J.M. is a Pew Scholar in the Biomedical Sciences. We thank Andrew Capaldi and Scott Rifkin for helpful discussions, I. Ispolatov for helpful suggestions on an earlier version of this manuscript, and Ben Wilson for help with programming and MySQL.

References

  1. Fraser HB, Hirsh AE, Giaever G, Kumm J, Eisen MB: Noise Minimization in Eukaryotic Gene Expression.

    PLoS Biol 2004, 2:e137. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Lehner B: Selection to minimise noise in living systems and its implications for the evolution of gene expression.

    Mol Syst Biol 2008.

    4.

    OpenURL

  3. Li J, Min R, Vizeacoumar FJ, Jin K, Xin X, Zhang Z: Exploiting the determinants of stochastic gene expression in Saccharomyces cerevisiae for genome-wide prediction of expression noise.

    P Natl Acad Sci USA 2010, 107:10472-10477. Publisher Full Text OpenURL

  4. Huang H, Jedynak BM, Bader JS: Where have all the interactions gone?

    Estimating the coverage of two-hybrid protein interaction maps. PLoS Comput Biol 2007, 3:e214. OpenURL

  5. Yu X, Ivanic J, Memisević V, Wallqvist A, Reifman J: Categorizing biases in high-confidence high-throughput protein-protein interaction data sets.

    Mol Cell Proteomics 2011, 10:M111.012500. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. de Silva E, Thorne T, Ingram P, Agrafioti I, Swire J, Wiuf C, Stumpf MPH: The effects of incomplete protein interaction data on structural and evolutionary inferences.

    BMC Biol 2006, 4:39. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  7. Kelly WP, Stumpf MPH: Assessing coverage of protein interaction data using capture-recapture models.

    Bull Math Biol 2012, 74:356-374. PubMed Abstract | Publisher Full Text OpenURL

  8. Gavin A-C, Aloy P, Grandi P, Krause R, Boesche M, Marzioch M, Rau C, Jensen LJ, Bastuck S, Dümpelfeld B, Edelmann A, Heurtier M-A, Hoffman V, Hoefert C, Klein K, Hudak M, Michon A-M, Schelder M, Schirle M, Remor M, Rudi T, Hooper S, Bauer A, Bouwmeester T, Casari G, Drewes G, Neubauer G, Rick JM, Kuster B, Bork P, Russell RB, Superti-Furga G: Proteome survey reveals modularity of the yeast cell machinery.

    Nature 2006, 440:631-636. PubMed Abstract | Publisher Full Text OpenURL

  9. Zhang J, Maslov S, Shakhnovich EI: Constraints imposed by non-functional protein-protein interactions on gene expression and proteome size.

    Mol Syst Biol 2008, 4:210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Deeds EJ, Ashenberg O, Shakhnovich EI: A simple physical model for scaling in protein-protein interaction networks.

    P Natl Acad Sci USA 2006, 103:311-316. Publisher Full Text OpenURL

  11. Lehner B: Conflict between Noise and Plasticity in Yeast.

    PLoS Genet 2010, 6:e1001185. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Dong D, Shao X, Deng N, Zhang Z: Gene expression variations are predictive for stochastic noise.

    Nucleic Acids Res 2010, 39:403-413. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. 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.

    Nature 2006, 441:840-846. PubMed Abstract | Publisher Full Text OpenURL

  14. Batada NN, Hurst LD: Evolution of chromosome organization driven by selection for reduced gene expression noise.

    Nat Genet 2007, 39:945-949. PubMed Abstract | Publisher Full Text OpenURL

  15. Tirosh I, Barkai N, Verstrepen KJ: Promoter architecture and the evolvability of gene expression.

    J Biol 2009, 8:95. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Choi JK, Kim Y-J: Intrinsic variability of gene expression encoded in nucleosome positioning sequences.

    Nat Genet 2009, 41:498-503. PubMed Abstract | Publisher Full Text OpenURL

  17. Landry CR, Lemos B, Rifkin SA, Dickinson WJ, Hartl DL: Genetic Properties Influencing the Evolvability of Gene Expression.

    Science 2007, 317:118-121. PubMed Abstract | Publisher Full Text OpenURL

  18. Tirosh I, Weinberger A, Carmi M, Barkai N: A genetic signature of interspecies variations in gene expression.

    Nat Genet 2006, 38:830-834. PubMed Abstract | Publisher Full Text OpenURL

  19. Tirosh I, Barkai N: Two strategies for gene regulation by promoter nucleosomes.

    Genome Res 2008, 18:1084-1091. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Zhou L, Ma X, Sun F: The effects of protein interactions, gene essentiality and regulatory regions on expression variation.

    BMC Syst Biol 2008, 2:54. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  21. Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics.

    Nat Rev Genet 2009, 10:57-63. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Ivanic J, Yu X, Wallqvist A, Reifman J: Influence of Protein Abundance on High-Throughput Protein-Protein Interaction Detection.

    PLoS One 2009, 4:e5815. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Drummond DA, Wilke CO: Mistranslation-Induced Protein Misfolding as a Dominant Constraint on Coding-Sequence Evolution.

    Cell 2008, 134:341-352. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Heo M, Maslov S, Shakhnovich E: Topology of protein interaction network shapes protein abundances and strengths of their functional and nonspecific interactions.

    P Natl Acad Sci USA 2011, 108:4258-4263. Publisher Full Text OpenURL

  25. Ispolatov I, Yuryev A, Maslov S: Binding properties and evolution of homodimers in protein-protein interaction networks.

    Nucleic Acids Res 2005, 33:3629-3635. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Deeds EJ, Ashenberg O, Gerardin J, Shakhnovich EI: Robust protein protein interactions in crowded cellular environments.

    P Natl Acad Sci USA 2007, 104:14952-14957. Publisher Full Text OpenURL

  27. Yang J-R, Liao B-Y, Zhuang S-M, Zhang J: PNAS Plus: Protein misinteraction avoidance causes highly expressed proteins to evolve slowly.

    P Natl Acad Sci USA 2012, 109:E831-E840. Publisher Full Text OpenURL

  28. Vavouri T, Semple JI, Garcia-Verdugo R, Lehner B: Intrinsic Protein Disorder and Interaction Promiscuity Are Widely Associated with Dosage Sensitivity.

    Cell 2009, 138:198-208. PubMed Abstract | Publisher Full Text OpenURL

  29. Ghim C-M, Almaas E: Genetic noise control via protein oligomerization.

    BMC Syst Biol 2008, 2:94. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  30. Bundschuh R, Hayot F, Jayaprakash C: The Role of Dimerization in Noise Reduction of Simple Genetic Networks.

    J Theor Biol 2003, 220:261-269. PubMed Abstract | Publisher Full Text OpenURL

  31. Maslov S, Ispolatov I: Propagation of large concentration changes in reversible protein-binding networks.

    P Natl Acad Sci USA 2007, 104:13655-13660. Publisher Full Text OpenURL

  32. Ferrell JE Jr: Tripping the switch fantastic: how a protein kinase cascade can convert graded inputs into switch-like outputs.

    Trends Biochem Sci 1996, 21:460-466. PubMed Abstract | Publisher Full Text OpenURL

  33. Ferrell JE Jr: Q&A: Cooperativity.

    J Biol 2009, 8:53. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. Masel J, Siegal ML: Robustness: mechanisms and consequences.

    Trends Genet 2009, 25:395-403. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Amoutzias GD, Robertson DL, Van de Peer Y, Oliver SG: Choose your partners: dimerization in eukaryotic transcription factors.

    Trends Biochem Sci 2008, 33:220-229. PubMed Abstract | Publisher Full Text OpenURL

  36. Ghaemmaghami S, Huh W-K, Bower K, Howson RW, Belle A, Dephoure N, O’Shea EK, Weissman JS: Global analysis of protein expression in yeast.

    Nature 2003, 425:737-741. PubMed Abstract | Publisher Full Text OpenURL

  37. Basehoar AD, Zanton SJ, Pugh BF: Identification and Distinct Regulation of Yeast TATA Box-Containing Genes.

    Cell 2004, 116:699-709. PubMed Abstract | Publisher Full Text OpenURL

  38. Mewes HW: MIPS: analysis and annotation of proteins from whole genomes in 2005.

    Nucleic Acids Res 2006, 34:D169-D172. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. SGD Expression Connectionhttp://downloads.yeastgenome.org/expression/microarray/archive/Expression_connection_data/ webcite.

  40. Cherry JM, Hong EL, Amundsen C, Balakrishnan R, Binkley G, Chan ET, Christie KR, Costanzo MC, Dwight SS, Engel SR, Fisk DG, Hirschman JE, Hitz BC, Karra K, Krieger CJ, Miyasato SR, Nash RS, Park J, Skrzypek MS, Simison M, Weng S, Wong ED: Saccharomyces Genome Database: the genomics resource of budding yeast.

    Nucleic Acids Res 2012, 40:D700-705. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Roberts CJ, Nelson B, Marton MJ, Stoughton R, Meyer MR, Bennett HA, He YD, Dai H, Walker WL, Hughes TR, Tyers M, Boone C, Friend SH: Signaling and circuitry of multiple MAPK pathways revealed by a matrix of global gene expression profiles.

    Science 2000, 287:873-880. PubMed Abstract | Publisher Full Text OpenURL

  42. Roberts GG, Hudson AP: Transcriptome profiling of Saccharomyces cerevisiae during a transition from fermentative to glycerol-based respiratory growth reveals extensive metabolic and structural remodeling.

    Mol Genet Genomics 2006, 276:170-186. PubMed Abstract | Publisher Full Text OpenURL

  43. Lyons TJ, Gasch AP, Gaither LA, Botstein D, Brown PO, Eide DJ: Genome-wide characterization of the Zap1p zinc-responsive regulon in yeast.

    P Natl Acad Sci USA 2000, 97:7957-7962. Publisher Full Text OpenURL

  44. Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes.

    Mol Biol Cell 2000, 11:4241-4257. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Ogawa N, DeRisi J, Brown PO: New components of a system for phosphate accumulation and polyphosphate metabolism in Saccharomyces cerevisiae revealed by genomic expression analysis.

    Mol Biol Cell 2000, 11:4309-4321. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Chu S, DeRisi J, Eisen M, Mulholland J, Botstein D, Brown PO, Herskowitz I: The transcriptional program of sporulation in budding yeast.

    Science 1998, 282:699-705. PubMed Abstract | Publisher Full Text OpenURL

  47. Yoshimoto H, Saltsman K, Gasch AP, Li HX, Ogawa N, Botstein D, Brown PO, Cyert MS: Genome-wide analysis of gene expression regulated by the calcineurin/Crz1p signaling pathway in Saccharomyces cerevisiae.

    J Biol Chem 2002, 277:31079-31088. PubMed Abstract | Publisher Full Text OpenURL

  48. Haugen AC, Kelley R, Collins JB, Tucker CJ, Deng C, Afshari CA, Brown JM, Ideker T, Van Houten B: Integrating phenotypic and expression profiles to map arsenic-response networks.

    Genome Biol 2004, 5:R95. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  49. Lai L-C, Kosorukoff AL, Burke PV, Kwast KE: Metabolic-state-dependent remodeling of the transcriptome in response to anoxia and subsequent reoxygenation in Saccharomyces cerevisiae.

    Eukaryot Cell 2006, 5:1468-1489. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale.

    Science 1997, 278:680-686. PubMed Abstract | Publisher Full Text OpenURL

  51. Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization.

    Mol Biol Cell 1998, 9:3273-3297. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Stark C, Breitkreutz B-J, Chatr-Aryamontri A, Boucher L, Oughtred R, Livstone MS, Nixon J, Van Auken K, Wang X, Shi X, Reguly T, Rust JM, Winter A, Dolinski K, Tyers M: The BioGRID Interaction Database: 2011 update.

    Nucleic Acids Res 2011, 39:D698-704. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  53. Stark C, Breitkreutz B-J, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets.

    Nucleic Acids Res 2006, 34:D535-539. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Salwinski L, Miller CS, Smith AJ, Pettit FK, Bowie JU, Eisenberg D: The Database of Interacting Proteins: 2004 update.

    Nucleic Acids Res 2004, 32:D449-451. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. Xenarios I, Rice DW, Salwinski L, Baron MK, Marcotte EM, Eisenberg D: DIP: the database of interacting proteins.

    Nucleic Acids Res 2000, 28:289-291. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Kerrien S, Aranda B, Breuza L, Bridge A, Broackes-Carter F, Chen C, Duesbury M, Dumousseau M, Feuermann M, Hinz U, Jandrasits C, Jimenez RC, Khadake J, Mahadevan U, Masson P, Pedruzzi I, Pfeiffenberger E, Porras P, Raghunath A, Roechert B, Orchard S, Hermjakob H: The IntAct molecular interaction database in 2012.

    Nucleic Acids Res 2011, 40:D841-D846. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  57. Hermjakob H: IntAct: an open source molecular interaction database.

    Nucleic Acids Res 2004, 32:452D-455. Publisher Full Text OpenURL

  58. Licata L, Briganti L, Peluso D, Perfetto L, Iannuccelli M, Galeota E, Sacco F, Palma A, Nardozza AP, Santonico E, Castagnoli L, Cesareni G: MINT, the molecular interaction database: 2012 update.

    Nucleic Acids Res 2012, 40:D857-861. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Chatr-aryamontri A, Ceol A, Palazzi LM, Nardelli G, Schneider MV, Castagnoli L, Cesareni G: MINT: the Molecular INTeraction database.

    Nucleic Acids Res 2007, 35:D572-D574. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL