- Research article
- Open Access
CpG islands under selective pressure are enriched with H3K4me3, H3K27ac and H3K36me3 histone modifications
BMC Evolutionary Biology volume 13, Article number: 145 (2013)
Histone modification is an epigenetic mechanism that influences gene regulation in eukaryotes. In particular, histone modifications in CpG islands (CGIs) are associated with different chromatin states and with transcription activity. Changes in gene expression play a crucial role in adaptation and evolution.
In this paper, we have studied, using a computational biology approach, the relationship between histone modifications in CGIs and selective pressure in Homo sapiens. We considered three histone modifications: histone H3 lysine 4 trimethylation (H3K4me3), histone H3 lysine 27 acetylation (H3K27ac) and histone H3 lysine 36 trimethylation (H3K36me3), and we used the publicly available genomic-scale histone modification data of thirteen human cell lines. To define regions under selective pressure, we used three distinct signatures that mark selective events from different evolutionary periods. We found that CGIs under selective pressure showed significant enrichments for histone modifications.
Our result suggests that, CGIs that have undergone selective events are characterized by epigenetic signatures, in particular, histone modifications that are distinct from CGIs with no evidence of selection.
CpG islands (CGIs) are unmethylated segments of a genome that have an increased level of CpG dinucleotides and a high GC content [1, 2]. In the human genome, most CGIs are either inside or close to the promoter regions of genes. Typically these CGIs occur at or close to transcription start sites (TSSs) . It is well established that CpG sites in promoter CGIs are undermethylated in expressed genes, while hypermethylation of promoter CpG sites is associated with gene silencing . Others CGIs that are distant from known TSSs have been found in intergenic, 3’ and intragenic regions .
There is an extensive literature demonstrating that structural modifications to chromatin, along with CGI methylation, contribute to the functional output of related genes . The N-terminal tails of histone proteins can be modified covalently by small molecules (for example, phosphorylation, acetylation, methylation) and by macromolecules (for example, ubiquitination, sumoylation etc.). The precise environment of the CGI chromatin that controls gene regulation is not definitively established. The general understanding is that by altering the state of the CGI chromatin, histone modification can regulate access of the transcription machinery to particular DNA sequences . Of all the possible histone modifications, methylation of the lysine or arginine residues has received the main attention. These modifications can activate or repress the associated genes depending on which lysine or arginine residues are methylated . Methylation of histone H3 at lysine 9 (H3K9) or lysine 27 (H3K27) is considered to be a repressive mark . In contrast, H3K4me3, perhaps the best established epigenetic marker, is robustly associated with activation of transcription . In mammals, the trimethylation of H3K4 can be catalyzed by different histone methyltransferases, such as MLL1 or ASH1L [10, 11]. The majority of H3K4me3 sites overlap with the 5’ ends of annotated human genes  and several studies have reported the inverse correlation between two epigenetic marks, DNA methylation and H3K4me3 [13, 14]. The H3K4me3 mark also plays a crucial role in mammalian development , and its alteration has been found to be associated with cancer and other diseases [16–18]. In addition, both H3K27ac and H3K36me3, which are known as a promoter mark  and a gene body mark , respectively, are associated with transcriptional activation [19, 21].
Alterations in gene regulation are thought to play an important role both in adaptation and evolution . A recent report proposed that differences in gene expression levels among primates are associated with the changes in H3K4me3 . Moreover, another recent study identified human-specific changes in H3K4me3 levels at TSSs and related regulatory sequences in comparison with chimpanzees and macaques . Besides, the Encyclopedia of DNA Elements (ENCODE) project is studying different functional elements of human genome including regions of histone modifications. In particular they assayed chromosomal locations for 12 histone modification in 46 different cell types . In a previous study , we demonstrated that CGIs under selective pressure are hypomethylated compared to the CGIs in other regions of the genome. In this study, we explored the relationship between selective pressure signatures of and histone modification (H3K4me3, H3K27ac and H3K36me3) enrichment in CGIs. We used the genome-wide histone modification data of thirteen human cell lines produced by the ENCODE consortium [27, 28]. To define regions under selective pressure we used three distinct methods  that are able to detect both recent and ancient selective pressure events .
We analyzed thirteen cell lines from the ENCODE/Broad Institute, derived from nine normal and four cancer tissues, respectively. A list of features for each considered cell line is presented in Additional file 1: Table S1. For each cell line, we downloaded histone modification data for H3K4me3, H3K27ac and H3K36me3 marks.
We used the “Peaks Signal” (PS), representing regions of statically significant enrichment of a specific histone modification (see Materials and Methods). We downloaded genomic coordinates of 27718 unique CGIs defined according to criteria described in the University of California Santa Cruz Genome Browser (UCSC GB) (http://genome.ucsc.edu/) (see Materials and Methods). For each cell line, we estimated the number of CGIs containing at least one PS of histone modification and found, on average, 15478, 11903 and 10182 CGIs containing PSs of H3K4me3, H3K27ac and H3K36me3, respectively.
To identify genomic regions that may have undergone selective pressure we used three different approaches that are sensitive to selective pressure events that occurred in distinct evolutionary epochs.
The first method uses the per-continent “integrated Haplotype Score” (iHS)  and marks recent positive selection (see Materials and Methods). Using the iHS we identified 586 genomic regions that have putatively undergone recent selective pressure. We denoted these regions as “high iHS regions” (HIRs). Within the HIRs regions we found 2545 CGIs.
The second approach is based on a comparison between Homo sapiens and Neanderthal genomes (see Materials and Methods). The selective sweep scan score (S score) was used to identify regions of the human genome with a strong signal for depletion of Neanderthal-derived alleles. This score, when negative, may indicate an episode of positive selection in early humans . We found 212 genomic regions with a significant negative score (5% lowest S regions, hereafter denoted as 5LSRs) containing 348 CGIs.
In the third approach, we looked for sequences that were conserved across ten primate genomes. These sequences are the so-called “Conserved Elements” (CEs) (see Materials and Methods) and they mark ancient selective pressure events. We downloaded 725627 CEs and used them to search for CGIs that contain CEs . We identified 13288 unique CGIs that contained at least one CE.
We, then, computed the fraction of CGIs containing histone modification marks that show signatures of natural selection (HIRs, CEs and 5LSRs), and compared it with an analogous quantity computed for CGIs shown to have no signals of selective pressure. The presence of a possible enrichment/diminishment, defined as the ratio of the percentages of the above two groups, was assessed by means of a hypergeometric test (see Material and Methods).
We found a significant enrichment of H3K4me3 and H3K27ac markers for all three signatures of selection in almost all cell lines (Figure 1, Figure 2 and Additional file 2: Table S2) while for H3K36me3 the enrichment reached significance only for the CE signature (Figure 3 and Additional file 2: Table S2).
In order to understand if the results are due to the same CGIs identified by the three different methods, we estimated the overlaps among the corresponding CGIs lists (see Additional file 3). As shown in the figure, the different sets of CGIs identified are not essentially included one in the other.
We investigated whether or not these differences were dependent on the position of the CGIs in various genomic regions. To do this we followed the same approach described by Medvedeva et al.  dividing the CGIs into four groups according to their positions with respect to genes: at the 5’ end of a gene, in the intragenic region, at the 3’ end of a gene, and in the intergenic region. Results of this analysis are presented as additional information (see Additional file 4: Table S3) and summarized below.
Analysis of 5’ CGIs demonstrated the same significant enrichment pattern as seen in the overall analysis with significant enrichment of H3K4me3 and H3K27ac (Additional file 5 and Additional file 6), in almost all cell lines for all signatures of selection, and significant enrichment of H3K36me3 in all cell lines for the CE signature only (Additional file 7).
Both intragenic and 3’ CGIs were significantly enriched for H3K36me3 in the majority of cells lines (nine and ten out of the thirteen, respectively) for the CE signature (Additional file 8 and Additional file 9), while analysis of the other markers did not reach significance in almost all other cases (Additional file 10, Additional file 11, Additional file 12 and Additional file 13).
Finally, regarding intergenic CGIs we found a significant enrichment in all cell lines for all considered markers for the CE signature and in twelve out of thirteen cell lines for H3K4me3 in the HIR signature (Additional file 14, Additional file 15 and Additional file 16).
Evolutionary forces analysis
Two major evolutionary forces result in enriched CpG content: one is based on low levels of DNA methylation and, consequentially, deamination; and the other is biased gene conversion (BGC), which acts to repair TG mismatches caused by the deamination of methyl-cytosine . According to the role that these two forces play in CGI maintenance, CGIs can be classified as hypo-deaminated CGIs or BGC CGIs. We examined whether or not the relationship that we found between selective pressure and histone mark enrichment was present in both classes.
We found that both hypo-deaminated and BCG CGIs showed an enrichment of all markers in the CE signature in all cell lines, while only BGG CGIs showed significant enrichment of H3K4me3 and H3K27ac in the HIR signature in all cell lines (Additional file 17, Additional file 18, Additional file 19, Additional file 20, Additional file 21, Additional file 22 and Additional file 23: Table S4).
Expression divergence between humans and chimps and histone modifications
Cain et al.  have identified genes whose expression levels differ between species. In particular, by using their approach we have classified 1888 genes as Differentially Expressed (DE) between humans and chimpanzees, at a FDR of 0.050 (see Materials and Methods) among the whole set of 12559 genes considered in their study . Inside this class, and by applying a GREAT analysis  we associated corresponding sets of genes to CE CGIs (7436), HIR CGIs (1190) and 5SLR CGIs (214) (see Materials and Methods). As shown in Additional file 24: Table S5, the CE CGI genes are significantly enriched for DE, while HIR CGI genes, even if they do not match our conservative threshold for statistical significance of 0.001, show a p-value = 0.008. No significant signal is present in 5SLR CGI genes. In the Additional file 25: Table S6 we presented the histone modification enrichments referred to the different classes of genes labeled by their selective signal (present or not), and by their possible membership to the class of DE. As clear from the table, we found a significant enrichment of histone modifications in DE + CE genes, with respect to DE – CE class (defined as the set of elements in DE not belonging to CE), and this is independent of the histone modification chosen. The same is not true when we compare the enrichment in DE + CE genes, with respect to CE – DE. This suggests that CE plays the central role for the enrichment. The same is true for HIR if we only consider H3K4me3, whereas we get at most a marginal statistical significance if we take H3K27ac and H3K36me3. No significance at all is found for 5SLR.
As discussed in Materials and Methods, we assumed the conservative approach to define as PS genes the ones associated to the list of PS CGIs that were simultaneously present in all cell lines. In order to understand if our results depend on such choice we have also considered the other extreme case, namely labeling as PS genes the ones associated to a PS CGIs in at least one cell line. This approach confirms the previous finding (data not shown).
In this study, we investigated the hypothesis that CGIs under selective pressure are enriched with histone modifications that are associated with gene activation. To do this, we analyzed data from thirteen human cell lines for three well-known histone modifications (H3K4me3, H3K27ac and H3K36me3) to explore their relationship with both recent and ancient events of selective pressure.
H3K4me3 and H3K27ac are epigenetic marks that are generally associated with gene activation [9, 19] while H3K36me3 is associated with transcriptional elongation . Moreover, H3K4me3 and H3K27ac are evolutionarily conserved among species  and negatively correlated with DNA methylation [13, 14, 36, 37]. Also H3K36me3 in exons is found to be conserved between human and mouse .
Using the entire set of human CGIs we found that the CGIs associated with signatures of selective pressure were significantly enriched with H3K4me3 and H3K27ac in almost all considered cell lines. H3K36me3, on the other hand, showed a significant enrichment in global CGIs only in CE regions: this could be due to a small sample size effect (Additional file 2). These findings support a previous study in which we found that CGIs located in regions under selective pressure are more protected from DNA methylation compared the CGIs in other genomic regions . In the same study, we found that CGIs under selective pressure show a lower SNP content as well. Here we checked two parameters more: C+G content and CGIs’ length, and we found that CE CGIs have a slight but statistically higher G+C content compared with the remaining CGIs (mean = 0.689 vs. 0.683, t-test p < 2.2 10-16), and (adopting the classification described by Elango and Soojin ) an enrichment of long (>2000 bp) CGIs (Fisher’s Exact Test p < 2.2 10-16). We did not find significant differences in length and G+C content for HIR and 5LSR CGIs.
When we divided CGIs according to their positions with respect to the genes, we found that the statistical differences between CGIs with and without signatures of selective pressure were clearest for CGIs located in the 5’ regions for both H3K4me3 and H3K27ac. This result is intriguing in the light of the well established evidence that CGIs at the 5’ ends of genes are involved mainly in the control of gene expression . It is also possible that the small sample size led to a lack of statistical confidence in the results for CGIs in other positions. We noticed a different behavior for H3K36me3. H3K36me3 was the only mark to be enriched in 3’ and intragenic CGIs in CE regions for majority of cell lines; this finding is again intriguing considering that H3K36me3 is reported to be a gene body mark [8, 20]. In a recent study, H3K36me3 mark was found to be significantly associated with alternative splicing . It is well known that, alternative splicing is a key reason for protein diversity in higher eukaryotes . It has been a fundamental question in evolutionary study, how species having the similar repertoires of protein-coding genes differ strikingly at the phenotypic level. A very recent study proposed a link between alternative splicing and species-specific phenotypic differences among vertebrate species .
Two different evolutionary processes, namely hypo-deamination and BCG, are involved in the generation and maintenance of CGIs . The majority of hypo-deaminated CGIs are usually unmethylated while most BGC CGIs are constitutively methylated and clustered in subtelomeric regions. We found H3K4me3, H3K27ac and H3K36me3 enrichment in CGIs in CE regions, independently of the evolutionary process involved in their generation. Since CGIs belonging to these two groups differ in their DNA methylation levels, our finding seems to suggest that the difference we found was quite independent of the DNA methylation status.
The impact of natural selection on functional elements in human genome is also addressed in the last report from ENCODE project . In that case the authors focused their attention mainly on the relationship between negative selection and a subset of functional elements but they did not specifically address histone modifications. Positive selection, on the other hand, was addressed in a recent work by Vernot et al.  who studied the impact of this kind of selective pressure on DNase I peaks.
Cain et al  identified a list of genes that were differentially expressed between humans and other primates. They proposed that epigenetic changes could be, at least in part, involved in these differences. When we compared this list with the list of genes associated to the CGIs, we found a significant enrichment of differentially expressed genes in the CE CGI genes list, while the HIR CGI genes showed a p-value (0.008) near our conservative threshold for statistical significance (0.001). This seems to suggest that genes transcriptionally different between species are more likely located near CGIs with signature of selective pressure. In literature it is growing the evidence that epigenetic mechanisms provide a significant source of phenotypic variation that, in turn, can cause evolutionary novelty and potentially influence adaptation and evolution. Although the exact evolutionary significance of our results need further experiments to be completely defined, our data seem to support the idea of a close connection between adaptation, evolution and epigenetic mechanisms.
It has been hypothesized that CGIs are fundamental regulatory structures that have evolved under selection in genomes where DNA methylation plays a regulatory role [24, 45]. In particular, CGIs act as a platform where chromatin modifications and additional signaling help to define the functional output of the respective genes.
The present analysis concerning H3K4me3, H3K27ac and H3K36me3 enrichment in CGIs under selective pressure, supplements the findings of a previous paper by Cocozza et al. . In that study, we demonstrated a DNA undermethylation of CGIs under selective pressure. It is well established that a complex, perhaps bidirectional, crosstalk exists between DNA methylation and histone modification  suggesting that these two epigenetic mechanisms are, at some extent, dependent one each other. The overall picture emerging from the two studies is that CGIs under selective pressure seem to share definite epigenetic features.
To our knowledge, the present study is the first report addressing the relationship between histone modifications and natural selection and the overall framework emerging from our analyses support the hypothesis that CGIs that have experienced selection could be characterized by distinct epigenetic signatures.
Analyzing thirteen human cell lines, we found H3K4me3, H3K27ac and H3K36me3 enrichment in the CGIs that experienced selective events. Further studies using other epigenetic marks could help to clarify the relation between epigenetic modification and selective pressure in human genome.
CGIs coordinates were downloaded from the “CpgIslandExt” track of the UCSC GB (http://genome.ucsc.edu/). The CGIs in this track were predicted by searching the human genome assembly (GRCh37/hg19) sequence, scoring each dinucleotide and identifying maximally scoring segments. In this dataset, a CpG island was defined according to the following criteria: i) GC content of 50% or greater, ii) length of at least 200 bp, and iii) observed CpG / expected CpG ratio greater than 0.6. The CGI set that we obtained consisted of 27718 CGIs (this excluded the CGIs in the data related to the alternative haplotype sequences).
5’, intragenic, 3’, and intergenic CGIs
We used the classification system that was described previously by Medvedeva et al.  in which the CGIs were classified according to their locations. Thus, the CGIs were classified into four classes:
5’ CGIs - located in the 5’ flank region (3 kb upstream the TSS), the 5’ UTR-exon, the 5’ UTR-intron, the initial coding exon or the initial intron.
Intragenic CGIs are located in the internal exons and introns.
3’ CGIs are located in the final exon, the final introns, the 3’ UTR-exon or in the 3’ UTR-intron.
Intergenic CGIs are located at least 3 kb upstream or downstream from any known gene.
Hypo-deaminated and biased gene conversion (BGC) CGIs
Two sets of CGIs were described by Cohen et al.  using a new parameter-rich evolutionary model in combination with high resolution DNA methylation data to study the origin of the CpG repertoire in primate genomes (marmoset, rhesus, orangutan, chimp and human). Following a clustering analysis, they observed that most CGIs were constitutively unmethylated and underwent slow C-to-T deamination. They denoted this group as hypo-deaminated CGIs. In contrast, another class of CGI was constitutively methylated with a rapid deamination rate and was termed as BGC CGIs. For our analysis, we considered the 9091 hypo-deaminated and 4782 BGC CGIs from the UCSC CGIs sample.
Histone modification data
The histone modification (H3K4me3, H3K27ac and H3K36me3) data for thirteen human cell lines (HUVEC, Monocytes-CD14+_RO01746 (CD14+), HMEC, HSMM, HSMMtube, NH-A, NHDF-Ad, NHEK, NHLF, K562, HeLa-S3, HepG2 and Dnd41) were downloaded from the “Broad histone” track of the UCSC GB. This track contains genome-wide histone modification data of different cell lines, generated using ChIP-seq high-throughput sequencing as a part of the ENCODE project . In this study, we used the “Peaks Signal” (PS), which identifies discrete intervals of ChIP-seq fragment enrichment. In particular, we considered the CGIs in our sample that contained at least one PS.
Integrated haplotype score (iHS)
The iHS belongs to the Extended Haplotype Homozygosity statistic “family”  and is a marker of recent positive selection . The iHS measures the decay of identity, as a function of distance, of haplotypes that carry a specified “core” allele. We downloaded the iHS normalized values from the “HGDP iHS” track of the UCSC GB. The scores were calculated using SNPs genotyped in 1043 individual taken from 53 populations worldwide by the Human Genome Diversity Project in collaboration with the Centre d’Etude du Polymorphisme Humain (HGDP-CEPH). The 53 populations were divided into seven continental groups: Africa (Bantu populations only), Middle East, Europe, South Asia, East Asia, Oceania and the Americas. For each population group, the iHS was calculated and then normalized . Per-SNP iHSs were smoothed in windows of 31 SNPs, centered on each SNP. The final score is -log10 of the proportion of smoothed scores higher than each SNP’s smoothed score. For our analysis, we used the Batch Coordinate Conversion (liftOver) utility (UCSC GB) to convert the genome coordinates from assembly NCBI36/hg18 to assembly GRCh37/hg19. We scanned the normalized iHSs across the whole genome and selected the genomic intervals where the iHS was ≥ 2. After these regions were identified, we extended their boundaries to the nearest loci where the iHS exactly vanished.
Selective sweep scan (S): the 5% lowest S scores
The S score is based on a comparison between Homo sapiens DNA and Neanderthal DNA . We downloaded the regions with S scores from the “5% Lowest S” track of the UCSC GB and denoted them as “5LSRs” (5% lowest S regions). Green et. al. identified polymorphic sites among five modern human genomes and determined the ancestral or derived state of each SNP. The states of the human alleles were used to estimate the expected number of derived alleles in Neanderthal in a 100000-base window around each SNP. The S scores were used to compare the observed number of Neanderthal alleles to the expected number in each window. A positive S score indicates more derived alleles in Neanderthal than expected given the frequency of derived alleles in human; a negative S score, on the other hand, indicates fewer derived alleles in Neanderthal, which might suggest positive selection in the human lineage after divergence from Neanderthal and before divergence in human populations. The 5LSRs represent the regions in the 5% lower percentile of the S score.
Conserved elements (CEs)
CEs are sequences in the genome that are conserved across species . Conserved regions have a reduced rate of evolution compared to the expected rate under neutral drift. The CEs used in this study were downloaded from the “Conservation (cons46way)” track of the UCSC GB. This track shows measurements of evolutionary conserved elements using two phylogenetic methods, phastCons and phyloP. The CEs used in this study were predicted using ten primates, Homo sapiens (reference species), Pan troglodytes, Gorilla gorilla, Pongo pygmaeus abelii, Macaca mulatta, Papio hamadryas, Callithrix jacchus, Tarsier syrichta, Microcebus murinus and Otolemur garnettii.
Gene expression data
The complete list of the 12559 genes expressed in lymphoblastoid cell lines LBEG (LymphoBlastoid Expressed Genes) and studied in  has been obtained from supplementary data (FileS2.xls) available on “Genetics” journal web site. From this list we selected the subset of 1888 genes differently expressed (DE genes) between H. sapiens and P. troglodytes using an FDR cut-off of 0.050.
We associated each group of CGIs (CE, HIR, 5SLR) with their putative target genes trough Genomic Regions Enrichment of Annotations Tool (GREAT) by using the default association rule  obtaining: 10867 CE CGIs genes, 1726 HIR CGIs genes and 275 5SLR CGIs genes. We then restricted each set of genes to the intersection with LBEG, retaining 7436 CE CGIs genes, 1190 HIR CGIs genes and 214 5SLR genes.
We tested the enrichment of DE in each set of genes associated with the considered signatures of selection by using Fisher’s Exact Test.
To be conservative, we have intersected the list of PS CGIs identified by considering each cell line and each particular histone modification. Hence, the previous GREAT analysis allows labelling as PS the associated genes. We have tested with a Fisher’s Exact Test the enrichment of PS genes inside the classes obtained by overlapping selective signatures and DE.
We used a hypergeometric-based approach to test the null hypothesis that the possible enrichment of H3K4me3, H3K27ac and H3K36me3 is independent of the presence of signals of natural selection. In particular we considered: k, the observed number of CGIs containing both PSs and signatures of selective pressure, as the number of success in the sample; n, the number of CGIs characterized by signatures of selective pressure only, as the sample size; M, the total number of CGIs with PS, as the number of successes in the population; and N, the total number of CGIs, as the population size (see Additional file 2: Table S2, Additional file 4: Table S3 and Additional file 23: Table S4). For statistical significance we set the threshold for the Bonferroni corrected p-value at 10-3. All the statistical analyses were performed with R ver. 2.14.2 (R Foundation for Statistical Computing, Vienna, Austria; http://www.r-project.org/).
Histone H3 lysine 4 trimethylation
Histone H3 lysine 27 acetylation
Histone H3 lysine 36 trimethylation
Transcription start sites
Integrated haplotype score
High iHS regions
- S score:
Selective sweep scan score
% lowest S regions
Biased gene conversion
Bird A, Taggart M, Frommer M, Miller OJ, Macleod D: A fraction of the mouse genome that is derived from islands of nonmethylated, CpG-rich DNA. Cell. 1985, 40 (1): 91-99. 10.1016/0092-8674(85)90312-5.
Illingworth RS, Bird AP: CpG islands-“a rough guide”. FEBS Lett. 2009, 583 (11): 1713-1720. 10.1016/j.febslet.2009.04.012.
Larsen F, Gundersen G, Lopez R, Prydz H: CpG islands as gene markers in the human genome. Genomics. 1992, 13 (4): 1095-1107. 10.1016/0888-7543(92)90024-M.
Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr ARW, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP: Orphan CpG islands identify numerous conserved promoters in the mammalian genome. PLoS Genet. 2010, 6 (9): e1001134-10.1371/journal.pgen.1001134.
Jones PA: Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012, 13 (7): 484-492. 10.1038/nrg3230.
Medvedeva YA, Fridman MV, Oparina NJ, Malko DB, Ermakova EO, Kulakovskiy IV, Heinzel A, Makeev VJ: Intergenic, gene terminal, and intragenic CpG islands in the human genome. BMC Genomics. 2010, 11: 48-10.1186/1471-2164-11-48.
Blackledge NP, Klose R: CpG island chromatin: a platform for gene regulation. Epigenetics: official journal of the DNA Methylation Society. 2011, 6 (2): 147-10.4161/epi.6.2.13640.
Li B, Carey M, Workman JL: The role of chromatin during transcription. Cell. 2007, 128 (4): 707-719. 10.1016/j.cell.2007.01.015.
Kouzarides T: Chromatin modifications and their function. Cell. 2007, 128 (4): 693-705. 10.1016/j.cell.2007.02.005.
Dou Y, Milne TA, Ruthenburg AJ, Lee S, Lee JW, Verdine GL, Allis CD, Roeder RG: Regulation of MLL1 H3K4 methyltransferase activity by its core components. Nat Struct Mol Biol. 2006, 13 (8): 713-719. 10.1038/nsmb1128.
Gregory GD, Vakoc CR, Rozovskaia T, Zheng X, Patel S, Nakamura T, Canaani E, Blobel GA: Mammalian ASH1L is a histone methyltransferase that occupies the transcribed region of active genes. Mol Cell Biol. 2007, 27 (24): 8466-10.1128/MCB.00993-07.
Guenther MG, Levine SS, Boyer LA, Jaenisch R, Young RA: A chromatin landmark and transcription initiation at most promoters in human cells. Cell. 2007, 130 (1): 77-88. 10.1016/j.cell.2007.05.042.
Balasubramanian D, Akhtar-Zaidi B, Song L, Bartels CF, Veigl M, Beard L, Myeroff L, Guda K, Lutterbaugh J, Willis J, et al: H3K4me3 inversely correlates with DNA methylation at a large class of non-CpG-island containing start sites. Genome Med. 2012, 4 (5): 47-10.1186/gm346.
Okitsu CY, Hsieh CL: DNA methylation dictates histone H3K4 methylation. Mol Cell Biol. 2007, 27 (7): 2746-10.1128/MCB.02291-06.
Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K: A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006, 125 (2): 315-326. 10.1016/j.cell.2006.02.041.
Ke XS, Qu Y, Rostad K, Li WC, Lin B, Halvorsen OJ, Haukaas SA, Jonassen I, Petersen K, Goldfinger N: Genome-wide profiling of histone h3 lysine 4 and lysine 27 trimethylation reveals an epigenetic signature in prostate carcinogenesis. PLoS One. 2009, 4 (3): e4687-10.1371/journal.pone.0004687.
Sandgren J, Andersson R, Rada-Iglesias A, Enroth S, Akerstrom G, Dumanski JP, Komorowski J, Westin G, Wadelius C: Integrative epigenomic and genomic analysis of malignant pheochromocytoma. Exp Mol Med. 2010, 42 (7): 484-10.3858/emm.2010.42.7.050.
Kaneda R, Takada S, Yamashita Y, Choi YL, Nonaka-Sarukawa M, Soda M, Misawa Y, Isomura T, Shimada K, Mano H: Genome-wide histone methylation profile for heart failure. Genes Cells. 2009, 14 (1): 69-77. 10.1111/j.1365-2443.2008.01252.x.
Wang Z, Zang C, Rosenfeld JA, Schones DE, Barski A, Cuddapah S, Cui K, Roh TY, Peng W, Zhang MQ, et al: Combinatorial patterns of histone acetylations and methylations in the human genome. Nat Genet. 2008, 40 (7): 897-903. 10.1038/ng.154.
Barski A, Cuddapah S, Cui K, Roh TY, Schones DE, Wang Z, Wei G, Chepelev I, Zhao K: High-resolution profiling of histone methylations in the human genome. Cell. 2007, 129 (4): 823-837. 10.1016/j.cell.2007.05.009.
Wagner EJ, Carpenter PB: Understanding the language of Lys36 methylation at histone H3. Nat Rev Mol Cell Biol. 2012, 13 (2): 115-126. 10.1038/nrm3274.
Britten RJ, Davidson EH: Gene regulation for higher cells: a theory. Science. 1969, 165 (891): 349-357.
Cain CE, Blekhman R, Marioni JC, Gilad Y: Gene expression differences among primates are associated with changes in a histone epigenetic modification. Genetics. 2011, 187 (4): 1225-1234. 10.1534/genetics.110.126177.
Shulha HP, Crisci JL, Reshetov D, Tushir JS, Cheung I, Bharadwaj R, Chou HJ, Houston IB, Peter CJ, Mitchell AC: Human-Specific Histone Methylation Signatures at Transcription Start Sites in Prefrontal Neurons. PLoS Biol. 2012, 10 (11): e1001427-10.1371/journal.pbio.1001427.
The-ENCODE-Project-Consortium: An integrated encyclopedia of DNA elements in the human genome. Nature. 2012, 489 (7414): 57-74. 10.1038/nature11247.
Cocozza S, Akhtar MM, Miele G, Monticelli A: CpG Islands Undermethylation in Human Genomic Regions under Selective Pressure. PLoS One. 2011, 6 (8): e23156-10.1371/journal.pone.0023156.
Celniker SE, Dillon LA, Gerstein MB, Gunsalus KC, Henikoff S, Karpen GH, Kellis M, Lai EC, Lieb JD, MacAlpine DM, et al: Unlocking the secrets of the genome. Nature. 2009, 459 (7249): 927-930. 10.1038/459927a.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, et al: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473 (7345): 43-49. 10.1038/nature09906.
Nielsen R: Molecular signatures of natural selection. Annu Rev Genet. 2005, 39: 197-218. 10.1146/annurev.genet.39.073003.112420.
Voight BF, Kudaravalli S, Wen X, Pritchard JK: A map of recent positive selection in the human genome. PLoS Biol. 2006, 4 (3): e72-10.1371/journal.pbio.0040072.
Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz MH, et al: A draft sequence of the Neandertal genome. Science. 2010, 328 (5979): 710-722. 10.1126/science.1188021.
Siepel A, Bejerano G, Pedersen JS, Hinrichs AS, Hou M, Rosenbloom K, Clawson H, Spieth J, Hillier LW, Richards S, et al: Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. Genome Res. 2005, 15 (8): 1034-1050. 10.1101/gr.3715005.
Cohen NM, Kenigsberg E, Tanay A: Primate CpG islands are maintained by heterogeneous evolutionary regimes involving minimal selection. Cell. 2011, 145 (5): 773-10.1016/j.cell.2011.04.024.
McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G: GREAT improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010, 28 (5): 495-501. 10.1038/nbt.1630.
Woo YH, Li WH: Evolutionary conservation of histone modifications in mammals. Mol Biol Evol. 2012, 29 (7): 1757-1767. 10.1093/molbev/mss022.
Thomson JP, Skene PJ, Selfridge J, Clouaire T, Guy J, Webb S, Kerr ARW, Deaton A, Andrews R, James KD: CpG islands influence chromatin structure via the CpG-binding protein Cfp1. Nature. 2010, 464 (7291): 1082-1086. 10.1038/nature08924.
Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, Degner JF, Gilad Y, Pritchard JK: DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol. 2011, 12 (1): R10-10.1186/gb-2011-12-1-r10.
Kolasinska-Zwierz P, Down T, Latorre I, Liu T, Liu XS, Ahringer J: Differential chromatin marking of introns and expressed exons by H3K36me3. Nat Genet. 2009, 41 (3): 376-381. 10.1038/ng.322.
Elango N, Soojin VY: Functional relevance of CpG island length for regulation of gene expression. Genetics. 2011, 187 (4): 1077-1083. 10.1534/genetics.110.126094.
Ruthenburg AJ, Allis CD, Wysocka J: Methylation of lysine 4 on histone H3: intricacy of writing and reading a single epigenetic mark. Mol Cell. 2007, 25 (1): 15-30. 10.1016/j.molcel.2006.12.014.
Zhou Y, Lu Y, Tian W: Epigenetic features are significantly associated with alternative splicing. BMC Genomics. 2012, 13 (1): 123-10.1186/1471-2164-13-123.
Luco RF, Allo M, Schor IE, Kornblihtt AR, Misteli T: Epigenetics in alternative pre-mRNA splicing. Cell. 2012, 144 (1): 16-26.
Barbosa-Morais NL, Irimia M, Pan Q, Xiong HY, Gueroussov S, Lee LJ, Slobodeniuc V, Kutter C, Watt S, Colak R, et al: The evolutionary landscape of alternative splicing in vertebrate species. Science. 2012, 338 (6114): 1587-1593. 10.1126/science.1230612.
Vernot B, Stergachis AB, Maurano MT, Vierstra J, Neph S, Thurman RE, Stamatoyannopoulos JA, Akey JM: Personal and population genomics of human regulatory variation. Genome Res. 2012, 22 (9): 1689-1697. 10.1101/gr.134890.111.
Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes Dev. 2011, 25 (10): 1010-10.1101/gad.2037511.
Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, Gabriel SB, Platko JV, Patterson NJ, McDonald GJ, et al: Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002, 419 (6909): 832-837. 10.1038/nature01140.
Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A: Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010, 20 (1): 110-121. 10.1101/gr.097857.109.
This work was supported by Institutional funds by University “Federico II”, Naples and POR Campania FSE 2007-2013, Project CREME. and partially supported by The Epigenomic Flagship Project- Epigen, CNR. MMA and GS are the recipients of a fellowship by Doctorate of Computational Biology and Bioinformatics, University “Federico II”, Naples.
The authors declare that they have no competing interests.
MMA, SC, AM and GM conceived the study and participated in its design and coordination. MMA, GS and SC performed the statistical analysis. MMA, GM, GS and AM helped to draft the manuscript. All authors read and approved the final manuscript.
Giovanni Scala contributed equally to this work.
Electronic supplementary material
Additional file 2: Table S2: Raw data used to calculate CGIs enriched with H3K4me3, H3K27ac and H3K36me3 in different cell lines. (XLS 46 KB)
Additional file 3: Euler diagram showing the overlaps among CGIs localized in the regions under selective pressure detected by the three methods used.(PDF 41 KB)
Additional file 4: Table S3: Raw data used to calculate CGIs enriched with H3K4me3, H3K27ac and H3K36me3 at 5’, intragenic, 3’ and intergenic locations in different cell lines. (XLS 112 KB)
Additional file 5: Enrichment of H3K4me3 modification in 5’ CpG islands under selective pressure. Black bars represent the fraction of CGIs containing histone modification marks within regions that show signatures of natural selection (HIRs, CEs and 5LSRs). Grey bars represent the fraction of CGIs containing histone modification marks within regions that do not show signatures of selective events. The X-axis indicates the analyzed cell lines. An asterisk (*) above a bar indicates a statistically non-significant difference. (PDF 565 KB)
Additional file 6: Same notation as Additional file 5. (PDF 544 KB)
Additional file 7: Same notation as Additional file 5. (PDF 559 KB)
Additional file 8: Same notation as Additional file 5. (PDF 571 KB)
Additional file 9: Same notation as Additional file 5. (PDF 215 KB)
Additional file 10: Same notation as Additional file 5. (PDF 521 KB)
Additional file 11: Same notation as Additional file 5. (PDF 516 KB)
Additional file 12: Same notation as Additional file 5. (PDF 522 KB)
Additional file 13: Same notation as Additional file 5. (PDF 525 KB)
Additional file 14: Same notation as Additional file 5. (PDF 548 KB)
Additional file 15: Same notation as Additional file 5. (PDF 527 KB)
Additional file 16: Same notation as Additional file 5. (PDF 510 KB)
Additional file 17: Same notation as Additional file 5. (PDF 605 KB)
Additional file 18: Same notation as Additional file 5. (PDF 533 KB)
Additional file 19: Same notation as Additional file 5. (PDF 588 KB)
Additional file 20: Same notation as Additional file 5. (PDF 520 KB)
Additional file 21: Same notation as Additional file 5. (PDF 576 KB)
Additional file 22: Same notation as Additional file 5. (PDF 559 KB)
Additional file 23: Table S4: Raw data used to calculate CGIs enriched with H3K4me3, H3K27ac and H3K36me3 in different cell lines according to the CGIs evolutionary model. (XLS 44 KB)
Additional file 24: Table S5: Lists number of PS genes that are present in each class defined by the possible presence of each signature of selection (“+” stands for presence, “-” stands for absence) and by the membership to DE set (“+” stands for belonging, “-” stands for not belonging). The exact Fisher’s test p-values are reported, highlighting the statistical significant ones. (XLS 9 KB)
Additional file 25: Table S6: Lists percentages of PS genes for classes defined by the possible presence of each signature of selection (“+” stands for presence, “-” stands for absence) and by the membership to DE set (“+” stands for belonging, “-” stands for not belonging). For each selection signature we give the p-value for the percentages of PS genes in case of DE + compared with DE - (p-values are reported on the same row), and in case of DE + considering the classes of genes under selection or not under selection (p-values are reported on the same column). The p-values are obtained by means of an exact Fisher’s test, highlighting the statistical significant values. (XLS 12 KB)
About this article
Cite this article
Akhtar, M.M., Scala, G., Cocozza, S. et al. CpG islands under selective pressure are enriched with H3K4me3, H3K27ac and H3K36me3 histone modifications. BMC Evol Biol 13, 145 (2013). https://doi.org/10.1186/1471-2148-13-145
- Selective pressure
- CpG islands