Genetic Control of Expression and Splicing in Developing Human Brain Informs Disease Mechanisms.

Tissue-specific regulatory regions harbor substantial genetic risk for disease. Because brain development is a critical epoch for neuropsychiatric disease susceptibility, we characterized the genetic control of the transcriptome in 201 mid-gestational human brains, identifying 7,962 expression quantitative trait loci (eQTL) and 4,635 spliceQTL (sQTL), including several thousand prenatal-specific regulatory regions. We show that significant genetic liability for neuropsychiatric disease lies within prenatal eQTL and sQTL. Integration of eQTL and sQTL with genome-wide association studies (GWAS) via transcriptome-wide association identified dozens of novel candidate risk genes, highlighting shared and stage-specific mechanisms in schizophrenia (SCZ). Gene network analysis revealed that SCZ and autism spectrum disorder (ASD) affect distinct developmental gene co-expression modules. Yet, in each disorder, common and rare genetic variation converges within modules, which in ASD implicates superficial cortical neurons. More broadly, these data, available as a web browser and our analyses, demonstrate the genetic mechanisms by which developmental events have a widespread influence on adult anatomical and behavioral phenotypes.

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.

Brain-trait-associated variants impact cell-type-specific gene regulation during neurogenesis.

Interpretation of the function of non-coding risk loci for neuropsychiatric disorders and brain-relevant traits via gene expression and alternative splicing quantitative trait locus (e/sQTL) analyses is generally performed in bulk post-mortem adult tissue. However, genetic risk loci are enriched in regulatory elements active during neocortical differentiation, and regulatory effects of risk variants may be masked by heterogeneity in bulk tissue. Here, we map e/sQTLs, and allele-specific expression in cultured cells representing two major developmental stages, primary human neural progenitors (n = 85) and their sorted neuronal progeny (n = 74), identifying numerous loci not detected in either bulk developing cortical wall or adult cortex. Using colocalization and genetic imputation via transcriptome-wide association, we uncover cell-type-specific regulatory mechanisms underlying risk for brain-relevant traits that are active during neocortical differentiation. Specifically, we identified a progenitor-specific eQTL for CENPW co-localized with common variant associations for cortical surface area and educational attainment.

Author Correction: Genome-wide changes in lncRNA, splicing, and regional gene expression patterns in autism.

Change history: In this Letter, the labels for splicing events A3SS and A5SS were swapped in column D of Supplementary Table 3a and b. This has been corrected online.

Chromosome conformation elucidates regulatory relationships in developing human brain.

Three-dimensional physical interactions within chromosomes dynamically regulate gene expression in a tissue-specific manner. However, the 3D organization of chromosomes during human brain development and its role in regulating gene networks dysregulated in neurodevelopmental disorders, such as autism or schizophrenia, are unknown. Here we generate high-resolution 3D maps of chromatin contacts during human corticogenesis, permitting large-scale annotation of previously uncharacterized regulatory relationships relevant to the evolution of human cognition and disease. Our analyses identify hundreds of genes that physically interact with enhancers gained on the human lineage, many of which are under purifying selection and associated with human cognitive function. We integrate chromatin contacts with non-coding variants identified in schizophrenia genome-wide association studies (GWAS), highlighting multiple candidate schizophrenia risk genes and pathways, including transcription factors involved in neurogenesis, and cholinergic signalling molecules, several of which are supported by independent expression quantitative trait loci and gene expression analyses. Genome editing in human neural progenitors suggests that one of these distal schizophrenia GWAS loci regulates FOXG1 expression, supporting its potential role as a schizophrenia risk gene. This work provides a framework for understanding the effect of non-coding regulatory elements on human brain development and the evolution of cognition, and highlights novel mechanisms underlying neuropsychiatric disorders.

FoxO6 regulates memory consolidation and synaptic function.

The FoxO family of transcription factors is known to slow aging downstream from the insulin/IGF (insulin-like growth factor) signaling pathway. The most recently discovered FoxO isoform in mammals, FoxO6, is highly enriched in the adult hippocampus. However, the importance of FoxO factors in cognition is largely unknown. Here we generated mice lacking FoxO6 and found that these mice display normal learning but impaired memory consolidation in contextual fear conditioning and novel object recognition. Using stereotactic injection of viruses into the hippocampus of adult wild-type mice, we found that FoxO6 activity in the adult hippocampus is required for memory consolidation. Genome-wide approaches revealed that FoxO6 regulates a program of genes involved in synaptic function upon learning in the hippocampus. Consistently, FoxO6 deficiency results in decreased dendritic spine density in hippocampal neurons in vitro and in vivo. Thus, FoxO6 may promote memory consolidation by regulating a program coordinating neuronal connectivity in the hippocampus, which could have important implications for physiological and pathological age-dependent decline in memory.

A FOXO-Pak1 transcriptional pathway controls neuronal polarity.

Neuronal polarity is essential for normal brain development and function. However, cell-intrinsic mechanisms that govern the establishment of neuronal polarity remain to be identified. Here, we report that knockdown of endogenous FOXO proteins in hippocampal and cerebellar granule neurons, including in the rat cerebellar cortex in vivo, reveals a requirement for the FOXO transcription factors in the establishment of neuronal polarity. The FOXO transcription factors, including the brain-enriched protein FOXO6, play a critical role in axo-dendritic polarization of undifferentiated neurites, and hence in a switch from unpolarized to polarized neuronal morphology. We also identify the gene encoding the protein kinase Pak1, which acts locally in neuronal processes to induce polarity, as a critical direct target gene of the FOXO transcription factors. Knockdown of endogenous Pak1 phenocopies the effect of FOXO knockdown on neuronal polarity. Importantly, exogenous expression of Pak1 in the background of FOXO knockdown in both primary neurons and postnatal rat pups in vivo restores the polarized morphology of neurons. These findings define the FOXO proteins and Pak1 as components of a cell-intrinsic transcriptional pathway that orchestrates neuronal polarity, thus identifying a novel function for the FOXO transcription factors in a unique aspect of neural development.

A novel Hap1-Tsc1 interaction regulates neuronal mTORC1 signaling and morphogenesis in the brain.

Tuberous sclerosis complex (TSC) is a leading genetic cause of autism. The TSC proteins Tsc1 and Tsc2 control the mTORC1 signaling pathway in diverse cells, but how the mTORC1 pathway is specifically regulated in neurons remains to be elucidated. Here, using an interaction proteomics approach in neural cells including neurons, we uncover the brain-enriched protein huntingtin-associated protein 1 (Hap1) as a novel functional partner of Tsc1. Knockdown of Hap1 promotes specification of supernumerary axons in primary hippocampal neurons and profoundly impairs the positioning of pyramidal neurons in the mouse hippocampus in vivo. The Hap1 knockdown-induced phenotypes in primary neurons and in vivo recapitulate the phenotypes induced by Tsc1 knockdown. We also find that Hap1 knockdown in hippocampal neurons induces the downregulation of Tsc1 and stimulates the activity of mTORC1, as reflected by phosphorylation of the ribosomal protein S6. Inhibition of mTORC1 activity suppresses the Hap1 knockdown-induced polarity phenotype in hippocampal neurons. Collectively, these findings define a novel link between Hap1 and Tsc1 that regulates neuronal mTORC1 signaling and neuronal morphogenesis, with implications for our understanding of developmental disorders of cognition.

A conserved PTEN/FOXO pathway regulates neuronal morphology during C. elegans development.

The phosphatidylinositol 3-kinase (PI3K) signaling pathway is a conserved signal transduction cascade that is fundamental for the correct development of the nervous system. The major negative regulator of PI3K signaling is the lipid phosphatase DAF-18/PTEN, which can modulate PI3K pathway activity during neurodevelopment. Here, we identify a novel role for DAF-18 in promoting neurite outgrowth during development in Caenorhabditis elegans. We find that DAF-18 modulates the PI3K signaling pathway to activate DAF-16/FOXO and promote developmental neurite outgrowth. This activity of DAF-16 in promoting outgrowth is isoform-specific, being effected by the daf-16b isoform but not the daf-16a or daf-16d/f isoform. We also demonstrate that the capacity of DAF-16/FOXO in regulating neuron morphology is conserved in mammalian neurons. These data provide a novel mechanism by which the conserved PI3K signaling pathway regulates neuronal cell morphology during development through FOXO.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms.

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

MicroRNA-eQTLs in the developing human neocortex link miR-4707-3p expression to brain size.

Expression quantitative trait loci (eQTL) data have proven important for linking non-coding loci to protein-coding genes. But eQTL studies rarely measure microRNAs (miRNAs), small non-coding RNAs known to play a role in human brain development and neurogenesis. Here, we performed small-RNA sequencing across 212 mid-gestation human neocortical tissue samples, measured 907 expressed miRNAs, discovering 111 of which were novel, and identified 85 local-miRNA-eQTLs. Colocalization of miRNA-eQTLs with GWAS summary statistics yielded one robust colocalization of miR-4707-3p expression with educational attainment and brain size phenotypes, where the miRNA expression increasing allele was associated with decreased brain size. Exogenous expression of miR-4707-3p in primary human neural progenitor cells decreased expression of predicted targets and increased cell proliferation, indicating miR-4707-3p modulates progenitor gene regulation and cell fate decisions. Integrating miRNA-eQTLs with existing GWAS yielded evidence of a miRNA that may influence human brain size and function via modulation of neocortical brain development.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms.

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders, yet the role of cell-type-specific splicing or transcript-isoform diversity during human brain development has not been systematically investigated. Here, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 unique isoforms, of which 72.6% are novel (unannotated in Gencode-v33), and uncovered a substantial contribution of transcript-isoform diversity, regulated by RNA binding proteins, in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to re-prioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders. One-Sentence Summary: A cell-specific atlas of gene isoform expression helps shape our understanding of brain development and disease. Structured Abstract: INTRODUCTION: The development of the human brain is regulated by precise molecular and genetic mechanisms driving spatio-temporal and cell-type-specific transcript expression programs. Alternative splicing, a major mechanism increasing transcript diversity, is highly prevalent in the human brain, influences many aspects of brain development, and has strong links to neuropsychiatric disorders. Despite this, the cell-type-specific transcript-isoform diversity of the developing human brain has not been systematically investigated.RATIONALE: Understanding splicing patterns and isoform diversity across the developing neocortex has translational relevance and can elucidate genetic risk mechanisms in neurodevelopmental disorders. However, short-read sequencing, the prevalent technology for transcriptome profiling, is not well suited to capturing alternative splicing and isoform diversity. To address this, we employed third-generation long-read sequencing, which enables capture and sequencing of complete individual RNA molecules, to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution.RESULTS: We profiled microdissected GZ and CP regions of post-conception week (PCW) 15-17 human neocortex in bulk and at single-cell resolution across six subjects using high-fidelity long-read sequencing (PacBio IsoSeq). We identified 214,516 unique isoforms, of which 72.6% were novel (unannotated in Gencode), and >7,000 novel exons, expanding the proteome by 92,422 putative proteoforms. We uncovered thousands of isoform switches during cortical neurogenesis predicted to impact RNA regulatory domains or protein structure and implicating previously uncharacterized RNA-binding proteins in cellular identity and neuropsychiatric disease. At the single-cell level, early-stage excitatory neurons exhibited the greatest isoform diversity, and isoform-centric single-cell clustering led to the identification of previously uncharacterized cell states. We systematically assessed the contribution of transcriptomic features, and localized cell and spatio-temporal transcript expression signatures across neuropsychiatric disorders, revealing predominant enrichments in dynamic isoform expression and utilization patterns and that the number and complexity of isoforms per gene is strongly predictive of disease. Leveraging this resource, we re-prioritized thousands of rare de novo risk variants associated with autism spectrum disorders (ASD), intellectual disability (ID), and neurodevelopmental disorders (NDDs), more broadly, to potentially more severe consequences and revealed a larger proportion of cryptic splice variants with the expanded transcriptome annotation provided in this study.CONCLUSION: Our study offers a comprehensive landscape of isoform diversity in the human neocortex during development. This extensive cataloging of novel isoforms and splicing events sheds light on the underlying mechanisms of neurodevelopmental disorders and presents an opportunity to explore rare genetic variants linked to these conditions. The implications of our findings extend beyond fundamental neuroscience, as they provide crucial insights into the molecular basis of developmental brain disorders and pave the way for targeted therapeutic interventions. To facilitate exploration of this dataset we developed an online portal ( https://sciso.gandallab.org/ ).

Degradation of Id2 by the anaphase-promoting complex couples cell cycle exit and axonal growth.

In the developing nervous system, Id2 (inhibitor of DNA binding 2, also known as inhibitor of differentiation 2) enhances cell proliferation, promotes tumour progression and inhibits the activity of neurogenic basic helix-loop-helix (bHLH) transcription factors. The anaphase promoting complex/cyclosome and its activator Cdh1 (APC/C(Cdh1)) restrains axonal growth but the targets of APC/C(Cdh1) in neurons are unknown. Id2 and other members of the Id family are very unstable proteins that are eliminated as cells enter the quiescent state, but how they are targeted for degradation has remained elusive. Here we show that Id2 interacts with the core subunits of APC/C and Cdh1 in primary neurons. APC/C(Cdh1) targets Id2 for degradation through a destruction box motif (D box) that is conserved in Id1 and Id4. Depletion of Cdh1 stabilizes Id proteins in neurons, whereas Id2 D-box mutants are impaired for Cdh1 binding and remain stable in cells that exit from the cell cycle and contain active APC/C(Cdh1). Mutants of the Id2 D box enhance axonal growth in cerebellar granule neurons in vitro and in the context of the cerebellar cortex, and overcome the myelin inhibitory signals for growth. Conversely, activation of bHLH transcription factors induces a cluster of genes with potent axonal inhibitory functions including the gene coding for the Nogo receptor, a key transducer of myelin inhibition. Degradation of Id2 in neurons permits the accumulation of the Nogo receptor, thereby linking APC/C(Cdh1) activity with bHLH target genes for the inhibition of axonal growth. These findings indicate that deregulated Id activity might be useful to reprogramme quiescent neurons into the axonal growth mode.

Defining the evolutionary and gene regulatory logic of vertebrate neuronal diversity.

How is the exquisite cellular complexity of the vertebrate brain achieved? In this issue of Neuron, Closser et al. (2022) reveal that an expanded neuronal gene regulatory landscape may drive evolutionary cellular diversification by providing complex context and cell-specific control of effector genes.

A JIP3-regulated GSK3ß/DCX signaling pathway restricts axon branching.

Axon branching plays a critical role in establishing the accurate patterning of neuronal circuits in the brain. However, the mechanisms that control axon branching remain poorly understood. Here we report that knockdown of the brain-enriched signaling protein JNK-interacting protein 3 (JIP3) triggers exuberant axon branching and self-contact in primary granule neurons of the rat cerebellar cortex. JIP3 knockdown in cerebellar slices and in postnatal rat pups in vivo leads to the formation of ectopic branches in granule neuron parallel fiber axons in the cerebellar cortex. We also find that JIP3 restriction of axon branching is mediated by the protein kinase glycogen synthase kinase 3ß (GSK3ß). JIP3 knockdown induces the downregulation of GSK3ß in neurons, and GSK3ß knockdown phenocopies the effect of JIP3 knockdown on axon branching and self-contact. Finally, we establish doublecortin (DCX) as a novel substrate of GSK3ß in the control of axon branching and self-contact. GSK3ß phosphorylates DCX at the distinct site of Ser327 and thereby contributes to DCX function in the restriction of axon branching. Together, our data define a JIP3-regulated GSK3ß/DCX signaling pathway that restricts axon branching in the mammalian brain. These findings may have important implications for our understanding of neuronal circuitry during development, as well as the pathogenesis of neurodevelopmental disorders of cognition.

Transcriptomic Insight Into the Polygenic Mechanisms Underlying Psychiatric Disorders.

Over the past decade, large-scale genetic studies have successfully identified hundreds of genetic variants robustly associated with risk for psychiatric disorders. However, mechanistic insight and clinical translation continue to lag the pace of risk variant identification, hindered by the sheer number of targets and their predominant noncoding localization, as well as pervasive pleiotropy and incomplete penetrance. Successful next steps require identification of "causal" genetic variants and their proximal biological consequences; placing variants within biologically defined functional contexts, reflecting specific molecular pathways, cell types, circuits, and developmental windows; and characterizing the downstream, convergent neurobiological impact of polygenicity within an individual. Here, we discuss opportunities and challenges of high-throughput transcriptomic profiling in the human brain, and how transcriptomic approaches can help pinpoint mechanisms underlying genetic risk for psychiatric disorders at a scale necessary to tackle daunting levels of polygenicity. These include transcriptome-wide association studies for risk gene prioritization through integration of genome-wide association studies with expression quantitative trait loci. We outline transcriptomic results that inform our understanding of the brain-level molecular pathology of psychiatric disorders, including autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia. Finally, we discuss systems-level approaches for integration of distinct genetic, genomic, and phenotypic levels, including combining spatially resolved gene expression and human neuroimaging maps. Results highlight the importance of understanding gene expression (dys)regulation across human brain development as a major contributor to psychiatric disease pathogenesis, from common variants acting as expression quantitative trait loci to rare variants enriched for gene expression regulatory pathways.

Brain gene co-expression networks link complement signaling with convergent synaptic pathology in schizophrenia.

The most significant common variant association for schizophrenia (SCZ) reflects increased expression of the complement component 4A (C4A). Yet, it remains unclear how C4A interacts with other SCZ risk genes or whether the complement system more broadly is implicated in SCZ pathogenesis. Here, we integrate several existing, large-scale genetic and transcriptomic datasets to interrogate the functional role of the complement system and C4A in the human brain. Unexpectedly, we find no significant genetic enrichment among known complement system genes for SCZ. Conversely, brain co-expression network analyses using C4A as a seed gene reveal that genes downregulated when C4A expression increases exhibit strong and specific genetic enrichment for SCZ risk. This convergent genomic signal reflects synaptic processes, is sexually dimorphic and most prominent in frontal cortical brain regions, and is accentuated by smoking. Overall, these results indicate that synaptic pathways-rather than the complement system-are the driving force conferring SCZ risk.

Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation.

Common genetic risk for neuropsychiatric disorders is enriched in regulatory elements active during cortical neurogenesis. However, it remains poorly understood as to how these variants influence gene regulation. To model the functional impact of common genetic variation on the noncoding genome during human cortical development, we performed the assay for transposase accessible chromatin using sequencing (ATAC-seq) and analyzed chromatin accessibility quantitative trait loci (QTL) in cultured human neural progenitor cells and their differentiated neuronal progeny from 87 donors. We identified significant genetic effects on 988/1,839 neuron/progenitor regulatory elements, with highly cell-type and temporally specific effects. A subset (roughly 30%) of chromatin accessibility-QTL were also associated with changes in gene expression. Motif-disrupting alleles of transcriptional activators generally led to decreases in chromatin accessibility, whereas motif-disrupting alleles of repressors led to increases in chromatin accessibility. By integrating cell-type-specific chromatin accessibility-QTL and brain-relevant genome-wide association data, we were able to fine-map and identify regulatory mechanisms underlying noncoding neuropsychiatric disorder risk loci.

Author Correction: Transcriptomic and cellular decoding of regional brain vulnerability to neurogenetic disorders.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.