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.

Critical test of the assumption that the hypothalamic entopeduncular nucleus of rodents is homologous with the primate internal pallidum.

The globus pallidus (GP) of primates is divided conventionally into distinct internal and external parts. The literature repeats since 1930 the opinion that the homolog of the primate internal pallidum in rodents is the hypothalamic entopeduncular nucleus (embedded within fiber tracts of the cerebral peduncle). To test this idea, we explored its historic fundaments, checked the development and genoarchitecture of mouse entopeduncular and pallidal neurons, and examined relevant comparative connectivity data. We found that the extratelencephalic mouse entopeduncular structure consists of four different components arrayed along a dorsoventral sequence in the alar hypothalamus. The ventral entopeduncular nucleus (EPV), with GABAergic neurons expressing Dlx5&6 and Nkx2-1, lies within the hypothalamic peduncular subparaventricular area. Three other formations-the dorsal entopeduncular nucleus (EPD), the prereticular entopeduncular nucleus (EPPRt ), and the preeminential entopeduncular nucleus (EPPEm )-lie within the overlying paraventricular area, under the subpallium. EPD contains glutamatergic neurons expressing Tbr1, Otp, and Pax6. The EPPRt has GABAergic cells expressing Isl1 and Meis2, whereas the EPPEm population expresses Foxg1 and may be glutamatergic. Genoarchitectonic observations on relevant areas of the mouse pallidal/diagonal subpallium suggest that the GP of rodents is constituted as in primates by two adjacent but molecularly and hodologically differentiable telencephalic portions (both expressing Foxg1). These and other reported data oppose the notion that the rodent extratelencephalic entopeduncular nucleus is homologous to the primate internal pallidum. We suggest instead that all mammals, including rodents, have dual subpallial GP components, whereas primates probably also have a comparable set of hypothalamic entopeduncular nuclei. Remarkably, there is close similarity in some gene expression properties of the telencephalic internal GP and the hypothalamic EPV. This apparently underlies their notable functional analogy, sharing GABAergic neurons and thalamopetal connectivity.

The transcription regulator Lmo3 is required for the development of medial ganglionic eminence derived neurons in the external globus pallidus.

The external globus pallidus (GPe) is an essential component of the basal ganglia, a group of subcortical nuclei that are involved in control of action. Changes in the firing of GPe neurons are associated with both passive and active body movements. Aberrant activity of GPe neurons has been linked to motor symptoms of a variety of movement disorders, such as Parkinson's Disease, Huntington's disease and dystonia. Recent studies have helped delineate functionally distinct subtypes of GABAergic GPe projection neurons. However, not much is known about specific molecular mechanisms underlying the development of GPe neuronal subtypes. We show that the transcriptional regulator Lmo3 is required for the development of medial ganglionic eminence derived Nkx2.1+ and PV+ GPe neurons, but not lateral ganglionic eminence derived FoxP2+ neurons. As a consequence of the reduction in PV+ neurons, Lmo3-null mice have a reduced GPe input to the subthalamic nucleus.

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.

Machine learning dissection of human accelerated regions in primate neurodevelopment.

Using machine learning (ML), we interrogated the function of all human-chimpanzee variants in 2,645 human accelerated regions (HARs), finding 43% of HARs have variants with large opposing effects on chromatin state and 14% on neurodevelopmental enhancer activity. This pattern, consistent with compensatory evolution, was confirmed using massively parallel reporter assays in chimpanzee and human neural progenitor cells. The species-specific enhancer activity of HARs was accurately predicted from the presence and absence of transcription factor footprints in each species. Despite these striking cis effects, activity of a given HAR sequence was nearly identical in human and chimpanzee cells. This suggests that HARs did not evolve to compensate for changes in the trans environment but instead altered their ability to bind factors present in both species. Thus, ML prioritized variants with functional effects on human neurodevelopment and revealed an unexpected reason why HARs may have evolved so rapidly.

Dorsoventral Arrangement of Lateral Hypothalamus Populations in the Mouse Hypothalamus: a Prosomeric Genoarchitectonic Analysis.

The lateral hypothalamus (LH) has a heterogeneous cytoarchitectonic organization that has not been elucidated in detail. In this work, we analyzed within the framework of the prosomeric model the differential expression pattern of 59 molecular markers along the ventrodorsal dimension of the medial forebrain bundle in the mouse, considering basal and alar plate subregions of the LH. We found five basal (LH1-LH5) and four alar (LH6-LH9) molecularly distinct sectors of the LH with neuronal cell groups that correlate in topography with previously postulated alar and basal hypothalamic progenitor domains. Most peptidergic populations were restricted to one of these LH sectors though some may have dispersed into a neighboring sector. For instance, histaminergic Hdc-positive neurons were mostly contained within the basal LH3, Nts (neurotensin)- and Tac2 (tachykinin 2)-expressing cells lie strictly within LH4, Hcrt (hypocretin/orexin)-positive and Pmch (pro-melanin-concentrating hormone)-positive neurons appeared within separate LH5 subdivisions, Pnoc (prepronociceptin)-expressing cells were mainly restricted to LH6, and Sst (somatostatin)-positive cells were identified within the LH7 sector. The alar LH9 sector, a component of the Foxg1-positive telencephalo-opto-hypothalamic border region, selectively contained Satb2-expressing cells. Published studies of rodent LH subdivisions have not described the observed pattern. Our genoarchitectonic map should aid in systematic approaches to elucidate LH connectivity and function.

Tsc1 Loss in VIP-Lineage Cortical Interneurons Results in More VIP+ Interneurons and Enhanced Excitability.

The mammalian target of rapamycin (mTOR) signaling pathway is a powerful regulator of cell proliferation, growth, synapse maintenance and cell fate. While intensely studied for its role in cancer, the role of mTOR signaling is just beginning to be uncovered in specific cell types that are implicated in neurodevelopmental disorders. Previously, loss of the Tsc1 gene, which results in hyperactive mTOR, was shown to affect the function and molecular properties of GABAergic cortical interneurons (CINs) derived from the medial ganglionic eminence. To assess if other important classes of CINs could be impacted by mTOR dysfunction, we deleted Tsc1 in a caudal ganglionic eminence-derived interneuron group, the vasoactive intestinal peptide (VIP)+ subtype, whose activity disinhibits local circuits. Tsc1 mutant VIP+ CINs reduced their pattern of apoptosis from postnatal days 15-20, resulting in increased VIP+ CINs. The mutant CINs exhibited synaptic and electrophysiological properties that could contribute to the high rate of seizure activity in humans that harbor Tsc1 mutations.

DLX1 and the NuRD complex cooperate in enhancer decommissioning and transcriptional repression.

In the developing subpallium, the fate decision between neurons and glia is driven by expression of Dlx1/2 or Olig1/2, respectively, two sets of transcription factors with a mutually repressive relationship. The mechanism by which Dlx1/2 repress progenitor and oligodendrocyte fate, while promoting transcription of genes needed for differentiation, is not fully understood. We identified a motif within DLX1 that binds RBBP4, a NuRD complex subunit. ChIP-seq studies of genomic occupancy of DLX1 and six different members of the NuRD complex show that DLX1 and NuRD colocalize to putative regulatory elements enriched near other transcription factor genes. Loss of Dlx1/2 leads to dysregulation of genome accessibility at putative regulatory elements near genes repressed by Dlx1/2, including Olig2. Consequently, heterozygosity of Dlx1/2 and Rbbp4 leads to an increase in the production of OLIG2+ cells. These findings highlight the importance of the interplay between transcription factors and chromatin remodelers in regulating cell-fate decisions.

Trans-Seq maps a selective mammalian retinotectal synapse instructed by Nephronectin.

The mouse visual system serves as an accessible model to understand mammalian circuit wiring. Despite rich knowledge in retinal circuits, the long-range connectivity map from distinct retinal ganglion cell (RGC) types to diverse brain neuron types remains unknown. In this study, we developed an integrated approach, called Trans-Seq, to map RGCs to superior collicular (SC) circuits. Trans-Seq combines a fluorescent anterograde trans-synaptic tracer, consisting of codon-optimized wheat germ agglutinin fused to mCherry, with single-cell RNA sequencing. We used Trans-Seq to classify SC neuron types innervated by genetically defined RGC types and predicted a neuronal pair from αRGCs to Nephronectin-positive wide-field neurons (NPWFs). We validated this connection using genetic labeling, electrophysiology and retrograde tracing. We then used transcriptomic data from Trans-Seq to identify Nephronectin as a determinant for selective synaptic choice from αRGC to NPWFs via binding to Integrin α8ß1. The Trans-Seq approach can be broadly applied for post-synaptic circuit discovery from genetically defined pre-synaptic neurons.

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

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.

LiCl treatment leads to long-term restoration of spine maturation and synaptogenesis in adult Tbr1 mutants.

BACKGROUND: Tbr1 encodes a T-box transcription factor and is considered a high confidence autism spectrum disorder (ASD) gene. Tbr1 is expressed in the postmitotic excitatory neurons of the deep neocortical layers 5 and 6. Postnatally and neonatally, Tbr1 conditional mutants (CKOs) have immature dendritic spines and reduced synaptic density. However, an understanding of Tbr1's function in the adult mouse brain remains elusive. METHODS: We used conditional mutagenesis to interrogate Tbr1's function in cortical layers 5 and 6 of the adult mouse cortex. RESULTS: Adult Tbr1 CKO mutants have dendritic spine and synaptic deficits as well as reduced frequency of mEPSCs and mIPSCs. LiCl, a WNT signaling agonist, robustly rescues the dendritic spine maturation, synaptic defects, and excitatory and inhibitory synaptic transmission deficits. CONCLUSIONS: LiCl treatment could be used as a therapeutic approach for some cases of ASD with deficits in synaptic transmission.

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.

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.

Autism risk gene POGZ promotes chromatin accessibility and expression of clustered synaptic genes.

Deleterious genetic variants in POGZ, which encodes the chromatin regulator Pogo Transposable Element with ZNF Domain protein, are strongly associated with autism spectrum disorder (ASD). Although it is a high-confidence ASD risk gene, the neurodevelopmental functions of POGZ remain unclear. Here we reveal the genomic binding of POGZ in the developing forebrain at euchromatic loci and gene regulatory elements (REs). We profile chromatin accessibility and gene expression in Pogz-/- mice and show that POGZ promotes the active chromatin state and transcription of clustered synaptic genes. We further demonstrate that POGZ forms a nuclear complex and co-occupies loci with ADNP, another high-confidence ASD risk gene, and provide evidence that POGZ regulates other neurodevelopmental disorder risk genes as well. Our results reveal a neurodevelopmental function of an ASD risk gene and identify molecular targets that may elucidate its function in ASD.

Developmental dynamics of voltage-gated sodium channel isoform expression in the human and mouse brain.

BACKGROUND: Genetic variants in the voltage-gated sodium channels SCN1A, SCN2A, SCN3A, and SCN8A are leading causes of epilepsy, developmental delay, and autism spectrum disorder. The mRNA splicing patterns of all four genes vary across development in the rodent brain, including mutually exclusive copies of the fifth protein-coding exon detected in the neonate (5N) and adult (5A). A second pair of mutually exclusive exons is reported in SCN8A only (18N and 18A). We aimed to quantify the expression of individual exons in the developing human brain. METHODS: RNA-seq data from 783 human brain samples across development were analyzed to estimate exon-level expression. Developmental changes in exon utilization were validated by assessing intron splicing. Exon expression was also estimated in RNA-seq data from 58 developing mouse neocortical samples. RESULTS: In the mature human neocortex, exon 5A is consistently expressed at least 4-fold higher than exon 5N in all four genes. For SCN2A, SCN3A, and SCN8A, a brain-wide synchronized 5N to 5A transition occurs between 24 post-conceptual weeks (2nd trimester) and 6 years of age. In mice, the equivalent 5N to 5A transition begins at or before embryonic day 15.5. In SCN8A, over 90% of transcripts in the mature human cortex include exon 18A. Early in fetal development, most transcripts include 18N or skip both 18N and 18A, with a transition to 18A inclusion occurring from 13 post-conceptual weeks to 6 months of age. No other protein-coding exons showed comparably dynamic developmental trajectories. CONCLUSIONS: Exon usage in SCN1A, SCN2A, SCN3A, and SCN8A changes dramatically during human brain development. These splice isoforms, which alter the biophysical properties of the encoded channels, may account for some of the observed phenotypic differences across development and between specific variants. Manipulation of the proportion of splicing isoforms at appropriate stages of development may act as a therapeutic strategy for specific mutations or even epilepsy in general.

COUP-TFI specifies the medial entorhinal cortex identity and induces differential cell adhesion to determine the integrity of its boundary with neocortex.

Development of cortical regions with precise, sharp, and regular boundaries is essential for physiological function. However, little is known of the mechanisms ensuring these features. Here, we show that determination of the boundary between neocortex and medial entorhinal cortex (MEC), two abutting cortical regions generated from the same progenitor lineage, relies on COUP-TFI (chicken ovalbumin upstream promoter-transcription factor I), a patterning transcription factor with graded expression in cortical progenitors. In contrast with the classical paradigm, we found that increased COUP-TFI expression expands MEC, creating protrusions and disconnected ectopic tissue. We further developed a mathematical model that predicts that neuronal specification and differential cell affinity contribute to the emergence of an instability region and boundary sharpness. Correspondingly, we demonstrated that high expression of COUP-TFI induces MEC cell fate and protocadherin 19 expression. Thus, we conclude that a sharp boundary requires a subtle interplay between patterning transcription factors and differential cell affinity.

Ultraconserved enhancer function does not require perfect sequence conservation.

Ultraconserved enhancer sequences show perfect conservation between human and rodent genomes, suggesting that their functions are highly sensitive to mutation. However, current models of enhancer function do not sufficiently explain this extreme evolutionary constraint. We subjected 23 ultraconserved enhancers to different levels of mutagenesis, collectively introducing 1,547 mutations, and examined their activities in transgenic mouse reporter assays. Overall, we find that the regulatory properties of ultraconserved enhancers are robust to mutation. Upon mutagenesis, nearly all (19/23, 83%) still functioned as enhancers at one developmental stage, as did most of those tested again later in development (5/9, 56%). Replacement of endogenous enhancers with mutated alleles in mice corroborated results of transgenic assays, including the functional resilience of ultraconserved enhancers to mutation. Our findings show that the currently known activities of ultraconserved enhancers do not necessarily require the perfect conservation observed in evolution and suggest that additional regulatory or other functions contribute to their sequence constraint.