Svitlana V. Bach

Major depressive disorder (MDD) and bipolar disorder (BD) exhibit both overlapping and distinct clinical features, including sex differences in prevalence and symptomatology, but how these similarities and differences are reflected at the molecular and cellular level in the human brain remains poorly understood. Here, we applied spatial transcriptomics to the human dorsolateral prefrontal cortex from individuals with MDD, BD, and neurotypical controls to resolve molecular programs associated with diagnosis and biological sex. Our analysis revealed sex- and diagnosis-specific differences across white matter, the neurovascular unit (NVU), and neuronal compartments. Within white matter, we identified female-specific survival-response programs shared in both MDD and BD, and divergent myelin-related gene expression patterns across disorders ( CLDN11 , CNDP1 ). Within the NVU, transcriptional evidence of impaired vascular integrity was accompanied by astrocytic response profiles in females with MDD and males with BD, with inflammatory responses further evident in females with MDD ( IFITM3 , CEBPD , IL1R1 ). In donors with BD of both sexes and in females with MDD, inhibitory neurons exhibited impaired neuropeptide expression ( SST , CORT , CRH ). Integration with donor genotype information highlighted that genetic risk for mood disorders influences MAPK3 expression, implicating altered ERK signaling, which was further supported by disorder-associated changes in Ras/MEK/ERK pathway components in neurons ( MAPK3, DUSP6 , RASD1 , ELK1 ). Together, these findings reveal coordinated, cell type- and spatial domain-specific molecular programs that distinguish shared and disorder-specific transcriptomic changes across mood disorders in a sex-dependent manner.

The habenula (Hb) is a small epithalamic midline structure critical for reward processing and mood regulation that is implicated in many complex brain disorders, including depression and addiction. The molecular anatomy of the Hb has been well-characterized in model organisms, but little is known about its organization in the human brain, including the transcriptomic and epigenetic signatures of spatially-organized lateral (LHb) and medial (MHb) cell types. Here, we generated a high resolution spatiomolecular map of the human Hb using both multiomic single nucleus sequencing (snMultiome) and spatially-resolved transcriptomics. We measured gene expression and chromatin accessibility across 4 MHb and 5 LHb neuronal cell types as defined by transcriptomic conservation with rodents and zebrafish. We identified cell type-specific molecular specializations related to GABAergic and glutamatergic signaling as well as regulatory relationships between open chromatin regions, transcription factors, and target genes associated with depression and addiction. High resolution spatial transcriptomics identified the anatomical location of heterogeneous LHb subpopulations, including two putative inhibitory populations in the LHb. We further characterized cell-type relationships between Hb neurons and glial subtypes and predicted ligand-receptor interactions between LHb and MHb neurons. Finally, we delineated Hb cell types and spatial domains enriched in genes associated with genetic risk for psychiatric disorders, identifying the strongest enrichment in the MHb. Collectively, the generated data and analyses resolve the spatiomolecular organization of the human Hb and identify molecular-genetic associations with neuropsychiatric disease.

Spatial transcriptomics has transformed molecular characterization of the human brain, but most studies profile individual two-dimensional tissue sections and therefore cannot capture molecular organization across 3D neuroanatomical axes. Here, we developed a spatial genomics framework to reconstruct the 3D cellular architecture of the human nucleus accumbens (NAc) by densely sampling and aligning serial sections across its anterior-posterior (AP) extent. Combining single-cell Xenium profiling with near-transcriptome-wide VisiumHD, we mapped cell type composition and topography across AP, mediolateral and dorsoventral axes. Neuronal populations varied in relative abundance across the AP axis while maintaining characteristic spatial topographies. These differences were particularly pronounced among D1-islands, a topographically organized structure within the NAc in which transcriptionally distinct subpopulations of DRD1 -expressing medium spiny neurons (MSNs) aggregate into discrete cellular formations. D1-islands comprised two molecularly distinct spatial domains with conserved counterparts across species. Both domains varied in relative abundance across the AP axis, but demonstrated distinct 3D distribution; one was present throughout the NAc, whereas the other, corresponding to islands of Calleja (ICj), was more restricted to intermediate-to-posterior levels. Local spatial analyses revealed region- and AP-dependent cellular relationships, including segregation of the two island populations despite their anatomical proximity. Together, these findings establish a framework for three dimensional spatial genomics of human brain tissue and reveal molecular and cellular organization not captured by individual tissue sections.

The amygdala is a central hub for emotional learning that is dysregulated across numerous psychiatric disorders, yet its molecular architecture in humans remains poorly defined. We here present a spatial transcriptomic atlas of the human amygdala from nine neurotypical donors. Integrating Visium, Xenium, and VisiumHD technologies, we profiled over 1.1 million cells across 13 spatial domains spanning the basolateral complex, central, medial, and cortical nuclei, as well as the intercalated islands, each defined by distinct marker genes and gene co-expression networks. By integrating snRNA-seq reference atlases, we found that each subdivision of the basolateral subnucleus contains distinct excitatory neuron classes. We additionally found three transcriptionally distinct populations of intercalated neuron types, two of which form spatially distinct islands neighboring the basolateral complex, and were able to refine cell type diversity across the central nucleus and related amygdalostriatal transition areas. Finally, we localized psychiatric genetic risk across the amygdala, revealing both broad neuronal enrichment and subnuclear specificity. Together, this atlas provides a foundational resource and framework for accelerating cross-species comparisons and human disease-focused investigations.

Summary: The dorsolateral prefrontal cortex (dlPFC) is central to cognitive dysfunction in schizophrenia (SCZ), yet how molecular changes are organized across cortex remains unclear. Here, we applied complementary spatial transcriptomic approaches spanning laminar domains, microenvironments, and cell types in postmortem human dlPFC. At the laminar level, SCZ-associated transcriptional changes were strongest in glia-enriched domains (layer 1/meninges and white matter), including down-regulation of microglia-associated genes, whereas genetic risk localized to neuronal-rich gray matter. We next analyzed SCZ-linked microenvironments including neuropil, neuronal, perineuronal net, and vascular compartments. Within these, neuronal and synaptic transcriptional changes were most prominent in neuropil, showing down-regulation of activity-dependent synaptic genes and inhibitory neuron markers. At the cellular level, these signals reflected intrinsic alterations localized to cell types. Across analyses, patterns converged on altered BDNF-TrkB signaling and inhibitory circuit dysfunction. Together, our findings highlight spatial scale as a key determinant in resolving neuronal and non-neuronal aspects of SCZ-associated biology.

The entorhinal cortex (ERC) is implicated in early progression of Alzheimer’s disease (AD). Here we investigated the impact of established biological risk factors for AD, including APOE genotype (E2 versus E4 alleles), sex, and ancestry, on gene expression in the human ERC. We generated paired spatially-resolved transcriptomics (SRT) and single-nucleus RNA sequencing data (snRNA-seq) in postmortem human ERC tissue from middle aged brain donors with no history of AD. APOE -dependent changes in gene expression predominantly mapped to a transcriptionally-defined oligodendrocyte subtype, which varied substantially with ancestry, and suggested differences in oligodendrocyte differentiation and myelination. Integration of SRT and snRNA-seq data identified a common gene expression signature associated with APOE genotype, which we localized to the same oligodendrocyte subtype and a white matter spatial domain. This suggests that AD risk in ERC may be associated with disrupted oligodendrocyte function, potentially contributing to future neurodegeneration. Lay Summary: Alzheimer’s disease (AD) is a neurodegenerative disorder that accounts for 60-80% dementia cases. Apolipoprotein E ( APOE) genotype is the strongest genetic risk factor for AD, and the entorhinal cortex (ERC) is a brain region implicated in its earliest progression. Our study investigated how APOE genotype impacts gene expression in the ERC. We identified genotype-dependent effects on oligodendrocytes with different transcriptional profiles related to maturation that may help explain how APOE genotype mediates its effects on AD risk.

The locus coeruleus (LC) is a small noradrenergic nucleus in the dorsal pons that sends projections across the brain regulating sleep, arousal, attention, stress responses, and some forms of cognition. LC neurons show pathology in the earliest stages of Alzheimer’s disease (AD), including age-related accumulation of hyperphosphorylated tau (pTau) and accelerated loss of neuromelanin (NM) pigmentation. NM-sensitive neuroimaging of the LC predicts previous cognitive decline, clinical severity, and future AD progression. While these findings suggest that the LC plays an etiologic role in AD, the molecular landscape of the LC prior to clinical manifestation of sporadic AD remains largely uncharacterized. This information is critical for developing interventions that preserve LC integrity and function. We performed spatially-resolved transcriptomics on 85 sections of human postmortem LC from N =33 neurotypical middle-aged donors, balanced for epidemiologic AD risk factors including sex, African or European ancestry, and APOE genotype (carriers of the E4/risk or E2/protective alleles). Comparing across APOE genotypes, we find astrocytic gene expression differences proximal to LC neurons. Associating NM content with local gene expression, we show that higher overall APOE gene expression correlates with reduced NM content and an enrichment of NM-associated genes in aging pathways. Unexpectedly, we find enriched LC expression of cholesterol synthesis genes, alongside evidence for lipid synthesis gene regulatory network activity in NM-containing LC specifically, revealing a potential intersection between intrinsic lipid metabolism in LC neurons, NM, and the role of APOE-mediated lipid biology in AD. Together, these data illuminate the molecular features of the human LC at spatial resolution with unprecedented sampling depth, revealing how AD risk factors and NM content influence resilience and susceptibility of this critical brain nucleus to pathology accumulation and degeneration.

The amygdaloid complex mediates learning, memory, and emotions. Understanding cellular and anatomical features specialized in the primate amygdala versus other mammals requires a systematic, anatomically resolved molecular analysis of neuron types. We analyzed five nuclear subdivisions of the primate amygdala with single-nucleus RNA sequencing in macaques, baboons, and humans to examine gene expression profiles for excitatory and inhibitory neurons. Integrated analyses across species identified diverse subtypes of glutamatergic and GABAergic neurons that are highly conserved across primates. Compositional analyses revealed that subdivisions of the primate basolateral complex contain distinct classes of glutamatergic neurons and divergent gene expression profiles for parvalbumin and somatostatin GABAergic neurons. Referencing primate neuron types to transcriptomic atlases of the murine amygdala highlighted primate-specific specializations for glutamatergic neurons and identified which neuron types are conserved across mammals and amygdalar subdivisions. Understanding the molecular heterogeneity of anatomically resolved amygdalar neuron types provides a cellular framework for improving models of how amygdalar circuits contribute to cognition and mental health.

The nucleus accumbens (NAc) is a key component of the mesolimbic dopamine system that critically regulates many behaviors related to reward and motivation. The NAc is implicated in several neuropsychiatric disorders, including major depressive disorder, schizophrenia, and substance use disorders. Rodent studies have identified spatial organization of heterogeneous medium spiny neuron (MSN) subtypes across the NAc core and shell, but the extent to which this cellular diversity and spatial organization is conserved in the human brain remains unclear. Here, we generated a spatiomolecular atlas of NAc cell types and spatial domains by integrating spatial transcriptomics and single-nucleus RNA sequencing data from postmortem NAc tissue from 10 neurotypical adult donors. We identified 20 transcriptionally unique cell populations and 8 spatial domains, including specialized D1 islands composed of distinct dopamine receptor 1 (DRD1) MSN subtypes, which were enriched for OPRM1 . In contrast to a discrete core vs. shell division, we observed continuous spatial gradients of gene expression across MSN domains, suggesting a more complex organization of DRD1 and DRD2 MSNs. Cross-species comparisons demonstrated conservation of MSN subtypes and spatial features between human, rodent, and nonhuman primate NAc. Genetic enrichment analysis with stratified linkage disequilibrium score regression revealed specific spatial domains associated with risk for psychiatric and addiction-related traits. To investigate this further, we spatially mapped ligand-receptor interactions involving neuropsychiatric risk genes. Finally, we leveraged existing rodent NAc data to identify drug-responsive transcriptional programs and predict their spatial distribution in the human NAc. Collectively, we provide a spatiomolecular framework for understanding the human NAc and its relevance to neuropsychiatric disease.