ScISOr-ATAC reveals convergent and divergent splicing and chromatin specificities between matched cell types across cortical regions, evolution, and in Alzheimer's Disease.

Multimodal measurements have become widespread in genomics, however measuring open chromatin accessibility and splicing simultaneously in frozen brain tissues remains unconquered. Hence, we devised Single-Cell-ISOform-RNA sequencing coupled with the Assay-for-Transposase-Accessible-Chromatin (ScISOr-ATAC). We utilized ScISOr-ATAC to assess whether chromatin and splicing alterations in the brain convergently affect the same cell types or divergently different ones. We applied ScISOr-ATAC to three major conditions: comparing (i) the Rhesus macaque (Macaca mulatta) prefrontal cortex (PFC) and visual cortex (VIS), (ii) cross species divergence of Rhesus macaque versus human PFC, as well as (iii) dysregulation in Alzheimer's disease in human PFC. We found that among cortical-layer biased excitatory neuron subtypes, splicing is highly brain-region specific for L3-5/L6 IT_RORB neurons, moderately specific in L2-3 IT_CUX2.RORB neurons and unspecific in L2-3 IT_CUX2 neurons. In contrast, at the chromatin level, L2-3 IT_CUX2.RORB neurons show the highest brain-region specificity compared to other subtypes. Likewise, when comparing human and macaque PFC, strong evolutionary divergence on one molecular modality does not necessarily imply strong such divergence on another molecular level in the same cell type. Finally, in Alzheimer's disease, oligodendrocytes show convergently high dysregulation in both chromatin and splicing. However, chromatin and splicing dysregulation most strongly affect distinct oligodendrocyte subtypes. Overall, these results indicate that chromatin and splicing can show convergent or divergent results depending on the performed comparison, justifying the need for their concurrent measurement to investigate complex systems. Taken together, ScISOr-ATAC allows for the characterization of single-cell splicing and chromatin patterns and the comparison of sample groups in frozen brain samples.

Neural and metabolic dysregulation in PMM2-deficient human in vitro neural models.

Phosphomannomutase 2-congenital disorder of glycosylation (PMM2-CDG) is a rare inborn error of metabolism caused by deficiency of the PMM2 enzyme, which leads to impaired protein glycosylation. While the disorder presents with primarily neurological symptoms, there is limited knowledge about the specific brain-related changes caused by PMM2 deficiency. Here, we demonstrate aberrant neural activity in 2D neuronal networks from PMM2-CDG individuals. Utilizing multi-omics datasets from 3D human cortical organoids (hCOs) derived from PMM2-CDG individuals, we identify widespread decreases in protein glycosylation, highlighting impaired glycosylation as a key pathological feature of PMM2-CDG, as well as impaired mitochondrial structure and abnormal glucose metabolism in PMM2-deficient hCOs, indicating disturbances in energy metabolism. Correlation between PMM2 enzymatic activity in hCOs and symptom severity suggests that the level of PMM2 enzyme function directly influences neurological manifestations. These findings enhance our understanding of specific brain-related perturbations associated with PMM2-CDG, offering insights into the underlying mechanisms and potential directions for therapeutic interventions.

Gliomas: a reflection of temporal gliogenic principles.

The hijacking of early developmental programs is a canonical feature of gliomas where neoplastic cells resemble neurodevelopmental lineages and possess mechanisms of stem cell resilience. Given these parallels, uncovering how and when in developmental time gliomagenesis intersects with normal trajectories can greatly inform our understanding of tumor biology. Here, we review how elapsing time impacts the developmental principles of astrocyte (AS) and oligodendrocyte (OL) lineages, and how these same temporal programs are replicated, distorted, or circumvented in pathological settings such as gliomas. Additionally, we discuss how normal gliogenic processes can inform our understanding of the temporal progression of gliomagenesis, including when in developmental time gliomas originate, thrive, and can be pushed towards upon therapeutic coercion.

Cross-species analysis identifies mitochondrial dysregulation as a functional consequence of the schizophrenia-associated 3q29 deletion.

The 1.6-megabase deletion at chromosome 3q29 (3q29Del) is the strongest identified genetic risk factor for schizophrenia, but the effects of this variant on neurodevelopment are not well understood. We interrogated the developing neural transcriptome in two experimental model systems with complementary advantages: isogenic human cortical organoids and isocortex from the 3q29Del mouse model. We profiled transcriptomes from isogenic cortical organoids that were aged for 2 and 12 months, as well as perinatal mouse isocortex, all at single-cell resolution. Systematic pathway analysis implicated dysregulation of mitochondrial function and energy metabolism. These molecular signatures were supported by analysis of oxidative phosphorylation protein complex expression in mouse brain and assays of mitochondrial function in engineered cell lines, which revealed a lack of metabolic flexibility and a contribution of the 3q29 gene PAK2. Together, these data indicate that metabolic disruption is associated with 3q29Del and is conserved across species.

Identification of ligand-receptor pairs that drive human astrocyte development.

Extrinsic signaling between diverse cell types is crucial for nervous system development. Ligand binding is a key driver of developmental processes. Nevertheless, it remains a significant challenge to disentangle which and how extrinsic signals act cooperatively to affect changes in recipient cells. In the developing human brain, cortical progenitors transition from neurogenesis to gliogenesis in a stereotyped sequence that is in part influenced by extrinsic ligands. Here we used published transcriptomic data to identify and functionally test five ligand-receptor pairs that synergistically drive human astrogenesis. We validate the synergistic contributions of TGFß2, NLGN1, TSLP, DKK1 and BMP4 ligands on astrocyte development in both hCOs and primary fetal tissue. We confirm that the cooperative capabilities of these five ligands are greater than their individual capacities. Additionally, we discovered that their combinatorial effects converge in part on the mTORC1 signaling pathway, resulting in transcriptomic and morphological features of astrocyte development. Our data-driven framework can leverage single-cell and bulk genomic data to generate and test functional hypotheses surrounding cell-cell communication regulating neurodevelopmental processes.

GPR37L1 controls maturation and organization of cortical astrocytes during development.

Astrocyte maturation is crucial to proper brain development and function. This maturation process includes the ramification of astrocytic morphology and the establishment of astrocytic domains. While this process has been well-studied, the mechanisms by which astrocyte maturation is initiated are not well understood. GPR37L1 is an astrocyte-specific G protein-coupled receptor (GPCR) that is predominantly expressed in mature astrocytes and has been linked to the modulation of seizure susceptibility in both humans and mice. To investigate the role of GPR37L1 in astrocyte biology, RNA-seq analyses were performed on astrocytes immunopanned from P7 Gpr37L1-/- knockout (L1KO) mouse cortex and compared to those from wild-type (WT) mouse cortex. These RNA-seq studies revealed that pathways involved in central nervous system development were altered and that L1KO cortical astrocytes express lower amounts of mature astrocytic genes compared to WT astrocytes. Immunohistochemical studies of astrocytes from L1KO mouse brain revealed that these astrocytes exhibit overall shorter total process length, and are also less complex and spaced further apart from each other in the mouse cortex. This work sheds light on how GPR37L1 regulates cellular processes involved in the control of astrocyte biology and maturation.

Tyrosinase-induced neuromelanin accumulation triggers rapid dysregulation and degeneration of the mouse locus coeruleus.

Motor symptoms in Parkinson's disease (PD) are caused by degeneration of dopamine (DA) neurons of the substantia nigra (SN), while early non-motor symptoms such as anxiety and sleep disturbances are likely mediated by dysfunction of locus coeruleus (LC) norepinephrine (NE) neurons. The LC develops α-synuclein pathology prior to SN DA neurons in PD, and later undergoes degeneration, but the mechanisms responsible for its vulnerability are unknown. The SN and LC are the only structures in the brain that produces appreciable amounts of neuromelanin (NM), a dark cytoplasmic pigment. It has been proposed that NM initially plays a protective role by sequestering toxic catecholamine metabolites and heavy metals, but may become harmful during aging and PD as they overwhelm cellular machinery and are released during neurodegeneration. Rodents do not naturally produce NM, limiting the study of causal relationships between NM and PD-associated LC pathology. Adapting a viral-mediated approach for expression of human tyrosinase, the enzyme responsible for peripheral melanin production, we successfully promoted pigmentation in mouse LC neurons that recapitulates key features of endogenous NM found in primates, including eumelanin and pheomelanin, lipid droplets, and a double-membrane encasement. Pigment expression results in mild neurodegeneration, reduced NE levels, transcriptional changes, and novelty-induced anxiety phenotypes as early as 1-week post-injection. By 6-weeks, NM accumulation is associated with severe LC neurodegeneration and a robust neuroinflammatory response. These phenotypes are reminiscent of LC dysfunction in PD, validating this model for studying the consequences of pigment accumulation in the LC as it relates to neurodegenerative disease.

Cross-species transcriptomic analysis identifies mitochondrial dysregulation as a functional consequence of the schizophrenia-associated 3q29 deletion.

Recent advances in the genetics of schizophrenia (SCZ) have identified rare variants that confer high disease risk, including a 1.6 Mb deletion at chromosome 3q29 with a staggeringly large effect size (O.R. > 40). Understanding the impact of the 3q29 deletion (3q29Del) on the developing CNS may therefore lead to insights about the pathobiology of schizophrenia. To gain clues about the molecular and cellular perturbations caused by the 3q29 deletion, we interrogated transcriptomic effects in two experimental model systems with complementary advantages: isogenic human forebrain cortical organoids and isocortex from the 3q29Del mouse model. We first created isogenic lines by engineering the full 3q29Del into an induced pluripotent stem cell line from a neurotypical individual. We profiled transcriptomes from isogenic cortical organoids that were aged for 2 months and 12 months, as well as day p7 perinatal mouse isocortex, all at single cell resolution. Differential expression analysis by genotype in each cell-type cluster revealed that more than half of the differentially expressed genes identified in mouse cortex were also differentially expressed in human cortical organoids, and strong correlations were observed in mouse-human differential gene expression across most major cell-types. We systematically filtered differentially expressed genes to identify changes occurring in both model systems. Pathway analysis on this filtered gene set implicated dysregulation of mitochondrial function and energy metabolism, although the direction of the effect was dependent on developmental timepoint. Transcriptomic changes were validated at the protein level by analysis of oxidative phosphorylation protein complexes in mouse brain tissue. Assays of mitochondrial function in human heterologous cells further confirmed robust mitochondrial dysregulation in 3q29Del cells, and these effects are partially recapitulated by ablation of the 3q29Del gene PAK2 . Taken together these data indicate that metabolic disruption is associated with 3q29Del and is conserved across species. These results converge with data from other rare SCZ-associated variants as well as idiopathic schizophrenia, suggesting that mitochondrial dysfunction may be a significant but overlooked contributing factor to the development of psychotic disorders. This cross-species scRNA-seq analysis of the SCZ-associated 3q29 deletion reveals that this copy number variant may produce early and persistent changes in cellular metabolism that are relevant to human neurodevelopment.

The Neurotoxin DSP-4 Dysregulates the Locus Coeruleus-Norepinephrine System and Recapitulates Molecular and Behavioral Aspects of Prodromal Neurodegenerative Disease.

The noradrenergic locus coeruleus (LC) is among the earliest sites of tau and α-synuclein pathology in Alzheimer's disease (AD) and Parkinson's disease (PD), respectively. The onset of these pathologies coincides with loss of noradrenergic fibers in LC target regions and the emergence of prodromal symptoms including sleep disturbances and anxiety. Paradoxically, these prodromal symptoms are indicative of a noradrenergic hyperactivity phenotype, rather than the predicted loss of norepinephrine (NE) transmission following LC damage, suggesting the engagement of complex compensatory mechanisms. Because current therapeutic efforts are targeting early disease, interest in the LC has grown, and it is critical to identify the links between pathology and dysfunction. We employed the LC-specific neurotoxin N-(2-chloroethyl)-N-ethyl-2-bromobenzylamine (DSP-4), which preferentially damages LC axons, to model early changes in the LC-NE system pertinent to AD and PD in male and female mice. DSP-4 (two doses of 50 mg/kg, one week apart) induced LC axon degeneration, triggered neuroinflammation and oxidative stress, and reduced tissue NE levels. There was no LC cell death or changes to LC firing, but transcriptomics revealed reduced expression of genes that define noradrenergic identity and other changes relevant to neurodegenerative disease. Despite the dramatic loss of LC fibers, NE turnover and signaling were elevated in terminal regions and were associated with anxiogenic phenotypes in multiple behavioral tests. These results represent a comprehensive analysis of how the LC-NE system responds to axon/terminal damage reminiscent of early AD and PD at the molecular, cellular, systems, and behavioral levels, and provides potential mechanisms underlying prodromal neuropsychiatric symptoms.

An RNA-sequencing transcriptome of the rodent Schwann cell response to peripheral nerve injury.

BACKGROUND: The important contribution of glia to mechanisms of injury and repair of the nervous system is increasingly recognized. In stark contrast to the central nervous system (CNS), the peripheral nervous system (PNS) has a remarkable capacity for regeneration after injury. Schwann cells are recognized as key contributors to PNS regeneration, but the molecular underpinnings of the Schwann cell response to injury and how they interact with the inflammatory response remain incompletely understood. METHODS: We completed bulk RNA-sequencing of Schwann cells purified acutely using immunopanning from the naïve and injured rodent sciatic nerve at 3, 5, and 7 days post-injury. We used qRT-PCR and in situ hybridization to assess cell purity and probe dataset integrity. Finally, we used bioinformatic analysis to probe Schwann cell-specific injury-induced modulation of cellular pathways. RESULTS: Our data confirm Schwann cell purity and validate RNAseq dataset integrity. Bioinformatic analysis identifies discrete modules of genes that follow distinct patterns of regulation in the 1st days after injury and their corresponding molecular pathways. These findings enable improved differentiation of myeloid and glial components of neuroinflammation after peripheral nerve injury and highlight novel molecular aspects of the Schwann cell injury response such as acute downregulation of the AGE/RAGE pathway and of secreted molecules Sparcl1 and Sema5a. CONCLUSIONS: We provide a helpful resource for further deciphering the Schwann cell injury response and a depth of transcriptional data that can complement the findings of recent single cell sequencing approaches. As more data become available on the response of CNS glia to injury, we anticipate that this dataset will provide a valuable platform for understanding key differences in the PNS and CNS glial responses to injury and for designing approaches to ameliorate CNS regeneration.

Single-nuclei isoform RNA sequencing unlocks barcoded exon connectivity in frozen brain tissue.

Single-nuclei RNA sequencing characterizes cell types at the gene level. However, compared to single-cell approaches, many single-nuclei cDNAs are purely intronic, lack barcodes and hinder the study of isoforms. Here we present single-nuclei isoform RNA sequencing (SnISOr-Seq). Using microfluidics, PCR-based artifact removal, target enrichment and long-read sequencing, SnISOr-Seq increased barcoded, exon-spanning long reads 7.5-fold compared to naive long-read single-nuclei sequencing. We applied SnISOr-Seq to adult human frontal cortex and found that exons associated with autism exhibit coordinated and highly cell-type-specific inclusion. We found two distinct combination patterns: those distinguishing neural cell types, enriched in TSS-exon, exon-polyadenylation-site and non-adjacent exon pairs, and those with multiple configurations within one cell type, enriched in adjacent exon pairs. Finally, we observed that human-specific exons are almost as tightly coordinated as conserved exons, implying that coordination can be rapidly established during evolution. SnISOr-Seq enables cell-type-specific long-read isoform analysis in human brain and in any frozen or hard-to-dissociate sample.

Growing Glia: Cultivating Human Stem Cell Models of Gliogenesis in Health and Disease.

Glia are present in all organisms with a central nervous system but considerably differ in their diversity, functions, and numbers. Coordinated efforts across many model systems have contributed to our understanding of glial-glial and neuron-glial interactions during nervous system development and disease, but human glia exhibit prominent species-specific attributes. Limited access to primary samples at critical developmental timepoints constrains our ability to assess glial contributions in human tissues. This challenge has been addressed throughout the past decade via advancements in human stem cell differentiation protocols that now offer the ability to model human astrocytes, oligodendrocytes, and microglia. Here, we review the use of novel 2D cell culture protocols, 3D organoid models, and bioengineered systems derived from human stem cells to study human glial development and the role of glia in neurodevelopmental disorders.

Elevated preoptic brain activity in zebrafish glial glycine transporter mutants is linked to lethargy-like behaviors and delayed emergence from anesthesia.

Delayed emergence from anesthesia was previously reported in a case study of a child with Glycine Encephalopathy. To investigate the neural basis of this delayed emergence, we developed a zebrafish glial glycine transporter (glyt1 - / -) mutant model. We compared locomotor behaviors; dose-response curves for tricaine, ketamine, and 2,6-diisopropylphenol (propofol); time to emergence from these anesthetics; and time to emergence from propofol after craniotomy in glyt1-/- mutants and their siblings. To identify differentially active brain regions in glyt1-/- mutants, we used pERK immunohistochemistry as a proxy for brain-wide neuronal activity. We show that glyt1-/- mutants initiated normal bouts of movement less frequently indicating lethargy-like behaviors. Despite similar anesthesia dose-response curves, glyt1-/- mutants took over twice as long as their siblings to emerge from ketamine or propofol, mimicking findings from the human case study. Reducing glycine levels rescued timely emergence in glyt1-/- mutants, pointing to a causal role for elevated glycine. Brain-wide pERK staining showed elevated activity in hypnotic brain regions in glyt1-/- mutants under baseline conditions and a delay in sensorimotor integration during emergence from anesthesia. Our study links elevated activity in preoptic brain regions and reduced sensorimotor integration to lethargy-like behaviors and delayed emergence from propofol in glyt1-/- mutants.

A spatially resolved brain region- and cell type-specific isoform atlas of the postnatal mouse brain.

Splicing varies across brain regions, but the single-cell resolution of regional variation is unclear. We present a single-cell investigation of differential isoform expression (DIE) between brain regions using single-cell long-read sequencing in mouse hippocampus and prefrontal cortex in 45 cell types at postnatal day 7 ( www.isoformAtlas.com ). Isoform tests for DIE show better performance than exon tests. We detect hundreds of DIE events traceable to cell types, often corresponding to functionally distinct protein isoforms. Mostly, one cell type is responsible for brain-region specific DIE. However, for fewer genes, multiple cell types influence DIE. Thus, regional identity can, although rarely, override cell-type specificity. Cell types indigenous to one anatomic structure display distinctive DIE, e.g. the choroid plexus epithelium manifests distinct transcription-start-site usage. Spatial transcriptomics and long-read sequencing yield a spatially resolved splicing map. Our methods quantify isoform expression with cell-type and spatial resolution and it contributes to further our understanding of how the brain integrates molecular and cellular complexity.

3D Bioprinting of Neural Tissues.

The human nervous system is a remarkably complex physiological network that is inherently challenging to study because of obstacles to acquiring primary samples. Animal models offer powerful alternatives to study nervous system development, diseases, and regenerative processes, however, they are unable to address some species-specific features of the human nervous system. In vitro models of the human nervous system have expanded in prevalence and sophistication, but still require further advances to better recapitulate microenvironmental and cellular features. The field of neural tissue engineering (TE) is rapidly adopting new technologies that enable scientists to precisely control in vitro culture conditions and to better model nervous system formation, function, and repair. 3D bioprinting is one of the major TE technologies that utilizes biocompatible hydrogels to create precisely patterned scaffolds, designed to enhance cellular responses. This review focuses on the applications of 3D bioprinting in the field of neural TE. Important design parameters are considered when bioprinting neural stem cells are discussed. The emergence of various bioprinted in vitro platforms are also reviewed for developmental and disease modeling and drug screening applications within the central and peripheral nervous systems, as well as their use as implants for in vivo regenerative therapies.

Astrocyte-to-astrocyte contact and a positive feedback loop of growth factor signaling regulate astrocyte maturation.

Astrocytes are critical for the development and function of the central nervous system. In developing brains, immature astrocytes undergo morphological, molecular, cellular, and functional changes as they mature. Although the mechanisms that regulate the maturation of other major cell types in the central nervous system such as neurons and oligodendrocytes have been extensively studied, little is known about the cellular and molecular mechanisms that control astrocyte maturation. Here, we identified molecular markers of astrocyte maturation and established an in vitro assay for studying the mechanisms of astrocyte maturation. Maturing astrocytes in vitro exhibit similar molecular changes and represent multiple molecular subtypes of astrocytes found in vivo. Using this system, we found that astrocyte-to-astrocyte contact strongly promotes astrocyte maturation. In addition, secreted signals from microglia, oligodendrocyte precursor cells, or endothelial cells affect a small subset of astrocyte genes but do not consistently change astrocyte maturation. To identify molecular mechanisms underlying astrocyte maturation, we treated maturing astrocytes with molecules that affect the function of tumor-associated genes. We found that a positive feedback loop of heparin-binding epidermal growth factor-like growth factor (HBEGF) and epidermal growth factor receptor (EGFR) signaling regulates astrocytes maturation. Furthermore, HBEGF, EGFR, and tumor protein 53 (TP53) affect the expression of genes important for cilium development, the circadian clock, and synapse function. These results revealed cellular and molecular mechanisms underlying astrocytes maturation with implications for the understanding of glioblastoma.

Differentiation and maturation of oligodendrocytes in human three-dimensional neural cultures.

Investigating human oligodendrogenesis and the interaction of oligodendrocytes with neurons and astrocytes would accelerate our understanding of the mechanisms underlying white matter disorders. However, this is challenging because of the limited accessibility of functional human brain tissue. Here, we developed a new differentiation method of human induced pluripotent stem cells to generate three-dimensional brain organoids that contain oligodendrocytes as well as neurons and astrocytes, called human oligodendrocyte spheroids. We found that oligodendrocyte lineage cells derived in human oligodendrocyte spheroids transitioned through developmental stages similar to primary human oligodendrocytes and that the migration of oligodendrocyte lineage cells and their susceptibility to lysolecithin exposure could be captured by live imaging. Moreover, their morphology changed as they matured over time in vitro and started myelinating neurons. We anticipate that this method can be used to study oligodendrocyte development, myelination, and interactions with other major cell types in the CNS.

Single-cell isoform RNA sequencing characterizes isoforms in thousands of cerebellar cells.

Full-length RNA sequencing (RNA-Seq) has been applied to bulk tissue, cell lines and sorted cells to characterize transcriptomes, but applying this technology to single cells has proven to be difficult, with less than ten single-cell transcriptomes having been analyzed thus far. Although single splicing events have been described for ≤200 single cells with statistical confidence, full-length mRNA analyses for hundreds of cells have not been reported. Single-cell short-read 3' sequencing enables the identification of cellular subtypes, but full-length mRNA isoforms for these cell types cannot be profiled. We developed a method that starts with bulk tissue and identifies single-cell types and their full-length RNA isoforms without fluorescence-activated cell sorting. Using single-cell isoform RNA-Seq (ScISOr-Seq), we identified RNA isoforms in neurons, astrocytes, microglia, and cell subtypes such as Purkinje and Granule cells, and cell-type-specific combination patterns of distant splice sites. We used ScISOr-Seq to improve genome annotation in mouse Gencode version 10 by determining the cell-type-specific expression of 18,173 known and 16,872 novel isoforms.

Generation and assembly of human brain region-specific three-dimensional cultures.

The ability to generate region-specific three-dimensional (3D) models to study human brain development offers great promise for understanding the nervous system in both healthy individuals and patients. In this protocol, we describe how to generate and assemble subdomain-specific forebrain spheroids, also known as brain region-specific organoids, from human pluripotent stem cells (hPSCs). We describe how to pattern the neural spheroids toward either a dorsal forebrain or a ventral forebrain fate, establishing human cortical spheroids (hCSs) and human subpallial spheroids (hSSs), respectively. We also describe how to combine the neural spheroids in vitro to assemble forebrain assembloids that recapitulate the interactions of glutamatergic and GABAergic neurons seen in vivo. Astrocytes are also present in the human forebrain-specific spheroids, and these undergo maturation when the forebrain spheroids are cultured long term. The initial generation of neural spheroids from hPSCs occurs in <1 week, with regional patterning occurring over the subsequent 5 weeks. After the maturation stage, brain region-specific spheroids are amenable to a variety of assays, including live-cell imaging, calcium dynamics, electrophysiology, cell purification, single-cell transcriptomics, and immunohistochemistry studies. Once generated, forebrain spheroids can also be matured for >24 months in culture.