The inability of neurons to regenerate long axons within the CNS is a major impediment to improving outcome after spinal cord injury, stroke, and other CNS insults. Recent advances have uncovered an intrinsic program that involves coordinate regulation by multiple transcription factors that can be manipulated to enhance growth in the peripheral nervous system. Here, we use a systems genomics approach to characterize regulatory relationships of regeneration-associated transcription factors, identifying RE1-Silencing Transcription Factor (REST; Neuron-Restrictive Silencer Factor, NRSF) as a predicted upstream suppressor of a pro-regenerative gene program associated with axon regeneration in the CNS. We validate our predictions using multiple paradigms, showing that mature mice bearing cell type-specific deletions of REST or expressing dominant-negative mutant REST show improved regeneration of the corticospinal tract and optic nerve after spinal cord injury and optic nerve crush, which is accompanied by upregulation of regeneration-associated genes in cortical motor neurons and retinal ganglion cells, respectively. These analyses identify a role for REST as an upstream suppressor of the intrinsic regenerative program in the CNS and demonstrate the utility of a systems biology approach involving integrative genomics and bio-informatics to prioritize hypotheses relevant to CNS repair.

Regulatory programs governing neuronal death and axon regeneration in neurodegenerative diseases remain poorly understood. In adult mice, optic nerve crush (ONC) injury by severing retinal ganglion cell (RGC) axons results in massive RGC death and regenerative failure. We performed an in vivo CRISPR-Cas9-based genome-wide screen of 1,893 transcription factors (TFs) to seek repressors of RGC survival and axon regeneration following ONC. In parallel, we profiled the epigenetic and transcriptional landscapes of injured RGCs by ATAC-seq and RNA-seq to identify injury-responsive TFs and their targets. These analyses converged on four TFs as critical survival regulators, of which ATF3/CHOP preferentially regulate pathways activated by cytokines and innate immunity and ATF4/C/EBPγ regulate pathways engaged by intrinsic neuronal stressors. Manipulation of these TFs protects RGCs in a glaucoma model. Our results reveal core transcription programs that transform an initial axonal insult into a degenerative process and suggest novel strategies for treating neurodegenerative diseases.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

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