BACKGROUND
gene expression changes in the cell allow for reprogramming of cellular behavior depending on the extracellular conditions.
global gene expression profiling of cells has become a routine procedure in biology, and extensive work has been done in the recent years studying gene expression programs under various conditions  <cit> .
in addition, many aspects of gene expression behavior at the dna and chromatin level have also been identified  <cit> .
although these studies yielded much insight into the regulation of gene expression under the specific conditions studied, we do not have a clear understanding of global patterns in gene expression regulation in human cells in response to extracellular stimuli.
some notable studies addressing functional aspects of gene expression regulation at a systems level have been performed in yeast  <cit> , however, an analysis of general trends in the gene expression response of human cells to extracellular cues and of their functional consequences on the regulation of human cell behavior has not been performed.
we undertook a functional analysis of global trends in the expression variance of human genes in response to extracellular cues.
expression variance of a gene can be defined as the frequency and magnitude of change in its mrna levels in response to changing extracellular conditions and can be thought of as regulatability of a gene at the mrna level.
first, we report that human genes display a wide spectrum of ev under physiological conditions, with some genes showing very little variation in their mrna levels, while some have extremely variable expression across a wide range of conditions.
the ev pattern of genes strongly correlates with their promoter architecture, such that genes with lowest ev have open promoters with constitutive rna polymerase occupancy, while those with highest ev have closed promoters with little or no rna polymerase occupancy.
then, we show that this pattern of ev under physiological conditions reflects positioning of genes in the hierarchy of cell signaling, such that the most highly regulated genes are located at the apical parts of signaling hierarchy and are generally functionally more specialized.
finally, we discuss implications of these findings on our understanding of the generic mechanisms of regulation of cell behavior as it relates to restructuring of the intracellular protein interactome.
this study uncovers some of the basic principles of transcriptional response in human cells and expands our understanding of conditional gene expression at the protein network level suggested by earlier studies on tissue-specific gene expression  <cit> .
RESULTS
calculating expression variance of human genes
in order to calculate global patterns of expression variance  of human genes, we used the extensive collection of human cancer tissue microarrays of the expression project for oncology  of the international genomics consortium , which contains expression microarray profiles of  <dig> tumor tissue samples.
this dataset contains samples dissected from diverse tissues with various types of cancer with different characteristics and treatments, and therefore spans a wide spectrum of cellular environments.
we calculated evs for each human gene by taking statistical variance of normalized expression levels of each gene across the whole dataset.
normalization of expression levels of each gene was done by first normalizing the samples in the dataset by quantile normalization  <cit> , so that each sample in the dataset has an identical distribution, and then dividing the expression level of the gene in each sample by the median of its expression level across all samples, which resulted in a measure of deviation of the expression levels of each gene from their median.
representative plots of expression profiles of some low and high ev genes are shown in figure  <dig>  varying levels of expression variation is evident between high and low ev genes.
expression levels of genes with low ev seem to be relatively stable regardless of extracellular conditions, while some genes seem to have extremely variable expression pattern across many different conditions.
this suggests that some genes may be preferentially regulated during cellular adaptation to its environment, while some genes are generally not regulated at a transcriptional level.
it is possible that the ev values simply reflect basal tissue specific expression variations of genes and not the variability of their expression under different cellular conditions.
in order to test this, we calculated tissue-specific evs of genes using only samples in the expo dataset collected from breast, lung or colon, thereby obtaining evs of genes for each individual tissue type.
if the evs reflect tissue-specific expression variations of genes, there should not be high correlation between tissue-specific evs.
however, there is a high correlation between breast and lung tissue-specific evs .
similarly high correlation was also observed between ovary and lung tissue evs , indicating that the ev values identified here represent gene-specific variations of mrna levels.
ev reflects gene regulation under varying extracellular conditions
to further confirm that our ev values reflect cellular response to varying extracellular conditions rather than being an artifact of tissue samples, we compiled an independent collection of microarray gene expression datasets from  <dig> different studies measuring responses of cultured human cells to various receptor ligands  .
we normalized each microarray sample in the dataset by their corresponding controls  and discarded the control samples, so that each sample in the ck dataset reflects fold-changes in response to the corresponding stimulus.
therefore, this dataset contains measurements of gene expression change in various human cell lines under  <dig> different stimulation conditions.
the expression variance of genes calculated using the ck dataset is in a significantly high agreement with the ev values calculated using the expo dataset of  <dig> tissue samples .
since evck values reflect fold changes of gene expression upon a large number of different stimulations, our observation indicates that evexpo values reflect true expression variations of genes within cells under different extracellular environments, rather than an artifact of tissue- or cell type-specific expressional variations.
ev is not an artifact of mrna abundance
total mrna expression levels of genes are extremely variable , and this can substantially contribute to the variability of genes between different conditions.
indeed, ev of genes has a significant negative correlation with their average expression levels across the whole expo dataset , so that genes with low mrna abundance are more likely to have variable expression.
therefore, it is possible that our observations above and below simply reflect the correlations of total expression levels of gene mrnas rather than their variability.
in order to test this, we calculated partial correlation between evexpo and evck having controlled for average mrna expression levels of genes and find that the correlation strength between these ev values calculated from different datasets is still significantly strong , indicating that the observed ev is not an artifact of mrna abundance.
in order to confirm this observation, we selected genes with similar average mrna levels , and tested if the correlation between evexpo and evck is still high.
indeed, although the correlation of total mrna levels with either evck or evexpo is lost , the correlation between evck and evexpo is still significantly high , which strongly suggests that expression variability is an intrinsic characteristic of genes rather than an artifact of their total mrna abundance.
importantly, the correlations we present below are also reproducible having controlled for the total mrna levels of genes .
ev reflects rna polymerase ii promoter occupancy
next, we asked if evs of genes correlates with the pattern of their promoter activities.
kim et al  conducted a comprehensive study mapping active promoters across the human genome and identified  <dig> classes of genes based on their expression and rna polymerase ii pre-initiation complex  promoter occupancy  <cit> .
the first class  of genes had pic occupancy and increased histone acetylation in their promoters and were actively transcribed.
the second class  had pic occupancy, however were not actively transcribed.
the third class of genes were actively transcribed although no pic could be detected, while the fourth class  had no pic occupancy or detectable transcript levels and had reduced histone acetylation at their promoters  <cit> .
we find that there is a high concordance between ev values and these gene classes , which was also reproduced with evs of the ck dataset .
low-ev genes mostly belong to classes i and iii, reflecting constitutive and high promoter activity, whereas high-ev genes mostly belong to classes ii and iv, which are mostly expressed in a condition-specific manner.
this is concordant with the invariant constitutive expression pattern of low ev genes and highly variant condition-specific expression pattern of high ev genes.
these observations suggest that low ev genes are generally highly active and abundantly transcribed, while high ev genes are transcribed in a condition-specific manner.
functional distinction of genes based on ev
we have previously shown organization of genes into separate modules based on their expression variation in yeast  <cit> .
we wanted to determine if expression variation in human genes has a functional significance similar to that observed in yeast.
in order to answer this question, we constructed a comprehensive network of human genes based on their functional similarity, where each interaction is between two genes sharing a significant functional annotation from either gene ontology  <cit>  or kegg  <cit>  .
then, we tested whether subnetworks of genes with specific functional associations segregated based on their ev.
in order to gain a comprehensive view of gene-gene association preferences in the fun-net based on their ev, we binned genes into  <dig> bins based on their ev and calculated interaction preferences between each bin pair in the fun-net.
as expected, the heatmap of interaction preferences shows a clear clustering of low and high ev genes into separate functional categories .
this is not an artifact of the network connectivity, as this pattern is not observed in a network where node positions have been randomly shuffled .
similarly analysis using different bin sizes did not significantly alter the outcome .
a similar high correlation and interaction preference pattern, albeit weaker, is observed when protein-protein interaction network is used for gene-gene interactions instead of the fun-net .
these observations show that human genes can be functionally separated based on their ev patterns.
low overall association of genes with low and high ev genes in either network suggests that the cellular functions performed by the low and high ev genes are distinct, similar to what we have shown for yeast.
next, in order to see which cellular functions are represented by the high and low ev genes, we calculated relative enrichment of the top and bottom  <dig> genes within the ev distribution for specific go functional categories.
figure 2c shows pie-charts of most enriched  functional categories in  <dig> genes with lowest or highest ev.
genes encoding for cellular functions pertaining to cellular homeostasis: mrna transcription and processing, protein synthesis and proteasomal protein degradation, are the most significantly enriched functional categories among genes with lowest ev .
however, genes exhibiting highest ev are mainly composed of genes encoding proteins in the extracellular space, including extracellular matrix  components, growth factors and extracellular proteases.
a similar pattern is identified using the evck values , where the values reflect fold inductions of genes within the same cell line in response to a treatment.
therefore, the differential enrichment of high and low ev genes for, respectively, extracellular space and intracellular homeostasis genes reflects biological pattern of cellular response to extracellular conditions.
ev of signaling genes reflects their role in the signaling hierarchy
based on the observations above, we reasoned that extracellular ligands for cellular transmembrane receptors may be more variable than their receptors, meaning that cells are more likely to modulate the expression levels of secreted factors rather than their receptors in response to extracellular cues.
we compared evs of genes annotated as "receptor binding" , "growth factor activity"  or "cytokine activity"  and "extracellular space"   to those annotated as "transmembrane receptor activity" , "receptor activity"  and "plasma membrane"   .
although evs of both classes are significantly higher when compared to the rest of genes, evs of the sf list are significantly higher than those of the gr list .
we wanted to determine if evs of genes involved in signal transduction reflect the hierarchical position of the corresponding signaling molecules within the signaling network.
in order to answer this question, we compiled a comperehensive signaling network from online databases , and defined  <dig> levels of signaling hierarchy based on the positions of the signaling molecules .
the first level, growth factor modulators , are secreted molecules that modulate the activities of receptor-binding secreted factors.
this class includes genes such as sfrp <dig> , mmp <dig> , igfbp <dig>  and ltbp <dig> .
the next two levels are secreted factors  and receptors , explained above.
receptor substrates  are molecules immediately downstream of receptors , such as g-proteins , receptor-associated kinases , and adaptor proteins  among others; and the next class are molecules that mediate signal transduction downstream of rs  .
strikingly, ev patterns of these levels display a gradient, with the gm level being the most variable, and rs being the least variable among these hierarchy levels .
this pattern is also reproduced with evck values .
this suggests that transcriptional regulation of intracellular signaling pathways mostly happens at the level of secreted growth factor modulators and growth factors, while signaling molecules immediately downstream of signaling receptors seem to be the least transcriptionally modulated.
accordingly, the rs and rs <dig> levels are mostly found in class i through iii of genes based on their pic occupancy, while gm, sf and gr levels are in class iv .
interestingly, class ii, which represents genes with pic occupancy but no detectable transcription, is enriched for intracellular signaling genes in the rs level, not secreted factors, although this class has significantly high ev .
this may indicate that class ii contains condition-specific intracellular signaling genes, while classes i and iii are enriched for constitutively expressed intracellular signaling genes.
indeed, genes with class ii promoters contain high ev genes of the rs level, while classes i and iii contain the low ev genes of the rs level .
importantly, these observations suggest not only that genes coding for extra- and intra-cellular proteins can be distinguished based on their promoter architecture, but also that promoters of intracellular proteins among themselves are distinguished based on whether they are constitutive or condition-specific.
in addition, while genes for condition-specific extracellular proteins are located within densely packed hypo-acetylated regions, condition-specific intracellular genes have relatively open promoters with a pre-assembled pic.
this suggests that regulation of transcription of genes coding for extracellular proteins may be fundamentally different from those coding for intracellular signaling genes .
ev reflects functional centrality
next, we asked whether the correlation of ev patterns of signaling genes with their positioning within the signaling hierarchy has a biological significance in terms of regulation of signaling pathways.
secreted factors are very specific in terms of their signaling targets  therefore variation of their levels through transcriptional regulation may provide a high level of specificity in the regulation of signaling pathways.
however, the receptor substrates are utilized by many signaling pathways , and therefore variation in their expression can lead to major rearrangements in the signaling architecture of the cell.
therefore, it is possible that transcriptional regulation of signaling in response to extracellular cues involves highly selective activation/inhibition of specific pathways, rather than involving large rearrangements of the signaling network.
in order to test this, we compared total number of protein-protein interactions within each hierarchical class.
since the functions of signaling molecules mainly involve protein-protein interactions, the total number of protein-protein interactions of a signaling protein may provide an estimate of the number of different processes/functions that it can be involved in.
indeed, the rs class has the most overall number of interactions among all the classes .
in order to see if lower ev of most central proteins is a general trend in the intracellular protein interaction network, we correlated total number of interactions of the rs and rs <dig> proteins with their evs.
there is a significant negative correlation of evs of rs and rs <dig> proteins and their number of protein interactions , which supports our hypothesis that differential expressional variation of genes within the signaling hierarchy can at least in part be explained by their functional centrality in the signaling network.
moreover, the rs and rs <dig> level proteins with high ev  have significantly less number of protein-protein interactions than those with low ev  .
in addition, the intracellular signaling genes with class ii promoters, which are mainly condition-specific , have significantly less number of protein-protein interactions than those with class i promoters , which are mainly constitutive .
in order to get a view of the layout of proteins with different evs within the signaling network, we plotted the signaling network of proteins of the rs and rs <dig> levels with high and low evs.
this network in figure  <dig> shows a clear differential distribution of high and low ev proteins.
while proteins with low ev are mainly located in the central dense regions of the network, those with high ev are mainly located at the periphery and generally have sparser connections.
these observations support our hypothesis that condition-specific genes encode proteins with less central roles in the signaling network.
the finding that genes with lowest evs mostly comprise genes involved in the cellular homeostatic processes , which can be regarded as the most central cellular processes, also adds to the hypothesis that mrna levels of genes with highest functional centrality in the cell are modulated less and the most variable genes are those encoding more specialized regulatory proteins.
discussion &
CONCLUSIONS
expressional variation of human genes
computational studies in yeast combining large-scale gene expression data with protein interaction networks have revealed high level of modularity in the network with respect to transcriptional regulation  <cit> .
however, with the exception of some recent studies  <cit> , such studies with human data have not been performed.
here, we report a study of global patterns in the expressional variation of human genes across a wide spectrum of conditions, and the functional significance of ev with respect to the regulation of signaling network architecture.
our findings were reproduced using two independent data compendiums, suggesting that these observations reflect true biological relationships.
in addition, since variations in mrna levels of genes have been shown to be in a relatively high agreement with corresponding variations in protein levels  <cit> , the patterns of ev discovered in this study give insight into the patterns of regulation of signaling networks in response to extracellular stimuli.
our results show that human genes are extremely variable in the extent of regulation of their mrna levels.
while some genes' mrna levels are highly variable across many conditions, some show very tight expression patterns with very little variation.
as expected, genes with lowest ev are those involved in cellular "housekeeping" functions, such as mrna synthesis and processing as well as protein synthesis and degradation.
in agreement with prior data about condition-specific genes  <cit> , genes with high ev mainly have "covered" promoters with reduced histone acetylation and no rna polymerase pre-initiation complex  occupancy, while genes with low ev have high pic occupancy and increased histone acetylation in their promoters.
transcriptional regulation of intracellular and extracellular proteins
our analyses correlating previous classification of genes into  <dig> distinct classes of promoters by kim et al   <cit>  revealed that there is a high concordance of ev values with their promoter architectures.
low ev genes are abundantly and actively transcribed, while high ev genes are generally not active.
most interestingly, high ev genes coding for intracellular signaling proteins have acetylated promoters with pre-assembled pic, while high ev genes coding for extracellular proteins have hypo-acetylated promoters without pic.
importantly, this may imply that the regulation of gene expression for extracellular proteins involves chromatin remodeling and pic assembly, while that for intracellular proteins occurs at the level of rna polymerase ii elongation, rather than pic assembly and chromatin remodeling.
it has been reported that promoter-proximal pausing of the rna polymerase ii and its subsequent release for elongation is a major mechanism of regulation of human gene expression  <cit> , which suggests that this mechanism may be employed for the regulation of intracellular proteins.
it should be noted that the study of kim et al was performed on human fibroblasts, and therefore it could be argued that the classification of genes into distinct promoter classes may be specific to fibroblasts, despite the observed high correlation of the ev patterns with these classes.
we find that ev of genes is highly similar between different tissues  and different conditions , suggesting that a common pool of condition-specific genes may exist, selective modulation of which may drive cell adaptation.
similarly, genes with lowest ev are primarily those with housekeeping functions in the cell, and are therefore likely to be expressed in all cell types.
therefore, it is likely that the overall chromatin architecture of most human promoters is also largely conserved between different cell types, and tissue and cell type-specific promoters may constitute a relative minority.
this hypothesis is not far-fetched, as another recent study analyzing chromatin architecture around gene promoters in a number of different human cell lines reported more than 70% similarity in observed positioning of nucleosomes in promoters of different cell types  <cit> .
regulation of signaling network architecture
the observation that genes regulated most in response to extracellular stimuli are secreted factors and their receptors implies that regulation of cell behavior mostly involves modulation of the composition of the extracellular environment.
even the intracellular signaling proteins with high ev seem to be mainly those with specialized roles in the regulation of signaling and with fewer number of functional interactions.
this indicates that the repertoire of the extracellular space and of their receptors mostly determines cell behavior, while the intracellular signaling hubs are mainly common for different cell types/conditions.
since in a scale-free network, such as the protein-protein interaction network  <cit> , highly connected hubs play an important role in determining the overall architecture  <cit> , our findings may suggest that the overall architecture of the signaling network is relatively stable across different conditions.
therefore, regulation of cell signaling during cell adaptation is mainly at the level of signaling inputs at the extracellular space, and minor highly specific rearrangements within the intracellular network.
this in turn suggests that relatively same signaling network architecture allows for integration of various inputs to elicit a variety of cell fates, reminiscent of a multifunctional electronic circuit.
a relatively stable network architecture where the hubs are involved in multiple processes may be evolutionarily more advantageous over a highly dynamic network architecture where hubs are condition-specific.
it is interesting to note that similar conclusions have been drawn from recent studies on tissue-specific genes, where it was reported that tissue-specific proteins are enriched for extracellular proteins  <cit> , and another reporting that tissue-specific proteins generally have less number of protein interactions  <cit> .
therefore, it is possible that the regulatory principles in response to diverse external stimuli uncovered in this study also apply to tissue-specific modulation of cell behavior.
it can be argued that quantitating protein-protein interactions to show relative centrality of proteins may introduce artifacts of historically more studied proteins.
however, we suggest it is a fair assumption that the distribution of well-studied proteins across the ev spectrum is relatively uniform so as to allow for the detection of statistically significant patterns.

 gene expression changes in the cell allow for reprogramming of cellular behavior depending on the extracellular conditions. global gene expression profiling of cells has become a routine procedure in biology, and extensive work has been done in the recent years studying gene expression programs under various conditions. in addition, many aspects of gene expression behavior at the dna and chromatin level have also been identified. although these studies yielded much insight into the regulation of gene expression under the specific conditions studied, we do not have a clear understanding of global patterns in gene expression regulation in human cells in response to extracellular stimuli. some notable studies addressing functional aspects of gene expression regulation at a systems level have been performed in yeast, however, an analysis of general trends in the gene expression response of human cells to extracellular cues and of their functional consequences on the regulation of human cell behavior has not been performed. we undertook a functional analysis of global trends in the expression variance of human genes in response to extracellular cues. expression variance of a gene can be defined as the frequency and magnitude of change in its mrna levels in response to changing extracellular conditions and can be thought of as regulatability of a gene at the mrna level. first, we report that human genes display a wide spectrum of ev under physiological conditions, with some genes showing very little variation in their mrna levels, while some have extremely variable expression across a wide range of conditions. the ev pattern of genes strongly correlates with their promoter architecture, such that genes with lowest ev have open promoters with constitutive rna polymerase occupancy, while those with highest ev have closed promoters with little or no rna polymerase occupancy. then, we show that this pattern of ev under physiological conditions reflects positioning of genes in the hierarchy of cell signaling, such that the most highly regulated genes are located at the apical parts of signaling hierarchy and are generally functionally more specialized. finally, we discuss implications of these findings on our understanding of the generic mechanisms of regulation of cell behavior as it relates to restructuring of the intracellular protein interactome. this study uncovers some of the basic principles of transcriptional response in human cells and expands our understanding of conditional gene expression at the protein network level suggested by earlier studies on tissue-specific gene expression.  calculating expression variance of human genes in order to calculate global patterns of expression variance of human genes, we used the extensive collection of human cancer tissue microarrays of the expression project for oncology of the international genomics consortium, which contains expression microarray profiles of tumor tissue samples. this dataset contains samples dissected from diverse tissues with various types of cancer with different characteristics and treatments, and therefore spans a wide spectrum of cellular environments. we calculated evs for each human gene by taking statistical variance of normalized expression levels of each gene across the whole dataset. normalization of expression levels of each gene was done by first normalizing the samples in the dataset by quantile normalization, so that each sample in the dataset has an identical distribution, and then dividing the expression level of the gene in each sample by the median of its expression level across all samples, which resulted in a measure of deviation of the expression levels of each gene from their median. representative plots of expression profiles of some low and high ev genes are shown in figure varying levels of expression variation is evident between high and low ev genes. expression levels of genes with low ev seem to be relatively stable regardless of extracellular conditions, while some genes seem to have extremely variable expression pattern across many different conditions. this suggests that some genes may be preferentially regulated during cellular adaptation to its environment, while some genes are generally not regulated at a transcriptional level. it is possible that the ev values simply reflect basal tissue specific expression variations of genes and not the variability of their expression under different cellular conditions. in order to test this, we calculated tissue-specific evs of genes using only samples in the expo dataset collected from breast, lung or colon, thereby obtaining evs of genes for each individual tissue type. if the evs reflect tissue-specific expression variations of genes, there should not be high correlation between tissue-specific evs. however, there is a high correlation between breast and lung tissue-specific evs. similarly high correlation was also observed between ovary and lung tissue evs, indicating that the ev values identified here represent gene-specific variations of mrna levels. ev reflects gene regulation under varying extracellular conditions to further confirm that our ev values reflect cellular response to varying extracellular conditions rather than being an artifact of tissue samples, we compiled an independent collection of microarray gene expression datasets from different studies measuring responses of cultured human cells to various receptor ligands. we normalized each microarray sample in the dataset by their corresponding controls and discarded the control samples, so that each sample in the ck dataset reflects fold-changes in response to the corresponding stimulus. therefore, this dataset contains measurements of gene expression change in various human cell lines under different stimulation conditions. the expression variance of genes calculated using the ck dataset is in a significantly high agreement with the ev values calculated using the expo dataset of tissue samples. since evck values reflect fold changes of gene expression upon a large number of different stimulations, our observation indicates that evexpo values reflect true expression variations of genes within cells under different extracellular environments, rather than an artifact of tissue- or cell type-specific expressional variations. ev is not an artifact of mrna abundance total mrna expression levels of genes are extremely variable, and this can substantially contribute to the variability of genes between different conditions. indeed, ev of genes has a significant negative correlation with their average expression levels across the whole expo dataset, so that genes with low mrna abundance are more likely to have variable expression. therefore, it is possible that our observations above and below simply reflect the correlations of total expression levels of gene mrnas rather than their variability. in order to test this, we calculated partial correlation between evexpo and evck having controlled for average mrna expression levels of genes and find that the correlation strength between these ev values calculated from different datasets is still significantly strong, indicating that the observed ev is not an artifact of mrna abundance. in order to confirm this observation, we selected genes with similar average mrna levels, and tested if the correlation between evexpo and evck is still high. indeed, although the correlation of total mrna levels with either evck or evexpo is lost, the correlation between evck and evexpo is still significantly high, which strongly suggests that expression variability is an intrinsic characteristic of genes rather than an artifact of their total mrna abundance. importantly, the correlations we present below are also reproducible having controlled for the total mrna levels of genes. ev reflects rna polymerase ii promoter occupancy next, we asked if evs of genes correlates with the pattern of their promoter activities. kim et al conducted a comprehensive study mapping active promoters across the human genome and identified classes of genes based on their expression and rna polymerase ii pre-initiation complex promoter occupancy. the first class of genes had pic occupancy and increased histone acetylation in their promoters and were actively transcribed. the second class had pic occupancy, however were not actively transcribed. the third class of genes were actively transcribed although no pic could be detected, while the fourth class had no pic occupancy or detectable transcript levels and had reduced histone acetylation at their promoters. we find that there is a high concordance between ev values and these gene classes, which was also reproduced with evs of the ck dataset. low-ev genes mostly belong to classes i and iii, reflecting constitutive and high promoter activity, whereas high-ev genes mostly belong to classes ii and iv, which are mostly expressed in a condition-specific manner. this is concordant with the invariant constitutive expression pattern of low ev genes and highly variant condition-specific expression pattern of high ev genes. these observations suggest that low ev genes are generally highly active and abundantly transcribed, while high ev genes are transcribed in a condition-specific manner. functional distinction of genes based on ev we have previously shown organization of genes into separate modules based on their expression variation in yeast. we wanted to determine if expression variation in human genes has a functional significance similar to that observed in yeast. in order to answer this question, we constructed a comprehensive network of human genes based on their functional similarity, where each interaction is between two genes sharing a significant functional annotation from either gene ontology or kegg. then, we tested whether subnetworks of genes with specific functional associations segregated based on their ev. in order to gain a comprehensive view of gene-gene association preferences in the fun-net based on their ev, we binned genes into bins based on their ev and calculated interaction preferences between each bin pair in the fun-net. as expected, the heatmap of interaction preferences shows a clear clustering of low and high ev genes into separate functional categories. this is not an artifact of the network connectivity, as this pattern is not observed in a network where node positions have been randomly shuffled. similarly analysis using different bin sizes did not significantly alter the outcome. a similar high correlation and interaction preference pattern, albeit weaker, is observed when protein-protein interaction network is used for gene-gene interactions instead of the fun-net. these observations show that human genes can be functionally separated based on their ev patterns. low overall association of genes with low and high ev genes in either network suggests that the cellular functions performed by the low and high ev genes are distinct, similar to what we have shown for yeast. next, in order to see which cellular functions are represented by the high and low ev genes, we calculated relative enrichment of the top and bottom genes within the ev distribution for specific go functional categories. figure 2c shows pie-charts of most enriched functional categories in genes with lowest or highest ev. genes encoding for cellular functions pertaining to cellular homeostasis: mrna transcription and processing, protein synthesis and proteasomal protein degradation, are the most significantly enriched functional categories among genes with lowest ev. however, genes exhibiting highest ev are mainly composed of genes encoding proteins in the extracellular space, including extracellular matrix components, growth factors and extracellular proteases. a similar pattern is identified using the evck values, where the values reflect fold inductions of genes within the same cell line in response to a treatment. therefore, the differential enrichment of high and low ev genes for, respectively, extracellular space and intracellular homeostasis genes reflects biological pattern of cellular response to extracellular conditions. ev of signaling genes reflects their role in the signaling hierarchy based on the observations above, we reasoned that extracellular ligands for cellular transmembrane receptors may be more variable than their receptors, meaning that cells are more likely to modulate the expression levels of secreted factors rather than their receptors in response to extracellular cues. we compared evs of genes annotated as "receptor binding", "growth factor activity" or "cytokine activity" and "extracellular space" to those annotated as "transmembrane receptor activity", "receptor activity" and "plasma membrane". although evs of both classes are significantly higher when compared to the rest of genes, evs of the sf list are significantly higher than those of the gr list. we wanted to determine if evs of genes involved in signal transduction reflect the hierarchical position of the corresponding signaling molecules within the signaling network. in order to answer this question, we compiled a comperehensive signaling network from online databases, and defined levels of signaling hierarchy based on the positions of the signaling molecules. the first level, growth factor modulators, are secreted molecules that modulate the activities of receptor-binding secreted factors. this class includes genes such as sfrp, mmp, igfbp and ltbp. the next two levels are secreted factors and receptors, explained above. receptor substrates are molecules immediately downstream of receptors, such as g-proteins, receptor-associated kinases, and adaptor proteins among others; and the next class are molecules that mediate signal transduction downstream of rs. strikingly, ev patterns of these levels display a gradient, with the gm level being the most variable, and rs being the least variable among these hierarchy levels. this pattern is also reproduced with evck values. this suggests that transcriptional regulation of intracellular signaling pathways mostly happens at the level of secreted growth factor modulators and growth factors, while signaling molecules immediately downstream of signaling receptors seem to be the least transcriptionally modulated. accordingly, the rs and rs levels are mostly found in class i through iii of genes based on their pic occupancy, while gm, sf and gr levels are in class iv. interestingly, class ii, which represents genes with pic occupancy but no detectable transcription, is enriched for intracellular signaling genes in the rs level, not secreted factors, although this class has significantly high ev. this may indicate that class ii contains condition-specific intracellular signaling genes, while classes i and iii are enriched for constitutively expressed intracellular signaling genes. indeed, genes with class ii promoters contain high ev genes of the rs level, while classes i and iii contain the low ev genes of the rs level. importantly, these observations suggest not only that genes coding for extra- and intra-cellular proteins can be distinguished based on their promoter architecture, but also that promoters of intracellular proteins among themselves are distinguished based on whether they are constitutive or condition-specific. in addition, while genes for condition-specific extracellular proteins are located within densely packed hypo-acetylated regions, condition-specific intracellular genes have relatively open promoters with a pre-assembled pic. this suggests that regulation of transcription of genes coding for extracellular proteins may be fundamentally different from those coding for intracellular signaling genes. ev reflects functional centrality next, we asked whether the correlation of ev patterns of signaling genes with their positioning within the signaling hierarchy has a biological significance in terms of regulation of signaling pathways. secreted factors are very specific in terms of their signaling targets therefore variation of their levels through transcriptional regulation may provide a high level of specificity in the regulation of signaling pathways. however, the receptor substrates are utilized by many signaling pathways, and therefore variation in their expression can lead to major rearrangements in the signaling architecture of the cell. therefore, it is possible that transcriptional regulation of signaling in response to extracellular cues involves highly selective activation/inhibition of specific pathways, rather than involving large rearrangements of the signaling network. in order to test this, we compared total number of protein-protein interactions within each hierarchical class. since the functions of signaling molecules mainly involve protein-protein interactions, the total number of protein-protein interactions of a signaling protein may provide an estimate of the number of different processes/functions that it can be involved in. indeed, the rs class has the most overall number of interactions among all the classes. in order to see if lower ev of most central proteins is a general trend in the intracellular protein interaction network, we correlated total number of interactions of the rs and rs proteins with their evs. there is a significant negative correlation of evs of rs and rs proteins and their number of protein interactions, which supports our hypothesis that differential expressional variation of genes within the signaling hierarchy can at least in part be explained by their functional centrality in the signaling network. moreover, the rs and rs level proteins with high ev have significantly less number of protein-protein interactions than those with low ev. in addition, the intracellular signaling genes with class ii promoters, which are mainly condition-specific, have significantly less number of protein-protein interactions than those with class i promoters, which are mainly constitutive. in order to get a view of the layout of proteins with different evs within the signaling network, we plotted the signaling network of proteins of the rs and rs levels with high and low evs. this network in figure shows a clear differential distribution of high and low ev proteins. while proteins with low ev are mainly located in the central dense regions of the network, those with high ev are mainly located at the periphery and generally have sparser connections. these observations support our hypothesis that condition-specific genes encode proteins with less central roles in the signaling network. the finding that genes with lowest evs mostly comprise genes involved in the cellular homeostatic processes, which can be regarded as the most central cellular processes, also adds to the hypothesis that mrna levels of genes with highest functional centrality in the cell are modulated less and the most variable genes are those encoding more specialized regulatory proteins. discussion &  expressional variation of human genes computational studies in yeast combining large-scale gene expression data with protein interaction networks have revealed high level of modularity in the network with respect to transcriptional regulation. however, with the exception of some recent studies, such studies with human data have not been performed. here, we report a study of global patterns in the expressional variation of human genes across a wide spectrum of conditions, and the functional significance of ev with respect to the regulation of signaling network architecture. our findings were reproduced using two independent data compendiums, suggesting that these observations reflect true biological relationships. in addition, since variations in mrna levels of genes have been shown to be in a relatively high agreement with corresponding variations in protein levels, the patterns of ev discovered in this study give insight into the patterns of regulation of signaling networks in response to extracellular stimuli. our results show that human genes are extremely variable in the extent of regulation of their mrna levels. while some genes' mrna levels are highly variable across many conditions, some show very tight expression patterns with very little variation. as expected, genes with lowest ev are those involved in cellular "housekeeping" functions, such as mrna synthesis and processing as well as protein synthesis and degradation. in agreement with prior data about condition-specific genes, genes with high ev mainly have "covered" promoters with reduced histone acetylation and no rna polymerase pre-initiation complex occupancy, while genes with low ev have high pic occupancy and increased histone acetylation in their promoters. transcriptional regulation of intracellular and extracellular proteins our analyses correlating previous classification of genes into distinct classes of promoters by kim et al revealed that there is a high concordance of ev values with their promoter architectures. low ev genes are abundantly and actively transcribed, while high ev genes are generally not active. most interestingly, high ev genes coding for intracellular signaling proteins have acetylated promoters with pre-assembled pic, while high ev genes coding for extracellular proteins have hypo-acetylated promoters without pic. importantly, this may imply that the regulation of gene expression for extracellular proteins involves chromatin remodeling and pic assembly, while that for intracellular proteins occurs at the level of rna polymerase ii elongation, rather than pic assembly and chromatin remodeling. it has been reported that promoter-proximal pausing of the rna polymerase ii and its subsequent release for elongation is a major mechanism of regulation of human gene expression, which suggests that this mechanism may be employed for the regulation of intracellular proteins. it should be noted that the study of kim et al was performed on human fibroblasts, and therefore it could be argued that the classification of genes into distinct promoter classes may be specific to fibroblasts, despite the observed high correlation of the ev patterns with these classes. we find that ev of genes is highly similar between different tissues and different conditions, suggesting that a common pool of condition-specific genes may exist, selective modulation of which may drive cell adaptation. similarly, genes with lowest ev are primarily those with housekeeping functions in the cell, and are therefore likely to be expressed in all cell types. therefore, it is likely that the overall chromatin architecture of most human promoters is also largely conserved between different cell types, and tissue and cell type-specific promoters may constitute a relative minority. this hypothesis is not far-fetched, as another recent study analyzing chromatin architecture around gene promoters in a number of different human cell lines reported more than 70% similarity in observed positioning of nucleosomes in promoters of different cell types. regulation of signaling network architecture the observation that genes regulated most in response to extracellular stimuli are secreted factors and their receptors implies that regulation of cell behavior mostly involves modulation of the composition of the extracellular environment. even the intracellular signaling proteins with high ev seem to be mainly those with specialized roles in the regulation of signaling and with fewer number of functional interactions. this indicates that the repertoire of the extracellular space and of their receptors mostly determines cell behavior, while the intracellular signaling hubs are mainly common for different cell types/conditions. since in a scale-free network, such as the protein-protein interaction network, highly connected hubs play an important role in determining the overall architecture, our findings may suggest that the overall architecture of the signaling network is relatively stable across different conditions. therefore, regulation of cell signaling during cell adaptation is mainly at the level of signaling inputs at the extracellular space, and minor highly specific rearrangements within the intracellular network. this in turn suggests that relatively same signaling network architecture allows for integration of various inputs to elicit a variety of cell fates, reminiscent of a multifunctional electronic circuit. a relatively stable network architecture where the hubs are involved in multiple processes may be evolutionarily more advantageous over a highly dynamic network architecture where hubs are condition-specific. it is interesting to note that similar conclusions have been drawn from recent studies on tissue-specific genes, where it was reported that tissue-specific proteins are enriched for extracellular proteins, and another reporting that tissue-specific proteins generally have less number of protein interactions. therefore, it is possible that the regulatory principles in response to diverse external stimuli uncovered in this study also apply to tissue-specific modulation of cell behavior. it can be argued that quantitating protein-protein interactions to show relative centrality of proteins may introduce artifacts of historically more studied proteins. however, we suggest it is a fair assumption that the distribution of well-studied proteins across the ev spectrum is relatively uniform so as to allow for the detection of statistically significant patterns. 
we report that human genes are highly diverse in terms of their ev; while some genes have highly variable expression pattern, some seem to be relatively ubiquitously expressed across a wide range of conditions. although much has been learnt about gene expression programs under specific conditions, global patterns of expressional variation of genes and their relationship to cellular functions and physiological responses is poorly understood. our analysis shows that this pattern of ev reflects functional centrality: proteins with highly specific signaling functions are modulated more frequently than those with highly central functions in the network, which is also consistent with previous studies on tissue-specific gene expression. to understand global principles of relationship between transcriptional regulation of human genes and their functions, we have leveraged large-scale datasets of human gene expression measurements across a wide spectrum of cell conditions. our analyses suggest a generic systems level mechanism of regulation of the cellular signaling network at the transcriptional level. the wide spectrum of gene ev strongly correlates with the positioning of proteins within the signaling network hierarchy, such that, secreted extracellular receptor ligands and membrane receptors have the highest ev, and intracellular signaling proteins have the lowest ev in the genome. interestingly, these patterns of ev along the signaling network hierarchy have significant correlations with promoter architectures of respective genes.
 understanding organizational principles of cellular networks is one of the central goals of systems biology. although much has been learnt about gene expression programs under specific conditions, global patterns of expressional variation of genes and their relationship to cellular functions and physiological responses is poorly understood.  to understand global principles of relationship between transcriptional regulation of human genes and their functions, we have leveraged large-scale datasets of human gene expression measurements across a wide spectrum of cell conditions. we report that human genes are highly diverse in terms of their ev; while some genes have highly variable expression pattern, some seem to be relatively ubiquitously expressed across a wide range of conditions. the wide spectrum of gene ev strongly correlates with the positioning of proteins within the signaling network hierarchy, such that, secreted extracellular receptor ligands and membrane receptors have the highest ev, and intracellular signaling proteins have the lowest ev in the genome. our analysis shows that this pattern of ev reflects functional centrality: proteins with highly specific signaling functions are modulated more frequently than those with highly central functions in the network, which is also consistent with previous studies on tissue-specific gene expression. interestingly, these patterns of ev along the signaling network hierarchy have significant correlations with promoter architectures of respective genes.  our analyses suggest a generic systems level mechanism of regulation of the cellular signaling network at the transcriptional level. 