- Open Access
DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines
Genome Biologyvolume 12, Article number: R10 (2011)
The Erratum to this article has been published in Genome Biology 2011 12:405
DNA methylation is an essential epigenetic mechanism involved in gene regulation and disease, but little is known about the mechanisms underlying inter-individual variation in methylation profiles. Here we measured methylation levels at 22,290 CpG dinucleotides in lymphoblastoid cell lines from 77 HapMap Yoruba individuals, for which genome-wide gene expression and genotype data were also available.
Association analyses of methylation levels with more than three million common single nucleotide polymorphisms (SNPs) identified 180 CpG-sites in 173 genes that were associated with nearby SNPs (putatively in cis, usually within 5 kb) at a false discovery rate of 10%. The most intriguing trans signal was obtained for SNP rs10876043 in the disco-interacting protein 2 homolog B gene (DIP2B, previously postulated to play a role in DNA methylation), that had a genome-wide significant association with the first principal component of patterns of methylation; however, we found only modest signal of trans-acting associations overall. As expected, we found significant negative correlations between promoter methylation and gene expression levels measured by RNA-sequencing across genes. Finally, there was a significant overlap of SNPs that were associated with both methylation and gene expression levels.
Our results demonstrate a strong genetic component to inter-individual variation in DNA methylation profiles. Furthermore, there was an enrichment of SNPs that affect both methylation and gene expression, providing evidence for shared mechanisms in a fraction of genes.
DNA methylation plays an important regulatory role in eukaryotic genomes. Alterations in methylation can affect transcription and phenotypic variation , but the source of variation in DNA methylation itself remains poorly understood. Substantial evidence of inter-individual variation in DNA methylation exists with age [2, 3], tissue [4, 5], and species . In mammals, DNA methylation is mediated by DNA methyltransferases (DNMTs) that are responsible for de novo methylation and maintenance of methylation patterns during replication. Genes involved in the synthesis of methylation and in DNA demethylation can also affect methylation variation. For example, mutations in DNMT3L  and MTHFR  associate with global DNA hypo-methylation in human blood. These changes occur at a genome-wide level and are distinct from genetic variants that impact DNA methylation variability in targeted genomic regions, for example, genetic polymorphisms associated with differential methylation in the H19/IGF2 locus .
Recent evidence suggests a dependence of DNA methylation on local sequence content [10–12]. A strong genetic effect is supported by studies of methylation patterns in families  and in twins , but stochastic and environmental factors are also likely to play an important role [2, 14]. Recent work indicates that genetic variation may have a substantial impact on local methylation patterns [5, 15–18], but neither the extent to which methylation is affected by genetic variation, nor the mechanisms are yet clear. Furthermore, the degree to which variation in DNA methylation underlies variation in gene expression across individuals remains unknown.
DNA methylation has long been considered a key regulator of gene expression. The genetic basis of gene expression has been investigated across tissues  and populations . Both lines of evidence suggest genetic variants associated with gene expression variation are located predominantly near transcription start sites. However, not much is known about the precise mechanisms by which genetic variants modify gene-expression. Combining genetic, epigenetic, and gene expression data can inform the underlying relationship between these processes, but such studies are rare on a genome-wide scale. Two recent studies have examined the link between DNA methylation and expression in human brain samples [5, 18]. Both studies identified substantial numbers of quantitative trait loci underlying each type of phenotype, but few examples of individual loci driving variation in both methylation and expression.
To better understand the role of genetic variation in controlling DNA methylation variation, and its resulting effects on gene expression variation, we studied DNA promoter methylation across the genome in 77 human lymphoblastoid cell lines (LCLs) from the HapMap collection. These cell lines represent a unique resource as they have been densely genotyped by the HapMap Project , and are now being genome-sequenced by the 1,000 Genomes Project. In addition, these cell lines have been studied by numerous groups studying variation in gene expression using microarrays [20, 22] and RNA sequencing [23, 24], as well as smaller studies of variation in chromatin accessibility and PolII binding [25, 26]. Finally, one of the HapMap cell lines is now being intensely studied by the ENCODE Project . This convergence of diverse types of genome-wide data from the same cell lines should ultimately enable a clearer understanding of the mechanisms by which genetic variation impacts gene regulation.
Characteristics of DNA promoter methylation patterns
To study inter-individual variation in methylation profiles we measured methylation levels across the genome in 77 lymphoblastoid cell lines (LCLs) derived from unrelated individuals from the HapMap Yoruba (YRI) collection. For these samples we also had publicly available genotypes , as well as estimates of gene expression levels from RNA-sequencing in 69 of the 77 samples . Methylation profiling was performed in duplicate using the Illumina HumanMethylation27 DNA Analysis BeadChip assay, which is based on genotyping of bisulfite-converted genomic DNA at individual CpG-sites to provide a quantitative measure of DNA methylation. The Illumina array includes probes that target 27,578 CpG-sites. However, we limited analyses to probes that mapped uniquely to the genome and did not contain known sequence variation, leaving us with a data set of 22,290 CpG-sites in the promoter regions of 13,236 genes (see Methods). Following hybridization, methylation levels were estimated as the ratio of intensity signal obtained from the methylated allele over the sum of methylated and unmethylated allele intensity signals. Methylation levels were quantile-normalized  across two replicates. We tested for correlations with potential confounding variables that could affect methylation levels in LCLs , such as LCL cell growth rate, copy numbers of Epstein-Barr virus, and other measures of biological variation (see Additional file 1) that were available for 60 of the individuals in our study ; these did not significantly explain variation in the methylation levels in our sample (Figure S1 in Additional file 1). However, we observed an influence of HapMap Phase (samples from Phase 1/2 vs 3) on the distribution of the first principal component loadings in the autosomal data, suggesting that the first methylation principal component may in part capture technical variation potentially related to LCL culture. In the downstream association mapping analyses, we applied a correction using principal component analysis regressing the first three principal components to account for unmeasured confounders and increase power to detect quantitative trait loci.
Global patterns of methylation
Distinct patterns of methylation were observed for CpG-sites located on the autosomes, X-chromosome, and in the vicinity of imprinted genes (Figure 1a). The majority (71.4%) of autosomal CpG-sites were primarily unmethylated (observed fraction of methylation <0.3), 15.6% were hemi-methylated (fraction of methylation was between 0.3 and 0.7), and 13% were methylated. As expected, these patterns were consistent with previously observed lower levels of methylation near promoters relative to genome-wide levels [4, 31]. We did not find evidence for sex-specific autosomal methylation patterns, consistent with a previous report . In contrast, CpG-sites on the X-chromosome exhibited highly significant sex-specific differences (Figure S2) with hemi-methylated patterns in females that were consistent with X-chromosome inactivation. A similar hemi-methylation peak was observed for CpG-sites located near the transcription start sites (TSSs) of known autosomal imprinted genes in the entire sample.
We observed a previously reported  drop in methylation levels for CpG-sites located within 1 kb of TSSs (Figure 1b). Promoter methylation levels have been reported to vary with respect to CpG islands . We found that although distance to the CpG island (CGI) border  (including CpG shores ) did not significantly affect methylation levels, CpG-sites located in CGIs were under-methylated and less variable (Wilcoxon rank-sum test P < 2.2 × 10-16) compared to sites outside of CGIs (Figure 1, Figure S3 in Additional file 1).
Methylation is often found to be correlated across genomic regions at the scale of 1-2 kb [4, 35]. We investigated whether the correlation between autosomal methylation levels (co-methylation) depended on the distance between CpG-sites. We observed that methylation levels at probes located in close proximity (up to 2 kb apart) were highly correlated (Figure 1c), indicating that variation in methylation levels between individuals is correlated within cell type. Figure 1c also shows that pairs of CpG-sites that were both within a CGI showed greater evidence for co-methylation than pairs of CpG sites for which at least one was outside the CGI, controlling for distance, implying differential regulation of DNA methylation for CpGs inside and outside of CGIs .
DNA methylation correlates with transcription and histone modifications
Methylation has long been implicated in the regulation of gene expression. To examine the role of methylation in gene expression variation, we compared methylation levels to estimates of gene expression based on RNA-sequencing (Figure 2a). Within individuals, we found a significant negative correlation between methylation and gene expression levels (Figure S4 in Additional file 1) across 11,657 genes (mean rank correlation r = -0.454). We divided the genes into quartiles from high to low gene expression and observed that the drop in methylation levels near to the TSS (Figure 1b) was only seen in highly expressed genes (Figure 2b). We also asked whether variation in methylation levels across individuals correlates with variation in gene expression levels. Comparisons at the gene level across 69 individuals indicated a modest but significant excess of negatively correlated genes (permutation P < 0.0001).
DNA methylation is thought to interact with histone modifications during the regulation of gene-expression [36, 37]. We compared methylation levels in our sample with histone modification ChIP-seq data from the ENCODE project in one of the CEPH HapMap LCLs (GM12878). We found strong negative correlations between DNA methylation levels and the presence of histone marks that target active genes (Figure 1d; Figures S3 and S5 in Additional file 1). For example, DNA methylation was low in H3K27ac peaks, which are indicative of enhancers , have previously been positively correlated with transcription levels  and negatively correlated with DNA methylation levels . Similarly, the transcription marks H3K4me3 and H3K9ac were both negatively correlated with DNA methylation levels. We also observed lower methylation levels in transcription factor binding sites predicted by the CENTIPEDE algorithm, using cell-type specific data including DNase1 sequencing reads , consistent with the expectation that the absence of methylation is important for transcription factor binding.
Genome-wide association of DNA methylation with SNP genotypes
We next assessed whether genetic variation contributes to inter-individual variation in DNA methylation levels. We first tested whether any SNPs were associated with overall patterns of DNA methylation, as measured by principal component analysis (see Methods). The most interesting signal was obtained for SNP rs10876043, which had a genome-wide significant association with variation in the first principal component of methylation (P = 4.5 × 10-9), and which also showed a modest association with average genome-wide methylation levels (P = 4.0 × 10-5) (Table S1 in Additional file 1). This SNP lies within the intron of the gene DIP2B, which contains a DMAP1-binding domain, and has been previously proposed to play a role in DNA methylation .
Associations in trans
After assessing the possibility that SNPs can have genome-wide effects on overall methylation patterns, we next transformed the methylation data by regressing out the first three principal components (see Methods), as we have previously found that this procedure can greatly reduce noise in the data and improve quantitative trait locus (QTL) mapping  (see also [42, 43]). At a genome-wide false discovery rate (FDR) of 10% (P = 2.1 × 10-10) methylation levels at 37 CpG-sites showed evidence for association with SNP genotypes (Table S2 in Additional file 1). The majority of these CpG-sites (27 of 37) were putative cis association signals, that is, the most significant SNP was within 50 kb of the measured CpG site (Figure S6 in Additional file 1). We observed a modest enrichment of distal associations (putative trans associations) that was primarily due to signals in 10 CpG-sites (Figure S7 in Additional file 1). We then examined distal association at SNPs that had previously been implicated in methylation (Table S3 in Additional file 1) and found a significant proximal association between SNP rs8075575, which is 150 kb from gene ZBTB4 that binds methylated DNA, and methylation at probe cg24181591 in gene EIF5A that encodes a translation initiation factor. Three previously reported  significant distal associations were also observed for SNP rs7225527 (38 kb from gene RHBDL3) and methylation at probe cg17704839 in gene UBL5 that encodes ubiquitin-like protein, and for SNPs rs2638971 (106 kb from gene DDX11) and rs17804971 (49 kb from gene DDX12) and methylation at probe cg18906795 in gene RANBP6, which may function in nuclear protein import as a nuclear transport receptor. Associations were also seen at SNPs located 165 kb from the gene encoding methyl-binding protein MBD2, 22 kb from the methyltransferase gene DNMT1, 192 kb from the methyltransferase gene DNMT3B, and at three SNPs with previous evidence for association but to different regions  (Figure S8 in Additional file 1). Overall however, we obtained relatively weak evidence for associations in trans and weak to moderate enrichment of trans association signals at more relaxed significance thresholds in candidate regions of interest.
Associations in cis
Since the majority of the genome-wide association signals were proximal to the corresponding CpG-sites, we next focused on association testing for SNPs within 50 kb of each CpG-site (Figure 3). At a genome-wide FDR of 10% (P = 2.0 × 10-5) there were 180 CpG-sites with cis methylation quantitative trait loci (meQTLs). The strongest association signal (P = 8.0 × 10-18) was obtained at SNP rs2187102 with probe cg27519424 in gene HLCS, which is thought to be involved in gene-regulation by mediating histone biotinylation . The proportion of variance explained by meQTLs for normalized methylation data ranged between 22% and 63%. If mechanisms affecting DNA methylation generally act over distances of up to approximately 2 kb (Figure 1c), then SNPs impacting methylation should be detected as meQTLs at multiple nearby CpG-sites. We observed that SNPs associated with methylation were also enriched for association with additional CpG-sites within 2 kb of the best-associated CpG-site with the most-significant P-value (Figure 3b), suggesting that a single genetic variant often affects methylation at numerous nearby CpG-sites.
Genetic variation has previously been associated with methylation at specific imprinted regions . The 180 CpG-sites with meQTLs in our data were nearest to the TSSs of 173 genes, of which two-MEST and CPA4, were known to be imprinted genes. Previous observations suggested that eQTL and imprinting effects can be sex-specific , raising the possibility that some of the meQTLs may act in a sex-dependent manner. However, we did not find compelling genome-wide significant sex-specific cis meQTL effects (see Additional file 1). Of the 180 associations of CpG-sites with proximal meQTLs, 27 were previously reported in human brain samples .
Little is known about the biological mechanisms that may underlie meQTL effects. To this end we applied a Bayesian hierarchical model  to test for enrichment of meQTLs in transcription factor binding sites, in histone modification categories, and in the vicinity of the associated probes. We found that SNPs located nearest to the probe, and specifically in the 5 kb immediately surrounding the probe, were significantly enriched for meQTLs (Figure 3c). Transcription factor binding sites, including CTCF-binding sites, showed a modest but non-significant enrichment for meQTLs (Figure S9 in Additional file 1).
Methylation QTLs are enriched for expression QTLs
Finally, we examined the overlap in regulatory variation that affects both methylation and gene expression levels using RNA-sequencing data . We hypothesized that since DNA methylation can regulate gene expression, then variants that affect methylation should often have consequent effects on gene expression. The first way that we looked at this was to take the set of 180 SNPs that are meQTLs at FDR <10% (taking only the most significant SNP for each meQTL). We then tested each of these SNPs for association with expression levels of nearby genes (Figure 4a, red points). There is a clear enrichment of association with expression levels compared to the null hypothesis (black line) and compared to sets of control SNPs that are matched in terms of allele frequency and distance-to-probe distributions (black dots).
One example of a SNP, rs8133082, that is both a meQTL and eQTL for the gene C21orf56 is illustrated in Figure 5. When we regress out methylation, this completely removes the association of this SNP with gene expression (Figure 5a, b, c, d). We validated the methylation assay findings at C21orf56 by bisulfite sequencing the methylation probe region in eight samples in our study, four from each homozygote genotype class for the SNP (Figure 5f). The two methylation probes at C21orf56 both had cis meQTLs and overlapped the likely promoter region as indicated by histone modification data (Figure 5e), suggesting that genetic variation may affect the chromatin structure in this region. C21orf56 appears to modulate the response of human LCLs to alkylating agents, and may act as a genomic predictor for inter-individual differences in response to DNA damaging agents .
To examine further the overlap between eQTLs and meQTLs, we re-analyzed the eQTL data by incorporating methylation as a gene-specific covariate. If variation in methylation underlies variation in gene-expression, we expect to observe a drop in the number of eQTLs in the methylation-residual gene expression data. At an FDR of 10% (P = 2.5 × 10-5) there were 484 original eQTLs and 463 methylation-residual eQTLs, where 439 eQTLs overlapped, 45 eQTLs were present only in the original data, and 24 new eQTLs were present only in the methylation-residuals (Figure 4b). Interestingly, the SNPs that were eQTLs for the 45 genes with reduced signals in the methylation-residuals were enriched for significant methylation associations (Figure S10 in Additional file 1), suggesting that these are true underlying meQTLs, where genetic variation affects methylation, which in turn regulates gene expression [5, 18]. In summary our results indicate a significant enrichment of SNPs that affect both methylation and gene expression, suggesting a shared mechanism (for example, that increased DNA methylation might drive lower gene expression). However the number of genes that show such a signal is a modest fraction of the total number of meQTLs.
We report association between DNA methylation with genetic and gene expression variation at a genome-wide level. We have identified methylation QTLs genome-wide, the majority of which act over very short distances, namely less than 5 kb. Furthermore, methylation patterns generally covary within individuals over distances of approximately 2 kb and in conjunction with this, meQTLs frequently affect multiple neighboring CpG sites. Our findings are consistent with previous methylation associations [5, 16, 18], familial aggregation [13, 14], correlation with local sequence , allele-specific methylation [15, 17], and effects of histone modifications . Little is known about the biological mechanisms that underlie meQTL effects, however, this is one important route to identify how genetic variation affects gene regulation.
We find an overall enrichment of significant associations of genetic variants with methylation CpG-sites, which is consistent with the results from two recent reports examining genome-wide methylation QTLs in human brain samples [5, 18]. Overall, the number of genome-wide significant meQTLs varies across the three studies, which is likely due to differences in sample sizes, differences in multiple testing corrections and definition of cis intervals, and the presence of large tissue-specific differences in DNA methylation with tissue-specific meQTLs. In general, power to detect meQTLs will depend on many factors including sample size, genome-wide coverage of genetic variation, genome-wide coverage of methylation variation, and the effect size of the genetic variants associated with methylation variation in the tissue of interest.
Additionally, our analyses are based on Epstein-Barr virus transformed lymphoblastoid cell lines. The choice of cell type will affect the observed genome-wide DNA methylation patterns, and in particular, high-passage LCLs may exhibit methylation alterations over time . Sun et al. , for example, investigated genome-wide differences in DNA methylation between LCLs and peripheral blood cells (PBCs), and identified 3,723 autosomal DNA methylation sites that had significantly different methylation patterns across cell types. In that respect, it is expected that a subset of our results reflect LCL-specific events. We have tested potential confounding variables that could affect methylation levels specifically in LCLs , but do not observe significant effects of these on overall DNA methylation patterns in our data. However, variation in methylation are slightly different in HapMap Phase 1/2 samples compared to HapMap Phase 3 samples, suggesting that technical variation related to LCL culture may influence DNA methylation. We took this into account when performing all downstream methylation QTL analyses, and our analyses of the uncorrected methylation patterns are consistent with the results of previous studies in primary cells [4, 31, 35].
We obtained interesting results from the trans analysis highlighting several loci with potential long-range effects on DNA methylation. Furthermore, an intriguing association of a SNP within the intron of DIP2B, which contains a DMAP1-binding domain, with the first principal component of autosomal methylation patterns suggests novel genome-wide effects on methylation variability. However, we do not observe a strong effect of polymorphisms in many of the candidate methylation regulatory genes on overall patterns of methylation or on specific probes. The sample size used in the study limits our power to detect trans signals, rendering these analyses more difficult to interpret. In general, the moderate sample sizes used in all three genome-wide methylation studies to date do not allow for the detection of subtle effects of genetic variants on methylation variation and correspondingly the majority of methylation sites assayed across all studies remains unexplained by the GWAS analyses. However, the findings indicate that genetic regulation of methylation is as complex as expression or phenotypic variation.
Relating genetic variation to both DNA methylation and gene expression variation reveals complex patterns. We observe significant overlap between meQTLs and eQTLs for cis regulatory variants. These findings were obtained when we both focus exclusively on meQTL SNPs (Figure 4a) and when we compare the genome-wide meQTL results for all SNPs classified as eQTLs in the hierarchical model framework (Figure S9 in Additional file 1). The observations indicate evidence for shared regulatory mechanisms in a fraction of genes. However, in the re-analyses of the eQTL data taking into account DNA methylation, in only 10% of eQTLs was the genetic effect of the SNP on expression affected by controlling for methylation, suggesting that variation in methylation accounts for only a small fraction of variation in gene expression levels. There may be several explanation for this. First, the coverage of the methylation array provides a relatively low resolution snapshot of the genome-wide DNA methylation patterns. Second, steady state gene expression levels (as measured by RNA-sequencing) are controlled by many other factors in addition to DNA methylation, such as transcription factor binding, chromatin state including histone marks and nucleosome positioning, and regulation by small RNAs. Finally, our study sample size provides modest power, both for eQTL and meQTL mapping. However, compared to previous studies addressing this issue [5, 18], we find more convincing evidence for meQTL and eQTL overlap. For example, Zhang et al.  found ten cases where genetic variants associated with both methylation and expression, but they only examined gene expression data for fewer than 100 genes in these comparisons in a subset of the sample, while Gibbs et al.  found that approximately 5% of SNPs in their study were significant as both meQTLs and eQTLs. Also, Gibbs et al.  find proportionally similar number of QTLs for methylation and gene expression, while we find more eQTLs. A potential explanation for the greater overlap obtained in our data is that our study examines one cell type in comparison to heterogeneous cell-types in human brain tissue samples used in both other studies [5, 18].
Characterizing the genetic control of methylation and its association to the regulation of gene expression is an important area for research, critical to our understanding of how complex living systems are regulated. Our study has the potential to help disease mapping studies, by informing the phenotypic consequences of this variation. Altogether, of the 173 genes with proximal meQTLs in our study, eighteen genes were previously reported to be differentially methylated in cancer, in other diseases, or across multiple tissues (see Table S4 in Additional file 1). Furthermore, thirty of the meQTL associations reported in our study were also observed in human brain samples . These findings provide a framework to help the interpretation of GWAS findings and improve our understanding of the underlying biology in multiple complex phenotypes.
Our results, together with recent findings of heritable allele-specific chromatin modification [25, 47] and transcription factor binding [26, 49] demonstrate a strong genetic component to inter-individual variation in epigenetic and chromatin signature, with likely downstream transcriptional and phenotypic consequences. Importantly, we found an enrichment for SNPs that affect both methylation and gene expression, implying a single causal mechanism by which one SNP may affect both processes, although such shared QTLs represent a minority of both meQTLs and eQTLs. Our data also have implications for the functional interpretation of mechanisms underlying association of genetic variants with disease.
Materials and methods
DNA was extracted from lymphoblastoid cell lines from 77 individuals from the Yoruba (YRI) population from the International HapMap project (60 HapMap Phase 1/2 and 17 HapMap Phase 3 individuals). Lymphoblastoid cell lines were previously established by Epstein-Barr Virus transformation of peripheral blood mononuclear cells using phytohemagluttinin. We obtained the transformed cell lines from the Coriell Cell Repositories. Methylation data were obtained using the Illumina HumanMethylation27 DNA Analysis BeadChip assay. Methylation estimates were assayed using two technical replicates per individual and methylation levels were quantile normalized across replicates . At each CpG-site the methylation level is presented as β, which is the fraction of signal obtained from the methylated beads over the sum of methylated and unmethylated bead signals. We considered different approaches to normalizing values across replicates, as well as using the log of the ratio of methylated to unmethylated signal instead of β, and found the results robust to normalization procedure, measure of methylation, and across technical replicates (see Additional file 1). The methylation data are publicly available  and have been submitted to the NCBI Gene Expression Omnibus  under accession no. [GSE26133].
We mapped the 27,578 Illumina probes to the human genome sequence (hg18) using BLAT  and MAQ . We selected 26,690 probes that unambiguously mapped to single locations in the human genome at a sequence identity of 100%, discarding probes that mapped to multiple locations with up to two mismatches. We excluded a further 4,400 probes that contained sequence variants, including 3,960 probes with SNPs (from the 1,000 genomes project , July 2009 release, YRI population) and 440 probes which overlapped copy number variants . This resulted in a final set of 22,290 probes (21,289 autosomal probes) that were used in all further analyses. The 22,290 probes were nearest to the TSSs of 13,236 Ensembl genes, of which 12,901 genes had at least one methylation CpG-site within 2 kb of the TSS.
Bisulfite sequencing was performed in the C21orf56 region for eight individuals. DNA was bisulfite-converted using the EZ DNA Methylation-Gold Kit (Zymo Research). PCR amplification was performed using primers designed around CpG-site cg07747299 from the HumanMethylation27 array and the nearest CpG island in the region (using Methyl Primer Express from Applied Biosystems) for a total of 411 bp amplified in the 5' UTR of the C21orf56 gene. PCR products were sequenced and cytosine peak heights compared to overall peak height were called using 4Peaks Software.
Gene expression data
RNA-sequencing data were obtained for LCLs from 69 individuals in our study from . The methylation and RNA-sequencing data were obtained from the same cultures of the LCLs. RNA-sequencing gene expression values are presented as the number of GC-corrected reads mapping to a gene in an individual, divided by the length of the gene. In the methylation to gene expression comparisons we split genes into quantiles based on the mean gene expression per gene. For the eQTL analyses, RNA-sequencing data were corrected and normalized exactly as in . Of the 22,683 genes in the original study, 10,167 autosomal genes had both gene expression counts and methylation CpG-sites within 2 kb of the TSS.
HapMap release 27 genotype data were obtained for 3.8 million autosomal SNPs in HapMap (combined Phase 1/2 and 3). Missing genotypes were imputed by BIMBAM  using the posterior mean genotype. Non-polymorphic SNPs were excluded, reducing the set to 3,035,566 autosomal SNPs for association analyses.
Spearman rank correlations were used to assess co-methylation between probes and to compare methylation and gene expression. We used 10,000 permutations of the gene expression to methylation assignments to assess the enrichment of negatively and positively correlated genes in the 25% and 5% tails within genes. Wilcoxon rank-sum tests were used to compare probe means and variances for subsets of probes.
Genome-wide association was performed using the methylation values at each CpG-site as phenotypes and three million autosomal SNP genotypes. We used least squares linear regression with a single-locus additive effects model, where we estimated the effect of the minor SNP allele on the increase in methylation levels. Prior to the association analyses, we normalized the methylation values at each CpG-site to N(0, 1) and applied a correction using principal component analysis regressing the first three principal components to account for unmeasured confounders following similar approaches to reduce expression heterogeneity in gene expression experiments [24, 42, 43] (see Additional file 1). Sex-specific analyses were performed using sex as a covariate and assessing the significance of the sex by additive-QTL interaction term.
We assessed the enrichment of association at SNPs and probes that were previously reported to be associated with methylation [7, 8, 15–18] and at SNPs within 200 kb of genes known to affect DNA methylation (Table S3 in Additional file 1). We also compared genetic variation to normalized variation in the principal components loadings for the autosomal methylation data (see Additional file 1). Results from the 180 cis meQTLs are available online .
We performed genome-wide permutations to assess the significance of the genome-wide association results in the least-square linear regressions. We permuted the methylation values for the 21,289 autosomal probes (phenotypes), performed genome-wide association on the 21,289 permuted and normalized phenotypes, and repeated this procedure for 10 (cis-analyses) or 1 (trans-analyses) replicates selecting the best signal per probe per replicate. Results are presented at an FDR of 10%, meaning that an estimated 10% of the meQTLs are false positives. Results for additional FDR thresholds are shown in Additional file 1. FDR was calculated as the fraction of significant hits in the permuted versus the observed data at a given P-value threshold. The association analyses and FDR calculations were performed for all autosomal principal components and CpG-sites in the methylation data, and for all autosomal genes in the RNA-sequencing data.
We fitted a Bayesian hierarchical model  to test whether meQTLs were over-represented in transcription factor binding sites, histone-modifications, and with respect to distance to the probe. We extended the model to fit the methylation data, where the reference point was the location of the methylation probe. Each annotation category that we examined was included in the model while accounting for distance effects.
Genome annotation data were obtained from UCSC (hg18). Histone modification data were obtained from ChIP-seq reads from the ENCODE project (Bernstein lab) for GM12878 for seven histone marks. Histone modification categories were based on estimated peaks in the read-depth distribution (see Additional file 1).
Transcription factor binding site locations were estimated using the algorithm CENTIPEDE [40, 57]. For the results presented here, CENTIPEDE started by identifying all matches in the genome to a large number of transcription factor binding motifs obtained from the TRANSFAC and JASPAR databases. It then estimated which potential binding sites are actually occupied by transcription factors in LCLs, by incorporating input data from sequence conservation, location with respect to nearby genes, and cell-specific experimental data, including DNaseI data. We used 1,136,620 non-overlapping sites from 751 transcription factor motif matches that overlapped 1,913 CpG-sites.
Centre d'Etude du Polymorphisme Humain
chromatin immunoprecipitation followed by sequencing
disco-interacting protein 2 homolog B gene
expression quantitative trait locus
false discovery rate
lymphoblastoid cell line
methylation quantitative trait locus
quantitative trait locus
single nucleotide polymorphism
transcription start site
University of California Santa Cruz genome browser
Murrell A, Heeson S, Cooper WN, Douglas E, Apostolidou S, Moore GE, Maher ER, Reik W: An association between variants in the IGF2 gene and Beckwith-Wiedemann syndrome: interaction between genotype and epigenotype. Hum Mol Genet. 2004, 13: 247–255. 10.1093/hmg/ddh013.
Rakyan VK, Down TA, Maslau S, Andrew T, Yang TP, Beyan H, Whittaker P, McCann OT, Finer S, Valdes AM, Leslie RD, Deloukas P, Spector TD: Human aging-associated DNA hypermethylation occurs preferentially at bivalent chromatin domains. Genome Res. 2010, 20: 434–439. 10.1101/gr.103101.109.
Teschendorff AE, Menon U, Gentry-Maharaj A, Ramus SJ, Weisenberger DJ, Shen H, Campan M, Noushmehr H, Bell CG, Maxwell AP, Savage DA, Mueller-Holzner E, Marth C, Kocjan G, Gayther SA, Jones A, Beck S, Wagner W, Laird PW, Jacobs IJ, Widschwendter M: Age-dependent DNA methylation of genes that are suppressed in stem cells is a hallmark of cancer. Genome Res. 2010, 20: 440–446. 10.1101/gr.103606.109.
Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, Burton J, Cox TV, Davies R, Down TA, Haefliger C, Horton R, Howe K, Jackson DK, Kunde J, Koenig C, Liddle J, Niblett D, Otto T, Pettett R, Seemann S, Thompson C, West T, Rogers J, Olek A, Berlin K, Beck S: DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006, 38: 1378–1385. 10.1038/ng1909.
Gibbs JR, van der Brug MP, Hernandez DG, Traynor BJ, Nalls MA, Lai SL, Arepalli S, Dillman A, Rafferty IP, Troncoso J, Johnson R, Zielke HR, Ferrucci L, Longo DL, Cookson MR, Singleton AB: Abundant quantitative trait Loci exist for DNA methylation and gene expression in human brain. PLoS Genet. 2010, 6: e1000952-10.1371/journal.pgen.1000952.
Enard W, Fassbender A, Model F, Adorján P, Pääbo S, Olek A: Differences in DNA methylation patterns between humans and chimpanzees. Curr Biol. 2004, 14: R148–149.
El-Maarri O, Kareta MS, Mikeska T, Becker T, Diaz-Lacava A, Junen J, Nüsgen N, Behne F, Wienker T, Waha A, Oldenburg J, Chédin F: A systematic search for DNA methyltransferase polymorphisms reveals a rare DNMT3L variant associated with subtelomeric hypomethylation. Hum Mol Genet. 2009, 18: 1755–1768. 10.1093/hmg/ddp088.
Friso S, Girelli D, Trabetti E, Olivieri O, Guarini P, Pignatti PF, Corrocher R, Choi SW: The MTHFR 1298A > C polymorphism and genomic DNA methylation in human lymphocytes. Cancer Epidemiol Biomarkers Prev. 2005, 14: 938–943. 10.1158/1055-9965.EPI-04-0601.
Heijmans BT, Kremer D, Tobi EW, Boomsma DI, Slagboom PE: Heritable rather than age-related environmental and stochastic factors dominate variation in DNA methylation of the human IGF2/H19 locus. Hum Mol Genet. 2007, 16: 547–554. 10.1093/hmg/ddm010.
Bock C, Paulsen M, Tierling S, Mikeska T, Lengauer T, Walter J: CpG island methylation in human lymphocytes is highly correlated with DNA sequence, repeats, and predicted DNA structure. PLoS Genet. 2006, 2: e26-10.1371/journal.pgen.0020026.
Bhasin M, Zhang H, Reinherz EL, Reche PA: Prediction of methylated CpGs in DNA sequences using a support vector machine. FEBS Lett. 2005, 579: 4302–4308. 10.1016/j.febslet.2005.07.002.
Handa V, Jeltsch A: Profound flanking sequence preference of Dnmt3a and Dnmt3b mammalian DNA methyltransferases shape the human epigenome. J Mol Biol. 2005, 348: 1103–1112. 10.1016/j.jmb.2005.02.044.
Bjornsson HT, Sigurdsson MI, Fallin MD, Irizarry RA, Aspelund T, Cui H, Yu W, Rongione MA, Ekström TJ, Harris TB, Launer LJ, Eiriksdottir G, Leppert MF, Sapienza C, Gudnason V, Feinberg AP: Intra-individual change over time in DNA methylation with familial clustering. JAMA. 2008, 299: 2877–2883. 10.1001/jama.299.24.2877.
Kaminsky ZA, Tang T, Wang SC, Ptak C, Oh GHT, Wong AHC, Feldcamp LA, Virtanen C, Halfvarson J, Tysk C, McRae AF, Visscher PM, Montgomery GW, Gottesman II, Martin NG, Petronis A: DNA methylation profiles in monozygotic and dizygotic twins. Nat Genet. 2009, 41: 240–245. 10.1038/ng.286.
Kerkel K, Spadola A, Yuan E, Kosek J, Jiang L, Hod E, Li K, Murty VV, Schupf N, Vilain E, Morris M, Haghighi F, Tycko B: Genomic surveys by methylation-sensitive SNP analysis identify sequence-dependent allele-specific DNA methylation. Nat Genet. 2008, 40: 904–908. 10.1038/ng.174.
Boks MP, Derks EM, Weisenberger DJ, Strengman E, Janson E, Sommer IE, Kahn RS, Ophoff RA: The relationship of DNA methylation with age, gender and genotype in twins and healthy controls. PLoS One. 2009, 4: e6767-10.1371/journal.pone.0006767.
Schalkwyk LC, Meaburn EL, Smith R, Dempster EL, Jeffries AR, Davies MN, Plomin R, Mill J: Allelic skewing of DNA methylation is widespread across the genome. Am J Hum Genet. 2010, 86: 196–212. 10.1016/j.ajhg.2010.01.014.
Zhang D, Cheng L, Badner JA, Chen C, Chen Q, Luo W, Craig DW, Redman M, Gershon ES, Liu C: Genetic control of individual differences in gene-specific methylation in human brain. Am J Hum Genet. 2010, 86: 411–419. 10.1016/j.ajhg.2010.02.005.
Dimas AS, Deutsch S, Stranger BE, Montgomery SB, Borel C, Attar-Cohen H, Ingle C, Beazley C, Gutierrez Arcelus M, Sekowska M, Gagnebin M, Nisbett J, Deloukas P, Dermitzakis ET, Antonarakis SE: Common regulatory variation impacts gene expression in a cell type-dependent manner. Science. 2009, 325: 1246–1250. 10.1126/science.1174148.
Stranger BE, Nica AC, Forrest MS, Dimas A, Bird CP, Beazley C, Ingle CE, Dunning M, Flicek P, Koller D, Montgomery S, Tavaré S, Deloukas P, Dermitzakis ET: Population genomics of human gene expression. Nat Genet. 2007, 39: 1217–1224. 10.1038/ng2142.
International HapMap Consortium, Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H, Hu W, Li C, Lin W, Liu S, Pan H, Tang X, Wang J, Wang W, Yu J, Zhang B, Zhang Q, Zhao H, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449: 851–861. 10.1038/nature06258.
Veyrieras JB, Kudaravalli S, Kim SY, Dermitzakis ET, Gilad Y, Stephens M, Pritchard JK: High-resolution mapping of expression-QTLs yields insight into human gene regulation. PLoS Genet. 2008, 4: e1000214-10.1371/journal.pgen.1000214.
Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Lach RP, Ingle C, Nisbett J, Guigo R, Dermitzakis ET: Transcriptome genetics using second generation sequencing in a Caucasian population. Nature. 2010, 464: 773–777. 10.1038/nature08903.
Pickrell JK, Marioni JC, Pai AA, Degner JF, Engelhardt BE, Nkadori E, Veyrieras JB, Stephens M, Gilad Y, Pritchard JK: Understanding mechanisms underlying human gene expression variation with RNA sequencing. Nature. 2010, 464: 768–772. 10.1038/nature08872.
McDaniell R, Lee BK, Song L, Liu Z, Boyle AP, Erdos MR, Scott LJ, Morken MA, Kucera KS, Battenhouse A, Keefe D, Collins FS, Willard HF, Lieb JD, Furey TS, Crawford GE, Iyer VR, Birney E: Heritable individual-specific and allele-specific chromatin signatures in humans. Science. 2010, 328: 235–239. 10.1126/science.1184655.
Kasowski M, Grubert F, Heffelfinger C, Hariharan M, Asabere A, Waszak SM, Habegger L, Rozowsky J, Shi M, Urban AE, Hong MY, Karczewski KJ, Huber W, Weissman SM, Gerstein MB, Korbel JO, Snyder M: Variation in transcription factor binding among humans. Science. 2010, 328: 232–235. 10.1126/science.1183621.
ENCODE Project Consortium, 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, Fiegler H, et al: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007, 447: 799–816. 10.1038/nature05874.
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. 10.1093/bioinformatics/19.2.185.
Grafodatskaya D, Choufani S, Ferreira JC, Butcher DT, Lou Y, Zhao C, Scherer SW, Weksberg R: EBV transformation and cell culturing destabilizes DNA methylation in human lymphoblastoid cell lines. Genomics. 2010, 95: 73–83. 10.1016/j.ygeno.2009.12.001.
Choy E, Yelensky R, Bonakdar S, Plenge RM, Saxena R, De Jager PL, Shaw SY, Wolfish CS, Slavik JM, Cotsapas C, Rivas M, Dermitzakis ET, Cahir-McFarland E, Kieff E, Hafler D, Daly MJ, Altshuler D: Genetic analysis of human traits in vitro: drug response and gene expression in lymphoblastoid cell lines. PLoS Genet. 2008, 4: e1000287-10.1371/journal.pgen.1000287.
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. 10.1038/nature08514.
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. 10.1038/ng1990.
Gardiner-Garden M, Frommer M: CpG islands in vertebrate genomes. J Mol Biol. 1987, 196: 261–282. 10.1016/0022-2836(87)90689-9.
Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo-and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178–186. 10.1038/ng.298.
Ball MP, Li JB, Gao Y, Lee JH, LeProust EM, Park IH, Xie B, Daley GQ, Church GM: Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat Biotechnol. 2009, 27: 361–368. 10.1038/nbt.1533.
Cedar H, Bergman Y: Linking DNA methylation and histone modification: patterns and paradigms. Nat Rev Genet. 2009, 10: 295–304. 10.1038/nrg2540.
Thomson JP, Skene PJ, Selfridge J, Clouaire T, Guy J, Webb S, Kerr ARW, 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. 10.1038/nature08924.
Heintzman ND, Hon GC, Hawkins RD, Kheradpour P, Stark A, Harp LF, Ye Z, Lee LK, Stuart RK, Ching CW, Ching KA, Antosiewicz-Bourget JE, Liu H, Zhang X, Green RD, Lobanenkov VV, Stewart R, Thomson JA, Crawford GE, Kellis M, Ren B: Histone modifications at human enhancers reflect global cell-type-specific gene expression. Nature. 2009, 459: 108–112. 10.1038/nature07829.
Kurdistani SK, Tavazoie S, Grunstein M: Mapping global histone acetylation patterns to gene expression. Cell. 2004, 117: 721–733. 10.1016/j.cell.2004.05.023.
Pique-Regi R, Degner JF, Pai AA, Gaffney DJ, Gilad Y, Pritchard JK: Accurate inference of transcription factor binding from DNA sequence and chromatin accessibility data. Genome Res. 2011,
Winnepenninckx B, Debacker K, Ramsay J, Smeets D, Smits A, FitzPatrick DR, Kooy RF: CGG-repeat expansion in the DIP2B gene is associated with the fragile site FRA12A on chromosome 12q13.1. Am J Hum Genet. 2007, 80: 221–231. 10.1086/510800.
Leek JT, Storey JD: Capturing heterogeneity in gene expression studies by surrogate variable analysis. PLoS Genet. 2007, 3: 1724–1735. 10.1371/journal.pgen.0030161.
Kang HM, Ye C, Eskin E: Accurate discovery of expression quantitative trait loci under confounding from spurious and genuine regulatory hotspots. Genetics. 2008, 180: 1909–1925. 10.1534/genetics.108.094201.
Zempleni J, Chew YC, Bao B, Pestinger V, Wijeratne SSK: Repression of transposable elements by histone biotinylation. J Nutr. 2009, 139: 2389–2392. 10.3945/jn.109.111856.
Ober C, Loisel DA, Gilad Y: Sex-specific genetic architecture of human disease. Nat Rev Genet. 2008, 9: 911–922. 10.1038/nrg2415.
Fry RC, Svensson JP, Valiathan C, Wang E, Hogan BJ, Bhattacharya S, Bugni JM, Whittaker CA, Samson LD: Genomic predictors of interindividual differences in response to DNA damaging agents. Genes Dev. 2008, 22: 2621–2626. 10.1101/gad.1688508.
Kadota M, Yang HH, Hu N, Wang C, Hu Y, Taylor PR, Buetow KH, Lee MP: Allele-specific chromatin immunoprecipitation studies show genetic influence on chromatin state in human genome. PLoS Genet. 2007, 3: e81-10.1371/journal.pgen.0030081.
Sun YV, Turner ST, Smith JA, Hammond PI, Lazarus A, Van De Rostyne JL, Cunningham JM, Kardia SLR: Comparison of the DNA methylation profiles of human peripheral blood cells and transformed B-lymphocytes. Hum Genet. 2010, 127: 651–658. 10.1007/s00439-010-0810-y.
Zheng W, Zhao H, Mancera E, Steinmetz LM, Snyder M: Genetic analysis of variation in transcription factor binding in yeast. Nature. 2010, 464: 1187–1191. 10.1038/nature08934.
Complete methylation data and results. [https://doi.org/eqtl.uchicago.edu/]
NCBI Gene Expression Omnibus. [https://doi.org/www.ncbi.nlm.nih.gov/geo/]
Kent WJ: BLAT-the BLAST-like alignment tool. Genome Res. 2002, 12: 656–664.
Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18: 1851–1858. 10.1101/gr.078212.108.
The 1000 genomes project. [https://doi.org/www.1000genomes.org/]
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, Aerts J, Andrews TD, Barnes C, Campbell P, Fitzgerald T, Hu M, Ihm CH, Kristiansson K, Macarthur DG, Macdonald JR, Onyiah I, Pang AWC, Robson S, Stirrups K, Valsesia A, Walter K, Wei J, Wellcome Trust Case Control Consortium, Tyler-Smith C, Carter NP, Lee C, Scherer SW, Hurles ME: Origins and functional impact of copy number variation in the human genome. Nature. 2010, 464: 704–712. 10.1038/nature08516.
Guan Y, Stephens M: Practical issues in imputation-based association mapping. PLoS Genet. 2008, 4: e1000279-10.1371/journal.pgen.1000279.
Servin B, Stephens M: Imputation-based analysis of association studies: candidate regions and quantitative traits. PLoS Genet. 2007, 3: e114-10.1371/journal.pgen.0030114.
Devlin AM, Singh R, Wade RE, Innis SM, Bottiglieri T, Lentz SR: Hypermethylation of Fads2 and altered hepatic fatty acid and phospholipid metabolism in mice with hyperhomocysteinemia. J Biol Chem. 2007, 282: 37082–37090. 10.1074/jbc.M704256200.
Gómez E, Caamaño JN, Bermejo-Alvarez P, Díez C, Muñoz M, Martín D, Carrocera S, Gutiérrez-Adán A: Gene expression in early expanded parthenogenetic and in vitro fertilized bovine blastocysts. J Reprod Dev. 2009, 55: 607–614.
Sandell LL, Guan XJ, Ingram R, Tilghman SM: Gatm, a creatine synthesis enzyme, is imprinted in mouse placenta. Proc Natl Acad Sci USA. 2003, 100: 4622–4627. 10.1073/pnas.0230424100.
Kim M, Patel B, Schroeder KE, Raza A, Dejong J: Organization and transcriptional output of a novel mRNA-like piRNA gene (mpiR) located on mouse chromosome 10. RNA. 2008, 14: 1005–1011. 10.1261/rna.974608.
Gius D, Cui H, Bradbury CM, Cook J, Smart DK, Zhao S, Young L, Brandenburg SA, Hu Y, Bisht KS, Ho AS, Mattson D, Sun L, Munson PJ, Chuang EY, Mitchell JB, Feinberg AP: Distinct effects on gene expression of chemical and genetic manipulation of the cancer epigenome revealed by a multimodality approach. Cancer Cell. 2004, 6: 361–371. 10.1016/j.ccr.2004.08.029.
Sun L, Huang L, Nguyen P, Bisht KS, Bar-Sela G, Ho AS, Bradbury CM, Yu W, Cui H, Lee S, Trepel JB, Feinberg AP, Gius D: DNA methyltransferase 1 and 3B activate BAG-1 expression via recruitment of CTCFL/BORIS and modulation of promoter histone methylation. Cancer Res. 2008, 68: 2726–2735. 10.1158/0008-5472.CAN-07-6654.
Morison IM, Paton CJ, Cleverley SD: The imprinted gene and parent-of-origin effect database. Nucleic Acids Res. 2001, 29: 275–276. 10.1093/nar/29.1.275.
Kong A, Steinthorsdottir V, Masson G, Thorleifsson G, Sulem P, Besenbacher S, Jonasdottir A, Sigurdsson A, Kristinsson KT, Jonasdottir A, Frigge ML, Gylfason A, Olason PI, Gudjonsson SA, Sverrisson S, Stacey SN, Sigurgeirsson B, Benediktsdottir KR, Sigurdsson H, Jonsson T, Benediktsson R, Olafsson JH, Johannsson OT, Hreidarsson AB, Sigurdsson G, DIAGRAM Consortium, Ferguson-Smith AC, Gudbjartsson DF, Thorsteinsdottir U, Stefansson K: Parental origin of sequence variants associated with complex diseases. Nature. 2009, 462: 868–874. 10.1038/nature08625.
We thank Joseph deYoung (UCLA Southern California Genotyping Consortium) for performing the Illumina methylation assays. We thank the anonymous reviewers for helpful comments. We thank Matthew Stephens, Anna di Rienzo, Barbara Engelhardt, Jean-Baptiste Veyrieras, Yongtao Guan, Kevin Bullaughey, Gorka Alkorta-Aranburu, and members of the Pritchard, Przeworski, and Stephens labs for helpful discussions. We acknowledge the ENCODE Project for providing publicly-available histone modification and DNase data (collected by the Bernstein and Crawford labs). JTB is supported by a Sir Henry Wellcome postdoctoral fellowship. RPR is supported by the Chicago Fellows Program. AAP is supported by an American Heart Association predoctoral fellowship. This work was supported by the Howard Hughes Medical Institute, and by grants from the National Institutes of Health (Genetics and Regulation Training T 532 GM007197-34 support for JFD and AAP; RO1 MH084703-01 to JKPr; and GM077959 to YG).
The authors declare that they have no competing interests.
JTB, JKPr, and YG wrote the paper and interpreted the results. JKPr and YG designed the study. JTB analyzed the data. AAP performed bisulfite sequencing and sample preparation. JKPi mapped and processed the RNA-sequencing data, and helped with the analyses. DJG mapped and processed the histone modification data. RP-R and JFD provided estimates for the transcription factor binding sites. All authors read and approved the final manuscript.
An erratum to this article is available at http://dx.doi.org/10.1186/gb-2011-12-6-405.