Layer-specific chromatin accessibility landscapes reveal regulatory networks in adult mouse visual cortex

A suite of enhancer AAVs and transgenic mouse lines for genetic access to cortical cell types

The mammalian cortex is comprised of cells with different morphological, physiological, and molecular properties that can be classified according to shared properties into cell types. Defining the contribution of each cell type to the computational and cognitive processes that are guided by the cortex is essential for understanding its function in health and disease. We use transcriptomic and epigenomic cortical cell type taxonomies from mice and humans to define marker genes and enhancers, and to build genetic tools for cortical cell types. Here, we present a large toolkit for selective targeting of cortical populations, including mouse transgenic lines and recombinant adeno-associated virus (AAV) vectors containing genomic enhancers. We report evaluation of fifteen new transgenic driver lines and over 680 different enhancer AAVs covering all major subclasses of cortical cells, with many achieving a high degree of specificity, comparable with existing transgenic lines. We find that the transgenic lines based on marker genes can provide exceptional specificity and completeness of cell type labeling, but frequently require generation of a triple-transgenic cross for best usability/specificity. On the other hand, enhancer AAVs are easy to screen and use, and can be easily modified to express diverse cargo, such as recombinases. However, their use depends on many factors, such as viral titer and route of administration. The tools reported here as well as the scaled process of tool creation provide an unprecedented resource that should enable diverse experimental strategies towards understanding mammalian cortex and brain function.

Molecular mechanisms reconstruction from single-cell multi-omics data with HuMMuS

MOTIVATION

The molecular identity of a cell results from a complex interplay between heterogeneous molecular layers. Recent advances in single-cell sequencing technologies have opened the possibility to measure such molecular layers of regulation.

RESULTS

Here, we present HuMMuS, a new method for inferring regulatory mechanisms from single-cell multi-omics data. Differently from the state-of-the-art, HuMMuS captures cooperation between biological macromolecules and can easily include additional layers of molecular regulation. We benchmarked HuMMuS with respect to the state-of-the-art on both paired and unpaired multi-omics datasets. Our results proved the improvements provided by HuMMuS in terms of transcription factor (TF) targets, TF binding motifs and regulatory regions prediction. Finally, once applied to snmC-seq, scATAC-seq and scRNA-seq data from mouse brain cortex, HuMMuS enabled to accurately cluster scRNA profiles and to identify potential driver TFs.

AVAILABILITY AND IMPLEMENTATION

HuMMuS is available at https://github.com/cantinilab/HuMMuS.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

During gene regulation, DNA accessibility is thought to limit the availability of transcription factor (TF) binding sites, while TFs can increase DNA accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events in the modulation of gene expression remain unknown for the vast majority of genes. We utilized deeply sequenced ATAC-Seq data and site-specific knock-in reporter genes to investigate the relationship between the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of Cebpa during macrophage-neutrophil differentiation. While the enhancers upregulate reporter expression during the earliest stages of differentiation, there is little corresponding increase in their total accessibility. Conversely, total accessibility peaks during the last stages of differentiation without any increase in enhancer activity. The accessibility of positions neighboring C/EBP-family TF binding sites, which indicates TF occupancy, does increase significantly during early differentiation, showing that the early upregulation of enhancer activity is driven by TF binding. These results imply that a generalized increase in DNA accessibility is not sufficient, and binding by enhancer-specific TFs is necessary, for the upregulation of gene expression. Additionally, high-coverage ATAC-Seq combined with time-series expression data can infer the sequence of regulatory events at binding-site resolution.

An enhancer-AAV approach selectively targeting dentate granule cells of the mouse hippocampus

The mammalian brain contains a diverse array of cell types, including dozens of neuronal subtypes with distinct anatomical and functional characteristics. The brain leverages these neuron-type specializations to perform diverse circuit operations and thus execute different behaviors properly. Through the use of Cre lines, access to specific neuron types has improved over past decades. Despite their extraordinary utility, development and cross-breeding of Cre lines is time consuming and expensive, presenting a significant barrier to entry for investigators. Furthermore, cell-based therapeutics developed in Cre mice are not clinically translatable. Recently, several adeno-associated virus (AAV) vectors utilizing neuron-type-specific regulatory transcriptional sequences (enhancer-AAVs) were developed that overcome these limitations. Using a publicly available RNA sequencing (RNA-seq) dataset, we evaluated the potential of several candidate enhancers for neuron-type-specific targeting in the hippocampus. Here, we demonstrate that a previously identified enhancer-AAV selectively targets dentate granule cells over other excitatory neuron types in the hippocampus of wild-type adult mice.

Multi-omic profiling of the developing human cerebral cortex at the single-cell level

The cellular complexity of the human brain is established via dynamic changes in gene expression throughout development that is mediated, in part, by the spatiotemporal activity of cis-regulatory elements (CREs). We simultaneously profiled gene expression and chromatin accessibility in 45,549 cortical nuclei across six broad developmental time points from fetus to adult. We identified cell type-specific domains in which chromatin accessibility is highly correlated with gene expression. Differentiation pseudotime trajectory analysis indicates that chromatin accessibility at CREs precedes transcription and that dynamic changes in chromatin structure play a critical role in neuronal lineage commitment. In addition, we mapped cell type-specific and temporally specific genetic loci implicated in neuropsychiatric traits, including schizophrenia and bipolar disorder. Together, our results describe the complex regulation of cell composition at critical stages in lineage determination and shed light on the impact of spatiotemporal alterations in gene expression on neuropsychiatric disease.

CUT homeobox genes: transcriptional regulation of neuronal specification and beyond

CUT homeobox genes represent a captivating gene class fulfilling critical functions in the development and maintenance of multiple cell types across a wide range of organisms. They belong to the larger group of homeobox genes, which encode transcription factors responsible for regulating gene expression patterns during development. CUT homeobox genes exhibit two distinct and conserved DNA binding domains, a homeodomain accompanied by one or more CUT domains. Numerous studies have shown the involvement of CUT homeobox genes in diverse developmental processes such as body axis formation, organogenesis, tissue patterning and neuronal specification. They govern these processes by exerting control over gene expression through their transcriptional regulatory activities, which they accomplish by a combination of classic and unconventional interactions with the DNA. Intriguingly, apart from their roles as transcriptional regulators, they also serve as accessory factors in DNA repair pathways through protein-protein interactions. They are highly conserved across species, highlighting their fundamental importance in developmental biology. Remarkably, evolutionary analysis has revealed that CUT homeobox genes have experienced an extraordinary degree of rearrangements and diversification compared to other classes of homeobox genes, including the emergence of a novel gene family in vertebrates. Investigating the functions and regulatory networks of CUT homeobox genes provides significant understanding into the molecular mechanisms underlying embryonic development and tissue homeostasis. Furthermore, aberrant expression or mutations in CUT homeobox genes have been associated with various human diseases, highlighting their relevance beyond developmental processes. This review will overview the well known roles of CUT homeobox genes in nervous system development, as well as their functions in other tissues across phylogeny.

A novel enhancer-AAV approach selectively targeting dentate granule cells

The mammalian brain contains the most diverse array of cell types of any organ, including dozens of neuronal subtypes with distinct anatomical and functional characteristics. The brain leverages these neuron-type-specializations to perform diverse circuit operations and thus execute different behaviors properly. Through the use of Cre lines, access to specific neuron types has steadily improved over past decades. Despite their extraordinary utility, development and cross-breeding of Cre lines is time-consuming and expensive, presenting a significant barrier to entry for many investigators. Furthermore, cell-based therapeutics developed in Cre mice are not clinically translatable. Recently, several AAV vectors utilizing neuron-type-specific regulatory transcriptional sequences (enhancer-AAVs) were developed which overcome these limitations. Using a publicly available RNAseq dataset, we evaluated the potential of several candidate enhancers for neuron-type-specific targeting in the hippocampus. Here we identified a promising enhancer-AAV for targeting dentate granule cells and validated its selectivity in wild-type adult mice.

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.

The contributions of DNA accessibility and transcription factor occupancy to enhancer activity during cellular differentiation

The upregulation of gene expression by enhancers depends upon the interplay between the binding of sequence-specific transcription factors (TFs) and DNA accessibility. DNA accessibility is thought to limit the ability of TFs to bind to their sites, while TFs can increase accessibility to recruit additional factors that upregulate gene expression. Given this interplay, the causative regulatory events underlying the modulation of gene expression during cellular differentiation remain unknown for the vast majority of genes. We investigated the binding-site resolution dynamics of DNA accessibility and the expression dynamics of the enhancers of an important neutrophil gene, Cebpa, during macrophage-neutrophil differentiation. Reporter genes were integrated in a site-specific manner in PUER cells, which are progenitors that can be differentiated into neutrophils or macrophages in vitro by activating the pan-leukocyte TF PU.1. Time series data show that two enhancers upregulate reporter expression during the first 48 hours of neutrophil differentiation. Surprisingly, there is little or no increase in the total accessibility, measured by ATAC-Seq, of the enhancers during the same time period. Conversely, total accessibility peaks 96 hrs after PU.1 activation-consistent with its role as a pioneer-but the enhancers do not upregulate gene expression. Combining deeply sequenced ATAC-Seq data with a new bias-correction method allowed the profiling of accessibility at single-nucleotide resolution and revealed protected regions in the enhancers that match all previously characterized TF binding sites and ChIP-Seq data. Although the accessibility of most positions does not change during early differentiation, that of positions neighboring TF binding sites, an indicator of TF occupancy, did increase significantly. The localized accessibility changes are limited to nucleotides neighboring C/EBP-family TF binding sites, showing that the upregulation of enhancer activity during early differentiation is driven by C/EBP-family TF binding. These results show that increasing the total accessibility of enhancers is not sufficient for upregulating their activity and other events such as TF binding are necessary for upregulation. Also, TF binding can cause upregulation without a perceptible increase in total accessibility. Finally, this study demonstrates the feasibility of comprehensively mapping individual TF binding sites as footprints using high coverage ATAC-Seq and inferring the sequence of events in gene regulation by combining with time-series gene expression data.

Remembering through the genome: the role of chromatin states in brain functions and diseases

Chromatin is the physical substrate of the genome that carries the DNA sequence and ensures its proper functions and regulation in the cell nucleus. While a lot is known about the dynamics of chromatin during programmed cellular processes such as development, the role of chromatin in experience-dependent functions remains not well defined. Accumulating evidence suggests that in brain cells, environmental stimuli can trigger long-lasting changes in chromatin structure and tri-dimensional (3D) organization that can influence future transcriptional programs. This review describes recent findings suggesting that chromatin plays an important role in cellular memory, particularly in the maintenance of traces of prior activity in the brain. Inspired by findings in immune and epithelial cells, we discuss the underlying mechanisms and the implications for experience-dependent transcriptional regulation in health and disease. We conclude by presenting a holistic view of chromatin as potential molecular substrate for the integration and assimilation of environmental information that may constitute a conceptual basis for future research.

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.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

Reorganization of postmitotic neuronal chromatin accessibility for maturation of serotonergic identity

Assembly of transcriptomes encoding unique neuronal identities requires selective accessibility of transcription factors to cis-regulatory sequences in nucleosome-embedded postmitotic chromatin. Yet, the mechanisms controlling postmitotic neuronal chromatin accessibility are poorly understood. Here, we show that unique distal enhancers define the Pet1 neuron lineage that generates serotonin (5-HT) neurons in mice. Heterogeneous single-cell chromatin landscapes are established early in postmitotic Pet1 neurons and reveal the putative regulatory programs driving Pet1 neuron subtype identities. Distal enhancer accessibility is highly dynamic as Pet1 neurons mature, suggesting the existence of regulatory factors that reorganize postmitotic neuronal chromatin. We find that Pet1 and Lmx1b control chromatin accessibility to select Pet1-lineage-specific enhancers for 5-HT neurotransmission. Additionally, these factors are required to maintain chromatin accessibility during early maturation suggesting that postmitotic neuronal open chromatin is unstable and requires continuous regulatory input. Together, our findings reveal postmitotic transcription factors that reorganize accessible chromatin for neuron specialization.

Genetic Regulation of Vertebrate Forebrain Development by Homeobox Genes

Forebrain development in vertebrates is regulated by transcription factors encoded by homeobox, bHLH and forkhead gene families throughout the progressive and overlapping stages of neural induction and patterning, regional specification and generation of neurons and glia from central nervous system (CNS) progenitor cells. Moreover, cell fate decisions, differentiation and migration of these committed CNS progenitors are controlled by the gene regulatory networks that are regulated by various homeodomain-containing transcription factors, including but not limited to those of the Pax (paired), Nkx, Otx (orthodenticle), Gsx/Gsh (genetic screened), and Dlx (distal-less) homeobox gene families. This comprehensive review outlines the integral role of key homeobox transcription factors and their target genes on forebrain development, focused primarily on the telencephalon. Furthermore, links of these transcription factors to human diseases, such as neurodevelopmental disorders and brain tumors are provided.

Analysis of long and short enhancers in melanoma cell states

Understanding how enhancers drive cell-type specificity and efficiently identifying them is essential for the development of innovative therapeutic strategies. In melanoma, the melanocytic (MEL) and the mesenchymal-like (MES) states present themselves with different responses to therapy, making the identification of specific enhancers highly relevant. Using massively parallel reporter assays (MPRAs) in a panel of patient-derived melanoma lines (MM lines), we set to identify and decipher melanoma enhancers by first focusing on regions with state-specific H3K27 acetylation close to differentially expressed genes. An in-depth evaluation of those regions was then pursued by investigating the activity of overlapping ATAC-seq peaks along with a full tiling of the acetylated regions with 190 bp sequences. Activity was observed in more than 60% of the selected regions, and we were able to precisely locate the active enhancers within ATAC-seq peaks. Comparison of sequence content with activity, using the deep learning model DeepMEL2, revealed that AP-1 alone is responsible for the MES enhancer activity. In contrast, SOX10 and MITF both influence MEL enhancer function with SOX10 being required to achieve high levels of activity. Overall, our MPRAs shed light on the relationship between long and short sequences in terms of their sequence content, enhancer activity, and specificity across melanoma cell states.

Comprehensive understanding of Tn5 insertion preference improves transcription regulatory element identification

Tn5 transposase, which can efficiently tagment the genome, has been widely adopted as a molecular tool in next-generation sequencing, from short-read sequencing to more complex methods such as assay for transposase-accessible chromatin using sequencing (ATAC-seq). Here, we systematically map Tn5 insertion characteristics across several model organisms, finding critical parameters that affect its insertion. On naked genomic DNA, we found that Tn5 insertion is not uniformly distributed or random. To uncover drivers of these biases, we used a machine learning framework, which revealed that DNA shape cooperatively works with DNA motif to affect Tn5 insertion preference. These intrinsic insertion preferences can be modeled using nucleotide dependence information from DNA sequences, and we developed a computational pipeline to correct for these biases in ATAC-seq data. Using our pipeline, we show that bias correction improves the overall performance of ATAC-seq peak detection, recovering many potential false-negative peaks. Furthermore, we found that these peaks are bound by transcription factors, underscoring the biological relevance of capturing this additional information. These findings highlight the benefits of an improved understanding and precise correction of Tn5 insertion preference.

Cell-type-specific Hox regulatory strategies orchestrate tissue identity

Hox proteins are homeodomain transcription factors that diversify serially homologous segments along the animal body axis, as revealed by the classic bithorax phenotype of Drosophila melanogaster, in which mutations in Ultrabithorax (Ubx) transform the third thoracic segment into the likeness of the second thoracic segment. To specify segment identity, we show that Ubx both increases and decreases chromatin accessibility, coinciding with its dual role as both an activator and repressor of transcription. However, the choice of transcriptional activity executed by Ubx is spatially regulated and depends on the availability of cofactors, with Ubx acting as a repressor in some populations and as an activator in others. Ubx-mediated changes to chromatin accessibility positively and negatively affect the binding of Scalloped (Sd), a transcription factor that is required for appendage development in both segments. These findings illustrate how a single Hox protein can modify complex gene regulatory networks to transform the identity of an entire tissue.

Conserved and species-specific chromatin remodeling and regulatory dynamics during mouse and chicken limb bud development

Chromatin remodeling and genomic alterations impact spatio-temporal regulation of gene expression, which is central to embryonic development. The analysis of mouse and chicken limb development provides important insights into the morphoregulatory mechanisms, however little is known about the regulatory differences underlying their morphological divergence. Here, we identify the underlying shared and species-specific epigenomic and genomic variations. In mouse forelimb buds, we observe striking synchrony between the temporal dynamics of chromatin accessibility and gene expression, while their divergence in chicken wing buds uncovers species-specific regulatory heterochrony. In silico mapping of transcription factor binding sites and computational footprinting establishes the developmental time-restricted transcription factor-DNA interactions. Finally, the construction of target gene networks for HAND2 and GLI3 transcriptional regulators reveals both conserved and species-specific interactions. Our analysis reveals the impact of genome evolution on the regulatory interactions orchestrating vertebrate limb bud morphogenesis and provides a molecular framework for comparative Evo-Devo studies.

The vertebrate limb bud is a paradigm to uncover the fundamental mechanisms that govern embryogenesis and evolutionary diversification. Here the authors compare mouse and chicken limb bud development to study the impact of genome evolution on conserved and divergent gene regulatory interactions.

Consequences of NMDA receptor deficiency can be rescued in the adult brain

N-methyl-D-aspartate receptors (NMDARs) are required to shape activity-dependent connections in the developing and adult brain. Impaired NMDAR signalling through genetic or environmental insults causes a constellation of neurodevelopmental disorders that manifest as intellectual disability, epilepsy, autism, or schizophrenia. It is not clear whether the developmental impacts of NMDAR dysfunction can be overcome by interventions in adulthood. This question is paramount for neurodevelopmental disorders arising from mutations that occur in the GRIN genes, which encode NMDAR subunits, and the broader set of mutations that disrupt NMDAR function. We developed a mouse model where a congenital loss-of-function allele of Grin1 can be restored to wild type by gene editing with Cre recombinase. Rescue of NMDARs in adult mice yields surprisingly robust improvements in cognitive functions, including those that are refractory to treatment with current medications. These results suggest that neurodevelopmental disorders arising from NMDAR deficiency can be effectively treated in adults.

Chromatin accessibility in neuropsychiatric disorders

Epigenetic mechanisms have recently emerged as critical regulators of brain function in health and disease. By controlling the accessibility and the expression of specific genes, these pathways can mediate transient and long-lasting changes in neuronal function in both physiological and pathological contexts. Due to the extreme complexity of the epigenetic regulatory landscape, emerging methods that directly assay chromatin accessibility are of particular interest. Here, I review recent insights gained on open and closed chromatin states in the brain, with emphasis on neuropsychiatric disorders. These advances generated an invaluable wealth of information that can help us better understand gene regulation in the brain and its impairments that contribute to the development of disease.

Enhancer viruses for combinatorial cell-subclass-specific labeling

Rapid cell type identification by new genomic single-cell analysis methods has not been met with efficient experimental access to these cell types. To facilitate access to specific neural populations in mouse cortex, we collected chromatin accessibility data from individual cells and identified enhancers specific for cell subclasses and types. When cloned into recombinant adeno-associated viruses (AAVs) and delivered to the brain, these enhancers drive transgene expression in specific cortical cell subclasses. We extensively characterized several enhancer AAVs to show that they label different projection neuron subclasses as well as a homologous neuron subclass in human cortical slices. We also show how coupling enhancer viruses expressing recombinases to a newly generated transgenic mouse, Ai213, enables strong labeling of three different neuronal classes/subclasses in the brain of a single transgenic animal. This approach combines unprecedented flexibility with specificity for investigation of cell types in the mouse brain and beyond.

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.

Cell type-specific chromatin accessibility analysis in the mouse and human brain

The Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) is becoming popular in the neuroscience field where chromatin regulation is thought to be involved in neurodevelopment, activity-dependent gene regulation, hormonal and environmental responses, and pathophysiology of neuropsychiatric disorders. The advantages of using ATAC-seq include a small amount of material needed, fast protocol, and the ability to capture a range of gene regulatory elements with a single assay. With increasing interest in chromatin research, it is an imperative to have feasible, reliable assays that are compatible with a range of neuroscience study designs. Here we tested three protocols for neuronal chromatin accessibility analysis, including a varying brain tissue freezing method followed by fluorescence-activated nuclei sorting (FANS) and ATAC-seq. Our study shows that the cryopreservation method impacts the number of open chromatin regions identified from frozen brain tissue using ATAC-seq. However, we show that all protocols generate consistent and robust data and enable the identification of functional regulatory elements in neuronal cells. Our study implies that the broad biological interpretation of chromatin accessibility data is not significantly affected by the freezing condition. We also reveal additional challenges of doing chromatin analysis on post-mortem human brain tissue. Overall, ATAC-seq coupled with FANS is a powerful method to capture cell-type-specific chromatin accessibility information in mouse and human brain. Our study provides alternative brain preservation methods that generate high-quality ATAC-seq data while fitting in different study designs, and further encourages the use of this method to uncover the role of epigenetic (dys)regulation in the brain.

Spatially mapped single-cell chromatin accessibility

High-throughput single-cell epigenomic assays can resolve cell type heterogeneity in complex tissues, however, spatial orientation is lost. Here, we present single-cell combinatorial indexing on Microbiopsies Assigned to Positions for the Assay for Transposase Accessible Chromatin, or sciMAP-ATAC, as a method for highly scalable, spatially resolved, single-cell profiling of chromatin states. sciMAP-ATAC produces data of equivalent quality to non-spatial sci-ATAC and retains the positional information of each cell within a 214 micron cubic region, with up to hundreds of tracked positions in a single experiment. We apply sciMAP-ATAC to assess cortical lamination in the adult mouse primary somatosensory cortex and in the human primary visual cortex, where we produce spatial trajectories and integrate our data with non-spatial single-nucleus RNA and other chromatin accessibility single-cell datasets. Finally, we characterize the spatially progressive nature of cerebral ischemic infarction in the mouse brain using a model of transient middle cerebral artery occlusion.

Spatial orientation of cells in an interconnected network is lost in high-throughput single-cell epigenomic assays. Here the authors present sciMAP-ATAC to produce spatially resolved single-cell ATAC-seq data.

Identification of determinants of differential chromatin accessibility through a massively parallel genome-integrated reporter assay

A key mechanism in cellular regulation is the ability of the transcriptional machinery to physically access DNA. Transcription factors interact with DNA to alter the accessibility of chromatin, which enables changes to gene expression during development or disease or as a response to environmental stimuli. However, the regulation of DNA accessibility via the recruitment of transcription factors is difficult to study in the context of the native genome because every genomic site is distinct in multiple ways. Here we introduce the multiplexed integrated accessibility assay (MIAA), an assay that measures chromatin accessibility of synthetic oligonucleotide sequence libraries integrated into a controlled genomic context with low native accessibility. We apply MIAA to measure the effects of sequence motifs on cell type-specific accessibility between mouse embryonic stem cells and embryonic stem cell-derived definitive endoderm cells, screening 7905 distinct DNA sequences. MIAA recapitulates differential accessibility patterns of 100-nt sequences derived from natively differential genomic regions, identifying E-box motifs common to epithelial-mesenchymal transition driver transcription factors in stem cell-specific accessible regions that become repressed in endoderm. We show that a single binding motif for a key regulatory transcription factor is sufficient to open chromatin, and classify sets of stem cell-specific, endoderm-specific, and shared accessibility-modifying transcription factor motifs. We also show that overexpression of two definitive endoderm transcription factors, T and Foxa2, results in changes to accessibility in DNA sequences containing their respective DNA-binding motifs and identify preferential motif arrangements that influence accessibility.

Cell-Type Specificity of Genomic Imprinting in Cerebral Cortex

In mammalian genomes, a subset of genes is regulated by genomic imprinting, resulting in silencing of one parental allele. Imprinting is essential for cerebral cortex development, but prevalence and functional impact in individual cells is unclear. Here, we determined allelic expression in cortical cell types and established a quantitative platform to interrogate imprinting in single cells. We created cells with uniparental chromosome disomy (UPD) containing two copies of either the maternal or the paternal chromosome; hence, imprinted genes will be 2-fold overexpressed or not expressed. By genetic labeling of UPD, we determined cellular phenotypes and transcriptional responses to deregulated imprinted gene expression at unprecedented single-cell resolution. We discovered an unexpected degree of cell-type specificity and a novel function of imprinting in the regulation of cortical astrocyte survival. More generally, our results suggest functional relevance of imprinted gene expression in glial astrocyte lineage and thus for generating cortical cell-type diversity.

Transcription factor expression defines subclasses of developing projection neurons highly similar to single-cell RNA-seq subtypes

We are only just beginning to catalog the vast diversity of cell types in the cerebral cortex. Such categorization is a first step toward understanding how diversification relates to function. All cortical projection neurons arise from a uniform pool of progenitor cells that lines the ventricles of the forebrain. It is still unclear how these progenitor cells generate the more than 50 unique types of mature cortical projection neurons defined by their distinct gene-expression profiles. Moreover, exactly how and when neurons diversify their function during development is unknown. Here we relate gene expression and chromatin accessibility of two subclasses of projection neurons with divergent morphological and functional features as they develop in the mouse brain between embryonic day 13 and postnatal day 5 in order to identify transcriptional networks that diversify neuron cell fate. We compare these gene-expression profiles with published profiles of single cells isolated from similar populations and establish that layer-defined cell classes encompass cell subtypes and developmental trajectories identified using single-cell sequencing. Given the depth of our sequencing, we identify groups of transcription factors with particularly dense subclass-specific regulation and subclass-enriched transcription factor binding motifs. We also describe transcription factor-adjacent long noncoding RNAs that define each subclass and validate the function of Myt1l in balancing the ratio of the two subclasses in vitro. Our multidimensional approach supports an evolving model of progressive restriction of cell fate competence through inherited transcriptional identities.

Integrative analysis of in vivo recording with single-cell RNA-seq data reveals molecular properties of light-sensitive neurons in mouse V1

Vision formation is classically based on projections from retinal ganglion cells (RGC) to the lateral geniculate nucleus (LGN) and the primary visual cortex (V1). Neurons in the mouse V1 are tuned to light stimuli. Although the cellular information of the retina and the LGN has been widely studied, the transcriptome profiles of single light-stimulated neuron in V1 remain unknown. In our study, in vivo calcium imaging and whole-cell electrophysiological patch-clamp recording were utilized to identify 53 individual cells from layer 2/3 of V1 as light-sensitive (LS) or non-light-sensitive (NS) by single-cell light-evoked calcium evaluation and action potential spiking. The contents of each cell after functional tests were aspirated in vivo through a patch-clamp pipette for mRNA sequencing. Moreover, the three-dimensional (3-D) morphological characterizations of the neurons were reconstructed in a live mouse after the whole-cell recordings. Our sequencing results indicated that V1 neurons with a high expression of genes related to transmission regulation, such as Rtn4r and Rgs7, and genes involved in membrane transport, such as Na⁺/K⁺ ATPase and NMDA-type glutamatergic receptors, preferentially responded to light stimulation. Furthermore, an antagonist that blocks Rtn4r signals could inactivate the neuronal responses to light stimulation in live mice. In conclusion, our findings of the vivo-seq analysis indicate the key role of the strength of synaptic transmission possesses neurons in V1 of light sensory.

Leveraging mouse chromatin data for heritability enrichment informs common disease architecture and reveals cortical layer contributions to schizophrenia

Genome-wide association studies have implicated thousands of noncoding variants across common human phenotypes. However, they cannot directly inform the cellular context in which disease-associated variants act. Here, we use open chromatin profiles from discrete mouse cell populations to address this challenge. We applied stratified linkage disequilibrium score regression and evaluated heritability enrichment in 64 genome-wide association studies, emphasizing schizophrenia. We provide evidence that mouse-derived human open chromatin profiles can serve as powerful proxies for difficult to obtain human cell populations, facilitating the illumination of common disease heritability enrichment across an array of human phenotypes. We demonstrate that signatures from discrete subpopulations of cortical excitatory and inhibitory neurons are significantly enriched for schizophrenia heritability with maximal enrichment in cortical layer V excitatory neurons. We also show that differences between schizophrenia and bipolar disorder are concentrated in excitatory neurons in cortical layers II-III, IV, and V, as well as the dentate gyrus. Finally, we leverage these data to fine-map variants in 177 schizophrenia loci nominating variants in 104/177. We integrate these data with transcription factor binding site, chromatin interaction, and validated enhancer data, placing variants in the cellular context where they may modulate risk.

Epigenetic regulation for acquiring glial identity by neural stem cells during cortical development

The nervous system consists of several hundred neuronal subtypes and glial cells that show specific gene expression and are generated from common ancestors, neural stem cells (NSCs). As the experimental techniques and molecular tools to analyze epigenetics and chromatin structures are rapidly advancing, the comprehensive events and genome-wide states of DNA methylation, histone modifications, and chromatin accessibility in developing NSCs are gradually being unveiled. Here, we review recent advances in elucidating the role of epigenetic and chromatin regulation in NSCs, especially focusing on the acquisition of glial identity and how epigenetic regulation enables the temporal regulation of NSCs during murine cortical development.

Poor Concordance of Floxed Sequence Recombination in Single Neural Stem Cells: Implications for Cell Autonomous Studies

To manipulate target gene function in specific adult cell populations, tamoxifen (TAM)-dependent CreER^T2 is widely used to drive inducible, site-specific recombination of loxP flanked sequences. In studies of cell autonomous target gene function, it is common practice to combine these CreER^T2-lox systems with a ubiquitously expressed stop-floxed fluorescent reporter gene to identify single cells supposedly undergoing target gene recombination. Here, we studied the reliability of using Cre-induced recombination of one gene to predict recombination in another gene at the single-cell level in adult hippocampal neural stem and progenitor cells (NSPCs). Using both probabilistic predictions in a generic experimental paradigm, as well as a mouse model with two separate stop-floxed reporters plus a Nestin promoter-driven CreER^T2, we found that, in individual cells, recombination of one gene was a poor predictor of recombination in another. This poor concordance in floxed sequence recombination across genes suggests that use of stop-floxed reporters to investigate cell autonomous gene function may not be universally reliable and could lead to false conclusions.

Development and Arealization of the Cerebral Cortex

Adult cortical areas consist of specialized cell types and circuits that support unique higher-order cognitive functions. How this regional diversity develops from an initially uniform neuroepithelium has been the subject of decades of seminal research, and emerging technologies, including single-cell transcriptomics, provide a new perspective on area-specific molecular diversity. Here, we review the early developmental processes that underlie cortical arealization, including both cortex intrinsic and extrinsic mechanisms as embodied by the protomap and protocortex hypotheses, respectively. We propose an integrated model of serial homology whereby intrinsic genetic programs and local factors establish early transcriptomic differences between excitatory neurons destined to give rise to broad "proto-regions," and activity-dependent mechanisms lead to progressive refinement and formation of sharp boundaries between functional areas. Finally, we explore the potential of these basic developmental processes to inform our understanding of the emergence of functional neural networks and circuit abnormalities in neurodevelopmental disorders.

Chromatin Decondensation by FOXP2 Promotes Human Neuron Maturation and Expression of Neurodevelopmental Disease Genes

Forkhead box P2 (FOXP2) is a transcription factor expressed in the human brain that peaks during fetal development, and disruption in its ability to regulate downstream target genes leads to vulnerability to neurodevelopmental disorders. However, the mechanisms by which FOXP2 exerts regulatory control over targets during neuronal maturation have not been fully elucidated. Here, we use genome-wide chromatin accessibility assays and transcriptome-wide expression analyses in differentiating human neurons to show that FOXP2 represses proliferation-promoting genes in a DNA-binding-dependent manner. In contrast, FOXP2 and its cofactors, NFIA and NFIB, activate neuronal maturation genes in a manner that does not require FOXP2 to interact with DNA directly. Moreover, comparisons with expression data from the developing human brain suggest that FOXP2 and NFIA- or NFIB-dependent chromatin alterations drive maturation of excitatory cortical neurons. Thus, FOXP2 and its NFI cofactors may be specifically important for the development of cortical circuits underlying neurodevelopmental disorders.

cisTopic: cis-regulatory topic modeling on single-cell ATAC-seq data

We present cisTopic, a probabilistic framework used to simultaneously discover coaccessible enhancers and stable cell states from sparse single-cell epigenomics data ( http://github.com/aertslab/cistopic ). Using a compendium of single-cell ATAC-seq datasets from differentiating hematopoietic cells, brain and transcription factor perturbations, we demonstrate that topic modeling can be exploited for robust identification of cell types, enhancers and relevant transcription factors. cisTopic provides insight into the mechanisms underlying regulatory heterogeneity in cell populations.

Combined analysis of dissimilar promoter accessibility and gene expression profiles identifies tissue-specific genes and actively repressed networks

BACKGROUND

The assay for transposase-accessible chromatin (ATAC-seq) is a powerful method to examine chromatin accessibility. While many studies have reported a positive correlation between gene expression and promoter accessibility, few have investigated the genes that deviate from this trend. In this study, we aimed to understand the relationship between gene expression and promoter accessibility in multiple cell types while also identifying gene regulatory networks in the placenta, an understudied organ that is critical for a successful pregnancy.

RESULTS

We started by assaying the open chromatin landscape in the mid-gestation placenta, when the fetal vasculature has started developing. After incorporating transcriptomic data generated in the placenta at the same time point, we grouped genes based on their expression levels and ATAC-seq promoter coverage. We found that the genes with the strongest correlation (high expression and high coverage) are likely involved in housekeeping functions, whereas tissue-specific genes were highly expressed and had only medium-low coverage. We also predicted that genes with medium-low expression and high promoter coverage were actively repressed. Within this group, we extracted a protein-protein interaction network enriched for neuronal functions, likely preventing the cells from adopting a neuronal fate. We further confirmed that a repressive histone mark is bound to the promoters of genes in this network. Finally, we ran our pipeline using ATAC-seq and RNA-seq data generated in ten additional cell types. We again found that genes with the strongest correlation are enriched for housekeeping functions and that genes with medium-low promoter coverage and high expression are more likely to be tissue-specific. These results demonstrate that only two data types, both of which require relatively low starting material to generate and are becoming more commonly available, can be integrated to understand multiple aspects of gene regulation.

CONCLUSIONS

Within the placenta, we identified an active placenta-specific gene network as well as a repressed neuronal network. Beyond the placenta, we demonstrate that ATAC-seq data and RNA-seq data can be integrated to identify tissue-specific genes and actively repressed gene networks in multiple cell types.

SYNGAP1 heterozygosity disrupts sensory processing by reducing touch-related activity within somatosensory cortex circuits

In addition to cognitive impairments, neurodevelopmental disorders often result in sensory processing deficits. However, the biological mechanisms that underlie impaired sensory processing associated with neurodevelopmental disorders are generally understudied and poorly understood. We found that SYNGAP1 haploinsufficiency in humans, which causes a sporadic neurodevelopmental disorder defined by cognitive impairment, autistic features, and epilepsy, also leads to deficits in tactile-related sensory processing. In vivo neurophysiological analysis in Syngap1 mouse models revealed that upper-lamina neurons in somatosensory cortex weakly encode information related to touch. This was caused by reduced synaptic connectivity and impaired intrinsic excitability within upper-lamina somatosensory cortex neurons. These results were unexpected, given that Syngap1 heterozygosity is known to cause circuit hyperexcitability in brain areas more directly linked to cognitive functions. Thus, Syngap1 heterozygosity causes a range of circuit-specific pathologies, including reduced activity within cortical neurons required for touch processing, which may contribute to sensory phenotypes observed in patients.

Developmental Chromatin Restriction of Pro-Growth Gene Networks Acts as an Epigenetic Barrier to Axon Regeneration in Cortical Neurons

Axon regeneration in the central nervous system is prevented in part by a developmental decline in the intrinsic regenerative ability of maturing neurons. This loss of axon growth ability likely reflects widespread changes in gene expression, but the mechanisms that drive this shift remain unclear. Chromatin accessibility has emerged as a key regulatory mechanism in other cellular contexts, raising the possibility that chromatin structure may contribute to the age-dependent loss of regenerative potential. Here we establish an integrated bioinformatic pipeline that combines analysis of developmentally dynamic gene networks with transcription factor regulation and genome-wide maps of chromatin accessibility. When applied to the developing cortex, this pipeline detected overall closure of chromatin in sub-networks of genes associated with axon growth. We next analyzed mature CNS neurons that were supplied with various pro-regenerative transcription factors. Unlike prior results with SOX11 and KLF7, here we found that neither JUN nor an activated form of STAT3 promoted substantial corticospinal tract regeneration. Correspondingly, chromatin accessibility in JUN or STAT3 target genes was substantially lower than in predicted targets of SOX11 and KLF7. Finally, we used the pipeline to predict pioneer factors that could potentially relieve chromatin constraints at growth-associated loci. Overall this integrated analysis substantiates the hypothesis that dynamic chromatin accessibility contributes to the developmental decline in axon growth ability and influences the efficacy of pro-regenerative interventions in the adult, while also pointing toward selected pioneer factors as high-priority candidates for future combinatorial experiments. © 2018 Wiley Periodicals, Inc. Develop Neurobiol 00: 000-000, 2018.

Marked Diversity of Unique Cortical Enhancers Enables Neuron-Specific Tools by Enhancer-Driven Gene Expression

Understanding neural circuit function requires individually addressing their component parts: specific neuronal cell types. However, not only do the precise genetic mechanisms specifying neuronal cell types remain obscure, access to these neuronal cell types by transgenic techniques also remains elusive. Whereas most genes are expressed in the brain, the vast majority are expressed in many different kinds of neurons, suggesting that promoters alone are not sufficiently specific to distinguish cell types. However, there are orders of magnitude more distal genetic cis-regulatory elements controlling transcription (i.e., enhancers), so we screened for enhancer activity in microdissected samples of mouse cortical subregions. This identified thousands of novel putative enhancers, many unique to particular cortical subregions. Pronuclear injection of expression constructs containing such region-specific enhancers resulted in transgenic lines driving expression in distinct sets of cells specifically in the targeted cortical subregions, even though the parent gene's promoter was relatively non-specific. These data showcase the promise of utilizing the genetic mechanisms underlying the specification of diverse neuronal cell types for the development of genetic tools potentially capable of targeting any neuronal circuit of interest, an approach we call enhancer-driven gene expression (EDGE).

Specialized Subpopulations of Deep-Layer Pyramidal Neurons in the Neocortex: Bridging Cellular Properties to Functional Consequences

Neocortical pyramidal neurons with somata in layers 5 and 6 are among the most visually striking and enigmatic neurons in the brain. These deep-layer pyramidal neurons (DLPNs) integrate a plethora of cortical and extracortical synaptic inputs along their impressive dendritic arbors. The pattern of cortical output to both local and long-distance targets is sculpted by the unique physiological properties of specific DLPN subpopulations. Here we revisit two broad DLPN subpopulations: those that send their axons within the telencephalon (intratelencephalic neurons) and those that project to additional target areas outside the telencephalon (extratelencephalic neurons). While neuroscientists across many subdisciplines have characterized the intrinsic and synaptic physiological properties of DLPN subpopulations, our increasing ability to selectively target and manipulate these output neuron subtypes advances our understanding of their distinct functional contributions. This Viewpoints article summarizes our current knowledge about DLPNs and highlights recent work elucidating the functional differences between DLPN subpopulations.

Gene regulatory mechanisms underlying sex differences in brain development and psychiatric disease

The sexual differentiation of the mammalian nervous system requires the precise coordination of the temporal and spatial regulation of gene expression in diverse cell types. Sex hormones act at multiple developmental time points to specify sex-typical differentiation during embryonic and early development and to coordinate subsequent responses to gonadal hormones later in life by establishing sex-typical patterns of epigenetic modifications across the genome. Thus, mutations associated with neuropsychiatric conditions may result in sexually dimorphic symptoms by acting on different neural substrates or chromatin landscapes in males and females. Finally, as stress hormone signaling may directly alter the molecular machinery that interacts with sex hormone receptors to regulate gene expression, the contribution of chronic stress to the pathogenesis or presentation of mental illness may be additionally different between the sexes. Here, we review the mechanisms that contribute to sexual differentiation in the mammalian nervous system and consider some of the implications of these processes for sex differences in neuropsychiatric conditions.

Single-nucleus analysis of accessible chromatin in developing mouse forebrain reveals cell-type-specific transcriptional regulation

Analysis of chromatin accessibility can reveal transcriptional regulatory sequences, but heterogeneity of primary tissues poses a significant challenge in mapping the precise chromatin landscape in specific cell types. Here we report single-nucleus ATAC-seq, a combinatorial barcoding-assisted single-cell assay for transposase-accessible chromatin that is optimized for use on flash-frozen primary tissue samples. We apply this technique to the mouse forebrain through eight developmental stages. Through analysis of more than 15,000 nuclei, we identify 20 distinct cell populations corresponding to major neuronal and non-neuronal cell types. We further define cell-type-specific transcriptional regulatory sequences, infer potential master transcriptional regulators and delineate developmental changes in forebrain cellular composition. Our results provide insight into the molecular and cellular dynamics that underlie forebrain development in the mouse and establish technical and analytical frameworks that are broadly applicable to other heterogeneous tissues.

The Dynamic Landscape of Open Chromatin during Human Cortical Neurogenesis

Non-coding regions comprise most of the human genome and harbor a significant fraction of risk alleles for neuropsychiatric diseases, yet their functions remain poorly defined. We created a high-resolution map of non-coding elements involved in human cortical neurogenesis by contrasting chromatin accessibility and gene expression in the germinal zone and cortical plate of the developing cerebral cortex. We link distal regulatory elements (DREs) to their cognate gene(s) together with chromatin interaction data and show that target genes of human-gained enhancers (HGEs) regulate cortical neurogenesis and are enriched in outer radial glia, a cell type linked to human cortical evolution. We experimentally validate the regulatory effects of predicted enhancers for FGFR2 and EOMES. We observe that common genetic variants associated with educational attainment, risk for neuropsychiatric disease, and intracranial volume are enriched within regulatory elements involved in cortical neurogenesis, demonstrating the importance of this early developmental process for adult human cognitive function.

The BRAIN Initiative Cell Census Consortium: Lessons Learned toward Generating a Comprehensive Brain Cell Atlas

A comprehensive characterization of neuronal cell types, their distributions, and patterns of connectivity is critical for understanding the properties of neural circuits and how they generate behaviors. Here we review the experiences of the BRAIN Initiative Cell Census Consortium, ten pilot projects funded by the U.S. BRAIN Initiative, in developing, validating, and scaling up emerging genomic and anatomical mapping technologies for creating a complete inventory of neuronal cell types and their connections in multiple species and during development. These projects lay the foundation for a larger and longer-term effort to generate whole-brain cell atlases in species including mice and humans.

An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues

We present Omni-ATAC, an improved ATAC-seq protocol for chromatin accessibility profiling that works across multiple applications with substantial improvement of signal-to-background ratio and information content. The Omni-ATAC protocol generates chromatin accessibility profiles from archival frozen tissue samples and 50-μm sections, revealing the activities of disease-associated DNA elements in distinct human brain structures. The Omni-ATAC protocol enables the interrogation of personal regulomes in tissue context and translational studies.