A high-resolution enhancer atlas of the developing telencephalon

A massively parallel reporter assay library to screen short synthetic promoters in mammalian cells

Cellular responses to stimuli underpin discoveries in drug development, synthetic biology, and general life sciences. We introduce a library comprising 6144 synthetic promoters, each shorter than 250 bp, designed as transcriptional readouts of cellular stimulus responses in massively parallel reporter assay format. This library facilitates precise detection and amplification of transcriptional activity from our promoters, enabling the systematic development of tunable reporters with dynamic ranges of 50-100 fold. Our library proved functional in numerous cell lines and responsive to a variety of stimuli, including metabolites, mitogens, toxins, and pharmaceutical agents, generating robust and scalable reporters effective in screening assays, biomarkers, and synthetic circuits attuned to endogenous cellular activities. Particularly valuable in therapeutic development, our library excels in capturing candidate reporters to signals mediated by drug targets, a feature we illustrate across nine diverse G-protein coupled receptors (GPCRs), critical targets in drug development. We detail how this tool isolates and defines discrete signaling pathways associated with specific GPCRs, elucidating their transcriptional signatures. With its ease of implementation, broad utility, publicly available data, and comprehensive documentation, our library will be beneficial in synthetic biology, cellular engineering, ligand exploration, and drug development.

Synthetic Promoters in Gene Therapy: Design Approaches, Features and Applications

Gene therapy is a promising approach to the treatment of various inherited diseases, but its development is complicated by a number of limitations of the natural promoters used. The currently used strong ubiquitous natural promoters do not allow for the specificity of expression, while natural tissue-specific promoters have lowactivity. These limitations of natural promoters can be addressed by creating new synthetic promoters that achieve high levels of tissue-specific target gene expression. This review discusses recent advances in the development of synthetic promoters that provide a more precise regulation of gene expression. Approaches to the design of synthetic promoters are reviewed, including manual design and bioinformatic methods using machine learning. Examples of successful applications of synthetic promoters in the therapy of hereditary diseases and cancer are presented, as well as prospects for their clinical use.

Combinatorial transcription factor binding encodes cis-regulatory wiring of mouse forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to cis-regulatory elements, orchestrating transcriptional programs. Although studies of chromatin state and chromosomal interactions have demonstrated dynamic neurodevelopmental cis-regulatory landscapes, parallel understanding of TF interactions lags. To elucidate combinatorial TF binding driving mouse basal ganglia development, we integrated chromatin immunoprecipitation sequencing (ChIP-seq) for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and gene expression data, and functional enhancer assays. We identified sets of putative regulatory elements with shared TF binding (TF-pRE modules) that orchestrate distinct processes of GABAergic neurogenesis and suppress other cell fates. The majority of pREs were bound by one or two TFs; however, a small proportion were extensively bound. These sequences had exceptional evolutionary conservation and motif density, complex chromosomal interactions, and activity as in vivo enhancers. Our results provide insights into the combinatorial TF-pRE interactions that activate and repress expression programs during telencephalon neurogenesis and demonstrate the value of TF binding toward modeling developmental transcriptional wiring.

Nr2f1 enhancers have distinct functions in controlling Nr2f1 expression during cortical development

There is evidence that transcription factor (TF) encoding genes, which temporally control development in multiple cell types, can have tens of enhancers that regulate their expression. The NR2F1 TF developmentally promotes caudal and ventral cortical regional fates. Here, we epigenomically compared the activity of Nr2f1's enhancers during mouse cortical development with their activity in a transgenic assay. We identified at least six that are likely to be important in prenatal cortical development, with three harboring de novo mutants identified in ASD individuals. We chose to study the function of two of the most robust enhancers by deleting them singly or together. We found that they have distinct and overlapping functions in driving Nr2f1's regional and laminar expression in the developing cortex. Thus, these two enhancers, probably in combination with the others that we defined epigenetically, precisely tune Nr2f1's regional, cell type, and temporal expression during corticogenesis.

Neuronal enhancers fine-tune adaptive circuit plasticity

Neuronal activity-regulated gene expression plays a crucial role in sculpting neural circuits that underpin adaptive brain function. Transcriptional enhancers are now recognized as key components of gene regulation that orchestrate spatiotemporally precise patterns of gene transcription. We propose that the dynamics of enhancer activation uniquely position these genomic elements to finely tune activity-dependent cellular plasticity. Enhancer specificity and modularity can be exploited to gain selective genetic access to specific cell states, and the precise modulation of target gene expression within restricted cellular contexts enabled by targeted enhancer manipulation allows for fine-grained evaluation of gene function. Mounting evidence also suggests that enduring stimulus-induced changes in enhancer states can modify target gene activation upon restimulation, thereby contributing to a form of cell-wide metaplasticity. We advocate for focused exploration of activity-dependent enhancer function to gain new insight into the mechanisms underlying brain plasticity and cognitive dysfunction.

Mutagenesis Sensitivity Mapping of Human Enhancers In Vivo

Distant-acting enhancers are central to human development. However, our limited understanding of their functional sequence features prevents the interpretation of enhancer mutations in disease. Here, we determined the functional sensitivity to mutagenesis of human developmental enhancers in vivo. Focusing on seven enhancers active in the developing brain, heart, limb and face, we created over 1700 transgenic mice for over 260 mutagenized enhancer alleles. Systematic mutation of 12-basepair blocks collectively altered each sequence feature in each enhancer at least once. We show that 69% of all blocks are required for normal in vivo activity, with mutations more commonly resulting in loss (60%) than in gain (9%) of function. Using predictive modeling, we annotated critical nucleotides at base-pair resolution. The vast majority of motifs predicted by these machine learning models (88%) coincided with changes to in vivo function, and the models showed considerable sensitivity, identifying 59% of all functional blocks. Taken together, our results reveal that human enhancers contain a high density of sequence features required for their normal in vivo function and provide a rich resource for further exploration of human enhancer logic.

Does a Vertebrate Morphotype of Pallial Subdivisions Really Exist?

BACKGROUND

Comparative neuroanatomists have long sought to determine which part of the pallium in nonmammals is homologous to the mammalian neocortex. A number of similar connectivity patterns across species have led to the idea that the basic organization of the vertebrate brain is relatively conserved; thus, efforts of the last decades have been focused on determining a vertebrate "morphotype" - a model comprising the characteristics believed to have been present in the last common ancestor of all vertebrates.

SUMMARY

The endeavor to determine the vertebrate morphotype has been riddled with controversies due to the extensive morphological diversity of the pallium among vertebrate taxa. Nonetheless, most proposed scenarios of pallial homology are variants of a common theme where the vertebrate pallium is subdivided into subdivisions homologous to the hippocampus, neocortex, piriform cortex, and amygdala, in a one-to-one manner. We review the rationales of major propositions of pallial homology and identify the source of the discrepancies behind different hypotheses. We consider that a source of discrepancies is the prevailing assumption that there is a single "morphotype of the pallial subdivisions" throughout vertebrates. Instead, pallial subdivisions present in different taxa probably evolved independently in each lineage.

KEY MESSAGES

We encounter discrepancies when we search for a single morphotype of subdivisions across vertebrates. These discrepancies can be resolved by considering that several subdivisions within the pallium were established after the divergence of the different lineages. The differences of pallial organization are especially remarkable between actinopterygians (including teleost fishes) and other vertebrates. Thus, the prevailing notion of a simple one-to-one homology between the mammalian and teleost pallia needs to be reconsidered.

Five autism-associated transcriptional regulators target shared loci proximal to brain-expressed genes

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.

DLX genes and proteins in mammalian forebrain development

The vertebrate Dlx gene family encode homeobox transcription factors that are related to the Drosophila Distal-less (Dll) gene and are crucial for development. Over the last ∼35 years detailed information has accrued about the redundant and unique expression and function of the six mammalian Dlx family genes. DLX proteins interact with general transcriptional regulators, and co-bind with other transcription factors to enhancer elements with highly specific activity in the developing forebrain. Integration of the genetic and biochemical data has yielded a foundation for a gene regulatory network governing the differentiation of forebrain GABAergic neurons. In this Primer, we describe the discovery of vertebrate Dlx genes and their crucial roles in embryonic development. We largely focus on the role of Dlx family genes in mammalian forebrain development revealed through studies in mice. Finally, we highlight questions that remain unanswered regarding vertebrate Dlx genes despite over 30 years of research.

Using a comprehensive atlas and predictive models to reveal the complexity and evolution of brain-active regulatory elements

Most genetic variants associated with psychiatric disorders are located in noncoding regions of the genome. To investigate their functional implications, we integrate epigenetic data from the PsychENCODE Consortium and other published sources to construct a comprehensive atlas of candidate brain cis-regulatory elements. Using deep learning, we model these elements' sequence syntax and predict how binding sites for lineage-specific transcription factors contribute to cell type-specific gene regulation in various types of glia and neurons. The elements' evolutionary history suggests that new regulatory information in the brain emerges primarily via smaller sequence mutations within conserved mammalian elements rather than entirely new human- or primate-specific sequences. However, primate-specific candidate elements, particularly those active during fetal brain development and in excitatory neurons and astrocytes, are implicated in the heritability of brain-related human traits. Additionally, we introduce PsychSCREEN, a web-based platform offering interactive visualization of PsychENCODE-generated genetic and epigenetic data from diverse brain cell types in individuals with psychiatric disorders and healthy controls.

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.

Non-CG DNA methylation and MeCP2 stabilize repeated tuning of long genes that distinguish closely related neuron types

The extraordinary diversity of neuron types in the mammalian brain is delineated at the highest resolution by subtle gene expression differences that may require specialized molecular mechanisms to be maintained. Neurons uniquely express the longest genes in the genome and utilize neuron-enriched non-CG DNA methylation (mCA) together with the Rett syndrome protein, MeCP2, to control gene expression, but the function of these unique gene structures and machinery in regulating finely resolved neuron type-specific gene programs has not been explored. Here, we employ epigenomic and spatial transcriptomic analyses to discover a major role for mCA and MeCP2 in maintaining neuron type-specific gene programs at the finest scale of cellular resolution. We uncover differential susceptibility to MeCP2 loss in neuronal populations depending on global mCA levels and dissect methylation patterns and intragenic enhancer repression that drive overlapping and distinct gene regulation between neuron types. Strikingly, we show that mCA and MeCP2 regulate genes that are repeatedly tuned to differentiate neuron types at the highest cellular resolution, including spatially resolved, vision-dependent gene programs in the visual cortex. These repeatedly tuned genes display genomic characteristics, including long length, numerous intragenic enhancers, and enrichment for mCA, that predispose them to regulation by MeCP2. Thus, long gene regulation by the MeCP2 pathway maintains differential gene expression between closely-related neurons to facilitate the exceptional cellular diversity in the complex mammalian brain.

EnhancerTracker: Comparing cell-type-specific enhancer activity of DNA sequence triplets via an ensemble of deep convolutional neural networks

Motivation

Transcriptional enhancers - unlike promoters - are unrestrained by distance or strand orientation with respect to their target genes, making their computational identification a challenge. Further, there are insufficient numbers of confirmed enhancers for many cell types, preventing robust training of machine-learning-based models for enhancer prediction for such cell types.

Results

We present EnhancerTracker , a novel tool that leverages an ensemble of deep separable convolutional neural networks to identify cell-type-specific enhancers with the need of only two confirmed enhancers. EnhancerTracker is trained, validated, and tested on 52,789 putative enhancers obtained from the FANTOM5 Project and control sequences derived from the human genome. Unlike available tools, which accept one sequence at a time, the input to our tool is three sequences; the first two are enhancers active in the same cell type. EnhancerTracker outputs 1 if the third sequence is an enhancer active in the same cell type(s) where the first two enhancers are active. It outputs 0 otherwise. On a held-out set (15%), EnhancerTracker achieved an accuracy of 64%, a specificity of 93%, a recall of 35%, a precision of 84%, and an F1 score of 49%.

Availability and implementation

https://github.com/BioinformaticsToolsmith/EnhancerTracker.

Contact

hani.girgis@tamuk.edu.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

Joint epigenome profiling reveals cell-type-specific gene regulatory programmes in human cortical organoids

Gene expression is regulated by multiple epigenetic mechanisms, which are coordinated in development and disease. However, current multiomics methods are frequently limited to one or two modalities at a time, making it challenging to obtain a comprehensive gene regulatory signature. Here, we describe a method-3D genome, RNA, accessibility and methylation sequencing (3DRAM-seq)-that simultaneously interrogates spatial genome organization, chromatin accessibility and DNA methylation genome-wide and at high resolution. We combine 3DRAM-seq with immunoFACS and RNA sequencing in cortical organoids to map the cell-type-specific regulatory landscape of human neural development across multiple epigenetic layers. Finally, we apply a massively parallel reporter assay to profile cell-type-specific enhancer activity in organoids and to functionally assess the role of key transcription factors for human enhancer activation and function. More broadly, 3DRAM-seq can be used to profile the multimodal epigenetic landscape in rare cell types and different tissues.

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 ~7,000 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.

Combinatorial transcription factor binding encodes cis-regulatory wiring of forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to genomic cis-regulatory elements (cREs), orchestrating transcription programs. While studies of chromatin state and chromosomal interactions have revealed dynamic neurodevelopmental cRE landscapes, parallel understanding of the underlying TF binding lags. To elucidate the combinatorial TF-cRE interactions driving mouse basal ganglia development, we integrated ChIP-seq for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and transcriptional state, and transgenic enhancer assays. We identified TF-cREs modules with distinct chromatin features and enhancer activity that have complementary roles driving GABAergic neurogenesis and suppressing other developmental fates. While the majority of distal cREs were bound by one or two TFs, a small proportion were extensively bound, and these enhancers also exhibited exceptional evolutionary conservation, motif density, and complex chromosomal interactions. Our results provide new insights into how modules of combinatorial TF-cRE interactions activate and repress developmental expression programs and demonstrate the value of TF binding data in modeling gene regulatory wiring.

Discovery and Validation of Context-Dependent Synthetic Mammalian Promoters

Cellular transcription enables cells to adapt to various stimuli and maintain homeostasis. Transcription factors bind to transcription response elements (TREs) in gene promoters, initiating transcription. Synthetic promoters, derived from natural TREs, can be engineered to control exogenous gene expression using endogenous transcription machinery. This technology has found extensive use in biological research for applications including reporter gene assays, biomarker development, and programming synthetic circuits in living cells. However, a reliable and precise method for selecting minimally-sized synthetic promoters with desired background, amplitude, and stimulation response profiles has been elusive. In this study, we introduce a massively parallel reporter assay library containing 6184 synthetic promoters, each less than 250 bp in length. This comprehensive library allows for rapid identification of promoters with optimal transcriptional output parameters across multiple cell lines and stimuli. We showcase this library's utility to identify promoters activated in unique cell types, and in response to metabolites, mitogens, cellular toxins, and agonism of both aminergic and non-aminergic GPCRs. We further show these promoters can be used in luciferase reporter assays, eliciting 50-100 fold dynamic ranges in response to stimuli. Our platform is effective, easily implemented, and provides a solution for selecting short-length promoters with precise performance for a multitude of applications.

An unbiased AAV-STARR-seq screen revealing the enhancer activity map of genomic regions in the mouse brain in vivo

Enhancers are important cis-regulatory elements controlling cell-type specific expression patterns of genes. Furthermore, combinations of enhancers and minimal promoters are utilized to construct small, artificial promoters for gene delivery vectors. Large-scale functional screening methodology to construct genomic maps of enhancer activities has been successfully established in cultured cell lines, however, not yet applied to terminally differentiated cells and tissues in a living animal. Here, we transposed the Self-Transcribing Active Regulatory Region Sequencing (STARR-seq) technique to the mouse brain using adeno-associated-viruses (AAV) for the delivery of a highly complex screening library tiling entire genomic regions and covering in total 3 Mb of the mouse genome. We identified 483 sequences with enhancer activity, including sequences that were not predicted by DNA accessibility or histone marks. Characterizing the expression patterns of fluorescent reporters controlled by nine candidate sequences, we observed differential expression patterns also in sparse cell types. Together, our study provides an entry point for the unbiased study of enhancer activities in organisms during health and disease.

Advances in AAV technology for delivering genetically encoded cargo to the nonhuman primate nervous system

Modern neuroscience approaches including optogenetics, calcium imaging, and other genetic manipulations have facilitated our ability to dissect specific circuits in rodent models to study their role in neurological disease. These approaches regularly use viral vectors to deliver genetic cargo (e.g., opsins) to specific tissues and genetically-engineered rodents to achieve cell-type specificity. However, the translatability of these rodent models, cross-species validation of identified targets, and translational efficacy of potential therapeutics in larger animal models like nonhuman primates remains difficult due to the lack of efficient primate viral vectors. A refined understanding of the nonhuman primate nervous system promises to deliver insights that can guide the development of treatments for neurological and neurodegenerative diseases. Here, we outline recent advances in the development of adeno-associated viral vectors for optimized use in nonhuman primates. These tools promise to help open new avenues for study in translational neuroscience and further our understanding of the primate brain.

Evolutionary relevance of single nucleotide variants within the forebrain exclusive human accelerated enhancer regions

Background

Human accelerated regions (HARs) are short conserved genomic sequences that have acquired significantly more nucleotide substitutions than expected in the human lineage after divergence from chimpanzees. The fast evolution of HARs may reflect their roles in the origin of human-specific traits. A recent study has reported positively-selected single nucleotide variants (SNVs) within brain-exclusive human accelerated enhancers (BE-HAEs) hs1210 (forebrain), hs563 (hindbrain) and hs304 (midbrain/forebrain). By including data from archaic hominins, these SNVs were shown to be Homo sapiens-specific, residing within transcriptional factors binding sites (TFBSs) for SOX2 (hs1210), RUNX1/3 (hs563), and FOS/JUND (hs304). Although these findings suggest that the predicted modifications in TFBSs may have some role in present-day brain structure, work is required to verify the extent to which these changes translate into functional variation.

Results

To start to fill this gap, we investigate the SOX2 SNV, with both forebrain expression and strong signal of positive selection in humans. We demonstrate that the HMG box of SOX2 binds in vitro with Homo sapiens-specific derived A-allele and ancestral T-allele carrying DNA sites in BE-HAE hs1210. Molecular docking and simulation analysis indicated highly favourable binding of HMG box with derived A-allele containing DNA site when compared to site carrying ancestral T-allele.

Conclusion

These results suggest that adoptive changes in TF affinity within BE-HAE hs1210 and other HAR enhancers in the evolutionary history of Homo sapiens might.

have brought about changes in gene expression patterns and have functional consequences on forebrain formation and evolution.

Methods

The present study employ electrophoretic mobility shift assays (EMSA) and molecular docking and molecular dynamics simulations approaches.

K27M in canonical and noncanonical H3 variants occurs in distinct oligodendroglial cell lineages in brain midline gliomas

Canonical (H3.1/H3.2) and noncanonical (H3.3) histone 3 K27M-mutant gliomas have unique spatiotemporal distributions, partner alterations and molecular profiles. The contribution of the cell of origin to these differences has been challenging to uncouple from the oncogenic reprogramming induced by the mutation. Here, we perform an integrated analysis of 116 tumors, including single-cell transcriptome and chromatin accessibility, 3D chromatin architecture and epigenomic profiles, and show that K27M-mutant gliomas faithfully maintain chromatin configuration at developmental genes consistent with anatomically distinct oligodendrocyte precursor cells (OPCs). H3.3K27M thalamic gliomas map to prosomere 2-derived lineages. In turn, H3.1K27M ACVR1-mutant pontine gliomas uniformly mirror early ventral NKX6-1⁺/SHH-dependent brainstem OPCs, whereas H3.3K27M gliomas frequently resemble dorsal PAX3⁺/BMP-dependent progenitors. Our data suggest a context-specific vulnerability in H3.1K27M-mutant SHH-dependent ventral OPCs, which rely on acquisition of ACVR1 mutations to drive aberrant BMP signaling required for oncogenesis. The unifying action of K27M mutations is to restrict H3K27me3 at PRC2 landing sites, whereas other epigenetic changes are mainly contingent on the cell of origin chromatin state and cycling rate.

A multi-omics approach to visualize early neuronal differentiation from hESCs in 4D

Neuronal differentiation of pluripotent stem cells is an established method to study physiology, disease, and medication safety. However, the sequence of events in human neuronal differentiation and the ability of in vitro models to recapitulate early brain development are poorly understood. We developed a protocol optimized for the study of early human brain development and neuropharmacological applications. We comprehensively characterized gene expression and epigenetic profiles at four timepoints, because the cells differentiate from embryonic stem cells towards a heterogeneous population of progenitors, immature and mature neurons bearing telencephalic signatures. A multi-omics roadmap of neuronal differentiation, combined with searchable interactive gene analysis tools, allows for extensive exploration of early neuronal development and the effect of medications.

Spatially resolved epigenomic profiling of single cells in complex tissues

The recent development of spatial omics methods has enabled single-cell profiling of the transcriptome and 3D genome organization with high spatial resolution. Expanding the repertoire of spatial omics tools, a spatially resolved single-cell epigenomics method will accelerate understanding of the spatial regulation of cell and tissue functions. Here, we report a method for spatially resolved epigenomic profiling of single cells using in situ tagmentation and transcription followed by multiplexed imaging. We demonstrated the ability to profile histone modifications marking active promoters, putative enhancers, and silent promoters in individual cells, and generated high-resolution spatial atlas of hundreds of active promoters and putative enhancers in embryonic and adult mouse brains. Our results suggested putative promoter-enhancer pairs and enhancer hubs regulating developmentally important genes. We envision this approach will be generally applicable to spatial profiling of epigenetic modifications and DNA-binding proteins, advancing our understanding of how gene expression is spatiotemporally regulated by the epigenome.

A map of cis-regulatory modules and constituent transcription factor binding sites in 80% of the mouse genome

BACKGROUND

Mouse is probably the most important model organism to study mammal biology and human diseases. A better understanding of the mouse genome will help understand the human genome, biology and diseases. However, despite the recent progress, the characterization of the regulatory sequences in the mouse genome is still far from complete, limiting its use to understand the regulatory sequences in the human genome.

RESULTS

Here, by integrating binding peaks in ~ 9,000 transcription factor (TF) ChIP-seq datasets that cover 79.9% of the mouse mappable genome using an efficient pipeline, we were able to partition these binding peak-covered genome regions into a cis-regulatory module (CRM) candidate (CRMC) set and a non-CRMC set. The CRMCs contain 912,197 putative CRMs and 38,554,729 TF binding sites (TFBSs) islands, covering 55.5% and 24.4% of the mappable genome, respectively. The CRMCs tend to be under strong evolutionary constraints, indicating that they are likely cis-regulatory; while the non-CRMCs are largely selectively neutral, indicating that they are unlikely cis-regulatory. Based on evolutionary profiles of the genome positions, we further estimated that 63.8% and 27.4% of the mouse genome might code for CRMs and TFBSs, respectively.

CONCLUSIONS

Validation using experimental data suggests that at least most of the CRMCs are authentic. Thus, this unprecedentedly comprehensive map of CRMs and TFBSs can be a good resource to guide experimental studies of regulatory genomes in mice and humans.

Accurate prediction of functional states of cis-regulatory modules reveals common epigenetic rules in humans and mice

BACKGROUND

Predicting cis-regulatory modules (CRMs) in a genome and their functional states in various cell/tissue types of the organism are two related challenging computational tasks. Most current methods attempt to simultaneously achieve both using data of multiple epigenetic marks in a cell/tissue type. Though conceptually attractive, they suffer high false discovery rates and limited applications. To fill the gaps, we proposed a two-step strategy to first predict a map of CRMs in the genome, and then predict functional states of all the CRMs in various cell/tissue types of the organism. We have recently developed an algorithm for the first step that was able to more accurately and completely predict CRMs in a genome than existing methods by integrating numerous transcription factor ChIP-seq datasets in the organism. Here, we presented machine-learning methods for the second step.

RESULTS

We showed that functional states in a cell/tissue type of all the CRMs in the genome could be accurately predicted using data of only 1~4 epigenetic marks by a variety of machine-learning classifiers. Our predictions are substantially more accurate than the best achieved so far. Interestingly, a model trained on a cell/tissue type in humans can accurately predict functional states of CRMs in different cell/tissue types of humans as well as of mice, and vice versa. Therefore, epigenetic code that defines functional states of CRMs in various cell/tissue types is universal at least in humans and mice. Moreover, we found that from tens to hundreds of thousands of CRMs were active in a human and mouse cell/tissue type, and up to 99.98% of them were reutilized in different cell/tissue types, while as small as 0.02% of them were unique to a cell/tissue type that might define the cell/tissue type.

CONCLUSIONS

Our two-step approach can accurately predict functional states in any cell/tissue type of all the CRMs in the genome using data of only 1~4 epigenetic marks. Our approach is also more cost-effective than existing methods that typically use data of more epigenetic marks. Our results suggest common epigenetic rules for defining functional states of CRMs in various cell/tissue types in humans and mice.

A Ctnnb1 enhancer regulates neocortical neurogenesis by controlling the abundance of intermediate progenitors

β-catenin-dependent canonical Wnt signaling plays a plethora of roles in neocortex (Ncx) development, but its function in regulating the abundance of intermediate progenitors (IPs) is elusive. Here we identified neCtnnb1, an evolutionarily conserved cis-regulatory element with typical enhancer features in developing Ncx. neCtnnb1 locates 55 kilobase upstream of and spatially close to the promoter of Ctnnb1, the gene encoding β-catenin. CRISPR/Cas9-mediated activation or interference of the neCtnnb1 locus enhanced or inhibited transcription of Ctnnb1. neCtnnb1 drove transcription predominantly in the subventricular zone of developing Ncx. Knock-out of neCtnnb1 in mice resulted in compromised expression of Ctnnb1 and the Wnt reporter in developing Ncx. Importantly, knock-out of neCtnnb1 lead to reduced production and transit-amplification of IPs, which subsequently generated fewer upper-layer Ncx projection neurons (PNs). In contrast, enhancing the canonical Wnt signaling by stabilizing β-catenin in neCtnnb1-active cells promoted the production of IPs and upper-layer Ncx PNs. ASH2L was identified as the key trans-acting factor that associates with neCtnnb1 and Ctnnb1’s promoter to maintain Ctnnb1’s transcription in both mouse and human Ncx progenitors. These findings advance understanding of transcriptional regulation of Ctnnb1, and provide insights into mechanisms underlying Ncx expansion during development.

Analysis of somatic mutations in 131 human brains reveals aging-associated hypermutability

We analyzed 131 human brains (44 neurotypical, 19 with Tourette syndrome, 9 with schizophrenia, and 59 with autism) for somatic mutations after whole genome sequencing to a depth of more than 200×. Typically, brains had 20 to 60 detectable single-nucleotide mutations, but ~6% of brains harbored hundreds of somatic mutations. Hypermutability was associated with age and damaging mutations in genes implicated in cancers and, in some brains, reflected in vivo clonal expansions. Somatic duplications, likely arising during development, were found in ~5% of normal and diseased brains, reflecting background mutagenesis. Brains with autism were associated with mutations creating putative transcription factor binding motifs in enhancer-like regions in the developing brain. The top-ranked affected motifs corresponded to MEIS (myeloid ectopic viral integration site) transcription factors, suggesting a potential link between their involvement in gene regulation and autism.

PCRMS: a database of predicted cis-regulatory modules and constituent transcription factor binding sites in genomes

More accurate and more complete predictions of cis-regulatory modules (CRMs) and constituent transcription factor (TF) binding sites (TFBSs) in genomes can facilitate characterizing functions of regulatory sequences. Here, we developed a database predicted cis-regulatory modules (PCRMS) (https://cci-bioinfo.uncc.edu) that stores highly accurate and unprecedentedly complete maps of predicted CRMs and TFBSs in the human and mouse genomes. The web interface allows the user to browse CRMs and TFBSs in an organism, find the closest CRMs to a gene, search CRMs around a gene and find all TFBSs of a TF. PCRMS can be a useful resource for the research community to characterize regulatory genomes. Database URL: https://cci-bioinfo.uncc.edu/.

Single cell enhancer activity distinguishes GABAergic and cholinergic lineages in embryonic mouse basal ganglia

During brain development, neurons are generated by spatially and temporally distinct processes that remain to be fully characterized. The ganglionic eminences (GEs) in the embryonic subpallium give rise to GABAergic and cholinergic neuron lineages that form the basal ganglia or migrate to the cerebral cortex. Beyond a limited set of canonical RNA markers, the transcriptional states of GE progenitors and immature neurons cells remain poorly defined. We combine enhancer labeling, single-cell transcriptomics using transcription factor-anchored clustering, and integration with in situ hybridization data to distinguish emerging neuronal populations in embryonic mouse basal ganglia. Our results demonstrate the specificity of enhancer-based labeling at single-cell resolution and reveal developmental origins and specification processes of critical neuronal lineages.

Enhancers integrate transcription factor signaling pathways that drive cell fate specification in the developing brain. We paired enhancer labeling and single-cell RNA-sequencing (scRNA-seq) to delineate and distinguish specification of neuronal lineages in mouse medial, lateral, and caudal ganglionic eminences (MGE, LGE, and CGE) at embryonic day (E)11.5. We show that scRNA-seq clustering using transcription factors improves resolution of regional and developmental populations, and that enhancer activities identify specific and overlapping GE-derived neuronal populations. First, we mapped the activities of seven evolutionarily conserved brain enhancers at single-cell resolution in vivo, finding that the selected enhancers had diverse activities in specific progenitor and neuronal populations across the GEs. We then applied enhancer-based labeling, scRNA-seq, and analysis of in situ hybridization data to distinguish transcriptionally distinct and spatially defined subtypes of MGE-derived GABAergic and cholinergic projection neurons and interneurons. Our results map developmental origins and specification paths underlying neurogenesis in the embryonic basal ganglia and showcase the power of scRNA-seq combined with enhancer-based labeling to resolve the complex paths of neuronal specification underlying mouse brain development.

AAV Deployment of Enhancer-Based Expression Constructs In Vivo in Mouse Brain

Enhancers are binding platforms for a diverse array of transcription factors that drive specific expression patterns of tissue- and cell-type-specific genes. Multiple means of assessing non-coding DNA and various chromatin states have proven useful in predicting the presence of enhancer sequences in the genome, but validating the activity of these sequences and finding the organs and developmental stages they are active in is a labor-intensive process. Recent advances in adeno-associated virus (AAV) vectors have enabled the widespread delivery of transgenes to mouse tissues, enabling in vivo enhancer testing without necessitating a transgenic animal. This protocol shows how a reporter construct that expresses EGFP under the control of a minimal promoter, which does not drive significant expression on its own, can be used to study the activity patterns of candidate enhancer sequences in the mouse brain. An AAV-packaged reporter construct is delivered to the mouse brain and incubated for 1-4 weeks, after which the animal is sacrificed, and brain sections are observed under a microscope. EGFP appears in cells in which the tested enhancer is sufficient to initiate gene expression, pinpointing the location and developmental stage in which the enhancer is active in the brain. Standard cloning methods, low-cost AAV packaging, and expanding AAV serotypes and methods for in vivo delivery and standard imaging readout make this an accessible approach for the study of how gene expression is regulated in the brain.

Further Delineation of Duplications of ARX Locus Detected in Male Patients with Varying Degrees of Intellectual Disability

The X-linked gene encoding aristaless-related homeobox (ARX) is a bi-functional transcription factor capable of activating or repressing gene transcription, whose mutations have been found in a wide spectrum of neurodevelopmental disorders (NDDs); these include cortical malformations, paediatric epilepsy, intellectual disability (ID) and autism. In addition to point mutations, duplications of the ARX locus have been detected in male patients with ID. These rearrangements include telencephalon ultraconserved enhancers, whose structural alterations can interfere with the control of ARX expression in the developing brain. Here, we review the structural features of 15 gain copy-number variants (CNVs) of the ARX locus found in patients presenting wide-ranging phenotypic variations including ID, speech delay, hypotonia and psychiatric abnormalities. We also report on a further novel Xp21.3 duplication detected in a male patient with moderate ID and carrying a fully duplicated copy of the ARX locus and the ultraconserved enhancers. As consequences of this rearrangement, the patient-derived lymphoblastoid cell line shows abnormal activity of the ARX-KDM5C-SYN1 regulatory axis. Moreover, the three-dimensional (3D) structure of the Arx locus, both in mouse embryonic stem cells and cortical neurons, provides new insight for the functional consequences of ARX duplications. Finally, by comparing the clinical features of the 16 CNVs affecting the ARX locus, we conclude that—depending on the involvement of tissue-specific enhancers—the ARX duplications are ID-associated risk CNVs with variable expressivity and penetrance.

Perfect and imperfect views of ultraconserved sequences

Across the human genome, there are nearly 500 'ultraconserved' elements: regions of at least 200 contiguous nucleotides that are perfectly conserved in both the mouse and rat genomes. Remarkably, the majority of these sequences are non-coding, and many can function as enhancers that activate tissue-specific gene expression during embryonic development. From their first description more than 15 years ago, their extreme conservation has both fascinated and perplexed researchers in genomics and evolutionary biology. The intrigue around ultraconserved elements only grew with the observation that they are dispensable for viability. Here, we review recent progress towards understanding the general importance and the specific functions of ultraconserved sequences in mammalian development and human disease and discuss possible explanations for their extreme conservation.

Dlx1/2-dependent expression of Meis2 promotes neuronal fate determination in the mammalian striatum

The striatum is a central regulator of behavior and motor function through the actions of D1 and D2 medium-sized spiny neurons (MSNs), which arise from a common lateral ganglionic eminence (LGE) progenitor. The molecular mechanisms of cell fate specification of these two neuronal subtypes are incompletely understood. Here, we found that deletion of murine Meis2, which is highly expressed in the LGE and derivatives, led to a large reduction in striatal MSNs due to a block in their differentiation. Meis2 directly binds to the Zfp503 and Six3 promoters and is required for their expression and specification of D1 and D2 MSNs, respectively. Finally, Meis2 expression is regulated by Dlx1/2 at least partially through the enhancer hs599 in the LGE subventricular zone. Overall, our findings define a pathway in the LGE whereby Dlx1/2 drives expression of Meis2, which subsequently promotes the fate determination of striatal D1 and D2 MSNs via Zfp503 and Six3.

Multimodal profiling of the transcriptional regulatory landscape of the developing mouse cortex identifies Neurog2 as a key epigenome remodeler

How multiple epigenetic layers and transcription factors (TFs) interact to facilitate brain development is largely unknown. Here, to systematically map the regulatory landscape of neural differentiation in the mouse neocortex, we profiled gene expression and chromatin accessibility in single cells and integrated these data with measurements of enhancer activity, DNA methylation and three-dimensional genome architecture in purified cell populations. This allowed us to identify thousands of new enhancers, their predicted target genes and the temporal relationships between enhancer activation, epigenome remodeling and gene expression. We characterize specific neuronal transcription factors associated with extensive and frequently coordinated changes across multiple epigenetic modalities. In addition, we functionally demonstrate a new role for Neurog2 in directly mediating enhancer activity, DNA demethylation, increasing chromatin accessibility and facilitating chromatin looping in vivo. Our work provides a global view of the gene regulatory logic of lineage specification in the cerebral cortex.

By profiling multiple epigenetic layers and enhancer activity in vivo, the authors show a widespread remodeling of the regulatory landscape during mouse cortical development and identify Neurog2 as a key transcription factor driving this process.

Common Genetic Variants Contribute to Risk of Transposition of the Great Arteries

RATIONALE

Dextro-transposition of the great arteries (D-TGA) is a severe congenital heart defect which affects approximately 1 in 4,000 live births. While there are several reports of D-TGA patients with rare variants in individual genes, the majority of D-TGA cases remain genetically elusive. Familial recurrence patterns and the observation that most cases with D-TGA are sporadic suggest a polygenic inheritance for the disorder, yet this remains unexplored.

OBJECTIVE

We sought to study the role of common single nucleotide polymorphisms (SNPs) in risk for D-TGA.

METHODS AND RESULTS

We conducted a genome-wide association study in an international set of 1,237 patients with D-TGA and identified a genome-wide significant susceptibility locus on chromosome 3p14.3, which was subsequently replicated in an independent case-control set (rs56219800, meta-analysis P=8.6x10-10, OR=0.69 per C allele). SNP-based heritability analysis showed that 25% of variance in susceptibility to D-TGA may be explained by common variants. A genome-wide polygenic risk score derived from the discovery set was significantly associated to D-TGA in the replication set (P=4x10-5). The genome-wide significant locus (3p14.3) co-localizes with a putative regulatory element that interacts with the promoter of WNT5A, which encodes the Wnt Family Member 5A protein known for its role in cardiac development in mice. We show that this element drives reporter gene activity in the developing heart of mice and zebrafish and is bound by the developmental transcription factor TBX20. We further demonstrate that TBX20 attenuates Wnt5a expression levels in the developing mouse heart.

CONCLUSIONS

This work provides support for a polygenic architecture in D-TGA and identifies a susceptibility locus on chromosome 3p14.3 near WNT5A. Genomic and functional data support a causal role of WNT5A at the locus.

An expansion of the non-coding genome and its regulatory potential underlies vertebrate neuronal diversity

Proper assembly and function of the nervous system requires the generation of a uniquely diverse population of neurons expressing a cell-type-specific combination of effector genes that collectively define neuronal morphology, connectivity, and function. How countless partially overlapping but cell-type-specific patterns of gene expression are controlled at the genomic level remains poorly understood. Here we show that neuronal genes are associated with highly complex gene regulatory systems composed of independent cell-type- and cell-stage-specific regulatory elements that reside in expanded non-coding genomic domains. Mapping enhancer-promoter interactions revealed that motor neuron enhancers are broadly distributed across the large chromatin domains. This distributed regulatory architecture is not a unique property of motor neurons but is employed throughout the nervous system. The number of regulatory elements increased dramatically during the transition from invertebrates to vertebrates, suggesting that acquisition of new enhancers might be a fundamental process underlying the evolutionary increase in cellular complexity.

Characterization of Mammalian In Vivo Enhancers Using Mouse Transgenesis and CRISPR Genome Editing

Embryonic morphogenesis is strictly dependent on tight spatiotemporal control of developmental gene expression, which is typically achieved through the concerted activity of multiple enhancers driving cell type-specific expression of a target gene. Mammalian genomes are organized in topologically associated domains, providing a preferred environment and framework for interactions between transcriptional enhancers and gene promoters. While epigenomic profiling and three-dimensional chromatin conformation capture have significantly increased the accuracy of identifying enhancers, assessment of subregional enhancer activities via transgenic reporter assays in mice remains the gold standard for assigning enhancer activity in vivo. Once this activity is defined, the ideal method to explore the functional necessity of a transcriptional enhancer and its contribution to target gene dosage and morphological or physiological processes is deletion of the enhancer sequence from the mouse genome. Here we present detailed protocols for efficient introduction of enhancer-reporter transgenes and CRISPR-mediated genomic deletions into the mouse genome, including a step-by-step guide for pronuclear microinjection of fertilized mouse eggs. We provide instructions for the assembly and genomic integration of enhancer-reporter cassettes that have been used for validation of thousands of putative enhancer sequences accessible through the VISTA enhancer browser, including a recently published method for robust site-directed transgenesis at the H11 safe-harbor locus. Together, these methods enable rapid and large-scale assessment of enhancer activities and sequence variants in mice, which is essential to understand mammalian genome function and genetic diseases.

Transcriptional network orchestrating regional patterning of cortical progenitors

We uncovered a transcription factor (TF) network that regulates cortical regional patterning in radial glial stem cells. Screening the expression of hundreds of TFs in the developing mouse cortex identified 38 TFs that are expressed in gradients in the ventricular zone (VZ). We tested whether their cortical expression was altered in mutant mice with known patterning defects (Emx2, Nr2f1, and Pax6), which enabled us to define a cortical regionalization TF network (CRTFN). To identify genomic programming underlying this network, we performed TF ChIP-seq and chromatin-looping conformation to identify enhancer-gene interactions. To map enhancers involved in regional patterning of cortical progenitors, we performed assays for epigenomic marks and DNA accessibility in VZ cells purified from wild-type and patterning mutant mice. This integrated approach has identified a CRTFN and VZ enhancers involved in cortical regional patterning in the mouse.

Cortical Cartography: Mapping Arealization Using Single-Cell Omics Technology

The cerebral cortex derives its cognitive power from a modular network of specialized areas processing a multitude of information. The assembly and organization of these regions is vital for human behavior and perception, as evidenced by the prevalence of area-specific phenotypes that manifest in neurodevelopmental and psychiatric disorders. Generations of scientists have examined the architecture of the human cortex, but efforts to capture the gene networks which drive arealization have been hampered by the lack of tractable models of human neurodevelopment. Advancements in "omics" technologies, imaging, and computational power have enabled exciting breakthroughs into the molecular and structural characteristics of cortical areas, including transcriptomic, epigenomic, metabolomic, and proteomic profiles of mammalian models. Here we review the single-omics atlases that have shaped our current understanding of cortical areas, and their potential to fuel a new era of multi-omic single-cell endeavors to interrogate both the developing and adult human cortex.

Single-cell epigenomics reveals mechanisms of human cortical development

During mammalian development, differences in chromatin state coincide with cellular differentiation and reflect changes in the gene regulatory landscape1. In the developing brain, cell fate specification and topographic identity are important for defining cell identity2 and confer selective vulnerabilities to neurodevelopmental disorders3. Here, to identify cell-type-specific chromatin accessibility patterns in the developing human brain, we used a single-cell assay for transposase accessibility by sequencing (scATAC-seq) in primary tissue samples from the human forebrain. We applied unbiased analyses to identify genomic loci that undergo extensive cell-type- and brain-region-specific changes in accessibility during neurogenesis, and an integrative analysis to predict cell-type-specific candidate regulatory elements. We found that cerebral organoids recapitulate most putative cell-type-specific enhancer accessibility patterns but lack many cell-type-specific open chromatin regions that are found in vivo. Systematic comparison of chromatin accessibility across brain regions revealed unexpected diversity among neural progenitor cells in the cerebral cortex and implicated retinoic acid signalling in the specification of neuronal lineage identity in the prefrontal cortex. Together, our results reveal the important contribution of chromatin state to the emerging patterns of cell type diversity and cell fate specification and provide a blueprint for evaluating the fidelity and robustness of cerebral organoids as a model for cortical development.

A statistical framework for predicting critical regions of p53-dependent enhancers

P53 is the 'guardian of the genome' and is responsible for regulating cell cycle and apoptosis. The genomic p53 binding regions, where activating transcriptional factors and cofactors like p300 simultaneously bind, are called 'p53-dependent enhancers', which play an important role in tumorigenesis. Current experimental assays generally provide a broad peak of each enhancer element, leaving our knowledge about critical enhancer regions (CERs) limited. Under the inspiration of enhancer dissection by CRISPR-Cas9 screen library on genome-wide p53 binding sites, here we introduce a statistical framework called 'Computational CRISPR Strategy' (CCS), to predict whether a given DNA fragment will be a p53-dependent CER by employing 7-mer as feature extractions along with random forest as the regressor. When training on a p53 CRISPR enhancer dataset, CCS not only accurately fitted the top-ranked enriched single guide RNAs (sgRNAs) but also successfully reproduced two known CERs that were validated by experiments. When applying it to an independent testing dataset on a tilling of a 2K-b genomic region of CRISPR-deCDKN1A-Lib, the trained model shows great generalizability by identifying a CER containing five top-ranked sgRNAs. A feature importance analysis further indicates that top-ranked 7-mers are mapped onto informative TF motifs including POU5F1 and SOX5, which are differentially enriched in p53-dependent CERs and are potential factors to make a general p53 binding site to form a p53-dependent CER, providing the interpretability of the trained model. Our results demonstrate that CCS is an alternative way of the CRISPR experiment to screen the genome for mapping p53-dependent CERs.

The impact of (ab)normal maternal environment on cortical development

The cortex in the mammalian brain is the most complex brain region that integrates sensory information and coordinates motor and cognitive processes. To perform such functions, the cortex contains multiple subtypes of neurons that are generated during embryogenesis. Newly born neurons migrate to their proper location in the cortex, grow axons and dendrites, and form neuronal circuits. These developmental processes in the fetal brain are regulated to a large extent by a great variety of factors derived from the mother - starting from simple nutrients as building blocks and ending with hormones. Thus, when the normal maternal environment is disturbed due to maternal infection, stress, malnutrition, or toxic substances, it might have a profound impact on cortical development and the offspring can develop a variety of neurodevelopmental disorders. Here we first describe the major developmental processes which generate neuronal diversity in the cortex. We then review our knowledge of how most common maternal insults affect cortical development, perturb neuronal circuits, and lead to neurodevelopmental disorders. We further present a concept of selective vulnerability of cortical neuronal subtypes to maternal-derived insults, where the vulnerability of cortical neurons and their progenitors to an insult depends on the time (developmental period), place (location in the developing brain), and type (unique features of a cell type and an insult). Finally, we provide evidence for the existence of selective vulnerability during cortical development and identify the most vulnerable neuronal types, stages of differentiation, and developmental time for major maternal-derived insults.

Functional and structural basis of extreme conservation in vertebrate 5' untranslated regions

The lack of knowledge about extreme conservation in genomes remains a major gap in our understanding of the evolution of gene regulation. Here, we reveal an unexpected role of extremely conserved 5' untranslated regions (UTRs) in noncanonical translational regulation that is linked to the emergence of essential developmental features in vertebrate species. Endogenous deletion of conserved elements within these 5' UTRs decreased gene expression, and extremely conserved 5' UTRs possess cis-regulatory elements that promote cell-type-specific regulation of translation. We further developed in-cell mutate-and-map (icM²), a new methodology that maps RNA structure inside cells. Using icM², we determined that an extremely conserved 5' UTR encodes multiple alternative structures and that each single nucleotide within the conserved element maintains the balance of alternative structures important to control the dynamic range of protein expression. These results explain how extreme sequence conservation can lead to RNA-level biological functions encoded in the untranslated regions of vertebrate genomes.

Functional enhancer elements drive subclass-selective expression from mouse to primate neocortex

Viral genetic tools that target specific brain cell types could transform basic neuroscience and targeted gene therapy. Here, we use comparative open chromatin analysis to identify thousands of human-neocortical-subclass-specific putative enhancers from across the genome to control gene expression in adeno-associated virus (AAV) vectors. The cellular specificity of reporter expression from enhancer-AAVs is established by molecular profiling after systemic AAV delivery in mouse. Over 30% of enhancer-AAVs produce specific expression in the targeted subclass, including both excitatory and inhibitory subclasses. We present a collection of Parvalbumin (PVALB) enhancer-AAVs that show highly enriched expression not only in cortical PVALB cells but also in some subcortical PVALB populations. Five vectors maintain PVALB-enriched expression in primate neocortex. These results demonstrate how genome-wide open chromatin data mining and cross-species AAV validation can be used to create the next generation of non-species-restricted viral genetic tools.

A flexible repertoire of transcription factor binding sites and a diversity threshold determines enhancer activity in embryonic stem cells

Transcriptional enhancers are critical for development and phenotype evolution and are often mutated in disease contexts; however, even in well-studied cell types, the sequence code conferring enhancer activity remains unknown. To examine the enhancer regulatory code for pluripotent stem cells, we identified genomic regions with conserved binding of multiple transcription factors in mouse and human embryonic stem cells (ESCs). Examination of these regions revealed that they contain on average 12.6 conserved transcription factor binding site (TFBS) sequences. Enriched TFBSs are a diverse repertoire of 70 different sequences representing the binding sequences of both known and novel ESC regulators. Using a diverse set of TFBSs from this repertoire was sufficient to construct short synthetic enhancers with activity comparable to native enhancers. Site-directed mutagenesis of conserved TFBSs in endogenous enhancers or TFBS deletion from synthetic sequences revealed a requirement for 10 or more different TFBSs. Furthermore, specific TFBSs, including the POU5F1:SOX2 comotif, are dispensable, despite cobinding the POU5F1 (also known as OCT4), SOX2, and NANOG master regulators of pluripotency. These findings reveal that a TFBS sequence diversity threshold overrides the need for optimized regulatory grammar and individual TFBSs that recruit specific master regulators.

Heart Enhancers: Development and Disease Control at a Distance

Bound by lineage-determining transcription factors and signaling effectors, enhancers play essential roles in controlling spatiotemporal gene expression profiles during development, homeostasis and disease. Recent synergistic advances in functional genomic technologies, combined with the developmental biology toolbox, have resulted in unprecedented genome-wide annotation of heart enhancers and their target genes. Starting with early studies of vertebrate heart enhancers and ending with state-of-the-art genome-wide enhancer discovery and testing, we will review how studying heart enhancers in metazoan species has helped inform our understanding of cardiac development and disease.

Exploring the biological role of postzygotic and germinal de novo mutations in ASD

De novo mutations (DNMs), including germinal and postzygotic mutations (PZMs), are a strong source of causality for Autism Spectrum Disorder (ASD). However, the biological processes involved behind them remain unexplored. Our aim was to detect DNMs (germinal and PZMs) in a Spanish ASD cohort (360 trios) and to explore their role across different biological hierarchies (gene, biological pathway, cell and brain areas) using bioinformatic approaches. For the majority of the analysis, a combined ASD cohort (N = 2171 trios) was created using previously published data by the Autism Sequencing Consortium (ASC). New plausible candidate genes for ASD such as FMR1 and NFIA were found. In addition, genes harboring PZMs were significantly enriched for miR-137 targets in comparison with germinal DNMs that were enriched in GO terms related to synaptic transmission. The expression pattern of genes with PZMs was restricted to early mid-fetal cortex. In contrast, the analysis of genes with germinal DNMs revealed a spatio-temporal window from early to mid-fetal development stages, with expression in the amygdala, cerebellum, cortex and striatum. These results provide evidence of the pathogenic role of PZMs and suggest the existence of distinct mechanisms between PZMs and germinal DNMs that are influencing ASD risk.

Parallel functional testing identifies enhancers active in early postnatal mouse brain

Enhancers are cis-regulatory elements that play critical regulatory roles in modulating developmental transcription programs and driving cell-type-specific and context-dependent gene expression in the brain. The development of massively parallel reporter assays (MPRAs) has enabled high-throughput functional screening of candidate DNA sequences for enhancer activity. Tissue-specific screening of in vivo enhancer function at scale has the potential to greatly expand our understanding of the role of non-coding sequences in development, evolution, and disease. Here, we adapted a self-transcribing regulatory element MPRA strategy for delivery to early postnatal mouse brain via recombinant adeno-associated virus (rAAV). We identified and validated putative enhancers capable of driving reporter gene expression in mouse forebrain, including regulatory elements within an intronic CACNA1C linkage disequilibrium block associated with risk in neuropsychiatric disorder genetic studies. Paired screening and single enhancer in vivo functional testing, as we show here, represents a powerful approach towards characterizing regulatory activity of enhancers and understanding how enhancer sequences organize gene expression in the brain.