Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type–specific gene regulation in adult human brains.

Single-cell technologies measure unique cellular signatures but are typically limited to a single modality. Computational approaches allow the fusion of diverse single-cell data types, but their efficacy is difficult to validate in the absence of authentic multi-omic measurements. To comprehensively assess the molecular phenotypes of single cells, we devised single-nucleus methylcytosine, chromatin accessibility, and transcriptome sequencing (snmCAT-seq) and applied it to postmortem human frontal cortex tissue. We developed a cross-validation approach using multi-modal information to validate fine-grained cell types and assessed the effectiveness of computational data fusion methods. Correlation analysis in individual cells revealed distinct relations between methylation and gene expression. Our integrative approach enabled joint analyses of the methylome, transcriptome, chromatin accessibility, and conformation for 63 human cortical cell types. We reconstructed regulatory lineages for cortical cell populations and found specific enrichment of genetic risk for neuropsychiatric traits, enabling the prediction of cell types that are associated with diseases.

In this report, we investigated the pathophysiology of a novel hypertension-associated mitochondrial tRNA^Ala 5655A → G (m.5655A → G) mutation. The destabilization of a highly conserved base pairing (A1-U72) at the aminoacyl acceptor stem by an m.5655A → G mutation altered the tRNA^Ala function. An in vitro processing analysis showed that the m.5655A → G mutation reduced the efficiency of tRNA^Ala precursor 5′ end cleavage catalyzed by RNase P. By using cybrids constructed by transferring mitochondria from lymphoblastoid cell lines derived from a Chinese family into mitochondrial DNA (mtDNA)-less (ρ^o) cells, we showed a 41% reduction in the steady-state level of tRNA^Ala in mutant cybrids. The mutation caused an improperly aminoacylated tRNA^Ala, as suggested by aberrantly aminoacylated tRNA^Ala and slower electrophoretic mobility of mutated tRNA. A failure in tRNA^Ala metabolism contributed to variable reductions in six mtDNA-encoded polypeptides in mutant cells, ranging from 21% to 37.5%, with an average of a 29.1% reduction, compared to levels of the controls. The impaired translation caused reduced activities of mitochondrial respiration chains. Furthermore, marked decreases in the levels of mitochondrial ATP and membrane potential were observed in mutant cells. These caused increases in the production of reactive oxygen species in the mutant cybrids. The data provide evidence for the association of the tRNA^Ala 5655A → G mutation with hypertension.

The primary motor cortex (M1) is essential for voluntary fine motor control and is functionally conserved across mammals. Using high-throughput transcriptomic and epigenomic profiling of over 450,000 single nuclei in human, marmoset monkey, and mouse, we demonstrate a broadly conserved cellular makeup of this region, whose similarity mirrors evolutionary distance and is consistent between the transcriptome and epigenome. The core conserved molecular identity of neuronal and non-neuronal types allowed the generation of a cross-species consensus cell type classification and inference of conserved cell type properties across species. Despite overall conservation, many species specializations were apparent, including differences in cell type proportions, gene expression, DNA methylation, and chromatin state. Few cell type marker genes were conserved across species, providing a short list of candidate genes and regulatory mechanisms responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allowed the Patch-seq identification of layer 5 (L5) corticospinal Betz cells in non-human primate and human and characterization of their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell type diversity in M1 across mammals and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

The pathophysiology underlying spiral ganglion cell defect–induced deafness remains elusive. Using the whole exome sequencing approach, in combination with functional assays and a mouse disease model, we identified the potentially novel deafness-causative MAP1B gene encoding a highly conserved microtubule-associated protein. Three novel heterozygous MAP1B mutations (c.4198A>G, p.1400S>G; c.2768T>C, p.923I>T; c.5512T>C, p.1838F>L) were cosegregated with autosomal dominant inheritance of nonsyndromic sensorineural hearing loss in 3 unrelated Chinese families. Here, we show that MAP1B is highly expressed in the spiral ganglion neurons in the mouse cochlea. Using otic sensory neuron–like cells, generated by pluripotent stem cells from patients carrying the MAP1B mutation and control subject, we demonstrated that the p.1400S>G mutation caused the reduced levels and deficient phosphorylation of MAP1B, which are involved in the microtubule stability and dynamics. Strikingly, otic sensory neuron–like cells exhibited disturbed dynamics of microtubules, axonal elongation, and defects in electrophysiological properties. Dysfunctions of these derived otic sensory neuron–like cells were rescued by genetically correcting MAP1B mutation using CRISPR/Cas9 technology. Involvement of MAP1B in hearing was confirmed by audiometric evaluation of Map1b heterozygous KO mice. These mutant mice displayed late-onset progressive sensorineural hearing loss that was more pronounced in the high frequencies. The spiral ganglion neurons isolated from Map1b mutant mice exhibited the deficient phosphorylation and disturbed dynamics of microtubules. Map1b deficiency yielded defects in the morphology and electrophysiology of spiral ganglion neurons, but it did not affect the morphologies of cochlea in mice. Therefore, our data demonstrate that dysfunctions of spiral ganglion neurons induced by MAP1B deficiency caused hearing loss.

Single cell transcriptomics has transformed the characterization of brain cell identity by providing quantitative molecular signatures for large, unbiased samples of brain cell populations. With the proliferation of taxonomies based on individual datasets, a major challenge is to integrate and validate results toward defining biologically meaningful cell types. We used a battery of single-cell transcriptome and epigenome measurements generated by the BRAIN Initiative Cell Census Network (BICCN) to comprehensively assess the molecular signatures of cell types in the mouse primary motor cortex (MOp). We further developed computational and statistical methods to integrate these multimodal data and quantitatively validate the reproducibility of the cell types. The reference atlas, based on more than 600,000 high quality single-cell or -nucleus samples assayed by six molecular modalities, is a comprehensive molecular account of the diverse neuronal and non-neuronal cell types in MOp. Collectively, our study indicates that the mouse primary motor cortex contains over 55 neuronal cell types that are highly replicable across analysis methods, sequencing technologies, and modalities. We find many concordant multimodal markers for each cell type, as well as thousands of genes and gene regulatory elements with discrepant transcriptomic and epigenomic signatures. These data highlight the complex molecular regulation of brain cell types and will directly enable design of reagents to target specific MOp cell types for functional analysis.