Single-nucleus transcriptome analysis reveals cell-type-specific molecular signatures across reward circuitry in the human brain

Genetic Implication of Specific Glutamatergic Neurons of the Prefrontal Cortex in the Pathophysiology of Schizophrenia

Background

The prefrontal cortex (PFC) has been strongly implicated in the pathophysiology of schizophrenia. Here, we combined high-resolution single-nuclei RNA sequencing data from the human PFC with large-scale genomic data for schizophrenia to identify constituent cell populations likely to mediate genetic liability to the disorder.

Methods

Gene expression specificity values were calculated from a single-nuclei RNA sequencing dataset comprising 84 cell populations from the human PFC, spanning gestation to adulthood. Enrichment of schizophrenia common variant liability and burden of rare protein-truncating coding variants were tested in genes with high expression specificity for each cell type. We also explored schizophrenia common variant associations in relation to gene expression across the developmental trajectory of implicated neurons.

Results

Common risk variation for schizophrenia was prominently enriched in genes with high expression specificity for a population of mature layer 4 glutamatergic neurons emerging in infancy. Common variant liability to schizophrenia increased along the developmental trajectory of this neuronal population. Fine-mapped genes at schizophrenia genome-wide association study risk loci had significantly higher expression specificity than other genes in these neurons and in a population of layer 5/6 glutamatergic neurons. People with schizophrenia had a higher rate of rare protein-truncating coding variants in genes expressed by cells of the PFC than control individuals, but no cell population was significantly enriched above this background rate.

Conclusions

We identified a population of layer 4 glutamatergic PFC neurons likely to be particularly affected by common variant genetic risk for schizophrenia, which may contribute to disturbances in thalamocortical connectivity in the condition.

NPAS4 supports cocaine-conditioned cues in rodents by controlling the cell type-specific activation balance in the nucleus accumbens

Powerful associations that link drugs of abuse with cues in the drug-paired environment often serve as prepotent relapse triggers. Drug-associated contexts and cues activate ensembles of nucleus accumbens (NAc) neurons, including D1-class medium spiny neurons (MSNs) that typically promote, and D2-class MSNs that typically oppose, drug seeking. We found that in mice, cocaine conditioning upregulated transiently the activity-regulated transcription factor, Neuronal PAS Domain Protein 4 (NPAS4), in a small subset of NAc neurons. The NPAS4+ NAc ensemble was required for cocaine conditioned place preference. We also observed that NPAS4 functions within NAc D2-, but not D1-, MSNs to support cocaine-context associations and cue-induced cocaine, but not sucrose, seeking. Together, our data show that the NPAS4+ ensemble of NAc neurons is essential for cocaine-context associations in mice, and that NPAS4 itself functions in NAc D2-MSNs to support cocaine-context associations by suppressing drug-induced counteradaptations that oppose relapse-related behaviour.

Multi-omic analysis of Huntington's disease reveals a compensatory astrocyte state

The mechanisms underlying the selective regional vulnerability to neurodegeneration in Huntington's disease (HD) have not been fully defined. To explore the role of astrocytes in this phenomenon, we used single-nucleus and bulk RNAseq, lipidomics, HTT gene CAG repeat-length measurements, and multiplexed immunofluorescence on HD and control post-mortem brains. We identified genes that correlated with CAG repeat length, which were enriched in astrocyte genes, and lipidomic signatures that implicated poly-unsaturated fatty acids in sensitizing neurons to cell death. Because astrocytes play essential roles in lipid metabolism, we explored the heterogeneity of astrocytic states in both protoplasmic and fibrous-like (CD44+) astrocytes. Significantly, one protoplasmic astrocyte state showed high levels of metallothioneins and was correlated with the selective vulnerability of distinct striatal neuronal populations. When modeled in vitro, this state improved the viability of HD-patient-derived spiny projection neurons. Our findings uncover key roles of astrocytic states in protecting against neurodegeneration in HD.

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.

A brain cell atlas integrating single-cell transcriptomes across human brain regions

While single-cell technologies have greatly advanced our comprehension of human brain cell types and functions, studies including large numbers of donors and multiple brain regions are needed to extend our understanding of brain cell heterogeneity. Integrating atlas-level single-cell data presents a chance to reveal rare cell types and cellular heterogeneity across brain regions. Here we present the Brain Cell Atlas, a comprehensive reference atlas of brain cells, by assembling single-cell data from 70 human and 103 mouse studies of the brain throughout major developmental stages across brain regions, covering over 26.3 million cells or nuclei from both healthy and diseased tissues. Using machine-learning based algorithms, the Brain Cell Atlas provides a consensus cell type annotation, and it showcases the identification of putative neural progenitor cells and a cell subpopulation of PCDH9high microglia in the human brain. We demonstrate the gene regulatory difference of PCDH9high microglia between hippocampus and prefrontal cortex and elucidate the cell-cell communication network. The Brain Cell Atlas presents an atlas-level integrative resource for comparing brain cells in different environments and conditions within the Human Cell Atlas.

Interneuron diversity in the human dorsal striatum

Deciphering the striatal interneuron diversity is key to understanding the basal ganglia circuit and to untangling the complex neurological and psychiatric diseases affecting this brain structure. We performed snRNA-seq and spatial transcriptomics of postmortem human caudate nucleus and putamen samples to elucidate the diversity and abundance of interneuron populations and their inherent transcriptional structure in the human dorsal striatum. We propose a comprehensive taxonomy of striatal interneurons with eight main classes and fourteen subclasses, providing their full transcriptomic identity and spatial expression profile as well as additional quantitative FISH validation for specific populations. We have also delineated the correspondence of our taxonomy with previous standardized classifications and shown the main transcriptomic and class abundance differences between caudate nucleus and putamen. Notably, based on key functional genes such as ion channels and synaptic receptors, we found matching known mouse interneuron populations for the most abundant populations, the recently described PTHLH and TAC3 interneurons. Finally, we were able to integrate other published datasets with ours, supporting the generalizability of this harmonized taxonomy.

Comparative basolateral amygdala connectomics reveals dissociable single-neuron projection patterns to frontal cortex in macaques and mice

Basolateral amygdala (BLA) is a key hub for affect in the brain,1,2,3 and dysfunction within this area contributes to a host of psychiatric disorders.4,5 BLA is extensively and reciprocally interconnected with frontal cortex,6,7,8,9 and some aspects of its function are evolutionarily conserved across rodents, anthropoid primates, and humans.10 Neuron density in BLA is substantially lower in primates compared to murine rodents,11 and frontal cortex (FC) is dramatically expanded in primates, particularly the more anterior granular and dysgranular areas.12,13,14 Yet, how these anatomical differences influence the projection patterns of single BLA neurons to frontal cortex across rodents and primates is unknown. Using a barcoded connectomic approach, we assessed the single BLA neuron connections to frontal cortex in mice and macaques. We found that BLA neurons are more likely to project to multiple distinct parts of FC in mice than in macaques. Further, while single BLA neuron projections to nucleus accumbens were similarly organized in mice and macaques, BLA-FC connections differed substantially. Notably, BLA connections to subcallosal anterior cingulate cortex (scACC) in macaques were least likely to branch to other medial frontal cortex areas compared to perigenual ACC (pgACC). This pattern of connections was reversed in the mouse homologues of these areas, infralimbic and prelimbic cortex (IL and PL), mirroring functional differences between rodents and non-human primates. Taken together, these results indicate that BLA connections to FC are not linearly scaled from mice to macaques and instead the organization of single-neuron BLA connections is distinct between these species.

Transcriptomic analysis of the human habenula in schizophrenia

Pathophysiology of many neuropsychiatric disorders, including schizophrenia (SCZD), is linked to habenula (Hb) function. While pharmacotherapies and deep brain stimulation targeting the Hb are emerging as promising therapeutic treatments, little is known about the cell type-specific transcriptomic organization of the human Hb or how it is altered in SCZD. Here we define the molecular neuroanatomy of the human Hb and identify transcriptomic changes in individuals with SCZD compared to neurotypical controls. Utilizing Hb-enriched postmortem human brain tissue, we performed single nucleus RNA-sequencing (snRNA-seq; n=7 neurotypical donors) and identified 17 molecularly defined Hb cell types across 16,437 nuclei, including 3 medial and 7 lateral Hb populations, several of which were conserved between rodents and humans. Single molecule fluorescent in situ hybridization (smFISH; n=3 neurotypical donors) validated snRNA-seq Hb cell types and mapped their spatial locations. Bulk RNA-sequencing and cell type deconvolution in Hb-enriched tissue from 35 individuals with SCZD and 33 neurotypical controls yielded 45 SCZD-associated differentially expressed genes (DEGs, FDR < 0.05), with 32 (71%) unique to Hb-enriched tissue. eQTL analysis identified 717 independent SNP-gene pairs (FDR < 0.05), where either the SNP is a SCZD risk variant (16 pairs) or the gene is a SCZD DEG (7 pairs). eQTL and SCZD risk colocalization analysis identified 16 colocalized genes. These results identify topographically organized cell types with distinct molecular signatures in the human Hb and demonstrate unique genetic changes associated with SCZD, thereby providing novel molecular insights into the role of Hb in neuropsychiatric disorders.

One Sentence Summary

Transcriptomic analysis of the human habenula and identification of molecular changes associated with schizophrenia risk and illness state.

Nucleus Accumbens Corticotropin-Releasing Hormone Neurons Projecting to the Bed Nucleus of the Stria Terminalis Promote Wakefulness and Positive Affective State

The nucleus accumbens (NAc) plays an important role in various emotional and motivational behaviors that rely on heightened wakefulness. However, the neural mechanisms underlying the relationship between arousal and emotion regulation in NAc remain unclear. Here, we investigated the roles of a specific subset of inhibitory corticotropin-releasing hormone neurons in the NAc (NAcCRH) in regulating arousal and emotional behaviors in mice. We found an increased activity of NAcCRH neurons during wakefulness and rewarding stimulation. Activation of NAcCRH neurons converts NREM or REM sleep to wakefulness, while inhibition of these neurons attenuates wakefulness. Remarkably, activation of NAcCRH neurons induces a place preference response (PPR) and decreased basal anxiety level, whereas their inactivation induces a place aversion response and anxious state. NAcCRH neurons are identified as the major NAc projection neurons to the bed nucleus of the stria terminalis (BNST). Furthermore, activation of the NAcCRH-BNST pathway similarly induced wakefulness and positive emotional behaviors. Taken together, we identified a basal forebrain CRH pathway that promotes the arousal associated with positive affective states.

Bulk and Single-cell Transcriptomic Brain Data Identify Overlapping Processes and Cell-types with Human AUD and Mammalian Models of Alcohol Use

This study explores the neurobiological underpinnings of alcohol use disorder (AUD) by integrating bulk and single-cell transcriptomic data from humans, primates, and mice across three brain regions associated with addiction (i.e., prefrontal cortex (PFC), nucleus accumbens (NAc), and central amygdala (CeA)). We compared AUD RNA expression and cell-type abundance from 92 human brain to data from 53 primates and 90 mice engaged in diverse alcohol use paradigms. The findings revealed significant and reproducible correlations between human AUD and mammalian models of alcohol use that vary by tissue, species, and behavioral paradigm. The strongest correlations occurred between primate and mouse models of binge drinking (i.e., high drinking in the dark). Certain primate models demonstrated that the brain RNA correlations with human alcohol use disorder (AUD) were approximately 40% as strong as the correlations observed within human samples themselves. By integrating single-cell transcriptomic data, this study observed decreased oligodendrocyte proportions in the PFC and NAc of human AUD with similar trends in animal models. Gene co-expression network analyses revealed conserved systems associated with human AUD and animal models of heavy/binge alcohol consumption. Gene co-expression networks were enriched for pathways related to inflammation, myelination, and synaptic plasticity and the genes within them accounted for ~20% of the heritability in human alcohol consumption. Identified hub genes were associated with relevant traits (e.g., impulsivity, motivation) in humans and mice. This study sheds light on conserved biological entities underlying AUD and chronic alcohol use, providing insights into the cellular, genetic, and neuromolecular basis across species.

Microglial function interacts with the environment to affect sex-specific depression risk

There is a two-fold higher incidence of depression in females compared to men with recent studies suggesting a role for microglia in conferring this sex-dependent depression risk. In this study we investigated the nature of this relation. Using GWAS enrichment, gene-set enrichment analysis and Mendelian randomization, we found minimal evidence for a direct relation between genes functionally related to microglia and sex-dependent genetic risk for depression. We then used expression quantitative trait loci and single nucleus RNA-sequencing resources to generate polygenic scores (PGS) representative of individual variation in microglial function in the adult (UK Biobank; N = 54753-72682) and fetal (ALSPAC; N = 1452) periods. The adult microglial PGS moderated the association between BMI (UK Biobank; beta = 0.001, 95 %CI 0.0009 to 0.003, P = 7.74E-6) and financial insecurity (UK Biobank; beta = 0.001, 95 %CI 0.005 to 0.015, P = 2E-4) with depressive symptoms in females. The fetal microglia PGS moderated the association between maternal prenatal depressive symptoms and offspring depressive symptoms at 24 years in females (ALSPAC; beta = 0.04, 95 %CI 0.004 to 0.07, P = 0.03). We found no evidence for an interaction between the microglial PGS and depression risk factors in males. Our results illustrate a role for microglial function in the conferral of sex-dependent depression risk following exposure to a depression risk factor.

Quantification and statistical modeling of droplet-based single-nucleus RNA-sequencing data

In complex tissues containing cells that are difficult to dissociate, single-nucleus RNA-sequencing (snRNA-seq) has become the preferred experimental technology over single-cell RNA-sequencing (scRNA-seq) to measure gene expression. To accurately model these data in downstream analyses, previous work has shown that droplet-based scRNA-seq data are not zero-inflated, but whether droplet-based snRNA-seq data follow the same probability distributions has not been systematically evaluated. Using pseudonegative control data from nuclei in mouse cortex sequenced with the 10x Genomics Chromium system and mouse kidney sequenced with the DropSeq system, we found that droplet-based snRNA-seq data follow a negative binomial distribution, suggesting that parametric statistical models applied to scRNA-seq are transferable to snRNA-seq. Furthermore, we found that the quantification choices in adapting quantification mapping strategies from scRNA-seq to snRNA-seq can play a significant role in downstream analyses and biological interpretation. In particular, reference transcriptomes that do not include intronic regions result in significantly smaller library sizes and incongruous cell type classifications. We also confirmed the presence of a gene length bias in snRNA-seq data, which we show is present in both exonic and intronic reads, and investigate potential causes for the bias.

Deciphering brain cellular and behavioral mechanisms: Insights from single-cell and spatial RNA sequencing

The brain is a complex computing system composed of a multitude of interacting neurons. The computational outputs of this system determine the behavior and perception of every individual. Each brain cell expresses thousands of genes that dictate the cell's function and physiological properties. Therefore, deciphering the molecular expression of each cell is of great significance for understanding its characteristics and role in brain function. Additionally, the positional information of each cell can provide crucial insights into their involvement in local brain circuits. In this review, we briefly overview the principles of single-cell RNA sequencing and spatial transcriptomics, the potential issues and challenges in their data processing, and their applications in brain research. We further outline several promising directions in neuroscience that could be integrated with single-cell RNA sequencing, including neurodevelopment, the identification of novel brain microstructures, cognition and behavior, neuronal cell positioning, molecules and cells related to advanced brain functions, sleep-wake cycles/circadian rhythms, and computational modeling of brain function. We believe that the deep integration of these directions with single-cell and spatial RNA sequencing can contribute significantly to understanding the roles of individual cells or cell types in these specific functions, thereby making important contributions to addressing critical questions in those fields. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA in Disease and Development > RNA in Development RNA in Disease and Development > RNA in Disease.

Single-cell transcriptomic reveals a cell atlas and diversity of chicken amygdala responded to social hierarchy

Amygdala serves as a highly cellular, heterogeneous brain region containing excitatory and inhibitory neurons and is involved in the dopamine and serotoninergic neuron systems. An increasing number of studies have revealed the underpinned mechanism mediating social hierarchy in mammal and vertebrate, however, there are rare studies conducted on how amygdala on social hierarchy in poultry. In this study, we conducted food competition tests and determined the social hierarchy of the rooster. We performed cross-species analysis with mammalian amygdala, and found that cell types of human and rhesus monkeys were more closely related and that of chickens were more distant. We identified 26 clusters and divided them into 10 main clusters, of which GABAergic and glutamatergic neurons were associated with social behaviors. In conclusion, our results provide to serve the developmental studies of the amygdala neuron system and new insights into the underpinned mechanism of social hierarchy in roosters.

The Advancement and Application of the Single-Cell Transcriptome in Biological and Medical Research

Single-cell RNA sequencing technology (scRNA-seq) has been steadily developing since its inception in 2009. Unlike bulk RNA-seq, scRNA-seq identifies the heterogeneity of tissue cells and reveals gene expression changes in individual cells at the microscopic level. Here, we review the development of scRNA-seq, which has gone through iterations of reverse transcription, in vitro transcription, smart-seq, drop-seq, 10 × Genomics, and spatial single-cell transcriptome technologies. The technology of 10 × Genomics has been widely applied in medicine and biology, producing rich research results. Furthermore, this review presents a summary of the analytical process for single-cell transcriptome data and its integration with other omics analyses, including genomes, epigenomes, proteomes, and metabolomics. The single-cell transcriptome has a wide range of applications in biology and medicine. This review analyzes the applications of scRNA-seq in cancer, stem cell research, developmental biology, microbiology, and other fields. In essence, scRNA-seq provides a means of elucidating gene expression patterns in single cells, thereby offering a valuable tool for scientific research. Nevertheless, the current single-cell transcriptome technology is still imperfect, and this review identifies its shortcomings and anticipates future developments. The objective of this review is to facilitate a deeper comprehension of scRNA-seq technology and its applications in biological and medical research, as well as to identify avenues for its future development in alignment with practical needs.

Reelin marks cocaine-activated striatal ensembles, promotes neuronal excitability, and regulates cocaine reward

Drugs of abuse activate defined neuronal ensembles in brain reward structures such as the nucleus accumbens (NAc), which are thought to promote the enduring synaptic, circuit, and behavioral consequences of drug exposure. While the molecular and cellular effects arising from experience with drugs like cocaine are increasingly well understood, the mechanisms that sculpt NAc ensemble participation are largely unknown. Here, we leveraged unbiased single-nucleus transcriptional profiling to identify expression of the secreted glycoprotein Reelin (encoded by the Reln gene) as a marker of cocaine-activated neuronal ensembles within the rat NAc. Multiplexed in situ detection confirmed selective expression of the immediate early gene Fos in Reln+ neurons after cocaine experience, and also revealed enrichment of Reln mRNA in Drd1 + medium spiny neurons (MSNs) in both the rat and human brain. Using a novel CRISPR interference strategy enabling selective Reln knockdown in the adult NAc, we observed altered expression of genes linked to calcium signaling, emergence of a transcriptional trajectory consistent with loss of cocaine sensitivity, and a striking decrease in MSN intrinsic excitability. At the behavioral level, loss of Reln prevented cocaine locomotor sensitization, abolished cocaine place preference memory, and decreased cocaine self-administration behavior. Together, these results identify Reelin as a critical mechanistic link between ensemble participation and cocaine-induced behavioral adaptations.

The single-cell opioid responses in the context of HIV (SCORCH) consortium

Substance use disorders (SUD) and drug addiction are major threats to public health, impacting not only the millions of individuals struggling with SUD, but also surrounding families and communities. One of the seminal challenges in treating and studying addiction in human populations is the high prevalence of co-morbid conditions, including an increased risk of contracting a human immunodeficiency virus (HIV) infection. Of the ~15 million people who inject drugs globally, 17% are persons with HIV. Conversely, HIV is a risk factor for SUD because chronic pain syndromes, often encountered in persons with HIV, can lead to an increased use of opioid pain medications that in turn can increase the risk for opioid addiction. We hypothesize that SUD and HIV exert shared effects on brain cell types, including adaptations related to neuroplasticity, neurodegeneration, and neuroinflammation. Basic research is needed to refine our understanding of these affected cell types and adaptations. Studying the effects of SUD in the context of HIV at the single-cell level represents a compelling strategy to understand the reciprocal interactions among both conditions, made feasible by the availability of large, extensively-phenotyped human brain tissue collections that have been amassed by the Neuro-HIV research community. In addition, sophisticated animal models that have been developed for both conditions provide a means to precisely evaluate specific exposures and stages of disease. We propose that single-cell genomics is a uniquely powerful technology to characterize the effects of SUD and HIV in the brain, integrating data from human cohorts and animal models. We have formed the Single-Cell Opioid Responses in the Context of HIV (SCORCH) consortium to carry out this strategy.

Genetic Implication of Prenatal GABAergic and Cholinergic Neuron Development in Susceptibility to Schizophrenia

BACKGROUND

The ganglionic eminences (GE) are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine-releasing neurons of the forebrain. Given the evidence for GABAergic, cholinergic, and neurodevelopmental disturbances in schizophrenia, we tested the potential involvement of GE neuron development in mediating genetic risk for the condition.

STUDY DESIGN

We combined data from a recent large-scale genome-wide association study of schizophrenia with single-cell RNA sequencing data from the human GE to test the enrichment of schizophrenia risk variation in genes with high expression specificity for developing GE cell populations. We additionally performed the single nuclei Assay for Transposase-Accessible Chromatin with Sequencing (snATAC-Seq) to map potential regulatory genomic regions operating in individual cell populations of the human GE, using these to test for enrichment of schizophrenia common genetic variant liability and to functionally annotate non-coding variants-associated with the disorder.

STUDY RESULTS

Schizophrenia common variant liability was enriched in genes with high expression specificity for developing neuron populations that are predicted to form dopamine D1 and D2 receptor-expressing GABAergic medium spiny neurons of the striatum, cortical somatostatin-positive GABAergic interneurons, calretinin-positive GABAergic neurons, and cholinergic neurons. Consistent with these findings, schizophrenia genetic risk was concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the GE.

CONCLUSIONS

Our study implicates prenatal development of specific populations of GABAergic and cholinergic neurons in later susceptibility to schizophrenia, and provides a map of predicted regulatory genomic elements operating in cells of the GE.

Circuit-Wide Gene Network Analysis Reveals Sex-Specific Roles for Phosphodiesterase 1b in Cocaine Addiction

Cocaine use disorder is a significant public health issue without an effective pharmacological treatment. Successful treatments are hindered in part by an incomplete understanding of the molecular mechanisms that underlie long-lasting maladaptive plasticity and addiction-like behaviors. Here, we leverage a large RNA sequencing dataset to generate gene coexpression networks across six interconnected regions of the brain's reward circuitry from mice that underwent saline or cocaine self-administration. We identify phosphodiesterase 1b (Pde1b), a Ca2+/calmodulin-dependent enzyme that increases cAMP and cGMP hydrolysis, as a central hub gene within a nucleus accumbens (NAc) gene module that was bioinformatically associated with addiction-like behavior. Chronic cocaine exposure increases Pde1b expression in NAc D2 medium spiny neurons (MSNs) in male but not female mice. Viral-mediated Pde1b overexpression in NAc reduces cocaine self-administration in female rats but increases seeking in both sexes. In female mice, overexpressing Pde1b in D1 MSNs attenuates the locomotor response to cocaine, with the opposite effect in D2 MSNs. Overexpressing Pde1b in D1/D2 MSNs had no effect on the locomotor response to cocaine in male mice. At the electrophysiological level, Pde1b overexpression reduces sEPSC frequency in D1 MSNs and regulates the excitability of NAc MSNs. Lastly, Pde1b overexpression significantly reduced the number of differentially expressed genes (DEGs) in NAc following chronic cocaine, with discordant effects on gene transcription between sexes. Together, we identify novel gene modules across the brain's reward circuitry associated with addiction-like behavior and explore the role of Pde1b in regulating the molecular, cellular, and behavioral responses to cocaine.

Analysis of gene expression in the postmortem brain of neurotypical Black Americans reveals contributions of genetic ancestry

Ancestral differences in genomic variation affect the regulation of gene expression; however, most gene expression studies have been limited to European ancestry samples or adjusted to identify ancestry-independent associations. Here, we instead examined the impact of genetic ancestry on gene expression and DNA methylation in the postmortem brain tissue of admixed Black American neurotypical individuals to identify ancestry-dependent and ancestry-independent contributions. Ancestry-associated differentially expressed genes (DEGs), transcripts and gene networks, while notably not implicating neurons, are enriched for genes related to the immune response and vascular tissue and explain up to 26% of heritability for ischemic stroke, 27% of heritability for Parkinson disease and 30% of heritability for Alzheimer's disease. Ancestry-associated DEGs also show general enrichment for the heritability of diverse immune-related traits but depletion for psychiatric-related traits. We also compared Black and non-Hispanic white Americans, confirming most ancestry-associated DEGs. Our results delineate the extent to which genetic ancestry affects differences in gene expression in the human brain and the implications for brain illness risk.

Chst9 marks a spatially and transcriptionally unique population of Oprm1-expressing neurons in the nucleus accumbens

Opioids produce addictive, analgesic, and euphoric effects via actions at mu opioid receptors (μORs). The μOR is encoded by the Oprm1 gene and is expressed in multiple brain regions that regulate reward and motivation, such as the nucleus accumbens (NAc). Oprm1 expression in NAc medium spiny neurons (MSNs) mediates opioid place preference, seeking, and consumption. However, recent single nucleus RNA sequencing (snRNA-seq) studies have revealed that multiple subpopulations of NAc neurons express Oprm1 mRNA, making it unclear which populations mediate diverse behaviors resulting from μOR activation. Using published snRNA-seq datasets from the rat NAc, we identified a novel population of MSNs that express the highest levels of Oprm1 of any NAc cell type. Here, we show that this population is selectively marked by expression of Chst9, a gene encoding a carbohydrate sulfotransferase. Notably, Chst9+ neurons exhibited more abundant expression of Oprm1 as compared to other cell types, and formed discrete cellular clusters along the medial and ventral borders of the NAc shell subregion. Moreover, CHST9 mRNA was also found to mark specific MSN populations in published human and primate snRNA-seq studies, indicating that this unique population may be conserved across species. Together, these results identify a spatially and transcriptionally distinct NAc neuron population characterized by the expression of Chst9. The abundant expression of Oprm1 in this population and the conservation of these cells across species suggests that they may play a key functional role in opioid response and identify this subpopulation as a target for further investigation.

Omics Approaches to Investigate the Pathogenesis of Suicide

Suicide is the second leading cause of death in U.S. adolescents and young adults and is generally associated with a psychiatric disorder. Suicidal behavior has a complex etiology and pathogenesis. Moderate heritability suggests genetic causes. Associations between childhood and recent life adversity indicate contributions from epigenetic factors. Genomic contributions to suicide pathogenesis remain largely unknown. This article is based on a workshop held to design strategies to identify molecular drivers of suicide neurobiology that would be putative new treatment targets. The panel determined that while bulk tissue studies provide comprehensive information, single-nucleus approaches that identify cell type-specific changes are needed. While single-nuclei techniques lack information on cytoplasm, processes, spines, and synapses, spatial multiomic technologies on intact tissue detect cell alterations specific to brain tissue layers and subregions. Because suicide has genetic and environmental drivers, multiomic approaches that combine cell type-specific epigenome, transcriptome, and proteome provide a more complete picture of pathogenesis. To determine the direction of effect of suicide risk gene variants on RNA and protein expression and how these interact with epigenetic marks, single-nuclei and spatial multiomics quantitative trait loci maps should be integrated with whole-genome sequencing and genome-wide association databases. The workshop concluded with a recommendation for the formation of an international suicide biology consortium that will bring together brain banks and investigators with expertise in cutting-edge omics technologies to delineate the biology of suicide and identify novel potential treatment targets to be tested in cellular and animal models for drug and biomarker discovery to guide suicide prevention.

A single-nucleus transcriptomic atlas of medium spiny neurons in the rat nucleus accumbens

Neural processing of rewarding stimuli involves several distinct regions, including the nucleus accumbens (NAc). The majority of NAc neurons are GABAergic projection neurons known as medium spiny neurons (MSNs). MSNs are broadly defined by dopamine receptor expression, but evidence suggests that a wider array of subtypes exist. To study MSN heterogeneity, we analyzed single-nucleus RNA sequencing data from the largest available rat NAc dataset. Analysis of 48,040 NAc MSN nuclei identified major populations belonging to the striosome and matrix compartments. Integration with mouse and human data indicated consistency across species and disease-relevance scoring using genome-wide association study results revealed potentially differential roles for MSN populations in substance use disorders. Additional high-resolution clustering identified 34 transcriptomically distinct subtypes of MSNs definable by a limited number of marker genes. Together, these data demonstrate the diversity of MSNs in the NAc and provide a basis for more targeted genetic manipulation of specific populations.

Lipidome atlas of the adult human brain

Lipids are the most abundant but poorly explored components of the human brain. Here, we present a lipidome map of the human brain comprising 75 regions, including 52 neocortical ones. The lipidome composition varies greatly among the brain regions, affecting 93% of the 419 analyzed lipids. These differences reflect the brain's structural characteristics, such as myelin content (345 lipids) and cell type composition (353 lipids), but also functional traits: functional connectivity (76 lipids) and information processing hierarchy (60 lipids). Combining lipid composition and mRNA expression data further enhances functional connectivity association. Biochemically, lipids linked with structural and functional brain features display distinct lipid class distribution, unsaturation extent, and prevalence of omega-3 and omega-6 fatty acid residues. We verified our conclusions by parallel analysis of three adult macaque brains, targeted analysis of 216 lipids, mass spectrometry imaging, and lipidome assessment of sorted murine neurons.

Systems biology dissection of PTSD and MDD across brain regions, cell types, and blood

The molecular pathology of stress-related disorders remains elusive. Our brain multiregion, multiomic study of posttraumatic stress disorder (PTSD) and major depressive disorder (MDD) included the central nucleus of the amygdala, hippocampal dentate gyrus, and medial prefrontal cortex (mPFC). Genes and exons within the mPFC carried most disease signals replicated across two independent cohorts. Pathways pointed to immune function, neuronal and synaptic regulation, and stress hormones. Multiomic factor and gene network analyses provided the underlying genomic structure. Single nucleus RNA sequencing in dorsolateral PFC revealed dysregulated (stress-related) signals in neuronal and non-neuronal cell types. Analyses of brain-blood intersections in >50,000 UK Biobank participants were conducted along with fine-mapping of the results of PTSD and MDD genome-wide association studies to distinguish risk from disease processes. Our data suggest shared and distinct molecular pathology in both disorders and propose potential therapeutic targets and biomarkers.

A data-driven single-cell and spatial transcriptomic map of the human prefrontal cortex

The molecular organization of the human neocortex historically has been studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally defined spatial domains that move beyond classic cytoarchitecture. We used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex. Integration with paired single-nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we mapped the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains.

Dopamine signaling enriched striatal gene set predicts striatal dopamine synthesis and physiological activity in vivo

The polygenic architecture of schizophrenia implicates several molecular pathways involved in synaptic function. However, it is unclear how polygenic risk funnels through these pathways to translate into syndromic illness. Using tensor decomposition, we analyze gene co-expression in the caudate nucleus, hippocampus, and dorsolateral prefrontal cortex of post-mortem brain samples from 358 individuals. We identify a set of genes predominantly expressed in the caudate nucleus and associated with both clinical state and genetic risk for schizophrenia that shows dopaminergic selectivity. A higher polygenic risk score for schizophrenia parsed by this set of genes predicts greater dopamine synthesis in the striatum and greater striatal activation during reward anticipation. These results translate dopamine-linked genetic risk variation into in vivo neurochemical and hemodynamic phenotypes in the striatum that have long been implicated in the pathophysiology of schizophrenia.

An integrated single-nucleus and spatial transcriptomics atlas reveals the molecular landscape of the human hippocampus

The hippocampus contains many unique cell types, which serve the structure's specialized functions, including learning, memory and cognition. These cells have distinct spatial topography, morphology, physiology, and connectivity, highlighting the need for transcriptome-wide profiling strategies that retain cytoarchitectural organization. Here, we generated spatially-resolved transcriptomics (SRT) and single-nucleus RNA-sequencing (snRNA-seq) data from adjacent tissue sections of the anterior human hippocampus across ten adult neurotypical donors. We defined molecular profiles for hippocampal cell types and spatial domains. Using non-negative matrix factorization and transfer learning, we integrated these data to define gene expression patterns within the snRNA-seq data and infer the expression of these patterns in the SRT data. With this approach, we leveraged existing rodent datasets that feature information on circuit connectivity and neural activity induction to make predictions about axonal projection targets and likelihood of ensemble recruitment in spatially-defined cellular populations of the human hippocampus. Finally, we integrated genome-wide association studies with transcriptomic data to identify enrichment of genetic components for neurodevelopmental, neuropsychiatric, and neurodegenerative disorders across cell types, spatial domains, and gene expression patterns of the human hippocampus. To make this comprehensive molecular atlas accessible to the scientific community, both raw and processed data are freely available, including through interactive web applications.

Cell-type and sex-specific rhythmic gene expression in the nucleus accumbens

Circadian rhythms are critical for human health and are highly conserved across species. Disruptions in these rhythms contribute to many diseases, including psychiatric disorders. Previous results suggest that circadian genes modulate behavior through specific cell types in the nucleus accumbens (NAc), particularly dopamine D1-expressing medium spiny neurons (MSNs). However, diurnal rhythms in transcript expression have not been investigated in NAc MSNs. In this study we identified and characterized rhythmic transcripts in D1- and D2-expressing neurons and compared rhythmicity results to homogenate as well as astrocyte samples taken from the NAc of male and female mice. We find that all cell types have transcripts with diurnal rhythms and that top rhythmic transcripts are largely core clock genes, which peak at approximately the same time of day in each cell type and sex. While clock-controlled rhythmic transcripts are enriched for protein regulation pathways across cell type, cell signaling and signal transduction related processes are most commonly enriched in MSNs. In contrast to core clock genes, these clock-controlled rhythmic transcripts tend to reach their peak in expression about 2-h later in females than males, suggesting diurnal rhythms in reward may be delayed in females. We also find sex differences in pathway enrichment for rhythmic transcripts peaking at different times of day. Protein folding and immune responses are enriched in transcripts that peak in the dark phase, while metabolic processes are primarily enriched in transcripts that peak in the light phase. Importantly, we also find that several classic markers used to categorize MSNs are rhythmic in the NAc. This is critical since the use of rhythmic markers could lead to over- or under-enrichment of targeted cell types depending on the time at which they are sampled. This study greatly expands our knowledge of how individual cell types contribute to rhythms in the NAc.

Transcriptomic characterization of human lateral septum neurons reveals conserved and divergent marker genes across species

The lateral septum (LS) is a midline, subcortical structure, which regulates social behaviors that are frequently impaired in neurodevelopmental disorders including schizophrenia and autism spectrum disorder. Mouse studies have identified neuronal populations within the LS that express a variety of molecular markers, including vasopressin receptor, oxytocin receptor, and corticotropin releasing hormone receptor, that control specific facets of social behavior. Despite its critical role in the regulation of social behavior and notable gene expression patterns, comprehensive molecular profiling of the human LS has not been performed. Here, we conducted single nucleus RNA-sequencing (snRNA-seq) to generate the first transcriptomic profiles of the human LS using postmortem human brain tissue samples from 3 neurotypical donors. Our analysis identified 4 transcriptionally distinct neuronal cell types within the human LS that are enriched for TRPC4, the gene encoding Trp-related protein 4. Differential expression analysis revealed a distinct LS neuronal cell type that is enriched for OPRM1, the gene encoding the μ-opioid receptor. Leveraging recently collected mouse LS snRNA-seq datasets, we also conducted a cross-species analysis. Our results demonstrate that TRPC4 enrichment in the LS is highly conserved between human and mouse, while FREM2, which encodes FRAS1 related extracellular matrix protein 2, is enriched only in the human LS. Together, these results highlight transcriptional heterogeneity of the human LS, and identify robust marker genes for the human LS.

Differentially Altered Metabolic Pathways in the Amygdala of Subjects with Schizophrenia, Bipolar Disorder and Major Depressive Disorder

Background and hypothesis

A growing number of studies implicate a key role for metabolic processes in psychiatric disorders. Recent studies suggest that ketogenic diet may be therapeutically effective for subgroups of people with schizophrenia (SCZ), bipolar disorder (BPD) and possibly major depressive disorder (MDD). Despite this promise, there is currently limited information regarding brain energy metabolism pathways across these disorders, limiting our understanding of how brain metabolic pathways are altered and who may benefit from ketogenic diets. We conducted gene expression profiling on the amygdala, a key region involved in in the regulation of mood and appetitive behaviors, to test the hypothesis that amygdala metabolic pathways are differentially altered between these disorders.

Study Design

We used a cohort of subjects diagnosed with SCZ, BPD or MDD, and non-psychiatrically ill control subjects (n=15/group), together with our bioinformatic 3-pod analysis consisting of full transcriptome pathway analysis, targeted pathway analysis, leading-edge gene analysis and iLINCS perturbagen analysis.

Study Results

We identified differential expression of metabolic pathways in each disorder. Subjects with SCZ displayed downregulation of mitochondrial respiration and nucleotide metabolism pathways. In comparison, we observed upregulation of mitochondrial respiration pathways in subjects with MDD, while subjects with BPD displayed enrichment of pathways involved in carbohydrate metabolism. Several pathways associated with brain metabolism including immune system processes and calcium ion transport were also differentially altered between diagnosis groups.

Conclusion

Our findings suggest metabolic pathways are differentially altered in the amygdala in these disorders, which may impact approaches for therapeutic strategies.

Drugs of abuse hijack a mesolimbic pathway that processes homeostatic need

Drugs of abuse are thought to promote addiction in part by "hijacking" brain reward systems, but the underlying mechanisms remain undefined. Using whole-brain FOS mapping and in vivo single-neuron calcium imaging, we found that drugs of abuse augment dopaminoceptive ensemble activity in the nucleus accumbens (NAc) and disorganize overlapping ensemble responses to natural rewards in a cell type-specific manner. Combining FOS-Seq, CRISPR-perturbation, and single-nucleus RNA sequencing, we identified Rheb as a molecular substrate that regulates cell type-specific signal transduction in NAc while enabling drugs to suppress natural reward consumption. Mapping NAc-projecting regions activated by drugs of abuse revealed input-specific effects on natural reward consumption. These findings characterize the dynamic, molecular and circuit basis of a common reward pathway, wherein drugs of abuse interfere with the fulfillment of innate needs.

lute: estimating the cell composition of heterogeneous tissue with varying cell sizes using gene expression

Relative cell type fraction estimates in bulk RNA-sequencing data are important to control for cell composition differences across heterogenous tissue samples. Current computational tools estimate relative RNA abundances rather than cell type proportions in tissues with varying cell sizes, leading to biased estimates. We present lute, a computational tool to accurately deconvolute cell types with varying sizes. Our software wraps existing deconvolution algorithms in a standardized framework. Using simulated and real datasets, we demonstrate how lute adjusts for differences in cell sizes to improve the accuracy of cell composition. Software is available from https://bioconductor.org/packages/lute.

Morphological heterogeneity of neurons in the human central amygdaloid nucleus

The central amygdaloid nucleus (CeA) has an ancient phylogenetic development and functions relevant for animal survival. Local cells receive intrinsic amygdaloidal information that codes emotional stimuli of fear, integrate them, and send cortical and subcortical output projections that prompt rapid visceral and social behavior responses. We aimed to describe the morphology of the neurons that compose the human CeA (N = 8 adult men). Cells within CeA coronal borders were identified using the thionine staining and were further analyzed using the "single-section" Golgi method followed by open-source software procedures for two-dimensional and three-dimensional image reconstructions. Our results evidenced varied neuronal cell body features, number and thickness of primary shafts, dendritic branching patterns, and density and shape of dendritic spines. Based on these criteria, we propose the existence of 12 morphologically different spiny neurons in the human CeA and discuss the variability in the dendritic architecture within cellular types, including likely interneurons. Some dendritic shafts were long and straight, displayed few collaterals, and had planar radiation within the coronal neuropil volume. Most of the sampled neurons showed a few to moderate density of small stubby/wide spines. Long spines (thin and mushroom) were observed occasionally. These novel data address the synaptic processing and plasticity in the human CeA. Our morphological description can be combined with further transcriptomic, immunohistochemical, and electrophysiological/connectional approaches. It serves also to investigate how neurons are altered in neurological and psychiatric disorders with hindered emotional perception, in anxiety, following atrophy in schizophrenia, and along different stages of Alzheimer's disease.

Blood extracellular vesicles carrying brain-specific mRNAs are potential biomarkers for detecting gene expression changes in the female brain

The absence of non-invasive tests that can monitor the status of the brain is a major obstacle for psychiatric care. In order to address this need, we assessed the feasibility of using tissue-specific gene expression to determine the origin of extracellular vesicle (EV) mRNAs in peripheral blood. Using the placenta as a model, we discovered that 26 messenger RNAs that are specifically expressed in the placenta are present in EVs circulating in maternal blood. Twenty-three of these transcripts were either exclusively or highly expressed in maternal blood during pregnancy only and not in the postpartum period, verifying the feasibility of using tissue-specific gene expression to infer the tissue of origin for EV mRNAs. Using the same bioinformatic approach, which provides better specificity than isolating L1 cell-adhesion molecule containing EVs, we discovered that 181 mRNAs that are specifically expressed in the female brain are also present in EVs circulating in maternal blood. Gene set enrichment analysis revealed that these transcripts, which are involved in synaptic functions and myelination, are enriched for genes implicated in mood disorders, schizophrenia, and substance use disorders. The EV mRNA levels of 13 of these female brain-specific transcripts are associated with postpartum depression (adjusted p-vals = 3 × 10-5 to 0.08), raising the possibility that they can be used to infer the state of the brain. In order to determine the extent to which EV mRNAs reflect transcription in the brain, we compared mRNAs isolated from cells and EVs in an iPSC-derived brain microphysiological system differentiated for 3 and 9 weeks. We discovered that, although cellular and extracellular mRNA levels are not identical, they do correlate, and it is possible to extrapolate cellular RNA expression changes in the brain via EV mRNA levels. Our findings bring EV mRNAs to the forefront of peripheral biomarker development efforts in psychiatric diseases by demonstrating the feasibility of inferring transcriptional changes in the brain via blood EV mRNA levels.

Multi-slice spatial transcriptome domain analysis with SpaDo

With the rapid advancements in spatial transcriptome sequencing, multiple tissue slices are now available, enabling the integration and interpretation of spatial cellular landscapes. Herein, we introduce SpaDo, a tool for multi-slice spatial domain analysis, including modules for multi-slice spatial domain detection, reference-based annotation, and multiple slice clustering at both single-cell and spot resolutions. We demonstrate SpaDo's effectiveness with over 40 multi-slice spatial transcriptome datasets from 7 sequencing platforms. Our findings highlight SpaDo's potential to reveal novel biological insights in multi-slice spatial transcriptomes.

Thyroid hormone action in adult neurogliogenic niches: the known and unknown

Over the last decades, thyroid hormones (THs) signaling has been established as a key signaling cue for the proper maintenance of brain functions in adult mammals, including humans. One of the most fascinating roles of THs in the mature mammalian brain is their ability to regulate adult neurogliogenic processes. In this respect, THs control the generation of new neuronal and glial progenitors from neural stem cells (NSCs) as well as their final differentiation and maturation programs. In this review, we summarize current knowledge on the cellular organization of adult rodent neurogliogenic niches encompassing well-established niches in the subventricular zone (SVZ) lining the lateral ventricles, the hippocampal subgranular zone (SGZ), and the hypothalamus, but also less characterized niches in the striatum and the cerebral cortex. We then discuss critical questions regarding how THs availability is regulated in the respective niches in rodents and larger mammals as well as how modulating THs availability in those niches interferes with lineage decision and progression at the molecular, cellular, and functional levels. Based on those alterations, we explore the novel therapeutic avenues aiming at harnessing THs regulatory influences on neurogliogenic output to stimulate repair processes by influencing the generation of either new neurons (i.e. Alzheimer's, Parkinson's diseases), oligodendrocytes (multiple sclerosis) or both (stroke). Finally, we point out future challenges, which will shape research in this exciting field in the upcoming years.

Drug targeting in psychiatric disorders - how to overcome the loss in translation?

In spite of major efforts and investment in development of psychiatric drugs, many clinical trials have failed in recent decades, and clinicians still prescribe drugs that were discovered many years ago. Although multiple reasons have been discussed for the drug development deadlock, we focus here on one of the major possible biological reasons: differences between the characteristics of drug targets in preclinical models and the corresponding targets in patients. Importantly, based on technological advances in single-cell analysis, we propose here a framework for the use of available and newly emerging knowledge from single-cell and spatial omics studies to evaluate and potentially improve the translational predictivity of preclinical models before commencing preclinical and, in particular, clinical studies. We believe that these recommendations will improve preclinical models and the ability to assess drugs in clinical trials, reducing failure rates in expensive late-stage trials and ultimately benefitting psychiatric drug discovery and development.

Single-Nucleus RNA-Seq Characterizes the Cell Types Along the Neuronal Lineage in the Adult Human Subependymal Zone and Reveals Reduced Oligodendrocyte Progenitor Abundance with Age

The subependymal zone (SEZ), also known as the subventricular zone (SVZ), constitutes a neurogenic niche that persists during postnatal life. In humans, the neurogenic potential of the SEZ declines after the first year of life. However, studies discovering markers of stem and progenitor cells highlight the neurogenic capacity of progenitors in the adult human SEZ, with increased neurogenic activity occurring under pathological conditions. In the present study, the complete cellular niche of the adult human SEZ was characterized by single-nucleus RNA sequencing, and compared between four youth (age 16-22) and four middle-aged adults (age 44-53). We identified 11 cellular clusters including clusters expressing marker genes for neural stem cells (NSCs), neuroblasts, immature neurons, and oligodendrocyte progenitor cells. The relative abundance of NSC and neuroblast clusters did not differ between the two age groups, indicating that the pool of SEZ NSCs does not decline in this age range. The relative abundance of oligodendrocyte progenitors and microglia decreased in middle-age, indicating that the cellular composition of human SEZ is remodeled between youth and adulthood. The expression of genes related to nervous system development was higher across different cell types, including NSCs, in youth as compared with middle-age. These transcriptional changes suggest ongoing central nervous system plasticity in the SEZ in youth, which declined in middle-age.

Single nuclei transcriptomics in human and non-human primate striatum in opioid use disorder

In brain, the striatum is a heterogenous region involved in reward and goal-directed behaviors. Striatal dysfunction is linked to psychiatric disorders, including opioid use disorder (OUD). Striatal subregions are divided based on neuroanatomy, each with unique roles in OUD. In OUD, the dorsal striatum is involved in altered reward processing, formation of habits, and development of negative affect during withdrawal. Using single nuclei RNA-sequencing, we identified both canonical (e.g., dopamine receptor subtype) and less abundant cell populations (e.g., interneurons) in human dorsal striatum. Pathways related to neurodegeneration, interferon response, and DNA damage were significantly enriched in striatal neurons of individuals with OUD. DNA damage markers were also elevated in striatal neurons of opioid-exposed rhesus macaques. Sex-specific molecular differences in glial cell subtypes associated with chronic stress were found in OUD, particularly female individuals. Together, we describe different cell types in human dorsal striatum and identify cell type-specific alterations in OUD.

Spatial transcriptomics and single-nucleus RNA sequencing reveal a transcriptomic atlas of adult human spinal cord

Despite the recognized importance of the spinal cord in sensory processing, motor behaviors, and neural diseases, the underlying organization of neuronal clusters and their spatial location remain elusive. Recently, several studies have attempted to define the neuronal types and functional heterogeneity in the spinal cord using single-cell or single-nucleus RNA sequencing in animal models or developing humans. However, molecular evidence of cellular heterogeneity in the adult human spinal cord is limited. Here, we classified spinal cord neurons into 21 subclusters and determined their distribution from nine human donors using single-nucleus RNA sequencing and spatial transcriptomics. Moreover, we compared the human findings with previously published single-nucleus data of the adult mouse spinal cord, which revealed an overall similarity in the neuronal composition of the spinal cord between the two species while simultaneously highlighting some degree of heterogeneity. Additionally, we examined the sex differences in the spinal neuronal subclusters. Several genes, such as SCN10A and HCN1, showed sex differences in motor neurons. Finally, we classified human dorsal root ganglia (DRG) neurons using spatial transcriptomics and explored the putative interactions between DRG and spinal cord neuronal subclusters. In summary, these results illustrate the complexity and diversity of spinal neurons in humans and provide an important resource for future research to explore the molecular mechanisms underlying spinal cord physiology and diseases.

The gene expression landscape of the human locus coeruleus revealed by single-nucleus and spatially-resolved transcriptomics

Norepinephrine (NE) neurons in the locus coeruleus (LC) make long-range projections throughout the central nervous system, playing critical roles in arousal and mood, as well as various components of cognition including attention, learning, and memory. The LC-NE system is also implicated in multiple neurological and neuropsychiatric disorders. Importantly, LC-NE neurons are highly sensitive to degeneration in both Alzheimer's and Parkinson's disease. Despite the clinical importance of the brain region and the prominent role of LC-NE neurons in a variety of brain and behavioral functions, a detailed molecular characterization of the LC is lacking. Here, we used a combination of spatially-resolved transcriptomics and single-nucleus RNA-sequencing to characterize the molecular landscape of the LC region and the transcriptomic profile of LC-NE neurons in the human brain. We provide a freely accessible resource of these data in web-accessible and downloadable formats.

TrkB-dependent regulation of molecular signaling across septal cell types

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning, and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Identifying novel gene dysregulation associated with opioid overdose death: A meta-analysis of differential gene expression in human prefrontal cortex

Only recently have human postmortem brain studies of differential gene expression (DGE) associated with opioid overdose death (OOD) been published; sample sizes from these studies have been modest (N = 40-153). To increase statistical power to identify OOD-associated genes, we leveraged human prefrontal cortex RNAseq data from four independent OOD studies and conducted a transcriptome-wide DGE meta-analysis (N = 285). Using a unified gene expression data processing and analysis framework across studies, we meta-analyzed 20 098 genes and found 335 significant differentially expressed genes (DEGs) by OOD status (false discovery rate < 0.05). Of these, 66 DEGs were among the list of 303 genes reported as OOD-associated in prior prefrontal cortex molecular studies, including genes/gene families (e.g., OPRK1, NPAS4, DUSP, EGR). The remaining 269 DEGs were not previously reported (e.g., NR4A2, SYT1, HCRTR2, BDNF). There was little evidence of genetic drivers for the observed differences in gene expression between opioid addiction cases and controls. Enrichment analyses for the DEGs across molecular pathway and biological process databases highlight an interconnected set of genes and pathways from orexin and tyrosine kinase receptors through MEK/ERK/MAPK signaling to affect neuronal plasticity.

Cell-type brain-region specific changes in prefrontal cortex of a mouse model of alcohol dependence

The prefrontal cortex is a crucial regulator of alcohol drinking, and dependence, and other behavioral phenotypes associated with AUD. Comprehensive identification of cell-type specific transcriptomic changes in alcohol dependence will improve our understanding of mechanisms underlying the excessive alcohol use associated with alcohol dependence and will refine targets for therapeutic development. We performed single nucleus RNA sequencing (snRNA-seq) and Visium spatial gene expression profiling on the medial prefrontal cortex (mPFC) obtained from C57BL/6 J mice exposed to the two-bottle choice-chronic intermittent ethanol (CIE) vapor exposure (2BC-CIE, defined as dependent group) paradigm which models phenotypes of alcohol dependence including escalation of alcohol drinking. Gene co-expression network analysis and differential expression analysis identified highly dysregulated co-expression networks in multiple cell types. Dysregulated modules and their hub genes suggest novel understudied targets for studying molecular mechanisms contributing to the alcohol dependence state. A subtype of inhibitory neurons was the most alcohol-sensitive cell type and contained a downregulated gene co-expression module; the hub gene for this module is Cpa6, a gene previously identified by GWAS to be associated with excessive alcohol consumption. We identified an astrocytic Gpc5 module significantly upregulated in the alcohol-dependent group. To our knowledge, there are no studies linking Cpa6 and Gpc5 to the alcohol-dependent phenotype. We also identified neuroinflammation related gene expression changes in multiple cell types, specifically enriched in microglia, further implicating neuroinflammation in the escalation of alcohol drinking. Here, we present a comprehensive atlas of cell-type specific alcohol dependence mediated gene expression changes in the mPFC and identify novel cell type-specific targets implicated in alcohol dependence.

Sex differences in bipolar disorder: The dorsolateral prefrontal cortex as an etiopathogenic region

Bipolar disorder (BD) is worldwide a prevalent mental illness and a leading risk factor for suicide. Over the past three decades, it has been discovered that sex differences exist throughout the entire panorama of BD, but the etiologic regions and mechanisms that generate such differences remain poorly characterized. Available evidence indicates that the dorsolateral prefrontal cortex (DLPFC), a critical region that controls higher-order cognitive processing and mood, exhibits biological disparities between male and female patients with psychiatric disorders, which are highly correlated with the co-occurrence of psychotic symptoms. This review addresses the sex differences in BD concerning epidemiology, cognitive impairments, clinical manifestations, neuroimaging, and laboratory abnormalities. It also provides strong evidence linking DLPFC to the etiopathogenesis of these sex differences. We emphasize the importance of identifying gene signatures using human brain transcriptomics, which can depict sexually different variations, explain sex-biased symptomatic features, and provide novel targets for sex-specific therapeutics.

Automatic cell-type harmonization and integration across Human Cell Atlas datasets

Harmonizing cell types across the single-cell community and assembling them into a common framework is central to building a standardized Human Cell Atlas. Here, we present CellHint, a predictive clustering tree-based tool to resolve cell-type differences in annotation resolution and technical biases across datasets. CellHint accurately quantifies cell-cell transcriptomic similarities and places cell types into a relationship graph that hierarchically defines shared and unique cell subtypes. Application to multiple immune datasets recapitulates expert-curated annotations. CellHint also reveals underexplored relationships between healthy and diseased lung cell states in eight diseases. Furthermore, we present a workflow for fast cross-dataset integration guided by harmonized cell types and cell hierarchy, which uncovers underappreciated cell types in adult human hippocampus. Finally, we apply CellHint to 12 tissues from 38 datasets, providing a deeply curated cross-tissue database with ∼3.7 million cells and various machine learning models for automatic cell annotation across human tissues.

Genetic implication of prenatal GABAergic and cholinergic neuron development in susceptibility to schizophrenia

Background

The ganglionic eminences are fetal-specific structures that give rise to gamma-aminobutyric acid (GABA)- and acetylcholine- releasing neurons of the forebrain. Given evidence for GABAergic and cholinergic disturbances in schizophrenia, as well as an early neurodevelopmental component to the disorder, we tested the potential involvement of developing cells of the ganglionic eminences in mediating genetic risk for the condition.

Study Design

We combined data from a recent large-scale genome-wide association study of schizophrenia with single cell RNA sequencing data from the human ganglionic eminences to test enrichment of schizophrenia risk variation in genes with high expression specificity for particular developing cell populations within these structures. We additionally performed the single nuclei Assay for Transposase-Accessible Chromatin with Sequencing (snATAC-Seq) to map potential regulatory genomic regions operating in individual cell populations of the human ganglionic eminences, using these to additionally test for enrichment of schizophrenia common genetic variant liability and to functionally annotate non-coding variants associated with the disorder.

Study Results

Schizophrenia common variant liability was enriched in genes with high expression specificity for developing neuron populations that are predicted to form dopamine D1 and D2 receptor expressing GABAergic medium spiny neurons of the striatum, cortical somatostatin-positive GABAergic interneurons, calretinin-positive GABAergic neurons and cholinergic neurons. Consistent with these findings, schizophrenia genetic risk was also concentrated in predicted regulatory genomic sequence mapped in developing neuronal populations of the ganglionic eminences.

Conclusions

Our study provides evidence for a role of prenatal GABAergic and cholinergic neuron development in later susceptibility to schizophrenia.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single-cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding disease pathologies. However, several experimental and computational challenges impede transcriptomics-based deconvolution approaches using single-cell/nucleus RNA-seq reference atlases. Cells from the brain and blood have substantially different sizes, total mRNA, and transcriptional activities, and existing approaches may quantify total mRNA instead of cell type proportions. Further, standards are lacking for the use of cell reference atlases and integrative analyses of single-cell and spatial transcriptomics data. We discuss how to approach these key challenges with orthogonal "gold standard" datasets for evaluating deconvolution methods.