Skip to main content

Human-specific epigenetic variation in the immunological Leukotriene B4 Receptor (LTB4R/BLT1) implicated in common inflammatory diseases



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 [3]. 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 [6] with close relatives such as the chimpanzee [7], rhesus macaque [8] and other primates [912] have explored these alterations at the sequence level. Accelerated regions have highlighted potentially modified neurological [13, 14] and anatomic [15] pathways, human-specific duplications have been implicated in enhanced neuronal migration [16] and human-specific deletions have identified lost regulatory regions [17]. Comparative epigenomics has only recently started to be explored in a number of tissues [1822]. 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 [25] 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 [26]. 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) [27]; or associated with CpG density changes, due to gain or loss of species-specific CpGs [28], impacting on CpG island (CpGi) strength [2931]. The major proportional components of peripheral blood cells are similar between all three species [3234]; 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 [35].

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)) [37] will have then left not only a genetic [38] 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 [39].


DNA methylomes

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 [40]. 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 [32], has very similar ratios in chimpanzee (male neutrophils 59%, lymphocytes 28%; female neutrophils 49%, lymphocytes 41% [33]) and rhesus macaque (neutrophils approximately 67% and lymphocytes 26 to 35% [34]). 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 [41] 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) [42]. Sequence quality control was performed using FASTQC [43]. MeDUSA utilized BWA (v0.5.8) [44] 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 [45] 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) [46]. 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 [47] and BayesPeak v1.8 [48]. 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 [49] with minMatch = 0.7. This parameter is appropriate due the recent shared evolutionary history of these species [50]. Subsequently, human-specific differentially methylated regions (s-DMRs) were located through a series of intersection analyses utilizing the BEDTools [51] 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 [52].

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 [53]. Significance was calculated by Genomic Hyperbrowser [54] 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 [53] 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 [51] 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 [49] with the BEDTools [51] FastaFromBed function. The transcription factor affinity prediction tool TRAP [55] (multiple sequences) was implemented with the 904 TRANSFAC [56] 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 [55] 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.

Statistical analysis

All other statistical analyses were performed in the R environment [57]. 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. [19] 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. [21] 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.

Real-time PCR

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 [40]. 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 [60] Automated MeDIP-seq was performed as previously described [41], and data were processed using the MeDUSA pipeline [42].

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 [61]. The chimpanzee and macaque results can also be viewed on the human sequence having been converted using the liftOver utility [49] 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 [47] peak-calling algorithm, but with the additional conservative requirement that regions must also have an overlap with a BayesPeak v1.8 [48] 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 [49] 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% [62]) 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 [65] 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[66].

Figure 1
figure 1

Genome-wide view of human s-DMRs. Pooled DNA from uncultured whole blood cell samples including both sexes (all 60% male) were analyzed for each species. Methylated fragments via MeDIP-seq were aligned to the appropriate species’ genome; then peaks were called within these, using two peak-calling algorithms (MACS v1.4.1 [47] and BayesPeak v1.8 [48]). Triangulation reciprocal LiftOver [49] comparison then identified human s-DMRs. Outer ring, human chromosomes; blue ring, hypermethylated human s-DMRs (15,858); yellow ring, hypomethylated human s-DMRs (22,758); green ring, CNV density; inner ring, gene density.

Genome feature annotation enrichment

We compared the location of human s-DMRs to genomic functional annotation, as defined by the ChromHMM segmentation analysis [67], using the HapMap line interrogated by Encode (B-lymphocyte GM12878 data). We calculated significance via Epiexplorer [53] and the Genomic Hyperbrowser [54]. 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 [36]. Increased levels of s-DMRs in these transcription-associated regions may be indicative of the potential role methylation variability plays in gene body methylation [68] 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. [73] (data not shown, χ2P < 2.2 × 10-16).

Figure 2
figure 2

Proportional genomic annotation coverage of s-DMRs compared with HapMap B-lymphocyte EBV GM12878 ChromHMM data[68]. (A) Hypomethylated s-DMRs; (B) hypermethylated s-DMRs. Various genomic annotations were significantly enriched, as defined by this segmentation analysis, calculated via the Genomic Hyperbrowser [54] (P-value overlap MC * < 0.05, ** < 0.005).

Gene enrichment

Gene enrichment analysis of these s-DMR regions was performed using the Genomic Regions Enrichment of Annotations Tool (GREAT) [74]. 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 [75]. 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 [76] are enriched in these s-DMR regions by approximately 1.40× compared to the genome average (data not shown; χ2P = 5.09 × 10-3).

Repetitive elements

The hypermethylated s-DMR set was found to be increased within repetitive elements, specifically the SINE group (Figure 3; via Epiexplorer [53]). 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 [77], 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 AluS2P < 2.2 × 10-16), but also showed a slight increase in the youngest and most active AluY2P = 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 [78]. Primate-specific open chromatin regions (DNAse I hypersensitivity sites) are enriched for transposable elements (approximately 63%) [79] and methylation changes in trans-species experiments have also indicated their regulatory potential [80].

Figure 3
figure 3

Subcategorization of repeat element increase in hypermethylated s-DMRs via Epiexplorer in comparison with a reshuffled control set[53](medium overlap ≥10%). Increased hyper s-DMRs within the SINE group are identified, which comprises predominately Alu s. LINE, long interspersed nucleotide element; LTR, long terminal repeat; RC, rolling circle; snRNA, short nuclear RNA; srpRNA, signal recognition particle RNA.

CpG islands

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).

Figure 4
figure 4

Human s-DMRs in CpG islands versus CpG island shores. Number of hypo- and hypermethylated s-DMRs for island (77 and 45) and shore (821 and 431), respectively, is corrected for proportion of genome size (genomic space for islands = 23.8 Mb and for shores = 89.5 Mb). Directly comparing island and shore s-DMRs with regard to possible locations of hypomethylated s-DMRs (all co-locating chimpanzee and rhesus macaque peaks in these regions), χ2P = 2.90 × 10-29, and hypermethylated s-DMRs (from all location of human peaks in these regions), χ2P = 1.93 × 10-8; combined χ2P = 1.80 × 10-32. Therefore, human-specific peaks are more likely than non-human-specific peaks to reside within CpGi shores.

Sequence similarity of s-DMRs in CpG islands

Sequence divergence could be contributing to the s-DMRs within these CpGis [27]. 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.

Figure 5
figure 5

Sequence divergence between human and chimpanzee within peak regions. There was no significant difference (Wilcoxon P > 0.05 for both) between s-DMRs in islands and all chimpanzee and rhesus macaque (C&M) combined peaks or all human peaks in islands.

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 [27]. 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 [55] and calculated the binding prediction of 904 TFBSs (TRANSFAC vertebrate V$motifs [56]) 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 [27]. This trans-species effect is consistent with the disruption of CTCF occupancy being linked to increased methylation, recently observed across numerous human cell types [83], as well as the different role of CTCF binding regions within open and closed chromatin [84]. 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.

Figure 6
figure 6

Change in transcription factor motif binding prediction within s-DMRs between primates calculated via TRAP [[55]] with TRANSFAC motifs[56]. (A,B) Difference in binding prediction (total corrected Benjamini-Hochberg -log10P-value) between human and chimpanzee (y-axis) and human and rhesus macaque (x-axis) for each motif within the total set of hypomethylated (A) and hypermethylated (B) CpGi s-DMRs. Known TFBSs with MDR effects are highlighted in color (SP1 in red, CTCF in blue, RFX motif family in green). Both CTCF motifs show a consistent increase in the hypomethylated s-DMRs, as well as a consistent decrease in the hypermethylated DMRs, with respect to human. The MeCP2 motif is identified as a strongly increased outlier in the hypermethylated s-DMRs (orange).

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 [89].

As highly expressed genes are associated with low gene promoter methylation and high gene body methylation [90], 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) [91] 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).

Figure 7
figure 7

Comparative DNA methylation of LTB4R visualized in the UCSC browser. Human hypomethylated s-DMRs (yellow) are shown in the promoter CpGi (CpG: 99) of LTB4R (major isoform LTB4R-001 outlined in red). Methylation scale is in reads per millions (RPM) for each species from MeDIP-seq (human, light blue; chimpanzee, orange; rhesus macaque, olive green). As well as reduced promoter methylation, larger gene-body methylation, which is related to higher expression [69], was also seen in human compared with the other species over the sole exon (approximately 1.29-fold stronger peak MAC P-value over gene body CpGi (CpG:76)). In this complex locus the promoter of the major isoform of LTB4R (highlighted with a red rectangle) also co-locates with the gene body of the low-specificity receptor LTB4R2 and CIDEB. LTB4R has strong expression in all blood subtypes, particularly the myeloid lineage, including monocytes (Additional file 6).

Table 1 Hypomethylated 5′ human promoter with primate differential exon methylation

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. [19] 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. [92], 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.

Figure 8
figure 8

Results of HpaII digestion of LTB4R followed by sequencing for promoter and gene body CpGs in human (Hs) and chimpanzee (Pt) (Martin et al .[19]via GEO). Methylation within the (A) promoter (Prom = CpG:99 in Figure 7) and (B) gene body (Body = CpG:76 in Figure 7) CpG islands of LTB4R from reverse scores (1 - p(U)) of all included HpaII MethylSeq sites analyzed by MetMap6 in purified neutrophils. These data replicate the significant difference identified in MeDIP peripheral blood. Four human and four chimpanzee samples had average methylation of 18.4% and 70.3% in the promoter (Prom), and 75.9% and 62.0% in the exonic gene body (Body) orthologous CpGis (Wilcoxon P < 0.05 for both).

Primate comparative chromatin data were available from B-cell lymphoblastoid cell lines [21]. 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.

Figure 9
figure 9

Human-specific H3K4me3 enrichment in the LTB4R promoter from ChIP-seq data derived from B cells (lymphoblastoid cell line; data from Cain et al. [[21]] via GEO). These data identified significant activating peaks in all three human samples, but the signal was not strong enough for any peaks to be called in all three chimpanzee and all three rhesus macaque samples analyzed.

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. [93]. 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 [94] 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).

Figure 10
figure 10

Data from the MARMAL-AID Human 450k Methylation array repository[94]for the LTB4R promoter and gene body CpG islands. These data are derived from 1,665 whole blood samples from healthy and non-cancer disease subjects from multiple experiments. This also showed a consistent low average level of human LTB4R promoter methylation (mean 29.4%, standard deviation 7.4%) and high average LTB4R gene body methylation (mean 87.9%, standard deviation 2.3%). Human LTB4R methylation whole blood 1,665 samples.

Methylation of this LTB4R promoter has previously been established by Kato et al. [95] 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 [96], including asthma [97], inflammatory arthritis [98], atherosclerosis [99, 100], inflammatory bowel disease [101], and psoriasis [102]. It has also previously been postulated to play a role in HIV infection [103]. 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 [104], and early effector T-cell recruitment [105], acting as a potent non-chemokine pathway for cytotoxic effector cell traffic [106]. 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 [27], although whether one motif would be expected to lead to such a considerable difference in methylation is not certain. Lienert et al. [27] 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).

Figure 11
figure 11

Expression and signaling of LTB4R (BLT1) in PBMCs isolated from human and rhesus macaque peripheral blood. (A) Real-time RT-PCR analysis of LTB4R mRNA in isolated human and rhesus macaque PBMCs normalized to the reference gene (18s); results from seven different human donors and from pooled rhesus macaque PBMCs run as four separate experiments; mean ± standard error of the mean. (B) Human and rhesus macaque PBMCs were stimulated with indicated concentrations of LTB4 and intracellular calcium mobilizations were recorded. Results are expressed as the ratio of stimulated over basal (S/B) peak calcium fluxes obtained from three different human donors and pooled rhesus macaque PBMCs run as four experiments analyzed simultaneously, mean ± standard error of the mean. (C) Human PBMCs were pre-incubated for 10 minutes with different concentrations of BLT1 inhibitor LY22398 and calcium mobilization was analyzed in response to LTB4 (300 nmol/L). Data are presented as mean ± standard error of the mean percentages of maximum response to LTB4 (N = 3).

LTB4 signaling through the LTB4R/BLT1 receptor induces calcium mobilization in cells known to express the receptor [58]. 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 [110]. 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 [111]. 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 [52] and genetic loss of human-specific regulatory DNA has previously highlighted regulation change as a human evolutionary mechanism [17].

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 [114]. DNA binding proteins localizing within these boundary regions, or shores, may help define these methylation states [18]. Thus, these species-specific modifications may be acting in similar fashion.

Pathogens have imposed a strong selective influence on the human genome [115]. 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 [118]. 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 [119].

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 [17]. Furthermore, a recently published comparative array analysis in peripheral blood from Hernando-Herraez et al. [119], 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) [97], atherosclerosis, Alzheimer’s disease [120], obesity-related inflammation [121], and inflammatory bowel disease [101]. LTB4R is involved in amplifying T-cell recruitment [122] and a heightened T-cell response in human compared with chimpanzee has previously been identified [123]. 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 [124] and furthermore that LTB4R expression on T cells is regulated by inflammation and may link innate and adaptive immune responses [122]. 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 [125], and the MECP2 gene, which is mutated in Rett syndrome [126].

The effect of genetic mutation in particular TFBSs on the methylation state of CpGis has recently been refined by Lienert et al. [27]. 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.



base pair


chromatin immunoprecipitation


CpG island


differentially methylated region


false discovery rate


Gene Expression Omnibus


leukocyte subtype-related differentially methylated CpG position


Monte Carlo


methylation determining region


methylated DNA immunoprecipitation


million years ago


peripheral blood mononuclear cell


polymerase chain reaction


species-specific DMR


short interspersed nucleotide element


single-nucleotide polymorphism


transcription factor binding site


University of California Santa Cruz.


  1. Eder W, Ege MJ, von Mutius E: The asthma epidemic. N Engl J Med. 2006, 355: 2226-2235.

    CAS  PubMed  Google Scholar 

  2. 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.

    CAS  PubMed  Google Scholar 

  3. 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.

    CAS  PubMed  Google Scholar 

  4. Feinberg AP: Phenotypic plasticity and the epigenetics of human disease. Nature. 2007, 447: 433-440.

    CAS  PubMed  Google Scholar 

  5. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. 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.

    CAS  PubMed  Google Scholar 

  7. The Chimpanzee Sequencing and Analysis Consortium: Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005, 437: 69-87.

    Google Scholar 

  8. 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.

    CAS  PubMed  Google Scholar 

  9. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  10. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. 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.

    PubMed Central  PubMed  Google Scholar 

  12. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  13. 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.

    CAS  PubMed  Google Scholar 

  14. Prabhakar S, Noonan JP, Paabo S, Rubin EM: Accelerated evolution of conserved noncoding sequences in humans. Science. 2006, 314: 786-

    CAS  PubMed  Google Scholar 

  15. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  16. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  17. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  18. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  19. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  20. 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-

    PubMed Central  CAS  PubMed  Google Scholar 

  21. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  22. 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-

    PubMed Central  CAS  PubMed  Google Scholar 

  23. Feil R, Fraga MF: Epigenetics and the environment: emerging patterns and implications. Nat Rev Genet. 2011, 13: 97-109.

    Google Scholar 

  24. Eichten S, Borevitz J: Epigenomics: Methylation’s mark on inheritance. Nature. 2013, 495: 181-182.

    CAS  PubMed  Google Scholar 

  25. Harris RA, Rogers J, Milosavljevic A: Human-specific changes of genome structure detected by genomic triangulation. Science. 2007, 316: 235-237.

    CAS  PubMed  Google Scholar 

  26. Richards EJ: Inherited epigenetic variation–revisiting soft inheritance. Nat Rev Genet. 2006, 7: 395-401.

    CAS  PubMed  Google Scholar 

  27. 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.

    CAS  PubMed  Google Scholar 

  28. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  29. 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.

    CAS  PubMed  Google Scholar 

  30. Bock C, Walter J, Paulsen M, Lengauer T: CpG island mapping by epigenome prediction. PLoS Comput Biol. 2007, 3: e110-

    PubMed Central  PubMed  Google Scholar 

  31. Cohen NM, Kenigsberg E, Tanay A: Primate CpG islands are maintained by heterogeneous evolutionary regimes involving minimal selection. Cell. 2011, 145: 773-786.

    CAS  PubMed  Google Scholar 

  32. 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.

    CAS  PubMed  Google Scholar 

  33. 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.

    CAS  PubMed  Google Scholar 

  34. 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.

    CAS  PubMed  Google Scholar 

  35. Watson JA, Watson CJ, McCann A, Baugh J: Epigenetics, the epicenter of the hypoxic response. Epigenetics. 2010, 5: 293-296.

    CAS  PubMed  Google Scholar 

  36. Bock C: Analysing and interpreting DNA methylation data. Nat Rev Genet. 2012, 13: 705-719.

    CAS  PubMed  Google Scholar 

  37. Olson MV, Varki A: Sequencing the chimpanzee genome: insights into human evolution and disease. Nat Rev Genet. 2003, 4: 20-28.

    CAS  PubMed  Google Scholar 

  38. 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-

    PubMed Central  PubMed  Google Scholar 

  39. Varki A, Geschwind DH, Eichler EE: Explaining human uniqueness: genome interactions with environment, behaviour and culture. Nat Rev Genet. 2008, 9: 749-763.

    PubMed Central  CAS  PubMed  Google Scholar 

  40. The American Society of Primatologists: Principles for the ethical treatment of non-human primates. []

    Google Scholar 

  41. Butcher LM, Beck S: AutoMeDIP-seq: a high-throughput, whole genome, DNA methylation assay. Methods. 2010, 52: 223-231.

    PubMed Central  CAS  PubMed  Google Scholar 

  42. 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.

    CAS  Google Scholar 

  43. FastQC. [], []

  44. Li H, Durbin R: Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009, 25: 1754-1760.

    PubMed Central  CAS  PubMed  Google Scholar 

  45. 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-

    Google Scholar 

  46. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. 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-

    PubMed Central  PubMed  Google Scholar 

  48. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  49. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  50. Bardet AF, He Q, Zeitlinger J, Stark A: A computational pipeline for comparative ChIP-seq analyses. Nat Protoc. 2012, 7: 45-61.

    CAS  Google Scholar 

  51. Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26: 841-842.

    PubMed Central  CAS  PubMed  Google Scholar 

  52. 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.

    Google Scholar 

  53. 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-

    PubMed Central  PubMed  Google Scholar 

  54. 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-

    Google Scholar 

  55. 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-

    Google Scholar 

  56. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  57. Ihaka R, Gentleman R: R: a language for data analysis and graphics. J Comput Graph Stat. 1996, 5: 299-314.

    Google Scholar 

  58. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  59. 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.

    CAS  PubMed  Google Scholar 

  60. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  61. TriPrimate Methylome. [], []

  62. 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.

    PubMed Central  PubMed  Google Scholar 

  63. 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-

    Google Scholar 

  64. 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-

    Google Scholar 

  65. 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.

    CAS  PubMed  Google Scholar 

  66. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  67. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  68. Hellman A, Chess A: Gene body-specific methylation on the active X chromosome. Science. 2007, 315: 1141-1143.

    CAS  PubMed  Google Scholar 

  69. 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.

    CAS  PubMed  Google Scholar 

  70. Jones PA: Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012, 13: 484-492.

    CAS  PubMed  Google Scholar 

  71. Spitz F, Furlong EE: Transcription factors: from enhancer binding to developmental control. Nat Rev Genet. 2012, 13: 613-626.

    CAS  PubMed  Google Scholar 

  72. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  73. 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.

    CAS  PubMed  Google Scholar 

  74. 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.

    CAS  PubMed  Google Scholar 

  75. The Center for Academic Research and Training in Anthropogeny: The Matrix of Comparative Anthropogeny (MOCA). []

    Google Scholar 

  76. The National Human Genome Research Institute (NHGRI): A catalog of published genome-wide association studies. [], []

    Google Scholar 

  77. 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-

    Google Scholar 

  78. 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

    CAS  PubMed  Google Scholar 

  79. Jacques PE, Jeyakani J, Bourque G: The majority of primate-specific regulatory sequences are derived from transposable elements. PLoS Genet. 2013, 9: e1003504-

    PubMed Central  CAS  PubMed  Google Scholar 

  80. 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.

    Google Scholar 

  81. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  82. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  83. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  84. Ernst J, Kellis M: Interplay between chromatin state, regulator binding, and regulatory motifs in six human cell types. Genome Res. 2013, 23: 1142-1154.

    PubMed Central  CAS  PubMed  Google Scholar 

  85. 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.

    CAS  PubMed  Google Scholar 

  86. 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-

    PubMed Central  PubMed  Google Scholar 

  87. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  88. McGraw CM, Samaco RC, Zoghbi HY: Adult neural function requires MeCP2. Science. 2011, 333: 186-

    PubMed Central  CAS  PubMed  Google Scholar 

  89. Siegfried Z, Eden S, Mendelsohn M, Feng X, Tsuberi BZ, Cedar H: DNA methylation represses transcription in vivo. Nat Genet. 1999, 22: 203-206.

    CAS  PubMed  Google Scholar 

  90. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  91. Tager AM, Luster AD: BLT1 and BLT2: the leukotriene B(4) receptors. Prostaglandins Leukot Essent Fatty Acids. 2003, 69: 123-134.

    CAS  PubMed  Google Scholar 

  92. 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-

    PubMed Central  PubMed  Google Scholar 

  93. 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.

    Google Scholar 

  94. Lowe R, Rakyan VK: Marmal-aid - a database for infinium HumanMethylation450. BMC Bioinforma. 2013, 14: 359-10.1186/1471-2105-14-359.

    Google Scholar 

  95. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  96. Peters-Golden M, Henderson WR: Leukotrienes. N Engl J Med. 1841–1854, 2007: 357-

    Google Scholar 

  97. Hallstrand TS, Henderson WR: An update on the role of leukotrienes in asthma. Curr Opin Allergy Clin Immunol. 2010, 10: 60-66.

    PubMed Central  CAS  PubMed  Google Scholar 

  98. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  99. Jala VR, Haribabu B: Leukotrienes and atherosclerosis: new roles for old mediators. Trends Immunol. 2004, 25: 315-322.

    CAS  PubMed  Google Scholar 

  100. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  101. 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.

    PubMed  Google Scholar 

  102. 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-

    Google Scholar 

  103. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  104. 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.

    PubMed  Google Scholar 

  105. 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.

    CAS  PubMed  Google Scholar 

  106. 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.

    CAS  PubMed  Google Scholar 

  107. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  108. 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.

    PubMed Central  PubMed  Google Scholar 

  109. Somel M, Liu X, Khaitovich P: Human brain evolution: transcripts, metabolites and their regulators. Nat Rev Neurosci. 2013, 14: 112-127.

    CAS  PubMed  Google Scholar 

  110. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  111. 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.

    CAS  PubMed  Google Scholar 

  112. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  113. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  114. Barreiro LB, Quintana-Murci L: From evolutionary genetics to human immunology: how selection shapes host defence genes. Nat Rev Genet. 2010, 11: 17-30.

    CAS  PubMed  Google Scholar 

  115. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  116. Of men, not mice. Nat Med. 2013, 19: 379-379.

    Google Scholar 

  117. 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.

    CAS  PubMed  Google Scholar 

  118. 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.

    CAS  PubMed  Google Scholar 

  119. 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-

    PubMed Central  CAS  PubMed  Google Scholar 

  120. 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.

    Google Scholar 

  121. Osborn O, Olefsky JM: The cellular and signaling networks linking the immune system and metabolism in disease. Nat Med. 2012, 18: 363-374.

    CAS  PubMed  Google Scholar 

  122. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  123. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  124. 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.

    PubMed Central  PubMed  Google Scholar 

  125. 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.

    PubMed Central  CAS  PubMed  Google Scholar 

  126. Bienvenu T, Chelly J: Molecular genetics of rett syndrome: when DNA methylation goes unrecognized. Nat Rev Genet. 2006, 7: 415-426.

    CAS  PubMed  Google Scholar 

Download references


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.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christopher G Bell.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

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 1: Table S1: Genomic Regions Enrichment of Annotations Tool (GREAT) s-DMR binomial analysis. (XLSX 113 KB)


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 [77]. Hypomethylated s-DMRs also show this pattern, with an increase in AluS2P < 2.2 × 10-16), but also an increase in the youngest and most active AluY2P = 1.15 × 10-4). (TIFF 13 MB)

Additional file 3: Table S2: CpG island and CpG island shore s-DMRs. (XLSX 45 KB)

Additional file 4: Table S3: Transcription Factor Affinity Prediction (TRAP) motif results across species. (XLSX 18 KB)

Additional file 5: Table S4: Gene body methylation of 5′ promoter s-DMR genes. (XLSX 64 KB)

Additional file 6: Figure S2: BioGPS - human LTB4R expression profile. (TIFF 1 MB)


Additional file 7: Figure S3: Change in transcription factor motif binding prediction within s-DMRs between primates calculated via TRAP [55] for all 904 TRANSFAC motifs [56] 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

Rights and permissions

Reprints and permissions

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).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: