Skip to main content

Human evolution: the non-coding revolution


What made us human? Gene expression changes clearly played a significant part in human evolution, but pinpointing the causal regulatory mutations is hard. Comparative genomics enabled the identification of human accelerated regions (HARs) and other human-specific genome sequences. The major challenge in the past decade has been to link diverged sequences to uniquely human biology. This review discusses approaches to this problem, progress made at the molecular level, and prospects for moving towards genetic causes for uniquely human biology.

Post-genomic challenges for determining uniquely human biology

When the human genome was first sequenced [1, 2], the big question was “how many genes do we have?” Most people guessed too high. Sequencing our closest living relative, the chimpanzee [3, 4], begged the question “which genes are different?” Here the answer was predicted a century before and supported by King and Wilson’s 1975 discovery that certain blood proteins have very few amino acid differences between human and chimpanzee [5, 6]. We now know that the vast majority of all genomic changes that happened since the human–chimpanzee ancestor are in non-coding regions, consistent with King and Wilson’s hypothesis that regulatory changes drove the differences between our species. In hindsight, the importance of gene regulation in human evolution is logical. There are many more DNA bases in regulatory regions than in protein-coding genes, making them a larger target for evolutionary innovation. Furthermore, genes frequently function in many different contexts, and this pleiotropy constrains their evolution compared to regulatory elements, which tend to be more modular [7]. Thus, regulatory sequences have great potential to be drivers of human evolution.

The challenge in the post-genomic era has been to determine which of the millions of human-specific non-coding sequence differences are responsible for the unique aspects of our biology. This is a hard problem for many reasons. First, the non-coding genome is vast, requiring methods to prioritize the mutations that matter. The neutral theory of molecular evolution, coupled with redundancy in biological networks, suggests that many human-specific DNA changes had little effect on our biology. Second, we know much less about how sequence determines function of regulatory elements compared to protein or RNA genes. Hence it is difficult to predict the molecular, cellular, and organismal consequences of human-specific regulatory mutations. Furthermore, most uniquely human traits are complex, and there is no doubt that they are encoded by a combination of mutations in different genomic loci. Finally, because gene regulation has diverged significantly between primates and model organisms such as mice, zebrafish or flies, it is hard to test hypotheses about the functional effects of regulatory mutations. In this review, we discuss advances to address these barriers with an emphasis on linking sequence to function, complementing other recent papers that explore genetic and regulatory changes in human evolution [8,9,10,11,12,13].

Discovering the fastest evolving regions in the human genome

Single nucleotide changes can have functional consequences, but currently these are difficult to predict in non-coding regions where small mutations are frequently tolerated and the function of a particular nucleotide is rarely known. Hence, human evolutionary genetics has mostly focused on genome regions with many human-specific differences (reviewed in [14,15,16,17]). Human accelerated regions (HARs) are short, evolutionarily conserved DNA sequences that have acquired significantly more DNA substitutions than expected in the human lineage since divergence from chimpanzees. A number of studies applied different tests to identify HARs either genome-wide [18, 19] or with protein-coding sequences specifically removed from the analysis [20,21,22,23]. We and the other authors of these studies had a common aim: to identify regulatory elements with human-specific activity (Fig. 1). These analyses started with regions conserved across non-human mammals in order to enrich for functional elements [24,25,26,27] and to increase power to detect acceleration (Box 1). Then they used various methods to identify a subset of conserved elements that accumulated human-specific changes. Differences in analysis choices and available data over time (for example, species in alignments, methods used to identify conserved elements, tests for acceleration, bioinformatics filters to remove artifacts) resulted in only modest overlap between the HARs identified in different studies despite their common aim (Fig. 1). These studies also differed in whether they specifically tested for positive selection compared to a neutral model or simply identified acceleration, which could be due to a variety of evolutionary processes including mutation rate increases or loss of constraint (Box 1). Collectively nearly 3000 non-coding HARs have been identified to date [14], representing a pool of candidates that can be searched for regulatory regions with human-specific activity.

Fig. 1.
figure 1

Identification of human accelerated elements. Top: the four different approaches used to identify human accelerated regions. Some key differences include (i) the conserved elements used as candidates to identify HARs (which depend on multiple sequence alignments, methods to detect conservation, and whether human was masked in the alignments), (ii) bioinformatics filters that aim to restrict to non-coding elements and/or remove assembly or alignment artifacts, and (iii) tests used to detect acceleration. Bottom: overlap of the different datasets of human accelerated regions. Abbreviations: ANC accelerated conserved non-coding sequences [20]; HACNS human accelerated conserved non-coding sequences [23]; HTBE human terminal branch elements [21]. HARs include the original HARs [19] and the second generation HARs or 2xHARs [100]

Insertions, deletions, duplications, and rearrangements—collectively known as structural variants (SVs)—contribute many more nucleotides to the genetic difference between humans and chimpanzees than do HARs and other single nucleotide variants. Indeed, the first differences detected between the human and chimp genomes were large SVs discovered prior to genome sequencing (for example, chromosome 2 fusion, inversions on chromosomes 1 and 18 [28]). Primates have accumulated SVs at an accelerated rate compared to other mammals, and the human genome contains numerous lineage-specific SVs [29, 30]. These SVs contribute significantly to the emergence of new genes and gene families (reviewed in [31]). Most HAR analyses have specifically filtered out SVs due to the difficulty of accurately assembling and aligning these regions. But HARs are in fact common in recent segmental duplications [20], and we found them to be enriched near duplicated genes [32]. In addition, SVs can change gene expression and phenotypes by associating non-coding elements with genes they did not previously regulate (‘enhancer hijacking’, see below) [33, 34]. Thus, further exploration of the regulatory consequences of human-specific SVs is clearly needed.

Beyond conservation: combining acceleration tests and epigenetic marks

Sequence conservation is a useful tool for predicting which non-coding substitutions will be functional. But many regulatory elements are not conserved [35, 36], and conservation does not indicate when and where a regulatory element is active. Since the initial discovery of HARs and human-specific SVs, understanding of the proteins and epigenetic marks found at particular categories of regulatory element has improved significantly, as have genome-wide functional genomics assays such as chromatin immunoprecipitation sequencing (ChIP-seq), ATAC-seq, and RNA-seq. Using these tools it is now possible to generate genome-wide maps of predicted regulatory regions across many cell types and developmental stages. We predicted enhancers in the human genome and tested 29 predictions overlapping HARs with in vivo reporter assays (see below), discovering that many of the previously identified HARs are enhancers active in the developing embryo [37].

Functional genomics data are also very useful for loosening or omitting the requirement that HARs be conserved in other species. Two recent investigations performed evolutionary analyses of putative human regulatory elements and identified many examples of human sequence acceleration missed by studies that relied on deep sequence conservation (for example, greater than 95% sequence identity across mammals or vertebrate phastCons elements). Gittelman et al. performed comprehensive evolutionary and population genetics analyses on open chromatin (specifically, DNase I hypersensitive sites (DHSs)) from 130 cell types to discover regulatory DNA with evidence of adaptive evolution in humans [38]. They discovered 524 DHSs with sequences that are conserved in non-human primates but show accelerated nucleotide substitution rates in the human lineage (haDHS). By functionally characterizing selected haDHSs using transgenic reporters and luciferase assays in different cell types, the authors identified several where the human sequence changes result in a gain of enhancer function (Table 1; Fig. 2).

Table 1 Selected examples of human-specific regulatory elements tested in functional assays that showed gain or loss of function
Fig. 2.
figure 2

Strategies to identify human-specific enhancers. a Methodology used to identify human-accelerated DNSase I hypersensitive sites (haDHSs) in [38]. Black bars, nucleotides that differ from the human sequence; blue bars, sites where all species differ from humans; dotted red lines, DHS and neutral sequence. b Five haDHSs showed different activity with human versus chimpanzee sequences in SK-N-MC cells. c Differentially active haDHS12/DAR12 (asterisk) overlaps previously identified HACNS219 [23] and is located near RNF145, a ring finger gene involved in cellular cholesterol metabolism [101]. d Brains of species studied in [43] (approximately to scale, colors label regions). Tree indicates approximate timing of splits between lineages (millions of years ago). e Workflow and (bottom) fraction of H3K27ac peaks that is differentially enriched (DE) between species per brain region. f Left: Percentage of HARs within conserved regulatory elements (CREs). Right: Most of the 240 HAR-containing CREs that align across species were not DE in human versus macaque and chimpanzee. g H3K27ac tracks for mouse, macaque, chimpanzee, and human cerebellum. Gray, shared enhancers; purple, DE enhancers higher in human versus macaque; red, HAR87. Enhancers gained in mouse or lost in primates (yellow) may compensate for the enhancer gained in primates (light blue) upstream of the CADM1 gene. Abbreviations: CB cerebellum, PFC prefrontal cortex, PcGm precentral gyrus, OP occipital pole, WM white matter, CN caudate nucleus, TN thalamic nucleus, Put putamen. Reproduced partially from [38, 43], with permission

Using a similar approach, Dong et al. compared substitution rates of DHSs to nearby ancestral repeats (LINE1 or LINE2 elements, assumed to be neutrally evolving) and identified 3538 accelerated DHSs (ace-DHS) [39]. Notably, only 17 accelerated DHSs overlap between these two studies, likely due to Gittelman et al. requiring conservation in primates so that they tested only 113,577 DHSs versus 808,943 in Dong et al.’s analysis, which had no such conservation filter. In addition, these two approaches used different data to estimate neutral rates for their acceleration tests. Regardless, it is clear that starting with DHSs rather than conserved elements reveals novel human-specific regulatory elements: the majority of haDHSs (454/524) and ace-DHSs (3520/3538) were not previously identified as HARs. It is important, however, to validate the regulatory functions of DHSs and other biochemically active genomic regions experimentally, especially if they are not conserved [36].

Comparative identification of regulatory regions

An alternative approach to identify human-specific regulatory elements is to perform functional genomics assays such as open chromatin or protein binding in multiple species and find regions with differential evidence of regulatory activity. Some of the first evidence that human-specific regulatory elements may be common came from the ENCODE project, which found that a high proportion (~50%) of the biochemically active regions identified in the human genome were not conserved across mammals [35, 40]. Several recent studies built on this observation and compared human epigenomic profiles (such as open chromatin, histone modifications) to those of primates and other species using either primary tissues or cell lines from various developmental stages (see below). To infer that regulatory marks are human-specific requires epigenomic profiles from chimpanzees (or bonobos), which is often impossible due to limited availability of material and ethical considerations. One potential solution is to investigate post hoc if the regions with epigenetic signatures that differ between human and rodents or monkeys intersect with HARs or harbor sequence changes that are unique to the human genome.

Cotney and colleagues performed ChIP-seq for histone H3 lysine 27 acetylation (H3K27ac, a mark of active enhancers) in human, rhesus macaque, and mouse embryonic limb at four different developmental stages from bud to digit separation, defined as orthologous based on morphology and HOXD gene expression [41]. They found 2175 promoters and 2915 putative enhancers with significantly higher H3K27ac signal in human at one or more time points. Sixteen of these regions overlap with HARs, but most show no evidence of human accelerated substitution rates or unusual diversity patterns across modern humans [41].

Using a similar approach, Reilly et al. mapped active promoters and enhancers using di-methylation of histone H2 lysine 27 (H3K27me2, a promoter mark) and H3K27ac, respectively, during three stages of human, rhesus macaque, and mouse cortex development [42]. Comparing signal for these epigenetic marks across species, the authors predicted 2855 promoters and 8996 enhancers unique to human samples (that is, with higher signal in human than other species). The majority of these regions show no evidence of human accelerated substitution rates, although 48 of them overlap HARs [42].

Vermunt and colleagues assayed H3K27ac in eight anatomical subdivisions of the adult brain (cerebellum, caudate nucleus, thalamic nuclei, putamen, white matter, precentral gyrus, prefrontal cortex, occipital pole) of human, chimpanzee and rhesus macaque [43]. The authors found that 14–43% of putative enhancers and 3–10% of promoters had differential signal between human and macaque brain regions. However, only a very small fraction of these showed a similar difference between human and chimpanzee. While 284 H3K27ac regions from this study overlap with HARs, only 32 of these have significantly higher or lower H3K27ac signal in human versus chimpanzee or macaque. One particularly interesting example is 2xHAR87/HACNS548, which interacts with the promoter of the CADM1 gene and shows steadily increasing H3K27ac signal from macaque to human (Fig. 2).

Prescott et al. compared transcription factor (TF) and co-activator binding, histone modifications, and chromatin accessibility genome-wide in human and chimpanzee cranial neural crest cells (CNCCs) to pinpoint putative enhancers [44]. They derived CNCCs from pluripotent stem cells using an in vitro protocol in which specification, migration, and maturation are recapitulated in the dish [45, 46], an approach that allows unavailable chimpanzee cell types to be derived in the lab. The authors predicted ~ 1800 regulatory elements with different levels of H3K27ac in human versus chimpanzee, suggestive of differential enhancer activity. A few of these overlap HARs (three of the top 1000; Table 1). The authors found that the variance in H3K27ac signal scales proportionally with human–chimp sequence divergence. We note that sequence differences are quite low overall (three to six substitutions per 500 bp) even at regions with the biggest human–chimpanzee differences in H3K27ac binding [44]. These results show that regions with one or two sequence changes can have different enhancer activity (measured by transgenic assays using luciferase or lacZ), which may make identifying causal sequence changes easier than in HARs with five to ten changes. Prescott et al. speculate that the enhancer activity differences they detect are due to changes in consensus transcription factor binding motifs, such as one called Coordinator. However, since enhancer assays are not highly quantitative (see below) and cannot identify upstream causes of differences in reporter gene activity, further functional experiments are important to validate the expression differences and test the hypothesis that transcription factor binding changes played a causal role in these regions.

Overall, these studies discovered many genome sequences with human-specific epigenomic signatures, which are exciting catalogs of putative regulatory elements for future studies. Since many of these regions are not deeply conserved across mammals, they complement previously identified HARs. However, most of the elements show little or no evidence of human-specific sequence change, suggesting that the causal mutations are in trans or that the catalogs contain many false positives due to the many challenges of cross-species functional genomics, such as antibody affinity, low samples sizes, and matching cell composition and developmental stage. As Prescott et al. note, regions with fewer sequence changes than a typical HAR can show enhancer activity differences [44], though we emphasize that it is challenging to discriminate these from the many cases where an equivalent number of sequence changes causes no activity difference. More functional validation of predicted differences in regulatory activity (for example, with reporter assays or CRISPR/Cas9 screens) will help to determine the false positive rate. Meanwhile, we should be careful when interpreting chromatin immunoprecipitation data since protein–DNA binding events cannot be directly translated into regulatory function, as has been discussed elsewhere [47,48,49]. Another challenge is evaluating the net effect of all genetic and epigenetic changes within a locus on species-specific gene expression (reviewed in [49]), which is particularly difficult when regulatory elements function synergistically or redundantly (Fig. 3 and below).

Fig. 3.
figure 3

Evolutionary mechanisms at the level of gene regulatory regions. a An example of accelerated sequence evolution affecting one enhancer in a locus and leading to gain of an expression domain in the developing forebrain. We characterized one such gain of function HAR in Kamm et al. [57]. b Shadow enhancers are multiple enhancers that direct a similar gene expression pattern and thus overlap in function. As such, they can act cooperatively to confer robustness in different physiological situations [77, 78]. The example shown is based on the work of Lam et al. [78], where two enhancers direct expression to the Arcuate Nucleus and both must be deleted in mice to produce a dramatic change in expression and phenotype. c Enhancer turnover is when one enhancer disappears and a new one appears in the same regulatory region, replacing the lost function. In the example shown by Domené et al. [72], turnover resulted in no net change in expression or phenotype

Transgenic approaches to study the functional impact of non-coding evolution

One approach to investigating the functional effects of human-specific non-coding sequences is to introduce them into model organisms. Despite differences in anatomy and genetic networks across species [50], transgenic and mutant mouse models have been useful for characterizing human-specific genes, such as ARHGAP11B [51] and SRGAP2C [52], and human-specific gene variants, such as FOXP2 [53, 54]. Predicted enhancers can also be validated and tested for effects of sequence differences using transient or stable reporter assays in mouse and zebrafish. In general, in transient transgenic enhancer assays (where the F0 generation is analyzed) enhancer function can be evaluated in only one or a small number of embryonic stages, and not in adult tissues. Although we found that many HARs appear to function during development [10], a substantial fraction likely cannot be characterized with this approach. Generation of stable transgenic lines (where the F1 generation is analyzed) could facilitate the analysis of many more developmental stages and adult expression. However, the low throughput of transgenic experiments has limited their utility for characterizing large catalogs of human-specific regulatory elements.

Despite these caveats, a number of human-specific non-coding sequences have been functionally characterized as enhancers in vivo with reporter assays in transgenic animal models. These include two conserved elements deleted in humans (hCONDELs) [55], nearly 70 HARs (summarized in [10]), including nine haDHSs, and an additional seven haDHSs that do not overlap HARs [38]. Reporter assays in transgenic animals can also be used to test the hypothesis that human-specific mutations in non-coding elements altered their enhancer activity, and indeed several HARs and haDHSs show expression differences between constructs carrying the chimpanzee versus human sequence [37,57,58,, 5659] (Table 1). It is important to note that transient transgenics can capture gains and losses of enhancer activity in specific tissues or cell types (expression patterns), but they are not quantitative (expression level depends on number and location of random genomic integrations) and hence are unable to capture changes in activity levels. These studies aiming to analyze the function of human-specific non-coding sequences and then to study them comparatively with transgenic model organism enhancer assays are adding valuable information about the functional impact of human-specific DNA changes. Looking ahead, it will be important to supplement these with additional approaches to link enhancer activity differences to gene expression and phenotypes.

Massively parallel reporter assays (MPRAs) show promise in cell lines and may become practical in whole animals [60,61,62]. In this technique, a library of candidate regulatory enhancer DNA sequences is cloned into a reporter construct containing a unique DNA barcode that will be transcribed if the enhancer is active. RNA-sequencing of the DNA barcodes enables quantitative measurements of enhancer activity, including effects of individual nucleotide variants [62]. Two potential limitations of this technique are length of the candidate enhancers (<200 bp with current DNA synthesis methods) and the fact that they are usually evaluated outside of their genome context. Knock-in strategies where regulatory element variants are replaced in the analyzed locus (such as via CRISPR/Cas9 genome editing) could address this issue. Low throughput transgenic and mutant model animals can complement both techniques. A combination of methods can together enable tests of hypotheses about the effects of human-specific variants on molecular and organismal phenotypes.

Human variation links accelerated regions to phenotypes

Another way to associate human-specific non-coding regions with traits is to investigate associations between any polymorphisms they harbor and phenotypic variability in humans. The polymorphic sites will be relatively recent mutations that are largely distinct from differences between human and chimpanzee reference genomes. But they may nonetheless shed light on the general function of the accelerated regions in which they occur. For example, gene copy number variants in the human-specific pericentric inversion of chromosome 1 have been associated with human developmental and neurogenetic diseases [15,64,, 52, 6365], and we showed that HARs that are conserved in primates but not across all mammals are associated with schizophrenia [66]. As many more human genomes are sequenced, population genetic variation in and near to HARs will also help researchers to test hypotheses about the evolutionary forces that created and maintain HARs (Box 1, Fig. 3).

Leveraging this approach, Doan et al. sequenced HARs in individuals with autism spectrum disorder (ASD) and unaffected controls [60]. ASD was associated with a significant 43% excess of rare biallelic variants in HARs. The authors showed that 29% of these variants alter enhancer activity in primary mouse neurospheres using a custom MPRA. They further identified rare homozygous mutations with active regulatory marks near neurodevelopmental and disease-associated genes in ASD patients who lack causal coding mutations. Using luciferase reporter assays and transgenic mice, the authors show that one particular ASD-associated HAR variant (HACNS426) that was previously shown to interact with the neuronal morphology gene Cux1 [67,68,69,70] is an enhancer (Fig. 4). The authors report that transgenic mice carrying the mutant allele fused to the CUX1 promoter and to GFP show elevated expression of the reporter gene in the developing brain compared to the wild-type allele. But it is important to remember that transgenic enhancer assays cannot capture quantitative expression differences (see above).

Fig. 4
figure 4

Testing HARs and HAR sequence variants for enhancer activity with reporter assays. a Example of a transient transgenic reporter assay to test a HAR (HACNS426) for enhancer activity in mouse embryos. The experiment compares enhancer activity of the major allele (G) to that of the autism-associated minor allele (A) that is never homozygous in healthy controls and is predicted to change transcription factor binding. Top: constructs carrying each of the two HACNS426 alleles fused to the human CUX1 promoter and the GFP reporter gene are separately injected into single-cell mouse embryos. Bottom: to assay enhancer activity, brain slices from embryonic day E16.5 are stained for GFP. The authors observed differences in GFP expression with the G allele (above) versus A allele (below). The major strength of the approach is the spatial and cellular resolution of in vivo measurements, while weaknesses include not being highly quantitative, the use of mouse to compare human and chimpanzee variants, low throughput, and relatively high cost. Adapted with permission from [23]. The study also performed in vitro luciferase reporter assays and showed HANCS426 interacts with the dosage-sensitive CUX1 promoter. b Massively parallel reporter assays (MPRAs) enable thousands of reporter constructs to be tested as a library (top) in which each HAR variant is associated with a unique DNA barcode (such as 20-bp sequence). RNA sequencing of barcodes (bottom) provides a quantitative estimate of the activity of each HAR variant. MPRAs are high throughput, allowing thousands of HARs and variants thereof to be tested, and they are quantitative, enabling detection of single nucleotide differences with moderate effects on expression. Current weaknesses of the technology include being limited to HARs or HAR segments less than 200 bp and being restricted to testing in cell lines or mouse tail vein assays

The study of Doan et al. demonstrates the potential for human sequencing studies coupled with deep phenotyping to shed light on HAR function. Using associations between specific HARs and ASD, the authors conclude that many HARs are essential for normal development. Even though the ASD-associated mutations are different from the substitutions that happened during human evolution, they nonetheless help to annotate the function of the genomic region containing a HAR. Polymorphic deletions of HARs could similarly be used to annotate their functions. Doan et al. further hypothesize that human-specific mutations in HARs could have altered social and/or cognitive behavior. To test this hypothesis precisely with medical re-sequencing requires that specific human–chimp differences and combinations thereof be polymorphic in living humans or primates. As targeted sequencing (for example, with molecular inversion probes [71]) enables HARs to be sequenced in millions of individuals, it will become clear how many human–chimp differences are present in modern humans. One caveat is that the genomes of modern humans and other primates differ from those of ancestral populations, so phenotype associations of HAR mutations must be interpreted in the context of the modern genetic background. For this approach to be a useful path to elucidating HAR functions in humans, it will also be important to couple sequencing data with a very wide range of reliably measured phenotypes ranging from diseases to behaviors. This is an exciting and increasingly feasible direction.

Changing things so everything stays the same? Regulatory element turnover and compensatory evolution

To measure the effects of non-coding mutations on expression and phenotypes, we must move towards considering the whole regulatory region of a gene, as well as changes in trans (such as, expression and activity of upstream regulators). One reason for this is that regulatory elements without sufficient sequence similarity to be aligned across species can direct conserved expression patterns, a phenomenon known as ‘enhancer turnover’ [72,73,74] (Fig. 3). Another consideration is the fact that regulatory elements interact with each other to drive gene expression, which means that non-coding mutations can amplify or buffer one another even if they are not adjacent. One particular type of compensation across regulatory elements involves ‘shadow enhancers’, which are distinct sequences capable of guiding similar spatiotemporal expression patterns [75,76,77,78,79,80] (Fig. 3). This apparently redundant function may ensure robust, precise, and stable regulation [77, 78]. Thus, gene expression evolution cannot be easily predicted from a single human-specific regulatory element without evaluating the complete regulatory landscape of the locus.

Most human-specific regulatory elements have nonetheless been characterized one at a time, or even one mutation at a time. Supporting the idea that this limits our ability to predict expression divergence, Vermunt et al. showed that many loci contain both gains and losses of predicted enhancers (4.9–10.3% across brain regions of mice, macaques, chimpanzees, and humans) [43]. For example, the locus containing HAR87 harbors multiple other predicted enhancers with increased or decreased activity during human evolution, measured via chromatin capture coverage (Fig. 2). These other changes could compensate for or otherwise interact with changes in the activity of HAR87. To investigate this possibility, comparative genomics should be combined with functional genomics assays, especially those that do not depend on antibody affinities (for example, chromatin capture, open chromatin).

An emerging approach is to leverage high-throughput genetic screens that can probe combinations of mutations across a locus, for example, via CRISPR genome editing [81, 82]. Using a combination of CRISPR activation and inhibition, as well as knock-outs and knock-ins, it will be possible to study the effect of human–chimpanzee non-coding sequence differences on gene expression in cells or in models such as mice. These studies should investigate all the putative regulatory regions of a gene individually and in combination to truly decipher the net regulatory effects.

Thinking again: what are we missing in the study of genetic bases of human evolution?

The past decade has seen significant progress towards addressing the major hurdles of associating human-specific non-coding elements with traits that make our species unique. Functional genomics has been transformative in terms of annotating and prioritizing HARs, and it has also been used directly to identify novel human-specific regulatory elements that are not divergent enough in sequence to be detectable in genome-wide tests for accelerated sequence evolution. As the mechanisms of gene regulation are increasingly understood, the functional effects of human-specific non-coding mutations are becoming less mysterious. High-throughput functional assays, such as massively parallel reporter assays and CRISPR screens, are increasing the numbers of characterized human-specific regulatory elements by orders of magnitude. Meanwhile, targeted and whole genome sequencing of thousands of people is enabling rare variant association and linkage studies that connect non-coding elements to traits and may even reveal the phenotypes of individuals carrying specific ancestral versus derived haplotypes. Human evolutionary studies will always be challenged by our inability to test specific hypotheses about genetic changes in the correct organismal context (an ancestral human). However, we now have a much broader collection of complementary research tools to address this problem, including on one hand genome editing of stem cell-derived cell lines, organoids, human primary cells, and model organisms, and on the other hand deep phenotyping and sequencing of humans and other primates to capture natural variation.

Despite these exciting advances towards understanding the role of regulatory changes in human evolution, it is important to ask what hurdles still remain. One major barrier is our very limited ability to assay functions of regulatory elements other than enhancers. Genome editing is helping to address this bias by enabling researchers to introduce individual human mutations or to knock out single human-specific elements [81, 82]. By performing these experiments in primate cells, human cells, and model organisms, it will be feasible to characterize putative human-specific insulating and repressing elements in terms of their downstream effects on molecular and cellular phenotypes. The development of new high-throughput assays to read out regulatory functions other than enhancing gene expression (such as repression or insulation) will help as well.

Another missing piece in the story of human regulatory evolution is the vast sequence differences encoded within SVs, many of which are missing or incorrect in human and non-human primate reference genomes [83]. Even within the SVs that have been accurately compared across primates, human non-coding evolution is relatively unexplored, due largely to the technical challenges of functionally testing regions with duplications and complex genomic architectures. SVs have immense potential to alter regulatory elements and their interactions with target genes. For example, inserting or deleting topologically associating domain (TAD) boundaries, which are often conserved across cell types and species [84], can associate regulatory elements with new genes (termed ‘enhancer hijacking’) or insulate them from their ancestral gene targets [85, 86]. It will be exciting to explore the evolution of human gene expression in terms of 3D genome organization and its effects on regulatory interactions.

Outlook: what will it take to crack the code?

The next decade could be the one in which the regulatory code is cracked, opening the door to reading out the functional effects of non-coding changes that distinguish humans from other primates. To realize this goal, we should not only push to characterize the human–chimp differences for which we already have functional hypotheses. We must also continue to ask ‘what genomic regions are missing from our analyses?’ and ‘where do current models of gene regulation fail to explain divergently expressed genes?’


  1. 1.

    Consortium HGS. Finishing the euchromatic sequence of the human genome. Nature. 2004;431:931–45.

    Article  Google Scholar 

  2. 2.

    Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, et al. Initial sequencing and analysis of the human genome. Nature. 2001;409:860–921.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Consortium CSaA. Initial sequence of the chimpanzee genome and comparison with the human genome. Nature. 2005;437:69–87.

    Article  Google Scholar 

  4. 4.

    Gibbs RA, Rogers J, Katze MG, Bumgarner R, Weinstock GM, Mardis ER, et al. Evolutionary and biomedical insights from the rhesus macaque genome. Science. 2007;316:222–34.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    King MC, Wilson AC. Evolution at two levels in humans and chimpanzees. Science. 1975;188:107–16.

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Nuttall. Blood immunity and blood relationships. Cambridge, UK: Cambridge University Press; 1904.

    Google Scholar 

  7. 7.

    Carroll SB. Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution. Cell. 2008;134:25–36.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Enard W. Human evolution: enhancing the brain. Curr Biol. 2015;25:R421–3.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Franchini LF, Pollard KS. Can a few non-coding mutations make a human brain? Bioessays. 2015;37:1054–61.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Franchini LF, Pollard KS. Genomic approaches to studying human-specific developmental traits. Development. 2015;142:3100–12.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Reilly SK, Noonan JP. Evolution of gene regulation in humans. Annu Rev Genomics Hum Genet. 2016;17:45–67.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Silver DL. Genomic divergence and brain evolution: How regulatory DNA influences development of the cerebral cortex. Bioessays. 2016;38:162–71.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Enard W. The molecular basis of human brain evolution. Curr Biol. 2016;26:R1109–R17.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Hubisz MJ, Pollard KS. Exploring the genesis and functions of Human Accelerated Regions sheds light on their role in human evolution. Curr Opin Genet Dev. 2014;29:15–21.

    CAS  Article  PubMed  Google Scholar 

  15. 15.

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

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Preuss TM. Human brain evolution: from gene discovery to phenotype discovery. Proc Natl Acad Sci U S A. 2012;109 Suppl 1:10709–16.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Vallender EJ, Mekel-Bobrov N, Lahn BT. Genetic basis of human brain evolution. Trends Neurosci. 2008;31:637–44.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Pollard KS, Salama SR, Lambert N, Lambot MA, Coppens S, Pedersen JS, et al. An RNA gene expressed during cortical development evolved rapidly in humans. Nature. 2006;443:167–72.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Pollard KS, Salama SR, King B, Kern AD, Dreszer T, Katzman S, et al. Forces shaping the fastest evolving regions in the human genome. PLoS Genet. 2006;2, e168.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Bird CP, Stranger BE, Liu M, Thomas DJ, Ingle CE, Beazley C, et al. Fast-evolving noncoding sequences in the human genome. Genome Biol. 2007;8:R118.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Bush EC, Lahn BT. A genome-wide screen for noncoding elements important in primate evolution. BMC Evol Biol. 2008;8:17.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Lowe CB, Haussler D. 29 mammalian genomes reveal novel exaptations of mobile elements for likely regulatory functions in the human genome. PLoS One. 2012;7, e43128.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

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

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Pennacchio LA, Ahituv N, Moses AM, Prabhakar S, Nobrega MA, Shoukry M, et al. In vivo enhancer analysis of human conserved non-coding sequences. Nature. 2006;444:499–502.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Visel A, Blow MJ, Li Z, Zhang T, Akiyama JA, Holt A, et al. ChIP-seq accurately predicts tissue-specific activity of enhancers. Nature. 2009;457:854–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Visel A, Prabhakar S, Akiyama JA, Shoukry M, Lewis KD, Holt A, et al. Ultraconservation identifies a small subset of extremely constrained developmental enhancers. Nat Genet. 2008;40:158–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Dickel DE, Visel A, Pennacchio LA. Functional anatomy of distant-acting mammalian enhancers. Philos Trans R Soc Lond B Biol Sci. 2013;368:20120359.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. 28.

    Yunis JJ, Prakash O. The origin of man: a chromosomal pictorial legacy. Science. 1982;215:1525–30.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Marques-Bonet T, Eichler EE. The evolution of human segmental duplications and the core duplicon hypothesis. Cold Spring Harb Symp Quant Biol. 2009;74:355–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Sudmant PH, Huddleston J, Catacchio CR, Malig M, Hillier LW, Baker C, et al. Evolution and diversity of copy number variation in the great ape lineage. Genome Res. 2013;23:1373–82.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  31. 31.

    Dennis MY, Eichler EE. Human adaptation and evolution by segmental duplication. Curr Opin Genet Dev. 2016;41:44–52.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Kostka D, Hahn MW, Pollard KS. Noncoding sequences near duplicated genes evolve rapidly. Genome Biol Evol. 2010;2:518–33.

    Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Lupianez DG, Kraft K, Heinrich V, Krawitz P, Brancati F, Klopocki E, et al. Disruptions of topological chromatin domains cause pathogenic rewiring of gene-enhancer interactions. Cell. 2015;161:1012–25.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Symmons O, Pan L, Remeseiro S, Aktas T, Klein F, Huber W, et al. The Shh topological domain facilitates the action of remote enhancers by reducing the effects of genomic distances. Dev Cell. 2016;39:529–43.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Consortium EP. An integrated encyclopedia of DNA elements in the human genome. Nature. 2012;489:57–74.

    Article  Google Scholar 

  36. 36.

    Ponting CP. Biological function in the twilight zone of sequence conservation. BMC Biol. 2017;15:71.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Capra JA, Erwin GD, McKinsey G, Rubenstein JL, Pollard KS. Many human accelerated regions are developmental enhancers. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130025.

    Article  PubMed  PubMed Central  Google Scholar 

  38. 38.

    Gittelman RM, Hun E, Ay F, Madeoy J, Pennacchio L, Noble WS, et al. Comprehensive identification and analysis of human accelerated regulatory DNA. Genome Res. 2015;25:1245–55.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Dong X, Wang X, Zhang F, Tian W. Genome-wide identification of regulatory sequences undergoing accelerated evolution in the human genome. Mol Biol Evol. 2016;33:2565–75.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Consortium EP, Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, et al. Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007;447:799–816.

    Article  Google Scholar 

  41. 41.

    Cotney J, Leng J, Yin J, Reilly SK, DeMare LE, Emera D, et al. The evolution of lineage-specific regulatory activities in the human embryonic limb. Cell. 2013;154:185–96.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Reilly SK, Yin J, Ayoub AE, Emera D, Leng J, Cotney J, et al. Evolutionary genomics. Evolutionary changes in promoter and enhancer activity during human corticogenesis. Science. 2015;347:1155–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Vermunt MW, Tan SC, Castelijns B, Geeven G, Reinink P, de Bruijn E, et al. Epigenomic annotation of gene regulatory alterations during evolution of the primate brain. Nat Neurosci. 2016;19:494–503.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Prescott SL, Srinivasan R, Marchetto MC, Grishina I, Narvaiza I, Selleri L, et al. Enhancer divergence and cis-regulatory evolution in the human and chimp neural crest. Cell. 2015;163:68–83.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Bajpai R, Chen DA, Rada-Iglesias A, Zhang J, Xiong Y, Helms J, et al. CHD7 cooperates with PBAF to control multipotent neural crest formation. Nature. 2010;463:958–62.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Rada-Iglesias A, Bajpai R, Swigut T, Brugmann SA, Flynn RA, Wysocka J. A unique chromatin signature uncovers early developmental enhancers in humans. Nature. 2011;470:279–83.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Doolittle WF. Is junk DNA bunk? A critique of ENCODE. Proc Natl Acad Sci U S A. 2013;110:5294–300.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Graur D, Zheng Y, Price N, Azevedo RB, Zufall RA, Elhaik E. On the immortality of television sets: ‘function’ in the human genome according to the evolution-free gospel of ENCODE. Genome Biol Evol. 2013;5:578–90.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Sakabe NJ, Nobrega MA. Beyond the ENCODE project: using genomics and epigenomics strategies to study enhancer evolution. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130022.

    Article  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Rossant J. Mouse and human blastocyst-derived stem cells: vive les differences. Development. 2015;142:9–12.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Florio M, Albert M, Taverna E, Namba T, Brandl H, Lewitus E, et al. Human-specific gene ARHGAP11B promotes basal progenitor amplification and neocortex expansion. Science. 2015;347(6229):1465–70.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    Charrier C, Joshi K, Coutinho-Budd J, Kim JE, Lambert N, de Marchena J, et al. Inhibition of SRGAP2 function by its human-specific paralogs induces neoteny during spine maturation. Cell. 2012;149:923–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Schreiweis C, Bornschein U, Burguiere E, Kerimoglu C, Schreiter S, Dannemann M, et al. Humanized Foxp2 accelerates learning by enhancing transitions from declarative to procedural performance. Proc Natl Acad Sci U S A. 2014;111:14253–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Enard W, Gehre S, Hammerschmidt K, Holter SM, Blass T, Somel M, et al. A humanized version of Foxp2 affects cortico-basal ganglia circuits in mice. Cell. 2009;137:961–71.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    McLean CY, Reno PL, Pollen AA, Bassan AI, Capellini TD, Guenther C, et al. Human-specific loss of regulatory DNA and the evolution of human-specific traits. Nature. 2011;471:216–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  56. 56.

    Boyd JL, Skove SL, Rouanet JP, Pilaz LJ, Bepler T, Gordan R, et al. Human-chimpanzee differences in a FZD8 enhancer alter cell-cycle dynamics in the developing neocortex. Curr Biol. 2015;25:772–9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  57. 57.

    Kamm GB, Lopez-Leal R, Lorenzo JR, Franchini LF. A fast-evolving human NPAS3 enhancer gained reporter expression in the developing forebrain of transgenic mice. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130019.

    Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Kamm GB, Pisciottano F, Kliger R, Franchini LF. The developmental brain gene NPAS3 contains the largest number of accelerated regulatory sequences in the human genome. Mol Biol Evol. 2013;30:1088–102.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Prabhakar S, Visel A, Akiyama JA, Shoukry M, Lewis KD, Holt A, et al. Human-specific gain of function in a developmental enhancer. Science. 2008;321:1346–50.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  60. 60.

    Doan RN, Bae BI, Cubelos B, Chang C, Hossain AA, Al-Saad S, et al. Mutations in human accelerated regions disrupt cognition and social behavior. Cell. 2016;167:341–54. e12.

    CAS  Article  PubMed  Google Scholar 

  61. 61.

    Melnikov A, Murugan A, Zhang X, Tesileanu T, Wang L, Rogov P, et al. Systematic dissection and optimization of inducible enhancers in human cells using a massively parallel reporter assay. Nat Biotechnol. 2012;30:271–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  62. 62.

    Patwardhan RP, Hiatt JB, Witten DM, Kim MJ, Smith RP, May D, et al. Massively parallel functional dissection of mammalian enhancers in vivo. Nat Biotechnol. 2012;30:265–70.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  63. 63.

    Fortna A, Kim Y, MacLaren E, Marshall K, Hahn G, Meltesen L, et al. Lineage-specific gene duplication and loss in human and great ape evolution. PLoS Biol. 2004;2, E207.

    Article  PubMed  PubMed Central  Google Scholar 

  64. 64.

    Davis JM, Searles VB, Anderson N, Keeney J, Dumas L, Sikela JM. DUF1220 dosage is linearly associated with increasing severity of the three primary symptoms of autism. PLoS Genet. 2014;10, e1004241.

    Article  PubMed  PubMed Central  Google Scholar 

  65. 65.

    Dennis MY, Nuttle X, Sudmant PH, Antonacci F, Graves TA, Nefedov M, et al. Evolution of human-specific neural SRGAP2 genes by incomplete segmental duplication. Cell. 2012;149:912–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  66. 66.

    Xu K, Schadt EE, Pollard KS, Roussos P, Dudley JT. Genomic and network patterns of schizophrenia genetic variation in human evolutionary accelerated regions. Mol Biol Evol. 2015;32:1148–60.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  67. 67.

    Cubelos B, Sebastian-Serrano A, Beccari L, Calcagnotto ME, Cisneros E, Kim S, et al. Cux1 and Cux2 regulate dendritic branching, spine morphology, and synapses of the upper layer neurons of the cortex. Neuron. 2010;66:523–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Fullwood MJ, Han Y, Wei CL, Ruan X, Ruan Y. Chromatin interaction analysis using paired-end tag sequencing. Curr Protoc Mol Biol. 2010;Chapter 21:Unit 21 15 1–25.

  69. 69.

    Grueber WB, Jan LY, Jan YN. Different levels of the homeodomain protein cut regulate distinct dendrite branching patterns of Drosophila multidendritic neurons. Cell. 2003;112:805–18.

    CAS  Article  PubMed  Google Scholar 

  70. 70.

    Li G, Fullwood MJ, Xu H, Mulawadi FH, Velkov S, Vega V, et al. ChIA-PET tool for comprehensive chromatin interaction analysis with paired-end tag sequencing. Genome Biol. 2010;11:R22.

    Article  PubMed  PubMed Central  Google Scholar 

  71. 71.

    Cantsilieris S, Stessman HA, Shendure J, Eichler EE. Targeted Capture and high-throughput sequencing using molecular inversion probes (MIPs). Methods Mol Biol. 2017;1492:95–106.

    Article  PubMed  PubMed Central  Google Scholar 

  72. 72.

    Domene S, Bumaschny VF, de Souza FS, Franchini LF, Nasif S, Low MJ, et al. Enhancer turnover and conserved regulatory function in vertebrate evolution. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130027.

    Article  PubMed  PubMed Central  Google Scholar 

  73. 73.

    Harmston N, Baresic A, Lenhard B. The mystery of extreme non-coding conservation. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130021.

    Article  PubMed  PubMed Central  Google Scholar 

  74. 74.

    Maeso I, Irimia M, Tena JJ, Casares F, Gomez-Skarmeta JL. Deep conservation of cis-regulatory elements in metazoans. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130020.

    Article  PubMed  PubMed Central  Google Scholar 

  75. 75.

    de Souza FS, Franchini LF, Rubinstein M. Exaptation of transposable elements into novel cis-regulatory elements: is the evidence always strong? Mol Biol Evol. 2013;30:1239–51.

    Article  PubMed  PubMed Central  Google Scholar 

  76. 76.

    Franchini LF, Lopez-Leal R, Nasif S, Beati P, Gelman DM, Low MJ, et al. Convergent evolution of two mammalian neuronal enhancers by sequential exaptation of unrelated retroposons. Proc Natl Acad Sci U S A. 2011;108(37):15270–5.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  77. 77.

    Frankel N, Davis GK, Vargas D, Wang S, Payre F, Stern DL. Phenotypic robustness conferred by apparently redundant transcriptional enhancers. Nature. 2010;466:490–3.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  78. 78.

    Lam DD, de Souza FS, Nasif S, Yamashita M, Lopez-Leal R, Otero-Corchon V, et al. Partially redundant enhancers cooperatively maintain Mammalian pomc expression above a critical functional threshold. PLoS Genet. 2015;11, e1004935.

    Article  PubMed  PubMed Central  Google Scholar 

  79. 79.

    Levine M. Transcriptional enhancers in animal development and evolution. Curr Biol. 2010;20:R754–63.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  80. 80.

    Rubinstein M, de Souza FS. Evolution of transcriptional enhancers and animal diversity. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130017.

    Article  PubMed  PubMed Central  Google Scholar 

  81. 81.

    Gilbert LA, Horlbeck MA, Adamson B, Villalta JE, Chen Y, Whitehead EH, et al. Genome-scale CRISPR-mediated control of gene repression and activation. Cell. 2014;159:647–61.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  82. 82.

    Zheng Q, Cai X, Tan MH, Schaffert S, Arnold CP, Gong X, et al. Precise gene deletion and replacement using the CRISPR/Cas9 system in human cells. Biotechniques. 2014;57:115–24.

    CAS  Article  PubMed  Google Scholar 

  83. 83.

    Cheng Z, Ventura M, She X, Khaitovich P, Graves T, Osoegawa K, et al. A genome-wide comparison of recent chimpanzee and human segmental duplications. Nature. 2005;437:88–93.

    CAS  Article  PubMed  Google Scholar 

  84. 84.

    Vietri Rudan M, Barrington C, Henderson S, Ernst C, Odom DT, Tanay A, et al. Comparative Hi-C reveals that CTCF underlies evolution of chromosomal domain architecture. Cell Rep. 2015;10:1297–309.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  85. 85.

    Acemel RD, Maeso I, Gomez-Skarmeta JL. Topologically associated domains: a successful scaffold for the evolution of gene regulation in animals. Wiley Interdiscip Rev Dev Biol. 2017;6.

  86. 86.

    Maeso I, Acemel RD, Gomez-Skarmeta JL. Cis-regulatory landscapes in development and evolution. Curr Opin Genet Dev. 2017;43:17–22.

    CAS  Article  PubMed  Google Scholar 

  87. 87.

    Pollard KS, Hubisz MJ, Rosenbloom KR, Siepel A. Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res. 2010;20:110–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  88. 88.

    Kostka D, Hubisz MJ, Siepel A, Pollard KS. The role of GC-biased gene conversion in shaping the fastest evolving regions of the human genome. Mol Biol Evol. 2012;29:1047–57.

    CAS  Article  PubMed  Google Scholar 

  89. 89.

    Katzman S, Kern AD, Pollard KS, Salama SR, Haussler D. GC-biased evolution near human accelerated regions. PLoS Genet. 2010;6, e1000960.

    Article  PubMed  PubMed Central  Google Scholar 

  90. 90.

    Sumiyama K, Saitou N. Loss-of-function mutation in a repressor module of human-specifically activated enhancer HACNS1. Mol Biol Evol. 2011;28:3005–7.

    CAS  Article  PubMed  Google Scholar 

  91. 91.

    Sikela JM. The jewels of our genome: the search for the genomic changes underlying the evolutionarily unique capacities of the human brain. PLoS Genet. 2006;2, e80.

    Article  PubMed  PubMed Central  Google Scholar 

  92. 92.

    Varki A, Altheide TK. Comparing the human and chimpanzee genomes: searching for needles in a haystack. Genome Res. 2005;15:1746–58.

    CAS  Article  PubMed  Google Scholar 

  93. 93.

    Herculano-Houzel S. The remarkable, yet not extraordinary, human brain as a scaled-up primate brain and its associated cost. Proc Natl Acad Sci U S A. 2012;109 Suppl 1:10661–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  94. 94.

    Preuss TM. Primate brain evolution in phylogenetic context. In: Kaas J, editor. Evolution of nervous systems. 4. Oxford, UK: Elsevier; 2007. p. 1–34.

  95. 95.

    Semendeferi K, Lu A, Schenker N, Damasio H. Humans and great apes share a large frontal cortex. Nat Neurosci. 2002;5:272–6.

    CAS  Article  PubMed  Google Scholar 

  96. 96.

    Gabi M, Neves K, Masseron C, Ribeiro PF, Ventura-Antunes L, Torres L, et al. No relative expansion of the number of prefrontal neurons in primate and human evolution. Proc Natl Acad Sci U S A. 2016;113:9617–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  97. 97.

    Semendeferi K, Teffer K, Buxhoeveden DP, Park MS, Bludau S, Amunts K, et al. Spatial organization of neurons in the frontal pole sets humans apart from great apes. Cereb Cortex. 2011;21:1485–97.

    Article  PubMed  Google Scholar 

  98. 98.

    Bianchi S, Stimpson CD, Duka T, Larsen MD, Janssen WG, Collins Z, et al. Synaptogenesis and development of pyramidal neuron dendritic morphology in the chimpanzee neocortex resembles humans. Proc Natl Acad Sci U S A. 2013;110 Suppl 2:10395–401.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  99. 99.

    Spocter MA, Hopkins WD, Barks SK, Bianchi S, Hehmeyer AE, Anderson SM, et al. Neuropil distribution in the cerebral cortex differs between humans and chimpanzees. J Comp Neurol. 2012;520:2917–29.

    Article  PubMed  PubMed Central  Google Scholar 

  100. 100.

    Lindblad-Toh K, Garber M, Zuk O, Lin MF, Parker BJ, Washietl S, et al. A high-resolution map of human evolutionary constraint using 29 mammals. Nature. 2011;478:476–82.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  101. 101.

    Cook EC, Nelson JK, Sorrentino V, Koenis D, Moeton M, Scheij S, et al. Identification of the ER-resident E3 ubiquitin ligase RNF145 as a novel LXR-regulated gene. PLoS One. 2017;12, e0172721.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


This work was supported by grants from the Agencia Nacional de Promoción Cientıífica y Tecnológica (PICT2013-1642; PICT 2015–1726) to LFF and by Gladstone Institutes to KSP.

Author information



Corresponding author

Correspondence to Katherine S. Pollard.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Franchini, L.F., Pollard, K.S. Human evolution: the non-coding revolution. BMC Biol 15, 89 (2017).

Download citation

  • Published:

  • DOI:


  • Human Accelerated Regions (HARs)
  • DNase I-hypersensitive Sites (DHSs)
  • Enhancer Hijacking
  • Chimpanzee
  • H3K27ac Signal