- Open Access
Human-specific epigenetic variation in the immunological Leukotriene B4 Receptor (LTB4R/BLT1) implicated in common inflammatory diseases
Genome Medicine volume 6, Article number: 19 (2014)
Common human diseases are caused by the complex interplay of genetic susceptibility as well as environmental factors. Due to the environment’s influence on the epigenome, and therefore genome function, as well as conversely the genome’s facilitative effect on the epigenome, analysis of this level of regulation may increase our knowledge of disease pathogenesis.
In order to identify human-specific epigenetic influences, we have performed a novel genome-wide DNA methylation analysis comparing human, chimpanzee and rhesus macaque.
We have identified that the immunological Leukotriene B4 receptor (LTB4R, BLT1 receptor) is the most epigenetically divergent human gene in peripheral blood in comparison with other primates. This difference is due to the co-ordinated active state of human-specific hypomethylation in the promoter and human-specific increased gene body methylation. This gene is significant in innate immunity and the LTB4/LTB4R pathway is involved in the pathogenesis of the spectrum of human inflammatory diseases. This finding was confirmed by additional neutrophil-only DNA methylome and lymphoblastoid H3K4me3 chromatin comparative data. Additionally we show through functional analysis that this receptor has increased expression and a higher response to the LTB4 ligand in human versus rhesus macaque peripheral blood mononuclear cells. Genome-wide we also find human species-specific differentially methylated regions (human s-DMRs) are more prevalent in CpG island shores than within the islands themselves, and within the latter are associated with the CTCF motif.
This result further emphasises the exclusive nature of the human immunological system, its divergent adaptation even from very closely related primates, and the power of comparative epigenomics to identify and understand human uniqueness.
In the past half century there has been a dramatic increase in chronic inflammatory and metabolic common human diseases [1, 2]. This is too rapid a time frame to be due to changes in common genetic allele frequency; therefore, the environment of this increasingly urbanized human population is implicated . These environmental effects can influence the genome through the more malleable epigenome [4, 5]. Thus, by studying this level of regulation, we may identify the genes and pathways modified by this modern environment that are involved in the pathogenesis of these diseases.
Over a far longer evolutionary term, the environmental pressure experienced since the last common primate ancestor has been a major driver in the human phenotype we see today. Comparative genetic analyses of human  with close relatives such as the chimpanzee , rhesus macaque  and other primates [9–12] have explored these alterations at the sequence level. Accelerated regions have highlighted potentially modified neurological [13, 14] and anatomic  pathways, human-specific duplications have been implicated in enhanced neuronal migration  and human-specific deletions have identified lost regulatory regions . Comparative epigenomics has only recently started to be explored in a number of tissues [18–22]. This regulatory level has a plausible role as a rapid adaptive response mechanism to external change. It is therefore hypothesized to hold insights into recent environmental effects impacting on genomic activity [23, 24]. Furthermore, this regulatory change can occur gradually, minimizing effects on viability.
To fully understand human susceptibility to disease we require precise knowledge of what makes us unique. This will be the result of the integration of genetic and epigenetic differences, and their interplay with the environment. The aim of this study was to find the most substantial human species-specific DNA methylome variation in peripheral blood, a tissue type integral to immune responses, by triangulation analysis  between human, chimpanzee, and rhesus macaque. The epigenomic mark of DNA methylation is critical for development and is strongly associated with gene regulation. This human-specific DNA methylome variation will be driven by obligatory and facilitative genetic factors, as well as pure epigenetic effects . Genetic divergence in orthologous loci will contribute to methylome variation; by alterations in transcription factor binding sites (TFBSs) that act as methylation determining regions (MDRs) ; or associated with CpG density changes, due to gain or loss of species-specific CpGs , impacting on CpG island (CpGi) strength [29–31]. The major proportional components of peripheral blood cells are similar between all three species [32–34]; therefore, methylation differences affecting the major cell type fractions or multiple subtypes will be detected. Finally there is the possibility that changes in DNA methylation may be driven as a response to the different environmental stimuli encountered by the primates .
The observed methylation state may play an active role in expression, or may passively accumulate due to lack of transcription factor binding, but in either situation is still informative of robust variation in effectors and regulation. Thus, the possible identification of environmentally driven DNA methylation variation between the species could be highly informative in terms of the analysis of the pathophysiology of human disease, particularly if this was able to quantify modern lifestyle disease risk factors [4, 36]. The various environmental pressures, bottlenecks, and drift over the separate courses of these species since the last common ancestor of human and chimpanzee (approximately 5 to 6 million years ago (MYA))  will have then left not only a genetic  but also an epigenetic human-specific signature. Understanding these modifications will reveal insights into human-specific physiology and potentially human-specific responses to environmental change and vulnerabilities to disease .
DNA was extracted from peripheral blood of five chimpanzees (Pan troglodytes; three males, two females) and five rhesus macaques (Macaca mulatta; three males, two females). Samples were taken from captive individuals at Tierpark Nordhorn, Basel Zoo, Leipzig Zoo and at the German Primate Centre during routine health checks and not specifically for this study. Microsatellite analysis conducted at the German Primate Centre verified that respective individuals are not related. Sample collection adhered to the American Society of Primatologists' Principles for the Ethical Treatment of Non-Human Primates . Human DNA was derived from 10 anonymous healthy human subjects (also 60% male).
Analysis was performed in whole peripheral blood in all primates in order to identify significant outlier changes. The proportional makeup of cells present in the blood, which are predominately neutrophils (47 to 67%) and lymphocytes (28 to 37%) in human , has very similar ratios in chimpanzee (male neutrophils 59%, lymphocytes 28%; female neutrophils 49%, lymphocytes 41% ) and rhesus macaque (neutrophils approximately 67% and lymphocytes 26 to 35% ). Whilst additional populations of cells, as well as extra subtypes such as lymphocytes, T cells (CD4, CD8, and others), B cells, natural killer cells, and others, will each possess their unique subtle signatures, this is far below the resolution of this study as only large significant global changes were examined.
DNA samples were pooled for each species at equal concentration for each individual to obtain averaged methylomes and reduce individual genotypic polymorphism effects. Methylated DNA immunoprecipitation (MeDIP was then executed according to the Auto-MeDIP-seq protocol as described in Butcher and Beck  and sequenced on an Illumina GAIIx. We generated a data set of over 171 million uniquely mapped fragments (>342 million mapped paired-end reads). Of these, 40,797,356 were human, 64,610,346 were chimpanzee and 65,824,761 macaque. This was performed with paired end reads of 36 bp with average fragment sizes of 197 bp in human, 222 bp in chimpanzee, and 217 bp in macaque.
MeDIP-seq data were processed using MeDUSA v1.0 (Methylated DNA Utility for Sequence Analysis) . Sequence quality control was performed using FASTQC . MeDUSA utilized BWA (v0.5.8)  for alignment to reference genomes. Human was aligned to hg19, chimpanzee to panTro2 and macaque to rheMac2 obtained from University of California Santa Cruz (UCSC). Alignment was performed with default parameters. Following alignment a number of filtering steps were performed. Initially SAMtools v0.1.18  was used to remove reads with low alignment score (q < 10). Only those forming a correctly aligned pair were kept. A final filtering step removed potential PCR artifacts by discarding all but one read pair within groups of non-unique fragments. MeDIP specific quality control was performed using MEDIPS (v1.0) . Bigwig files representing normalized read depth (reads per million or rpm) were generated for viewing in the relevant genomes. The corresponding methylome data were deposited in Gene Expression Omnibus (GEO) under accession number GSE48942.
Differential methylation region analysis
Regions enriched for methylated reads were identified using both MACS v1.4.1  and BayesPeak v1.8 . Both programs were run using the alignment files obtained from their native genome (that is, prior to liftOver). MACS was run using the parameters -nomodel with bandwidth set to the alignment calculated fragment length (human = 197, chimpanzee = 220 and rhesus macaque = 217) and shiftsize to half the fragment length. A P-value threshold of 1 × 10-3 was used for peak selection. Bayes-Peak was run with default values in multicore mode. Having identified the regions enriched for MeDIP-seq signal, to enable comparative analysis of these peak regions, the MACS and BayesPeak output for chimpanzee and rhesus macaque was lifted to human hg19 using liftOver  with minMatch = 0.7. This parameter is appropriate due the recent shared evolutionary history of these species . Subsequently, human-specific differentially methylated regions (s-DMRs) were located through a series of intersection analyses utilizing the BEDTools  command intersectBed (v2.10) with a minimum overlap threshold of 1 bp.
Regions hypomethylated in human were defined as regions containing peaks in chimpanzee and rhesus macaque but not in human. Firstly, a single peak set was defined for both chimpanzee and macaque by excluding any MACS peaks that failed to intersect with a Bayes-Peak call. Secondly, shared peaks were identified between the chimpanzee and macaque peak set prior to intersecting with both the human MACS and human BayesPeak set. Hypomethylated s-DMRs were defined by the presence of a shared chimpanzee/rhesus macaque peak in the absence of a human peak from either software output. The chimpanzee MACS peak location lifted to human determined the coordinates of the s-DMR. A similar process was followed to identify hypermethylated s-DMRs in human. Any regions found in both MACS and BayesPeak in human, but that did not intersect any peak in chimpanzee or rhesus macaque, were isolated. Potentially, false positive hypermethylated regions could be generated, caused by an issue in the conversion from non-human primate genome to human genome. To remove such errors, all potential hypermethylated regions were lifted back to their non-human primate genome and then again to human. Only those peaks mapping back to the original location were maintained as hypermethylated s-DMRs for further analysis. We furthermore excluded all of the declared poor mappability regions as identified by ENCODE .
To benchmark the robustness of our peak calling method we compared the raw read counts for each species within their native genomes within the lifted s-DMRs coordinates. Normalizing reads between all three of the species to a human constant ((Reads/Total aligned individual species reads) × Total human reads), we showed that there were significant differences in methylation signal prior to liftOver. The comparisons between human and chimpanzee and human and rhesus macaque are significantly greater in the hypermethylated s-DMRs and are significantly less in the hypomethylated s-DMRs, versus the difference between chimpanzee and rhesus macaque in these locations (all human differences versus non-human differences pairwise comparisons; Wilcoxon tests all P < 2.2 × 10-16; data not shown).
s-DMR enrichment analysis
Feature enrichment data were visualized using Epiexplorer . Significance was calculated by Genomic Hyperbrowser  via hypothesis overlap testing with ChromHMM GM12878 (P-values were computed under the null model defined by the following preservation and randomization rules: preserve segments (T2), segment lengths and inter-segment gaps (T1); randomize positions (T1), Monte Carlo (MC) false discovery rate (FDR) threshold 0.005). Gene enrichment via GREAT 2.0.2 Region-Based Binomial Analysis was performed with the default Basal + extension parameters (constitutive 5.0 kb upstream and 1.0 kb downstream, up to 1000.0 kb max extension). Curated regulatory domains were included. For the hypermethylated s-DMR repeat enrichment, the control set was calculated within Epiexplorer  by reshuffling the genomic positions while retaining the overall number of regions and the distribution of region sizes.
CpG islands and shores
CpGi annotation was obtained from Ensembl (build 64). CpGi shores were defined as regions extending 2,000 bp upstream and downstream of each CpGi. The BEDTools  function intersectBed, using a minimum overlap of 0.5, was used to determine the number of s-DMRs within each feature.
Transcription factor binding sites
FASTA sequence for each of the s-DMRs was obtained from the reference sequence for human and the orthologous sequence obtained for chimpanzee and rhesus macaque via LiftOver  with the BEDTools  FastaFromBed function. The transcription factor affinity prediction tool TRAP  (multiple sequences) was implemented with the 904 TRANSFAC  motifs (transfac_2010.1 vertebrates), background model of human promoters, and Benjamini-Hochberg multiple test correction. DMRs greater than 5 kb were excluded from this analysis. The SP1 motifs investigated included SP1_01, SP1_Q6, SP1_Q6_01, SP1_Q4_01, SP1_Q2_01 and SP1_02; and the RFX family members included RFX1_01, RFX1_02, RFX_Q6 and EFC_Q6. The empirical P-value calculation for CTCF was derived by randomization with R (4 by 1,000×). A set of random hypomethylated and hypermethylated s-DMRs, from the total set of shared chimpanzee and rhesus macaque peaks, and all human peaks, respectively, that overlap CpGis were selected and TRAP  motif prediction scores were calculated for CTCF_01 and CTCF_02 in the human as well as orthologous chimpanzee and rhesus macaque sequences (via LiftOver and FastaFromBed). Those results that exceeded the observed human CTCF motif divergence were then calculated.
Exons and gene body methylation
Hypomethylated sDMRs located in CpGis were further filtered to obtain a subset associated with 5′ promoter regions. The exons associated with the selected promoters were identified and their methylation status was compared in order to detect inverse changes in gene body methylation. This was performed by identifying MACS peak regions in the human data within the exons and comparing the MACS score (the transformed P-value attributed to the peak) between species. If no peak was found within an exon for either chimpanzee or rhesus macaque, a score of 0 was given. The resulting regions were ranked according to differential MACS score.
All other statistical analyses were performed in the R environment . Chi-squared calculations for enrichment for hypomethylated s-DMRs, were compared with the combined chimpanzee and rhesus macaque peak set (94,799 peaks) and those for hypermethylated s-DMRs, with the total human peak set (133,494). For the genomic feature intersection calculations, the Bayes Peak confirmed MACS peak coordinates were used.
Additional epigenomic datasets
Neutrophil MethylSeq (digestion with HpaII followed by sequencing) data from Martin et al.  for four humans and four chimpanzees was accessed from GEO accession number GSE22376. Results were expressed as the probability p(U) that the site is unmethylated, and vary between 0 (methylated site) and 1 (unmethylated site). The probability of being methylated was then taken as 1 - p(U), with 0 = unmethylated and 1 = methylated. Values across the two human CpGis and chimpanzee CpGis were averaged across the island. Comparative H3K4me3 ChIP-seq data from Cain et al.  were accessed via GEO (GSE24111). These data were derived from B cells (lymphoblastoid cell line) from three humans, three chimpanzees and three rhesus macaques.
Functional analysis of Leukotriene B4 Receptor
LTB4 and LY293111 were purchased from Cayman Chemical (Ann Arbor, MI, USA). Calcium ionophore (A23187) was purchased from Sigma-Aldrich (Dorset, UK). All other materials were purchased from Life Technologies (Carlsbad, CA, USA) unless otherwise stated.
Isolation of peripheral blood mononuclear cells
Isolation of human PBMCs was performed by density centrifugation using Polymorphprep (Axis-Shield, Oslo, Norway) according to the manufacturer’s protocol. The study was approved by the Research Ethics Committee of Guy’s Hospital REC (no: 09/H804/077) and participants had provided written consent prior to any procedure. In brief, blood collected over 5.4 mmol/L EDTA was layered over an equivalent volume of Polymorphprep and centrifuged (500 g) for 35 minutes at 20°C. PBMCs collected were resuspended in 1:1 volume of RPMI:H2O to restore osmolarity. After incubation with RBC lysis buffer (eBiosciences, San Diego, CA, USA) to eliminate erythrocyte contamination, cells were then washed twice in phosphate-buffered saline and counted using a NucleoCounter (Chemometec, Allerød, Denmark). Isolation of rhesus macaque PBMCs was performed at the German Primate Centre by standard protocols for human PBMCs using Biocoll 1.077 g/ml (Biochrome, Cambridge, UK) with an extended centrifugation step (45 minutes) for separation of PBMCs. Separate PMBC aliquots in Trizol and in medium (HEPES-buffered RPMI + 10% fetal calf serum + Pen/Strep) were transported.
Total RNA was extracted from cells lysed in Trizol (Life Technologies) using QIAshredder columns and a RNeasy mini kit (QIAGEN, Venlo, Limburg, Netherlands) with TURBO DNase (Ambion, Life technologies, Carlsbad, CA, USA) treatment following the manufacturer’s protocol. RNA was quantified using a NanoDrop ND 1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA) and ND-1000 software version 3.2.0. Reverse transcription of RNA was carried out using RevertAid M-MuLV reverse transcriptase (Fermentas, Thermo Scientific, Waltham, MA, USA) and primed using random hexamers according to the manufacturer’s protocol. Relative mRNA expression levels were measured using TaqMan Gene Expression Master Mix (Applied Biosystems, Life technologies, Carlsbad, CA, USA) and the following TaqMan probe sets: LTB4R- Hs01938704_s1 and 18s RNA- 4319413E, matching both human and rhesus macaque sequences. RT-PCR was performed on a ViiA 7 Real-Time PCR system (Applied Biosystems) and analyzed using ViiA 7 software version 1.0.
Calcium mobilization assay
Calcium mobilization assays were conducted using a FLIPR calcium 4 assay kit (Molecular Devices, Sunnydale, CA, USA) as described previously [58, 59]. PBMCs (2 × 105/well) were plated into poly-D-lysine-coated 96-well plates in RPMI 1640 supplemented with 10 mmol/L HEPES. Cells were incubated for 1 hour with FLIPR loading buffer prior to addition of ligand and fluorescent intensity was measured at 37°C using a Flexstation 3 (Molecular Devices). Results were analyzed with SoftMax Pro Software (Molecular Devices).
Chimpanzee and rhesus macaque samples were taken from captive individuals at Tierpark Nordhorn, Basel Zoo, Leipzig Zoo and at the German Primate Centre during routine health checks and not specifically for this study. Sample collection adhered to the American Society of Primatologists' Principles for the Ethical Treatment of Non-Human Primates . For the isolation of human PBMC samples for functional analysis, the study was approved by the Research Ethics Committee of Guy’s Hospital (REC no: 09/H804/077) and participants had provided written consent prior to any procedure.
Primate DNA methylomes
The DNA methylomes of peripheral blood DNA samples were generated from pooled individuals in order to reduce the effects of individual genetic variability. These were of healthy unrelated primates, comprising 10 humans, 5 chimpanzees and 5 rhesus macaques, all with a 60% male split. Also, uncultured cells were used to avoid the incorporation of additional stochastic artifacts  Automated MeDIP-seq was performed as previously described , and data were processed using the MeDUSA pipeline .
The DNA methylomes are displayable in the individual genomes, human (GRCh37), chimpanzee (CGCS 2.1/PanTrog2) and rhesus macaque (MGSC Merged 1.0/rheMac2) in the UCSC genome browser in the context of existing annotation and are available at . The chimpanzee and macaque results can also be viewed on the human sequence having been converted using the liftOver utility  to enable an initial comparative view.
Human species-specific differentially methylated regions
To find loci with human s-DMRs we first identified robust enrichment peaks using the MACS v1.4.1  peak-calling algorithm, but with the additional conservative requirement that regions must also have an overlap with a BayesPeak v1.8  algorithm result (see Methods). This additional step reduced the number of potential peaks by up to approximately 8%. We then identified all human peaks where no chimpanzee or rhesus macaque peak was present in the orthologous location (with either algorithm) to define hypermethylated human s-DMRs, and hypomethylated human s-DMRs where no human peak was present (again with either algorithm) but where both chimpanzee and macaque peaks were located in common.
A total of 22,758 hypomethylated and 15,858 hypermethylated human s-DMRs were identified. All of these lie in strongly sequence-similar regions between the three primate genomes, but are also unique within each primate genome due to the requirement of these s-DMRs to be able to pass through consistent reciprocal liftOver  steps across the primate genomes. Therefore, whilst genetic influence on the methylome is strong and high sequence similarity exists, particularly between human and chimpanzee (approximately 98.6% ) and primates overall [63, 64], substantial numbers of s-DMRs were identified. Without the additional out-group of rhesus macaque the calculated human-specific set via only a comparison with chimpanzee would be 69.2% larger; thus, its inclusion increases our power to identify true human-specific modification. The number of hypomethylated human s-DMRs is moderately greater compared with hypermethylated s-DMRs, which is likely due to an increased loss within the human peak set, as these must pass through reciprocal liftOver to the two slightly less well-characterized genome sequences of chimpanzee and macaque. Consistent with this, peaks from human chromosome 21 lifted via the higher quality chimpanzee chromosome 21  retained the highest proportion of peaks of any chromosome (data not shown). These regions are visible as hypomethylated (yellow) and hypermethylated (blue) s-DMRs also via the above UCSC tracks link and are visualized genome-wide in Figure 1.
Genome feature annotation enrichment
We compared the location of human s-DMRs to genomic functional annotation, as defined by the ChromHMM segmentation analysis , using the HapMap line interrogated by Encode (B-lymphocyte GM12878 data). We calculated significance via Epiexplorer  and the Genomic Hyperbrowser . This showed a similar pattern in both hypo- and hypermethylated s-DMRs, with a significant skew towards weak promoters, gene transcriptional regions, and both strong and weak enhancers, but depletion within strong promoters (Figure 2A,B). This is comparable to the finding in tissue and reprogramming s-DMRs, with most DMRs in moderate CpG dense regions . Increased levels of s-DMRs in these transcription-associated regions may be indicative of the potential role methylation variability plays in gene body methylation  and splicing [69, 70]. The number of hypomethylated s-DMRs residing over exons is 3,444 and over introns is 9,848, and for hypermethylated s-DMRs these numbers are 1,042 for exons and 6,514 for introns. The more consistent and stronger results of hypomethylated s-DMRs in enhancer regions is consistent with studies showing these regions may be protected from DNA methylation [71, 72] and are more species-specific. Furthermore s-DMR genomic regions are enriched in comparison to the rest of the genome for the dynamic fraction of the DNA methylome as defined by Ziller et al.  (data not shown, χ2P < 2.2 × 10-16).
Gene enrichment analysis of these s-DMR regions was performed using the Genomic Regions Enrichment of Annotations Tool (GREAT) . A large number of significant biological, human phenotype and disease-related enrichments were identified by FDR-corrected region-based binomial analysis (Additional file 1). A smaller set of intriguing categories were found to be significant in both the hypo- and hypermethylated s-DMR enrichment categories. This included monocyte activation and regulation of prostaglandin-endoperoxide synthase activity for biological processes; hypotrichosis for human phenotype; and diverticulitis for disease ontology, which has not been observed in non-human primates to date . The evolutionary timescale difference between human and chimpanzee divergence (5 to 6 MYA) and human population-based polymorphisms (<1 MYA) is large. Therefore, many potential genetic disease susceptibilities will be fixed between the species. Nevertheless genome-wide association study SNPs  are enriched in these s-DMR regions by approximately 1.40× compared to the genome average (data not shown; χ2P = 5.09 × 10-3).
The hypermethylated s-DMR set was found to be increased within repetitive elements, specifically the SINE group (Figure 3; via Epiexplorer ). Examination of this result showed that these hypermethylated s-DMRs were disproportionally increased within the second oldest AluS subcategory (χ2P < 2.2 × 10-16), which still possesses mobilization ability , whilst being proportionally reduced within both the most ancient AluJ and youngest AluY categories (Additional file 2). Moreover, hypomethylated s-DMRs also showed this pattern, with an increase in AluS (χ2P < 2.2 × 10-16), but also showed a slight increase in the youngest and most active AluY (χ2P = 1.15 × 10-4). Transposable elements constitute nearly half of the primate genome and there is increasing evidence for their functional role in influencing expression and potential modulation in human pathology . Primate-specific open chromatin regions (DNAse I hypersensitivity sites) are enriched for transposable elements (approximately 63%)  and methylation changes in trans-species experiments have also indicated their regulatory potential .
s-DMRs are more prevalent in CpG shores than CpG islands
Previously identified tissue-, cancer- and reprogramming-specific DMRs [81, 82] have all been found to be more prevalent in moderate CpG dense regions surrounding CpGis, termed CpGi ‘shores’, than within the island themselves. Within the CpGis and CpGi shores, we identified 77 hypo- and 45 hypermethylated human s-DMRs, and 821 hypo- and 431 hypermethylated s-DMRs, respectively (Additional file 3). After accounting for the almost fourfold larger genome size taken up by shores (approximately 89.5 Mb) than islands (approximately 23.8 Mb), these s-DMRs were still more prevalent within shore regions (Figure 4; Wilcoxon test, both hypo- and hypermethylated, P < 2.2 × 10-16). Additionally, this was also significant when comparing island and shore s-DMRs with regard to all possible locations of hypomethylated s-DMRs (all co-locating chimpanzee and macaque peaks in these regions; χ2P = 4.80 × 10-29) and hypermethylated s-DMRs (from location of all human peaks in these regions; χ2P = 2.74 × 10-8) and for s-DMRs combined (χ2P = 2.94 × 10-32). Therefore, the proportion of s-DMRs in CpGi shores is far higher than in islands (14.0 compared to 5.13 per Mb of feature sequence). However, other features are even more concentrated with s-DMRs: exons (23.7/Mb), CTCF sites (26.8/Mb) and DNAse I hypersensitivity sites (40.3/Mb).
Sequence similarity of s-DMRs in CpG islands
Sequence divergence could be contributing to the s-DMRs within these CpGis . Therefore, to assess whether there was greater divergence between species within s-DMRs, we compared the sequence similarity between the s-DMRs and the total background set of peaks that these were identified from. However, no significant difference in sequence change between human and chimpanzee was seen between hypomethylated s-DMR peaks and the total combined non-primate peaks, or hypermethylated s-DMR peaks and the total set of human peaks within CpGis (Wilcoxon test, both non-significant, P > 0.05; Figure 5). The observed sequence divergence for all four subsets is at the expected global genomic level of approximately 99% sequence identity. This does not rule out sequence influence within these s-DMRs, but implies that more severe genetic differences would not have passed our reciprocal liftOver requirements and within these s-DMR regions sequence change is not more prevalent than the genome average. Therefore, only minor levels of base change may be contributing to this methylation variation.
Transcription factor binding site modification
Sequence alteration in certain TFBS motifs within CpGis have been shown to modify local methylation and thereby act as MDRs . With the availability of both human s-DMRs from our dataset and sequence data from these three species, we were able to interrogate these regions for potential MDRs. Isolated species-specific genetic changes significantly modifying TFBS motif strength were looked for. To do this we utilized the Transcription Factor Affinity Prediction (TRAP) tool  and calculated the binding prediction of 904 TFBSs (TRANSFAC vertebrate V$motifs ) within CpGi-related s-DMRs in human, as well as the orthologous chimpanzee and macaque sequences.
We performed a global comparison for the total set of hypomethylated and hypermethylated CpGi s-DMRs by calculating a single total binding P-value (Benjamini-Hochberg corrected) for each motif in each of the three species. This could be performed in these CpG dense, predominately hypomethylated and classical promoter regions as a slight MDR motif sequence variation would be detectable in these highly similar orthologous sequences. We visualized this by plotting in two dimensions the difference in P-values (-log10) between human and chimpanzee (on the y-axis) and human and macaque (on the x-axis) binding scores to identify consistent human directional change (Figure 6A,B). Both CTCF motifs (CTCF_01 and CTCF_02) showed consistently increased and decreased motif binding in human hypomethylated and hypermethylated s-DMRs, respectively (blue dots in Figure 6A,B; Additional file 4; empirical P = 0.046). That is, there was an increase in human CTCF binding compared to both chimpanzee and macaque within hypomethylated s-DMRs, and inversely a consistent decrease in CTCF binding in hypermethylated s-DMRs. The CTCF motif has been previously identified as an MDR . This trans-species effect is consistent with the disruption of CTCF occupancy being linked to increased methylation, recently observed across numerous human cell types , as well as the different role of CTCF binding regions within open and closed chromatin . However, other potential MDRs, such as SP1 motifs [85, 86] and members of the RFX family, did not show consistent human directional divergence in this global analysis. One of the most extreme outliers within the hypermethylated set was for increased human binding strength of a MeCP2 (MECP2_02) motif. Due to the critical role of this gene in development and neural function [87, 88], this may represent motif variation that is contributing genetically to epigenetic effects across multiple tissues.
Canonical 5′ promoter s-DMRs
From the total s-DMRs within CpGis, a smaller subset of 53 hypomethylated and 25 hypermethylated s-DMRs reside within the canonical 5′ promoter region, and as such would have a recognized strong potential association with the expression of these genes .
As highly expressed genes are associated with low gene promoter methylation and high gene body methylation , we ranked these loci by differential gene body methylation between human and averaged chimpanzee and rhesus macaque methylation levels in the consensus transcripts (CCDC) that arise from these promoters (Table S4A in Additional file 5). The top hypomethylated promoter with greatest inverse differential in gene body methylation was the major isoform of LTB4R (Leukotriene B4 Receptor; Figure 7). This 5′ CpGi in fact also possesses two hypomethylated s-DMRs, with one also overlapping the CpGi shore region. This promoter also overlaps the Leukotriene B4 receptor 2 gene (LTB4R2)  and is a bidirectional promoter for CIDEB transcribed in the reverse direction. As mentioned, the major isoform of LTB4R showed higher exonic methylation than both chimpanzee and rhesus macaque (MACS -log10P-value 2943.7, versus 2397.8 and 2180.8) than any of the other 5′ CpGi hypomethylated genes, with these other co-locating genes (LTB4R2 and CIDEB) showing the inverse result (Table 1). We excluded any build inconsistencies influencing this result by also re-aligning to additional available primate builds (CSAC 2.1.4/PanTrog4 and BGI CR_1.0/rheMac3) and these all showed the identical comparative methylation difference and peak calls within this LTB4R locus. The gene body methylation analysis within the set of hypermethylated promoter s-DMRs revealed the genes GIT1, ABCG4 and RUSC2 as the highest ranked. However, in this analysis there was not one strong outlier identified, as there was for the LTB4R result (Table S4B in Additional file 5).
Leukotriene B4 receptor
To investigate the strong comparative methylation difference we had identified in LTB4R, we examined the available CpGi-focused methylation-sensitive sequence in neutrophils from Martin et al.  for four human and four chimpanzees. This replicated these findings within the promoter and exonic CpGis with average methylation of 18.4% and 70.3% across the orthologous CpGi promoter of LTB4R, and 75.9% and 62.0% for the exonic CpGi in human and chimpanzee, respectively (Figure 8; Wilcoxon P <0.05 for both). This >50% difference in promoter methylation is a level of change only seen within humans in cancer tissue studies, as opposed to the very small variations identified to date in common non-malignant diseases. We also investigated the available human peripheral blood bisulfite-seq DNA methylome from Li et al. , which also supports the findings in human with a methylation level of 3.43% within the promoter CpGi and 82.22% in the gene-body exonic CpGi.
Primate comparative chromatin data were available from B-cell lymphoblastoid cell lines . This identified a significant activating H3K4me3 peak in the LTB4R promoter of three human samples, but no declared peak in three chimpanzees or three rhesus macaque samples (Figure 9). Thus, although the initial findings were identified in mixed peripheral blood, the strong effect is corroborated in two major blood components; neutrophils and B lymphocyte-derived.
Therefore, these results in individual cell types support that variation in cell type composition is not a potential confounding factor in this DNA methylome analysis from whole blood-derived DNA. However, we further excluded this by examining the 500 leukocyte subtype-related differentially methylated CpG positions (L-DMPs) identified by Houseman et al. . We find that none of these CpGs overlap with the entire LTB4R locus. Furthermore, just 5 of our 15,858 hypermethylated DMRs (0.032%) and 6 of our 22,758 hypomethylated DMRs (0.026%) show overlap with these L-DMPs, respectively. Therefore, only 11 or approximately 2.2% of these 500 L-DMPs co-localize with any of the s-DMRs identified, indicating that our results are not significantly enriched for cell type-associated methylation changes.
Furthermore, we examined intra-human variability, via the MARMAL-AID repository  of all performed 450k DNA methylation array experiments, which includes 1,665 whole blood samples from healthy and non-cancer disease subjects. These data, even derived from multiple studies with significant experimental/batch variation, also confirm consistently low human LTB4R promoter methylation (mean 29.4%, standard deviation 7.4% and high LTB4R gene body methylation (mean 87.9%, standard deviation 2.3%) (Figure 10).
Methylation of this LTB4R promoter has previously been established by Kato et al.  to reduce this gene's expression. LTB4R (aka BLT1 receptor) is a high affinity G-protein-coupled receptor for LTB4, a potent chemoattractant involved in inflammation and immune response in the eicosanoid signaling (leukotriene and prostaglandin) pathway. It has been implicated in the whole spectrum of inflammatory diseases , including asthma , inflammatory arthritis , atherosclerosis [99, 100], inflammatory bowel disease , and psoriasis . It has also previously been postulated to play a role in HIV infection . LTB4R is highly expressed across all blood cell types, particularly in myeloid cells such as neutrophils and monocytes (Additional file 6). The LTB4-LTB4R axis is involved in linking early immune system activation, neutrophil auto-signaling and swarming , and early effector T-cell recruitment , acting as a potent non-chemokine pathway for cytotoxic effector cell traffic . Interestingly, the Gene Ontology category for the ligand of this receptor, ‘Leukotriene Production involved in Inflammatory Response’ (GO:0002540), was significantly enriched in the GREAT analysis for hypomethylated s-DMRs (P = 2.76 × 10-5, FDR q = 2.77 × 10-4). This is due to s-DMRs associated with genes such as ALOX5 (5-LO) and ALOX5AP (FLAP) within this pathway. Furthermore, the gene encoding LTA4 hydrolase (LTA4H), which catalyzes the conversion of LTA4 to LTB4, also contains two intragenic s-DMRs.
Whilst CpG density is inversely associated with CpGi methylation state [28, 29], the human LTB4R promoter island is slightly less CpG dense than the other primates (18.7% in human, 19.1% in chimpanzee and 19.3% in rhesus macaque), so from this could be expected to be more, not less, methylated. The total sequence similarity is 98.3% between human and chimpanzee and 96.3% between human and rhesus macaque. We examined the comparative TRANSFAC TFBS motif analysis via TRAP for the LTB4R-associated s-DMRs. This revealed no significant difference in known MDR’s CTCF or SP1 motifs in this region. The only highly significant human-specific gain of binding affinity within this island was for Rfx1 (Regulatory factor X1, aka Enhancer factor C, V$ECF_Q6), with a significant human motif binding P-value of 0.042, and non-significant chimpanzee (P = 0.635) and rhesus macaque (P = 0.491) prediction (Additional file 7). As previously stated, Rfx2, another member of the RFX winged-helix transcription factor motifs, has been found to contribute to methylation variation in MDRs , although whether one motif would be expected to lead to such a considerable difference in methylation is not certain. Lienert et al.  identified an incremental additive effect of multiple MDR motifs required to modify methylation to this level, thus possibly implicating additional non-genetic or environmental factors in this locus. The RFX1 motif shows human-specific binding due to a central T present at chr14:24779978, which is the C of a CpG in all other primates. There is also no human common polymorphism at this site (dbSNP138). Furthermore, it is interesting to note that the modern human motif sequence is identical in both the archaic Denisovan hominin and the recent high coverage Altai Neanderthal sequences, thus dating this particular genetic change to at least prior to approximately 600,000 years ago [107, 108].
Different LTB4R (BLT1) mRNA expression and signaling in human and rhesus macaque cells
To investigate the potential effects of the different epigenetic states of the LTB4R (BLT1) promoter, we performed comparative functional analysis between human and rhesus macaque PBMCs for differential mRNA expression and signaling. We evaluated the mRNA expression of LTB4R and the 18s gene as control using TaqMan primers and probes that are complementary to the gene sequence in both species. Similar levels of internal control 18s mRNA were observed in human and rhesus macaque PBMCs, while 40-fold higher levels of LTB4R mRNA were identified in human cells (Figure 11A).
LTB4 signaling through the LTB4R/BLT1 receptor induces calcium mobilization in cells known to express the receptor . To determine whether the increased expression of LTB4R mRNA in human PBMCs resulted in functional expression of LTB4R/BLT1, human and rhesus macaque PBMCs were stimulated with different concentrations of LTB4 and intracellular calcium flux was measured (as described in Methods). Human and rhesus macaque PBMCs responded similarly to a non-specific calcium activator, the calcium ionophore A23187 (1 μmol/L). However, concentration-dependent responses to LTB4 were observed in human cells in contrast to rhesus macaque cells where no specific calcium mobilization was identified (Figure 11B).
The LTB4R/BLT1 selective antagonist LY223982 was used to confirm specificity of calcium responses and full inhibition of calcium mobilization in response to LTB4 was observed in human PBMCs (Figure 11C). These data support that the identified human-specific differential methylation pattern within the LTB4R locus may affect mRNA expression, signaling and function of the LTB4R/BLT1 receptor.
Human uniqueness has arisen due to the accumulation of genetic, environmental, behavioral and cultural changes [39, 109]. Subtle variation within all these factors, separately and in complex combinations, also contributes to disease susceptibility. Additionally, the more extreme differences between species contribute to species-specific disease variation . Higher prevalence of some common human diseases has been postulated, even whilst taking environmental differences into consideration, including Alzheimer’s disease, coronary artery disease and some cancers . Humans have also been observed, comparatively, to have an excessively heightened inflammatory response, with genetic contributions identified, including SIGLEC genes [112, 113].
By comparative tri-primate peripheral blood DNA methylation analysis we have identified strong DNA methylation variation between three primate species, human, chimpanzee and macaque, and located human-specific DMRs with potential regulatory consequences. Many regulatory loci are not constrained across mammalian evolution  and genetic loss of human-specific regulatory DNA has previously highlighted regulation change as a human evolutionary mechanism .
Both the hypomethylated and hypermethylated s-DMRs identified were skewed towards weak promoters, enhancers and transcribed regions. We also identified these DMRs to be more prevalent within moderately CpG dense CpGi shores than the islands themselves, consistent with tissue-, cancer- and reprogramming-specific changes [82, 83]. Promoter hypomethylated regions that are common across diverse cell types are often found to possess a central region that remains unmethylated with an outer tidal region that shifts in a lineage-specific manner modulating associated gene expression . DNA binding proteins localizing within these boundary regions, or shores, may help define these methylation states . Thus, these species-specific modifications may be acting in similar fashion.
Pathogens have imposed a strong selective influence on the human genome . Therefore, the unique evolutionary histories of these primates post-speciation would be expected to lead to unique immune systems. These strongly divergent differences have been recently exposed to be extreme within mammals, whereby the mouse was shown not to reproduce any of the patterns of gene expression induced by human inflammatory disease [116, 117]. This poor model of human inflammation was found to perform no better than random. Human immunological variation is of particular interest due to the extreme rise of inflammatory diseases within the past few generations . Combined genetic and new environmental factors may be driving this increase, with effects on immune development proposed, and epigenetic alterations are being explored to explain these immune system-related diseases .
We identified strong human-specific variation in LTB4R, important in the innate immune leukotriene and prostaglandin signaling pathway. Hypomethylated s-DMRs were identified covering both the promoter CpG island and shore region of this gene. In addition, inverse species-specific methylation was identified in the gene body of LTB4R. These significant DNA methylation differences identified in peripheral whole blood were replicated in neutrophil-only data, within both the promoter and the gene-body region. These data were also consistent with a H3K4me3 ChIP-seq study from primate lymphocyte cell-lines, with a strong active promoter signature present only in human. Comparative functional analysis also supported these large differences found across major blood cell types by showing human-specific increased expression and ligand response in peripheral blood mononuclear cells. These results are therefore consistent with the epigenetic ‘dimmer switch’ being dialed up in this human locus, in contrast to on/off human-specific regulatory deletions . Furthermore, a recently published comparative array analysis in peripheral blood from Hernando-Herraez et al. , including further primate species (bonobo chimpanzee, gorilla and orangutan), supports the human-specific epigenetic state of the LTB4R locus. In this study three available complementary CpG dinucleotide positions located within the LTB4R promoter possessed human-specific hypomethylation.
The LTB4-LTB4R (BLT1) pathway, because of its significant immunological role, has been implicated in the pathophysiology of a number of diseases, including asthma (including airway hyper-responsiveness, severe attacks and asthma exacerbations) , atherosclerosis, Alzheimer’s disease , obesity-related inflammation , and inflammatory bowel disease . LTB4R is involved in amplifying T-cell recruitment  and a heightened T-cell response in human compared with chimpanzee has previously been identified . It is also of interest to note that up-regulation of LTB4R expression has been identified in seven out of nine experiments in human via the Expression Atlas  and furthermore that LTB4R expression on T cells is regulated by inflammation and may link innate and adaptive immune responses . A number of further promoter associated s-DMRs would also be of interest for further follow-up, including MAPK15 (mitogen-activated protein kinase 15), involved in autophagy , and the MECP2 gene, which is mutated in Rett syndrome .
The effect of genetic mutation in particular TFBSs on the methylation state of CpGis has recently been refined by Lienert et al. . Our identification of human divergent CTCF and RFX1 motif changes are consistent with these MDR data. These thus shift the genetic set point of methylation, within these promoters, between these primates. Additional environmental factors may then lead to further subtle variation.
Furthermore, from this comparative epigenomic work we propose the ‘s-DMR hypothesis’ whereby regions of significant epigenetic difference between humans and other primates that are highly conserved at the sequence level may include loci subject to, or indicative of, human-specific environmental influence, potentially in association with pathogenic conditions. Thus, focusing on these regions may enable the identification of endemic human environmental imprints on the epigenome.
Comparison with our most closely related primate relatives enables insights into human-unique physiology. It is vital to identify and explain these differences, particularly as we readily extrapolate disease pathology from models in other species. Comparative epigenomics enables robust regulatory modifications to be pinpointed, and by integration with genetic data, a more complete functional picture can be elucidated. We have shown within blood tissue that the DNA methylome is not identically conserved between these primate species. These variations may be strongly driven by facilitative genetic means, but also potentially additional environmental factors. We have identified that the LTB4R gene is a significant human-specific epigenomic outlier, whilst containing minimal genetic differences. Thus, this finding reveals human-specific change in the innate immune system that may be a human-specific susceptibility factor, an inbuilt resultant common primate response to modern human conditions, or a combination of both, possibly contributing to the high level of these common diseases in the current environment.
differentially methylated region
false discovery rate
Gene Expression Omnibus
leukocyte subtype-related differentially methylated CpG position
methylation determining region
methylated DNA immunoprecipitation
million years ago
peripheral blood mononuclear cell
polymerase chain reaction
short interspersed nucleotide element
transcription factor binding site
University of California Santa Cruz.
Eder W, Ege MJ, von Mutius E: The asthma epidemic. N Engl J Med. 2006, 355: 2226-2235.
Olshansky SJ, Passaro DJ, Hershow RC, Layden J, Carnes BA, Brody J, Hayflick L, Butler RN, Allison DB, Ludwig DS: A potential decline in life expectancy in the United States in the 21st century. N Engl J Med. 2005, 352: 1138-1145.
Renz H, von Mutius E, Brandtzaeg P, Cookson WO, Autenrieth IB, Haller D: Gene-environment interactions in chronic inflammatory disease. Nat Immunol. 2011, 12: 273-277.
Feinberg AP: Phenotypic plasticity and the epigenetics of human disease. Nature. 2007, 447: 433-440.
Joubert BR, Haberg SE, Nilsen RM, Wang X, Vollset SE, Murphy SK, Huang Z, Hoyo C, Midttun O, Cupul-Uicab LA, Ueland PM, Wu MC, Nystad W, Bell DA, Peddada SD, London SJ: 450 K epigenome-wide scan identifies differential DNA methylation in newborns related to maternal smoking during pregnancy. Environ Health Perspect. 2012, 120: 1425-1431.
Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, Devon K, Dewar K, Doyle M, FitzHugh W, Funke R, Gage D, Harris K, Heaford A, Howland J, Kann L, Lehoczky J, LeVine R, McEwan P, McKernan K, Meldrim J, Mesirov JP, Miranda C, Morris W, Naylor J, Raymond C, Rosetti M, Santos R, Sheridan A, Sougnez C: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921.
The Chimpanzee Sequencing and Analysis Consortium: Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005, 437: 69-87.
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: Evolutionary and biomedical insights from the rhesus macaque genome. Science. 2007, 316: 222-234.
Scally A, Dutheil JY, Hillier LW, Jordan GE, Goodhead I, Herrero J, Hobolth A, Lappalainen T, Mailund T, Marques-Bonet T, McCarthy S, Montgomery SH, Schwalie PC, Tang YA, Ward MC, Xue Y, Yngvadottir B, Alkan C, Andersen LN, Ayub Q, Ball EV, Beal K, Bradley BJ, Chen Y, Clee CM, Fitzgerald S, Graves TA, Gu Y, Heath P, Heger A: Insights into hominid evolution from the gorilla genome sequence. Nature. 2012, 483: 169-175.
Locke DP, Hillier LW, Warren WC, Worley KC, Nazareth LV, Muzny DM, Yang SP, Wang Z, Chinwalla AT, Minx P, Mitreva M, Cook L, Delehaunty KD, Fronick C, Schmidt H, Fulton LA, Fulton RS, Nelson JO, Magrini V, Pohl C, Graves TA, Markovic C, Cree A, Dinh HH, Hume J, Kovar CL, Fowler GR, Lunter G, Meader S, Heger A: Comparative and demographic analysis of orang-utan genomes. Nature. 2011, 469: 529-533.
Prufer K, Munch K, Hellmann I, Akagi K, Miller JR, Walenz B, Koren S, Sutton G, Kodira C, Winer R, Knight JR, Mullikin JC, Meader SJ, Ponting CP, Lunter G, Higashino S, Hobolth A, Dutheil J, Karakoc E, Alkan C, Sajjadian S, Catacchio CR, Ventura M, Marques-Bonet T, Eichler EE, Andre C, Atencia R, Mugisha L, Junhold J, Patterson N: The bonobo genome compared with the chimpanzee and human genomes. Nature. 2012, 486: 527-531.
Lindblad-Toh K, Garber M, Zuk O, Lin MF, Parker BJ, Washietl S, Kheradpour P, Ernst J, Jordan G, Mauceli E, Ward LD, Lowe CB, Holloway AK, Clamp M, Gnerre S, Alfoldi J, Beal K, Chang J, Clawson H, Cuff J, Di Palma F, Fitzgerald S, Flicek P, Guttman M, Hubisz MJ, Jaffe DB, Jungreis I, Kent WJ, Kostka D, Lara M: A high-resolution map of human evolutionary constraint using 29 mammals. Nature. 2011, 478: 476-482.
Pollard KS, Salama SR, Lambert N, Lambot MA, Coppens S, Pedersen JS, Katzman S, King B, Onodera C, Siepel A, Kern AD, Dehay C, Igel H, Ares M, Vanderhaeghen P, Haussler D: An RNA gene expressed during cortical development evolved rapidly in humans. Nature. 2006, 443: 167-172.
Prabhakar S, Noonan JP, Paabo S, Rubin EM: Accelerated evolution of conserved noncoding sequences in humans. Science. 2006, 314: 786-
Prabhakar S, Visel A, Akiyama JA, Shoukry M, Lewis KD, Holt A, Plajzer-Frick I, Morrison H, Fitzpatrick DR, Afzal V, Pennacchio LA, Rubin EM, Noonan JP: Human-specific gain of function in a developmental enhancer. Science. 2008, 321: 1346-1350.
Dennis Megan Y, Nuttle X, Sudmant Peter H, Antonacci F, Graves Tina A, Nefedov M, Rosenfeld Jill A, Sajjadian S, Malig M, Kotkiewicz H, Curry CJ, Shafer S, Shaffer LG, de Jong PJ, Wilson RK, Eichler EE: Evolution of human-specific neural SRGAP2 genes by incomplete segmental duplication. Cell. 2012, 149: 912-922.
McLean CY, Reno PL, Pollen AA, Bassan AI, Capellini TD, Guenther C, Indjeian VB, Lim X, Menke DB, Schaar BT, Wenger AM, Bejerano G, Kingsley DM: Human-specific loss of regulatory DNA and the evolution of human-specific traits. Nature. 2011, 471: 216-219.
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.
Martin DI, Singer M, Dhahbi J, Mao G, Zhang L, Schroth GP, Pachter L, Boffelli D: Phyloepigenomic comparison of great apes reveals a correlation between somatic and germline methylation states. Genome Res. 2011, 21: 2049-2057.
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-
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: 1225-1234.
Shulha HP, Crisci JL, Reshetov D, Tushir JS, Cheung I, Bharadwaj R, Chou H-J, Houston IB, Peter CJ, Mitchell AC, Yao W-D, Myers RH, Chen J-F, Preuss TM, Rogaev EI, Jensen JD, Weng Z, Akbarian S: Human-specific histone methylation signatures at transcription start sites in prefrontal neurons. PLoS Biol. 2012, 10: e1001427-
Feil R, Fraga MF: Epigenetics and the environment: emerging patterns and implications. Nat Rev Genet. 2011, 13: 97-109.
Eichten S, Borevitz J: Epigenomics: Methylation’s mark on inheritance. Nature. 2013, 495: 181-182.
Harris RA, Rogers J, Milosavljevic A: Human-specific changes of genome structure detected by genomic triangulation. Science. 2007, 316: 235-237.
Richards EJ: Inherited epigenetic variation–revisiting soft inheritance. Nat Rev Genet. 2006, 7: 395-401.
Lienert F, Wirbelauer C, Som I, Dean A, Mohn F, Schubeler D: Identification of genetic elements that autonomously determine DNA methylation states. Nat Genet. 2011, 43: 1091-1097.
Bell CG, Wilson GA, Butcher LM, Roos C, Walter L, Beck S: Human-specific CpG 'beacons' identify loci associated with human-specific traits and disease. Epigenetics. 2012, 7: 1188-1199.
Weber M, Hellmann I, Stadler MB, Ramos L, Paabo S, Rebhan M, Schubeler D: Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome. Nat Genet. 2007, 39: 457-466.
Bock C, Walter J, Paulsen M, Lengauer T: CpG island mapping by epigenome prediction. PLoS Comput Biol. 2007, 3: e110-
Cohen NM, Kenigsberg E, Tanay A: Primate CpG islands are maintained by heterogeneous evolutionary regimes involving minimal selection. Cell. 2011, 145: 773-786.
Wakeman L, Al-Ismail S, Benton A, Beddall A, Gibbs A, Hartnell S, Morris K, Munro R: Robust, routine haematology reference ranges for healthy adults. Int J Lab Hematol. 2007, 29: 279-283.
Howell S, Hoffman K, Bartel L, Schwandt M, Morris J, Fritz J: Normal hematologic and serum clinical chemistry values for captive chimpanzees (Pan troglodytes). Comp Med. 2003, 53: 413-423.
Buchl SJ, Howard B: Hematologic and serum biochemical and electrolyte values in clinically normal domestically bred rhesus monkeys (Macaca mulatta) according to age, sex, and gravidity. Lab Anim Sci. 1997, 47: 528-533.
Watson JA, Watson CJ, McCann A, Baugh J: Epigenetics, the epicenter of the hypoxic response. Epigenetics. 2010, 5: 293-296.
Bock C: Analysing and interpreting DNA methylation data. Nat Rev Genet. 2012, 13: 705-719.
Olson MV, Varki A: Sequencing the chimpanzee genome: insights into human evolution and disease. Nat Rev Genet. 2003, 4: 20-28.
Hu HY, He L, Fominykh K, Yan Z, Guo S, Zhang X, Taylor MS, Tang L, Li J, Liu J, Liu J, Wang W, Yu H, Khaitovich P: Evolution of the human-specific microRNA miR-941. Nat Commun. 2012, 3: 1145-
Varki A, Geschwind DH, Eichler EE: Explaining human uniqueness: genome interactions with environment, behaviour and culture. Nat Rev Genet. 2008, 9: 749-763.
The American Society of Primatologists: Principles for the ethical treatment of non-human primates. [https://www.asp.org/society/resolutions/EthicalTreatmentOfNonHumanPrimates.cfm]
Butcher LM, Beck S: AutoMeDIP-seq: a high-throughput, whole genome, DNA methylation assay. Methods. 2010, 52: 223-231.
Wilson G, Dhami P, Feber A, Cortazar D, Suzuki Y, Schulz R, Schar P, Beck S: Resources for methylome analysis suitable for gene knockout studies of potential epigenome modifiers. Giga Sci. 2012, 1: 3-10.1186/2047-217X-1-3.
FastQC. [http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc], 
Li H, Durbin R: Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009, 25: 1754-1760.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/Map format and SAMtools. Bioinformatics. 2078–2079, 2009: 25-
Chavez L, Jozefczuk J, Grimm C, Dietrich J, Timmermann B, Lehrach H, Herwig R, Adjaye J: Computational analysis of genome-wide DNA methylation during the differentiation of human embryonic stem cells along the endodermal lineage. Genome Res. 2010, 20: 1441-1450.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-
Cairns J, Spyrou C, Stark R, Smith ML, Lynch AG, Tavare S: BayesPeak–an R package for analysing ChIP-seq data. Bioinformatics. 2011, 27: 713-714.
Fujita PA, Rhead B, Zweig AS, Hinrichs AS, Karolchik D, Cline MS, Goldman M, Barber GP, Clawson H, Coelho A, Diekhans M, Dreszer TR, Giardine BM, Harte RA, Hillman-Jackson J, Hsu F, Kirkup V, Kuhn RM, Learned K, Li CH, Meyer LR, Pohl A, Raney BJ, Rosenbloom KR, Smith KE, Haussler D, Kent WJ: The UCSC genome browser database: update 2011. Nucleic Acids Res. 2011, 39: D876-882.
Bardet AF, He Q, Zeitlinger J, Stark A: A computational pipeline for comparative ChIP-seq analyses. Nat Protoc. 2012, 7: 45-61.
Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842.
Bernstein BE, Birney E, Dunham I, Green ED, Gunter C, Snyder M: An integrated encyclopedia of DNA elements in the human genome. Nature. 2012, 489: 57-74.
Halachev K, Bast H, Albrecht F, Lengauer T, Bock C: EpiExplorer: live exploration and global analysis of large epigenomic datasets. Genome Biol. 2012, 13: R96-
Sandve GK, Gundersen S, Rydbeck H, Glad IK, Holden L, Holden M, Liestol K, Clancy T, Ferkingstad E, Johansen M, Nygaard V, Tostesen E, Frigessi A, Hovig E: The genomic HyperBrowser: inferential genomics at the sequence level. Genome Biol. 2010–2011, 2010: 11-
Thomas-Chollier M, Hufton A, Heinig M, O’Keeffe S, Masri NE, Roider HG, Manke T, Vingron M: Transcription factor binding predictions using TRAP for the analysis of ChIP-seq data and regulatory SNPs. Nat Protoc. 1860–1869, 2011: 6-
Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, Kloos DU, Land S, Lewicki-Potapov B, Michael H, Munch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31: 374-378.
Ihaka R, Gentleman R: R: a language for data analysis and graphics. J Comput Graph Stat. 1996, 5: 299-314.
Chen LY, Eberlein M, Alsaaty S, Martinez-Anton A, Barb J, Munson PJ, Danner RL, Liu Y, Logun C, Shelhamer JH, Woszczek G: Cooperative and redundant signaling of leukotriene B4 and leukotriene D4 in human monocytes. Allergy. 2011, 66: 1304-1311.
Parmentier CN, Fuerst E, McDonald J, Bowen H, Lee TH, Pease JE, Woszczek G, Cousins DJ: Human T(H)2 cells respond to cysteinyl leukotrienes through selective expression of cysteinyl leukotriene receptor 1. J Allergy Clin Immunol. 2012, 129: 1136-1142.
Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, Zhang X, Bernstein BE, Nusbaum C, Jaffe DB, Gnirke A, Jaenisch R, Lander ES: Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008, 454: 766-770.
TriPrimate Methylome. [http://www2.cancer.ucl.ac.uk/medicalgenomics/triPrimateMethylome], 
O’Bleness M, Searles VB, Varki A, Gagneux P, Sikela JM: Evolution of genetic and genomic features unique to the human lineage. Nat Rev Genet. 2012, 13: 853-866.
Paten B, Herrero J, Fitzgerald S, Beal K, Flicek P, Holmes I, Birney E: Genome-wide nucleotide-level mammalian ancestor reconstruction. Genome Res. 1829–1843, 2008: 18-
Paten B, Herrero J, Beal K, Fitzgerald S, Birney E: Enredo and pecan: genome-wide mammalian consistency-based multiple alignment with paralogs. Genome Res. 1814–1828, 2008: 18-
Watanabe H, Fujiyama A, Hattori M, Taylor TD, Toyoda A, Kuroki Y, Noguchi H, BenKahla A, Lehrach H, Sudbrak R, Kube M, Taenzer S, Galgoczy P, Platzer M, Scharfe M, Nordsiek G, Blocker H, Hellmann I, Khaitovich P, Paabo S, Reinhardt R, Zheng HJ, Zhang XL, Zhu GF, Wang BF, Fu G, Ren SX, Zhao GP, Chen Z, Lee YS: DNA sequence and comparative analysis of chimpanzee chromosome 22. Nature. 2004, 429: 382-388.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19: 1639-1645.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB, Zhang X, Wang L, Issner R, Coyne M, Ku M, Durham T, Kellis M, Bernstein BE: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473: 43-49.
Hellman A, Chess A: Gene body-specific methylation on the active X chromosome. Science. 2007, 315: 1141-1143.
Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, Oberdoerffer P, Sandberg R, Oberdoerffer S: CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. 2011, 479: 74-79.
Jones PA: Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012, 13: 484-492.
Spitz F, Furlong EE: Transcription factors: from enhancer binding to developmental control. Nat Rev Genet. 2012, 13: 613-626.
Xu J, Pope SD, Jazirehi AR, Attema JL, Papathanasiou P, Watts JA, Zaret KS, Weissman IL, Smale ST: Pioneer factor interactions and unmethylated CpG dinucleotides mark silent tissue-specific enhancers in embryonic stem cells. Proc Natl Acad Sci U S A. 2007, 104: 12377-12382.
Ziller MJ, Gu H, Muller F, Donaghey J, Tsai LTY, Kohlbacher O, De Jager PL, Rosen ED, Bennett DA, Bernstein BE, Gnirke A, Meissner A: Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013, 500: 477-481.
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: 495-501.
The Center for Academic Research and Training in Anthropogeny: The Matrix of Comparative Anthropogeny (MOCA). [http://carta.anthropogeny.org/moca/topics/diverticulosis]
The National Human Genome Research Institute (NHGRI): A catalog of published genome-wide association studies. [http://www.genome.gov/gwastudies/], 
Bennett EA, Keller H, Mills RE, Schmidt S, Moran JV, Weichenrieder O, Devine SE: Active Alu retrotransposons in the human genome. Genome Res. 1875–1883, 2008: 18-
Lamprecht B, Walter K, Kreher S, Kumar R, Hummel M, Lenze D, Kochert K, Bouhlel MA, Richter J, Soler E, Stadhouders R, Johrens K, Wurster KD, Callen DF, Harte MF, Giefing M, Barlow R, Stein H, Anagnostopoulos I, Janz M, Cockerill PN, Siebert R, Dorken B, Bonifer C, Mathas S: Derepression of an endogenous long terminal repeat activates the CSF1R proto-oncogene in human lymphoma. Nat Med. 2010, 16: 571-579. 571p following 579
Jacques PE, Jeyakani J, Bourque G: The majority of primate-specific regulatory sequences are derived from transposable elements. PLoS Genet. 2013, 9: e1003504-
Ward MC, Wilson MD, Barbosa-Morais NL, Schmidt D, Stark R, Pan Q, Schwalie PC, Menon S, Lukk M, Watt S, Thybert D, Kutter C, Kirschner K, Flicek P, Blencowe BJ, Odom DT: Latent regulatory potential of human-specific repetitive elements. Mol Cell. 2012, 12: 00944-00946.
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.
Bock C, Beerman I, Lien WH, Smith ZD, Gu H, Boyle P, Gnirke A, Fuchs E, Rossi DJ, Meissner A: DNA methylation dynamics during in vivo differentiation of blood and skin stem cells. Mol Cell. 2012, 47: 633-647.
Wang H, Maurano MT, Qu H, Varley KE, Gertz J, Pauli F, Lee K, Canfield T, Weaver M, Sandstrom R, Thurman RE, Kaul R, Myers RM, Stamatoyannopoulos JA: Widespread plasticity in CTCF occupancy linked to DNA methylation. Genome Res. 2012, 22: 1680-1688.
Ernst J, Kellis M: Interplay between chromatin state, regulator binding, and regulatory motifs in six human cell types. Genome Res. 2013, 23: 1142-1154.
Brandeis M, Frank D, Keshet I, Siegfried Z, Mendelsohn M, Nemes A, Temper V, Razin A, Cedar H: Sp1 elements protect a CpG island from de novo methylation. Nature. 1994, 371: 435-438.
Boumber YA, Kondo Y, Chen X, Shen L, Guo Y, Tellez C, Estecio MR, Ahmed S, Issa JP: An Sp1/Sp3 binding polymorphism confers methylation protection. PLoS Genet. 2008, 4: e1000162-
Chahrour M, Jung SY, Shaw C, Zhou X, Wong ST, Qin J, Zoghbi HY: MeCP2, a key contributor to neurological disease, activates and represses transcription. Science. 2008, 320: 1224-1229.
McGraw CM, Samaco RC, Zoghbi HY: Adult neural function requires MeCP2. Science. 2011, 333: 186-
Siegfried Z, Eden S, Mendelsohn M, Feng X, Tsuberi BZ, Cedar H: DNA methylation represses transcription in vivo. Nat Genet. 1999, 22: 203-206.
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.
Tager AM, Luster AD: BLT1 and BLT2: the leukotriene B(4) receptors. Prostaglandins Leukot Essent Fatty Acids. 2003, 69: 123-134.
Li Y, Zhu J, Tian G, Li N, Li Q, Ye M, Zheng H, Yu J, Wu H, Sun J, Zhang H, Chen Q, Luo R, Chen M, He Y, Jin X, Zhang Q, Yu C, Zhou G, Huang Y, Cao H, Zhou X, Guo S, Hu X, Li X, Kristiansen K, Bolund L, Xu J, Wang W, Yang H, Wang J, Li R, Beck S, Zhang X: The DNA methylome of human peripheral blood mononuclear cells. PLoS Biol. 2010, 8: e1000533-
Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, Wiencke JK, Kelsey KT: DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinforma. 2012, 13: 86-10.1186/1471-2105-13-86.
Lowe R, Rakyan VK: Marmal-aid - a database for infinium HumanMethylation450. BMC Bioinforma. 2013, 14: 359-10.1186/1471-2105-14-359.
Kato K, Yokomizo T, Izumi T, Shimizu T: Cell-specific transcriptional regulation of human leukotriene B(4) receptor gene. J Exp Med. 2000, 192: 413-420.
Peters-Golden M, Henderson WR: Leukotrienes. N Engl J Med. 1841–1854, 2007: 357-
Hallstrand TS, Henderson WR: An update on the role of leukotrienes in asthma. Curr Opin Allergy Clin Immunol. 2010, 10: 60-66.
Kim ND, Chou RC, Seung E, Tager AM, Luster AD: A unique requirement for the leukotriene B4 receptor BLT1 for neutrophil recruitment in inflammatory arthritis. J Exp Med. 2006, 203: 829-835.
Jala VR, Haribabu B: Leukotrienes and atherosclerosis: new roles for old mediators. Trends Immunol. 2004, 25: 315-322.
Hartiala J, Li D, Conti DV, Vikman S, Patel Y, Tang WH, Brennan ML, Newman JW, Stephensen CB, Armstrong P, Hazen SL, Allayee H: Genetic contribution of the leukotriene pathway to coronary artery disease. Hum Genet. 2011, 129: 617-627.
Nancey S, Boschetti G, Hacini F, Sardi F, Durand PY, Le Borgne M, Furhmann L, Flourie B, Kaiserlian D: Blockade of LTB(4)/BLT(1) pathway improves CD8(+) T-cell-mediated colitis. Inflamm Bowel Dis. 2011, 17: 279-288.
Hicks A, Monkarsh SP, Hoffman AF, Goodnow R: Leukotriene B4 receptor antagonists as therapeutics for inflammatory disease: preclinical and clinical developments. Expert Opin Investig Drugs. 1909–1920, 2007: 16-
Owman C, Garzino-Demo A, Cocchi F, Popovic M, Sabirsh A, Gallo RC: The leukotriene B4 receptor functions as a novel type of coreceptor mediating entry of primary HIV-1 isolates into CD4-positive cells. Proc Natl Acad Sci U S A. 1998, 95: 9530-9534.
Lammermann T, Afonso PV, Angermann BR, Wang JM, Kastenmuller W, Parent CA, Germain RN: Neutrophil swarms require LTB4 and integrins at sites of cell death in vivo. Nature. 2013, 498: 371-375.
Tager AM, Bromley SK, Medoff BD, Islam SA, Bercury SD, Friedrich EB, Carafone AD, Gerszten RE, Luster AD: Leukotriene B4 receptor BLT1 mediates early effector T cell recruitment. Nat Immunol. 2003, 4: 982-990.
Goodarzi K, Goodarzi M, Tager AM, Luster AD, von Andrian UH: Leukotriene B4 and BLT1 control cytotoxic effector T cell recruitment to inflamed tissues. Nat Immunol. 2003, 4: 965-973.
Meyer M, Kircher M, Gansauge MT, Li H, Racimo F, Mallick S, Schraiber JG, Jay F, Prufer K, de Filippo C, Sudmant PH, Alkan C, Fu Q, Do R, Rohland N, Tandon A, Siebauer M, Green RE, Bryc K, Briggs AW, Stenzel U, Dabney J, Shendure J, Kitzman J, Hammer MF, Shunkov MV, Derevianko AP, Patterson N, Andres AM, Eichler EE: A high-coverage genome sequence from an archaic denisovan individual. Science. 2012, 338: 222-226.
Prufer K, Racimo F, Patterson N, Jay F, Sankararaman S, Sawyer S, Heinze A, Renaud G, Sudmant PH, de Filippo C, Li H, Mallick S, Dannemann M, Fu Q, Kircher M, Kuhlwilm M, Lachmann M, Meyer M, Ongyerth M, Siebauer M, Theunert C, Tandon A, Moorjani P, Pickrell J, Mullikin JC, Vohr SH, Green RE, Hellmann I, Johnson PLF, Blanche H: The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014, 505: 43-49.
Somel M, Liu X, Khaitovich P: Human brain evolution: transcripts, metabolites and their regulators. Nat Rev Neurosci. 2013, 14: 112-127.
Finch CE: Evolution in health and medicine sackler colloquium: evolution of the human lifespan and diseases of aging: roles of infection, inflammation, and nutrition. Proc Natl Acad Sci U S A. 2010, 107: 1718-1724.
Varki NM, Strobert E, Dick EJ, Benirschke K, Varki A: Biomedical differences between human and nonhuman hominids: potential roles for uniquely human aspects of sialic acid biology. Annu Rev Pathol. 2011, 6: 365-393.
Wang X, Mitra N, Secundino I, Banda K, Cruz P, Padler-Karavani V, Verhagen A, Reid C, Lari M, Rizzi E, Balsamo C, Corti G, De Bellis G, Longo L, Beggs W, Caramelli D, Tishkoff SA, Hayakawa T, Green ED, Mullikin JC, Nizet V, Bui J, Varki A: Specific inactivation of two immunomodulatory SIGLEC genes during human evolution. Proc Natl Acad Sci U S A. 2012, 109: 9935-9940.
Hodges E, Molaro A, Dos Santos CO, Thekkat P, Song Q, Uren PJ, Park J, Butler J, Rafii S, McCombie WR, Smith AD, Hannon GJ: Directional DNA methylation changes and complex intermediate States accompany lineage specificity in the adult hematopoietic compartment. Mol Cell. 2011, 44: 17-28.
Barreiro LB, Quintana-Murci L: From evolutionary genetics to human immunology: how selection shapes host defence genes. Nat Rev Genet. 2010, 11: 17-30.
Seok J, Warren HS, Cuenca AG, Mindrinos MN, Baker HV, Xu W, Richards DR, McDonald-Smith GP, Gao H, Hennessy L, Finnerty CC, Lopez CM, Honari S, Moore EE, Minei JP, Cuschieri J, Bankey PE, Johnson JL, Sperry J, Nathens AB, Billiar TR, West MA, Jeschke MG, Klein MB, Gamelli RL, Gibran NS, Brownstein BH, Miller-Graziano C, Calvano SE, Mason PH: Genomic responses in mouse models poorly mimic human inflammatory diseases. Proc Natl Acad Sci U S A. 2013, 110: 3507-3512.
Of men, not mice. Nat Med. 2013, 19: 379-379.
Martino DJ, Prescott SL: Progress in understanding the epigenetic basis for immune development, immune function, and the rising incidence of allergic disease. Curr Allergy Asthma Rep. 2013, 13: 85-92.
Suarez-Alvarez B, Rodriguez RM, Fraga MF, López-Larrea C: DNA methylation: a promising landscape for immune system-related diseases. Trends Genet. 2012, 28: 506-514.
Hernando-Herraez I, Prado-Martinez J, Garg P, Fernandez-Callejo M, Heyn H, Hvilsom C, Navarro A, Esteller M, Sharp AJ, Marques-Bonet T: Dynamics of DNA methylation in recent human and great ape evolution. PLoS Genet. 2013, 9: e1003763-
Joshi YB, Di Meco A, Pratico D: Modulation of amyloid-beta production by leukotriene B4 via the gamma-secretase pathway. JAD. 2013, 38: 503-506.
Osborn O, Olefsky JM: The cellular and signaling networks linking the immune system and metabolism in disease. Nat Med. 2012, 18: 363-374.
Islam SA, Thomas SY, Hess C, Medoff BD, Means TK, Brander C, Lilly CM, Tager AM, Luster AD: The leukotriene B4 lipid chemoattractant receptor BLT1 defines antigen-primed T cells in humans. Blood. 2006, 107: 444-453.
Soto PC, Stein LL, Hurtado-Ziola N, Hedrick SM, Varki A: Relative over-reactivity of human versus chimpanzee lymphocytes: implications for the human diseases associated with immune activation. J Immunol. 2010, 184: 4185-4195.
Petryszak R, Burdett T, Fiorelli B, Fonseca NA, Gonzalez-Porta M, Hastings E, Huber W, Jupp S, Keays M, Kryvych N, McMurry J, Marioni JC, Malone J, Megy K, Rustici G, Tang AY, Taubert J, Williams E, Mannion O, Parkinson HE, Brazma A: Expression Atlas update - a database of gene and transcript expression from microarray- and sequencing-based functional genomics experiments. Nucleic Acids Res. 2013, 42: D926-932.
Colecchia D, Strambi A, Sanzone S, Iavarone C, Rossi M, Dall’Armi C, Piccioni F, Chiariello M, Verrotti Di Pianella A: MAPK15/ERK8 stimulates autophagy by interacting with LC3 and GABARAP proteins. Autophagy. 2012, 8: 1724-1740.
Bienvenu T, Chelly J: Molecular genetics of rett syndrome: when DNA methylation goes unrecognized. Nat Rev Genet. 2006, 7: 415-426.
We thank the staff of the Tierpark Nordhorn, Basel Zoo, Leipzig Zoo and the German Primate Center for providing blood samples of chimpanzees and rhesus macaques. This work was supported by the Wellcome Trust (99148), EU-FP7 projects EPIGENESYS (257082) and BLUEPRINT (282510) and a Royal Society Wolfson Research Merit Award (WM100023) to SB. The research based at Guy’s and St Thomas’ NHS Foundation Trust and King’s College London (GW, HRF) was supported by the National Institute for Health Research (NIHR) Biomedical Research Centre. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health.
The authors declare that they have no competing interests.
Conceived and designed the experiments: CGB, SB, GW. Performed the experiments: CGB, LMB, HRF, GW. Analyzed the data: CGB, GAW, GW. Contributed reagents/materials/analysis tools: CR, LW, GW, CGB, GAW, AF. Wrote the paper: CGB. All authors reviewed and approved the manuscript.
Electronic supplementary material
Additional file 2: Figure S1:Alu distribution of human s-DMRs. Hypermethylated s-DMRs are depleted within both the most ancient AluJ and youngest AluY categories, but were enriched within the second oldest AluS set (χ2P < 2.2 × 10-16), which still possesses mobilization ability . Hypomethylated s-DMRs also show this pattern, with an increase in AluS (χ2P < 2.2 × 10-16), but also an increase in the youngest and most active AluY (χ2P = 1.15 × 10-4). (TIFF 13 MB)
Additional file 7: Figure S3: Change in transcription factor motif binding prediction within s-DMRs between primates calculated via TRAP  for all 904 TRANSFAC motifs  within all individual hypomethylated s-DMR and hypermethylated s-DMR CpGi regions. The red lines indicate a magnitude increase or decrease in predicted motif binding (-log10P-value). In (A) the red circle in human out-lying positive motifs indicates the LTB4R - RFX (EFC_ Q6) motif. (TIFF 3 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Wilson, G.A., Butcher, L.M., Foster, H.R. et al. Human-specific epigenetic variation in the immunological Leukotriene B4 Receptor (LTB4R/BLT1) implicated in common inflammatory diseases. Genome Med 6, 19 (2014) doi:10.1186/gm536
- CTCF Binding
- Gene Body Methylation
- TFBS Motif
- CTCF Motif
- Calcium Mobilization Assay