A single dose of cocaine rewires the 3D genome structure of midbrain dopamine neurons

Midbrain dopamine neurons (DNs) respond to a first exposure to addictive drugs and play key roles in chronic drug usage 1-3 . As the synaptic and transcriptional changes that follow an acute cocaine exposure are mostly resolved within a few days 4,5 , the molecular changes that encode the long-term cellular memory of the exposure within DNs remain unknown. To investigate whether a single cocaine exposure induces long-term changes in the 3D genome structure of DNs, we applied Genome Architecture Mapping and single nucleus transcriptomic analyses in the mouse midbrain. We found extensive rewiring of 3D genome architecture at 24 hours past exposure which remains or worsens by 14 days, outlasting transcriptional responses. The cocaine-induced chromatin rewiring occurs at all genomic scales and affects genes with major roles in cocaine-induced synaptic changes. A single cocaine exposure triggers extensive long-lasting changes in chromatin condensation in post-synaptic and post-transcriptional regulatory genes, for example the unfolding of Rbfox1 which becomes most prominent 14 days post exposure. Finally, structurally remodeled genes are most expressed in a specific DN sub-type characterized by low expression of the dopamine auto-receptor Drd2 , a key feature of highly cocaine-sensitive cells. These results reveal an important role for long-lasting 3D genome remodelling in the cellular memory of a single cocaine exposure, providing new hypotheses for understanding the inception of drug addiction and 3D genome plasticity.

Heterogeneity of layer 4 in visual areas of rhesus macaque cortex

Recently, single-cell RNA-sequencing (scRNA-seq) has enabled unprecedented insights to the cellular landscape of the brains of many different species, among them the rhesus macaque as a key animal model. Building on previous, broader surveys of the macaque brain, we closely examined five immediately neighboring areas within the visual cortex of the rhesus macaque: V1, V2, V4, MT and TEO. To facilitate this, we first devised a novel pipeline for brain spatial archive - the BrainSPACE - which enabled robust archiving and sampling from the whole unfixed brain. SnRNA-sequencing of ~100,000 nuclei from visual areas V1 and V4 revealed conservation within the GABAergic neuron subtypes, while seven and one distinct principle neuron subtypes were detected in V1 and V4, respectively, all most likely located in layer 4. Moreover, using small molecule fluorescence in situ hybridization, we identified cell type density gradients across V1, V2, V4, MT, and TEO appearing to reflect the visual hierarchy. These findings demonstrate an association between the clear areal specializations among neighboring areas with the hierarchical levels within the visual cortex of the rhesus macaque.

A maturational shift in the frontal cortex synaptic transcriptional landscape underlies schizophrenia-relevant behavioural traits: A congenital rat model

Disruption of brain development early in life may underlie the neurobiology behind schizophrenia. We have reported more immature synaptic spines in the frontal cortex (FC) of adult Roman High-Avoidance (RHA-I) rats, a behavioural model displaying schizophrenia-like traits. Here, we performed a whole transcriptome analysis in the FC of 4 months old male RHA-I (n=8) and its counterpart, the Roman Low-Avoidance (RLA-I) (n=8). We identified 203 significant genes with overrepresentation of genes involved in synaptic function. Next, we performed a gene set enrichment analysis (GSEA) for genes co-expressed during neurodevelopment. Gene networks were obtained by weighted gene co-expression network analysis (WGCNA) of a transcriptomic dataset containing human FC during lifespan (n=269). Out of thirty-one functional gene networks, six were significantly enriched in the RHA-I. These were differentially regulated during infancy and enriched in biological ontologies related to myelination, synaptic function, and immune response. We validated differential gene expression in a new cohort of adolescent (<=2 months old) and young-adult (>=3 months old) RHA-I and RLA-I rats. The results confirmed overexpression of Gsn, Nt5cd1, Ppp1r1b, and Slc9a3r1 in young-adult RHA-I, while Cables1, a regulator of Cdk5 phosphorylation in actin regulation and involved in synaptic plasticity and maturation, was significantly downregulated in adolescent RHA-I. This age-related expression change was also observed for presynaptic components Snap25 and Snap29. Our results show a different maturational expression profile of synaptic components in the RHA-I strain, supporting a shift in FC maturation underlying schizophrenia-like behavioural traits and adding construct validity to this strain as a neurodevelopmental model.

A New Technical Approach for Cross-species Examination of Neuronal Wiring and Adult Neuron-glia Functions

Advances in single cell sequencing have enabled the identification of a large number of genes, expressed in many different cell types, and across a variety of model organisms. In particular, the nervous system harbors an immense number of interacting cell types, which are poorly characterized. Future loss- and gain-of-function experiments will be essential in determining how novel genes play critical roles in diverse cellular, as well as evolutionarily adapted, contexts. However, functional analysis across species is often hampered by technical limitations, in non-genetic animal systems. Here, we describe a new single plasmid system, misPiggy. The system is based around the hyperactive piggyBac transposon system, which combines stable genomic integration of transgenes (for long-term expression) with large cargo capacity. Taking full advantage of these characteristics, we engineered novel expression modules into misPiggy that allow for cell-type specific loss- and gain-of-gene function. These modules work widely across species from frog to ferret. As a proof of principle, we present a loss-of-function analysis of the neuronal receptor Deleted in Colorectal Cancer (DCC) in retinal ganglion cells (RGCs) of Xenopus tropicalis tadpoles. Single axon tracings of mosaic knock-out cells reveal a specific cell-intrinsic requirement of DCC, specifically in axonal arborization within the frog tectum, rather than retina-to-brain axon guidance. Furthermore, we report additional technical advances that enable temporal control of knock-down or gain-of-function analysis. We applied this to visualize and manipulate labeled neurons, astrocytes and other glial cells in the central nervous system (CNS) of mouse, rat and ferret. We propose that misPiggy will be a valuable tool for rapid, flexible and cost-effective screening of gene function across a variety of animal models.

Upper cortical layer-driven network impairment in schizophrenia

Schizophrenia is one of the most widespread and complex mental disorders. To characterize the impact of schizophrenia, we performed single-nucleus RNA sequencing (snRNA-seq) of >220,000 neurons from the dorsolateral prefrontal cortex of patients with schizophrenia and matched controls. In addition, >115,000 neurons were analyzed topographically by immunohistochemistry. Compositional analysis of snRNA-seq data revealed a reduction in abundance of GABAergic neurons and a concomitant increase in principal neurons, most pronounced for upper cortical layer subtypes, which was substantiated by histological analysis. Many neuronal subtypes showed extensive transcriptomic changes, the most marked in upper-layer GABAergic neurons, including down-regulation in energy metabolism and up-regulation in neurotransmission. Transcription factor network analysis demonstrated a developmental origin of transcriptomic changes. Last, Visium spatial transcriptomics further corroborated upper-layer neuron vulnerability in schizophrenia. Overall, our results point toward general network impairment within upper cortical layers as a core substrate associated with schizophrenia symptomatology.

Identification of epilepsy-associated neuronal subtypes and gene expression underlying epileptogenesis

Epilepsy is one of the most common neurological disorders, yet its pathophysiology is poorly understood due to the high complexity of affected neuronal circuits. To identify dysfunctional neuronal subtypes underlying seizure activity in the human brain, we have performed single-nucleus transcriptomics analysis of >110,000 neuronal transcriptomes derived from temporal cortex samples of multiple temporal lobe epilepsy and non-epileptic subjects. We found that the largest transcriptomic changes occur in distinct neuronal subtypes from several families of principal neurons (L5-6_Fezf2 and L2-3_Cux2) and GABAergic interneurons (Sst and Pvalb), whereas other subtypes in the same families were less affected. Furthermore, the subtypes with the largest epilepsy-related transcriptomic changes may belong to the same circuit, since we observed coordinated transcriptomic shifts across these subtypes. Glutamate signaling exhibited one of the strongest dysregulations in epilepsy, highlighted by layer-wise transcriptional changes in multiple glutamate receptor genes and strong upregulation of genes coding for AMPA receptor auxiliary subunits. Overall, our data reveal a neuronal subtype-specific molecular phenotype of epilepsy.

The pathophysiology of epilepsy is unclear. Here, the authors present single-nuclei transcriptomic profiling of human temporal lobe epilepsy from patients. They identified epilepsy-associated neuronal subtypes, and a panel of dysregulated genes, predicting neuronal circuits contributing to epilepsy.

Astrocyte layers in the mammalian cerebral cortex revealed by a single-cell in situ transcriptomic map

Although the cerebral cortex is organized into six excitatory neuronal layers, it is unclear whether glial cells show distinct layering. In the present study, we developed a high-content pipeline, the large-area spatial transcriptomic (LaST) map, which can quantify single-cell gene expression in situ. Screening 46 candidate genes for astrocyte diversity across the mouse cortex, we identified superficial, mid and deep astrocyte identities in gradient layer patterns that were distinct from those of neurons. Astrocyte layer features, established in the early postnatal cortex, mostly persisted in adult mouse and human cortex. Single-cell RNA sequencing and spatial reconstruction analysis further confirmed the presence of astrocyte layers in the adult cortex. Satb2 and Reeler mutations that shifted neuronal post-mitotic development were sufficient to alter glial layering, indicating an instructive role for neuronal cues. Finally, astrocyte layer patterns diverged between mouse cortical regions. These findings indicate that excitatory neurons and astrocytes are organized into distinct lineage-associated laminae.

Identification of region-specific astrocyte subtypes at single cell resolution

Astrocytes, a major cell type found throughout the central nervous system, have general roles in the modulation of synapse formation and synaptic transmission, blood-brain barrier formation, and regulation of blood flow, as well as metabolic support of other brain resident cells. Crucially, emerging evidence shows specific adaptations and astrocyte-encoded functions in regions, such as the spinal cord and cerebellum. To investigate the true extent of astrocyte molecular diversity across forebrain regions, we used single-cell RNA sequencing. Our analysis identifies five transcriptomically distinct astrocyte subtypes in adult mouse cortex and hippocampus. Validation of our data in situ reveals distinct spatial positioning of defined subtypes, reflecting the distribution of morphologically and physiologically distinct astrocyte populations. Our findings are evidence for specialized astrocyte subtypes between and within brain regions. The data are available through an online database (https://holt-sc.glialab.org/), providing a resource on which to base explorations of local astrocyte diversity and function in the brain.

An immunoaffinity-based method for isolating ultrapure adult astrocytes based on ATP1B2 targeting by the ACSA-2 antibody

Astrocytes are a major cell type in the mammalian CNS. Astrocytes are now known to play a number of essential roles in processes including synapse formation and function, as well as blood-brain barrier formation and control of cerebral blood flow. However, our understanding of the molecular mechanisms underlying astrocyte development and function is still rudimentary. This lack of knowledge is at least partly due to the lack of tools currently available for astrocyte biology. ACSA-2 is a commercially available antibody originally developed for the isolation of astrocytes from young postnatal mouse brain, using magnetic cell-sorting methods, but its utility in isolating cells from adult tissue has not yet been published. Using a modified protocol, we now show that this tool can also be used to isolate ultrapure astrocytes from the adult brain. Furthermore, using a variety of techniques (including single-cell sequencing, overexpression and knockdown assays, immunoblotting, and immunohistochemistry), we identify the ACSA-2 epitope for the first time as ATP1B2 and characterize its distribution in the CNS. Finally, we show that ATP1B2 is stably expressed in multiple models of CNS injury and disease. Hence, we show that the ACSA-2 antibody possesses the potential to be an extremely valuable tool for astrocyte research, allowing the purification and characterization of astrocytes (potentially including injury and disease models) without the need for any specialized and expensive equipment. In fact, our results suggest that ACSA-2 should be a first-choice method for astrocyte isolation and characterization.

Compensatory reduction of Cav3.1 expression in thalamocortical neurons of juvenile rats of WAG/Rij model of absence epilepsy

Absence seizures are the non-convulsive form of generalized epilepsy critically dependent on T-type calcium channels (Cav3) in thalamic neurons. In humans, absences accompany only childhood or adolescent epileptic syndromes--though in its polygenic rat models WAG/Rij and GAERS the opposite developmental pattern is observed. Hereby we address this issue by transcriptional and functional study of thalamic Cav3 in juvenile (i.e., free of seizures) rats of the absence-prone WAG/Rij strain and their coevals of the maternal Wistar strain. First, we measured the low voltage-activated (LVA) Ca(2+) current in freshly isolated thalamocortical neurons from laterodorsal nucleus of thalamus. The difference between current densities in control (12.9 ± 1.8pA/pF) and absence epilepsy (7.9 ± 1.8pA/pF) groups reached ∼ 39%. Second, we assessed the contribution of different T-channel isoforms into the reduction of Cav3-mediated current in WAG/Rij juveniles by means of RT PCR. The expression of all three LVA calcium channels was revealed with the prevalence of G and I isoforms. The expression level of G isoform (Cav3.1) was 35% smaller in WAG/Rij strain if compared to the control animals while that of H and I isoforms (Cav3.2 and Cav3.3, respectively) remained stable. The weakened expression of Cav3.1 in juveniles of WAG/Rij rats could represent a compensatory mechanism determining the pattern of the age dependency in the disease manifestation by this model of absence epilepsy.