Abstract
Although genome sequencing has identified numerous noncoding alterations between primate species, which of those are regulatory and potentially relevant to the evolution of the human brain is unclear. Here we annotated cis-regulatory elements (CREs) in the human, rhesus macaque and chimpanzee genomes using chromatin immunoprecipitation followed by sequencing (ChIP-seq) in different anatomical regions of the adult brain. We found high similarity in the genomic positioning of rhesus macaque and human CREs, suggesting that the majority of these elements were already present in a common ancestor 25 million years ago. Most of the observed regulatory changes between humans and rhesus macaques occurred before the ancestral separation of humans and chimpanzees, leaving a modest set of regulatory elements with predicted human specificity. Our data refine previous predictions and hypotheses on the consequences of genomic changes between primate species and allow the identification of regulatory alterations relevant to the evolution of the brain.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
$209.00 per year
only $17.42 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Accession codes
References
Azevedo, F.A. et al. Equal numbers of neuronal and nonneuronal cells make the human brain an isometrically scaled-up primate brain. J. Comp. Neurol. 513, 532–541 (2009).
Rakic, P. Evolution of the neocortex: a perspective from developmental biology. Nat. Rev. Neurosci. 10, 724–735 (2009).
Geschwind, D.H. & Rakic, P. Cortical evolution: judge the brain by its cover. Neuron 80, 633–647 (2013).
Rogers, J. & Gibbs, R.A. Comparative primate genomics: emerging patterns of genome content and dynamics. Nat. Rev. Genet. 15, 347–359 (2014).
O'Bleness, M., Searles, V.B., Varki, A., Gagneux, P. & Sikela, J.M. Evolution of genetic and genomic features unique to the human lineage. Nat. Rev. Genet. 13, 853–866 (2012).
King, M.C. & Wilson, A.C. Evolution at two levels in humans and chimpanzees. Science 188, 107–116 (1975).
Haygood, R., Babbitt, C.C., Fedrigo, O. & Wray, G.A. Contrasts between adaptive coding and noncoding changes during human evolution. Proc. Natl. Acad. Sci. USA 107, 7853–7857 (2010).
Konopka, G. et al. Human-specific transcriptional networks in the brain. Neuron 75, 601–617 (2012).
Shlyueva, D., Stampfel, G. & Stark, A. Transcriptional enhancers: from properties to genome-wide predictions. Nat. Rev. Genet. 15, 272–286 (2014).
Hong, J.W., Hendrix, D.A. & Levine, M.S. Shadow enhancers as a source of evolutionary novelty. Science 321, 1314 (2008).
Perry, M.W., Boettiger, A.N., Bothma, J.P. & Levine, M. Shadow enhancers foster robustness of Drosophila gastrulation. Curr. Biol. 20, 1562–1567 (2010).
Wittkopp, P.J. & Kalay, G. Cis-regulatory elements: molecular mechanisms and evolutionary processes underlying divergence. Nat. Rev. Genet. 13, 59–69 (2012).
Rubinstein, M. & de Souza, F.S. Evolution of transcriptional enhancers and animal diversity. Phil. Trans. R. Soc. Lond. B 368, 20130017 (2013).
Wray, G.A. The evolutionary significance of cis-regulatory mutations. Nat. Rev. Genet. 8, 206–216 (2007).
Schmidt, D. et al. Five-vertebrate ChIP-seq reveals the evolutionary dynamics of transcription factor binding. Science 328, 1036–1040 (2010).
Villar, D. et al. Enhancer evolution across 20 mammalian species. Cell 160, 554–566 (2015).
Bulger, M. & Groudine, M. Functional and mechanistic diversity of distal transcription enhancers. Cell 144, 327–339 (2011).
Villar, D., Flicek, P. & Odom, D.T. Evolution of transcription factor binding in metazoans-mechanisms and functional implications. Nat. Rev. Genet. 15, 221–233 (2014).
Cotney, J. et al. The evolution of lineage-specific regulatory activities in the human embryonic limb. Cell 154, 185–196 (2013).
Reilly, S.K. et al. Evolutionary genomics. Evolutionary changes in promoter and enhancer activity during human corticogenesis. Science 347, 1155–1159 (2015).
Visel, A. et al. A high-resolution enhancer atlas of the developing telencephalon. Cell 152, 895–908 (2013).
Creyghton, M.P. et al. Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc. Natl. Acad. Sci. USA 107, 21931–21936 (2010).
Prescott, S.L. et al. Enhancer divergence and cis-regulatory evolution in the human and chimp neural crest. Cell 163, 68–83 (2015).
Kundaje, A. et al. Integrative analysis of 111 reference human epigenomes. Nature 518, 317–330 (2015).
Rada-Iglesias, A. et al. A unique chromatin signature uncovers early developmental enhancers in humans. Nature 470, 279–283 (2011).
Vermunt, M.W. et al. Large-scale identification of coregulated enhancer networks in the adult human brain. Cell Rep. 9, 767–779 (2014).
Kumar, S. & Hedges, S.B. A molecular timescale for vertebrate evolution. Nature 392, 917–920 (1998).
Ernst, J. & Kellis, M. ChromHMM: automating chromatin-state discovery and characterization. Nat. Methods 9, 215–216 (2012).
Ernst, J. & Kellis, M. Discovery and characterization of chromatin states for systematic annotation of the human genome. Nat. Biotechnol. 28, 817–825 (2010).
Ernst, J. & Kellis, M. Large-scale imputation of epigenomic datasets for systematic annotation of diverse human tissues. Nat. Biotechnol. 33, 364–376 (2015).
Gibbs, R.A. et al. Evolutionary and biomedical insights from the rhesus macaque genome. Science 316, 222–234 (2007).
Chimpanzee Sequencing and Analysis Consortium. Initial sequence of the chimpanzee genome and comparison with the human genome. Nature 437, 69–87 (2005).
Brawand, D. et al. The evolution of gene expression levels in mammalian organs. Nature 478, 343–348 (2011).
Wang, Z. et al. Combinatorial patterns of histone acetylations and methylations in the human genome. Nat. Genet. 40, 897–903 (2008).
Bardet, A.F., He, Q., Zeitlinger, J. & Stark, A. A computational pipeline for comparative ChIP-seq analyses. Nat. Protoc. 7, 45–61 (2012).
Love, M.I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Hare, E.E., Peterson, B.K., Iyer, V.N., Meier, R. & Eisen, M.B. Sepsid even-skipped enhancers are functionally conserved in Drosophila despite lack of sequence conservation. PLoS Genet. 4, e1000106 (2008).
Fisher, S., Grice, E.A., Vinton, R.M., Bessling, S.L. & McCallion, A.S. Conservation of RET regulatory function from human to zebrafish without sequence similarity. Science 312, 276–279 (2006).
Andersen, B.B., Korbo, L. & Pakkenberg, B. A quantitative study of the human cerebellum with unbiased stereological techniques. J. Comp. Neurol. 326, 549–560 (1992).
Arnold, C.D. et al. Quantitative genome-wide enhancer activity maps for five Drosophila species show functional enhancer conservation and turnover during cis-regulatory evolution. Nat. Genet. 46, 685–692 (2014).
Vierstra, J. et al. Mouse regulatory DNA landscapes reveal global principles of cis-regulatory evolution. Science 346, 1007–1012 (2014).
Rahl, P.B. & Young, R.A. MYC and transcription elongation. Cold Spring Harb. Perspect. Med. 4, a020990 (2014).
Lin, C.Y. et al. Transcriptional amplification in tumor cells with elevated c-Myc. Cell 151, 56–67 (2012).
Prabhakar, S., Noonan, J.P., Paabo, S. & Rubin, E.M. Accelerated evolution of conserved noncoding sequences in humans. Science 314, 786 (2006).
Bird, C.P. et al. Fast-evolving noncoding sequences in the human genome. Genome Biol. 8, R118 (2007).
Lindblad-Toh, K. et al. A high-resolution map of human evolutionary constraint using 29 mammals. Nature 478, 476–482 (2011).
Fujita, E., Tanabe, Y., Imhof, B.A., Momoi, M.Y. & Momoi, T. Cadm1-expressing synapses on Purkinje cell dendrites are involved in mouse ultrasonic vocalization activity. PLoS One 7, e30151 (2012).
Boyd, J.L. et al. Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex. Curr. Biol. 25, 772–779 (2015).
Buenrostro, J.D. et al. Single-cell chromatin accessibility reveals principles of regulatory variation. Nature 523, 486–490 (2015).
Cusanovich, D.A. et al. Epigenetics. Multiplex single-cell profiling of chromatin accessibility by combinatorial cellular indexing. Science 348, 910–914 (2015).
Langmead, B. & Salzberg, S.L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).
Landt, S.G. et al. ChIP-seq guidelines and practices of the ENCODE and modENCODE consortia. Genome Res. 22, 1813–1831 (2012).
Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9, R137 (2008).
Auerbach, R.K. et al. Mapping accessible chromatin regions using Sono-Seq. Proc. Natl. Acad. Sci. USA 106, 14926–14931 (2009).
Bonn, S. et al. Tissue-specific analysis of chromatin state identifies temporal signatures of enhancer activity during embryonic development. Nat. Genet. 44, 148–156 (2012).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
Quinlan, A.R. & Hall, I.M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842 (2010).
R Development Core Team. R: a Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna, Austria, 2012).
van der Maarten, L.J.P. & Hinton, G.E. Visualizing high-dimensional data using t-SNE. J. Mach. Learn. Res. 9, 2579–2605 (2008).
Saldanha, A.J. Java Treeview--extensible visualization of microarray data. Bioinformatics 20, 3246–3248 (2004).
Akhtar-Zaidi, B. et al. Epigenomic enhancer profiling defines a signature of colon cancer. Science 336, 736–739 (2012).
Bernstein, B.E. et al. The NIH Roadmap Epigenomics Mapping Consortium. Nat. Biotechnol. 28, 1045–1048 (2010).
McLean, C.Y. et al. GREAT improves functional interpretation of cis-regulatory regions. Nat. Biotechnol. 28, 495–501 (2010).
Splinter, E., de Wit, E., van de Werken, H.J., Klous, P. & de Laat, W. Determining long-range chromatin interactions for selected genomic sites using 4C-seq technology: from fixation to computation. Methods 58, 221–230 (2012).
van de Werken, H.J. et al. 4C technology: protocols and data analysis. Methods Enzymol. 513, 89–112 (2012).
Acknowledgements
We thank S. Levine at the MIT BioMicro Center for managing the Solexa pipeline. This work was funded by the Dutch Royal Academy of Sciences (KNAW), the Dutch Cancer Foundation (KWF) grant HUBR2012-5392 (M.P.C.) and the Parkinson's foundation (Stichting ParkinsonFonds) (M.P.C.).
Author information
Authors and Affiliations
Consortia
Contributions
M.W.V. and M.P.C. conceived and designed the experiments. N.B.B. contributed human brain hemispheres. R.B., I.K. and S.P. provided nonhuman primate material from rhesus macaque and chimpanzee. M.W.V., B.C. and P.R. performed the experiments. M.W.V. and S.C.T. analyzed the data, and were supervised by M.P.C. G.G. analyzed 4C experiments, and was supervised by W.d.L. Sequencing at the Utrecht Sequencing Facility was performed by E.d.B., supervised by E.C. M.W.V. and M.P.C. wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Additional information
Netherlands Institute for Neuroscience, Amsterdam, the Netherlands.
Integrated supplementary information
Supplementary Figure 1 Annotation of active cis-regulatory elements in distinct brain regions.
(a) Schematic of the eight brain regions that were analyzed. Prefrontal Cortex (PFC), Precentral Gyrus (PcGm), Occipital Pole (OP), White Matter (WM), Caudate Nucleus (CN), Thalamic Nuclei (TN), Putamen (Put) and Cerebellum (CB). (b) Anatomical location of the eight regions that were sampled from two chimpanzee brains and three rhesus macaque brains. (c) Hierarchical clustering of 24 human brain samples based on normalized H3K27ac enrichment of 59,786 predicted human CREs. Color bar indicates main anatomical subdivisions. (d) PCA of 24 human brain samples on normalized H3K27ac enrichment of 59,786 human CREs.
Supplementary Figure 2 Comparison of H3K27ac enrichment in the human cortex with previously defined enhancer and promoter states.
In the upper panel, the pie chart illustrates the fraction of overlap between H3K27ac-enriched regions in the cortex and enhancer and/or promoter states defined in the dorsolateral prefrontal cortex. Bar chart further specifies exact categories that comprise more than 10%. The lower panel depicts the fraction of different chromatin states in the dorsolateral prefrontal cortex24 that is enriched for H3K27ac.
Supplementary Figure 3 Characterization of putative CREs in the brain.
(a) Selection strategy for human and rhesus macaque CREs with orthologues of comparable mappability on all three primate genomes. (b) Size distribution of all CREs in human and rhesus macaque. (c) Ratio of rhesus and human CRE length. The fraction that does not differ more than 25% in size between both species is indicated. (d) Enrichment analysis for genes within indicated modules for linked brain enhancers. (e) Percentage of genes within modules of a given tissue33 as indicated below the panels, that are linked by proximity to cortical or cerebellar brain-specific enhancers, as compared to 13 unrelated tissues. Fisher’s exact test was used to compare enrichment for brain enhancers between modules; *** P < 0.0005. (f) Venn diagram showing the overlap of putative enhancers between the main anatomical subdivisions of the brain. Brain-specific enhancers per anatomical subdivision were analyzed for linked gene ontologies using GREAT.
Supplementary Figure 4 Positional conservation of enhancers and promoters in the brain across primate evolution.
(a) Upper panels show the fraction of enhancers and promoters in cerebellum (CB) and prefrontal cortex (PFC) that is enriched, as defined by peak calling, in both human and rhesus macaque (blue box) or one species only (red box). Species-specific CREs are further subdivided using DESeq2 into significantly differentially enriched (DE higher in human, purple; DE lower in human, light brown) and not DE (grey). Lower panels show normalized RNA read counts, indicated for human and rhesus macaque by silhouettes, for genes close to DE and not DE enhancers. Dissimilarity between distributions was calculated using a Wilcoxon rank-sum test; *** P < 0.0005, ** P < 0.005, * P < 0.05, # P > 0.05. (b) Fraction of enhancers (left) and promoters (right) for each brain region that is positionally conserved between human and rhesus macaque (blue and grey). (c) Average PhastCons scores for all nucleotides within the indicated elements. (d) Fraction of reads within CREs that was mappable onto the other genome for (not) positionally conserved CREs. (e) Upper panels show the fraction of enhancers and promoters in CB and PFC that is enriched, as defined by peak calling, in both human and rhesus macaque (blue box) or one species only (red box). Shared CREs (blue) were further subdivided into DE, higher in human (purple), DE, lower in human (light brown) and not DE (grey) using DESeq2. Lower panels show normalized RNA read counts, indicated for human and rhesus macaque by the silhouettes, for genes close to (not) differentially enriched enhancers. Dissimilarity between distributions was calculated using a Wilcoxon rank-sum test; *** P < 0.0005, * P < 0.05, # P > 0.05.
Supplementary Figure 5 Differential H3K27ac enrichment corresponds to changes in gene expression.
(a) Heatmaps and boxplots for the eight brain regions illustrate (the ratio of) average read counts (normalized for number of reads in peaks) of replicates in human and rhesus macaque. Upper panels in purple depict DE enhancers, higher in human. Middle panels in light brown are DE enhancers, lower in human. Lower panels in grey represent enhancers that are not DE. (b) Percentage of human promoters per brain region that is (not) DE between human and rhesus macaque as defined by DESeq2 (≥ twofold change, FDR < 0.01). The number of analyzed promoters per subdivision is indicated in the color bar below the graph. (c) Boxplot showing the comparison of normalized RNA read counts from human and rhesus macaque brain samples as indicated by the silhouettes. Displayed are the closest genes to a (not) DE enhancer in the prefrontal cortex (PFC). Dissimilarity between distributions was calculated using a Wilcoxon rank-sum test; # P > 0.05. (d) Comparison of normalized RNA read counts for genes with (no) differential H3K27ac enrichment at their promoters (shown for CB and PFC). Dissimilarity between distributions was calculated using a Wilcoxon rank-sum test; *** P < 0.0005, ** P < 0.005, * P < 0.05, # P > 0.05. (e) Ratio of enhancer length between rhesus macaque and human for (not) DE enhancers. (f) Left panels show a RPM normalized ChIP-seq read distribution (axis limit 5 RPM) for H3K27ac at a genomic region spanning 115kb containing the DFNB31 gene. Eight tracks represent the different anatomical brain regions that were analyzed in both human (HS1) and rhesus macaque (RM1) as indicated by the silhouettes. Blue shade highlights the repurposed enhancer and promoter. The right panel displays RPKM normalized RNA sequencing reads for DFNB31 relative to the human cortex.
Supplementary Figure 6 Changes in active chromatin at predicted CREs occurred mainly before the human-chimpanzee separation.
(a) t-Distributed Stochastic Neighbor Embedding (t-SNE) analysis for all human, chimpanzee and rhesus macaque brain samples using normalized H3K27ac enrichment at 60,702 predicted human and rhesus CREs with orthologues on all three genomes. (b) Upper panels show scaled H3K27ac enrichment of (not) DE enhancers as indicated between human and rhesus macaque in prefrontal cortex (PFC). In purple, DE enhancers, higher in human. In light brown, DE enhancers, lower in human. In grey, not DE enhancers. Lower panels display normalized RNA read counts from human, rhesus macaque and chimpanzee brain samples for genes close to (not) DE enhancers. Dissimilarity between distributions was calculated using a Wilcoxon rank-sum test; * P < 0.05, # P > 0.05. (c) Fraction of (not) DE enhancers in cortical and subcortical samples defined between human and rhesus macaque with a twofold difference in normalized read counts (mean of replicates) between human and the indicated species.
Supplementary Figure 7 Newly introduced and specifically depleted CREs in humans.
(a) RPM normalized ChIP-seq read distribution (axis limit 5 RPM) for H3K27ac at a genomic region spanning ~173kb. Eight tracks represent the different anatomical brain regions that were analyzed in HS1, Ch1 and RM1 as indicated by the silhouettes. Colored boxes highlight CREs that are significantly higher in human (purple, 3 putative enhancers and 1 promoter), significantly lower in human (light brown, 1 enhancer) and gained in great apes (yellow, 1 enhancer in the subcortical structures). (b) Same analysis as (a) for a human-specific depletions. (c) Left graph shows distributions of human enhancers over the subdivisions of the brain. Fractions of shared enhancers between subdivisions are plotted in diagonal stripes that are colored according to the different subdivisions. This is shown for all predicted human enhancers and for newly introduced enhancers in comparison to rhesus macaque (random sampling, P < 0.001). Right graph shows fraction of new enhancers per (shared) anatomical subdivision of the original set of all enhancers. (d) Similar analysis as in (c) for human-specific depletions.
Supplementary Figure 8 Chromosome conformation capture for a CRE at accelerated DNA.
In green, enrichment for H3K27ac in the human cerebellum for a genomic region containing CADM1. In blue, circular chromosome conformation capture (4C-seq) using the genomic region containing HAR87 as a viewpoint (exact location indicated by green arrow) in human cerebellum. The second CADM1 enhancer is illustrated with a green arrow. Black boxes represent genes. Interaction peaks with both the second CADM1 enhancer and the CADM1 promoter are indicated using a red arrow.
Supplementary information
Supplementary Text and Figures
Supplementary Figures 1–8 (PDF 1743 kb)
Supplementary Table 1
Background information on the brain donors used in this study. (XLSX 29 kb)
Supplementary Table 2
GEO Sample IDs, anatomical location, alignment summary, number of called peaks, FRiP score and read length for each sample. (XLSX 26 kb)
Supplementary Table 3
List of the H3K27ac-enriched regions in human. (XLSX 3110 kb)
Supplementary Table 4
List of the H3K27ac-enriched regions in chimpanzee. (XLSX 3053 kb)
Supplementary Table 5
List of the H3K27ac-enriched regions in rhesus macaque. (XLSX 3202 kb)
Supplementary Table 6
List of the 60702 human and rhesus macaque CREs with mappability on the three genomes. (XLSX 6261 kb)
Supplementary Table 7
Differentially enriched CREs between human and rhesus macaque. Enhancers higher in human. (XLSX 1253 kb)
Supplementary Table 8
Differentially enriched CREs between human and rhesus macaque. Enhancers lower in human. (XLSX 1071 kb)
Supplementary Table 9
Differentially enriched CREs between human and rhesus macaque. Promoters higher in human. (XLSX 184 kb)
Supplementary Table 10
Differentially enriched CREs between human and rhesus macaque. Promoters lower in human. (XLSX 96 kb)
Supplementary Table 11
Gene Ontology output for genes close to a DE higher and DE lower in human enhancer. (XLSX 13 kb)
Supplementary Table 12
Differentially enriched CREs between human and chimpanzee. Enhancers higher in human. (XLSX 277 kb)
Supplementary Table 13
Differentially enriched CREs between human and chimpanzee. Enhancers lower in human. (XLSX 275 kb)
Supplementary Table 14
Differentially enriched CREs between human and chimpanzee. Promoters higher in human. (XLSX 33 kb)
Supplementary Table 15
Differentially enriched CREs between human and chimpanzee. Promoters lower in human. (XLSX 28 kb)
Supplementary Table 16
Gene Ontology output for genes close to DE enhancers of human versus rhesus and chimpanzee. (XLSX 14 kb)
Supplementary Table 17
List of newly enriched enhancers in human. (XLSX 109 kb)
Supplementary Table 18
List of depleted enhancers in human. (XLSX 97 kb)
Supplementary Table 19
List of newly enriched promoters in human. (XLSX 13 kb)
Supplementary Table 20
List of depleted promoters in human. (XLSX 12 kb)
Supplementary Table 21
List of CREs that overlap a human accelerated region, with specification of whether the region is differentially enriched. (XLSX 55 kb)
Supplementary Data 1
This file provides raw read counts for all CREs per species. (ZIP 8662 kb)
Supplementary Data 2
This file provides raw and normalized read counts for the 60702 CREs with comparable sequence content in the three primate genomes. (ZIP 34430 kb)
Rights and permissions
About this article
Cite this article
Vermunt, M., Tan, S., Castelijns, B. et al. Epigenomic annotation of gene regulatory alterations during evolution of the primate brain. Nat Neurosci 19, 494–503 (2016). https://doi.org/10.1038/nn.4229
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nn.4229
This article is cited by
-
Epigenetic regulation of human-specific gene expression in the prefrontal cortex
BMC Biology (2023)
-
Developmental mechanisms underlying the evolution of human cortical circuits
Nature Reviews Neuroscience (2023)
-
Current advances in primate genomics: novel approaches for understanding evolution and disease
Nature Reviews Genetics (2023)
-
Comparative epigenomics reveals the impact of ruminant-specific regulatory elements on complex traits
BMC Biology (2022)
-
Inferring mammalian tissue-specific regulatory conservation by predicting tissue-specific differences in open chromatin
BMC Genomics (2022)