Machine learning dissection of human accelerated regions in primate neurodevelopment

Large-scale genomic analysis of the domestic dog informs biological discovery

Recent advances in genomics, coupled with a unique population structure and remarkable levels of variation, have propelled the domestic dog to new levels as a system for understanding fundamental principles in mammalian biology. Central to this advance are more than 350 recognized breeds, each a closed population that has undergone selection for unique features. Genetic variation in the domestic dog is particularly well characterized compared with other domestic mammals, with almost 3000 high-coverage genomes publicly available. Importantly, as the number of sequenced genomes increases, new avenues for analysis are becoming available. Herein, we discuss recent discoveries in canine genomics regarding behavior, morphology, and disease susceptibility. We explore the limitations of current data sets for variant interpretation, tradeoffs between sequencing strategies, and the burgeoning role of long-read genomes for capturing structural variants. In addition, we consider how large-scale collections of whole-genome sequence data drive rare variant discovery and assess the geographic distribution of canine diversity, which identifies Asia as a major source of missing variation. Finally, we review recent comparative genomic analyses that will facilitate annotation of the noncoding genome in dogs.

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

Human Accelerated Regions regulate gene networks implicated in apical-to-basal neural progenitor fate transitions

The evolution of the human cerebral cortex involved modifications in the composition and proliferative potential of the neural stem cell (NSC) niche during brain development. Human Accelerated Regions (HARs) exhibit a significant excess of human-specific sequence changes and have been implicated in human brain evolution. Multiple studies support that HARs include neurodevelopmental enhancers with novel activities in humans, but their biological functions in NSCs have not been empirically assessed at scale. Here we conducted a direct-capture Perturb-seq screen repressing 180 neurodevelopmentally active HARs in human iPSC-derived NSCs with single-cell transcriptional readout. After profiling >188,000 NSCs, we identified a set of HAR perturbations with convergent transcriptional effects on gene networks involved in NSC apicobasal polarity, a cellular process whose precise regulation is critical to the developmental emergence of basal radial glia (bRG), a progenitor population that is expanded in humans. Across multiple HAR perturbations, we found convergent dysregulation of specific apicobasal polarity and adherens junction regulators, including PARD3, ABI2, SETD2 , and PCM1 . We found that the repression of one candidate from the screen, HAR181, as well as its target gene CADM1 , disrupted apical PARD3 localization and NSC rosette formation. Our findings reveal interconnected roles for HARs in NSC biology and cortical development and link specific HARs to processes implicated in human cortical expansion.

Genomic, molecular, and cellular divergence of the human brain

While many core biological processes are conserved across species, the human brain has evolved with unique capacities. Current understanding of the neurobiological mechanisms that endow human traits as well as associated vulnerabilities remains limited. However, emerging data have illuminated species divergence in DNA elements and genome organization, in molecular, morphological, and functional features of conserved neural cell types, as well as temporal differences in brain development. Here, we summarize recent data on unique features of the human brain and their complex implications for the study and treatment of brain diseases. We also consider key outstanding questions in the field and discuss the technologies and foundational knowledge that will be required to accelerate understanding of human neurobiology.

Advances in computational and experimental approaches for deciphering transcriptional regulatory networks: Understanding the roles of cis-regulatory elements is essential, and recent research utilizing MPRAs, STARR-seq, CRISPR-Cas9, and machine learning has yielded valuable insights

Understanding the influence of cis-regulatory elements on gene regulation poses numerous challenges given complexities stemming from variations in transcription factor (TF) binding, chromatin accessibility, structural constraints, and cell-type differences. This review discusses the role of gene regulatory networks in enhancing understanding of transcriptional regulation and covers construction methods ranging from expression-based approaches to supervised machine learning. Additionally, key experimental methods, including MPRAs and CRISPR-Cas9-based screening, which have significantly contributed to understanding TF binding preferences and cis-regulatory element functions, are explored. Lastly, the potential of machine learning and artificial intelligence to unravel cis-regulatory logic is analyzed. These computational advances have far-reaching implications for precision medicine, therapeutic target discovery, and the study of genetic variations in health and disease.

Accelerated evolution in the human lineage led to gain and loss of transcriptional enhancers in the RBFOX1 locus

A long-standing goal of evolutionary biology is to decode how changes in gene regulatory networks contribute to human-specific traits. Human accelerated regions (HARs) are prime candidates for driving gene regulatory modifications in human development. The RBFOX1 locus is densely populated with HARs, providing a set of potential regulatory elements that could have changed its expression in the human lineage. Here, we examined the role of RBFOX1-HARs using transgenic zebrafish reporter assays and identified 15 transcriptional enhancers that are active in the developing nervous system, 9 of which displayed differential activity between the human and chimpanzee sequences. The engineered loss of two selected RBFOX1-HARs in knockout mouse models modified Rbfox1 expression at specific developmental stages and tissues in the brain, influencing the expression and splicing of a high number of Rbfox1 target genes. Our results provided insight into the spatial and temporal changes in gene expression driven by RBFOX1-HARs.

CpG island turnover events predict evolutionary changes in enhancer activity

BACKGROUND

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution.

RESULTS

We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon.

CONCLUSIONS

Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

A molecular and cellular perspective on human brain evolution and tempo

The evolution of the modern human brain was accompanied by distinct molecular and cellular specializations, which underpin our diverse cognitive abilities but also increase our susceptibility to neurological diseases. These features, some specific to humans and others shared with related species, manifest during different stages of brain development. In this multi-stage process, neural stem cells proliferate to produce a large and diverse progenitor pool, giving rise to excitatory or inhibitory neurons that integrate into circuits during further maturation. This process unfolds over varying time scales across species and has progressively become slower in the human lineage, with differences in tempo correlating with differences in brain size, cell number and diversity, and connectivity. Here we introduce the terms 'bradychrony' and 'tachycrony' to describe slowed and accelerated developmental tempos, respectively. We review how recent technical advances across disciplines, including advanced engineering of in vitro models, functional comparative genetics and high-throughput single-cell profiling, are leading to a deeper understanding of how specializations of the human brain arise during bradychronic neurodevelopment. Emerging insights point to a central role for genetics, gene-regulatory networks, cellular innovations and developmental tempo, which together contribute to the establishment of human specializations during various stages of neurodevelopment and at different points in evolution.

Massively parallel reporter assays and mouse transgenic assays provide complementary information about neuronal enhancer activity

Genetic studies find hundreds of thousands of noncoding variants associated with psychiatric disorders. Massively parallel reporter assays (MPRAs) and in vivo transgenic mouse assays can be used to assay the impact of these variants. However, the relevance of MPRAs to in vivo function is unknown and transgenic assays suffer from low throughput. Here, we studied the utility of combining the two assays to study the impact of non-coding variants. We carried out an MPRA on over 50,000 sequences derived from enhancers validated in transgenic mouse assays and from multiple fetal neuronal ATAC-seq datasets. We also tested over 20,000 variants, including synthetic mutations in highly active neuronal enhancers and 177 common variants associated with psychiatric disorders. Variants with a high impact on MPRA activity were further tested in mice. We found a strong and specific correlation between MPRA and mouse neuronal enhancer activity including changes in neuronal enhancer activity in mouse embryos for variants with strong MPRA effects. Mouse assays also revealed pleiotropic variant effects that could not be observed in MPRA. Our work provides a large catalog of functional neuronal enhancers and variant effects and highlights the effectiveness of combining MPRAs and mouse transgenic assays.

A human-specific enhancer fine-tunes radial glia potency and corticogenesis

Humans evolved an extraordinarily expanded and complex cerebral cortex, associated with developmental and gene regulatory modifications 1-3 . Human accelerated regions (HARs) are highly conserved genomic sequences with human-specific nucleotide substitutions. Although there are thousands of annotated HARs, their functional contribution to human-specific cortical development is largely unknown 4,5 . HARE5 is a HAR transcriptional enhancer of the WNT signaling receptor Frizzled8 (FZD8) active during brain development 6 . Here, using genome-edited mouse and primate models, we demonstrate that human (Hs) HARE5 fine-tunes cortical development and connectivity by controlling the proliferative and neurogenic capacity of neural progenitor cells (NPCs). Hs-HARE5 knock-in mice have significantly enlarged neocortices containing more neurons. By measuring neural dynamics in vivo we show these anatomical features correlate with increased functional independence between cortical regions. To understand the underlying developmental mechanisms, we assess progenitor fate using live imaging, lineage analysis, and single-cell RNA sequencing. This reveals Hs-HARE5 modifies radial glial progenitor behavior, with increased self-renewal at early developmental stages followed by expanded neurogenic potential. We use genome-edited human and chimpanzee (Pt) NPCs and cortical organoids to assess the relative enhancer activity and function of Hs-HARE5 and Pt-HARE5. Using these orthogonal strategies we show four human-specific variants in HARE5 drive increased enhancer activity which promotes progenitor proliferation. These findings illustrate how small changes in regulatory DNA can directly impact critical signaling pathways and brain development. Our study uncovers new functions for HARs as key regulatory elements crucial for the expansion and complexity of the human cerebral cortex.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using assay for transposase-accessible chromatin coupled to self-transcribing active regulatory region (ATAC-STARR) sequencing. In addition to thousands of cis changes, we discover an unexpected number (∼10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie ∼37% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-in regulatory differences between species.

The growth factor EPIREGULIN promotes basal progenitor cell proliferation in the developing neocortex

Neocortex expansion during evolution is linked to higher numbers of neurons, which are thought to result from increased proliferative capacity and neurogenic potential of basal progenitor cells during development. Here, we show that EREG, encoding the growth factor EPIREGULIN, is expressed in the human developing neocortex and in gorilla cerebral organoids, but not in the mouse neocortex. Addition of EPIREGULIN to the mouse neocortex increases proliferation of basal progenitor cells, whereas EREG ablation in human cortical organoids reduces proliferation in the subventricular zone. Treatment of cortical organoids with EPIREGULIN promotes a further increase in proliferation of gorilla but not of human basal progenitor cells. EPIREGULIN competes with the epidermal growth factor (EGF) to promote proliferation, and inhibition of the EGF receptor abrogates the EPIREGULIN-mediated increase in basal progenitor cells. Finally, we identify putative cis-regulatory elements that may contribute to the observed inter-species differences in EREG expression. Our findings suggest that species-specific regulation of EPIREGULIN expression may contribute to the increased neocortex size of primates by providing a tunable pro-proliferative signal to basal progenitor cells in the subventricular zone.

Timing neurogenesis: a clock or an algorithm?

Emerging evidence supports the existence of dedicated molecular mechanisms under evolutionary selection to control time during neurogenesis. Here, we briefly review these mechanisms and discuss a potentially useful conceptual framework inspired by computer science to think about how these biological mechanisms operate during brain development and evolution.

Functional associations of evolutionarily recent human genes exhibit sensitivity to the 3D genome landscape and disease

Genome organization is intricately tied to regulating genes and associated cell fate decisions. In this study, we examine the positioning and functional significance of human genes, grouped by their evolutionary age, within the 3D organization of the genome. We reveal that genes of different evolutionary origin have distinct positioning relationships with both domains and loop anchors, and remarkably consistent relationships with boundaries across cell types. While the functional associations of each group of genes are primarily cell type-specific, such associations of conserved genes maintain greater stability across 3D genomic features and disease than recently evolved genes. Furthermore, the expression of these genes across various tissues follows an evolutionary progression, such that RNA levels increase from young genes to ancient genes. Thus, the distinct relationships of gene evolutionary age, function, and positioning within 3D genomic features contribute to tissue-specific gene regulation in development and disease.

Characterization of enhancer activity in early human neurodevelopment using Massively Parallel Reporter Assay (MPRA) and forebrain organoids

Regulation of gene expression through enhancers is one of the major processes shaping the structure and function of the human brain during development. High-throughput assays have predicted thousands of enhancers involved in neurodevelopment, and confirming their activity through orthogonal functional assays is crucial. Here, we utilized Massively Parallel Reporter Assays (MPRAs) in stem cells and forebrain organoids to evaluate the activity of ~ 7000 gene-linked enhancers previously identified in human fetal tissues and brain organoids. We used a Gaussian mixture model to evaluate the contribution of background noise in the measured activity signal to confirm the activity of ~ 35% of the tested enhancers, with most showing temporal-specific activity, suggesting their evolving role in neurodevelopment. The temporal specificity was further supported by the correlation of activity with gene expression. Our findings provide a valuable gene regulatory resource to the scientific community.

Development and evolution of the primate neocortex from a progenitor cell perspective

The generation of neurons in the developing neocortex is a major determinant of neocortex size. Crucially, the increase in cortical neuron numbers in the primate lineage, notably in the upper-layer neurons, contributes to increased cognitive abilities. Here, we review major evolutionary changes affecting the apical progenitors in the ventricular zone and focus on the key germinal zone constituting the foundation of neocortical neurogenesis in primates, the outer subventricular zone (OSVZ). We summarize characteristic features of the OSVZ and its key stem cell type, the basal (or outer) radial glia. Next, we concentrate on primate-specific and human-specific genes, expressed in OSVZ-progenitors, the ability of which to amplify these progenitors by targeting the regulation of the cell cycle ultimately underlies the evolutionary increase in upper-layer neurons. Finally, we address likely differences in neocortical development between present-day humans and Neanderthals that are based on human-specific amino acid substitutions in proteins operating in cortical progenitors.

The role of transposable elements in human evolution and methods for their functional analysis: current status and future perspectives

Transposable elements (TEs) are mobile DNA sequences that can insert themselves into various locations within the genome, causing mutations that may provide advantages or disadvantages to individuals and species. The insertion of TEs can result in genetic variation that may affect a wide range of human traits including genetic disorders. Understanding the role of TEs in human biology is crucial for both evolutionary and medical research. This review discusses the involvement of TEs in human traits and disease susceptibility, as well as methods for functional analysis of TEs.

Evolutionary innovation in conserved regulatory elements across the mammalian tree of life

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Evolution of enhancer sequences can alter target gene expression without causing detrimental misexpression in other contexts. It has long been thought that this modularity allows evolutionary changes in enhancers to escape pleiotropic constraints, which is especially important for evolutionary constrained developmental patterning genes. However, there is still little data supporting this hypothesis. Here we identified signatures of accelerated evolution in conserved enhancer elements across the mammalian phylogeny. We found that pleiotropic genes involved in gene regulatory and developmental processes were enriched for accelerated sequence evolution within their enhancer elements. These genes were associated with an excess number of enhancers compared to other genes, and due to this they exhibit a substantial degree of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. We studied one acceleration event in depth and found that its sequence evolution led to the emergence of a new enhancer activity domain that may be involved in the evolution of digit reduction in hoofed mammals. Our results provide tangible evidence that enhancer evolution has been a frequent contributor to modifications involving constrained developmental signaling genes in mammals.

Genetics of human brain development

Brain development in humans is achieved through precise spatiotemporal genetic control, the mechanisms of which remain largely elusive. Recently, integration of technological advances in human stem cell-based modelling with genome editing has emerged as a powerful platform to establish causative links between genotypes and phenotypes directly in the human system. Here, we review our current knowledge of complex genetic regulation of each key step of human brain development through the lens of evolutionary specialization and neurodevelopmental disorders and highlight the use of human stem cell-derived 2D cultures and 3D brain organoids to investigate human-enriched features and disease mechanisms. We also discuss opportunities and challenges of integrating new technologies to reveal the genetic architecture of human brain development and disorders.

Mapping cis- and trans-regulatory target genes of human-specific deletions

Deletion of functional sequence is predicted to represent a fundamental mechanism of molecular evolution1,2. Comparative genetic studies of primates2,3 have identified thousands of human-specific deletions (hDels), and the cis-regulatory potential of short (≤31 base pairs) hDels has been assessed using reporter assays4. However, how structural variant-sized (≥50 base pairs) hDels influence molecular and cellular processes in their native genomic contexts remains unexplored. Here, we design genome-scale libraries of single-guide RNAs targeting 7.2 megabases of sequence in 6,358 hDels and present a systematic CRISPR interference (CRISPRi) screening approach to identify hDels that modify cellular proliferation in chimpanzee pluripotent stem cells. By intersecting hDels with chromatin state features and performing single-cell CRISPRi (Perturb-seq) to identify their cis- and trans-regulatory target genes, we discovered 19 hDels controlling gene expression. We highlight two hDels, hDel_2247 and hDel_585, with tissue-specific activity in the liver and brain, respectively. Our findings reveal a molecular and cellular role for sequences lost in the human lineage and establish a framework for functionally interrogating human-specific genetic variants.

MPRAbase: A Massively Parallel Reporter Assay Database

Massively parallel reporter assays (MPRAs) represent a set of high-throughput technologies that measure the functional effects of thousands of sequences/variants on gene regulatory activity. There are several different variations of MPRA technology and they are used for numerous applications, including regulatory element discovery, variant effect measurement, saturation mutagenesis, synthetic regulatory element generation or characterization of evolutionary gene regulatory differences. Despite their many designs and uses, there is no comprehensive database that incorporates the results of these experiments. To address this, we developed MPRAbase, a manually curated database that currently harbors 129 experiments, encompassing 17,718,677 elements tested across 35 cell types and 4 organisms. The MPRAbase web interface ( http://www.mprabase.com ) serves as a centralized user-friendly repository to download existing MPRA data for independent analysis and is designed with the ability to allow researchers to share their published data for rapid dissemination to the community.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Decoding DNA sequence-driven evolution of the human brain epigenome at cellular resolution

DNA-based evolutionary comparisons of regulatory genomic elements enable insight into functional changes, overcoming tissue inaccessibility. Here, we harnessed adult and fetal cortex single-cell ATAC-seq datasets to uncover DNA substitutions specific to the human and human-ancestral lineages within apes. We found that fetal microglia identity is evolutionarily divergent in all lineages, whereas other cell types are conserved. Using multiomic datasets, we further identified genes linked to multiple lineage-divergent gene regulatory elements and implicated biological pathways associated with these divergent features. We also uncovered patterns of transcription factor binding site evolution across lineages and identified expansion of bHLH-PAS factor targets in human-hominin lineages, and MEF2 factor targets in the ape lineage. Finally, conserved features were more enriched in brain disease variants, whereas there was no distinct enrichment on the human lineage compared to its ancestral lineages. Our study identifies major evolutionary patterns in the human brain epigenome at cellular resolution.

CpG island turnover events predict evolutionary changes in enhancer activity

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Three-dimensional genome rewiring in loci with human accelerated regions

Human accelerated regions (HARs) are conserved genomic loci that evolved at an accelerated rate in the human lineage and may underlie human-specific traits. We generated HARs and chimpanzee accelerated regions with an automated pipeline and an alignment of 241 mammalian genomes. Combining deep learning with chromatin capture experiments in human and chimpanzee neural progenitor cells, we discovered a significant enrichment of HARs in topologically associating domains containing human-specific genomic variants that change three-dimensional (3D) genome organization. Differential gene expression between humans and chimpanzees at these loci suggests rewiring of regulatory interactions between HARs and neurodevelopmental genes. Thus, comparative genomics together with models of 3D genome folding revealed enhancer hijacking as an explanation for the rapid evolution of HARs.

Genetic architecture of human brain evolution

Comparative studies of hominids have long sought to identify mutational events that shaped the evolution of the human nervous system. However, functional genetic differences are outnumbered by millions of nearly neutral mutations, and the developmental mechanisms underlying human nervous system specializations are difficult to model and incompletely understood. Candidate-gene studies have attempted to map select human-specific genetic differences to neurodevelopmental functions, but it remains unclear how to contextualize the relative effects of genes that are investigated independently. Considering these limitations, we discuss scalable approaches for probing the functional contributions of human-specific genetic differences. We propose that a systems-level view will enable a more quantitative and integrative understanding of the genetic, molecular and cellular underpinnings of human nervous system evolution.

Genetic Variation in Transcription Factor Binding Sites

The interaction between transcription factors (TFs) and DNA is the core process that determines the state of a cell's transcriptome [...].

Massively parallel characterization of psychiatric disorder-associated and cell-type-specific regulatory elements in the developing human cortex

Nucleotide changes in gene regulatory elements are important determinants of neuronal development and disease. Using massively parallel reporter assays in primary human cells from mid-gestation cortex and cerebral organoids, we interrogated the cis-regulatory activity of 102,767 sequences, including differentially accessible cell-type specific regions in the developing cortex and single-nucleotide variants associated with psychiatric disorders. In primary cells, we identified 46,802 active enhancer sequences and 164 disorder-associated variants that significantly alter enhancer activity. Activity was comparable in organoids and primary cells, suggesting that organoids provide an adequate model for the developing cortex. Using deep learning, we decoded the sequence basis and upstream regulators of enhancer activity. This work establishes a comprehensive catalog of functional gene regulatory elements and variants in human neuronal development.

Developmental mechanisms underlying the evolution of human cortical circuits

The brain of modern humans has evolved remarkable computational abilities that enable higher cognitive functions. These capacities are tightly linked to an increase in the size and connectivity of the cerebral cortex, which is thought to have resulted from evolutionary changes in the mechanisms of cortical development. Convergent progress in evolutionary genomics, developmental biology and neuroscience has recently enabled the identification of genomic changes that act as human-specific modifiers of cortical development. These modifiers influence most aspects of corticogenesis, from the timing and complexity of cortical neurogenesis to synaptogenesis and the assembly of cortical circuits. Mutations of human-specific genetic modifiers of corticogenesis have started to be linked to neurodevelopmental disorders, providing evidence for their physiological relevance and suggesting potential relationships between the evolution of the human brain and its sensitivity to specific diseases.

Human gene regulatory evolution is driven by the divergence of regulatory element function in both cis and trans

Gene regulatory divergence between species can result from cis-acting local changes to regulatory element DNA sequences or global trans-acting changes to the regulatory environment. Understanding how these mechanisms drive regulatory evolution has been limited by challenges in identifying trans-acting changes. We present a comprehensive approach to directly identify cis- and trans-divergent regulatory elements between human and rhesus macaque lymphoblastoid cells using ATAC-STARR-seq. In addition to thousands of cis changes, we discover an unexpected number (~10,000) of trans changes and show that cis and trans elements exhibit distinct patterns of sequence divergence and function. We further identify differentially expressed transcription factors that underlie >50% of trans differences and trace how cis changes can produce cascades of trans changes. Overall, we find that most divergent elements (67%) experienced changes in both cis and trans, revealing a substantial role for trans divergence-alone and together with cis changes-to regulatory differences between species.

Transcription factor binding sites are frequently under accelerated evolution in primates

Recent comparative genomic studies have identified many human accelerated elements (HARs) with elevated substitution rates in the human lineage. However, it remains unknown to what extent transcription factor binding sites (TFBSs) are under accelerated evolution in humans and other primates. Here, we introduce two pooling-based phylogenetic methods with dramatically enhanced sensitivity to examine accelerated evolution in TFBSs. Using these new methods, we show that more than 6000 TFBSs annotated in the human genome have experienced accelerated evolution in Hominini, apes, and Old World monkeys. Although these TFBSs individually show relatively weak signals of accelerated evolution, they collectively are more abundant than HARs. Also, we show that accelerated evolution in Pol III binding sites may be driven by lineage-specific positive selection, whereas accelerated evolution in other TFBSs might be driven by nonadaptive evolutionary forces. Finally, the accelerated TFBSs are enriched around developmental genes, suggesting that accelerated evolution in TFBSs may drive the divergence of developmental processes between primates.

Characterization of De Novo Promoter Variants in Autism Spectrum Disorder with Massively Parallel Reporter Assays

Autism spectrum disorder (ASD) is a common, complex, and highly heritable condition with contributions from both common and rare genetic variations. While disruptive, rare variants in protein-coding regions clearly contribute to symptoms, the role of rare non-coding remains unclear. Variants in these regions, including promoters, can alter downstream RNA and protein quantity; however, the functional impacts of specific variants observed in ASD cohorts remain largely uncharacterized. Here, we analyzed 3600 de novo mutations in promoter regions previously identified by whole-genome sequencing of autistic probands and neurotypical siblings to test the hypothesis that mutations in cases have a greater functional impact than those in controls. We leveraged massively parallel reporter assays (MPRAs) to detect transcriptional consequences of these variants in neural progenitor cells and identified 165 functionally high confidence de novo variants (HcDNVs). While these HcDNVs are enriched for markers of active transcription, disruption to transcription factor binding sites, and open chromatin, we did not identify differences in functional impact based on ASD diagnostic status.