Email updates

Keep up to date with the latest news and content from Epigenetics & Chromatin and BioMed Central.

Open Access Highly Accessed Research

DNA methylation and transcriptional noise

Iksoo Huh1, Jia Zeng2, Taesung Park1 and Soojin V Yi2*

Author Affiliations

1 Department of Statistics, Bioinformatics and Biostatistics Laboratory, Interdisciplinary Program in Bioinformatics, Seoul National University, Seoul, 151-742, Korea

2 School of Biology, Institute of Bioengineering and Biosciences, Georgia Institute of Technology, 310 Ferst Drive, Atlanta, GA, 30332, USA

For all author emails, please log on.

Epigenetics & Chromatin 2013, 6:9  doi:10.1186/1756-8935-6-9


The electronic version of this article is the complete one and can be found online at: http://www.epigeneticsandchromatin.com/content/6/1/9


Received:1 June 2012
Accepted:5 April 2013
Published:26 April 2013

© 2013 Huh et al.; 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

DNA methylation is one of the most phylogenetically widespread epigenetic modifications of genomic DNA. In particular, DNA methylation of transcription units (‘gene bodies’) is highly conserved across diverse taxa. However, the functional role of gene body methylation is not yet fully understood. A long-standing hypothesis posits that gene body methylation reduces transcriptional noise associated with spurious transcription of genes. Despite the plausibility of this hypothesis, an explicit test of this hypothesis has not been performed until now.

Results

Using nucleotide-resolution data on genomic DNA methylation and abundant microarray data, here we investigate the relationship between DNA methylation and transcriptional noise. Transcriptional noise measured from microarrays scales down with expression abundance, confirming findings from single-cell studies. We show that gene body methylation is significantly negatively associated with transcriptional noise when examined in the context of other biological factors.

Conclusions

This finding supports the hypothesis that gene body methylation suppresses transcriptional noise. Heavy methylation of vertebrate genomes may have evolved as a global regulatory mechanism to control for transcriptional noise. In contrast, promoter methylation exhibits positive correlations with the level of transcriptional noise. We hypothesize that methylated promoters tend to undergo more frequent transcriptional bursts than those that avoid DNA methylation.

Keywords:
DNA methylation; Gene expression; Spurious transcription; Transcriptional noise

Background

DNA methylation at CpG dinucleotides is a key epigenetic modification in the human genome crucial for regulatory and developmental processes [1,2]. The degree of DNA methylation in the human genome is extensive: most CpG dinucleotides are methylated in most tissues and developmental stages examined [3-6]. In particular, transcription units, or so-called ‘gene bodies’, are even more heavily methylated than the surrounding intergenic regions [6-9].

The functional consequences of promoter methylation on chromatin configuration and transcriptional regulation are extensively documented (see, for example, [10-12]). There is also considerable evidence suggesting that DNA methylation suppresses proliferation of transposable elements (TEs) [13-15]. However, the role of gene body methylation remains largely unresolved. Recently, studies have begun to identify molecular consequences of gene body methylation. For example, gene body methylation affects pol II occupancy and histone modifications [16]. Differential levels of DNA methylation between different exons have been linked to differential inclusion and exclusion of specific exons in transcripts [17,18]. Gene body methylation may also occur as a byproduct of transcriptional processes [19]. Another possibility is that gene body methylation is simply an extension of methylation of TEs; many genes harbor TEs within their transcription units, and the main role of methylation is to suppress the proliferation of these TEs [15].

Nevertheless, the main role of gene body DNA methylation remains unresolved. In fact, it is considered as one of the most long-standing open questions regarding genomic DNA methylation [20-25]. This question is even more pertinent in light of evolutionary patterns of DNA methylation. Comparative DNA methylation studies indicate that gene body methylation is the most conserved, ancestral form of genomic DNA methylation [7,9,23,26]. Thus, elucidating the role of gene body DNA methylation may provide significant insights into the evolutionary divergence of genomic DNA methylation across taxa [9,23,26,27].

A long-standing hypothesis posits that gene body DNA methylation suppresses spurious transcription within coding regions. By doing so, gene body methylation can effectively reduce ‘transcriptional noise’ [27,28]. This hypothesis is based upon the well-accepted idea that DNA methylation is generally repressive [29]. Pervasive DNA methylation of gene bodies, and the consequent suppression of transcriptional noise, may have served as a key facilitator enabling the evolution of complex vertebrate genomes [27]. Moreover, recent studies have begun to indicate that epigenetic mechanisms are deeply implicated in regulation of gene expression variability [30-33].

However, a detailed analysis of the relationship between transcriptional noise and DNA methylation has been lacking until now, due in large part to technical difficulties. Here, capitalizing on the recent progress in genomics and epigenomics, we investigated the impact of DNA methylation on transcriptional noise, using data from the human genome. Our analyses provide, for the first time, unequivocal evidence supporting the role of gene body methylation to reducing transcriptional noise. Furthermore, we show that promoter DNA methylation is also highly significantly associated with transcriptional noise.

Results

Transcriptional noise is negatively correlated with expression abundance and associate with specific functions

Levels of gene expression vary between cells even with the same genetic materials and under the same biological conditions [34-36]. Understanding the nature and mechanism of such variability, which is commonly referred to as ‘transcriptional noise’, has manifold functional consequences [37]. Recently, there have been significant improvements in experimental methods to measure transcriptional noise, as well as in the theoretical understanding of transcriptional noise. These studies indicate that transcriptional noise may occur due to transcriptional bursting of promoters, as well as spurious transcription within coding sequences [38-41].

Transcriptional noise in multicellular organisms, such as mammals, cannot be easily dissected using experimental means. However, they can be approximated using abundant expression datasets, for example utilizing normalized variation among microarray assays between replicates of populations [42,43]. For example, Yin et al. [42] compared the transcriptional noise measured from microarrays to those measured from single-cell experiments. The two results correspond remarkably well [42]. Similar results were seen in another study, comparing expression variation among populations to experimentally measured transcriptional noise [43]. Following these approaches, in this study we approximated transcriptional noise of human genes as the coefficient of variation of transcriptional abundance, assayed between replicates of populations of the same tissue samples under normal conditions (see Methods).

There have been significant recent technical improvements in analysis of genomic DNA methylation. In particular, researchers have begun to generate whole-genome maps of DNA methylation at the nucleotide level, via whole-genome sequencing of bisulfite-converted genomic DNA [5,44,45]. This method quantifies the methylation level of each CpG dinucleotide across the whole genome, enabling us to discern gene body methylation levels for individual genes.

In this study, we analyzed DNA methylation and transcriptional noise of the prefrontal cortex (brain) and the peripheral blood mononuclear cells (blood). We chose these two tissues for the following reasons. First, we decided to analyze ‘normal’ tissues (as opposed to cell lines). While there exists vast information on transcriptional variation of cell lines, gene expression profiles of cell lines are known to have significantly diverged from those of normal tissues [46]. Consequently, we chose not to consider cell lines in the current study. Second, we chose tissues whose genome-wide methylation maps are currently available. Finally, large numbers of microarray data in the ‘control’ (as opposed to disease) conditions exist for these tissues, thereby enabling us to measure transcriptional noise with confidence. We used rigorous quality control processes to curate microarray data from these tissues (see Methods). The resulting data are from the same technical platforms, and exhibit high correlation levels among experiments (Additional files 1 and 2).

Additional file 1. Correlation between 12 brain microarray datasets used.

Format: PDF Size: 64KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 2. Correlation between blood microarray datasets used. For the interest of space, we only show 12 microarray datasets. The remaining data exhibit similarly high correspondence between datasets.

Format: PDF Size: 62KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

We examined whether the transcriptional noise calculated from these curated data exhibited similar properties to those identified from previous studies. For example, from studies of yeast, genes involved in protein synthesis exhibited lower noise compared to other genes [39,40]. At the same time, genes responding to environmental signals or stress genes showed particularly high levels of noise [39,40]. We found similar patterns in the transcriptional noise of human genes (Additional file 3). One of the most striking findings from previous studies was that transcriptional noise is approximately proportional to the expression abundance [39,40]. We observed the same scaling behavior in which the transcriptional noise was negatively associated with expression abundance in both human tissues studied (Figure 1). This observation indicates that the scaling of transcriptional noise to expression abundance is likely to be a common phenomenon across diverse taxa, underscoring common molecular mechanisms, such as random birth and death processes of mRNAs [39,40,47]. It has been also proposed that transcriptional noise is minimized for essential genes [48]. However, in our data, we did not observe enrichment of low-noise genes in essential genes (Additional file 4).

thumbnailFigure 1. Transcriptional noise and expression abundance are significantly negatively correlated in (A) brain, and (B) blood. Transcriptional noise is measured as the coefficient of variation of transcriptional abundance (see Methods section). The regression coefficients between these variables are −0.60 (P <0.001) and −0.55 (P <0.001) for brain and blood, respectively.

Additional file 3. GO enrichment analyses of genes exhibiting high or low transcriptional noise.

Format: DOC Size: 55KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 4. No enrichment of low noise genes according to gene essentiality.

Format: DOC Size: 92KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Gene body methylation and promoter methylation exhibit negative and positive associations with transcriptional noise

Our interest was in determining whether DNA methylation influences transcriptional noise. To do so, we needed to first account for the effect of expression abundance on both of these variables. This is because DNA methylation is intimately related to expression abundance [6,10,23,25], and gene expression abundance is correlated with transcriptional noise (Figure 1). In addition, other genomic variables, such as gene lengths, are also correlated with expression abundance [49,50].

Our goal was to explain the variation found in the levels of transcriptional noise using several explanatory (independent) variables. We used the following variables as explanatory variables: expression abundance, gene body methylation, promoter methylation, and gene lengths. We first examined the variance inflation factors (VIFs), which are indicators of multicolinearity among variables. None of the explanatory variables exhibited VIFs greater than 5. This demonstrated that we could assess individual contributions of each genomic trait without the influence of multicolinearity [51].

We found that, in both tissues, gene body methylation shows significant negative relations to transcriptional noise (Table 1). This is in accord with the hypothesis that gene body DNA methylation suppresses transcriptional noise [27]. As gene length increases, there may be more opportunities for spurious transcription. In other words, gene length may be positively correlated with transcriptional noise. According to our multiple linear regression analysis, however, the effect of gene length on transcriptional noise, while controlling for other factors, was negligible in the brain data, but significantly negative in the blood data (Table 1). Analyzing more tissue samples would clarify the effect of gene length on transcriptional noise. Interestingly, promoter methylation again exhibited strong positive relations with the transcriptional noise in a multiple linear regression setting (Table 1).

Table 1. Multiple linear regression models explaining variation of transcriptional noise in different tissues

In the above analyses, we analyzed gene body methylation levels after removing TEs. We also sought to include methylation of TEs specifically in our model, using the following method. We first estimated methylation levels of gene bodies that are identified as TEs according to RepeatMasker [52]. Then we included this methylation level of TEs found within each gene as a separate variable in a multiple linear regression setting. The length of TEs themselves within each gene could not be included in this model because they exhibited high VIFs (7.39 in brain, 6.62 in blood, respectively), and thus could cause multicolinearity problems. The results of this analysis, presented in Table 2, demonstrate that TE methylation is significantly negatively correlated with transcriptional noise. In other words, TE methylation may also contribute to reducing transcriptional noise. The regression coefficients of other variables are highly similar to those from Table 1, indicating that the effects of other variables are not highly influenced by the level of TE methylation.

Table 2. Multiple linear regression models explaining variation of transcriptional noise in different tissues

To attest that our results were not biased due to statistical outliers, we next performed robust regression analyses using the same explanatory variables. We used several available methods including quantile regression as well as a few well-known loss functions such as bisquare, and Hampel and Huber ([53-55]; see also Methods section). The results of these analyses (Table 3 and Additional file 5) were unanimously consistent with the previous results, indicating highly significant negative associations between the level of gene body DNA methylation and transcriptional noise, and highly significant positive associations between promoter DNA methylation and transcriptional noise. In conclusion, these analyses reveal that after controlling for other factors, gene body methylation and promoter methylation are negatively and positively correlated with transcriptional noise, respectively.

Additional file 5. Robust regression analysis using transcriptional noise as a response variable and other biological variables as explanatory variables.

Format: DOC Size: 126KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Table 3. Robust regression analyses (quantile regression for median) for the model used in Table 1

Accounting for technical noise and among individual variability of DNA methylation

One potential caveat of our approach is the presence of technical noise, or variation of gene expression caused by technical variation among experiments, on the level of gene expression variability. Our interest is in the biological variability of gene expression. As defined previously, we approximated ‘transcriptional noise’ as the coefficient of variation (CV) among the replicates of expression data, as used previously [42]. However, this measure of gene expression variability is a composite of biological noise, which is our main interest, plus technical variation among experiments. This is problematic because it is possible that technical noise might be confounded with biological noise. For example, technical variation among experiments is negatively correlated with the expression level of genes [56,57]. Thus, it is important to take into account the impact of technical noise in assessing the relationship between biological noise and DNA methylation.

To address this issue, we used a dataset on technical and biological replicates of blood gene expression. In this dataset, gene expression is measured in two sets of technical replicates across two biological experiments [58]. Using this data, we can decompose total variation of gene expression into ‘biological’ versus ‘technical’ variation. Specifically, for a specific gene using yij as the expression level of the jth technical replicate from the ith biological sample, decomposition of variance can be expressed as in Equation 1 below:

<a onClick="popup('http://www.epigeneticsandchromatin.com/content/6/1/9/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.epigeneticsandchromatin.com/content/6/1/9/mathml/M1">View MathML</a>

(1)

The left term represents the total sum of square in a gene; the first term on the right-hand side is the biological sum of squares and the second term is the technical sum of squares. Using this decomposition, we can then assess the effect of gene body methylation on the pure biological variation and on the technical variation, separately. In our first analysis, we used the biological sum of squares as the response variable, and examined the statistical effects of several predictor variables. The results of this analysis showed that gene body methylation has a significant negative effect on biological variability among samples (model 1 in Table 4). In the second analysis, we used the total sum of square as the response variable and the technical sum of square as an explanatory variable. The results from this analysis again indicated that the effect of gene body methylation on ‘biological’ transcriptional noise, after adjusting for the technical noise, is negative (model 2 in Table 4). Thus, both methods provide consistent support to our finding that gene body methylation is negatively correlated with biological variation of gene expression.

Table 4. Multiple linear regression models in which technical versus biological components of transcriptional noise are separately analyzed

Another source of variability that needs to be accounted for is variation of DNA methylation between individuals. To determine the influence of between-individual variability of DNA methylation on our results, we analyzed datasets on gene body DNA methylation from the brains of three individuals [6]. We constructed an augmented regression model, allowing the effect of gene body methylation to vary across individuals. We defined an index for each individual as an ‘individual factor’ and included it in the new model. In addition, we included interaction terms between individual factors and gene body methylation to this model. The results of these analyses (Table 5, and Additional file 6) indicate that between-individual variations of gene body methylation do not affect our findings.

Additional file 6. Multiple linear regression analyses incorporating individual factors indicate little effect of between-individual DNA methylation on transcriptional noise.

Format: DOC Size: 42KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Table 5. Regression analysis accounting for individual variation indicates little effect of between-individual variability of DNA methylation on transcriptional noise

Discussion

The human genome and other vertebrate genomes are heavily methylated in most tissues and developmental stages, a pattern referred to as ‘global’ DNA methylation [23]. This pattern is very different from what is observed in other animals and plants. In most invertebrates examined, DNA methylation is targeted to the transcription units (gene bodies) of a subset of genes [7,9,23]. Notably, gene body methylation appears to have existed well before the emergence of DNA methylation of promoters and TEs, as an ancestral form of DNA methylation in diverse animal and plant genomes [23,60,61].

Determining the role of gene body methylation is of much interest, and studies are revealing associations between gene body methylation and gene expression [9,21,62,63], transcript composition [17,18,64] and chromatin structures [16]. Nonetheless, the global role of gene body methylation remains unresolved. In this respect, two long-standing hypotheses stand out. The first hypothesis posits that gene body methylation reduces transcriptional noise [27]. Another hypothesis focuses on the impact of DNA methylation to suppress the proliferation of TEs [15]. Many TEs are found in gene bodies, thus methylation of TEs may have caused expansive methylation of gene bodies [15].

In this study we examined the predictions of these two hypotheses using whole genome methylation data and statistical methods. Because gene body methylation and transcriptional noise are both significantly correlated with expression abundance, it is important to analyze the impact of gene body methylation while considering the effect of expression abundance. We used several statistical methods to achieve this goal. We also examined the impact of noise due to technical variation among experiments, as well as between-individual variation of DNA methylation on our results. These analyses all indicate that gene body methylation, when viewed in the context of other biological factors, has a negative relationship with transcriptional noise.

Transcriptional noise is abundantly present in diverse taxa. The origin of transcriptional noise may be related to ‘transcriptional bursts’, referring to the phenomenon that transcription tends to occur in bursts [65-67]. Transcriptional noise also occurs due to transcription of non-canonical promoters within gene bodies, potentially due to the overabundance of RNA polymerase II in cellular environment [38]. Our results showing that more heavily methylated gene bodies exhibit less transcriptional noise are consistent with the idea that transcriptional noise is reduced by pervasive gene body methylation. Alternatively, the negative relationship between gene body DNA methylation and transcriptional noise may reflect an indirect association due to a third, yet unknown biological factor(s) that influence both variables.

The details of the actual underlying molecular mechanisms of such process are yet to be fully characterized. There are some well established epigenetic modifications of gene bodies are shown to directly suppress the initiation of non-canonical transcripts within coding sequences [68-70]. Emerging evidence indicate that gene body DNA methylation is likely to complement or function together with other epigenetic modifications to generate chromatin states that are repressive of the initiation and elongation of spurious transcripts. For example, DNA methylation of gene bodies reduces the efficiency of transcriptional elongation, by excluding RNA polymerase II occupancy and recruiting several repressive histone marks [16]. Gene body DNA methylation effectively excludes deposition of the histone variant H2A.Z, which tend to mark lowly expressed genes with high expression variability among tissues and biological conditions [71]. The identities of molecular components of the crosstalk between DNA methylation and histone modifications continue to be discovered (see, for example, [72]).

Interestingly, our analyses indicate that promoter DNA methylation is positively correlated with the level of transcriptional noise. The underlying molecular mechanism of this phenomenon is of great interest. One possibility is that this is related to the intrinsic susceptibility of specific promoters toward transcriptional bursting. In the simplest case, promoters appear to switch randomly between ‘ON’ and ‘OFF’ states with respect to the initiation of transcription [37,47]. Some promoters, however, remain perpetually in the ‘ON’ state (permissive to transcription) and do not exhibit bursting [73]. Such promoters exhibit less transcriptional variability compared to those undergoing switches between different transcriptional states [73]. In other words, the degree of transcriptional bursting likely varies between promoters according to their propensity toward different transcriptional states, leading to different levels of transcriptional noise among genes.

Given that there exists considerable evidence that unmethylated promoters can maintain a ‘permissive’ chromatin state [72,74], we hypothesize the following: promoters with lower level of DNA methylation are more likely to adopt and maintain a permissive transcriptional state (similar to the ‘ON’ state referred to above) and exhibit little transcriptional bursting. However, promoters that are more susceptible to DNA methylation may be more likely to undergo stochastic fluctuations between different states, facilitating transcriptional bursts, and as a consequence exhibit increased transcriptional noise. The actual molecular mechanisms underlying these processes are again likely to involve highly orchestrated interactions between DNA methylation and other epigenetic mechanisms: in particular, studies in yeast have revealed the role of nucleosome positioning in regulation of gene expression variability [31,33].

Reducing transcriptional noise is particularly important for genes that perform housekeeping functions and are therefore constantly expressed [28]. Indeed, methylation maps of distantly related animal genomes reveal that gene body methylation usually targets genes that function in ‘housekeeping’ cellular processes [26,28,75]. Thus, we hypothesize that gene body methylation functions as a primary mechanism to suppress transcriptional noise of essential housekeeping genes in diverse organisms. Gene body DNA methylation is the main mode of DNA methylation in many invertebrate species. Reducing transcriptional noise may serve as the primary function of DNA methylation in such genomes. Furthermore, the human genome is characterized by heavy transcription of non-coding regions [76,77]. Global methylation of the whole genome may have evolved as a molecular mechanism to reduce global transcriptional noise [27].

Moreover, we found that methylation of TEs within gene bodies also contributes to the suppression of transcriptional noise. Several studies now indicate that methylation of TEs may have evolved after the evolution of gene body methylation [23,61]. It will be interesting to determine whether the origin of TE methylation is related to its function to reduce intragenic transcriptional noise. Our study cannot provide a clear resolution to this question. Analyses of genomic methylation patterns of species straddling the invertebrate-vertebrate boundaries (near the origin of global DNA methylation) will be informative to determine the evolutionary sequences of these processes.

DNA methylation is known to vary among different tissues [5,6]. Given the potential role of gene body methylation in regulating transcriptional noise, it is possible that among-tissue variation of DNA methylation levels may be related to among-tissue variation of transcriptional noise. In our data, the prefrontal cortex (brain) exhibited higher methylation levels than blood (P <0.0001 by Mann–Whitney test, Figure 2). Since gene body DNA methylation is negatively correlated with transcriptional noise, we tested whether the brain exhibits lower noise compared to blood. Indeed, we found that prefrontal cortex samples (brain) exhibited significantly lower transcriptional noise compared to blood samples (Figure 2). Thus, regulation of transcriptional noise may be one mechanism determining tissue-specific or cell type-specific levels of gene body DNA methylation.

thumbnailFigure 2. Comparison of gene expression noise and DNA methylation between studied tissues (A) Comparison of mean transcriptional noise between the two tissues. The brain exhibits significantly lower transcriptional noise compared to blood (paired t test, P <0.001). (B) Methylation levels, however, are significantly higher in the brain compared to blood (paired t test, P <0.001). We only used genes for which methylation and transcriptional noise data exist in both tissues (total no. of genes = 3,644).

Conclusion

We explored the relationship between transcriptional noise and DNA methylation, using gene expression variability among different populations of cells as a proxy for transcriptional noise. Our analysis confirms the inverse relationship between gene expression abundance and transcriptional noise, while revealing novel relationships between DNA methylation and transcriptional noise. In particular gene body DNA methylation exhibits a negative correlation with transcriptional noise. This observation supports a longstanding hypothesis that gene body DNA methylation may reduce transcriptional noise. In light of evolutionary findings that gene body methylation is a widespread, conserved form of DNA methylation, the ancestral role of DNA methylation may have been related to the reduction of transcriptional noise. On the other hand, promoter DNA methylation is positively related to transcriptional noise, raising the possibility that epigenetic status of promoters may affect transcriptional bursts.

Methods

Data sources

Gene expression data was obtained from National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/ webcite) (Additional file 7). Because there are considerable technical variations between platforms, we restricted platforms to only the Affymetrix Human Genome U133 series. After quality control, we obtained a total of 52 datasets (12 datasets for the prefrontal cortex and 40 datasets for blood). Gene lengths were determined based upon the RefSeq annotation provided by the UCSC genome browser. Nucleotide resolution whole DNA methylation maps of the human prefrontal cortex (brain) were obtained from a recent study ([6], data available at NCBI Gene Omnibus under the record number GSE37202). DNA methylation maps of mature peripheral blood mononuclear cells were from Li et al. [44], generated using a similar method.

Additional file 7. List of microarray datasets used in this study.

Format: DOC Size: 85KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

DNA methylation

To obtain gene body methylation levels of non-repetitive portions of genes, we used the annotation of TEs from the RepeatMasker database (http://www.repeatmasker.org webcite). A custom Perl script was used to mask the TEs in gene bodies. For each mapped cytosine, the fractional methylation value was calculated as: total number of ‘C’ reads/(total number of ‘C’ reads + total number of ‘T’ reads), following previous studies [5,8,44]. We then calculated the fractional methylation level of each transcription unit, using the RefSeq database of hg18. Gene body methylation level for each gene was estimated as the mean fractional methylation value for all the mapped cytosines within each transcription unit. When alternative transcripts were present, we chose the longest transcript for each gene. The promoter methylation level for each gene was estimated as fractional methylation for regions spanning 1,500 bp upstream and 500 bp downstream of the transcription start site (TSS), similar to Zeng et al. [6].

Microarray data processing

Microarray raw data files were first processed using raw intensity using the MAS5.0 method [78]. Using other normalization methods provided similar results. We used the median probe intensities assigned to each gene as gene expression levels. We then analyzed correlation between pairwise samples, to assess similarities between datasets from the same tissue. Datasets within the same tissues exhibiting correlation coefficient greater than 0.8 are included in this study (Additional files 1 and 2). Quantile normalization using the R package ‘preprocesscore’ [79] was conducted within each tissue. Transcriptional noise was defined as the coefficient of variation (CV: standard deviation/mean) of transcriptional abundance within each tissue, following Yin et al. [42].

Multiple linear regression models of transcriptional noise

We performed multiple linear regression analyses to elucidate relationships between transcriptional noise and several biological factors (gene expression abundance, gene body methylation, promoter methylation, and gene lengths) simultaneously. CV and gene length were log transformed to improve normality. Our analyses indicated that the gender is not a significant variable and thus excluded from further analyses. We also examined the significance of the interaction terms between predictors. The results showed that the interaction terms were generally not significant and they were therefore removed from subsequent analyses.

Robust regression analysis was performed using various loss functions. We summarized the result of quantile regression in Table 3. We also used other well-known loss functions such as bisquare, Hampel and Huber (Additional file 6) [53-55]. All these approaches provided consistent results to those of the ordinary least squares method. Therefore, we conclude that the significance and magnitude of the explanatory variable effect is essential.

Functional enrichment analyses

Functional enrichment pattern of specific subsets of genes was assessed using the DAVID tools (http://david.abcc.ncifcrf.gov/ webcite) [80]. We used the list of genes included in our analyses as the background, and tested enrichments of specific gene ontology terms using the GO FAT annotation. We examined the mean transcriptional noise of genes in the two tissues and investigated the specific gene ontology terms for top 5% high transcriptional noise genes and 5% low transcriptional noise genes. A Benjamini multiple testing correction of the EASE score (a modified Fisher exact P value) was used to determine statistical significance of gene enrichment.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

TP and SVY designed the experiments. IH performed most of the statistical analyses. JZ and SVY generated and analyzed methylation data. TP and SVY coadvised the analyses. All authors read and approved the final manuscript.

Acknowledgements

This study was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (2012R1A3A2026438) to TP and by Georgia Tech Fund for Innovation in Research and Education (GT-FIRE) and NSF grants (MCB-0950896 and BCS-0751481) to SVY. We thank Hema Nagrajan for computational support.

References

  1. Jones PA: Functions of DNA methylation: islands, start sites, gene bodies and beyond.

    Nat Rev Genet 2012, 13:484-492. PubMed Abstract | Publisher Full Text OpenURL

  2. Jaenisch R, Bird A: Epigenetic regulation of gene expression: how the genome integrates intrinsic and environmental signals.

    Nat Genet 2003, 33:245-254. PubMed Abstract | Publisher Full Text OpenURL

  3. Gama-Sosa MA, Midgett RM, Slagel VA, Githens S, Kuo KC, Gehrke CW, Ehrlich M: Tissue-specific differences in DNA methylation in various mammals.

    Biochim Biophys Acta 1983, 740:212-219. PubMed Abstract | Publisher Full Text OpenURL

  4. Ehrlich M, Gama-Sosa MA, Huang L-H, Midgett RM, Kuo KC, McCune RA, Gehrke C: Amount and distribution of 5-methylcytosine in human DNA from different types of tissues or cells.

    Nucleic Acids Res 1982, 10:2709-2721. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, Low HM, Kin Sung KW, Rigoutsos I, Loring J, Wei C-L: Dynamic changes in the human methylome during differentiation.

    Genome Res 2010, 20:320-331. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Zeng J, Konopka G, Hunt BG, Preuss TM, Geschwind D, Yi SV: Divergent whole-genome methylation maps of human and chimpanzee brains reveal epigenetic basis of human regulatory evolution.

    Am J Hum Genet 2012, 91:455-465. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Feng S, Cokus SJ, Zhang X, Chen PY, Bostick M, Goll MG, Hetzel J, Jain J, Strauss SH, Halpern ME, Ukomadu C, Sadler KC, Pradhan S, Pellegrini M, Jacobsen SE: Conservation and divergence of methylation patterning in plants and animals.

    Proc Natl Acad Sci U S A 2010, 107:8689-8694. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM, Edsall L, Antosiewicz-Bourget J, Stewart R, Ruotti V, Millar AH, Thomson JA, Ren B, Ecker JR: Human DNA methylomes at base resolution show widespread epigenomic differences.

    Nature 2009, 462:315-322. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation.

    Science 2010, 328:916-919. PubMed Abstract | Publisher Full Text OpenURL

  10. Jones PA, Takai D: The role of DNA methylation in mammalian epigenetics.

    Science 2001, 293:1068-1070. PubMed Abstract | Publisher Full Text OpenURL

  11. Klose RJ, Bird AP: Genomic DNA methylation: the mark and its mediators.

    Trends Biochem Sci 2006, 31:89-97. PubMed Abstract | Publisher Full Text OpenURL

  12. Weber M, Hellmann I, Stadler MB, Ramos L, Pääbo S, Rebhan M, Schübeler D: Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome.

    Nat Genet 2007, 39:457-466. PubMed Abstract | Publisher Full Text OpenURL

  13. Miura A, Yonebayashi S, Watanabe K, Toyama T, Shimada H, Kakutani T: Mobilization of transposons by a mutation abolishing full DNA methylation in Arabidopsis.

    Nature 2001, 411:212-214. PubMed Abstract | Publisher Full Text OpenURL

  14. Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome.

    Nat Rev Genet 2007, 8:272-285. PubMed Abstract | Publisher Full Text OpenURL

  15. Yoder JA, Walsh CP, Bestor TH: Cytosine methylation and the ecology of intragenomic parasites.

    Trends Genet 1997, 13:335-340. PubMed Abstract | Publisher Full Text OpenURL

  16. Lorincz MC, Dickerson DR, Schmitt M, Groudine M: Intragenic DNA methylation alters chromatin structure and elongation efficiency in mammalian cells.

    Nat Struct Mol Biol 2004, 11:1068-1075. PubMed Abstract | Publisher Full Text OpenURL

  17. Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D’Souza C, Fouse SD, Johnson BE, Hong C, Nielsen C, Zhao Y, Turecki G, Delaney A, Varhol R, Thiessen N, Shchors K, Heine VM, Rowitch DH, Xing X, Fiore C, Schillebeeckx M, Jones SJ, Haussler D, Marra MA, Hirst M, Wang T, Costello JF: Conserved role of intragenic DNA methylation in regulating alternative promoters.

    Nature 2010, 466:253-257. PubMed Abstract | Publisher Full Text OpenURL

  18. Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, Oberdoerffer P, Sandberg R, Oberdoerffer S: CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing.

    Nature 2011, 479:74-79. PubMed Abstract | Publisher Full Text OpenURL

  19. Zilberman D, Gehring M, Tran RK, Ballinger T, Henikoff S: Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription.

    Nat Genet 2007, 39:61-69. PubMed Abstract | Publisher Full Text OpenURL

  20. Ball MP, Li JB, Gao Y, Lee J-H, LeProust EM, Park I-H, Xie B, Daley GQ, Church GM: Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells.

    Nat Biotech 2009, 27:361-368. Publisher Full Text OpenURL

  21. Hellman A, Chess A: Gene body-specific methylation on the active X chromosome.

    Science 2007, 315:1141-1143. PubMed Abstract | Publisher Full Text OpenURL

  22. Jones PA: The DNA methylation paradox.

    Trends Genet 1999, 15:34-37. PubMed Abstract | Publisher Full Text OpenURL

  23. Suzuki MM, Bird A: DNA methylation landscapes: provocative insights from epigenomics.

    Nat Rev Genet 2008, 9:465-476. PubMed Abstract | Publisher Full Text OpenURL

  24. Shenker N, Flanagan JM: Intragenic DNA methylation: implications of this epigenetic mechanism for cancer research.

    Br J Cancer 2012, 106:248-253. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Jjingo D, Conley AB, Yi SV, Lunyak VV, Jordan IK: On the presence and role of human gene-body DNA methylation.

    Oncotarget 2012, 3:462-474. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Sarda S, Zeng J, Hunt BG, Yi SV: The evolution of invertebrate gene body methylation.

    Mol Biol Evol 2012, 29:1907-1016. PubMed Abstract | Publisher Full Text OpenURL

  27. Bird A: Gene number, noise reduction and biological complexity.

    Trends Genet 1995, 11:94-100. PubMed Abstract | Publisher Full Text OpenURL

  28. Suzuki MM, Kerr ARW, De Sousa D, Bird A: CpG methylation is targeted to transcription units in an invertebrate genome.

    Genome Res 2007, 17:625-631. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Bird AP, Wolffe AP: Methylation-induced repression - belts, braces, and chromatin.

    Cell 1999, 99:451-454. PubMed Abstract | Publisher Full Text OpenURL

  30. Choi JK, Kim Y-J: Epigenetic regulation and the variability of gene expression.

    Nat Genet 2008, 40:141-147. PubMed Abstract | Publisher Full Text OpenURL

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

  32. Ha M, Ng DWK, Li W-H, Chen ZJ: Coordinated histone modifications are associated with gene expression variation within and between species.

    Genome Res 2011, 21:590-598. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

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

  34. Bahar R, Hartmann CH, Rodriguez KA, Denny AD, Busuttil RA, Dolle MET, Calder RB, Chisholm GB, Pollock BH, Klein CA, Vijg J: Increased cell-to-cell variation in gene expression in ageing mouse heart.

    Nature 2006, 441:1011-1014. PubMed Abstract | Publisher Full Text OpenURL

  35. Novick A, Weiner M: Enzyme induction as an all-or-none phenomenon.

    Proc Nat Acad Sci USA 1957, 43:553-566. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Ozbudak EM, Thattai M, Kurtser I, Grossman AD, van Oudenaarden A: Regulation of noise in the expression of a single gene.

    Nat Genet 2002, 31:69-73. PubMed Abstract | Publisher Full Text OpenURL

  37. Raj A, van Oudenaarden A: Nature, nurture, or chance: stochastic gene expression and its consequences.

    Cell 2008, 135:216-226. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Struhl K: Transcriptional noise and the fidelity of initiation by RNA polymerase II.

    Nat Struct Mol Biol 2007, 14:103-105. PubMed Abstract | Publisher Full Text OpenURL

  39. Bar-Even A, Paulsson J, Maheshri N, Carmi M, O’Shea E, Pilpel Y, Barkai N: Noise in protein expression scales with natural protein abundance.

    Nat Genet 2006, 38:636-643. PubMed Abstract | Publisher Full Text OpenURL

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

  41. Raser JM, O’Shea EK: Noise in gene expression: origins, consequences, and control.

    Science 2005, 309:2010-2013. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Yin S, Wang P, Deng W, Zheng H, Hu L, Hurst L, Kong X: Dosage compensation on the active X chromosome minimizes transcriptional noise of X-linked genes in mammals.

    Genome Biol 2009, 10:R74. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

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

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

  44. Li Y, Zhu J, Tian G, Li N, Li Q, Ye M, Zheng H, Yu J, Wu H, Sun J, Zhang H, Chen Q, Luo R, Chen M, He Y, Jin X, Zhang Q, Yu C, Zhou G, Sun J, Huang Y, Zheng H, Cao H, Zhou X, Guo S, Hu X, Li X, Kristiansen K, Bolund L, Xu J, Wang W, Yang H, Wang J, Li R, Beck S, Wang J, Zhang X: The DNA methylome of human peripheral blood mononuclear cells.

    PLoS Biol 2010, 8:e1000533. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Lister R, Ecker JR: Finding the fifth base: genome-wide sequencing of cytosine methylation.

    Genome Res 2009, 19:959-966. PubMed Abstract | Publisher Full Text OpenURL

  46. Lukk M, Kapushesky M, Nikkila J, Parkinson H, Goncalves A, Huber W, Ukkonen E, Brazma A: A global map of human gene expression.

    Nat Biotech 2010, 28:322-324. Publisher Full Text OpenURL

  47. Maheshri N, O’Shea EK: Living with noisy genes: how cells function reliably with inherent variability in gene expression.

    Annu Rev Biophys Biomol Struct 2007, 36:413-434. PubMed Abstract | Publisher Full Text OpenURL

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

  49. Park J, Xu K, Park T, Yi SV: What are the determinants of gene expression levels and breadths in the human genome?

    Hum Mol Genet 2012, 21:46-56. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  50. Park S, Choi S: Expression breadth and expression abundance behave differently in correlations with evolutionary rates.

    BMC Evol Biol 2010, 10:241. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  51. Montgomery DC, Peck EA, Vining GG: Introduction to Linear Regression Analysis. 4th edition. Hoboken, NJ: Wiley; 2006. OpenURL

  52. Smit AFA, Hubely R, Green P: RepeatMasker Open-3.0. [http://www.repeatmasker.org webcite]

  53. Hampel F: Contributions to the Theory of Robust Estimation. Berkeley, CA: University of California; 1968. OpenURL

  54. Huber PH: Robust estimation of a location parameter.

    Ann Math Stat 1964, 6:73-101. OpenURL

  55. Tukey JW: A survey of sampling from contaminated distributions. In Contributions to Probability and Statistics. Edited by Olkin I, Ghurye SG, Hoeffding W, Madow WG, Mann HB. Stanford, CA: Stanford University Press; 1960:448-485. OpenURL

  56. Novak JP, Sladek R, Hudson TJ: Characterization of variability in large-scale gene expression data: implications for study design.

    Genomics 2002, 79:104-113. PubMed Abstract | Publisher Full Text OpenURL

  57. Zeisel A, Amir A, Kostler WJ, Domany E: Intensity dependent estimation of noise in microarrays improves detection of differentially expressed genes.

    BMC Bioinformatics 2010, 11:400. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  58. Smiljanovic B, Grün JR, Biesen R, Schulte-Wrede U, Baumgrass R, Stuhlmüller B, Maslinski W, Hiepe F, Burmester GR, Radbruch A, Häupl T, Grützkau A: The multifaceted balance of TNF-α and type I/II interferon responses in SLE and RA: how monocytes manage the impact of cytokines.

    J Mol Med (Berl) 2012, 90:1295-1309. Publisher Full Text OpenURL

  59. Fox J, Monette G: Generalized collinearity diagnostics.

    J Am Stat Assoc 1992, 87:178-183. Publisher Full Text OpenURL

  60. Elango N, Yi SV: DNA methylation and structural and functional bimodality of vertebrate promoters.

    Mol Biol Evol 2008, 25:1602-1608. PubMed Abstract | Publisher Full Text OpenURL

  61. Simmen MW, Leitgeb S, Charlton J, Jones SJM, Harris BR, Clark VH, Bird A: Nonmethylated transposable elements and methylated genes in a chordate genome.

    Science 1999, 283:1164-1167. PubMed Abstract | Publisher Full Text OpenURL

  62. Aran D, Toperoff G, Rosenberg M, Hellman A: Replication timing-related and gene body-specific methylation of active human genes.

    Hum Mol Genet 2011, 20:670-680. PubMed Abstract | Publisher Full Text OpenURL

  63. Zeng J, Yi S: DNA methylation and genome evolution in honeybee: gene length, expression, functional enrichment co-vary with the evolutionary signature of DNA methylation.

    Genome Biol Evol 2010, 2:770-780. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Park J, Peng Z, Zeng J, Elango N, Park T, Wheeler D, Werren JH, Yi SV: Comparative analyses of DNA methylation and sequence evolution using Nasonia genomes.

    Mol Biol Evol 2011, 28:3345-3354. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  65. Golding I, Paulsson J, Zawilski SM, Cox EC: Real-time kinetics of gene activity in individual bacteria.

    Cell 2005, 123:1025-1036. PubMed Abstract | Publisher Full Text OpenURL

  66. Levsky JM, Shenoy SM, Pezo RC, Singer RH: Single-cell gene expression profiling.

    Science 2002, 297:836-840. PubMed Abstract | Publisher Full Text OpenURL

  67. Raj A, Peskin CS, Tranchina D, Vargas DY, Tyagi S: Stochastic mRNA synthesis in mammalian cells.

    PLoS Biol 2006, 4:e309. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  68. Carrozza MJ, Li B, Florens L, Suganuma T, Swanson SK, Lee KK, Shia W-J, Anderson S, Yates J, Washburn MP, Workman JL: Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription.

    Cell 2005, 123:581-592. PubMed Abstract | Publisher Full Text OpenURL

  69. Imbeault D, Gamar L, Rufiange A, Paquet E, Nourani A: The Rtt106 histone chaperone is functionally linked to transcription elongation and is involved in the regulation of spurious transcription from cryptic promoters in yeast.

    J Biol Chem 2008, 283:27350-27354. PubMed Abstract | Publisher Full Text OpenURL

  70. Cheung V, Chua G, Batada NN, Landry CR, Michnick SW, Hughes TR, Winston F: Chromatin- and transcription-related factors repress transcription from within coding regions throughout the Saccharomyces cerevisiae genome.

    PLoS Biol 2008, 6:e277. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  71. Coleman-Derr D, Zilberman D: Deposition of histone variant H2A.Z within gene bodies regulates responsive genes.

    PLoS Genet 2012, 8:e1002988. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  72. Thomson JP, Skene PJ, Selfridge J, Clouaire T, Guy J, Webb S, Kerr AR, Deaton A, Andrews R, James KD, Turner DJ, Illingworth R, Bird A: CpG islands influence chromatin structure via the CpG-binding protein Cfp1.

    Nature 2010, 464:1082-1086. PubMed Abstract | Publisher Full Text OpenURL

  73. Zenklusen D, Larson DR, Singer RH: Single-RNA counting reveals alternative modes of gene expression in yeast.

    Nat Struct Mol Biol 2008, 15:1263-1271. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  74. Illingworth RS, Bird AP: CpG islands - ‘A rough guide’.

    FEBS Lett 2009, 583:1713-1720. PubMed Abstract | Publisher Full Text OpenURL

  75. Hunt BG, Brisson JA, Yi SV, Goodisman MAD: Functional conservation of DNA methylation in the pea aphid and the honeybee.

    Genome Biol Evol 2010, 2:719-728. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  76. Birney E, Stamatoyannopoulos JA, Dutta A, Guigó R, Gingeras TR, Margulies EH, Weng Z, Snyder M, Dermitzakis ET, Thurman RE, Kuehn MS, Taylor CM, Neph S, Koch CM, Asthana S, Malhotra A, Adzhubei I, Greenbaum JA, Andrews RM, Flicek P, Boyle PJ, Cao H, Carter NP, Clelland GK, Davis S, Day N, Dhami P, Dillon SC, Dorschner MO, ENCODE Project Consortium: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project.

    Nature 2007, 447:799-816. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  77. Ponting CP, Belgard TG: Transcribed dark matter: meaning or myth?

    Hum Mol Genet 2010, 19:R162-R168. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Hubbell E, Liu W-M, Mei R: Robust estimators for expression analysis.

    Bioinformatics 2002, 18:1585-1592. PubMed Abstract | Publisher Full Text OpenURL

  79. Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias.

    Bioinformatics 2003, 19:185-193. PubMed Abstract | Publisher Full Text OpenURL

  80. Dennis G, Sherman B, Hosack D, Yang J, Gao W, Lane H, Lempicki R: DAVID: Database for Annotation, Visualization, and Integrated Discovery.

    Genome Biol 2003, 4:R60. BioMed Central Full Text OpenURL