A Single-Cell Transcriptomic Atlas of Human Neocortical Development during Mid-gestation

A multi-omic single-cell landscape of cellular diversification in the developing human cerebral cortex

The vast neuronal diversity in the human neocortex is vital for high-order brain functions, necessitating elucidation of the regulatory mechanisms underlying such unparalleled diversity. However, recent studies have yet to comprehensively reveal the diversity of neurons and the molecular logic of neocortical origin in humans at single-cell resolution through profiling transcriptomic or epigenomic landscapes, owing to the application of unimodal data alone to depict exceedingly heterogeneous populations of neurons. In this study, we generated a comprehensive compendium of the developing human neocortex by simultaneously profiling gene expression and open chromatin from the same cell. We computationally reconstructed the differentiation trajectories of excitatory projection neurons of cortical origin and inferred the regulatory logic governing lineage bifurcation decisions for neuronal diversification. We demonstrated that neuronal diversity arises from progenitor cell lineage specificity and postmitotic differentiation at distinct stages. Our data paves the way for understanding the primarily coordinated regulatory logic for neuronal diversification in the neocortex.

Molecular Lineages and Spatial Distributions of Subplate Neurons in the Human Fetal Cerebral Cortex

The expansion of neural progenitors and production of distinct neurons are crucial for architectural assembly and formation of connectivity in human brains. Subplate neurons (SPNs) are among the firstborn neurons in the human fetal cerebral cortex, and play a critical role in establishing intra- and extracortical connections. However, little is known about SPN origin and developmental lineages. In this study, spatial landscapes and molecular trajectories of SPNs in the human fetal cortices from gestational weeks (GW) 10 to 25 are created by performing spatial transcriptomics and single-cell RNA sequencing. Genes known to be evolutionarily human-specific and genes associated with extracellular matrices (ECMs) are found to maintain stable proportions of subplate neurons among other neuronal types. Enriched ECM gene expression in SPNs varies in distinct cortical regions, with the highest level in the frontal lobe of human fetal brains. This study reveals molecular origin and lineage specification of subplate neurons in the human fetal cerebral cortices, and highlights underpinnings of SPNs to cortical neurogenesis and early structural folding.

Transcriptional and cellular response of hiPSC-derived microglia-neural progenitor co-cultures exposed to IL-6

Elevated interleukin (IL-)6 levels during prenatal development have been linked to increased risk for neurodevelopmental disorders (NDD) in the offspring, but the mechanism remains unclear. Human-induced pluripotent stem cell (hiPSC) models offer a valuable tool to study the effects of IL-6 on features relevant for human neurodevelopment in vitro. We previously reported that hiPSC-derived microglia-like cells (MGLs) respond to IL-6, but neural progenitor cells (NPCs) in monoculture do not. Therefore, we investigated whether co-culturing hiPSC-derived MGLs with NPCs would trigger a cellular response to IL-6 stimulation via secreted factors from the MGLs. Using N=4 donor lines without psychiatric diagnosis, we first confirmed that NPCs can respond to IL-6 through trans-signalling when recombinant IL-6Ra is present, and that this response is dose-dependent. MGLs secreted soluble IL-6R, but at lower levels than found in vivo and below that needed to activate trans-signalling in NPCs. Whilst transcriptomic and secretome analysis confirmed that MGLs undergo substantial transcriptomic changes after IL-6 exposure and subsequently secrete a cytokine milieu, NPCs in co-culture with MGLs exhibited a minimal transcriptional response. Furthermore, there were no significant cell fate-acquisition changes when differentiated into post-mitotic cultures, nor alterations in synaptic densities in mature neurons. These findings highlight the need to investigate if trans-IL-6 signalling to NPCs is a relevant disease mechanism linking prenatal IL-6 exposure to increased risk for psychiatric disorders. Moreover, our findings underscore the importance of establishing more complex in vitro human models with diverse cell types, which may show cell-specific responses to microglia-released cytokines to fully understand how IL-6 exposure may influence human neurodevelopment.

A polarized FGF8 source specifies frontotemporal signatures in spatially oriented cell populations of cortical assembloids

Organoids generating major cortical cell types in distinct compartments are used to study cortical development, evolution and disorders. However, the lack of morphogen gradients imparting cortical positional information and topography in current systems hinders the investigation of complex phenotypes. Here, we engineer human cortical assembloids by fusing an organizer-like structure expressing fibroblast growth factor 8 (FGF8) with an elongated organoid to enable the controlled modulation of FGF8 signaling along the longitudinal organoid axis. These polarized cortical assembloids mount a position-dependent transcriptional program that in part matches the in vivo rostrocaudal gene expression patterns and that is lost upon mutation in the FGFR3 gene associated with temporal lobe malformations and intellectual disability. By producing spatially oriented cell populations with signatures related to frontal and temporal area identity within individual assembloids, this model recapitulates in part the early transcriptional divergence embedded in the protomap and enables the study of cortical area-relevant alterations underlying human disorders.

Comparative single-cell multiome identifies evolutionary changes in neural progenitor cells during primate brain development

Understanding the cellular and genetic mechanisms driving human-specific features of cortical development remains a challenge. We generated a cell-type resolved atlas of transcriptome and chromatin accessibility in the developing macaque and mouse prefrontal cortex (PFC). Comparing with published human data, our findings demonstrate that although the cortex cellular composition is overall conserved across species, progenitor cells show significant evolutionary divergence in cellular properties. Specifically, human neural progenitors exhibit extensive transcriptional rewiring in growth factor and extracellular matrix (ECM) pathways. Expression of the human-specific progenitor marker ITGA2 in the fetal mouse cortex increases the progenitor proliferation and the proportion of upper-layer neurons. These transcriptional divergences are primarily driven by altered activity in the distal regulatory elements. The chromatin regions with human-gained accessibility are enriched with human-specific sequence changes and polymorphisms linked to intelligence and neuropsychiatric disorders. Our results identify evolutionary changes in neural progenitors and putative gene regulatory mechanisms shaping primate brain evolution.

Neuronal SLC39A8 deficiency impairs cerebellar development by altering manganese homeostasis

Solute carrier family 39, member 8 (SLC39A8), is a transmembrane transporter that mediates the cellular uptake of zinc, iron, and manganese (Mn). Human genetic studies document the involvement of SLC39A8 in Mn homeostasis, brain development, and function. However, the role and pathophysiological mechanisms of SLC39A8 in the central nervous system remain elusive. We generated Slc39a8 neuron-specific knockout (Slc39a8-NSKO) mice to study SLC39A8 function in neurons. The Slc39a8-NSKO mice displayed markedly decreased Mn levels in the whole brain and brain regions, especially the cerebellum. Radiotracer studies using 54Mn revealed that Slc39a8-NSKO mice had impaired brain uptake of Mn. Slc39a8-NSKO cerebellums exhibited morphological defects and abnormal dendritic arborization of Purkinje cells. Reduced neurogenesis and increased apoptotic cell death occurred in the cerebellar external granular layer of Slc39a8-NSKO mice. Brain Mn deficiency in Slc39a8-NSKO mice was associated with motor dysfunction. Unbiased RNA-Seq analysis revealed downregulation of key pathways relevant to neurodevelopment and synaptic plasticity, including cAMP signaling pathway genes. We further demonstrated that Slc39a8 was required for the optimal transcriptional response to the cAMP-mediated signaling pathway. In summary, our study highlighted the essential roles of SLC39A8 in brain Mn uptake and cerebellum development and functions.

Integrated electrophysiological and genomic profiles of single cells reveal spiking tumor cells in human glioma

Prior studies have described the complex interplay that exists between glioma cells and neurons; however, the electrophysiological properties endogenous to glioma cells remain obscure. To address this, we employed Patch-sequencing (Patch-seq) on human glioma specimens and found that one-third of patched cells in IDH mutant (IDHmut) tumors demonstrate properties of both neurons and glia. To define these hybrid cells (HCs), which fire single, short action potentials, and discern if they are of tumoral origin, we developed the single cell rule association mining (SCRAM) computational tool to annotate each cell individually. SCRAM revealed that HCs possess select features of GABAergic neurons and oligodendrocyte precursor cells, and include both tumor and non-tumor cells. These studies characterize the combined electrophysiological and molecular properties of human glioma cells and describe a cell type in human glioma with unique electrophysiological and transcriptomic properties that may also exist in the non-tumor brain.

Temporally distinct 3D multi-omic dynamics in the developing human brain

The human hippocampus and prefrontal cortex play critical roles in learning and cognition1,2, yet the dynamic molecular characteristics of their development remain enigmatic. Here we investigated the epigenomic and three-dimensional chromatin conformational reorganization during the development of the hippocampus and prefrontal cortex, using more than 53,000 joint single-nucleus profiles of chromatin conformation and DNA methylation generated by single-nucleus methyl-3C sequencing (snm3C-seq3)3. The remodelling of DNA methylation is temporally separated from chromatin conformation dynamics. Using single-cell profiling and multimodal single-molecule imaging approaches, we have found that short-range chromatin interactions are enriched in neurons, whereas long-range interactions are enriched in glial cells and non-brain tissues. We reconstructed the regulatory programs of cell-type development and differentiation, finding putatively causal common variants for schizophrenia strongly overlapping with chromatin loop-connected, cell-type-specific regulatory regions. Our data provide multimodal resources for studying gene regulatory dynamics in brain development and demonstrate that single-cell three-dimensional multi-omics is a powerful approach for dissecting neuropsychiatric risk loci.

Superior temporal sulcus folding, functional network connectivity, and autistic-like traits in a non-clinical population

BACKGROUND

Autistic-like traits (ALT) are prevalent across the general population and might be linked to some facets of a broader autism spectrum disorder (ASD) phenotype. Recent studies suggest an association of these traits with both genetic and brain structural markers in non-autistic individuals, showing similar spatial location of findings observed in ASD and thus suggesting a potential neurobiological continuum.

METHODS

In this study, we first tested an association of ALTs (assessed with the AQ questionnaire) with cortical complexity, a cortical surface marker of early neurodevelopment, and then the association with disrupted functional connectivity. We analysed structural T1-weighted and resting-state functional MRI scans in 250 psychiatrically healthy individuals without a history of early developmental disorders, in a first step using the CAT12 toolbox for cortical complexity analysis and in a second step we used regional cortical complexity findings to apply the CONN toolbox for seed-based functional connectivity analysis.

RESULTS

Our findings show a significant negative correlation of both AQ total and AQ attention switching subscores with left superior temporal sulcus (STS) cortical folding complexity, with the former being significantly correlated with STS to left lateral occipital cortex connectivity, while the latter showed significant positive correlation of STS to left inferior/middle frontal gyrus connectivity (n = 233; all p < 0.05, FWE cluster-level corrected). Additional analyses also revealed a significant correlation of AQ attention to detail subscores with STS to left lateral occipital cortex connectivity.

LIMITATIONS

Phenotyping might affect association results (e.g. choice of inventories); in addition, our study was limited to subclinical expressions of autistic-like traits.

CONCLUSIONS

Our findings provide further evidence for biological correlates of ALT even in the absence of clinical ASD, while establishing a link between structural variation of early developmental origin and functional connectivity.

The Human Neural Cell Atlas of Zika Infection in developing human brain tissue: viral pathogenesis, innate immunity, and lineage reprogramming

Zika virus (ZIKV) infection during pregnancy can lead to fetal brain infection and developmental anomalies collectively known as congenital Zika syndrome (CZS). To define the molecular features underlying CZS in a relevant human cell model, we evaluated ZIKV infection and neurodevelopment in primary fetal brain explants and induced pluripotent stem cell-derived mixed neural cultures at single cell resolution. We identified astrocytes as key innate immune sentinel cells detecting ZIKV and producing IFN-β. In contrast, neural progenitor cells displayed impaired innate immunity and supported high levels of viral replication. ZIKV infection of neurons suppressed differentiation and synaptic signaling networks and programmed a molecular switch from neurogenesis to astrogliogenesis. We identified a universal ZIKV-driven cellular stress response linked to intrinsic apoptosis and regulated by IFN-β. These findings reveal how innate immune signaling intersects with ZIKV-driven perturbations in cellular function to influence CZS outcomes including neuron developmental dysfunction and apoptotic cell death.

Multi-omics identification of extracellular components of the fetal monkey and human neocortex

During development, precursor cells are continuously and intimately interacting with their extracellular environment, which guides their ability to generate functional tissues and organs. Much is known about the development of the neocortex in mammals. This information has largely been derived from histological analyses, heterochronic cell transplants, and genetic manipulations in mice, and to a lesser extent from transcriptomic and histological analyses in humans. However, these approaches have not led to a characterization of the extracellular composition of the developing neocortex in any species. Here, using a combination of single-cell transcriptomic analyses from published datasets, and our proteomics and immunohistofluorescence analyses, we provide a more comprehensive and unbiased picture of the early developing fetal neocortex in humans and non-human primates. Our findings provide a starting point for further hypothesis-driven studies on structural and signaling components in the developing cortex that had previously not been identified.

Single-cell m6A profiling in the mouse brain uncovers cell type-specific RNA methylomes and age-dependent differential methylation

N6-methyladenosine (m6A) is an abundant mRNA modification in the brain that has important roles in neurodevelopment and brain function. However, because of technical limitations, global profiling of m6A sites within the individual cell types that make up the brain has not been possible. Here, we develop a mouse model that enables transcriptome-wide m6A detection in any tissue of interest at single-cell resolution. We use these mice to map m6A across different brain regions and within single cells of the mouse cortex and discover a high degree of shared methylation across brain regions and cell types. However, we also identify a small number of differentially methylated mRNAs in neurons that encode important regulators of neuronal signaling, and we discover that microglia have lower levels of m6A than other cell types. Finally, we perform single-cell m6A mapping in aged mice and identify many transcripts with age-dependent changes in m6A.

RING1 missense variants reveal sensitivity of DNA damage repair to H2A monoubiquitination dosage during neurogenesis

Polycomb repressive complex 1 (PRC1) modifies chromatin through catalysis of histone H2A lysine 119 monoubiquitination (H2AK119ub1). RING1 and RNF2 interchangeably serve as the catalytic subunit within PRC1. Pathogenic missense variants in PRC1 core components reveal functions of these proteins that are obscured in knockout models. While Ring1a knockout models remain healthy, the microcephaly and neuropsychiatric phenotypes associated with a pathogenic RING1 missense variant implicate unappreciated functions. Using an in vitro model of neurodevelopment, we observe that RING1 contributes to the broad placement of H2AK119ub1, and that its targets overlap with those of RNF2. PRC1 complexes harboring hypomorphic RING1 bind target loci but do not catalyze H2AK119ub1, reducing H2AK119ub1 by preventing catalytically active complexes from accessing the locus. This results in delayed DNA damage repair and cell cycle progression in neural progenitor cells (NPCs). Conversely, reduced H2AK119ub1 due to hypomorphic RING1 does not generate differential expression that impacts NPC differentiation. In contrast, hypomorphic RNF2 generates a greater reduction in H2AK119ub1 that results in both delayed DNA repair and widespread transcriptional changes. These findings suggest that the DNA damage response is more sensitive to H2AK119ub1 dosage change than is regulation of gene expression.

Cross-site reproducibility of human cortical organoids reveals consistent cell type composition and architecture

While guided human cortical organoid (hCO) protocols reproducibly generate cortical cell types at one site, variability in hCO phenotypes across sites using a harmonized protocol has not yet been evaluated. To determine the cross-site reproducibility of hCO differentiation, three independent research groups assayed hCOs in multiple differentiation replicates from one induced pluripotent stem cell (iPSC) line using a harmonized miniaturized spinning bioreactor protocol across 3 months. hCOs were mostly cortical progenitor and neuronal cell types in reproducible proportions that were consistently organized in cortical wall-like buds. Cross-site differences were detected in hCO size and expression of metabolism and cellular stress genes. Variability in hCO phenotypes correlated with stem cell gene expression prior to differentiation and technical factors associated with seeding, suggesting iPSC quality and treatment are important for differentiation outcomes. Cross-site reproducibility of hCO cell type proportions and organization encourages future prospective meta-analytic studies modeling neurodevelopmental disorders in hCOs.

Genome-Wide Mendelian Randomization Identifies Ferroptosis-Related Drug Targets for Alzheimer's Disease

Background

Alzheimer's disease (AD) currently lacks effective disease-modifying treatments. Recent research suggests that ferroptosis could be a potential therapeutic target. Mendelian randomization (MR) is a widely used method for identifying novel therapeutic targets.

Objective

Employ genetic information to evaluate the causal impact of ferroptosis-related genes on the risk of AD.

Methods

564 ferroptosis-related genes were obtained from FerrDb. We derived genetic instrumental variables for these genes using four brain quantitative trait loci (QTL) and two blood QTL datasets. Summary-data-based Mendelian randomization (SMR) and two-sample MR methods were applied to estimate the causal effects of ferroptosis-related genes on AD. Using extern transcriptomic datasets and triple-transgenic mouse model of AD (3xTg-AD) to further validate the gene targets identified by the MR analysis.

Results

We identified 17 potential AD risk gene targets from GTEx, 13 from PsychENCODE, and 22 from BrainMeta (SMR p < 0.05 and HEIDI test p > 0.05). Six overlapping ferroptosis-related genes associated with AD were identified, which could serve as potential therapeutic targets (PEX10, CDC25A, EGFR, DLD, LIG3, and TRIB3). Additionally, we further pinpointed risk genes or proteins at the blood tissue and pQTL levels. Notably, EGFR demonstrated significant dysregulation in the extern transcriptomic datasets and 3xTg-AD models.

Conclusions

This study provides genetic evidence supporting the potential therapeutic benefits of targeting the six druggable genes for AD treatment, especially for EGFR (validated by transcriptome and 3xTg-AD), which could be useful for prioritizing AD drug development in the field of ferroptosis.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Oxygen-induced stress reveals context-specific gene regulatory effects in human brain organoids

The interaction between genetic variants and environmental stressors is key to understanding the mechanisms underlying neurological diseases. In this study, we used human brain organoids to explore how varying oxygen levels expose context-dependent gene regulatory effects. By subjecting a genetically diverse panel of 21 brain organoids to hypoxic and hyperoxic conditions, we identified thousands of gene regulatory changes that are undetectable under baseline conditions, with 1,745 trait-associated genes showing regulatory effects only in response to oxygen stress. To capture more nuanced transcriptional patterns, we employed topic modeling, which revealed context-specific gene regulation linked to dynamic cellular processes and environmental responses, offering a deeper understanding of how gene regulation is modulated in the brain. These findings underscore the importance of genotype-environment interactions in genetic studies of neurological disorders and provide new insights into the hidden regulatory mechanisms influenced by environmental factors in the brain.

Synaptic plasticity in human thalamocortical assembloids

Synaptic plasticities, such as long-term potentiation (LTP) and depression (LTD), tune synaptic efficacy and are essential for learning and memory. Current studies of synaptic plasticity in humans are limited by a lack of adequate human models. Here, we modeled the thalamocortical system by fusing human induced pluripotent stem cell-derived thalamic and cortical organoids. Single-nucleus RNA sequencing revealed that >80% of cells in thalamic organoids were glutamatergic neurons. When fused to form thalamocortical assembloids, thalamic and cortical organoids formed reciprocal long-range axonal projections and reciprocal synapses detectable by light and electron microscopy, respectively. Using whole-cell patch-clamp electrophysiology and two-photon imaging, we characterized glutamatergic synaptic transmission. Thalamocortical and corticothalamic synapses displayed short-term plasticity analogous to that in animal models. LTP and LTD were reliably induced at both synapses; however, their mechanisms differed from those previously described in rodents. Thus, thalamocortical assembloids provide a model system for exploring synaptic plasticity in human circuits.

Dynamic convergence of autism disorder risk genes across neurodevelopment

Over a hundred risk genes underlie risk for autism spectrum disorder (ASD) but the extent to which they converge on shared downstream targets to increase ASD risk is unknown. To test the hypothesis that cellular context impacts the nature of convergence, here we apply a pooled CRISPR approach to target 29 ASD loss-of-function genes in human induced pluripotent stem cell (hiPSC)-derived neural progenitor cells, glutamatergic neurons, and GABAergic neurons. Two distinct approaches (gene-level and network-level analyses) demonstrate that convergence is greatest in mature glutamatergic neurons. Convergent effects are dynamic, varying in strength, composition, and biological role between cell types, increasing with functional similarity of the ASD genes examined, and driven by cell-type-specific gene co-expression patterns. Stratification of ASD genes yield targeted drug predictions capable of reversing gene-specific convergent signatures in human cells and ASD-related behaviors in zebrafish. Altogether, convergent networks downstream of ASD risk genes represent novel points of individualized therapeutic intervention.

CellCover Captures Neural Stem Cell Progression in Mammalian Neocortical Development

Definition of cell classes across the tissues of living organisms is central in the analysis of growing atlases of single-cell RNA sequencing (scRNA-seq) data across biomedicine. Marker genes for cell classes are most often defined by differential expression (DE) methods that serially assess individual genes across landscapes of diverse cells. This serial approach has been extremely useful, but is limited because it ignores possible redundancy or complementarity across genes that can only be captured by analyzing multiple genes simultaneously. We aim to identify discriminating panels of genes. To efficiently explore the vast space of possible marker panels, leverage the large number of cells often sequenced, and overcome zero-inflation in scRNA-seq data, we propose viewing gene panel selection as a variation of the "minimal set-covering problem" in combinatorial optimization. We show that this new method, CellCover, captures cell-class-specific signals in the developing mouse neocortex that are distinct from those defined by DE methods. Transfer learning experiments across mouse, primate, and human data demonstrate that CellCover identifies markers of conserved cell classes in neurogenesis, as well as temporal progression in both progenitors and neurons. Exploring markers of human outer radial glia (oRG, or basal RG) across mammals, we show that transcriptomic elements of this key cell type in the expansion of the human cortex appeared in gliogenic precursors of the rodent before the full program emerged in the primate lineage. We have assembled the public datasets we use in this report at NeMO analytics where the expression of individual genes {NeMO Individual Genes} and marker gene panels can be freely explored {NeMO: Telley 3 Sets Covering Panels}, {NeMO: Telley 12 Sets Covering Panels}, and {NeMO: Sorted Brain Cell Covering Panels}. CellCover is available in {CellCover R} and {CellCover Python}.

A multi-layered integrative analysis reveals a cholesterol metabolic program in outer radial glia with implications for human brain evolution

The definition of molecular and cellular mechanisms contributing to brain ontogenetic trajectories is essential to investigate the evolution of our species. Yet their functional dissection at an appropriate level of granularity remains challenging. Capitalizing on recent efforts that have extensively profiled neural stem cells from the developing human cortex, we develop an integrative computational framework to perform trajectory inference and gene regulatory network reconstruction, (pseudo)time-informed non-negative matrix factorization for learning the dynamics of gene expression programs, and paleogenomic analysis for a higher-resolution mapping of derived regulatory variants in our species in comparison with our closest relatives. We provide evidence for cell type-specific regulation of gene expression programs during indirect neurogenesis. In particular, our analysis uncovers a key role for a cholesterol program in outer radial glia, regulated by zinc-finger transcription factor KLF6. A cartography of the regulatory landscape impacted by Homo sapiens-derived variants reveals signals of selection clustering around regulatory regions associated with GLI3, a well-known regulator of radial glial cell cycle, and impacting KLF6 regulation. Our study contributes to the evidence of significant changes in metabolic pathways in recent human brain evolution.

Cell type-specific dysregulation of gene expression due to Chd8 haploinsufficiency during mouse cortical development

Disruptive variants in the chromodomain helicase CHD8, which acts as a transcriptional regulator during neurodevelopment, are strongly associated with risk for autism spectrum disorder (ASD). Loss of CHD8 function is hypothesized to perturb gene regulatory networks in the developing brain, thereby contributing to ASD etiology. However, insight into the cell type-specific transcriptional effects of CHD8 loss of function remains limited. We used single-cell and single-nucleus RNA-sequencing to globally profile gene expression and identify dysregulated genes in the embryonic and juvenile wild type and Chd8 +/- mouse cortex, respectively. Chd8 and other ASD risk-associated genes showed a convergent expression trajectory that was largely conserved between the mouse and human developing cortex, increasing from the progenitor zones to the cortical plate. Genes associated with risk for neurodevelopmental disorders and genes involved in neuron projection development, chromatin remodeling, signaling, and migration were dysregulated in Chd8 +/- embryonic day (E) 12.5 radial glia. Genes implicated in synaptic organization and activity were dysregulated in Chd8 +/- postnatal day (P) 25 deep- and upper-layer excitatory cortical neurons, suggesting a delay in synaptic maturation or impaired synaptogenesis due to CHD8 loss of function. Our findings reveal a complex pattern of transcriptional dysregulation in Chd8 +/- developing cortex, potentially with distinct biological impacts on progenitors and maturing neurons in the excitatory neuronal lineage.

A single-cell transcriptomic landscape of cadmium-hindered brain development in mice

The effects of neurotoxicant cadmium (Cd) exposure on brain development have not been well elucidated. To investigate this, we have herein subjected pregnant mice to low-dose Cd throughout gestation. Using single-cell RNA sequencing (scRNA-seq), we explored the cellular responses in the embryonic brain to Cd exposure, and identified 18 distinct cell subpopulations that exhibited varied responses to Cd. Typically, Cd exposure impeded the development and maturation of cells in the brain, especially progenitor cells such as neural progenitor cells (NPCs) and oligodendrocyte progenitor cells (OPCs). It also caused significant cell subpopulation shifts in almost all the types of cells in the brain. Additionally, Cd exposure reduced the dendritic sophistication of cortical neurons in the offspring. Importantly, these changes led to aberrant Ca2+ activity in the cortex and neural behavior changes in mature offspring. These data contribute to our understanding of the effects and mechanisms of Cd exposure on brain development and highlight the importance of controlling environmental neurotoxicant exposure at the population level.

Age-dependent functional development pattern in neonatal brain: An fMRI-based brain entropy study

The relationship between brain entropy (BEN) and early brain development has been established through animal studies. However, it remains unclear whether the BEN can be used to identify age-dependent functional changes in human neonatal brains and the genetic underpinning of the new neuroimaging marker remains to be elucidated. In this study, we analyzed resting-state fMRI data from the Developing Human Connectome Project, including 280 infants who were scanned at 37.5-43.5 weeks postmenstrual age. The BEN maps were calculated for each subject, and a voxel-wise analysis was conducted using a general linear model to examine the effects of age, sex, and preterm birth on BEN. Additionally, we evaluated the correlation between regional BEN and gene expression levels. Our results demonstrated that the BEN in the sensorimotor-auditory and association cortices, along the 'S-A' axis, was significantly positively correlated with postnatal age (PNA), and negatively correlated with gestational age (GA), respectively. Meanwhile, the BEN in the right rolandic operculum correlated significantly with both GA and PNA. Preterm-born infants exhibited increased BEN values in widespread cortical areas, particularly in the visual-motor cortex, when compared to term-born infants. Moreover, we identified five BEN-related genes (DNAJC12, FIG4, STX12, CETN2, and IRF2BP2), which were involved in protein folding, synaptic vesicle transportation and cell division. These findings suggest that the fMRI-based BEN can serve as an indicator of age-dependent brain functional development in human neonates, which may be influenced by specific genes.

Development of the basic architecture of neocortical circuitry in the human fetus as revealed by the coupling spatiotemporal pattern of synaptogenesis along with microstructure and macroscale in vivo MR imaging

In humans, a quantifiable number of cortical synapses appears early in fetal life. In this paper, we present a bridge across different scales of resolution and the distribution of synapses across the transient cytoarchitectonic compartments: marginal zone (MZ), cortical plate (CP), subplate (SP), and in vivo MR images. The tissue of somatosensory cortex (7-26 postconceptional weeks (PCW)) was prepared for electron microscopy, and classified synapses with a determined subpial depth were used for creating histograms matched to the histological sections immunoreacted for synaptic markers and aligned to in vivo MR images (1.5 T) of corresponding fetal ages (maternal indication). Two time periods and laminar patterns of synaptogenesis were identified: an early and midfetal two-compartmental distribution (MZ and SP) and a late fetal three-compartmental distribution (CP synaptogenesis). During both periods, a voluminous, synapse-rich SP was visualized on the in vivo MR. Another novel finding concerns the phase of secondary expansion of the SP (13 PCW), where a quantifiable number of synapses appears in the upper SP. This lamina shows a T2 intermediate signal intensity below the low signal CP. In conclusion, the early fetal appearance of synapses shows early differentiation of putative genetic mechanisms underlying the synthesis, transport and assembly of synaptic proteins. "Pioneering" synapses are likely to play a morphogenetic role in constructing of fundamental circuitry architecture due to interaction between neurons. They underlie spontaneous, evoked, and resting state activity prior to ex utero experience. Synapses can also mediate genetic and environmental triggers, adversely altering the development of cortical circuitry and leading to neurodevelopmental disorders.

Molecular and cellular dynamics of the developing human neocortex at single-cell resolution

The development of the human neocortex is a highly dynamic process and involves complex cellular trajectories controlled by cell-type-specific gene regulation1. Here, we collected paired single-nucleus chromatin accessibility and transcriptome data from 38 human neocortical samples encompassing both the prefrontal cortex and primary visual cortex. These samples span five main developmental stages, ranging from the first trimester to adolescence. In parallel, we performed spatial transcriptomic analysis on a subset of the samples to illustrate spatial organization and intercellular communication. This atlas enables us to catalog cell type-, age-, and area-specific gene regulatory networks underlying neural differentiation. Moreover, combining single-cell profiling, progenitor purification, and lineage-tracing experiments, we have untangled the complex lineage relationships among progenitor subtypes during the transition from neurogenesis to gliogenesis in the human neocortex. We identified a tripotential intermediate progenitor subtype, termed Tri-IPC, responsible for the local production of GABAergic neurons, oligodendrocyte precursor cells, and astrocytes. Remarkably, most glioblastoma cells resemble Tri-IPCs at the transcriptomic level, suggesting that cancer cells hijack developmental processes to enhance growth and heterogeneity. Furthermore, by integrating our atlas data with large-scale GWAS data, we created a disease-risk map highlighting enriched ASD risk in second-trimester intratelencephalic projection neurons. Our study sheds light on the gene regulatory landscape and cellular dynamics of the developing human neocortex.

Examining the associations between microglia genetic capacity, early life exposures and white matter development at the level of the individual

There are inter-individual differences in susceptibility to the influence of early life experiences for which the underlying neurobiological mechanisms are poorly understood. Microglia play a role in environmental surveillance and may influence individual susceptibility to environmental factors. As an index of neurodevelopment, we estimated individual slopes of mean white matter fractional anisotropy (WM-FA) across three time-points (age 4.5, 6.0, and 7.5 years) for 351 participants. Individual variation in microglia reactivity was derived from an expression-based polygenic score(ePGS) comprised of Single Nucleotide Polymorphisms (SNPs) functionally related to the expression of microglia-enriched genes.A higher ePGS denotes an increased genetic capacity for the expression of microglia-related genes, and thus may confer a greater capacity to respond to the early environment and to influence brain development. We hypothesized that this ePGS would associate with the WM-FA index of neurodevelopment and moderate the influence of early environmental factors.Our findings show sex dependency, where a significant association between WM-FA and microglia ePGS was only obtained for females.We then examined associations with perinatal factors known to decrease (optimal birth outcomes and familial conditions) or increase (systemic inflammation) the risk for later mental health problems.In females, individuals with high microglia ePGS showed a negative association between systemic inflammation and WM-FA and a positive association between more advantageous environmental conditions and WM-FA. The microglia ePGS in females thus accounted for variations in the influence of the quality of the early environment on WM-FA.Finally, WM-FA slopes mediated the association of microglia ePGS with interpersonal problems and social hostility in females. Our findings suggest the genetic capacity for microglia function as a potential factor underlying differential susceptibility to early life exposuresthrough influences on neurodevelopment.

Microglial function interacts with the environment to affect sex-specific depression risk

There is a two-fold higher incidence of depression in females compared to men with recent studies suggesting a role for microglia in conferring this sex-dependent depression risk. In this study we investigated the nature of this relation. Using GWAS enrichment, gene-set enrichment analysis and Mendelian randomization, we found minimal evidence for a direct relation between genes functionally related to microglia and sex-dependent genetic risk for depression. We then used expression quantitative trait loci and single nucleus RNA-sequencing resources to generate polygenic scores (PGS) representative of individual variation in microglial function in the adult (UK Biobank; N = 54753-72682) and fetal (ALSPAC; N = 1452) periods. The adult microglial PGS moderated the association between BMI (UK Biobank; beta = 0.001, 95 %CI 0.0009 to 0.003, P = 7.74E-6) and financial insecurity (UK Biobank; beta = 0.001, 95 %CI 0.005 to 0.015, P = 2E-4) with depressive symptoms in females. The fetal microglia PGS moderated the association between maternal prenatal depressive symptoms and offspring depressive symptoms at 24 years in females (ALSPAC; beta = 0.04, 95 %CI 0.004 to 0.07, P = 0.03). We found no evidence for an interaction between the microglial PGS and depression risk factors in males. Our results illustrate a role for microglial function in the conferral of sex-dependent depression risk following exposure to a depression risk factor.

Brain Chimeroids reveal individual susceptibility to neurotoxic triggers

Interindividual genetic variation affects the susceptibility to and progression of many diseases1,2. However, efforts to study how individual human brains differ in normal development and disease phenotypes are limited by the paucity of faithful cellular human models, and the difficulty of scaling current systems to represent multiple people. Here we present human brain Chimeroids, a highly reproducible, multidonor human brain cortical organoid model generated by the co-development of cells from a panel of individual donors in a single organoid. By reaggregating cells from multiple single-donor organoids at the neural stem cell or neural progenitor cell stage, we generate Chimeroids in which each donor produces all cell lineages of the cerebral cortex, even when using pluripotent stem cell lines with notable growth biases. We used Chimeroids to investigate interindividual variation in the susceptibility to neurotoxic triggers that exhibit high clinical phenotypic variability: ethanol and the antiepileptic drug valproic acid. Individual donors varied in both the penetrance of the effect on target cell types, and the molecular phenotype within each affected cell type. Our results suggest that human genetic background may be an important mediator of neurotoxin susceptibility and introduce Chimeroids as a scalable system for high-throughput investigation of interindividual variation in processes of brain development and disease.

Five autism-associated transcriptional regulators target shared loci proximal to brain-expressed genes

Many autism spectrum disorder (ASD)-associated genes act as transcriptional regulators (TRs). Chromatin immunoprecipitation sequencing (ChIP-seq) was used to identify the regulatory targets of ARID1B, BCL11A, FOXP1, TBR1, and TCF7L2, ASD-associated TRs in the developing human and mouse cortex. These TRs shared substantial overlap in the binding sites, especially within open chromatin. The overlap within a promoter region, 1-2,000 bp upstream of the transcription start site, was highly predictive of brain-expressed genes. This signature was observed in 96 out of 102 ASD-associated genes. In vitro CRISPRi against ARID1B and TBR1 delineated downstream convergent biology in mouse cortical cultures. After 8 days, NeuN+ and CALB+ cells were decreased, GFAP+ cells were increased, and transcriptomic signatures correlated with the postmortem brain samples from individuals with ASD. We suggest that functional convergence across five ASD-associated TRs leads to shared neurodevelopmental outcomes of haploinsufficient disruption.

Deciphering the functional specialization of whole-brain spatiomolecular gradients in the adult brain

Cortical arealization arises during neurodevelopment from the confluence of molecular gradients representing patterned expression of morphogens and transcription factors. However, whether similar gradients are maintained in the adult brain remains unknown. Here, we uncover three axes of topographic variation in gene expression in the adult human brain that specifically capture previously identified rostral-caudal, dorsal-ventral, and medial-lateral axes of early developmental patterning. The interaction of these spatiomolecular gradients i) accurately reconstructs the position of brain tissue samples, ii) delineates known functional territories, and iii) can model the topographical variation of diverse cortical features. The spatiomolecular gradients are distinct from canonical cortical axes differentiating the primary sensory cortex from the association cortex, but radiate in parallel with the axes traversed by local field potentials along the cortex. We replicate all three molecular gradients in three independent human datasets as well as two nonhuman primate datasets and find that each gradient shows a distinct developmental trajectory across the lifespan. The gradients are composed of several well-known transcription factors (e.g., PAX6 and SIX3), and a small set of genes shared across gradients are strongly enriched for multiple diseases. Together, these results provide insight into the developmental sculpting of functionally distinct brain regions, governed by three robust transcriptomic axes embedded within brain parenchyma.

Gene regulatory landscape of cerebral cortex folding

Folding of the cerebral cortex is a key aspect of mammalian brain development and evolution, and defects are linked to severe neurological disorders. Primary folding occurs in highly stereotyped patterns that are predefined in the cortical germinal zones by a transcriptomic protomap. The gene regulatory landscape governing the emergence of this folding protomap remains unknown. We characterized the spatiotemporal dynamics of gene expression and active epigenetic landscape (H3K27ac) across prospective folds and fissures in ferret. Our results show that the transcriptomic protomap begins to emerge at early embryonic stages, and it involves cell-fate signaling pathways. The H3K27ac landscape reveals developmental cell-fate restriction and engages known developmental regulators, including the transcription factor Cux2. Manipulating Cux2 expression in cortical progenitors changed their proliferation and the folding pattern in ferret, caused by selective transcriptional changes as revealed by single-cell RNA sequencing analyses. Our findings highlight the key relevance of epigenetic mechanisms in defining the patterns of cerebral cortex folding.

Analyses of the autism-associated neuroligin-3 R451C mutation in human neurons reveal a gain-of-function synaptic mechanism

Mutations in many synaptic genes are associated with autism spectrum disorders (ASD), suggesting that synaptic dysfunction is a key driver of ASD pathogenesis. Among these mutations, the R451C substitution in the NLGN3 gene that encodes the postsynaptic adhesion molecule Neuroligin-3 is noteworthy because it was the first specific mutation linked to ASDs. In mice, the corresponding Nlgn3 R451C-knockin mutation recapitulates social interaction deficits of ASD patients and produces synaptic abnormalities, but the impact of the NLGN3 R451C mutation on human neurons has not been investigated. Here, we generated human knockin neurons with the NLGN3 R451C and NLGN3 null mutations. Strikingly, analyses of NLGN3 R451C-mutant neurons revealed that the R451C mutation decreased NLGN3 protein levels but enhanced the strength of excitatory synapses without affecting inhibitory synapses; meanwhile NLGN3 knockout neurons showed reduction in excitatory synaptic strengths. Moreover, overexpression of NLGN3 R451C recapitulated the synaptic enhancement in human neurons. Notably, the augmentation of excitatory transmission was confirmed in vivo with human neurons transplanted into mouse forebrain. Using single-cell RNA-seq experiments with co-cultured excitatory and inhibitory NLGN3 R451C-mutant neurons, we identified differentially expressed genes in relatively mature human neurons corresponding to synaptic gene expression networks. Moreover, gene ontology and enrichment analyses revealed convergent gene networks associated with ASDs and other mental disorders. Our findings suggest that the NLGN3 R451C mutation induces a gain-of-function enhancement in excitatory synaptic transmission that may contribute to the pathophysiology of ASD.

Functional neurogenomics in autism spectrum disorders: A decade of progress

Advances in autism spectrum disorder (ASD) genetics have identified many genetic causes, reflecting remarkable progress while at the same time identifying challenges such as heterogeneity and pleiotropy, which complicate attempts to connect genetic risk to mechanisms. High-throughput functional genomic approaches have yielded progress by defining a molecular pathology in the brain of individuals with ASD and in identifying convergent biological pathways through which risk genes are predicted to act. These studies indicate that ASD genetic risk converges in early brain development, primarily during the period of cortical neurogenesis. Over development, genetic perturbations in turn lead to broad neuronal signaling dysregulation, most prominent in glutamatergic cortical-cortical projecting neurons and somatostatin positive interneurons, which is accompanied by glial dyshomeostasis throughout the cerebral cortex. Connecting these developmental perturbations to disrupted neuronal and glial function in the postnatal brain is an important direction in current research. Coupling functional genomic approaches with advances in induced pluripotent stem cell-derived neural organoid development provides a promising avenue for connecting brain pathology to developmental mechanisms.

Autism gene variants disrupt enteric neuron migration and cause gastrointestinal dysmotility

The comorbidity of autism spectrum disorders and severe gastrointestinal symptoms is well-established, yet the molecular underpinnings remain unknown. The identification of high-confidence large-effect autism risk genes offers the opportunity to identify convergent, underlying biology by studying these genes in the context of the gastrointestinal system. Here we show that the expression of these genes is enriched in human prenatal gut neurons as well as their migratory progenitors, suggesting that the development and/or function of these neurons may be disrupted by autism-associated pathogenic variants, leading to gastrointestinal dysfunction. Here we document the prevalence of gastrointestinal issues in patients with large-effect variants in sixteen of these genes, highlighting dysmotility, consistent with potential enteric neuron dysfunction. Using the high-throughput diploid frog Xenopus tropicalis , we individually target five of these genes ( SYNGAP1, CHD8, SCN2A, CHD2 , and DYRK1A ) and observe disrupted enteric neuronal progenitor migration for each. More extensive analysis of DYRK1A reveals that perturbation causes gut dysmotility in vivo , which can be ameliorated by treatment with a selective serotonin reuptake inhibitor (escitalopram) or a serotonin receptor 6 agonist, identified by in vivo drug screening. This work suggests that atypical development of enteric neurons contributes to the gastrointestinal distress commonly seen in individuals with autism and that increasing serotonin signaling may be a productive therapeutic avenue.

Using a comprehensive atlas and predictive models to reveal the complexity and evolution of brain-active regulatory elements

Most genetic variants associated with psychiatric disorders are located in noncoding regions of the genome. To investigate their functional implications, we integrate epigenetic data from the PsychENCODE Consortium and other published sources to construct a comprehensive atlas of candidate brain cis-regulatory elements. Using deep learning, we model these elements' sequence syntax and predict how binding sites for lineage-specific transcription factors contribute to cell type-specific gene regulation in various types of glia and neurons. The elements' evolutionary history suggests that new regulatory information in the brain emerges primarily via smaller sequence mutations within conserved mammalian elements rather than entirely new human- or primate-specific sequences. However, primate-specific candidate elements, particularly those active during fetal brain development and in excitatory neurons and astrocytes, are implicated in the heritability of brain-related human traits. Additionally, we introduce PsychSCREEN, a web-based platform offering interactive visualization of PsychENCODE-generated genetic and epigenetic data from diverse brain cell types in individuals with psychiatric disorders and healthy controls.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Cross-ancestry atlas of gene, isoform, and splicing regulation in the developing human brain

Neuropsychiatric genome-wide association studies (GWASs), including those for autism spectrum disorder and schizophrenia, show strong enrichment for regulatory elements in the developing brain. However, prioritizing risk genes and mechanisms is challenging without a unified regulatory atlas. Across 672 diverse developing human brains, we identified 15,752 genes harboring gene, isoform, and/or splicing quantitative trait loci, mapping 3739 to cellular contexts. Gene expression heritability drops during development, likely reflecting both increasing cellular heterogeneity and the intrinsic properties of neuronal maturation. Isoform-level regulation, particularly in the second trimester, mediated the largest proportion of GWAS heritability. Through colocalization, we prioritized mechanisms for about 60% of GWAS loci across five disorders, exceeding adult brain findings. Finally, we contextualized results within gene and isoform coexpression networks, revealing the comprehensive landscape of transcriptome regulation in development and disease.

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.

CellWalker2: multi-omic discovery of hierarchical cell type relationships and their associations with genomic annotations

CellWalker2 is a graph diffusion-based method for single-cell genomics data integration. It extends the CellWalker model by incorporating hierarchical relationships between cell types, providing estimates of statistical significance, and adding data structures for analyzing multi-omics data so that gene expression and open chromatin can be jointly modeled. Our open-source software enables users to annotate cells using existing ontologies and to probabilistically match cell types between two or more contexts, including across species. CellWalker2 can also map genomic regions to cell ontologies, enabling precise annotation of elements derived from bulk data, such as enhancers, genetic variants, and sequence motifs. Through simulation studies, we show that CellWalker2 performs better than existing methods in cell type annotation and mapping. We then use data from the brain and immune system to demonstrate CellWalker2's ability to discover cell type-specific regulatory programs and both conserved and divergent cell type relationships in complex tissues.

Mutant IDH inhibitors induce lineage differentiation in IDH-mutant oligodendroglioma

A subset of patients with IDH-mutant glioma respond to inhibitors of mutant IDH (IDHi), yet the molecular underpinnings of such responses are not understood. Here, we profiled by single-cell or single-nucleus RNA-sequencing three IDH-mutant oligodendrogliomas from patients who derived clinical benefit from IDHi. Importantly, the tissues were sampled on-drug, four weeks from treatment initiation. We further integrate our findings with analysis of single-cell and bulk transcriptomes from independent cohorts and experimental models. We find that IDHi treatment induces a robust differentiation toward the astrocytic lineage, accompanied by a depletion of stem-like cells and a reduction of cell proliferation. Furthermore, mutations in NOTCH1 are associated with decreased astrocytic differentiation and may limit the response to IDHi. Our study highlights the differentiating potential of IDHi on the cellular hierarchies that drive oligodendrogliomas and suggests a genetic modifier that may improve patient stratification.

Integrated electrophysiological and genomic profiles of single cells reveal spiking tumor cells in human glioma

Prior studies have described the complex interplay that exists between glioma cells and neurons, however, the electrophysiological properties endogenous to tumor cells remain obscure. To address this, we employed Patch-sequencing on human glioma specimens and found that one third of patched cells in IDH mutant (IDH mut ) tumors demonstrate properties of both neurons and glia by firing single, short action potentials. To define these hybrid cells (HCs) and discern if they are tumor in origin, we developed a computational tool, Single Cell Rule Association Mining (SCRAM), to annotate each cell individually. SCRAM revealed that HCs represent tumor and non-tumor cells that feature GABAergic neuron and oligodendrocyte precursor cell signatures. These studies are the first to characterize the combined electrophysiological and molecular properties of human glioma cells and describe a new cell type in human glioma with unique electrophysiological and transcriptomic properties that are likely also present in the non-tumor mammalian brain.

De novo variants in SP9 cause a novel form of interneuronopathy characterized by intellectual disability, autism spectrum disorder, and epilepsy with variable expressivity

PURPOSE

Interneuronopathies are a group of neurodevelopmental disorders characterized by deficient migration and differentiation of gamma-aminobutyric acidergic interneurons resulting in a broad clinical spectrum, including autism spectrum disorders, early-onset epileptic encephalopathy, intellectual disability, and schizophrenic disorders. SP9 is a transcription factor belonging to the Krüppel-like factor and specificity protein family, the members of which harbor highly conserved DNA-binding domains. SP9 plays a central role in interneuron development and tangential migration, but it has not yet been implicated in a human neurodevelopmental disorder.

METHODS

Cases with SP9 variants were collected through international data-sharing networks. To address the specific impact of SP9 variants, in silico and in vitro assays were carried out.

RESULTS

De novo heterozygous variants in SP9 cause a novel form of interneuronopathy. SP9 missense variants affecting the glutamate 378 amino acid result in severe epileptic encephalopathy because of hypomorphic and neomorphic DNA-binding effects, whereas SP9 loss-of-function variants result in a milder phenotype with epilepsy, developmental delay, and autism spectrum disorder.

CONCLUSION

De novo heterozygous SP9 variants are responsible for a neurodevelopmental disease. Interestingly, variants located in conserved DNA-binding domains of KLF/SP family transcription factors may lead to neomorphic DNA-binding functions resulting in a combination of loss- and gain-of-function effects.

Functional synergy of a human-specific and an ape-specific metabolic regulator in human neocortex development

Metabolism has recently emerged as a major target of genes implicated in the evolutionary expansion of human neocortex. One such gene is the human-specific gene ARHGAP11B. During human neocortex development, ARHGAP11B increases the abundance of basal radial glia, key progenitors for neocortex expansion, by stimulating glutaminolysis (glutamine-to-glutamate-to-alpha-ketoglutarate) in mitochondria. Here we show that the ape-specific protein GLUD2 (glutamate dehydrogenase 2), which also operates in mitochondria and converts glutamate-to-αKG, enhances ARHGAP11B's ability to increase basal radial glia abundance. ARHGAP11B + GLUD2 double-transgenic bRG show increased production of aspartate, a metabolite essential for cell proliferation, from glutamate via alpha-ketoglutarate and the TCA cycle. Hence, during human evolution, a human-specific gene exploited the existence of another gene that emerged during ape evolution, to increase, via concerted changes in metabolism, progenitor abundance and neocortex size.

Spatiotemporal expression of thyroid hormone transporter MCT8 and THRA mRNA in human cerebral organoids recapitulating first trimester cortex development

Thyroid hormones (TH) play critical roles during nervous system development and patients carrying coding variants of MCT8 (monocarboxylate transporter 8) or THRA (thyroid hormone receptor alpha) present a spectrum of neurological phenotypes resulting from perturbed local TH action during early brain development. Recently, human cerebral organoids (hCOs) emerged as powerful in vitro tools for disease modelling recapitulating key aspects of early human cortex development. To begin exploring prospects of this model for thyroid research, we performed a detailed characterization of the spatiotemporal expression of MCT8 and THRA in developing hCOs. Immunostaining showed MCT8 membrane expression in neuronal progenitor cell types including early neuroepithelial cells, radial glia cells (RGCs), intermediate progenitors and outer RGCs. In addition, we detected robust MCT8 protein expression in deep layer and upper layer neurons. Spatiotemporal SLC16A2 mRNA expression, detected by fluorescent in situ hybridization (FISH), was highly concordant with MCT8 protein expression across cortical cell layers. FISH detected THRA mRNA expression already in neuroepithelium before the onset of neurogenesis. THRA mRNA expression remained low in the ventricular zone, increased in the subventricular zone whereas strong THRA expression was observed in excitatory neurons. In combination with a robust up-regulation of known T3 response genes following T3 treatment, these observations show that hCOs provide a promising and experimentally tractable model to probe local TH action during human cortical neurogenesis and eventually to model the consequences of impaired TH function for early cortex development.

Multimodal Nature of the Single-cell Primate Brain Atlas: Morphology, Transcriptome, Electrophysiology, and Connectivity

Primates exhibit complex brain structures that augment cognitive function. The neocortex fulfills high-cognitive functions through billions of connected neurons. These neurons have distinct transcriptomic, morphological, and electrophysiological properties, and their connectivity principles vary. These features endow the primate brain atlas with a multimodal nature. The recent integration of next-generation sequencing with modified patch-clamp techniques is revolutionizing the way to census the primate neocortex, enabling a multimodal neuronal atlas to be established in great detail: (1) single-cell/single-nucleus RNA-seq technology establishes high-throughput transcriptomic references, covering all major transcriptomic cell types; (2) patch-seq links the morphological and electrophysiological features to the transcriptomic reference; (3) multicell patch-clamp delineates the principles of local connectivity. Here, we review the applications of these technologies in the primate neocortex and discuss the current advances and tentative gaps for a comprehensive understanding of the primate neocortex.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Chronic hypoxia remodels the tumor microenvironment to support glioma stem cell growth

Cerebral organoids co-cultured with patient derived glioma stem cells (GLICOs) are an experimentally tractable research tool useful for investigating the role of the human brain tumor microenvironment in glioblastoma. Here we describe long-term GLICOs, a novel model in which COs are grown from embryonic stem cell cultures containing low levels of GSCs and tumor development is monitored over extended durations (ltGLICOs). Single-cell profiling of ltGLICOs revealed an unexpectedly long latency period prior to GSC expansion, and that normal organoid development was unimpaired by the presence of low numbers of GSCs. However, as organoids age they experience chronic hypoxia and oxidative stress which remodels the tumor microenvironment to promote GSC expansion. Receptor-ligand modelling identified astrocytes, which secreted various pro-tumorigenic ligands including FGF1, as the primary cell type for GSC crosstalk and single-cell multi-omic analysis revealed these astrocytes were under the control of ischemic regulatory networks. Functional validation confirmed hypoxia as a driver of pro-tumorigenic astrocytic ligand secretion and that GSC expansion was accelerated by pharmacological induction of oxidative stress. When controlled for genotype, the close association between glioma aggressiveness and patient age has very few proposed biological explanations. Our findings indicate that age-associated increases in cerebral vascular insufficiency and associated regional chronic cerebral hypoxia may contribute to this phenomenon.

Synaptic plasticity in human thalamocortical assembloids

Synaptic plasticities, such as long-term potentiation (LTP) and depression (LTD), tune synaptic efficacy and are essential for learning and memory. Current studies of synaptic plasticity in humans are limited by a lack of adequate human models. Here, we modeled the thalamocortical system by fusing human induced pluripotent stem cell-derived thalamic and cortical organoids. Single-nucleus RNA-sequencing revealed that most cells in mature thalamic organoids were glutamatergic neurons. When fused to form thalamocortical assembloids, thalamic and cortical organoids formed reciprocal long-range axonal projections and reciprocal synapses detectable by light and electron microscopy, respectively. Using whole-cell patch-clamp electrophysiology and two-photon imaging, we characterized glutamatergic synaptic transmission. Thalamocortical and corticothalamic synapses displayed short-term plasticity analogous to that in animal models. LTP and LTD were reliably induced at both synapses; however, their mechanisms differed from those previously described in rodents. Thus, thalamocortical assembloids provide a model system for exploring synaptic plasticity in human circuits.