NeuroD2 controls inhibitory circuit formation in the molecular layer of the cerebellum

A single-cell multi-omic atlas spanning the adult rhesus macaque brain

Cataloging the diverse cellular architecture of the primate brain is crucial for understanding cognition, behavior, and disease in humans. Here, we generated a brain-wide single-cell multimodal molecular atlas of the rhesus macaque brain. Together, we profiled 2.58 M transcriptomes and 1.59 M epigenomes from single nuclei sampled from 30 regions across the adult brain. Cell composition differed extensively across the brain, revealing cellular signatures of region-specific functions. We also identified 1.19 M candidate regulatory elements, many previously unidentified, allowing us to explore the landscape of cis-regulatory grammar and neurological disease risk in a cell type-specific manner. Altogether, this multi-omic atlas provides an open resource for investigating the evolution of the human brain and identifying novel targets for disease interventions.

Protective Effect of Dexmedetomidine against Hyperoxia-Damaged Cerebellar Neurodevelopment in the Juvenile Rat

Impaired cerebellar development of premature infants and the associated impairment of cerebellar functions in cognitive development could be crucial factors for neurodevelopmental disorders. Anesthetic- and hyperoxia-induced neurotoxicity of the immature brain can lead to learning and behavioral disorders. Dexmedetomidine (DEX), which is associated with neuroprotective properties, is increasingly being studied for off-label use in the NICU. For this purpose, six-day-old Wistar rats (P6) were exposed to hyperoxia (80% O₂) or normoxia (21% O₂) for 24 h after DEX (5 µg/kg, i.p.) or vehicle (0.9% NaCl) application. An initial detection in the immature rat cerebellum was performed after the termination of hyperoxia at P7 and then after recovery in room air at P9, P11, and P14. Hyperoxia reduced the proportion of Calb1+-Purkinje cells and affected the dendrite length at P7 and/or P9/P11. Proliferating Pax6+-granule progenitors remained reduced after hyperoxia and until P14. The expression of neurotrophins and neuronal transcription factors/markers of proliferation, migration, and survival were also reduced by oxidative stress in different manners. DEX demonstrated protective effects on hyperoxia-injured Purkinje cells, and DEX without hyperoxia modulated neuronal transcription in the short term without any effects at the cellular level. DEX protects hyperoxia-damaged Purkinje cells and appears to differentially affect cerebellar granular cell neurogenesis following oxidative stress.

TrkB-dependent EphrinA reverse signaling regulates callosal axon fasciculate growth downstream of Neurod2/6

Abnormal development of corpus callosum is relatively common and causes a broad spectrum of cognitive impairments in humans. We use acallosal Neurod2/6-deficient mice to study callosal axon guidance within the ipsilateral cerebral cortex. Initial callosal tracts form but fail to traverse the ipsilateral cingulum and are not attracted towards the midline in the absence of Neurod2/6. We show that the restoration of Ephrin-A4 (EfnA4) expression in the embryonic neocortex of Neurod2/6-deficient embryos is sufficient to partially rescue targeted callosal axon growth towards the midline. EfnA4 cannot directly mediate reverse signaling within outgrowing axons, but it forms co-receptor complexes with TrkB (Ntrk2). The ability of EfnA4 to rescue the guided growth of a subset of callosal axons in Neurod2/6-deficient mice is abolished by the co-expression of dominant negative TrkBK571N (kinase-dead) or TrkBY515F (SHC-binding deficient) variants, but not by TrkBY816F (PLCγ1-binding deficient). Additionally, EphA4 is repulsive to EfnA4-positive medially projecting axons in organotypic brain slice culture. Collectively, we suggest that EfnA4-mediated reverse signaling acts via TrkB-SHC and is required for ipsilateral callosal axon growth accuracy towards the midline downstream of Neurod family factors.

Arc Regulates Transcription of Genes for Plasticity, Excitability and Alzheimer’s Disease

The immediate early gene Arc is a master regulator of synaptic function and a critical determinant of memory consolidation. Here, we show that Arc interacts with dynamic chromatin and closely associates with histone markers for active enhancers and transcription in cultured rat hippocampal neurons. Both these histone modifications, H3K27Ac and H3K9Ac, have recently been shown to be upregulated in late-onset Alzheimer’s disease (AD). When Arc induction by pharmacological network activation was prevented using a short hairpin RNA, the expression profile was altered for over 1900 genes, which included genes associated with synaptic function, neuronal plasticity, intrinsic excitability, and signalling pathways. Interestingly, about 100 Arc-dependent genes are associated with the pathophysiology of AD. When endogenous Arc expression was induced in HEK293T cells, the transcription of many neuronal genes was increased, suggesting that Arc can control expression in the absence of activated signalling pathways. Taken together, these data establish Arc as a master regulator of neuronal activity-dependent gene expression and suggest that it plays a significant role in the pathophysiology of AD.

Dendritic Inhibition by Shh Signaling-Dependent Stellate Cell Pool Is Critical for Motor Learning

Cerebellar inhibitory interneurons are important regulators of neural circuit activity for diverse motor and nonmotor functions. The molecular layer interneurons (MLIs), consisting of basket cells (BCs) and stellate cells (SCs), provide dendritic and somatic inhibitory synapses onto Purkinje cells, respectively. They are sequentially generated in an inside-out pattern from Pax2⁺ immature interneurons, which migrate from the prospective white matter to the ML of the cortex. However, little is known about how MLI subtype identities and pool sizes are determined, nor are their contributions to motor learning well understood. Here, we show that GABAergic progenitors fated to generate both BCs and SCs respond to the Sonic hedgehog (Shh) signal. Conditional abrogation of Shh signaling of either sex inhibited proliferation of GABAergic progenitors and reduced the number of Pax2⁺ cells, whereas persistent Shh pathway activation increased their numbers. These changes, however, did not affect early born BC numbers but selectively altered the SC pool size. Moreover, genetic depletion of GABAergic progenitors when BCs are actively generated also resulted in a specific reduction of SCs, suggesting that the specification of MLI subtypes is independent of Shh signaling and their birth order and likely occurs after Pax2⁺ cells settle into their laminar positions in an inside-out sequence. Mutant mice with reduced SC numbers displayed decreased dendritic inhibitory synapses and neurotransmission onto Purkinje cells, resulting in an impaired acquisition of eyeblink conditioning. These findings also reveal an essential role of Shh signaling-dependent SCs in regulating inhibitory dendritic synapses and motor learning.SIGNIFICANCE STATEMENT The cerebellar circuit that enables fine motor learning involves MLIs of BCs and SCs, which provide dendritic and somatic inhibitory synapses onto Purkinje cells. Little is known about how their identities and numbers are determined, nor are their specific contributions to motor learning well understood. We show that MLI subtypes are specified independent of Shh signaling and their birth orders but appear to occur in their terminal laminar positions according to the inside-out sequence. This finding challenges the current view that MLI subtypes are specified sequentially at the progenitor level. We also demonstrate that dendritic inhibition by Shh signaling-dependent SC pool is necessary for motor learning.

Morphological pseudotime ordering and fate mapping reveal diversification of cerebellar inhibitory interneurons

Understanding how diverse neurons are assembled into circuits requires a framework for describing cell types and their developmental trajectories. Here we combine genetic fate-mapping, pseudotemporal profiling of morphogenesis, and dual morphology and RNA labeling to resolve the diversification of mouse cerebellar inhibitory interneurons. Molecular layer interneurons (MLIs) derive from a common progenitor population but comprise diverse dendritic-, somatic-, and axon initial segment-targeting interneurons. Using quantitative morphology from 79 mature MLIs, we identify two discrete morphological types and presence of extensive within-class heterogeneity. Pseudotime trajectory inference using 732 developmental morphologies indicate the emergence of distinct MLI types during migration, before reaching their final positions. By comparing MLI identities from morphological and transcriptomic signatures, we demonstrate the dissociation between these modalities and that subtype divergence can be resolved from axonal morphogenesis prior to marker gene expression. Our study illustrates the utility of applying single-cell methods to quantify morphology for defining neuronal diversification.

The Conflicting Role of Caffeine Supplementation on Hyperoxia-Induced Injury on the Cerebellar Granular Cell Neurogenesis of Newborn Rats

Preterm birth disrupts cerebellar development, which may be mediated by systemic oxidative stress that damages neuronal developmental stages. Impaired cerebellar neurogenesis affects several downstream targets important for cognition, emotion, and speech. In this study, we demonstrate that oxidative stress induced with high oxygen (80%) for three or five postnatal days (P3/P5) could significantly damage neurogenesis and proliferative capacity of granular cell precursor and Purkinje cells in rat pups. Reversal of cellular neuronal damage after recovery to room air (P15) was augmented by treatment with caffeine. However, downstream transcripts important for migration and differentiation of postmitotic granular cells were irreversibly reduced by hyperoxia, without rescue by caffeine. Protective effects of caffeine in the cerebellum were limited to neuronal survival but failed to restore important transcript signatures.

Comparative description of the mRNA expression profile of Na+ /K+ -ATPase isoforms in adult mouse nervous system

Mutations in genes encoding Na⁺ /K⁺ -ATPase α1, α2, and α3 subunits cause a wide range of disabling neurological disorders, and dysfunction of Na⁺ /K⁺ -ATPase may contribute to neuronal injury in stroke and dementia. To better understand the pathogenesis of these diseases, it is important to determine the expression patterns of the different Na⁺ /K⁺ -ATPase subunits within the brain and among specific cell types. Using two available scRNA-Seq databases from the adult mouse nervous system, we examined the mRNA expression patterns of the different isoforms of the Na⁺ /K⁺ -ATPase α, β and Fxyd subunits at the single-cell level among brain regions and various neuronal populations. We subsequently identified specific types of neurons enriched with transcripts for α1 and α3 isoforms and elaborated how α3-expressing neuronal populations govern cerebellar neuronal circuits. We further analyzed the co-expression network for α1 and α3 isoforms, highlighting the genes that positively correlated with α1 and α3 expression. The top 10 genes for α1 were Chn2, Hpcal1, Nrgn, Neurod1, Selm, Kcnc1, Snrk, Snap25, Ckb and Ccndbp1 and for α3 were Sorcs3, Eml5, Neurod2, Ckb, Tbc1d4, Ptprz1, Pvrl1, Kirrel3, Pvalb, and Asic2.

The Role of Neurod Genes in Brain Development, Function, and Disease

Transcriptional regulation is essential for the correct functioning of cells during development and in postnatal life. The basic Helix-loop-Helix (bHLH) superfamily of transcription factors is well conserved throughout evolution and plays critical roles in tissue development and tissue maintenance. A subgroup of this family, called neural lineage bHLH factors, is critical in the development and function of the central nervous system. In this review, we will focus on the function of one subgroup of neural lineage bHLH factors, the Neurod family. The Neurod family has four members: Neurod1, Neurod2, Neurod4, and Neurod6. Available evidence shows that these four factors are key during the development of the cerebral cortex but also in other regions of the central nervous system, such as the cerebellum, the brainstem, and the spinal cord. We will also discuss recent reports that link the dysfunction of these transcription factors to neurological disorders in humans.

Excitatory granule neuron precursors orchestrate laminar localization and differentiation of cerebellar inhibitory interneuron subtypes

GABAergic interneurons migrate long distances through stereotyped migration programs toward specific laminar positions. During their migration, GABAergic interneurons are morphologically alike but then differentiate into a rich array of interneuron subtypes critical for brain function. How interneuron subtypes acquire their final phenotypic traits remains largely unknown. Here, we show that cerebellar molecular layer GABAergic interneurons, derived from the same progenitor pool, use separate migration paths to reach their laminar position and differentiate into distinct basket cell (BC) and stellate cell (SC) GABAergic interneuron subtypes. Using two-photon live imaging, we find that SC final laminar position requires an extra step of tangential migration supported by a subpopulation of glutamatergic granule cells (GCs). Conditional depletion of GCs affects SC differentiation but does not affect BCs. Our results reveal how timely feedforward control of inhibitory interneuron migration path regulates their terminal differentiation and, thus, establishment of the local inhibitory circuit assembly.

PKN1 Is a Novel Regulator of Hippocampal GluA1 Levels

Alterations in the processes that control α-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) expression, assembly and trafficking are closely linked to psychiatric and neurodegenerative disorders. We have recently shown that the serine/threonine kinase Protein kinase N1 (PKN1) is a developmentally active regulator of cerebellar synaptic maturation by inhibiting AKT and the neurogenic transcription factor neurogenic differentiation factor-2 (NeuroD2). NeuroD2 is involved in glutamatergic synaptic maturation by regulating expression levels of various synaptic proteins. Here we aimed to study the effect of Pkn1 knockout on AKT phosphorylation and NeuroD2 levels in the hippocampus and the subsequent expression levels of the NeuroD2 targets and AMPAR subunits: glutamate receptor 1 (GluA1) and GluA2/3. We show that PKN1 is expressed throughout the hippocampus. Interestingly, not only postnatal but also adult hippocampal phospho-AKT and NeuroD2 levels were significantly elevated upon Pkn1 knockout. Postnatal and adult Pkn1 ^-/- hippocampi showed enhanced expression of the AMPAR subunit GluA1, particularly in area CA1. Surprisingly, GluA2/3 levels were not different between both genotypes. In addition to higher protein levels, we also found an enhanced GluA1 content in the membrane fraction of postnatal and adult Pkn1 ^-/- animals, while GluA2/3 levels remained unchanged. This points toward a very specific regulation of GluA1 expression and/or trafficking by the novel PKN1-AKT-NeuroD2 axis. Considering the important role of GluA1 in hippocampal development as well as the pathophysiology of several disorders, ranging from Alzheimer's, to depression and schizophrenia, our results validate PKN1 for future studies into neurological disorders related to altered AMPAR subunit expression in the hippocampus.

Origins, Development, and Compartmentation of the Granule Cells of the Cerebellum

Granule cells (GCs) are the most numerous cell type in the cerebellum and indeed, in the brain: at least 99% of all cerebellar neurons are granule cells. In this review article, we first consider the formation of the upper rhombic lip, from which all granule cell precursors arise, and the way by which the upper rhombic lip generates the external granular layer, a secondary germinal epithelium that serves to amplify the upper rhombic lip precursors. Next, we review the mechanisms by which postmitotic granule cells are generated in the external granular layer and migrate radially to settle in the granular layer. In addition, we review the evidence that far from being a homogeneous population, granule cells come in multiple phenotypes with distinct topographical distributions and consider ways in which the heterogeneity of granule cells might arise during development.

Expansion of NEUROD2 phenotypes to include developmental delay without seizures

De novo heterozygous variants in the brain-specific transcription factor Neuronal Differentiation Factor 2 (NEUROD2) have been recently associated with early-onset epileptic encephalopathy and developmental delay. Here, we report an adolescent with developmental delay without seizures who was found to have a novel de novo heterozygous NEUROD2 missense variant, p.(Leu163Pro). Functional testing using an in vivo assay of neuronal differentiation in Xenopus laevis tadpoles demonstrated that the patient variant of NEUROD2 displays minimal protein activity, strongly suggesting a loss of function effect. In contrast, a second rare NEUROD2 variant, p.(Ala235Thr), identified in an adolescent with developmental delay but lacking parental studies for inheritance, showed normal in vivo NEUROD2 activity. We thus provide clinical, genetic, and functional evidence that NEUROD2 variants can lead to developmental delay without accompanying early-onset seizures, and demonstrate how functional testing can complement genetic data when determining variant pathogenicity.

Disruption of NEUROD2 causes a neurodevelopmental syndrome with autistic features via cell-autonomous defects in forebrain glutamatergic neurons

While the transcription factor NEUROD2 has recently been associated with epilepsy, its precise role during nervous system development remains unclear. Using a multi-scale approach, we set out to understand how Neurod2 deletion affects the development of the cerebral cortex in mice. In Neurod2 KO embryos, cortical projection neurons over-migrated, thereby altering the final size and position of layers. In juvenile and adults, spine density and turnover were dysregulated in apical but not basal compartments in layer 5 neurons. Patch-clamp recordings in layer 5 neurons of juvenile mice revealed increased intrinsic excitability. Bulk RNA sequencing showed dysregulated expression of many genes associated with neuronal excitability and synaptic function, whose human orthologs were strongly associated with autism spectrum disorders (ASD). At the behavior level, Neurod2 KO mice displayed social interaction deficits, stereotypies, hyperactivity, and occasionally spontaneous seizures. Mice heterozygous for Neurod2 had similar defects, indicating that Neurod2 is haploinsufficient. Finally, specific deletion of Neurod2 in forebrain excitatory neurons recapitulated cellular and behavioral phenotypes found in constitutive KO mice, revealing the region-specific contribution of dysfunctional Neurod2 in symptoms. Informed by these neurobehavioral features in mouse mutants, we identified eleven patients from eight families with a neurodevelopmental disorder including intellectual disability and ASD associated with NEUROD2 pathogenic mutations. Our findings demonstrate crucial roles for Neurod2 in neocortical development, whose alterations can cause neurodevelopmental disorders including intellectual disability and ASD.

Identifying transcriptomic correlates of histology using deep learning

Linking phenotypes to specific gene expression profiles is an extremely important problem in biology, which has been approached mainly by correlation methods or, more fundamentally, by studying the effects of gene perturbations. However, genome-wide perturbations involve extensive experimental efforts, which may be prohibitive for certain organisms. On the other hand, the characterization of the various phenotypes frequently requires an expert's subjective interpretation, such as a histopathologist's description of tissue slide images in terms of complex visual features (e.g. 'acinar structures'). In this paper, we use Deep Learning to eliminate the inherent subjective nature of these visual histological features and link them to genomic data, thus establishing a more precisely quantifiable correlation between transcriptomes and phenotypes. Using a dataset of whole slide images with matching gene expression data from 39 normal tissue types, we first developed a Deep Learning tissue classifier with an accuracy of 94%. Then we searched for genes whose expression correlates with features inferred by the classifier and demonstrate that Deep Learning can automatically derive visual (phenotypical) features that are well correlated with the transcriptome and therefore biologically interpretable. As we are particularly concerned with interpretability and explainability of the inferred histological models, we also develop visualizations of the inferred features and compare them with gene expression patterns determined by immunohistochemistry. This can be viewed as a first step toward bridging the gap between the level of genes and the cellular organization of tissues.

Transcriptomic and epigenomic dynamics associated with development of human iPSC-derived GABAergic interneurons

GABAergic interneurons (GINs) are a heterogeneous population of inhibitory neurons that collectively contribute to the maintenance of normal neuronal excitability and network activity. Identification of the genetic regulatory elements and transcription factors that contribute toward GIN function may provide new insight into the pathways underlying proper GIN activity while also indicating potential therapeutic targets for GIN-associated disorders, such as schizophrenia and epilepsy. In this study, we examined the temporal changes in gene expression and chromatin accessibility during GIN development by performing transcriptomic and epigenomic analyses on human induced pluripotent stem cell-derived neurons at 22, 50 and 78 days (D) post-differentiation. We observed 13 221 differentially accessible regions (DARs) of chromatin that associate with temporal changes in gene expression at D78 and D50, relative to D22. We also classified families of transcription factors that are increasingly enriched at DARs during differentiation, indicating regulatory networks that likely drive GIN development. Collectively, these data provide a resource for examining the molecular networks regulating GIN functionality.