Neuronal subtypes and diversity revealed by single-nucleus RNA sequencing of the human brain

Calcium-permeable AMPA receptors govern PV neuron feature selectivity

The brain helps us survive by forming internal representations of the external world1,2. Excitatory cortical neurons are often precisely tuned to specific external stimuli3,4. However, inhibitory neurons, such as parvalbumin-positive (PV) interneurons, are generally less selective5. PV interneurons differ from excitatory neurons in their neurotransmitter receptor subtypes, including AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptors (AMPARs)6,7. Excitatory neurons express calcium-impermeable AMPARs that contain the GluA2 subunit (encoded by GRIA2), whereas PV interneurons express receptors that lack the GluA2 subunit and are calcium-permeable (CP-AMPARs). Here we demonstrate a causal relationship between CP-AMPAR expression and the low feature selectivity of PV interneurons. We find low expression stoichiometry of GRIA2 mRNA relative to other subunits in PV interneurons that is conserved across ferrets, rodents, marmosets and humans, and causes abundant CP-AMPAR expression. Replacing CP-AMPARs in PV interneurons with calcium-impermeable AMPARs increased their orientation selectivity in the visual cortex. Manipulations to induce sparse CP-AMPAR expression demonstrated that this increase was cell-autonomous and could occur with changes beyond development. Notably, excitatory-PV interneuron connectivity rates and unitary synaptic strength were unaltered by CP-AMPAR removal, which suggested that the selectivity of PV interneurons can be altered without markedly changing connectivity. In Gria2-knockout mice, in which all AMPARs are calcium-permeable, excitatory neurons showed significantly degraded orientation selectivity, which suggested that CP-AMPARs are sufficient to drive lower selectivity regardless of cell type. Moreover, hippocampal PV interneurons, which usually exhibit low spatial tuning, became more spatially selective after removing CP-AMPARs, which indicated that CP-AMPARs suppress the feature selectivity of PV interneurons independent of modality. These results reveal a new role of CP-AMPARs in maintaining low-selectivity sensory representation in PV interneurons and implicate a conserved molecular mechanism that distinguishes this cell type in the neocortex.

Research on the shared function of central neurons and breast cancer based on gene expression profile data mining: The role of EMID1 protein antibody expression

In recent years, the incidence of breast cancer has gradually increased, and the research on it has become a hot spot in the scientific community. Central neurons play an important role in breast cancer. This study aims to explore the application of gene expression profile data mining in the study of shared function between central neurons and breast cancer, and focuses on the expression of EMID1 protein antibody. The study collected biomedical images and gene expression profile data of breast cancer patients. Then, we use image processing and analysis technology to extract and analyze features of biomedical images to obtain quantitative features of breast cancer. Gene expression profile data were preprocessed and analyzed to obtain information about breast cancer related genes. Integrating and fusing biomedical images and gene expression profile data, and exploring the sharing function between central neurons and breast cancer through data mining algorithms and statistical analysis methods. The results showed that the expression of EMID1 protein was high in breast cancer tissues, and the expression pattern was similar to that of central neurons. Further functional studies have shown that EMID1 protein is involved in the regulation of proliferation and invasion of breast cancer cells. By regulating the expression level of EMID1 protein, we observed that the proliferation and invasion ability of breast cancer cells were significantly affected. The research results show that through the comprehensive analysis of biomedical images and gene expression profile data, we found the sharing function between central neurons and breast cancer. The central neuronal cell marker genes EMID1 and GREB1L may be used as key biomarkers to regulate the pathogenesis of breast cancer and affect the occurrence and development of breast cancer.

Visualized lateral flow assay for logic determination of co-existing viral RNA fragments

Different types of pathogenic viruses that have common transmission path can be co-infected, inducing distinct disease procession in comparison to that infection of one. Also, in the post COVID-19 time, more types of respiratory infectious virus are becoming prevalent and are concurrent. Those bring an urgent need for detection of co-existing viruses. Here, we propose a visualized lateral flow assay for logic determination of co-existing viral RNA fragments. In the presence of specific viral RNA inputs, DNAzyme is de-blocked according to defined logic, and catalyzes the hydrolysis of hairpin-structural substrate. One of cleaved substrates contains DNAzyme domain to realize dual signal amplification, which obtains copious of the other cleaved substrates. The cleaved substrates act as linking strands for bridging DNA-modified gold nanoparticles onto lateral flow strip to induce coloration on test line. "AND", "OR" and "INHIBIT" controlled lateral flow assays are respectively demonstrated for co-existing viral RNA detection, and the visual results can be obtained by the same kind of prepared strip, without need of re-fabricating strips according to logic systems. The work provides a flexible, convenient, visual and logic-processing strategy for simultaneous analysis of co-existing viruses.

Regional patterns of human cortex development correlate with underlying neurobiology

Human brain morphology undergoes complex changes over the lifespan. Despite recent progress in tracking brain development via normative models, current knowledge of underlying biological mechanisms is highly limited. We demonstrate that human cortical thickness development and aging trajectories unfold along patterns of molecular and cellular brain organization, traceable from population-level to individual developmental trajectories. During childhood and adolescence, cortex-wide spatial distributions of dopaminergic receptors, inhibitory neurons, glial cell populations, and brain-metabolic features explain up to 50% of the variance associated with a lifespan model of regional cortical thickness trajectories. In contrast, modeled cortical thickness change patterns during adulthood are best explained by cholinergic and glutamatergic neurotransmitter receptor and transporter distributions. These relationships are supported by developmental gene expression trajectories and translate to individual longitudinal data from over 8000 adolescents, explaining up to 59% of developmental change at cohort- and 18% at single-subject level. Integrating neurobiological brain atlases with normative modeling and population neuroimaging provides a biologically meaningful path to understand brain development and aging in living humans.

Transcriptomic Analysis Reveals the Heterogeneous Role of Conducting Films Upon Electrical Stimulation

Central nervous system (CNS) injuries and neurodegenerative diseases have markedly poor prognoses and can result in permanent dysfunction due to the general inability of CNS neurons to regenerate. Differentiation of transplanted stem cells has emerged as a therapeutic avenue to regenerate tissue architecture in damaged areas. Electrical stimulation is a promising approach for directing the differentiation outcomes and pattern of outgrowth of transplanted stem cells, however traditional inorganic bio-electrodes can induce adverse effects such as inflammation. This study demonstrates the implementation of two organic thin films, a polymer/reduced graphene oxide nanocomposite (P(rGO)) and PEDOT:PSS, that have favorable properties for implementation as conductive materials for electrical stimulation, as well as an inorganic indium tin oxide (ITO) conductive film. Transcriptomic analysis reveals that electrical stimulation improves neuronal differentiation of SH-SY5Y cells on all three films, with the greatest effect for P(rGO). Unique material- and electrical stimuli-mediated effects are observed, associated with differentiation, cell-substrate adhesion, and translation. The work demonstrates that P(rGO) and PEDOT:PSS are highly promising organic materials for the development of biocompatible, conductive scaffolds that will enhance electrically-aided stem cell therapeutics for CNS injuries and neurodegenerative diseases.

The macroscale routing mechanism of structural brain connectivity related to body mass index

Understanding the brain's mechanisms in individuals with obesity is important for managing body weight. Prior neuroimaging studies extensively investigated alterations in brain structure and function related to body mass index (BMI). However, how the network communication among the large-scale brain networks differs across BMI is underinvestigated. This study used diffusion magnetic resonance imaging of 290 young adults to identify links between BMI and brain network mechanisms. Navigation efficiency, a measure of network routing, was calculated from the structural connectivity computed using diffusion tractography. The sensory and frontoparietal networks indicated positive associations between navigation efficiency and BMI. The neurotransmitter association analysis identified that serotonergic and dopaminergic receptors, as well as opioid and norepinephrine systems, were related to BMI-related alterations in navigation efficiency. The transcriptomic analysis found that genes associated with network routing across BMI overlapped with genes enriched in excitatory and inhibitory neurons, specifically, gene enrichments related to synaptic transmission and neuron projection. Our findings suggest a valuable insight into understanding BMI-related alterations in brain network routing mechanisms and the potential underlying cellular biology, which might be used as a foundation for BMI-based weight management.

Accurate identification of genes associated with brain disorders by integrating heterogeneous genomic data into a Bayesian framework

BACKGROUND

Genome-wide association studies (GWAS) have revealed many brain disorder-associated SNPs residing in the noncoding genome, rendering it a challenge to decipher the underlying pathogenic mechanisms.

METHODS

Here, we present an unsupervised Bayesian framework to identify disease-associated genes by integrating risk SNPs with long-range chromatin interactions (iGOAT), including SNP-SNP interactions extracted from ∼500,000 patients and controls from the UK Biobank, and enhancer-promoter interactions derived from multiple brain cell types at different developmental stages.

FINDINGS

The application of iGOAT to three psychiatric disorders and three neurodegenerative/neurological diseases predicted sets of high-risk (HRGs) and low-risk (LRGs) genes for each disorder. The HRGs were enriched in drug targets, and exhibited higher expression during prenatal brain developmental stages than postnatal stages, indicating their potential to affect brain development at an early stage. The HRGs associated with Alzheimer's disease were found to share genetic architecture with schizophrenia, bipolar disorder and major depressive disorder according to gene co-expression module analysis and rare variants analysis. Comparisons of this method to the eQTL-based method, the TWAS-based method, and the gene-level GWAS method indicated that the genes identified by our method are more enriched in known brain disorder-related genes, and exhibited higher precision. Finally, the method predicted 205 risk genes not previously reported to be associated with any brain disorder, of which one top-risk gene, MLH1, was experimentally validated as being schizophrenia-associated.

INTERPRETATION

iGOAT can successfully leverage epigenomic data, phenotype-genotype associations, and protein-protein interactions to advance our understanding of brain disorders, thereby facilitating the development of new therapeutic approaches.

FUNDING

The work was funded by the National Key Research and Development Program of China (2024YFF1204902), the Natural Science Foundation of China (82371482), Guangzhou Science and Technology Research Plan (2023A03J0659) and Natural Science Foundation of Guangdong (2024A1515011363).

Integrative genomics approach identifies glial transcriptomic dysregulation and risk in the cortex of individuals with Alcohol Use Disorder

Alcohol use disorder (AUD) is a prevalent neuropsychiatric disorder that is a major global health concern, affecting millions of people worldwide. Past molecular studies of AUD used underpowered single cell analysis or bulk homogenates of postmortem brain tissue, which obscures gene expression changes in specific cell types. Here we performed single nuclei RNA-sequencing analysis of 73 post-mortem samples from individuals with AUD (N=36, N nuclei = 248,873) and neurotypical controls (N=37, N nuclei = 210,573) in both sexes across two institutional sites. We identified 32 clusters and found widespread cell type-specific transcriptomic changes across the cortex in AUD, particularly affecting glia. We found the greatest dysregulation in novel microglial and astrocytic subtypes that accounted for the majority of differential gene expression and co-expression modules linked to AUD. Analysis for cell type-specific enrichment of aggregate genetic risk for AUD identified subtypes of microglia and astrocytes as potential key players not only affected by but causally linked to the progression of AUD. These results highlight the importance of cell-type specific molecular changes in AUD and offer opportunities to identify novel targets for treatment.

Regional patterns of human cortex development correlate with underlying neurobiology

Human brain morphology undergoes complex changes over the lifespan. Despite recent progress in tracking brain development via normative models, current knowledge of underlying biological mechanisms is highly limited. We demonstrate that human cortical thickness development and aging trajectories unfold along patterns of molecular and cellular brain organization, traceable from population-level to individual developmental trajectories. During childhood and adolescence, cortex-wide spatial distributions of dopaminergic receptors, inhibitory neurons, glial cell populations, and brain-metabolic features explain up to 50% of variance associated with a lifespan model of regional cortical thickness trajectories. In contrast, modeled cortical thickness change patterns during adulthood are best explained by cholinergic and glutamatergic neurotransmitter receptor and transporter distributions. These relationships are supported by developmental gene expression trajectories and translate to individual longitudinal data from over 8,000 adolescents, explaining up to 59% of developmental change at cohort- and 18% at single-subject level. Integrating neurobiological brain atlases with normative modeling and population neuroimaging provides a biologically meaningful path to understand brain development and aging in living humans.

Histone variant H2BE enhances chromatin accessibility in neurons to promote synaptic gene expression and long-term memory

Histone proteins affect gene expression through multiple mechanisms, including through exchange with histone variants. Recent findings link histone variants to neurological disorders, yet few are well studied in the brain. Most notably, widely expressed variants of H2B remain elusive. We applied recently developed antibodies, biochemical assays, and sequencing approaches to reveal broad expression of the H2B variant H2BE and defined its role in regulating chromatin structure, neuronal transcription, and mouse behavior. We find that H2BE is enriched at promoters, and a single unique amino acid allows it to dramatically enhance chromatin accessibility. Further, we show that H2BE is critical for synaptic gene expression and long-term memory. Together, these data reveal a mechanism linking histone variants to chromatin accessibility, transcriptional regulation, neuronal function, and memory. This work further identifies a widely expressed H2B variant and uncovers a single histone amino acid with profound effects on genomic structure.

Integrative multi-omic analysis reveals conserved cell-projection deficits in human Down syndrome brains

Down syndrome (DS) is the most common genetic cause of cognitive disability. However, it is largely unclear how triplication of a small gene subset may impinge on diverse aspects of DS brain physiopathology. Here, we took a multi-omic approach and simultaneously analyzed by RNA-seq and proteomics the expression signatures of two diverse regions of human postmortem DS brains. We found that the overexpression of triplicated genes triggered global expression dysregulation, differentially affecting transcripts, miRNAs, and proteins involved in both known and novel biological candidate pathways. Among the latter, we observed an alteration in RNA splicing, specifically modulating the expression of genes involved in cytoskeleton and axonal dynamics in DS brains. Accordingly, we found an alteration in axonal polarization in neurons from DS human iPSCs and mice. Thus, our study provides an integrated multilayer expression database capable of identifying new potential targets to aid in designing future clinical interventions for DS.

Molecular Organization of Autonomic, Respiratory, and Spinally-Projecting Neurons in the Mouse Ventrolateral Medulla

The ventrolateral medulla (VLM) is a crucial region in the brain for visceral and somatic control, serving as a significant source of synaptic input to the spinal cord. Experimental studies have shown that gene expression in individual VLM neurons is predictive of their function. However, the molecular and cellular organization of the VLM has remained uncertain. This study aimed to create a comprehensive dataset of VLM cells using single-cell RNA sequencing in male and female mice. The dataset was enriched with targeted sequencing of spinally-projecting and adrenergic/noradrenergic VLM neurons. Based on differentially expressed genes, the resulting dataset of 114,805 VLM cells identifies 23 subtypes of neurons, excluding those in the inferior olive, and five subtypes of astrocytes. Spinally-projecting neurons were found to be abundant in seven subtypes of neurons, which were validated through in situ hybridization. These subtypes included adrenergic/noradrenergic neurons, serotonergic neurons, and neurons expressing gene markers associated with premotor neurons in the ventromedial medulla. Further analysis of adrenergic/noradrenergic neurons and serotonergic neurons identified nine and six subtypes, respectively, within each class of monoaminergic neurons. Marker genes that identify the neural network responsible for breathing were concentrated in two subtypes of neurons, delineated from each other by markers for excitatory and inhibitory neurons. These datasets are available for public download and for analysis with a user-friendly interface. Collectively, this study provides a fine-scale molecular identification of cells in the VLM, forming the foundation for a better understanding of the VLM's role in vital functions and motor control.

Critical reasoning on the co-expression module QTL in the dorsolateral prefrontal cortex

Expression quantitative trait locus (eQTL) analysis is a popular method of gaining insight into the function of regulatory variation. While cis-eQTL resources have been instrumental in linking genome-wide association study variants to gene function, complex trait heritability may be additionally mediated by other forms of gene regulation. Toward this end, novel eQTL methods leverage gene co-expression (module-QTL) to investigate joint regulation of gene modules by single genetic variants. Here we broadly define a "module-QTL" as the association of a genetic variant with a summary measure of gene co-expression. This approach aims to reduce the multiple testing burden of a trans-eQTL search through the consolidation of gene-based testing and provide biological context to eQTLs shared between genes. In this article we provide an in-depth examination of the co-expression module eQTL (module-QTL) through literature review, theoretical investigation, and real-data application of the module-QTL to three large prefrontal cortex genotype-RNA sequencing datasets. We find module-QTLs in our study that are disease associated and reproducible are not additionally informative beyond cis- or trans-eQTLs for module genes. Through comparison to prior studies, we highlight promises and limitations of the module-QTL across study designs and provide recommendations for further investigation of the module-QTL framework.

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.

Traumatic brain injury promotes neurogenesis at the cost of astrogliogenesis in the adult hippocampus of male mice

Traumatic brain injury (TBI) can result in long-lasting changes in hippocampal function. The changes induced by TBI on the hippocampus contribute to cognitive deficits. The adult hippocampus harbors neural stem cells (NSCs) that generate neurons (neurogenesis), and astrocytes (astrogliogenesis). While deregulation of hippocampal NSCs and neurogenesis have been observed after TBI, it is not known how TBI may affect hippocampal astrogliogenesis. Using a controlled cortical impact model of TBI in male mice, single cell RNA sequencing and spatial transcriptomics, we assessed how TBI affected hippocampal NSCs and the neuronal and astroglial lineages derived from them. We observe an increase in NSC-derived neuronal cells and a concomitant decrease in NSC-derived astrocytic cells, together with changes in gene expression and cell dysplasia within the dentate gyrus. Here, we show that TBI modifies NSC fate to promote neurogenesis at the cost of astrogliogenesis and identify specific cell populations as possible targets to counteract TBI-induced cellular changes in the adult hippocampus.

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.

pTDP-43 levels correlate with cell type specific molecular alterations in the prefrontal cortex of C9orf72 ALS/FTD patients

Repeat expansions in the C9orf72 gene are the most common genetic cause of amyotrophic lateral sclerosis and familial frontotemporal dementia (ALS/FTD). To identify molecular defects that take place in the dorsolateral frontal cortex of patients with C9orf72 ALS/FTD, we compared healthy controls with C9orf72 ALS/FTD donor samples staged based on the levels of cortical phosphorylated TAR DNA binding protein (pTDP-43), a neuropathological hallmark of disease progression. We identified distinct molecular changes in different cell types that take place during FTD development. Loss of neurosurveillance microglia and activation of the complement cascade take place early, when pTDP-43 aggregates are absent or very low, and become more pronounced in late stages, suggesting an initial involvement of microglia in disease progression. Reduction of layer 2-3 cortical projection neurons with high expression of CUX2/LAMP5 also occurs early, and the reduction becomes more pronounced as pTDP-43 accumulates. Several unique features were observed only in samples with high levels of pTDP-43, including global alteration of chromatin accessibility in oligodendrocytes, microglia, and astrocytes; higher ratios of premature oligodendrocytes; increased levels of the noncoding RNA NEAT1 in astrocytes and neurons, and higher amount of phosphorylated ribosomal protein S6. Our findings reveal previously unknown progressive functional changes in major cell types found in the frontal cortex of C9orf72 ALS/FTD patients that shed light on the mechanisms underlying the pathology of this disease.

scEM: A New Ensemble Framework for Predicting Cell Type Composition Based on scRNA-Seq Data

With the advent of single-cell RNA sequencing (scRNA-seq) technology, many scRNA-seq data have become available, providing an unprecedented opportunity to explore cellular composition and heterogeneity. Recently, many computational algorithms for predicting cell type composition have been developed, and these methods are typically evaluated on different datasets and performance metrics using diverse techniques. Consequently, the lack of comprehensive and standardized comparative analysis makes it difficult to gain a clear understanding of the strengths and weaknesses of these methods. To address this gap, we reviewed 20 cutting-edge unsupervised cell type identification methods and evaluated these methods comprehensively using 24 real scRNA-seq datasets of varying scales. In addition, we proposed a new ensemble cell-type identification method, named scEM, which learns the consensus similarity matrix by applying the entropy weight method to the four representative methods are selected. The Louvain algorithm is adopted to obtain the final classification of individual cells based on the consensus matrix. Extensive evaluation and comparison with 11 other similarity-based methods under real scRNA-seq datasets demonstrate that the newly developed ensemble algorithm scEM is effective in predicting cellular type composition.

Decoding frontotemporal and cell-type-specific vulnerabilities to neuropsychiatric disorders and psychoactive drugs

Frontotemporal lobe abnormalities are linked to neuropsychiatric disorders and cognition, but the role of cellular heterogeneity between temporal lobe (TL) and frontal lobe (FL) in the vulnerability to genetic risk factors remains to be elucidated. We integrated single-nucleus transcriptome analysis in 'fresh' human FL and TL with genetic susceptibility, gene dysregulation in neuropsychiatric disease and psychoactive drug response data. We show how intrinsic differences between TL and FL contribute to the vulnerability of specific cell types to both genetic risk factors and psychoactive drugs. Neuronal populations, specifically PVALB neurons, were most highly vulnerable to genetic risk factors for psychiatric disease. These psychiatric disease-associated genes were mostly upregulated in the TL, and dysregulated in the brain of patients with obsessive-compulsive disorder, bipolar disorder and schizophrenia. Among these genes, GRIN2A and SLC12A5, implicated in schizophrenia and bipolar disorder, were significantly upregulated in TL PVALB neurons and in psychiatric disease patients' brain. PVALB neurons from the TL were twofold more vulnerable to psychoactive drugs than to genetic risk factors, showing the influence and specificity of frontotemporal lobe differences on cell vulnerabilities. These studies provide a cell type resolved map of the impact of brain regional differences on cell type vulnerabilities in neuropsychiatric disorders.

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.

Advances in long-read single-cell transcriptomics

Long-read single-cell transcriptomics (scRNA-Seq) is revolutionizing the way we profile heterogeneity in disease. Traditional short-read scRNA-Seq methods are limited in their ability to provide complete transcript coverage, resolve isoforms, and identify novel transcripts. The scRNA-Seq protocols developed for long-read sequencing platforms overcome these limitations by enabling the characterization of full-length transcripts. Long-read scRNA-Seq techniques initially suffered from comparatively poor accuracy compared to short read scRNA-Seq. However, with improvements in accuracy, accessibility, and cost efficiency, long-reads are gaining popularity in the field of scRNA-Seq. This review details the advances in long-read scRNA-Seq, with an emphasis on library preparation protocols and downstream bioinformatics analysis tools.

Cell-type specific inference from bulk RNA-sequencing data by integrating single cell reference profiles via EPIC-unmix

Cell type specific (CTS) analysis is essential to reveal biological insights obscured in bulk tissue data. However, single-cell (sc) or single-nuclei (sn) resolution data are still cost-prohibitive for large-scale samples. Thus, computational methods to perform deconvolution from bulk tissue data are highly valuable. We here present EPIC-unmix, a novel two-step empirical Bayesian method integrating reference sc/sn RNA-seq data and bulk RNA-seq data from target samples to enhance the accuracy of CTS inference. We demonstrate through comprehensive simulations across three tissues that EPIC-unmix achieved 4.6% - 109.8% higher accuracy compared to alternative methods. By applying EPIC-unmix to human bulk brain RNA-seq data from the ROSMAP and MSBB cohorts, we identified multiple genes differentially expressed between Alzheimer's disease (AD) cases versus controls in a CTS manner, including 57.4% novel genes not identified using similar sample size sc/snRNA-seq data, indicating the power of our in-silico approach. Among the 6-69% overlapping, 83%-100% are in consistent direction with those from sc/snRNA-seq data, supporting the reliability of our findings. EPIC-unmix inferred CTS expression profiles similarly empowers CTS eQTL analysis. Among the novel eQTLs, we highlight a microglia eQTL for AD risk gene AP3B2, obscured in bulk and missed by sc/snRNA-seq based eQTL analysis. The variant resides in a microglia-specific cCRE, forming chromatin loop with AP3B2 promoter region in microglia. Taken together, we believe EPIC-unmix will be a valuable tool to enable more powerful CTS analysis.

Molecular cascades and cell type-specific signatures in ASD revealed by single-cell genomics

Genomic profiling in postmortem brain from autistic individuals has consistently revealed convergent molecular changes. What drives these changes and how they relate to genetic susceptibility in this complex condition are not well understood. We performed deep single-nucleus RNA sequencing (snRNA-seq) to examine cell composition and transcriptomics, identifying dysregulation of cell type-specific gene regulatory networks (GRNs) in autism spectrum disorder (ASD), which we corroborated using single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) and spatial transcriptomics. Transcriptomic changes were primarily cell type specific, involving multiple cell types, most prominently interhemispheric and callosal-projecting neurons, interneurons within superficial laminae, and distinct glial reactive states involving oligodendrocytes, microglia, and astrocytes. Autism-associated GRN drivers and their targets were enriched in rare and common genetic risk variants, connecting autism genetic susceptibility and cellular and circuit alterations in the human brain.

The impact of chronic pain on brain gene expression

Background

Chronic pain affects one fifth of American adults, contributing significant public health burden. Chronic pain mechanisms can be further understood through investigating brain gene expression.

Methods

We tested differentially expressed genes (DEGs) in chronic pain, migraine, lifetime fentanyl and oxymorphone use, and with chronic pain genetic risk in four brain regions (dACC, DLPFC, MeA, BLA) and imputed cell type expression data from 304 postmortem donors. We compared findings across traits and with independent transcriptomics resources, and performed gene-set enrichment.

Results

We identified two chronic pain DEGs: B4GALT and VEGFB in bulk dACC. We found over 2000 (primarily BLA microglia) chronic pain cell type DEGs. Findings were enriched for mouse microglia pain genes, and for hypoxia and immune response. Cross-trait DEG overlap was minimal.

Conclusions

Chronic pain-associated gene expression is heterogeneous across cell type, largely distinct from that in pain-related traits, and shows BLA microglia are a key cell type.

Advanced sequencing-based high-throughput and long-read single-cell transcriptome analysis

Cells are the fundamental building blocks of living systems, exhibiting significant heterogeneity. The transcriptome connects the cellular genotype and phenotype, and profiling single-cell transcriptomes is critical for uncovering distinct cell types, states, and the interplay between cells in development, health, and disease. Nevertheless, single-cell transcriptome analysis faces daunting challenges due to the low abundance and diverse nature of RNAs in individual cells, as well as their heterogeneous expression. The advent and continuous advancements of next-generation sequencing (NGS) and third-generation sequencing (TGS) technologies have solved these problems and facilitated the high-throughput, sensitive, full-length, and rapid profiling of single-cell RNAs. In this review, we provide a broad introduction to current methodologies for single-cell transcriptome sequencing. First, state-of-the-art advancements in high-throughput and full-length single-cell RNA sequencing (scRNA-seq) platforms using NGS are reviewed. Next, TGS-based long-read scRNA-seq methods are summarized. Finally, a brief conclusion and perspectives for comprehensive single-cell transcriptome analysis are discussed.

Characterizing dysregulations via cell-cell communications in Alzheimer's brains using single-cell transcriptomes

BACKGROUND

Alzheimer's disease (AD) is a devastating neurodegenerative disorder affecting 44 million people worldwide, leading to cognitive decline, memory loss, and significant impairment in daily functioning. The recent single-cell sequencing technology has revolutionized genetic and genomic resolution by enabling scientists to explore the diversity of gene expression patterns at the finest resolution. Most existing studies have solely focused on molecular perturbations within each cell, but cells live in microenvironments rather than in isolated entities. Here, we leveraged the large-scale and publicly available single-nucleus RNA sequencing in the human prefrontal cortex to investigate cell-to-cell communication in healthy brains and their perturbations in AD. We uniformly processed the snRNA-seq with strict QCs and labeled canonical cell types consistent with the definitions from the BRAIN Initiative Cell Census Network. From ligand and receptor gene expression, we built a high-confidence cell-to-cell communication network to investigate signaling differences between AD and healthy brains.

RESULTS

Specifically, we first performed broad communication pattern analyses to highlight that biologically related cell types in normal brains rely on largely overlapping signaling networks and that the AD brain exhibits the irregular inter-mixing of cell types and signaling pathways. Secondly, we performed a more focused cell-type-centric analysis and found that excitatory neurons in AD have significantly increased their communications to inhibitory neurons, while inhibitory neurons and other non-neuronal cells globally decreased theirs to all cells. Then, we delved deeper with a signaling-centric view, showing that canonical signaling pathways CSF, TGFβ, and CX3C are significantly dysregulated in their signaling to the cell type microglia/PVM and from endothelial to neuronal cells for the WNT pathway. Finally, after extracting 23 known AD risk genes, our intracellular communication analysis revealed a strong connection of extracellular ligand genes APP, APOE, and PSEN1 to intracellular AD risk genes TREM2, ABCA1, and APP in the communication from astrocytes and microglia to neurons.

CONCLUSIONS

In summary, with the novel advances in single-cell sequencing technologies, we show that cellular signaling is regulated in a cell-type-specific manner and that improper regulation of extracellular signaling genes is linked to intracellular risk genes, giving the mechanistic intra- and inter-cellular picture of AD.

Characterizing cell type specific transcriptional differences between the living and postmortem human brain

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

Molecular and cellular mechanisms of selective vulnerability in neurodegenerative diseases

The selective vulnerability of specific neuronal subtypes is a hallmark of neurodegenerative diseases. In this Review, I summarize our current understanding of the brain regions and cell types that are selectively vulnerable in different neurodegenerative diseases and describe the proposed underlying cell-autonomous and non-cell-autonomous mechanisms. I highlight how recent methodological innovations - including single-cell transcriptomics, CRISPR-based screens and human cell-based models of disease - are enabling new breakthroughs in our understanding of selective vulnerability. An understanding of the molecular mechanisms that determine selective vulnerability and resilience would shed light on the key processes that drive neurodegeneration and point to potential therapeutic strategies to protect vulnerable cell populations.

NKX2-2 based nuclei sorting on frozen human archival pancreas enables the enrichment of islet endocrine populations for single-nucleus RNA sequencing

BACKGROUND

Current approaches to profile the single-cell transcriptomics of human pancreatic endocrine cells almost exclusively rely on freshly isolated islets. However, human islets are limited in availability. Furthermore, the extensive processing steps during islet isolation and subsequent single cell dissolution might alter gene expressions. In this work, we report the development of a single-nucleus RNA sequencing (snRNA-seq) approach with targeted islet cell enrichment for endocrine-population focused transcriptomic profiling using frozen archival pancreatic tissues without islet isolation.

RESULTS

We cross-compared five nuclei isolation protocols and selected the citric acid method as the best strategy to isolate nuclei with high RNA integrity and low cytoplasmic contamination from frozen archival human pancreata. We innovated fluorescence-activated nuclei sorting based on the positive signal of NKX2-2 antibody to enrich nuclei of the endocrine population from the entire nuclei pool of the pancreas. Our sample preparation procedure generated high-quality single-nucleus gene-expression libraries while preserving the endocrine population diversity. In comparison with single-cell RNA sequencing (scRNA-seq) library generated with live cells from freshly isolated human islets, the snRNA-seq library displayed comparable endocrine cellular composition and cell type signature gene expression. However, between these two types of libraries, differential enrichments of transcripts belonging to different functional classes could be observed.

CONCLUSIONS

Our work fills a technological gap and helps to unleash frozen archival pancreatic tissues for molecular profiling targeting the endocrine population. This study opens doors to retrospective mappings of endocrine cell dynamics in pancreatic tissues of complex histopathology. We expect that our protocol is applicable to enrich nuclei for transcriptomics studies from various populations in different types of frozen archival tissues.

Mitochondrial complex I deficiency stratifies idiopathic Parkinson's disease

Idiopathic Parkinson's disease (iPD) is believed to have a heterogeneous pathophysiology, but molecular disease subtypes have not been identified. Here, we show that iPD can be stratified according to the severity of neuronal respiratory complex I (CI) deficiency, and identify two emerging disease subtypes with distinct molecular and clinical profiles. The CI deficient (CI-PD) subtype accounts for approximately a fourth of all cases, and is characterized by anatomically widespread neuronal CI deficiency, a distinct cell type-specific gene expression profile, increased load of neuronal mtDNA deletions, and a predilection for non-tremor dominant motor phenotypes. In contrast, the non-CI deficient (nCI-PD) subtype exhibits no evidence of mitochondrial impairment outside the dopaminergic substantia nigra and has a predilection for a tremor dominant phenotype. These findings constitute a step towards resolving the biological heterogeneity of iPD with implications for both mechanistic understanding and treatment strategies.

Antipsychotic-induced epigenomic reorganization in frontal cortex of individuals with schizophrenia

Genome-wide association studies have revealed >270 loci associated with schizophrenia risk, yet these genetic factors do not seem to be sufficient to fully explain the molecular determinants behind this psychiatric condition. Epigenetic marks such as post-translational histone modifications remain largely plastic during development and adulthood, allowing a dynamic impact of environmental factors, including antipsychotic medications, on access to genes and regulatory elements. However, few studies so far have profiled cell-specific genome-wide histone modifications in postmortem brain samples from schizophrenia subjects, or the effect of antipsychotic treatment on such epigenetic marks. Here, we conducted ChIP-seq analyses focusing on histone marks indicative of active enhancers (H3K27ac) and active promoters (H3K4me3), alongside RNA-seq, using frontal cortex samples from antipsychotic-free (AF) and antipsychotic-treated (AT) individuals with schizophrenia, as well as individually matched controls (n=58). Schizophrenia subjects exhibited thousands of neuronal and non-neuronal epigenetic differences at regions that included several susceptibility genetic loci, such as NRG1, DISC1, and DRD3. By analyzing the AF and AT cohorts separately, we identified schizophrenia-associated alterations in specific transcription factors, their regulatees, and epigenomic and transcriptomic features that were reversed by antipsychotic treatment; as well as those that represented a consequence of antipsychotic medication rather than a hallmark of schizophrenia in postmortem human brain samples. Notably, we also found that the effect of age on epigenomic landscapes was more pronounced in frontal cortex of AT-schizophrenics, as compared to AF-schizophrenics and controls. Together, these data provide important evidence of epigenetic alterations in the frontal cortex of individuals with schizophrenia, and remark for the first time on the impact of age and antipsychotic treatment on chromatin organization.

Dissecting the spatiotemporal diversity of adult neural stem cells

Adult stem cells are important for tissue turnover and regeneration. However, in most adult systems it remains elusive how stem cells assume different functional states and support spatially patterned tissue architecture. Here, we dissected the diversity of neural stem cells in the adult zebrafish brain, an organ that is characterized by pronounced zonation and high regenerative capacity. We combined single-cell transcriptomics of dissected brain regions with massively parallel lineage tracing and in vivo RNA metabolic labeling to analyze the regulation of neural stem cells in space and time. We detected a large diversity of neural stem cells, with some subtypes being restricted to a single brain region, while others were found globally across the brain. Global stem cell states are linked to neurogenic differentiation, with different states being involved in proliferative and non-proliferative differentiation. Our work reveals principles of adult stem cell organization and establishes a resource for the functional manipulation of neural stem cell subtypes.

A latent clinical-anatomical dimension relating metabolic syndrome to brain structure and cognition

The link between metabolic syndrome (MetS) and neurodegenerative as well as cerebrovascular conditions holds substantial implications for brain health in at-risk populations. This study elucidates the complex relationship between MetS and brain health by conducting a comprehensive examination of cardiometabolic risk factors, brain morphology, and cognitive function in 40,087 individuals. Multivariate, data-driven statistics identified a latent dimension linking more severe MetS to widespread brain morphological abnormalities, accounting for up to 71% of shared variance in the data. This dimension was replicable across sub-samples. In a mediation analysis, we could demonstrate that MetS-related brain morphological abnormalities mediated the link between MetS severity and cognitive performance in multiple domains. Employing imaging transcriptomics and connectomics, our results also suggest that MetS-related morphological abnormalities are linked to the regional cellular composition and macroscopic brain network organization. By leveraging extensive, multi-domain data combined with a dimensional stratification approach, our analysis provides profound insights into the association of MetS and brain health. These findings can inform effective therapeutic and risk mitigation strategies aimed at maintaining brain integrity.

Microfluidic platform for omics analysis on single cells with diverse morphology and size: A review

BACKGROUND

Microfluidic techniques have emerged as powerful tools in single-cell research, facilitating the exploration of omics information from individual cells. Cell morphology is crucial for gene expression and physiological processes. However, there is currently a lack of integrated analysis of morphology and single-cell omics information. A critical challenge remains: what platform technologies are the best option to decode omics data of cells that are complex in morphology and size?

RESULTS

This review highlights achievements in microfluidic-based single-cell omics and isolation of cells based on morphology, along with other cell sorting methods based on physical characteristics. Various microfluidic platforms for single-cell isolation are systematically presented, showcasing their diversity and adaptability. The discussion focuses on microfluidic devices tailored to the distinct single-cell isolation requirements in plants and animals, emphasizing the significance of considering cell morphology and cell size in optimizing single-cell omics strategies. Simultaneously, it explores the application of microfluidic single-cell sorting technologies to single-cell sequencing, aiming to effectively integrate information about cell shape and size.

SIGNIFICANCE AND NOVELTY

The novelty lies in presenting a comprehensive overview of recent accomplishments in microfluidic-based single-cell omics, emphasizing the integration of different microfluidic platforms and their implications for cell morphology-based isolation. By underscoring the pivotal role of the specialized morphology of different cells in single-cell research, this review provides robust support for delving deeper into the exploration of single-cell omics data.

Transcriptional cartography integrates multiscale biology of the human cortex

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Spatially Self-Organized Three-Dimensional Neural Concentroid as a Novel Reductionist Humanized Model to Study Neurovascular Development

Although human pluripotent stem cell (PSC)-derived brain organoids have enabled researchers to gain insight into human brain development and disease, these organoids contain solely ectodermal cells and are not vascularized as occurs during brain development. Here it is created less complex and more homogenous large neural constructs starting from PSC-derived neuroprogenitor cells (NPC), by fusing small NPC spheroids into so-called concentroids. Such concentroids consisted of a pro-angiogenic core, containing neuronal and outer radial glia cells, surrounded by an astroglia-dense outer layer. Incorporating PSC-derived endothelial cells (EC) around and/or in the concentroids promoted vascularization, accompanied by differential outgrowth and differentiation of neuronal and astroglia cells, as well as the development of ectodermal-derived pericyte-like mural cells co-localizing with EC networks. Single nucleus transcriptomic analysis revealed an enhanced neural cell subtype maturation and diversity in EC-containing concentroids, which better resemble the fetal human brain compared to classical organoids or NPC-only concentroids. This PSC-derived "vascularized" concentroid brain model will facilitate the study of neurovascular/blood-brain barrier development, neural cell migration, and the development of effective in vitro vascularization strategies of brain mimics.

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.

Single-nucleus transcriptomic analysis reveals the relationship between gene expression in oligodendrocyte lineage and major depressive disorder

BACKGROUND

Major depressive disorder (MDD) is a common mental illness that affects millions of people worldwide and imposes a heavy burden on individuals, families and society. Previous studies on MDD predominantly focused on neurons and employed bulk homogenates of brain tissues. This paper aims to decipher the relationship between oligodendrocyte lineage (OL) development and MDD at the single-cell resolution level.

METHODS

Here, we present the use of a guided regularized random forest (GRRF) algorithm to explore single-nucleus RNA sequencing profiles (GSE144136) of the OL at four developmental stages, which contains dorsolateral prefrontal cortex of 17 healthy controls (HC) and 17 MDD cases, generated by Nagy C et al. We prioritized and ordered differentially expressed genes (DEGs) based on Nagy et al., which could predominantly discriminate cells in the four developmental stages and two adjacent developmental stages of the OL. We further screened top-ranked genes that distinguished between HC and MDD in four developmental stages. Moreover, we estimated the performance of the GRRF model via the area under the curve value. Additionally, we validated the pivotal candidate gene Malat1 in animal models.

RESULTS

We found that, among the four developmental stages, the onset development of OL (OPC2) possesses the best predictive power for distinguishing HC and MDD, and long noncoding RNA MALAT1 has top-ranked importance value in candidate genes of four developmental stages. In addition, results of fluorescence in situ hybridization assay showed that Malat1 plays a critical role in the occurrence of depression.

CONCLUSIONS

Our work elucidates the mechanism of MDD from the perspective of OL development at the single-cell resolution level and provides novel insight into the occurrence of depression.

Quantifying the proportion of different cell types in the human cortex using DNA methylation profiles

BACKGROUND

Due to interindividual variation in the cellular composition of the human cortex, it is essential that covariates that capture these differences are included in epigenome-wide association studies using bulk tissue. As experimentally derived cell counts are often unavailable, computational solutions have been adopted to estimate the proportion of different cell types using DNA methylation data. Here, we validate and profile the use of an expanded reference DNA methylation dataset incorporating two neuronal and three glial cell subtypes for quantifying the cellular composition of the human cortex.

RESULTS

We tested eight reference panels containing different combinations of neuronal- and glial cell types and characterised their performance in deconvoluting cell proportions from computationally reconstructed or empirically derived human cortex DNA methylation data. Our analyses demonstrate that while these novel brain deconvolution models produce accurate estimates of cellular proportions from profiles generated on postnatal human cortex samples, they are not appropriate for the use in prenatal cortex or cerebellum tissue samples. Applying our models to an extensive collection of empirical datasets, we show that glial cells are twice as abundant as neuronal cells in the human cortex and identify significant associations between increased Alzheimer's disease neuropathology and the proportion of specific cell types including a decrease in NeuNNeg/SOX10Neg nuclei and an increase of NeuNNeg/SOX10Pos nuclei.

CONCLUSIONS

Our novel deconvolution models produce accurate estimates for cell proportions in the human cortex. These models are available as a resource to the community enabling the control of cellular heterogeneity in epigenetic studies of brain disorders performed on bulk cortex tissue.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Decreased CNNM2 expression in prefrontal cortex affects sensorimotor gating function, cognition, dendritic spine morphogenesis and risk of schizophrenia

Genome-wide association studies (GWASs) have reported multiple single nucleotide polymorphisms (SNPs) associated with schizophrenia, yet the underlying molecular mechanisms are largely unknown. In this study, we aimed to identify schizophrenia relevant genes showing alterations in mRNA and protein expression associated with risk SNPs at the 10q24.32-33 GWAS locus. We carried out the quantitative trait loci (QTL) and summary data-based Mendelian randomization (SMR) analyses, using the PsychENCODE dorsolateral prefrontal cortex (DLPFC) expression QTL (eQTL) database, as well as the ROSMAP and Banner DLPFC protein QTL (pQTL) datasets. The gene CNNM2 (encoding a magnesium transporter) at 10q24.32-33 was identified to be a robust schizophrenia risk gene, and was highly expressed in human neurons according to single cell RNA-seq (scRNA-seq) data. We further revealed that reduced Cnnm2 in the mPFC of mice led to impaired cognition and compromised sensorimotor gating function, and decreased Cnnm2 in primary cortical neurons altered dendritic spine morphogenesis, confirming the link between CNNM2 and endophenotypes of schizophrenia. Proteomics analyses showed that reduced Cnnm2 level changed expression of proteins associated with neuronal structure and function. Together, these results identify a robust gene in the pathogenesis of schizophrenia.

Downregulation of Nrp1 transcription promotes blood-brain barrier disruption following experimental cerebral ischemia-reperfusion

Disruption of the blood-brain barrier (BBB) following cerebral ischemia-reperfusion injury (CIRI) is a major factor in the pathophysiology of stroke. Endothelial cell-cell communication is essential for maintaining BBB integrity. By analyzing GSE227651 data, we found that a decrease in endothelial cell-cell communication mediated by Sema3/Nrp1 may be due to the downregulation of Nrp1 transcription, which could contribute to BBB breakdown after CIRI. We confirmed this hypothesis by using western blot analysis to show a reduction in Nrp1 protein levels in penumbra endothelial cells after CIRI in mice. We then overexpressed Nrp1 specifically in brain endothelial cells using adeno-associated virus in mice. Furthermore, Nrp1 overexpression had a protective effect on BBB integrity, as evidenced by a decrease in IgG and albumin leakage caused by CIRI in mice. Finally, we found that Nrp1 overexpression also reduced brain cell death and neurological deficits induced by cerebral ischemia-reperfusion in mice. Our findings suggest that Nrp1 downregulation may be a key factor in the breakdown of endothelial cell-cell communication and subsequent BBB disruption following CIRI. Targeting Nrp1-mediated pathways may be a promising approach for mitigating BBB damage and alleviating neurological consequences in stroke patients.

Neuronal Spike Shapes (NSS): A straightforward approach to investigate heterogeneity in neuronal excitability states

The mammalian brain exhibits a remarkable diversity of neurons, contributing to its intricate architecture and functional complexity. The analysis of multimodal single-cell datasets enables the investigation of cell types and states heterogeneity. In this study, we introduce the Neuronal Spike Shapes (NSS), a straightforward approach for the exploration of excitability states of neurons based on their Action Potential (AP) waveforms. The NSS method describes the AP waveform based on a triangular representation complemented by a set of derived electrophysiological (EP) features. To support this hypothesis, we validate the proposed approach on two datasets of murine cortical neurons, focusing it on GABAergic neurons. The validation process involves a combination of NSS-based clustering analysis, features exploration, Differential Expression (DE), and Gene Ontology (GO) enrichment analysis. Results show that the NSS-based analysis captures neuronal excitability states that possess biological relevance independently of cell subtype. In particular, Neuronal Spike Shapes (NSS) captures, among others, a well-characterized fast-spiking excitability state, supported by both electrophysiological and transcriptomic validation. Gene Ontology Enrichment Analysis reveals voltage-gated potassium (K+) channels as specific markers of the identified NSS partitions. This finding strongly corroborates the biological relevance of NSS partitions as excitability states, as the expression of voltage-gated K+ channels regulates the hyperpolarization phase of the AP, being directly implicated in the regulation of neuronal excitability.

The use of single-cell RNA-seq to study heterogeneity at varying levels of virus-host interactions

The outcome of viral infection depends on the diversity of the infecting viral population and the heterogeneity of the cell population that is infected. Until almost a decade ago, the study of these dynamic processes during viral infection was challenging and limited to certain targeted measurements. Presently, with the use of single-cell sequencing technology, the complex interface defined by the interactions of cells with infecting virus can now be studied across the breadth of the transcriptome in thousands of individual cells simultaneously. In this review, we will describe the use of single-cell RNA sequencing (scRNA-seq) to study the heterogeneity of viral infections, ranging from individual virions to the immune response between infected individuals. In addition, we highlight certain key experimental limitations and methodological decisions that are critical to analyzing scRNA-seq data at each scale.

Common genetic variation impacts stress response in the brain

To explain why individuals exposed to identical stressors experience divergent clinical outcomes, we determine how molecular encoding of stress modifies genetic risk for brain disorders. Analysis of post-mortem brain (n=304) revealed 8557 stress-interactive expression quantitative trait loci (eQTLs) that dysregulate expression of 915 eGenes in response to stress, and lie in stress-related transcription factor binding sites. Response to stress is robust across experimental paradigms: up to 50% of stress-interactive eGenes validate in glucocorticoid treated hiPSC-derived neurons (n=39 donors). Stress-interactive eGenes show brain region- and cell type-specificity, and, in post-mortem brain, implicate glial and endothelial mechanisms. Stress dysregulates long-term expression of disorder risk genes in a genotype-dependent manner; stress-interactive transcriptomic imputation uncovered 139 novel genes conferring brain disorder risk only in the context of traumatic stress. Molecular stress-encoding explains individualized responses to traumatic stress; incorporating trauma into genomic studies of brain disorders is likely to improve diagnosis, prognosis, and drug discovery.

FTY720 requires vitamin B12-TCN2-CD320 signaling in astrocytes to reduce disease in an animal model of multiple sclerosis

Vitamin B12 (B12) deficiency causes neurological manifestations resembling multiple sclerosis (MS); however, a molecular explanation for the similarity is unknown. FTY720 (fingolimod) is a sphingosine 1-phosphate (S1P) receptor modulator and sphingosine analog approved for MS therapy that can functionally antagonize S1P1. Here, we report that FTY720 suppresses neuroinflammation by functionally and physically regulating the B12 pathways. Genetic and pharmacological S1P1 inhibition upregulates a transcobalamin 2 (TCN2)-B12 receptor, CD320, in immediate-early astrocytes (ieAstrocytes; a c-Fos-activated astrocyte subset that tracks with experimental autoimmune encephalomyelitis [EAE] severity). CD320 is also reduced in MS plaques. Deficiency of CD320 or dietary B12 restriction worsens EAE and eliminates FTY720's efficacy while concomitantly downregulating type I interferon signaling. TCN2 functions as a chaperone for FTY720 and sphingosine, whose complex induces astrocytic CD320 internalization, suggesting a delivery mechanism of FTY720/sphingosine via the TCN2-CD320 pathway. Taken together, the B12-TCN2-CD320 pathway is essential for the mechanism of action of FTY720.

A latent clinical-anatomical dimension relating metabolic syndrome to brain structure and cognition

The link between metabolic syndrome (MetS) and neurodegenerative as well cerebrovascular conditions holds substantial implications for brain health in at-risk populations. This study elucidates the complex relationship between MetS and brain health by conducting a comprehensive examination of cardiometabolic risk factors, cortical morphology, and cognitive function in 40,087 individuals. Multivariate, data-driven statistics identified a latent dimension linking more severe MetS to widespread brain morphological abnormalities, accounting for up to 71% of shared variance in the data. This dimension was replicable across sub-samples. In a mediation analysis we could demonstrate that MetS-related brain morphological abnormalities mediated the link between MetS severity and cognitive performance in multiple domains. Employing imaging transcriptomics and connectomics, our results also suggest that MetS-related morphological abnormalities are linked to the regional cellular composition and macroscopic brain network organization. By leveraging extensive, multi-domain data combined with a dimensional stratification approach, our analysis provides profound insights into the association of MetS and brain health. These findings can inform effective therapeutic and risk mitigation strategies aimed at maintaining brain integrity.