Genome-wide DNA methylation analyses in the brain reveal four differentially methylated regions between humans and non-human primates
© Wang et al.; licensee BioMed Central Ltd. 2012
Received: 9 January 2012
Accepted: 3 August 2012
Published: 16 August 2012
The highly improved cognitive function is the most significant change in human evolutionary history. Recently, several large-scale studies reported the evolutionary roles of DNA methylation; however, the role of DNA methylation on brain evolution is largely unknown.
To test if DNA methylation has contributed to the evolution of human brain, with the use of MeDIP-Chip and SEQUENOM MassARRAY, we conducted a genome-wide analysis to identify differentially methylated regions (DMRs) in the brain between humans and rhesus macaques. We first identified a total of 150 candidate DMRs by the MeDIP-Chip method, among which 4 DMRs were confirmed by the MassARRAY analysis. All 4 DMRs are within or close to the CpG islands, and a MIR3 repeat element was identified in one DMR, but no repeat sequence was observed in the other 3 DMRs. For the 4 DMR genes, their proteins tend to be conserved and two genes have neural related functions. Bisulfite sequencing and phylogenetic comparison among human, chimpanzee, rhesus macaque and rat suggested several regions of lineage specific DNA methylation, including a human specific hypomethylated region in the promoter of K6IRS2 gene.
Our study provides a new angle of studying human brain evolution and understanding the evolutionary role of DNA methylation in the central nervous system. The results suggest that the patterns of DNA methylation in the brain are in general similar between humans and non-human primates, and only a few DMRs were identified.
The most significant difference between humans and non-human primates is the highly developed cognitive ability in humans. Though many efforts have been made to delineate the genetic basis of the cognitive difference, it remains unclear. Because of the great similarity of genome sequences between humans and non-human primates [1, 2], more and more studies have started to look at gene regulation differences in the brain [3, 4], which are thought to be much more important for human brain evolution. As epigenetic modifications play important roles in learning and memory [5–7], it is possible that epigenetic modifications have contributed to human brain evolution. However, this has rarely been addressed.
Comparison of DNA methylation between humans and non-human primates at the orthologous regions in the genome could be an informative attempt to reveal the evolutionary roles of DNA methylation in human brain evolution. Two small-scale studies that compared DNA methylation of brain between humans and non-human primates have been conducted. It was reported that DNA methylation differences between human and chimpanzee were predominantly observed in the brain . Another study reported a species-specific methylated region in an Alu-Sg1 repeat . Recently, several large-scale studies of DNA methylation comparison between humans and non-human primates in multiple tissues (brain was not included) and cell lines have also been conducted, and the differentially methylated regions were identified mostly in repeats and CpG islands [10–12]. But the situation of DNA methylation in the brain, which is intuitively more important for primates, is largely unknown.
We aim to identify genes that show differential DNA methylation levels in the prefrontal cortex (PFC) between humans and non-human primates. We first conducted a genome-wide methylation scan using the MeDIP-Chip technology, and then the candidate regions were subject to the analysis using SEQUENOM MassARRAY  for confirmation and further validation using independent samples. Finally, we identified four regions/genes differentially methylated between humans and rhesus macaques in PFC. Two genes have neural related functions. The phylogenetic comparison among human, chimpanzee, rhesus macaque and rat suggested that a DMR in the promoter of K6IRS2 was hypomethylated specifically in the human brain, not the other species.
General profile of genome-wide DNA methylation
Information of samples
We then tested if the DNA methylation level was negatively correlated with gene expression using the published gene expression data . As the coverage of the expression data of rhesus macaques was not high enough, we only tested the correlation in humans. We used a 540 bp sliding window to access the correlations between gene expression and the methylation levels of different locations of the 2,700bp promoter regions. As shown in Additional file 1 (Figure S1), all the five sliding windows showed significant negative correlation. Interestingly, the region of -40bp-500bp showed a non-monotonic relationship, which is consistent with the recent finding that the gene-body methylation has a non-monotonic relationship with gene expression .
Identification of candidate regions showing differential DNA methylation
To identify the candidate differentially methylated regions between humans and rhesus macaques, we used two different methods to analyze the MeDIP-Chip data: the NimbleScan and the Batman methods . The regions that were estimated having differential DNA methylation levels by both methods were selected as the candidate regions for further analysis. Around 42% of the peaks identified by NimbleScan (i.e. the peaks only observed in either of the two species) were confirmed by Batman (Additional file 2). As a preliminary screening step, we used a nominal P value cutoff of 0.05 to select the candidate regions. Totally, we obtained 150 candidate DMRs (Additional file 3), and most of them were located (or partly located) in the CpG islands (142 for humans and 141 for rhesus macaques) based on the UCSC genome annotation (http://genome.ucsc.edu/). This is consistent with the features of the MeDIP-Chip technology, which tends to find CpG methylation in the high CpG density regions . Gene ontology (GO) analysis for the 150 candidate DMRs was conducted using DAVID, however, no enrichment was observed.
It has been suggested that there is a high false discovery rate for the MeDIP-Chip method . And even the DMRs identified between two human embryonic stem cell lines using the more advanced MeDIP-seq technique do not have good validation . Hence, the 150 candidate DMRs need to be further tested by independent methods.
Validation of the candidate DMRs using MassARRAY
For the 150 candidate regions likely having differential DNA methylation levels between humans and rhesus macaques, we conducted further analysis in the same samples (3 human versus 3 rhesus samples) using the Sequenom MassARRAY method , which has much higher resolution and accuracy for DNA methylation compared to the MeDIP-Chip method and is suitable for candidate gene validation. Among the 150 candidate DMRs, a total of 118 DMRs were successfully amplified and obtained eligible data for more than one sample in both species. It should be noted that for the MassARRAY method, the minimum measuring unit was CpG unit, which are one or several consecutive CpG sites. A total of 948 eligible CpG units were obtained for the 118 DMRs, and 688 of them can be paired based on their orthology between humans and rhesus macaques. The other 260 eligible CpG units could not be paired, therefore, were excluded from further analyses.
The four DMRs (differentially methylated regions) validated by the MassARRAY analysis
Human PCR regionb
Rhesus PCR region
Validation of DMRs in additional independent samples
Corrected P valuea
P value (all)b
Differentially methylated CpG units between humans and rhesus macaques
Corrected P valuef
NA 4 5
13 12 NA
14 15 16
NA 1 NA
7 8 9 10
7 8 9 10
17 18 19
16 17 18
9 NA 10 11
8 9 10 NA
10 11 12 NA
10 11 NA 12
5 NA 6 NA
4 5 NA 6
3 4 5
NA 5 NA
18 19 NA
18 19 20
15 16 NA 17 18 19
8 9 10 NA 11 NA
45 46 47
43 NA 44
NA 3 NA
4 NA 5
Features and functional implication of the DMRs
For the 4 DMRs, DMR143 is located in the promoter of RNF32; DMR103 covers the promoter and the first exon of K6iRS2; DMR8 covers the first exon and extends to the first intron of ICAM1, and DMR13 is located in the first exon of ProSAPiP1. Based on the UCSC human genome annotation, DMR143 is entirely located in the CpG island; DMR103 and DMR8 have more than half of their full length located in the CpG islands, and DMR13 is located 2,992bp from the nearest CpG island.
We then checked if there were any repeat elements in the DMRs. We found nearly all the 4 DMRs had no repeat element, because we removed the non-unique probes of MeDIP-Chip array, which would remove most of probes targeted to repeat elements. However, there was an exception in rhesus macaques where we found a MIR3 repeat element that was entirely located in DMR8 (ICAM1 gene) based on the UCSC genome annotation.
To look into the potential functional significance of the differentially methylated regions and its implications for human evolution, we conducted functional annotation for the four genes located in the 4DMRs. Interestingly, two genes are involved in neural functions. ICAM1 (intercellular adhesion molecule 1) is a transmembrane glycoprotein, which belongs to the immunoglobulin superfamily, and has been reported associated with schizophrenia  and Alzheimer’s disease . ProSAPiP1 is a synaptic protein that interacts with the PDZ domain of ProSAP2/Shank3. For the other two genes, K6IRS2 encodes a type II keratin that is expressed in the inner root sheath of hair follicles , and the function of K6IRS2 is not known. RNF32 is a ring finger protein, and its function is unknown.
Next, we addressed the functional conservation of the 4 DMRs genes, and we studied their protein conservation using human versus rhesus macaque dN/dS ratios (ratio of nonsynonymous substitution rate vs. synonymous substitution rate). The dN/dS ratios of all 4 genes were smaller than one, suggesting functional conservation (0.14 for RNF32, 0.16 for K6IRS2, 0.13 for ProSAPiP1 and 0.49 for ICAM1). These results suggest that the identified DMRs are not likely caused by protein functional divergence between different species.
DNA methylation comparison in multiple species
In this study, we employed a two-step method to identify the differentially methylated regions and 4 regions were found to have different methylation levels in PFC between humans and rhesus macaques, which are informative candidates for further study of human brain evolution. However, there are also methodological limitations. The 150 candidate regions are all high CpG content regions because the NimbleGen method tends to find highly methylated regions with high CpG content. It was reported that high CpG content regions were of higher quality in MeDIP-Chip data . In studies of tissue specific DNA methylation regions, high CpG content regions such as the CpG islands showed less between individual variance and less differences between different tissues [28–30], and this may be the same case between different species. Though we cannot exclude the possibility that there are differentially methylated regions in low CpG content regions, the high CpG content regions, on which we focused in this study, are more likely to be functionally important [14, 15] and more likely to have definite between-species differences and less within-species variations, thus evolutionarily more significant. In addition, considering the multiple data filtering steps using different methods, the number of identified DMRs in this study was likely an underestimation of the between-species methylation divergence in the brain.
Repeat elements are important in DNA methylation. Farcas et al. reported a small region of ALU-Sg1 element that was differentially methylated between human and chimpanzee cerebral cortex. In this study, we found another case of SINE element differentially methylated between humans and rhesus macaques, implying the roles of repeat elements in the epigenetic evolution of primates. In the rhesus macaque genome, the MIR3 element located in DMR8 has a length of 98 bp and divergence of 34.4% to its consensus sequence. In DMR8, the human ortholog of this region have 6 substitutions and 1 deletion compared with rhesus macaque, and it is not annotated as repeat region by UCSC. MIR3 belongs to the MIR (Mammalian-wide interspersed repeats) family, which is one of the oldest tRNA-derived SINE (short interspersed element) elements, and they integrated into the host genomes before the radiation of mammals . In DMR8, three CpG sites (CpG12, CpG13 and CpG14) of rhesus macaques are located in MIR3. Interestingly, the three CpG sites are the only CpG sites that showed significant methylation differences between humans and rhesus macaques (Figure 3 and Table 4), suggesting that the MIR3 element has led to the differential methylation of DMR8 between humans and rhesus macaques.
Among the 4 validated DMRs, only one is located in the repeat region, suggesting that the differentially methylated regions between humans and nonhuman primates are not necessarily restricted to repeat regions and can be high CpG content regions, which have been previously thought to be less variable [28–30].
As DNA methylation in the promoter regions will repress gene expression, it would be informative to see whether the differential DNA methylations between humans and rhesus macaques would lead to expression differences. We analyzed the published gene expression data from human and rhesus macaque cerebral cortices . Among the four DMRs identified in our study, two DMR related genes (ICAM1andProSAPiP1) had eligible expression data and showed significant gene expression differences between human and rhesus macaque (P = 0.000757 and 0.00764 for ICAM1 and ProSAPiP1 respectively). However, for both genes, the species having higher DNA methylation levels also showed higher gene expression levels, which is inconsistent with the expected repression of gene expression by a higher level of DNA methylation in the promoter. However, these two DMRs also covered the first exons of the corresponding genes, and the gene-body methylation was reported to be positively correlated with gene expression in human cells . Hence, detailed analysis in the future is needed to reveal the influence of DMRs on gene expression.
It is difficult to address how the observed methylation differences have been formed during primate evolution. Trans-generational inheritance is crucial if epigenetic modification should play a role in evolution. Many studies have reported that epigenetic modifications can be inherited across generations [33–35]. But the inheritance of DNA methylation is not as stable as that of DNA sequence and the mechanisms underlying the trans-generational inheritance of epigenetic modification are not well understood.
Recently, it was reported that DNA methylation is correlated with DNA sequence variations [36–38] and non-coding RNAs [39–42], and these correlations might act as an indirect mechanism that explains the trans-generational inheritance of epigenetic modification. The DNA sequence substitutions during evolution may cause the emergence of new CpG sites and/or loss of the existing CpG sites, eventually leading to methylation divergence of specific genomic regions. It has been reported that the CpG-SNPs in the human genome have contributions to allele specific DNA methylation , implying that even at the population level, the methylation divergence can occur due to DNA sequence polymorphisms. Alternatively, DNA sequence differences could also result from DNA methylation differences between species [10–12].
Another possibility is that epigenetic modifications altered by environment can be transmitted to the next generation directly. Global epigenetic reprogramming including demethylation of DNA occurs in the mammalian primordial germ cells and in early embryos . Although global epigenetic reprogramming will restrict the trans-generational epigenetic inheritance, the erasure of DNA methylation modifications is not absolute [44, 45], suggesting the possibility of direct transmitting of DNA methylation modifications to the next generation. Additionally, at the whole genome level, the over-all similar DNA methylation patterns between humans and macaques cannot be explained by the similarity of living environment because they inhabit totally different environments, suggesting potentially vertical inheritance of methylation.
A lot of candidate DMRs obtained from the MeDIP-Chip analysis were not validated by MassARRAY and bisulfite sequencing. This discrepancy could be explained by two possible sources. First, for MeDIP-Chip analysis, we used different arrays for humans and rhesus macaques. Even though the validation for each species-specific array is good, the technical bias may still exist, as the array for rhesus macaques has not been rigorously tested. Second, besides detecting the genome-wide DNA methylation patterns, we focused on identifying differentially methylated regions between humans and non-human primates. Thus, a relaxed cut-off (nominal P smaller than 0.05 without multiple test correction) was applied for selecting the 150 candidates from the MeDIP-Chip data, likely resulting in a relatively high number of false-positives.
Our study is the first large-scale comparison of DNA methylation between human and non-human primates in the brain. We found several regions of lineage specific DNA methylation, indicating that DNA methylation plays potential roles in brain evolution, and the regions we reported are good candidates for further functional studies, and will shed light on the study of human brain evolution. The results suggest that the patterns of DNA methylation in the brain are in general similar between humans and non-human primates.
Frozen tissues were obtained from PFC of 5 humans, 7 rhesus macaques, two chimpanzees and two rats with no known neuronal diseases or drug abuse. Five of the 7 rhesus macaques and all 5 human samples were sex and age matched, following the reported criteria that the ages of humans are in general 3 times of the ages of rhesus macaques . These age-matched human and macaque samples (3 from each species) were used in the initial genome-wide scanning using the MeDIP-Chip method. The information of the samples is shown in Table 1. For the human subjects, informed consents were obtained from their relatives. The use of animal tissues follows internationally recognized guidelines. The research protocol was approved by the internal review board of Kunming Institute of Zoology, Chinese Academy of Sciences (reference number: SWYX-2010-002).
MeDIP-Chip (Methylated DNA immunoprecipitation)
Genomic DNA from PFC was extracted using DNeasy Blood & Tissue Kit (Qiagen) and fragmented by ultrasonic method following the manufacturer’s instructions. Immunoprecipitation of methylated DNA was performed using anti-5-methyl cytidine (mouse) and BiomagTM magnetic beads coupled anti-mouse IgG. Immunoprecipitated DNA was eluted and purified by phenol chloroform extraction. Then the Input and IP DNAs were labeled with Cy5- and Cy3-labeled random 9-mers respectively. The labeled human and rhesus macaque DNAs (Input and IP) were hybridized to NimbleGen HG18 RefSeq promoter arrays and the customer designed rhesus macaque DNA methylation array respectively. The scanning was performed with the Axon GenePix 4000B microarray scanner.
The NimbleGen HG18 RefSeq promoter array is a single array design containing all known well-characterized 18,028 RefSeq promoter regions (from about -2200bp to +500bp of the transcription start sites) totally covered by about 385,000 probes. The customer-designed rhesus macaque array was designed according to the commercial NimbleGen HG18 RefSeq promoter array. The regions covered on HG18 RefSeq promoter array were aligned to the rhesus macaque genome using BLAT software, and the alignments that have more than two hits were removed. Then the aligned rhesus macaque regions were merged if they overlapped, regions shorter than 1,700bp were removed because these regions were likely the false orthologous regions or those with large sequence divergence between the two species. The regions shorter than its corresponding human regions were extend equally on both sides to make the regions of both species equal length. Then the probes of rhesus macaque array were designed and synthesized according to these genomic regions. In all the aligned regions, only the Ensembl annotated one to one orthologous gene pairs were subjected to the following DNA methylation and comparison between humans and rhesus macaques.
MeDIP-Chip data analyses
The raw data were extracted by NimbleScan software. In the analyses of the genome-wide profiling of MeDIP-Chip data, the methylation levels of every sliding window of 540bp were calculated using Batman software . Before the analysis, the quantile normalization has been conducted using Limma package  implemented in R, and the probes that are not unique in the respective genomes were removed from the analysis.
Two different methods were used to obtain the list of candidate DMRs from the MeDIP-Chip data in order to reduce the false discovery rate. In the first method, we used the NimbleScan software developed by NimbleGen under its standard procedure and instruction to find the candidate DMRs. The regions were defined to be candidate DMRs only if they were reported methylated in all three samples of one species and none in the other species. Then, the quantitative methylation levels were calculated using Batman and compared between the two species for the regions that were reported to show different methylation status by NimbleScan method, and the differentially methylated regions showed large differences of CpG numbers between the two species were removed as they were likely due to the technical bias of NimbleScan method, which did not use CpG density to do the adjustment. Custom Perl scripts were used to accomplish the above procedures, and the statistics were conducted using R.
The MeDIP-Chip data for the human and rhesus macaque PFC samples can be downloaded from NCBI’s GEO (http:/http://www.ncbi.nlm.nih.gov/geo/) database (accession number: GSE27461).
We used the EpiTect Bisulfite Kits (Qiagen) to conduct the bisulfite conversion of DNA under the manufacturer's instructions, and the PCR primers were designed using Methyl Primer Express (v1.0, ABi Corp.). At least 10 clones were sequenced for each fragment by ABI 3130 sequencer after PCR, ligation and cloning. The clone sequences were then analyzed using BiQ Analyzer software  to get the methylation levels after manual checking.
SEQUENOM MassARRAY experiments and data analyses
Experiments were conducted under manufacturer’s instructions. Because of the difficulty of designing perfect eligible PCR primers with no CpG, for some regions, we allowed one CpG in the primers, and in rare cases, two CpGs were also allowed, and the CpGs in the primers had been removed from data analyses. As for the 4 DMRs found by MassARRAY method, none of them had CpG in its primers. For some regions that no eligible PCR primer could be designed in at least one species, we tested the proximal regions of both species instead.
Any unreliable CpG units (one or more CpG sites that measured as one unit in SEQUENOM MassARRAY technique) reported by EpiTYPER software (SEQUENOM) were removed from our analyses. To compare the methylation levels between humans and rhesus macaques, we used Clustal W software to align the bisulfite-converted sequences of the corresponding human and rhesus macaque regions. Then the orthologous pairs of CpG units of these two species were determined and compared respectively. The orthologous CpG unit pair was defined as the CpG unit pair that had at least one orthologous CpG site and all other CpG sites that did not have orthologous CpG site were due to the lack of CpG site on the counterpart species. Therefore, in some occasions, multiple short CpG units need to be merged in order to pair the longer CpG unit in the counterpart species. To calculate the average DNA methylation level of each region, we used only the CpG units that had corresponding orthologous CpG units in the counterpart species, but the methylation levels of the CpG units without orthologous CpG units were also shown (Figure 3).
The dN/dS values of the 4 DMR genes were obtained from the Ensembl database (http://www.ensembl.org/).
Gene expression analyses
To analyze the relationship between gene expression and DNA methylation and compare gene expression differences between human and rhesus macaque, we used the published gene expression data , and the human and rhesus macaque cerebral cortex samples with similar ages compared to our samples were used (NCBI’s GEO accession number: GSM289929, GSM289938, GSM289944, GSM289947, GSM289948, GSM289950, GSM289951, GSM289972, GSM289973, GSM289974, GSM289975, GSM289976 and GSM289977).
We thank Minati Singh for revising the manuscript, and Hui Zhang and Yanjiao Li for their technical assistance in the experiments. This study was supported by grants from the National 973 project of China (grant numbers, 2011CBA00401), the National Natural Science Foundation of China (grant numbers, 31130051), and the Natural Science Foundation of Yunnan Province of China (grant numbers, Y003721111).
- Varki A, Altheide TK: Comparing the human and chimpanzee genomes: searching for needles in a haystack. Genome Res. 2005, 15: 1746-1758. 10.1101/gr.3737405.PubMedView ArticleGoogle Scholar
- Gibbs RA, Rogers J, Katze MG, Bumgarner R, Weinstock GM, Mardis ER, Remington KA, Strausberg RL, Venter JC, Wilson RK, Batzer MA, Bustamante CD, Eichler EE, Hahn MW, Hardison RC, Makova KD, Miller W, Milosavljevic A, Palermo RE, Siepel A, Sikela JM, Attaway T, Bell S, Bernard KE, Buhay CJ, Chandrabose MN, Dao M, Davis C, Delehaunty KD, Ding Y, Dinh HH, Dugan-Rocha S, Fulton LA, Gabisi RA, Garner TT, Godfrey J, Hawes AC, Hernandez J, Hines S, Holder M, Hume J, Jhangiani SN, Joshi V, Khan ZM, Kirkness EF, Cree A, Fowler RG, Lee S, Lewis LR, Li Z, Liu YS, Moore SM, Muzny D, Nazareth LV, Ngo DN, Okwuonu GO, Pai G, Parker D, Paul HA, Pfannkoch C, Pohl CS, Gibbs RA, Rogers J, Katze MG, Bumgarner R, Weinstock GM, Mardis ER, Remington KA, Strausberg RL, Venter JC, Wilson RK, Batzer MA, Bustamante CD, Eichler EE, Hahn MW, Hardison RC, Makova KD, Miller W, Milosavljevic A, Palermo RE, Siepel A, Sikela JM, Attaway T, Bell S, Bernard KE, Buhay CJ, Chandrabose MN, Dao M, Davis C, Delehaunty KD, Ding Y, Dinh HH, Dugan-Rocha S, Fulton LA, Gabisi RA, Garner TT, Godfrey J, Hawes AC, Hernandez J, Hines S, Holder M, Hume J, Jhangiani SN, Joshi V, Khan ZM, Kirkness EF, Cree A, Fowler RG, Lee S, Lewis LR, Li Z, Liu YS, Moore SM, Muzny D, Nazareth LV, Ngo DN, Okwuonu GO, Pai G, Parker D, Paul HA, Pfannkoch C, Pohl CS, Rogers YH, Ruiz SJ, Sabo A, Santibanez J, Schneider BW, Smith SM, Sodergren E, Svatek AF, Utterback TR, Vattathil S, Warren W, White CS, Chinwalla AT, Feng Y, Halpern AL, Hillier LW, Huang X, Minx P, Nelson JO, Pepin KH, Qin X, Sutton GG, Venter E, Walenz BP, Wallis JW, Worley KC, Yang SP, Jones SM, Marra MA, Rocchi M, Schein JE, Baertsch R, Clarke L, Csuros M, Glasscock J, Harris RA, Havlak P, Jackson AR, Jiang H, Liu Y, Messina DN, Shen Y, Song HX, Wylie T, Zhang L, Birney E, Han K, Konkel MK, Lee J, Smit AF, Ullmer B, Wang H, Xing J, Burhans R, Cheng Z, Karro JE, Ma J, Raney B, She X, Cox MJ, Demuth JP, Dumas LJ, Han SG, Hopkins J, Karimpour-Fard A, Kim YH, Pollack JR, Vinar T, Addo-Quaye C, Degenhardt J, Denby A, Hubisz MJ, Indap A, Kosiol C, Lahn BT, Lawson HA, Marklein A, Nielsen R, Vallender EJ, Clark AG, Ferguson B, Hernandez RD, Hirani K, Kehrer-Sawatzki H, Kolb J, Patil S, Pu LL, Ren Y, Smith DG, Wheeler DA, Schenck I, Ball EV, Chen R, Cooper DN, Giardine B, Hsu F, Kent WJ, Lesk A, Nelson DL, O'Brien WE, Prufer K, Stenson PD, Wallace JC, Ke H, Liu XM, Wang P, Xiang AP, Yang F, Barber GP, Haussler D, Karolchik D, Kern AD, Kuhn RM, Smith KE, Zwieg AS: Evolutionary and biomedical insights from the rhesus macaque genome. Science. 2007, 316: 222-234.PubMedView ArticleGoogle Scholar
- Khaitovich P, Hellmann I, Enard W, Nowick K, Leinweber M, Franz H, Weiss G, Lachmann M, Paabo S: Parallel patterns of evolution in the genomes and transcriptomes of humans and chimpanzees. Science. 2005, 309: 1850-1854. 10.1126/science.1108296.PubMedView ArticleGoogle Scholar
- Caceres M, Lachuer J, Zapala MA, Redmond JC, Kudo L, Geschwind DH, Lockhart DJ, Preuss TM, Barlow C: Elevated gene expression levels distinguish human from non-human primate brains. Proc Natl Acad Sci USA. 2003, 100: 13030-13035. 10.1073/pnas.2135499100.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsankova N, Renthal W, Kumar A, Nestler EJ: Epigenetic regulation in psychiatric disorders. Nat Rev Neurosci. 2007, 8: 355-367.PubMedView ArticleGoogle Scholar
- Graff J, Mansuy IM: Epigenetic codes in cognition and behaviour. Behav Brain Res. 2008, 192: 70-87. 10.1016/j.bbr.2008.01.021.PubMedView ArticleGoogle Scholar
- Levenson JM, Sweatt JD: Epigenetic mechanisms in memory formation. Nat Rev Neurosci. 2005, 6: 108-118. 10.1038/nrn1604.PubMedView ArticleGoogle Scholar
- Enard W, Fassbender A, Model F, Adorjan P, Paabo S, Olek A: Differences in DNA methylation patterns between humans and chimpanzees. Curr Biol. 2004, 14: R148-R149.PubMedView ArticleGoogle Scholar
- Farcas R, Schneider E, Frauenknecht K, Kondova I, Bontrop R, Bohl J, Navarro B, Metzler M, Zischler H, Zechner U, Daser A, Haaf T: Differences in DNA methylation patterns and expression of the CCRK gene in human and nonhuman primate cortices. Mol Biol Evol. 2009, 26: 1379-1389. 10.1093/molbev/msp046.PubMedView ArticleGoogle Scholar
- Pai AA, Bell JT, Marioni JC, Pritchard JK, Gilad Y: A Genome-Wide Study of DNA Methylation Patterns and Gene Expression Levels in Multiple Human and Chimpanzee Tissues. PLoS Genet. 2011, 7: e1001316-10.1371/journal.pgen.1001316.PubMedPubMed CentralView ArticleGoogle Scholar
- Molaro A, Hodges E, Fang F, Song Q, McCombie WR, Hannon GJ, Smith AD: Sperm methylation profiles reveal features of epigenetic inheritance and evolution in primates. Cell. 2011, 146: 1029-1041. 10.1016/j.cell.2011.08.016.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin D, Singer M, Dhahbi J, Mao G, Zhang L, Schroth G: Pachter L. 2011, Phyloepigenomic comparison of great apes reveals a correlation between somatic and germline methylation states. Genome Res, Boffelli DGoogle Scholar
- Ehrich M, Nelson MR, Stanssens P, Zabeau M, Liloglou T, Xinarianos G, Cantor CR, Field JK, van den Boom D: Quantitative high-throughput analysis of DNA methylation patterns by base-specific cleavage and mass spectrometry. Proc Natl Acad Sci USA. 2005, 102: 15785-15790. 10.1073/pnas.0507816102.PubMedPubMed CentralView ArticleGoogle Scholar
- Fazzari MJ, Greally JM: Epigenomics: beyond CpG islands. Nat Rev Genet. 2004, 5: 446-455. 10.1038/nrg1349.PubMedView ArticleGoogle Scholar
- Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Graf S, Tomazou EM, Backdahl L, Johnson N, Herberth M, Howe KL, Jackson DK, Miretti MM, Fiegler H, Marioni JC, Birney E, Hubbard TJ, Carter NP, Tavare S, Beck S: An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Res. 2008, 18: 1518-1529. 10.1101/gr.077479.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Somel M, Franz H, Yan Z, Lorenc A, Guo S, Giger T, Kelso J, Nickel B, Dannemann M, Bahn S, Webster MJ, Weickert CS, Lachmann M, Paabo S, Khaitovich P: Transcriptional neoteny in the human brain. Proc Natl Acad Sci USA. 2009, 106: 5743-5748. 10.1073/pnas.0900544106.PubMedPubMed CentralView ArticleGoogle Scholar
- Jjingo D, Conley AB, Yi S, Lunyak VV, Jordan KI: On the presence and role of human gene-body DNA methylation. Oncotarget. 2012, 3: 462-474.PubMedPubMed CentralView ArticleGoogle Scholar
- Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, Graf S, Johnson N, Herrero J, Tomazou EM, Thorne NP, Backdahl L, Herberth M, Howe KL, Jackson DK, Miretti MM, Marioni JC, Birney E, Hubbard TJ, Durbin R, Tavare S, Beck S: A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis. Nat Biotechnol. 2008, 26: 779-785. 10.1038/nbt1414.PubMedPubMed CentralView ArticleGoogle Scholar
- Laird PW: Principles and challenges of genome-wide DNA methylation analysis. Nat Rev Genet. 2010, 11: 191-203.PubMedView ArticleGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37: 1-13. 10.1093/nar/gkn923.PubMedView ArticleGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Carvalho B, Wu H, Brandenburg SA, Jeddeloh JA, Wen B, Feinberg AP: Comprehensive high-throughput arrays for relative methylation (CHARM). Genome Res. 2008, 18: 780-790. 10.1101/gr.7301508.PubMedPubMed CentralView ArticleGoogle Scholar
- Bock C, Tomazou EM, Brinkman AB, Muller F, Simmer F, Gu H, Jager N, Gnirke A, Stunnenberg HG, Meissner A: Quantitative comparison of genome-wide DNA methylation mapping technologies. Nat Biotechnol. 2010, 28: 1106-1114. 10.1038/nbt.1681.PubMedPubMed CentralView ArticleGoogle Scholar
- Ehrich M, Turner J, Gibbs P, Lipton L, Giovanneti M, Cantor C, van den Boom D: Cytosine methylation profiling of cancer cell lines. Proc Natl Acad Sci USA. 2008, 105: 4844-4849. 10.1073/pnas.0712251105.PubMedPubMed CentralView ArticleGoogle Scholar
- Kronig H, Riedel M, Schwarz MJ, Strassnig M, Moller HJ, Ackenheil M, Muller N: ICAM G241A polymorphism and soluble ICAM-1 serum levels: evidence for an active immune process in schizophrenia. Neuroimmunomodulation. 2005, 12: 54-59. 10.1159/000082364.PubMedView ArticleGoogle Scholar
- Pola R, Flex A, Gaetani E, Santoliquido A, Serricchio M, Pola P, Bernabei R: Intercellular adhesion molecule-1 K469E gene polymorphism and Alzheimer's disease. Neurobiol Aging. 2003, 24: 385-387. 10.1016/S0197-4580(02)00087-8.PubMedView ArticleGoogle Scholar
- Wendholt D, Spilker C, Schmitt A, Dolnik A, Smalla KH, Proepper C, Bockmann J, Sobue K, Gundelfinger ED, Kreutz MR, Boeckers TM: ProSAP-interacting protein 1 (ProSAPiP1), a novel protein of the postsynaptic density that links the spine-associated Rap-Gap (SPAR) to the scaffolding protein ProSAP2/Shank3. J Biol Chem. 2006, 281: 13805-13816. 10.1074/jbc.M601101200.PubMedView ArticleGoogle Scholar
- Langbein L, Rogers MA, Praetzel S, Winter H, Schweizer J: K6irs1, K6irs2, K6irs3, and K6irs4 represent the inner-root-sheath-specific type II epithelial keratins of the human hair follicle. J Invest Dermatol. 2003, 120: 512-522. 10.1046/j.1523-1747.2003.12087.x.PubMedView ArticleGoogle Scholar
- Bock C, Walter J, Paulsen M, Lengauer T: Inter-individual variation of DNA methylation and its implications for large-scale epigenome mapping. Nucleic Acids Res. 2008, 36: e55-10.1093/nar/gkn122.PubMedPubMed CentralView ArticleGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.PubMedPubMed CentralView ArticleGoogle Scholar
- Yagi S, Hirabayashi K, Sato S, Li W, Takahashi Y, Hirakawa T, Wu G, Hattori N, Ohgane J, Tanaka S, Liu XS, Shiota K: DNA methylation profile of tissue-dependent and differentially methylated regions (T-DMRs) in mouse promoter regions demonstrating tissue-specific gene expression. Genome Res. 2008, 18: 1969-1978. 10.1101/gr.074070.107.PubMedPubMed CentralView ArticleGoogle Scholar
- Jurka J, Zietkiewicz E, Labuda D: Ubiquitous mammalian-wide interspersed repeats (MIRs) are molecular fossils from the mesozoic era. Nucleic Acids Res. 1995, 23: 170-175. 10.1093/nar/23.1.170.PubMedPubMed CentralView ArticleGoogle Scholar
- Ball MP, Li JB, Gao Y, Lee JH, LeProust EM, Park IH, Xie B, Daley GQ, Church GM: Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat Biotechnol. 2009, 27: 361-368. 10.1038/nbt.1533.PubMedPubMed CentralView ArticleGoogle Scholar
- Jirtle RL, Skinner MK: Environmental epigenomics and disease susceptibility. Nat Rev Genet. 2007, 8: 253-262. 10.1038/nrg2045.PubMedView ArticleGoogle Scholar
- Rakyan VK, Beck S: Epigenetic variation and inheritance in mammals. Curr Opin Genet Dev. 2006, 16: 573-577. 10.1016/j.gde.2006.09.002.PubMedView ArticleGoogle Scholar
- Richards EJ: Inherited epigenetic variation–revisiting soft inheritance. Nat Rev Genet. 2006, 7: 395-401.PubMedView ArticleGoogle Scholar
- Kerkel K, Spadola A, Yuan E, Kosek J, Jiang L, Hod E, Li K, Murty VV, Schupf N, Vilain E, Morris M, Haghighi F, Tycko B: Genomic surveys by methylation-sensitive SNP analysis identify sequence-dependent allele-specific DNA methylation. Nat Genet. 2008, 40: 904-908. 10.1038/ng.174.PubMedView ArticleGoogle Scholar
- Schalkwyk LC, Meaburn EL, Smith R, Dempster EL, Jeffries AR, Davies MN, Plomin R, Mill J: Allelic skewing of DNA methylation is widespread across the genome. Am J Hum Genet. 2010, 86: 196-212. 10.1016/j.ajhg.2010.01.014.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang D, Cheng L, Badner JA, Chen C, Chen Q, Luo W, Craig DW, Redman M, Gershon ES, Liu C: Genetic control of individual differences in gene-specific methylation in human brain. Am J Hum Genet. 2010, 86: 411-419. 10.1016/j.ajhg.2010.02.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Mathieu O, Bender J: RNA-directed DNA methylation. J Cell Sci. 2004, 117: 4881-4888.PubMedGoogle Scholar
- Zaratiegui M, Irvine DV, Martienssen RA: Noncoding RNAs and gene silencing. Cell. 2007, 128: 763-776. 10.1016/j.cell.2007.02.016.PubMedView ArticleGoogle Scholar
- Mattick JS, Amaral PP, Dinger ME, Mercer TR, Mehler MF: RNA regulation of epigenetic processes. Bioessays. 2009, 31: 51-59. 10.1002/bies.080099.PubMedView ArticleGoogle Scholar
- Ponting CP, Oliver PL, Reik W: Evolution and functions of long noncoding RNAs. Cell. 2009, 136: 629-641. 10.1016/j.cell.2009.02.006.PubMedView ArticleGoogle Scholar
- Shoemaker R, Deng J, Wang W, Zhang K: Allele-specific methylation is prevalent and is contributed by CpG-SNPs in the human genome. Genome Res. 2010, 20: 883-889. 10.1101/gr.104695.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Lane N, Dean W, Erhardt S, Hajkova P, Surani A, Walter J, Reik W: Resistance of IAPs to methylation reprogramming may provide a mechanism for epigenetic inheritance in the mouse. Genesis. 2003, 35: 88-93. 10.1002/gene.10168.PubMedView ArticleGoogle Scholar
- Silva AJ, White R: Inheritance of allelic blueprints for methylation patterns. Cell. 1988, 54: 145-152. 10.1016/0092-8674(88)90546-6.PubMedView ArticleGoogle Scholar
- Ottinger MA, Maltison JA, Zelinski MB, Wu JM, Kohama SG, Roth GS, Lane MA, Ingram DK: The rhesus macaque as a model of human aging and age-related disease. Handbook of models for human aging. Volume 1. Edited by: Conn PM. 2006, Elservier Inc, Burlinglon, 457-464.View ArticleGoogle Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- Bock C, Reither S, Mikeska T, Paulsen M, Walter J, Lengauer T: BiQ Analyzer: visualization and quality control for DNA methylation data from bisulfite sequencing. Bioinformatics. 2005, 21: 4067-4068. 10.1093/bioinformatics/bti652.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.