Single-cell sequencing technologies, including transcriptomic and epigenomic assays, are transforming our understanding of the cellular building blocks of neural circuits. By directly measuring multiple molecular signatures in thousands to millions of individual cells, single-cell sequencing methods can comprehensively characterize the diversity of brain cell types. These measurements uncover gene regulatory mechanisms that shape cellular identity and provide insight into developmental and evolutionary relationships between brain cell populations. Single-cell sequencing data can aid the design of tools for targeted functional studies of brain circuit components, linking molecular signatures with anatomy, connectivity, morphology, and physiology. Here, we discuss the fundamental principles of single-cell transcriptome and epigenome sequencing, integrative computational analysis of the data, and key applications in neuroscience.

Posttranscriptional adenosine-to-inosine modifications amplify the functionality of RNA molecules in the brain, yet the cellular and genetic regulation of RNA editing is poorly described. We quantify base-specific RNA editing across three major cell populations from the human prefrontal cortex: glutamatergic neurons, medial ganglionic eminence-derived GABAergic neurons, and oligodendrocytes. We identify more selective editing and hyper-editing in neurons relative to oligodendrocytes. RNA editing patterns are highly cell type-specific, with 189,229 cell type-associated sites. The cellular specificity for thousands of sites is confirmed by single nucleus RNA-sequencing. Importantly, cell type-associated sites are enriched in GTEx RNA-sequencing data, edited ~twentyfold higher than all other sites, and variation in RNA editing is largely explained by neuronal proportions in bulk brain tissue. Finally, we uncover 661,791 cis-editing quantitative trait loci across thirteen brain regions, including hundreds with cell type-associated features. These data reveal an expansive repertoire of highly regulated RNA editing sites across human brain cell types and provide a resolved atlas linking cell types to editing variation and genetic regulatory effects.

Brain function depends on interaction of diverse cell types whose gene expression and identity are defined, in part, by epigenetic mechanisms. Neuronal DNA contains two major epigenetic modifications, methylcytosine (mC) and hydroxymethylcytosine (hmC), yet their cell type-specific landscapes and relationship with gene expression are poorly understood. We report high-resolution (h)mC analyses, together with transcriptome and histone modification profiling, in three major cell types in human prefrontal cortex: glutamatergic excitatory neurons, medial ganglionic eminence-derived γ-aminobutyric acid (GABA)ergic inhibitory neurons, and oligodendrocytes. We detected a unique association between hmC and gene expression in inhibitory neurons that differed significantly from the pattern in excitatory neurons and oligodendrocytes. We also found that risk loci associated with neuropsychiatric diseases were enriched near regions of reduced hmC in excitatory neurons and reduced mC in inhibitory neurons. Our findings indicate differential roles for mC and hmC in regulation of gene expression in different brain cell types, with implications for the etiology of human brain diseases.

A repeat expansion in the C9orf72 (C9) gene is the most common genetic cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). Here we investigate single nucleus transcriptomics (snRNA-seq) and epigenomics (snATAC-seq) in postmortem motor and frontal cortices from C9-ALS, C9-FTD, and control donors. C9-ALS donors present pervasive alterations of gene expression with concordant changes in chromatin accessibility and histone modifications. The greatest alterations occur in upper and deep layer excitatory neurons, as well as in astrocytes. In neurons, the changes imply an increase in proteostasis, metabolism, and protein expression pathways, alongside a decrease in neuronal function. In astrocytes, the alterations suggest activation and structural remodeling. Conversely, C9-FTD donors have fewer high-quality neuronal nuclei in the frontal cortex and numerous gene expression changes in glial cells. These findings highlight a context-dependent molecular disruption in C9-ALS and C9-FTD, indicating unique effects across cell types, brain regions, and diseases.

Early life experience influences stress reactivity and mental health through effects on cognitive-emotional functions that are, in part, linked to gene expression in the dorsal and ventral hippocampus. The hippocampal dentate gyrus (DG) is a major site for experience-dependent plasticity associated with sustained transcriptional alterations, potentially mediated by epigenetic modifications. Here, we report comprehensive DNA methylome, hydroxymethylome and transcriptome data sets from mouse dorsal and ventral DG. We find genome-wide transcriptional and methylation differences between dorsal and ventral DG, including at key developmental transcriptional factors. Peripubertal environmental enrichment increases hippocampal volume and enhances dorsal DG-specific differences in gene expression. Enrichment also enhances dorsal-ventral differences in DNA methylation, including at binding sites of the transcription factor NeuroD1, a regulator of adult neurogenesis. These results indicate a dorsal-ventral asymmetry in transcription and methylation that parallels well-known functional and anatomical differences, and that may be enhanced by environmental enrichment.

The mammalian brain contains diverse neuronal types, yet we lack single-cell epigenomic assays that are able to identify and characterize them. DNA methylation is a stable epigenetic mark that distinguishes cell types and marks regulatory elements. We generated >6000 methylomes from single neuronal nuclei and used them to identify 16 mouse and 21 human neuronal subpopulations in the frontal cortex. CG and non-CG methylation exhibited cell type-specific distributions, and we identified regulatory elements with differential methylation across neuron types. Methylation signatures identified a layer 6 excitatory neuron subtype and a unique human parvalbumin-expressing inhibitory neuron subtype. We observed stronger cross-species conservation of regulatory elements in inhibitory neurons than in excitatory neurons. Single-nucleus methylomes expand the atlas of brain cell types and identify regulatory elements that drive conserved brain cell diversity.

Altered transcriptional and epigenetic regulation of brain cell types may contribute to cognitive changes with advanced age. Using single-nucleus multi-omic DNA methylation and transcriptome sequencing (snmCT-seq) in frontal cortex from young adult and aged donors, we found widespread age- and sex-related variation in specific neuron types. The proportion of inhibitory SST- and VIP-expressing neurons was reduced in aged donors. Excitatory neurons had more profound age-related changes in their gene expression and DNA methylation than inhibitory cells. Hundreds of genes involved in synaptic activity, including EGR1, were less expressed in aged adults. Genes located in subtelomeric regions increased their expression with age and correlated with reduced telomere length. We further mapped cell-type-specific sex differences in gene expression and X-inactivation escape genes. Multi-omic single-nucleus epigenomes and transcriptomes provide new insight into the effects of age and sex on human neurons.

Two epigenetic pathways of transcriptional repression, DNA methylation and polycomb repressive complex 2 (PRC2), are known to regulate neuronal development and function. However, their respective contributions to brain maturation are unknown. We found that conditional loss of the de novo DNA methyltransferase Dnmt3a in mouse excitatory neurons altered expression of synapse-related genes, stunted synapse maturation, and impaired working memory and social interest. At the genomic level, loss of Dnmt3a abolished postnatal accumulation of CG and non-CG DNA methylation, leaving adult neurons with an unmethylated, fetal-like epigenomic pattern at ~222,000 genomic regions. The PRC2-associated histone modification, H3K27me3, increased at many of these sites. Our data support a dynamic interaction between two fundamental modes of epigenetic repression during postnatal maturation of excitatory neurons, which together confer robustness on neuronal regulation.

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.