A Genome-Wide Study of DNA Methylation Patterns and Gene Expression Levels in Multiple Human and Chimpanzee Tissues
The modification of DNA by methylation is an important epigenetic mechanism that affects the spatial and temporal regulation of gene expression. Methylation patterns have been described in many contexts within and across a range of species. However, the extent to which changes in methylation might underlie inter-species differences in gene regulation, in particular between humans and other primates, has not yet been studied. To this end, we studied DNA methylation patterns in livers, hearts, and kidneys from multiple humans and chimpanzees, using tissue samples for which genome-wide gene expression data were also available. Using the multi-species gene expression and methylation data for 7,723 genes, we were able to study the role of promoter DNA methylation in the evolution of gene regulation across tissues and species. We found that inter-tissue methylation patterns are often conserved between humans and chimpanzees. However, we also found a large number of gene expression differences between species that might be explained, at least in part, by corresponding differences in methylation levels. In particular, we estimate that, in the tissues we studied, inter-species differences in promoter methylation might underlie as much as 12%–18% of differences in gene expression levels between humans and chimpanzees.
Published in the journal:
. PLoS Genet 7(2): e32767. doi:10.1371/journal.pgen.1001316
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1001316
Summary
The modification of DNA by methylation is an important epigenetic mechanism that affects the spatial and temporal regulation of gene expression. Methylation patterns have been described in many contexts within and across a range of species. However, the extent to which changes in methylation might underlie inter-species differences in gene regulation, in particular between humans and other primates, has not yet been studied. To this end, we studied DNA methylation patterns in livers, hearts, and kidneys from multiple humans and chimpanzees, using tissue samples for which genome-wide gene expression data were also available. Using the multi-species gene expression and methylation data for 7,723 genes, we were able to study the role of promoter DNA methylation in the evolution of gene regulation across tissues and species. We found that inter-tissue methylation patterns are often conserved between humans and chimpanzees. However, we also found a large number of gene expression differences between species that might be explained, at least in part, by corresponding differences in methylation levels. In particular, we estimate that, in the tissues we studied, inter-species differences in promoter methylation might underlie as much as 12%–18% of differences in gene expression levels between humans and chimpanzees.
Introduction
Changes in the regulation of gene expression levels have long been hypothesized to play an important role in primate evolution [1], [2]. To begin to address this hypothesis, a large number of studies have characterized gene expression differences across primates, in particular between humans and chimpanzees [3]–[9]. These studies have pointed to several classes of biological processes (such as transcriptional regulation, oxidative stress response, and a number of metabolic pathways), which might have evolved under natural selection in primates. In addition, in a few cases, comparative studies in primates have been able to draw strong connections between regulatory adaptations and ultimate physiological or anatomical phenotypes [10]–[15].
Despite the wealth of comparative gene expression data, there are many fewer studies of the mechanisms that underlie inter-primate differences in gene regulation (e.g., [12], [13], [16]–[18]). In particular, we know relatively little about the degree to which changes in epigenetic profiles might explain differences in gene expression levels between primates.
One of the most extensively studied epigenetic mechanisms is DNA methylation – an epigenetic modification that facilitates fine-tuned regulation of transcription rates [19], [20]. Spatial and temporal regulation of transcription by DNA methylation has been shown to play an important role in many contexts, including in female X-chromosome inactivation [21], [22], genomic imprinting [23], [24], and susceptibility to complex diseases in humans, especially cancers [25], [26]. Methylation is also essential for proper differentiation and development of mammalian tissues [27], [28]. For instance, the knockout of genes encoding for the DNA-methyl-transferase (DNMT) enzymes, which are responsible for de-novo methylation of DNA, results in embryonic lethality in mice [29], [30].
The causal relationship between changes in promoter DNA methylation and differences in gene regulation has been well established [28], [31]. It has been shown that hyper-methylation at promoter CpG islands typically results in decreased transcription of downstream genes [32]. When methylation is experimentally removed from promoter regions, transcription levels rise [33]. The specific mechanisms by which DNA methylation affects gene regulation are less clear, though DNA methylation is thought to interact with proteins (such as methyl-DNA binding proteins) that associate with histone modifications or the nucleosome in order to maintain a silenced chromatin state [28], [31], [34], [35]. Additionally, it has been proposed that the binding of the transcriptional machinery and enhancer-related transcription factors to methylated genomic regions is less frequent, resulting in decreased transcription levels or absolute gene silencing [28], [36].
Previous studies have typically described patterns of DNA methylation in a single or few tissues across species [26], [37]–[41] or in multiple tissues or developmental stages within a single organism [26], [27], [34], [42]–[45]. Comparative studies of DNA methylation across mammals have suggested that the role of DNA methylation in tissue-specific gene regulation is generally conserved. For example, after identifying Tissue-specific Differentially Methylated Regions (T-DMRs [42]), in heart, colon, kidney, testis, spleen, and muscle tissues in mice, Kitamura and colleagues were able to use the methylation status in orthologous human regions to distinguish between the corresponding human tissues [44]. Irizarry and colleagues [26], who studied genome-wide DNA methylation patterns in spleen, liver, and brain tissues from human and mouse, reported that 51% of T-DMRs are shared across both species. However, there also are a large number of potentially functional differences in methylation levels across species. In particular, in primates, Gama-Sosa and colleagues [39] found that relative methylation levels within tissues generally differ between species, with the exception of hyper-methylation in the brain and thymus, which were observed regardless of species. In addition, Enard and colleagues [38], who compared methylation profiles of 36 genes in livers, brains, and lymphocytes from humans and chimpanzees, reported significant inter-species methylation level differences in 22 of the 36 genes, in at least one tissue.
With few exceptions, however (e.g., [46]), comparative studies in primates have not explored the extent to which methylation differences between species might contribute to the genome-wide regulation of inter-species differences in gene expression levels. Towards this goal, we compared genome-wide gene expression levels and DNA methylation data in tissue samples from humans and chimpanzees.
Results
We characterized DNA promoter methylation across the genome in samples from heart, liver, and kidney tissues from both humans and chimpanzees, using two technical replicates from six individual samples of each tissue from each species (see Figure S1 for an illustration of the study design). Since genome-wide gene expression data were previously collected from the same tissue samples [8], we were able to study the relationships between DNA methylation and gene expression levels across tissues and species. The gene-specific expression level estimates and methylation profiles, for all samples, are provided in Table S1.
DNA methylation varies more across tissues than between humans and chimpanzees
We obtained methylation profiles from each sample (using two independent DNA extraction replicates) by using the Illumina HumanMethylation27 DNA Analysis BeadChip assay, which provides reproducible (Figure S2) quantitative estimates of methylation levels at 27,578 CpG-loci near transcription start sites. Since the 50 bp probes on the Illumina array were designed to interrogate human samples, we limited our analysis to probes that were a perfect sequence match to the chimpanzee genome. In addition, we only used probes that were associated with genes for which we had expression measurements across the three tissues [8]. Following these exclusion criteria, we retained 10,575 CpG site probes in the putative promoter regions of 7,723 genes (see Methods for more details). At each probe, DNA methylation levels were estimated using the Illumina-recommended β values, which are essentially estimates of the proportion of methylated DNA at each CpG site (see Methods).
We note that limiting our analysis to identical methylation probes in humans and chimpanzees resulted in a slight (0.5%) but significant decrease of the median sequence divergence estimates within 500 bp windows around the retained probes (Figure S3). As a result, it is possible that, in what follows, we slightly underestimate the proportion of inter-species differences in methylation levels. However, we confirmed that limiting our analysis to identical methylation probes in the two species did not result in a noticeable shift in the distribution of expression levels of the associated genes, nor in the proportion of observed differences in gene expression levels between the two species.
As a first step of our analysis, we examined patterns of promoter methylation across tissues and species. As expected [28], [31], we found a negative correlation between methylation and gene expression levels in each individual, whereby, regardless of tissue and species, the promoters of highly expressed genes tended to be lowly methylated while the promoters of lowly expressed genes were usually highly methylated (Figure 1A; Figure S4). We also confirmed that methylation patterns on the X-chromosome account for variation due to sex, regardless of species, as expected due to X-inactivation in mammalian females [21] (the first component of variance, corresponding to sex, accounts for 67% of the overall variation in the X-chromosome data; Figure 1B). Finally, we found that genes known to be imprinted in humans tend to show a similar hemi-methylation pattern in chimpanzees (permutation tests P<0.001; Figure 1C), suggesting that the imprinted status of this set of genes is conserved in the two species.
For the remainder of the analyses, we considered only the methylation data from autosomal probes. We observed that methylation patterns across different tissues and species were quite distinct (Figure 2; similar patterns for the expression data in Figure S5). The first component of variance for the autosomal probes, accounting for 69.3% of the overall variation in methylation, distinguished samples based on tissue, while the second principal component (accounting for 12.7% of the overall variation), separated the species. Overall, an average of 14.5% (range of 8.2–26.1%, depending on the pairwise comparison) of the assayed promoter CpG sites were differentially methylated between tissues within a species, while an average of 8.6% of the CpG sites (range of 3.4–13.5%, depending on the tissue) were differentially methylated between humans and chimpanzees (at FDR<0.001). Reassuringly, these patterns recapitulate previous observations in human and mouse [26], [44].
Methylation patterns in T-DMRs are often conserved
We identified regions with tissue-specific patterns of methylation (T-DMRs [26], [42]) by analyzing the data from each species separately (Figure 3). Specifically, we modeled the methylation data (namely, the β values) from each autosomal CpG site independently, using a linear mixed-effects model with a fixed effect for the tissue and a random effect to account for variation between individuals. We tested for differences in methylation levels between tissues by using likelihood ratio tests within the framework of the linear model (see Methods). Using this approach, we identified 1,578 and 1,401 T-DMRs in humans and chimpanzees, respectively (at an FDR<0.001; Figure 3A; Table S1).
Tissue-specific methylation profiles are of interest because they may underlie tissue-specific patterns of gene expression levels. To test this hypothesis, we calculated, separately for each species, Pearson correlation values between promoter methylation profiles and the corresponding gene expression levels, across the three tissues. If methylation was consistently used to silence tissue-specific gene expression across the genome, we would expect to observe an abundance of negative correlations between the estimates of methylation and gene expression levels. However, when we considered the data for all genes that were expressed in at least one tissue, we found no evidence for an enrichment of negative correlations between methylation and gene expression levels (Figure 3B, Figure S6; 48% and 49% of the correlation values were negative in human and chimpanzee, respectively). In contrast, when we restricted the analysis to species-specific T-DMRs, we found an enrichment of negative correlations between methylation and gene expression levels (Figure 3B; 64% and 67% of correlation values were negative in human and chimpanzee, respectively; Fisher's exact P<10−16). This result suggests that T-DMRs underlie a subset of gene expression differences across tissues, a notion that is consistent with the important role played by DNA methylation in tissue differentiation in a wide range of species [42].
We then focused on the subset of T-DMRs with the same methylation pattern in both species. We found that 18–26% (depending on the tissue) of loci classified as T-DMRs in either human or chimpanzee are shared between the two species (Figure 3A, Table S2), a highly significant overlap compared to that expected by chance alone (hypergeometric distribution P values across all pairwise tissue comparisons <10−16). Importantly, the observation of a significant overlap in T-DMRs across species is robust with respect to the statistical cutoff used to classify T-DMRs (0.001≤FDR≤0.05; Table S2). Interestingly, when we considered correlations of methylation and gene expression levels only at conserved T-DMRs, we found an even more pronounced enrichment of negative correlations (Figure 3B and 3C; 72% of the correlation values were negative, regardless of species; Fisher's exact for an enrichment of negative correlations: P<10−23), suggesting that conservation of T-DMRs often relates to functionally important tissue-specific patterns of gene regulation.
It is perhaps interesting to note that we did not find a difference in the correlation of methylation and expression levels between T-DMR CpG sites that are located within or outside an annotated CpG island (as defined by [47]; Figure S7).
When we examined the functional annotations of genes associated with species-specific T-DMRs as well as conserved T-DMRs (using gene ontology annotations), we found an expected enrichment of genes annotated as important in ‘developmental’ processes, regardless of tissue (P<5×10−3; FDR<0.3; Table S3), congruent with the importance of epigenetic modification in tissue differentiation. We also found enrichments of tissue-specific biological processes, such as genes associated with cardiac muscle cell differentiation processes among heart T-DMRs (P<5×10−3; FDR<0.3), genes associated with embryonic organ morphogenesis and embryonic organ development processes among kidney T-DMRs (P<5×10−4; FDR<0.05), and genes associated with blood coagulation and with the regulation of body fluid levels (putatively involved in homeostatic functions) among liver T-DMRs (P<10−5; FDR<6×10−3 and P<10−4; FDR<0.007, respectively). The enrichment of genes associated with both developmental and tissue-specific processes among genes associated with T-DMRs is consistent with previous observations [27], [42]. Furthermore, when we considered only conserved T-DMRs, we observed a significant under-representation of genes associated with nucleic-acid and primary metabolic processes in all three tissues studied (all P<5×10−3; FDR<0.01; Table S4). This result suggests that the epigenetically-mediated tissue-specific regulation of these core processes tends to be conserved between humans and chimpanzees.
Inter-species differences in methylation
We next focused on the relationships between inter-species differences in methylation profiles and differences in gene expression levels between humans and chimpanzees. To estimate the relative contribution of changes in DNA methylation to inter-species differences in gene expression levels, we used linear regression analysis to account for promoter methylation effects (per autosomal CpG site) before analyzing the gene expression data from both species. We analyzed methylation and gene expression data in each tissue using a linear model framework similar to the one described in Blekhman et al. 2008 [8]. We then compared the evidence supporting an inter-species difference in gene expression levels before and after correcting for methylation profiles (see Methods for more details).
For the majority of genes (78%, 82%, and 77% in liver, kidney, and heart, respectively; Figure 4A), the evidence for a difference in expression level between the species was similar, regardless of whether or not methylation status was taken into account. For a small subset of genes (1%, 3%, and 2% in liver, kidney, and heart, respectively), we did not find compelling evidence for a difference in expression level between the species using the uncorrected expression level data, but after correcting for methylation levels using regression analysis, we rejected the null hypothesis of no inter-species differences in gene expression level (at an FDR<0.01). This observation, however, is unlikely to be biologically meaningful, since it is expected by chance alone (by permutation analysis; P>0.434 for all tissues; Figure S8).
In contrast, in all three tissues, we found a significant enrichment of genes for which the evidence for inter-species differences in expression level was compelling (FDR<0.01) before, but not after we corrected for the methylation levels (21%, 15%, and 21% in liver, kidney, and heart, respectively, permutation analysis yields P<0.001 for all tissues; Figure 4B and 4C). Based on the expectation of such a pattern by chance alone (by permutations – see Methods for details), we estimated that, in the three tissues we studied, inter-species differences in promoter DNA methylation might underlie as much as 12–18% of differences in gene expression levels between humans and chimpanzees.
When we analyzed the data considering only the sets of genes that have negative correlations between methylation and gene expression levels (as expected if methylation is used to silence gene expression), we found that 8.1%, 7.6%, and 8.8% of interspecies differences in gene expression levels in liver, kidney, and heart, respectively, might be explained by corresponding methylation differences. The extent to which inter-species gene expression differences might be explained by methylation differences between the species was similar regardless of whether the methylated site was within or outside an annotated CpG islands (Figure S9).
Discussion
We explored the extent to which putatively functional DNA methylation differences between tissues are conserved in humans and chimpanzees, and estimated the relative contribution of inter-species changes in methylation levels to gene expression differences between the two species. To do so, we collected DNA methylation and gene expression data from frozen human and chimpanzee primary tissue samples. While we chose to work with tissues that are relatively homogenous with respect to their cellular composition, we could not measure the precise composition or choose to work with particular cell types, because the samples were frozen. Similarly, we could not stage the tissues or control the environment of the donor individuals because the samples were collected post mortem. These are limitations shared by nearly all comparative molecular studies of primary tissues from humans and other apes (see [8], [9] for more detailed discussions of the limitations associated with studying gene regulation in primate tissues).
The challenge is therefore to focus on patterns in the data that should be robust with respect to the aspects of the study design that could not be controlled. For example, it is reasonable to expect that differences in environment, staging, and cellular composition across samples will tend to increase variation of measurements within, and especially between species. For that reason, our analysis of conserved inter-tissue gene expression differences and tissue-specific methylation patterns is likely to be conservative. Indeed, because of our inability to minimize environmental differences across the donor individuals, it is likely that we are underestimating the proportion of conserved inter-tissue gene expression differences and conserved T-DMRs.
In turn, when we focus on inter-species differences in DNA methylation and gene expression levels, it is important to note that our study design does not allow us to distinguish between regulatory differences due to either heritable or environmental effects. Studies in model organisms typically do so by controlling the environments of all subjects, a restriction we cannot apply when studying primate tissues. However, we have previously shown that estimates of differences in gene regulation between humans and chimpanzees based on six randomly sampled individuals are stable [8], [9]. Regardless of the underlying mechanism, it is likely that the analysis of the data uncovered mostly steady-state inter-species regulatory differences. Thus, even if differences in environments underlie a subset of the observed regulatory differences between humans and chimpanzees, our previous work suggests that it is likely that, in most cases, we capture the effects of general environmental differences between the species, not just between the samples used.
DNA methylation and differences in gene expression across tissues
We found a substantial degree of conservation of tissue-specific methylated regions in human and chimpanzee. This observation is not surprising given that previous studies found a marked conservation of T-DMRs between human and mouse, which are much more distantly related [26], [41], [43], [44]. On the other hand, 7.0%, 21.6%, and 23.8% of the kidney, heart, and liver T-DMRs, respectively (identified in either species), were differentially methylated (in the relevant tissue) between humans and chimpanzees, while only 3.3%, 8.0%, and 11.8% of non-TDMRs in these three tissues were differentially methylated between the two species (P<10−10 for all pairwise comparisons).
The conservation of T-DMR profiles yet the generally faster rate of inter-species change in promoter methylation at T-DMRs compared to non-T-DMRs are intriguing. These observations are difficult to explain by technical or uncontrolled aspects of the study design, because it is unlikely that those confounding factors would affect methylation at T-DMRs differently than at non-T-DMRs. Instead, it is likely that the different patterns truly reflect a functional difference between methylation at T-DMRs and at non-T-DMR CpG sites (in the studied tissues).
Though there is substantial evidence that DNA methylation levels upstream of genes are often inversely correlated with gene expression levels [24], [28], [31], recent studies proposed that methylation of promoters may play only a relatively minor role in the regulation of tissue-specific gene expression [34]. In particular, Maunakea et al. [48] posited that methylation of gene body regions (in regions that putatively serve as alternative promoters) might have a greater influence on regulatory differences across tissues. While we cannot use our data to ask about the relative importance of different types and locations of epigenetic marks to tissue-specific gene regulation, our observations strongly imply that any such debate would benefit from further investigation into the evolution of epigenetic profiles. Indeed, in addition to a faster rate of evolutionary change of the methylation profiles in T-DMRs, we found evidence for an enrichment of inverse correlations between inter-tissue gene expression patterns and promoter methylation profiles at genes associated with T-DMRs, but not when we considered all genes (the latter observation is consistent with the findings of Weber et al. [34] and Maunakea et al. [48]). Our results, therefore, imply that tissue-specific promoter methylation patterns may play especially important roles in regulating gene expression. The data also suggest that altered methylation levels, primarily at these sites, may underlie regulatory differences between species.
DNA methylation and inter-species differences in gene regulation
We estimated that as much as 12–18% (depending on the tissue) of inter-species differences in gene expression levels might be explained, at least in part, by changes in DNA methylation patterns. It is important to note that this statement is based on the proposed mechanism by which DNA methylation affects the rate of transcription and overall levels of gene expression [28], [31]. Though we did not perform experiments from which causality can be directly deduced, a causal relationship between changes in DNA methylation and gene regulation is strongly supported by previous studies (e.g., [24], [28], [31]). When we only consider negative correlations between methylation and gene expression levels to be indicative of a putative causal relationship, 8–9% of inter-species differences in gene expression levels might be explained by corresponding changes in DNA methylation.
However, other mechanisms are also likely [34], [43]. While DNA methylation is typically considered a silencing mechanism, high levels of methylation may be causally linked to increased gene expression levels. For example, the methylation of a repressor site could prevent the binding of repressor transcription factors, or enhancer transcription factors could favor binding to a methylated site rather than to the unmethylated site [49]–[51]. The observation of a small enrichment of positive correlations between methylation and expression when only T-DMRs are considered provides additional support for these types of mechanisms. Thus, perhaps as much as 12–18% of differences in gene expression levels between humans and chimpanzees might be explained by inter-species changes in DNA methylation.
Either way, our results suggest that DNA methylation differences in promoter regions might account for, at most, a modest proportion of inter-primate differences in gene expression levels (we confirmed that our estimates do not rely on arbitrary choices of specific statistical cutoffs; Tables S2 and S5). Many inter-species differences in promoter methylation are not associated with gene expression differences between the species. One explanation for that observation may simply be that these methylation patterns are not regulatory or functional. An alternative, more interesting possibility to consider, is that a subset of genes whose regulation differed between species later acquired modifications in nearby DNA methylation patterns to accommodate (or even partially counteract) the original expression level changes.
Since we assayed methylation using a pre-designed microarray, changes in DNA methylation in un-assayed genomic regions might explain additional regulatory differences between the species. In particular, while our assay focused on methylation at promoter regions, it has been recently shown that as a class, gene-body methylation profiles might explain a larger proportion of variation in gene expression levels than methylation profiles at currently annotated promoters [26], [48]. With the advent of new sequencing technologies, it will soon be feasible to extend our comparative approach to characterize genome-wide patterns of methylation.
In summary, we have taken some of the first steps towards characterizing variation in one mechanism that affects gene expression differences between closely related primate species [16], [17]. In a broader context, DNA methylation is just one of many mechanisms that have been posited to regulate gene expression levels [28], [31], [52]. In that sense, our study is a step towards the ultimate goal of understanding the relative importance of changes in different regulatory mechanisms to human evolution. Our observations indicate that at least 82% of gene expression differences between humans and chimpanzees (in the three studied tissues and specific promoter CpG sites examined) are not likely to be explained by differences in promoter DNA methylation.
Methods
DNA methylation data
We collected methylation data from the same human and chimpanzee liver, kidney, and heart tissue samples used in Blekhman et al. 2008 [8] (Figure S1; see Table S6 for details on the samples). DNA was extracted from each sample (6 human and 6 chimpanzee samples from each of the three tissues) in two independent technical replicates using the QIAamp DNA Mini Kit (Qiagen) (with the exception of chimpanzee sample CK2, for which DNA was only available for one replicate – see Table S4). The methylation profile of each sample was assayed using the Illumina HumanMethylation27 DNA Analysis BeadChip, which assays methylation at 27,578 CpG sites. Methylation array data are deposited to the NCBI GEO database under the accession number GSE26033 (http://www.ncbi.nlm.nih.gov/projects/geo/query/acc.cgi?acc=GSE26033).
To facilitate an unbiased comparison of methylation and gene expression levels in the human and chimpanzee samples, we first mapped the 27,578 50-bp Illumina probes to the human genome sequence (hg18) using BLAT [53] and MAQ [54]. We retained only the 26,690 probes that unambiguously mapped to a single location in the human genome with a maximum of two mismatches. These probes were then associated with the nearest gene using Ensembl gene annotation, and we retained only the subset of probes associated with genes that were represented on the multi-species gene expression microarray used by Blekhman et al. 2008 [8]. This resulted in the retention of 19,849 probes, associated with 11,059 genes. Finally, since the Illumina array was designed based on human genomic sequence, we limited our analysis to probes that were a perfect sequence match to a single location in the chimpanzee genome, by mapping the remaining 19,849 probes to the chimpanzee genome (panTro2) using BLAT [53] and MAQ [54]. We retained 10,575 probes that mapped uniquely to the chimpanzee genome with no sequence mismatches. This step ensures that our relative methylation measurements are not biased due to the effect of sequence mismatches on hybridization intensities. The resulting set of 10,575 probes is associated with 7,723 genes, which are present on every chromosome in the genome except for the Y-chromosome (Figure S10). The majority (97%) of retained probes are located within 2 kb of an annotated transcription start site of the associated gene (Figure S11). We note that a similar screen for probes that were a perfect match to the genomes of human, chimpanzee, and rhesus macaque resulted in the retention of only 1,944 probes (associated with 1,715 genes). For that reason, we limited our current study to a comparison between human and chimpanzee samples.
All samples were hybridized to the Illumina HumanMethylation27 DNA Analysis BeadChip at the Southern California Genotyping Consortium facility following standard manufacturer's instructions. Basic quality checks were performed using Illumina's BeadStudio software. Of the 10,575 probes we considered as the final dataset, 299 had missing data for one or more individuals and were discarded in all subsequent analyses. This resulted in 9,911 autosomal probes (corresponding to 7,291 genes) and 365 probes on the X-chromosome (corresponding to 266 genes). Since the probes map to distinct CpG island regions, which can affect downstream gene expression independently, we treated methylation levels from each CpG probe as distinct data points in all subsequent analyses. We further classified each probe as being located confidently within a CpG island region or outside of a strict CpG island region using the CpG Islands track information downloaded from UCSC [47].
For each sample, the methylation status at a probed location was summarized as: where M and U denote the signal emitted from the beads assaying the methylated and unmethylated versions at each site, respectively. Due to the number of samples being interrogated, it was necessary to hybridize the samples in two balanced batches. We observed a small difference in the mean β-value between batches, and corrected for this difference by standardizing the means across batches. After this correction, there was no further evidence for a batch effect.
To further assess the quality of the data, we calculated pairwise correlations between the β-values for all hybridized samples (Figure S2). As expected, technical replicates (which were independent DNA extractions) were the most highly correlated (36 comparisons; median r = 0.99), followed by samples from the same tissue and species (396 comparisons; median r = 0.98), samples from the same tissue across species (432 comparisons; median r = 0.97), samples from different tissues from the same species (864 comparisons; median r = 0.95), and samples from different tissues and different species (864 comparisons; median r = 0.93).
To look for evidence of imprinting in both humans and chimpanzees, we focused on a set of 27 genes (associated with 90 methylation probes) known to be imprinted based on the Imprinted Gene Catalog (IGC) at http://igc.otago.ac.nz/. To assess whether the patterns of DNA methylation at these imprinted genes were likely to occur by chance, we compared the observed proportion of hemi-methylated sites (defined as 0.3<β<0.7) to the distribution obtained by analyzing methylation patterns in 1000 randomly chosen sets of 90 methylation probes, associated with an average of 27 genes (range 26–28).
Gene expression data
Measurements of gene expression levels for all samples in our study were previously described by Blekhman et al. (2008) [8]. These data are available at the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/) under series accession number GSE11560. In that study, a multi-species microarray was used to estimate gene expression levels in cDNA samples from humans, chimpanzees, and rhesus macaques. The multi-species array includes orthologous probes for 18,109 genes, thus facilitating comparisons of gene expression levels between species without the confounding effects of sequence mismatches on hybridization intensities [8]. Since our current study focused only on the human and chimpanzee gene expression data, we re-normalized the expression data using only the human and chimpanzee probes on the array, using the same modified quantile normalization approach described in Blekhman et al. (2008) [8]. All further analyses used these re-normalized gene expression estimates. When examining the relationships between gene expression and methylation levels, we limited our analyses to genes that were either expressed in at least one tissue (for inter-tissue comparisons within a species) or expressed in at least one species (for the inter-species comparisons within a tissue), using a conservative threshold for defining expression, based on the entire distribution of expression values (normalized expression value of 8; see Figure S14 in Blekhman et al. (2008) [8]).
Statistical analysis
All statistical analyses were performed using the R statistical framework (http://www.r-project.org).
Identifying tissue-differentially methylated regions (T-DMRs)
To identify T-DMRs, we modeled the methylation level of each CpG site separately within both humans and chimpanzees using a linear mixed-effects model. Specifically, for each of the 9,913 probes (associated with 7,291 genes) located on the autosomal chromosomes, if yijk represents the β value for technical replicate k (k = 1 or 2), for individual j (j = 1,…,6), from tissue i (i = heart, liver, or kidney), we assume that: (1)where: Here, αi represents the mean methylation value at a given site in tissue i. To account for correlation between samples of the same tissue from different individuals, a random effect, ρij, which follows a N(0,σ2rand) distribution, is also included in the model.
To determine whether a CpG site was likely to fall within a T-DMR, we assessed how well the model (1) fitted the data under various parameterizations of μijk. The three types of parameterizations considered are:
In the simplest model (H0), the region's methylation value is assumed to be constant across all three tissues, while in the second alternative (H2) the methylation value is allowed to differ between all three tissues. The first alternative (H1) models the situation where the methylation level at the site of interest is constant in the two non-target tissues but differs in the target tissue. All models are fitted using a restricted maximum likelihood (REML) framework, and the maximum likelihoods were calculated.
In this study, we are interested in identifying sites whose methylation levels are best modeled by H1. To find such sites, we first used a likelihood-ratio test statistic (with one degree of freedom) to exclude sites where H2 provides a better fit to the data than H1 (specifically, if the likelihood-ratio p-value was less than 0.05, we removed these sites from the analysis). H2 provides a better fit for 1220 and 886 (in humans and chimpanzees, respectively) of the total 9911 autosomal CpG sites. For the remaining positions, we examined whether there was significant evidence to reject H0 in favor of H1 using a likelihood-ratio test statistic (which we compared to a χ2 distribution with 1 degree of freedom). We corrected for multiple testing using the FDR approach of Storey and Tibshirani [55].
Gene ontology analysis for T-DMRs
We used GeneTrail (http://genetrail.bioinf.uni-sb.de) [56] to test for enrichments of functional annotations among different classes of T-DMRs. In all tests, we used a background set of genes that were present in our study and classified as expressed in at least one tissue (conditional on a normalized expression value of 8). The tests were performed using all GO categories and KEGG pathways. We calculated p-values using a Hyper-geometric distribution and report false discovery rates for each p-value.
A joint analysis of methylation and gene expression levels
To examine whether changes in gene expression levels between humans and chimpanzees (within each tissue) can be explained by inter-species differences in methylation levels, we extended the linear mixed-effects model framework described in Blekhman et al. (2008) [8] to include methylation as a covariate. However, since we have to correct the multi-species array data for probe-effects [8], it is difficult to interpret the methylation coefficient when it is added directly to the model, since it is confounded with the probe effects. Consequently, we used an alternative approach in which we used regression to correct for the methylation effect. Specifically, for each gene-tissue combination, we tested for differences in expression level between human and chimpanzee after regressing out the following effects:
-
Expression microarray probe effects only
-
Expression microarray probe effects and CpG-specific methylation levels
To do this, we used a fully parameterized model where gene expression probe effects, CpG-probe methylation values, and species effects were explanatory variables. Additionally, a random effect was used to account for variability between biological replicates. Specifically, if ysroi denotes the normalized log2 intensity expression value for individual i (i = 1,…,6), from species s (s = human or chimpanzee) measure at probe r (r = 1,…,7), which is derived from species o, we assume that: where: Here, μs denotes the species effect, πro is a fixed-effect representing the probe effect for each individual probe within a probe-set and the composition effect of species-specific orthologous probes, and κsro is a fixed-effect representing the attenuation of hybridization intensities due to sequence mismatches between species of RNA and a species-specific derived probe, which are different for each individual probe within a probe set (see [8] for more details). Additionally, γsi is a random effect (following a N(0,σ2rand) distribution) and βsi denotes the β value for the methylation probe of interest for individual i from species s. Upon fitting this model, using the lmer package within the R statistical framework, estimates of the parameters and the residuals were obtained. To obtain corrected measures of expression for each individual from each species, when probe and methylation effects are regressed out (scenario 2), we defined . When we only regressed out probe effects (scenario 1), the corrected values are defined as . In both of these scenarios, once the corrected data were obtained, we tested for differences in gene expression levels as follows. If, for each tissue-gene combination, xsik denotes the (corrected) level of expression for replicate k of individual i from species s, we modeled these data as follows:where:Here, αs is a species effect, and ρsi is a random individual effect. Subsequently, to test for inter-species differences in expression levels, we compare the following hypotheses:Here, the null model assumes equal expression level between the two species, and the alternative assumes different expression levels. Evidence against the null model was determined using a likelihood-ratio test statistic (compared against a chi-squared distribution with one degree of freedom). By performing this analysis independently for each CpG-gene combination in all tissues, we obtained a p-value indicating the strength of the evidence against the null hypothesis, before (under scenario 1 above) and after (under scenario 2 above) accounting for the region's DNA methylation status. By comparing these p-values, we were able to identify genes within each tissue where the difference in expression level between human and chimpanzee was likely explained by inter-species differences in DNA methylation.
To assess the statistical significance of our observations, we permuted the methylation values for a given gene across all individuals (maintaining replicate correlations, but allowing labels to permute across species classifications). Subsequently, we repeated the analysis described above to obtain an expected distribution of discrepancies between the methylation-corrected and uncorrected data. We performed 1000 permutations and p-values were calculated based on the number of times we observed as many or more discrepancies in the permuted compared to the real data.
In order to estimate the proportion of genes for which methylation differences might underlie gene expression differences, we treated the medians of the permutation distributions from each tissue as background levels. For each tissue, we then subtracted the background level from the observed proportion of genes with reduced evidence for inter-species differences in gene expression levels, once methylation was taken into account.
Supporting Information
Zdroje
1. KingM-C
WilsonAC
1975 Evolution at Two Levels in Humans and Chimpanzees. Science 188 107 116
2. BrittenRJ
DavidsonEH
1969 Gene Regulation for Higher Cells: A Theory. Science 349 357
3. EnardW
KhaitovichP
KloseJ
ZöllnerS
HeissigF
2002 Intra- and interspecific variation in primate gene expression patterns. Science 296 340 343
4. KhaitovichP
MuetzelB
SheX
LachmannM
HellmannI
2004 Regional patterns of gene expression in human and chimpanzee brains. Genome Research 14 1462 1473
5. KhaitovichP
HellmannI
EnardW
NowickK
LeinweberM
2005 Parallel patterns of evolution in the genomes and transcriptomes of humans and chimpanzees. Science 309 1850 1854
6. CáceresM
LachuerJ
ZapalaMA
RedmondJC
KudoL
2003 Elevated gene expression levels distinguish human from non-human primate brains. Proc Natl Acad Sci USA 100 13030 13035
7. KaramanMW
HouckML
ChemnickLG
NagpalS
ChawannakulD
2003 Comparative analysis of gene-expression patterns in human and African great ape cultured fibroblasts. Genome Research 13 1619 1630
8. BlekhmanR
OshlackA
ChabotAE
SmythGK
GiladY
2008 Gene regulation in primates evolves under tissue-specific selection pressures. PLoS Genet 4 e1000271 doi:10.1371/journal.pgen.1000271
9. BlekhmanR
MarioniJC
ZumboP
StephensM
GiladY
2010 Sex-specific and lineage-specific alternative splicing in primates. Genome Research 20 180 189
10. PrabhakarS
ViselA
AkiyamaJA
ShoukryM
LewisKD
2008 Human-specific gain of function in a developmental enhancer. Science 321 1346 1350
11. BabbittCC
SilvermanJS
HaygoodR
ReiningaJM
RockmanMV
2010 Multiple Functional Variants in cis Modulate PDYN Expression. Molecular Biology and Evolution 27 465 479
12. WarnerLR
BabbittCC
PrimusAE
SeversonTF
HaygoodR
2009 Functional consequences of genetic variation in primates on tyrosine hydroxylase (TH) expression in vitro. Brain Research 1288 1 8
13. LoiselDA
RockmanMV
WrayGA
AltmannJ
AlbertsSC
2006 Ancient polymorphism and functional variation in the primate MHC-DQA1 5′ cis-regulatory region. Proc Natl Acad Sci USA 103 16331 16336
14. RockmanMV
HahnMW
SoranzoN
ZimprichF
GoldsteinDB
2005 Ancient and recent positive selection transformed opioid cis-regulation in humans. Plos Biol 3 e387 doi:10.1371/journal.pbio.0030387
15. PollardKS
SalamaSR
LambertN
LambotM-A
CoppensS
2006 An RNA gene expressed during cortical development evolved rapidly in humans. Nature 443 167 172
16. BlekhmanR
OshlackA
GiladY
2009 Segmental Duplications Contribute to Gene Expression Differences Between Humans and Chimpanzees. Genetics 182 627 630
17. ChabotA
ShritRA
BlekhmanR
GiladY
2007 Using reporter gene assays to identify cis regulatory differences between humans and chimpanzees. Genetics 176 2069 2076
18. BabbittCC
FedrigoO
PfefferleAD
BoyleAP
HorvathJE
2010 Both Noncoding and Protein-Coding RNAs Contribute to Gene Expression Evolution in the Primate Brain. Genome Biology and Evolution 2 67 79
19. GorenA
SimchenG
FibachE
SzaboPE
TanimotoK
2006 Fine Tuning of Globin Gene Expression by DNA Methylation. PLoS ONE 1 e46 doi:10.1371/journal.pone.0000046
20. HeardE
DistecheCM
2006 Dosage compensation in mammals: fine-tuning the expression of the X chromosome. Genes & Development 20 1848 1867
21. HeardE
ClercP
AvnerP
1997 X-Chromosome Inactivation in Mammals. Annu Rev Genet 31 571 610
22. SadoT
FennerMH
TanSS
TamP
ShiodaT
2000 X inactivation in the mouse embryo deficient for Dnmt1: distinct effect of hypomethylation on imprinted and random X inactivation. Dev Biol 225 294 303
23. LiE
BeardC
JaenischR
1993 Role for DNA methylation in genomic imprinting. Nature 366 362 365
24. ReikW
2007 Stability and flexibility of epigenetic gene regulation in mammalian development. Nature 447 425 432
25. EggerG
LiangG
AparicioA
JonesPA
2004 Epigenetics in human disease and prospects for epigenetic therapy. Nature 429 457 463
26. IrizarryRA
Ladd-AcostaC
WenB
WuZ
MontanoC
2009 The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nature Genetics 41 178 186
27. IllingworthR
KerrA
DeSousaD
JorgensenH
EllisP
2008 A Novel CpG Island Set Identifies Tissue-Specific Methylation at Developmental Gene Loci. Plos Biol 6 e22 doi:10.1371/journal.pbio.0060022
28. JaenischR
BirdA
2003 Epigenetic regulation of gene expression: how the genome intergrates intrinsic and environmental signals. Nature Genetics 33 245 254
29. LiE
BestorTH
JaenischR
1992 Targeted Mutation of the DNA Methyltransferase Gene Results in Embryonic Lethality. Cell 69 915 926
30. OkanoM
BellDW
HaberDA
LiE
1999 DNA Methyltransferases Dnmt3a and Dnmt3b Are Essential for De Novo Methylation and Mammalian Development. Cell 99 247 257
31. MurrellA
RakyanVK
BeckS
2005 From genome to epigenome. Hum Mol Genet 14 Spec No 1 R3 R10
32. SteinR
RazinA
CedarH
1982 In vitro methylation of the hamster adenine phosphoribosyltransferase gene inhibits its expression in mouse L cells. Proc Natl Acad Sci USA 79 3418 3422
33. HansenRS
GartlerSM
1990 5-Azacytidine-induced reactivation of the human X chromosome- linked PGK1 gene is associated with a large region of cytosine demethylation in the 5′ CpG island. Proc Natl Acad Sci USA 87 4174 4178
34. WeberM
HellmannI
StadlerMB
RamosL
PääboS
2007 Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nature Genetics 39 457 466
35. ThomsonJP
SkenePJ
SelfridgeJ
ClouaireT
GuyJ
2010 CpG islands influence chromatin structure via the CpG-binding protein Cfp1. Nature 464 1082 1086
36. WattF
MolloyPL
1988 Cytosine methylation prevents binding to DNA of a HeLa cell transcription factor required for optimal expression of the adenovirus major late promoter. Genes & Development 2 1136 1143
37. FengS
CokusSJ
ZhangX
ChenP-Y
BostickM
2010 Conservation and divergence of methylation patterning in plants and animals. Proceedings of the National Academy of Sciences 107 8689 8694
38. EnardW
FassbenderA
ModelF
AdorjanP
PaaboS
2004 Differences in DNA methylation patterns between humans and chimpanzees. Current Biology 14 R148 R149
39. Gama-SosaMA
MidgettRM
SlagelVA
GithensS
KuoKC
1983 Tissue-specific differences in DNA methylation in various mammals. Biochimica et Biophysica Acta 740 212 219
40. ZemachA
McDanielIE
SilvaP
ZilbermanD
2010 Genome-Wide Evolutionary Analysis of Eukaryotic DNA Methylation. Science 1 7
41. IgarashiJ
MuroiS
KawashimaH
WangX
ShinojimaY
2009 Quantitative analysis of human tissue-specific differences in methylation. Biochemical and Biophysical Research Communications 376 658 664
42. RakyanVK
DownTA
ThorneNP
FlicekP
KuleshaE
2008 An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Research 18 1518 1529
43. EckhardtF
LewinJ
CorteseR
RakyanVK
AttwoodJ
2006 DNA methylation profiling of human chromosomes 6, 20 and 22. Nature Genetics 38 1378 1385
44. KitamuraE
IgarashiJ
MorohashiA
HidaN
OinumaT
2007 Analysis of tissue-specific differentially methylated regions (TDMs) in humans. Genomics 89 326 337
45. GibbsJR
van der BrugMP
HernandezDG
TraynorBJ
NallsMA
2010 Abundant quantitative trait Loci exist for DNA methylation and gene expression in human brain. PLoS Genet 6 e1000952 doi:10.1371/journal.pgen.1000952
46. FarcasR
SchneiderE
FrauenknechtK
KondovaI
BontropR
2009 Differences in DNA methylation patterns and expression of the CCRK gene in human and nonhuman primate cortices. Mol Biol Evol 26 1379 1389
47. Gardiner-GardenM
FrommerM
1987 CpG Islands in Vertebrate Genomes. J Mol Biol 196 261 282
48. MaunakeaAK
NagarajanRP
BilenkyM
BallingerTJ
D'souzaC
2010 Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature 466 253 257
49. RishiV
BhattacharyaP
ChatterjeeR
RozenbergJ
ZhaoJ
2010 CpG methylation of half-CRE sequences creates C/EBP{alpha} binding sites that activate some tissue-specific genes. Proc Natl Acad Sci USA 107 20311 20316
50. SunL
HuangL
NguyenP
BishtKS
Bar-SelaG
2008 DNA methyltransferase 1 and 3B activate BAG-1 expression via recruitment of CTCFL/BORIS and modulation of promoter histone methylation. Cancer Res 68 2726 2735
51. GiusD
CuiH
BradburyCM
CookJ
SmartDK
2004 Distinct effects on gene expression of chemical and genetic manipulation of the cancer epigenome revealed by a multimodality approach. Cancer Cell 6 361 371
52. BergerSL
KouzaridesT
ShiekhattarR
ShilatifardA
2009 An operational definition of epigenetics. Genes & Development 23 781 783
53. KentWJ
2002 BLAT--The BLAST-Like Alignment Tool. Genome Research 12 656 664
54. LiH
RuanJ
DurbinR
2008 Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Research 18 1851 1858
55. StoreyJD
TibshiraniR
2003 Statistical significance for genomewide studies. Proc Natl Acad Sci U S A 100 9440 9445
56. BackesC
KellerA
KuentzerJ
KneisslB
ComtesseN
2007 GeneTrail--advanced gene set enrichment analysis. Nucleic Acids Res 35 W186 192
Štítky
Genetika Reprodukční medicínaČlánek vyšel v časopise
PLOS Genetics
2011 Číslo 2
- Primární hyperoxalurie – aktuální možnosti diagnostiky a léčby
- Srdeční frekvence embrya může být faktorem užitečným v předpovídání výsledku IVF
- Akutní intermitentní porfyrie
- Vztah užívání alkoholu a mužské fertility
- Šanci na úspěšný průběh těhotenství snižují nevhodné hladiny progesteronu vznikající při umělém oplodnění
Nejčtenější v tomto čísle
- Meta-Analysis of Genome-Wide Association Studies in Celiac Disease and Rheumatoid Arthritis Identifies Fourteen Non-HLA Shared Loci
- MiRNA Control of Vegetative Phase Change in Trees
- The Cardiac Transcription Network Modulated by Gata4, Mef2a, Nkx2.5, Srf, Histone Modifications, and MicroRNAs
- Genome-Wide Transcript Profiling of Endosperm without Paternal Contribution Identifies Parent-of-Origin–Dependent Regulation of