Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex

Decoding the spatiotemporal regulation of transcription factors during human spinal cord development

The spinal cord is a crucial component of the central nervous system that facilitates sensory processing and motor performance. Despite its importance, the spatiotemporal codes underlying human spinal cord development have remained elusive. In this study, we have introduced an image-based single-cell transcription factor (TF) expression decoding spatial transcriptome method (TF-seqFISH) to investigate the spatial expression and regulation of TFs during human spinal cord development. By combining spatial transcriptomic data from TF-seqFISH and single-cell RNA-sequencing data, we uncovered the spatial distribution of neural progenitor cells characterized by combinatorial TFs along the dorsoventral axis, as well as the molecular and spatial features governing neuronal generation, migration, and differentiation along the mediolateral axis. Notably, we observed a sandwich-like organization of excitatory and inhibitory interneurons transiently appearing in the dorsal horns of the developing human spinal cord. In addition, we integrated data from 10× Visium to identify early and late waves of neurogenesis in the dorsal horn, revealing the formation of laminas in the dorsal horns. Our study also illuminated the spatial differences and molecular cues underlying motor neuron (MN) diversification, and the enrichment of Amyotrophic Lateral Sclerosis (ALS) risk genes in MNs and microglia. Interestingly, we detected disease-associated microglia (DAM)-like microglia groups in the developing human spinal cord, which are predicted to be vulnerable to ALS and engaged in the TYROBP causal network and response to unfolded proteins. These findings provide spatiotemporal transcriptomic resources on the developing human spinal cord and potential strategies for spinal cord injury repair and ALS treatment.

Temporal transcriptomic dynamics in developing macaque neocortex

Despite intense research on mice, the transcriptional regulation of neocortical neurogenesis remains limited in humans and non-human primates. Cortical development in rhesus macaque is known to recapitulate multiple facets of cortical development in humans, including the complex composition of neural stem cells and the thicker supragranular layer. To characterize temporal shifts in transcriptomic programming responsible for differentiation from stem cells to neurons, we sampled parietal lobes of rhesus macaque at E40, E50, E70, E80, and E90, spanning the full period of prenatal neurogenesis. Single-cell RNA sequencing produced a transcriptomic atlas of developing parietal lobe in rhesus macaque neocortex. Identification of distinct cell types and neural stem cells emerging in different developmental stages revealed a terminally bifurcating trajectory from stem cells to neurons. Notably, deep-layer neurons appear in the early stages of neurogenesis, while upper-layer neurons appear later. While these different lineages show overlap in their differentiation program, cell fates are determined post-mitotically. Trajectories analysis from ventricular radial glia (vRGs) to outer radial glia (oRGs) revealed dynamic gene expression profiles and identified differential activation of BMP, FGF, and WNT signaling pathways between vRGs and oRGs. These results provide a comprehensive overview of the temporal patterns of gene expression leading to different fates of radial glial progenitors during neocortex layer formation.

Variable expression of MECP2, CDKL5, and FMR1 in the human brain: Implications for gene restorative therapies

MECP2, CDKL5, and FMR1 are three X-linked neurodevelopmental genes associated with Rett, CDKL5-, and fragile-X syndrome, respectively. These syndromes are characterized by distinct constellations of severe cognitive and neurobehavioral anomalies, reflecting the broad but unique expression patterns of each of the genes in the brain. As these disorders are not thought to be neurodegenerative and may be reversible, a major goal has been to restore expression of the functional proteins in the patient's brain. Strategies have included gene therapy, gene editing, and selective Xi-reactivation methodologies. However, tissue penetration and overall delivery to various regions of the brain remain challenging for each strategy. Thus, gaining insights into how much restoration would be required and what regions/cell types in the brain must be targeted for meaningful physiological improvement would be valuable. As a step toward addressing these questions, here we perform a meta-analysis of single-cell transcriptomics data from the human brain across multiple developmental stages, in various brain regions, and in multiple donors. We observe a substantial degree of expression variability for MECP2, CDKL5, and FMR1 not only across cell types but also between donors. The wide range of expression may help define a therapeutic window, with the low end delineating a minimum level required to restore physiological function and the high end informing toxicology margin. Finally, the inter-cellular and inter-individual variability enable identification of co-varying genes and will facilitate future identification of biomarkers.

Development and evolution of the primate neocortex from a progenitor cell perspective

The generation of neurons in the developing neocortex is a major determinant of neocortex size. Crucially, the increase in cortical neuron numbers in the primate lineage, notably in the upper-layer neurons, contributes to increased cognitive abilities. Here, we review major evolutionary changes affecting the apical progenitors in the ventricular zone and focus on the key germinal zone constituting the foundation of neocortical neurogenesis in primates, the outer subventricular zone (OSVZ). We summarize characteristic features of the OSVZ and its key stem cell type, the basal (or outer) radial glia. Next, we concentrate on primate-specific and human-specific genes, expressed in OSVZ-progenitors, the ability of which to amplify these progenitors by targeting the regulation of the cell cycle ultimately underlies the evolutionary increase in upper-layer neurons. Finally, we address likely differences in neocortical development between present-day humans and Neanderthals that are based on human-specific amino acid substitutions in proteins operating in cortical progenitors.

Loss of Ezh2 in the medial ganglionic eminence alters interneuron fate, cell morphology and gene expression profiles

Introduction

Enhancer of zeste homolog 2 (Ezh2) is responsible for trimethylation of histone 3 at lysine 27 (H3K27me3), resulting in repression of gene expression. Here, we explore the role of Ezh2 in forebrain GABAergic interneuron development.

Methods

We removed Ezh2 in the MGE by generating Nkx2-1Cre;Ezh2 conditional knockout mice. We then characterized changes in MGE-derived interneuron fate and electrophysiological properties in juvenile mice, as well as alterations in gene expression, chromatin accessibility and histone modifications in the MGE.

Results

Loss of Ezh2 increases somatostatin-expressing (SST+) and decreases parvalbumin-expressing (PV+) interneurons in the forebrain. We observe fewer MGE-derived interneurons in the first postnatal week, indicating reduced interneuron production. Intrinsic electrophysiological properties in SST+ and PV+ interneurons are normal, but PV+ interneurons display increased axonal complexity in Ezh2 mutant mice. Single nuclei multiome analysis revealed differential gene expression patterns in the embryonic MGE that are predictive of these cell fate changes. Lastly, CUT&Tag analysis revealed that some genomic loci are particularly resistant or susceptible to shifts in H3K27me3 levels in the absence of Ezh2, indicating differential selectivity to epigenetic perturbation.

Discussion

Thus, loss of Ezh2 in the MGE alters interneuron fate, morphology, and gene expression and regulation. These findings have important implications for both normal development and potentially in disease etiologies.

Molecular signatures of cortical expansion in the human fetal brain

The third trimester of human gestation is characterised by rapid increases in brain volume and cortical surface area. A growing catalogue of cells in the prenatal brain has revealed remarkable molecular diversity across cortical areas.1,2 Despite this, little is known about how this translates into the patterns of differential cortical expansion observed in humans during the latter stages of gestation. Here we present a new resource, μBrain, to facilitate knowledge translation between molecular and anatomical descriptions of the prenatal developing brain. Built using generative artificial intelligence, μBrain is a three-dimensional cellular-resolution digital atlas combining publicly-available serial sections of the postmortem human brain at 21 weeks gestation3 with bulk tissue microarray data, sampled across 29 cortical regions and 5 transient tissue zones.4 Using μBrain, we evaluate the molecular signatures of preferentially-expanded cortical regions during human gestation, quantified in utero using magnetic resonance imaging (MRI). We find that differences in the rates of expansion across cortical areas during gestation respect anatomical and evolutionary boundaries between cortical types5 and are founded upon extended periods of upper-layer cortical neuron migration that continue beyond mid-gestation. We identify a set of genes that are upregulated from mid-gestation and highly expressed in rapidly expanding neocortex, which are implicated in genetic disorders with cognitive sequelae. Our findings demonstrate a spatial coupling between areal differences in the timing of neurogenesis and rates of expansion across the neocortical sheet during the prenatal epoch. The μBrain atlas is available from: https://garedaba.github.io/micro-brain/ and provides a new tool to comprehensively map early brain development across domains, model systems and resolution scales.

A foundational atlas of autism protein interactions reveals molecular convergence

Translating high-confidence (hc) autism spectrum disorder (ASD) genes into viable treatment targets remains elusive. We constructed a foundational protein-protein interaction (PPI) network in HEK293T cells involving 100 hcASD risk genes, revealing over 1,800 PPIs (87% novel). Interactors, expressed in the human brain and enriched for ASD but not schizophrenia genetic risk, converged on protein complexes involved in neurogenesis, tubulin biology, transcriptional regulation, and chromatin modification. A PPI map of 54 patient-derived missense variants identified differential physical interactions, and we leveraged AlphaFold-Multimer predictions to prioritize direct PPIs and specific variants for interrogation in Xenopus tropicalis and human forebrain organoids. A mutation in the transcription factor FOXP1 led to reconfiguration of DNA binding sites and altered development of deep cortical layer neurons in forebrain organoids. This work offers new insights into molecular mechanisms underlying ASD and describes a powerful platform to develop and test therapeutic strategies for many genetically-defined conditions.

IDHwt glioblastomas can be stratified by their transcriptional response to standard treatment, with implications for targeted therapy

BACKGROUND

Glioblastoma (GBM) brain tumors lacking IDH1 mutations (IDHwt) have the worst prognosis of all brain neoplasms. Patients receive surgery and chemoradiotherapy but tumors almost always fatally recur.

RESULTS

Using RNA sequencing data from 107 pairs of pre- and post-standard treatment locally recurrent IDHwt GBM tumors, we identify two responder subtypes based on longitudinal changes in gene expression. In two thirds of patients, a specific subset of genes is upregulated from primary to recurrence (Up responders), and in one third, the same genes are downregulated (Down responders), specifically in neoplastic cells. Characterization of the responder subtypes indicates subtype-specific adaptive treatment resistance mechanisms that are associated with distinct changes in the tumor microenvironment. In Up responders, recurrent tumors are enriched in quiescent proneural GBM stem cells and differentiated neoplastic cells, with increased interaction with the surrounding normal brain and neurotransmitter signaling, whereas Down responders commonly undergo mesenchymal transition. ChIP-sequencing data from longitudinal GBM tumors suggests that the observed transcriptional reprogramming could be driven by Polycomb-based chromatin remodeling rather than DNA methylation.

CONCLUSIONS

We show that the responder subtype is cancer-cell intrinsic, recapitulated in in vitro GBM cell models, and influenced by the presence of the tumor microenvironment. Stratifying GBM tumors by responder subtype may lead to more effective treatment.

Human fetal brain self-organizes into long-term expanding organoids

Human brain development involves an orchestrated, massive neural progenitor expansion while a multi-cellular tissue architecture is established. Continuously expanding organoids can be grown directly from multiple somatic tissues, yet to date, brain organoids can solely be established from pluripotent stem cells. Here, we show that healthy human fetal brain in vitro self-organizes into organoids (FeBOs), phenocopying aspects of in vivo cellular heterogeneity and complex organization. FeBOs can be expanded over long time periods. FeBO growth requires maintenance of tissue integrity, which ensures production of a tissue-like extracellular matrix (ECM) niche, ultimately endowing FeBO expansion. FeBO lines derived from different areas of the central nervous system (CNS), including dorsal and ventral forebrain, preserve their regional identity and allow to probe aspects of positional identity. Using CRISPR-Cas9, we showcase the generation of syngeneic mutant FeBO lines for the study of brain cancer. Taken together, FeBOs constitute a complementary CNS organoid platform.

On the ever-growing functional versatility of the CRISPR-Cas13 system

CRISPR-Cas systems evolved in prokaryotes to implement a powerful antiviral immune response as a result of sequence-specific targeting by ribonucleoproteins. One of such systems consists of an RNA-guided RNA endonuclease, known as CRISPR-Cas13. In very recent years, this system is being repurposed in different ways in order to decipher and engineer gene expression programmes. Here, we discuss the functional versatility of the CRISPR-Cas13 system, which includes the ability for RNA silencing, RNA editing, RNA tracking, nucleic acid detection and translation regulation. This functional palette makes the CRISPR-Cas13 system a relevant tool in the broad field of systems and synthetic biology.

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

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

iPSC-derived models of PACS1 syndrome reveal transcriptional and functional deficits in neuron activity

PACS1 syndrome is a neurodevelopmental disorder characterized by intellectual disability and distinct craniofacial abnormalities resulting from a de novo p.R203W variant in phosphofurin acidic cluster sorting protein 1 (PACS1). PACS1 is known to have functions in the endosomal pathway and nucleus, but how the p.R203W variant affects developing neurons is not fully understood. Here we differentiated stem cells towards neuronal models including cortical organoids to investigate the impact of the PACS1 syndrome-causing variant on neurodevelopment. While few deleterious effects were detected in PACS1(+/R203W) neural precursors, mature PACS1(+/R203W) glutamatergic neurons exhibited impaired expression of genes involved in synaptic signaling processes. Subsequent characterization of neural activity using calcium imaging and multielectrode arrays revealed the p.R203W PACS1 variant leads to a prolonged neuronal network burst duration mediated by an increased interspike interval. These findings demonstrate the impact of the PACS1 p.R203W variant on developing human neural tissue and uncover putative electrophysiological underpinnings of disease.

Patient derived model of UBA5-associated encephalopathy identifies defects in neurodevelopment and highlights potential therapies

UBA5 encodes for the E1 enzyme of the UFMylation cascade, which plays an essential role in ER homeostasis. The clinical phenotypes of UBA5-associated encephalopathy include developmental delays, epilepsy and intellectual disability. To date, there is no humanized neuronal model to study the cellular and molecular consequences of UBA5 pathogenic variants. We developed and characterized patient-derived cortical organoid cultures and identified defects in GABAergic interneuron development. We demonstrated aberrant neuronal firing and microcephaly phenotypes in patient-derived organoids. Mechanistically, we show that ER homeostasis is perturbed along with exacerbated unfolded protein response pathway in cells and organoids expressing UBA5 pathogenic variants. We also assessed two gene expression modalities that augmented UBA5 expression to rescue aberrant molecular and cellular phenotypes. Our study provides a novel humanized model that allows further investigations of UBA5 variants in the brain and highlights novel systemic approaches to alleviate cellular aberrations for this rare, developmental disorder.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Proteomic Changes in the Human Cerebrovasculature in Alzheimer's Disease and Related Tauopathies Linked to Peripheral Biomarkers in Plasma and Cerebrospinal Fluid

Dysfunction of the neurovascular unit stands as a significant pathological hallmark of Alzheimer's disease (AD) and age-related neurodegenerative diseases. Nevertheless, detecting vascular changes in the brain within bulk tissues has proven challenging, limiting our ability to characterize proteomic alterations from less abundant cell types. To address this challenge, we conducted quantitative proteomic analyses on both bulk brain tissues and cerebrovascular-enriched fractions from the same individuals, encompassing cognitively unimpaired control, progressive supranuclear palsy (PSP), and AD cases. Protein co-expression network analysis identified modules unique to the cerebrovascular fractions, specifically enriched with pericytes, endothelial cells, and smooth muscle cells. Many of these modules also exhibited significant correlations with amyloid plaques, cerebral amyloid angiopathy (CAA), and/or tau pathology in the brain. Notably, the protein products within AD genetic risk loci were found concentrated within modules unique to the vascular fractions, consistent with a role of cerebrovascular deficits in the etiology of AD. To prioritize peripheral AD biomarkers associated with vascular dysfunction, we assessed the overlap between differentially abundant proteins in AD cerebrospinal fluid (CSF) and plasma with a vascular-enriched network modules in the brain. This analysis highlighted matrisome proteins, SMOC1 and SMOC2, as being increased in CSF, plasma, and brain. Immunohistochemical analysis revealed SMOC1 deposition in both parenchymal plaques and CAA in the AD brain, whereas SMOC2 was predominantly localized to CAA. Collectively, these findings significantly enhance our understanding of the involvement of cerebrovascular abnormalities in AD, shedding light on potential biomarkers and molecular pathways associated with CAA and vascular dysfunction in neurodegenerative diseases.

STAB2: an updated spatio-temporal cell atlas of the human and mouse brain

The brain is constituted of heterogeneous types of neuronal and non-neuronal cells, which are organized into distinct anatomical regions, and show precise regulation of gene expression during development, aging and function. In the current database release, STAB2 provides a systematic cellular map of the human and mouse brain by integrating recently published large-scale single-cell and single-nucleus RNA-sequencing datasets from diverse regions and across lifespan. We applied a hierarchical strategy of unsupervised clustering on the integrated single-cell transcriptomic datasets to precisely annotate the cell types and subtypes in the human and mouse brain. Currently, STAB2 includes 71 and 61 different cell subtypes defined in the human and mouse brain, respectively. It covers 63 subregions and 15 developmental stages of human brain, and 38 subregions and 30 developmental stages of mouse brain, generating a comprehensive atlas for exploring spatiotemporal transcriptomic dynamics in the mammalian brain. We also augmented web interfaces for querying and visualizing the gene expression in specific cell types. STAB2 is freely available at https://mai.fudan.edu.cn/stab2.

Autism genes converge on microtubule biology and RNA-binding proteins during excitatory neurogenesis

Recent studies have identified over one hundred high-confidence (hc) autism spectrum disorder (ASD) genes. Systems biological and functional analyses on smaller subsets of these genes have consistently implicated excitatory neurogenesis. However, the extent to which the broader set of hcASD genes are involved in this process has not been explored systematically nor have the biological pathways underlying this convergence been identified. Here, we leveraged CROP-Seq to repress 87 hcASD genes in a human in vitro model of cortical neurogenesis. We identified 17 hcASD genes whose repression significantly alters developmental trajectory and results in a common cellular state characterized by disruptions in proliferation, differentiation, cell cycle, microtubule biology, and RNA-binding proteins (RBPs). We also characterized over 3,000 differentially expressed genes, 286 of which had expression profiles correlated with changes in developmental trajectory. Overall, we uncovered transcriptional disruptions downstream of hcASD gene perturbations, correlated these disruptions with distinct differentiation phenotypes, and reinforced neurogenesis, microtubule biology, and RBPs as convergent points of disruption in ASD.

Genetics of human brain development

Brain development in humans is achieved through precise spatiotemporal genetic control, the mechanisms of which remain largely elusive. Recently, integration of technological advances in human stem cell-based modelling with genome editing has emerged as a powerful platform to establish causative links between genotypes and phenotypes directly in the human system. Here, we review our current knowledge of complex genetic regulation of each key step of human brain development through the lens of evolutionary specialization and neurodevelopmental disorders and highlight the use of human stem cell-derived 2D cultures and 3D brain organoids to investigate human-enriched features and disease mechanisms. We also discuss opportunities and challenges of integrating new technologies to reveal the genetic architecture of human brain development and disorders.

Astrocytic β-catenin signaling via TCF7L2 regulates synapse development and social behavior

The Wnt/β-catenin pathway contains multiple high-confidence risk genes that are linked to neurodevelopmental disorders, including autism spectrum disorder. However, its ubiquitous roles across brain cell types and developmental stages have made it challenging to define its impact on neural circuit development and behavior. Here, we show that TCF7L2, which is a key transcriptional effector of the Wnt/β-catenin pathway, plays a cell-autonomous role in postnatal astrocyte maturation and impacts adult social behavior. TCF7L2 was the dominant Wnt effector that was expressed in both mouse and human astrocytes, with a peak during astrocyte maturation. The conditional knockout of Tcf7l2 in postnatal astrocytes led to an enlargement of astrocytes with defective tiling and gap junction coupling. These mice also exhibited an increase in the number of cortical excitatory and inhibitory synapses and a marked increase in social interaction by adulthood. These data reveal an astrocytic role for developmental Wnt/β-catenin signaling in restricting excitatory synapse numbers and regulating adult social behavior.

A human embryonic stem cell-based model reveals the cell of origin of FOXR2-activated CNS neuroblastoma

Background

FOXR2-activated central nervous system (CNS) neuroblastoma (CNS NB-FOXR2) is a recently identified subtype of brain tumor characterized by the elevated expression of the transcription factor FOXR2 mainly due to genomic rearrangements. However, the precise pathogenic mechanisms, including the cell type of origin, remain elusive.

Methods

A gene expression analysis of patient tumors was performed to identify putative cell types of origin. Based on this prediction, a new human embryonic stem cell-based model was developed to validate the origin and to examine the molecular and cellular mechanisms underlying the formation of CNS NB-FOXR2.

Results

Our data showed that CNS NB-FOXR2 tumors express a high level of lineage marker genes associated with the medial ganglionic eminence (MGE), a transient structure located in the developing ventral forebrain. Our model confirmed the cell-type-specific effect of FOXR2 on the proliferation and in vivo tumorigenicity. Additionally, we found that FOXR2 overexpression activated the MEK/ERK signaling pathway through a suppression of the endogenous RAS inhibitor DIRAS3. The MEK inhibitor trametinib suppressed the proliferation of FOXR2-expressing MGE progenitors more than nonexpressing cells.

Conclusions

Our study collectively demonstrates that MGE progenitors are the cell of origin of CNS NB-FOXR2 and that FOXR2 activates the MEK/ERK signaling pathway, providing a potential therapeutic target.

Concurrent impact of de novo mutations on cranial and cortical development in nonsyndromic craniosynostosis

OBJECTIVE

Nonsyndromic craniosynostosis (nsCS), characterized by premature cranial suture fusion, is considered a primary skull disorder in which impact on neurodevelopment, if present, results from the mechanical hindrance of brain growth. Despite surgical repair of the cranial defect, neurocognitive deficits persist in nearly half of affected children. Therefore, the authors performed a functional genomics analysis of nsCS to determine when, where, and in what cell types nsCS-associated genes converge during development.

METHODS

The authors integrated whole-exome sequencing data from 291 nsCS proband-parent trios with 29,803 single-cell transcriptomes of the prenatal and postnatal neurocranial complex to inform when, where, and in what cell types nsCS-mutated genes might exert their pathophysiological effects.

RESULTS

The authors found that nsCS-mutated genes converged in cranial osteoprogenitors and pial fibroblasts and their transcriptional networks that regulate both skull ossification and cerebral neurogenesis. Nonsyndromic CS-mutated genes also converged in inhibitory neurons and gene coexpression modules that overlapped with autism and other developmental disorders. Ligand-receptor cell-cell communication analysis uncovered crosstalk between suture osteoblasts and neurons via the nsCS-associated BMP, FGF, and noncanonical WNT signaling pathways.

CONCLUSIONS

These data implicate a concurrent impact of nsCS-associated de novo mutations on cranial morphogenesis and cortical development via cell- and non-cell-autonomous mechanisms in a developmental nexus of fetal osteoblasts, pial fibroblasts, and neurons. These results suggest that neurodevelopmental outcomes in nsCS patients may be driven more by mutational status than surgical technique.

Long and Short-Term Effect of mTOR Regulation on Cerebral Organoid Growth and Differentiations

BACKGROUND

The mammalian target of rapamycin (mTOR) signaling is critical for the maintenance and differentiation of neurogenesis, and conceivably for many other brain developmental processes. However, in vivo studies of mTOR functions in the brain are often hampered due to the essential role of the associated signaling in brain development.

METHODS

We monitored the long- and short-term effects of mTOR signaling regulation on cerebral organoids growth, differentiation and function using an mTOR inhibitor (everolimus) and an mTOR activator (MHY1485).

RESULTS

Short-term treatment with MHY1485 induced faster organoid growth and differentiation, while long-term treatment induced the maturation of cerebral organoids.

CONCLUSION

These data suggest that the optimal activity of mTOR is crucial in maintaining normal brain development, and its role is not confined to the early neurogenic phase of brain development.

Shaping the brain: The emergence of cortical structure and folding

The cerebral cortex-the brain's covering and largest region-has increased in size and complexity in humans and supports higher cognitive functions such as language and abstract thinking. There is a growing understanding of the human cerebral cortex, including the diversity and number of cell types that it contains, as well as of the developmental mechanisms that shape cortical structure and organization. In this review, we discuss recent progress in our understanding of molecular and cellular processes, as well as mechanical forces, that regulate the folding of the cerebral cortex. Advances in human genetics, coupled with experimental modeling in gyrencephalic species, have provided insights into the central role of cortical progenitors in the gyrification and evolutionary expansion of the cerebral cortex. These studies are essential for understanding the emergence of structural and functional organization during cortical development and the pathogenesis of neurodevelopmental disorders associated with cortical malformations.

Pathway centric analysis for single-cell RNA-seq and spatial transcriptomics data with GSDensity

Advances in single-cell technology have enabled molecular dissection of heterogeneous biospecimens at unprecedented scales and resolutions. Cluster-centric approaches are widely applied in analyzing single-cell data, however they have limited power in dissecting and interpreting highly heterogenous, dynamically evolving data. Here, we present GSDensity, a graph-modeling approach that allows users to obtain pathway-centric interpretation and dissection of single-cell and spatial transcriptomics (ST) data without performing clustering. Using pathway gene sets, we show that GSDensity can accurately detect biologically distinct cells and reveal novel cell-pathway associations ignored by existing methods. Moreover, GSDensity, combined with trajectory analysis can identify curated pathways that are active at various stages of mouse brain development. Finally, GSDensity can identify spatially relevant pathways in mouse brains and human tumors including those following high-order organizational patterns in the ST data. Particularly, we create a pan-cancer ST map revealing spatially relevant and recurrently active pathways across six different tumor types.

FEZ1 participates in human embryonic brain development by modulating neuronal progenitor subpopulation specification and migrations

Mutations in the human fasciculation and elongation protein zeta 1 (FEZ1) gene are found in schizophrenia and Jacobsen syndrome patients. Here, using human cerebral organoids (hCOs), we show that FEZ1 expression is turned on early during brain development and is detectable in both neuroprogenitor subtypes and immature neurons. FEZ1 deletion disrupts expression of neuronal and synaptic development genes. Using single-cell RNA sequencing, we detected abnormal expansion of homeodomain-only protein homeobox (HOPX)- outer radial glia (oRG), concurrent with a reduction of HOPX+ oRG, in FEZ1-null hCOs. HOPX- oRGs show higher cell mobility as compared to HOPX+ oRGs. Ectopic localization of neuroprogenitors to the outer layer is seen in FEZ1-null hCOs. Anomalous encroachment of TBR2+ intermediate progenitors into CTIP2+ deep layer neurons further indicated abnormalities in cortical layer formation these hCOs. Collectively, our findings highlight the involvement of FEZ1 in early cortical brain development and how it contributes to neurodevelopmental disorders.

Human brain organoid model of maternal immune activation identifies radial glia cells as selectively vulnerable

Maternal immune activation (MIA) during critical windows of gestation is correlated with long-term neurodevelopmental deficits in the offspring, including increased risk for autism spectrum disorder (ASD) in humans. Interleukin 6 (IL-6) derived from the gestational parent is one of the major molecular mediators by which MIA alters the developing brain. In this study, we establish a human three-dimensional (3D) in vitro model of MIA by treating induced pluripotent stem cell-derived dorsal forebrain organoids with a constitutively active form of IL-6, Hyper-IL-6. We validate our model by showing that dorsal forebrain organoids express the molecular machinery necessary for responding to Hyper-IL-6 and activate STAT signaling upon Hyper-IL-6 treatment. RNA sequencing analysis reveals the upregulation of major histocompatibility complex class I (MHCI) genes in response to Hyper-IL-6 exposure, which have been implicated with ASD. We find a small increase in the proportion of radial glia cells after Hyper-IL-6 treatment through immunohistochemistry and single-cell RNA-sequencing. We further show that radial glia cells are the cell type with the highest number of differentially expressed genes, and Hyper-IL-6 treatment leads to the downregulation of genes related to protein translation in line with a mouse model of MIA. Additionally, we identify differentially expressed genes not found in mouse models of MIA, which might drive species-specific responses to MIA. Finally, we show abnormal cortical layering as a long-term consequence of Hyper-IL-6 treatment. In summary, we establish a human 3D model of MIA, which can be used to study the cellular and molecular mechanisms underlying the increased risk for developing disorders such as ASD.

Non-synaptic function of the autism spectrum disorder-associated gene SYNGAP1 in cortical neurogenesis

Genes involved in synaptic function are enriched among those with autism spectrum disorder (ASD)-associated rare genetic variants. Dysregulated cortical neurogenesis has been implicated as a convergent mechanism in ASD pathophysiology, yet it remains unknown how 'synaptic' ASD risk genes contribute to these phenotypes, which arise before synaptogenesis. Here, we show that the synaptic Ras GTPase-activating (RASGAP) protein 1 (SYNGAP1, a top ASD risk gene) is expressed within the apical domain of human radial glia cells (hRGCs). In a human cortical organoid model of SYNGAP1 haploinsufficiency, we find dysregulated cytoskeletal dynamics that impair the scaffolding and division plane of hRGCs, resulting in disrupted lamination and accelerated maturation of cortical projection neurons. Additionally, we confirmed an imbalance in the ratio of progenitors to neurons in a mouse model of Syngap1 haploinsufficiency. Thus, SYNGAP1-related brain disorders may arise through non-synaptic mechanisms, highlighting the need to study genes associated with neurodevelopmental disorders (NDDs) in diverse human cell types and developmental stages.

GTF2I dosage regulates neuronal differentiation and social behavior in 7q11.23 neurodevelopmental disorders

Copy number variations at 7q11.23 cause neurodevelopmental disorders with shared and opposite manifestations. Deletion causes Williams-Beuren syndrome featuring hypersociability, while duplication causes 7q11.23 microduplication syndrome (7Dup), frequently exhibiting autism spectrum disorder (ASD). Converging evidence indicates GTF2I as key mediator of the cognitive-behavioral phenotypes, yet its role in cortical development and behavioral hallmarks remains largely unknown. We integrated proteomic and transcriptomic profiling of patient-derived cortical organoids, including longitudinally at single-cell resolution, to dissect 7q11.23 dosage-dependent and GTF2I-specific disease mechanisms. We observed dosage-dependent impaired dynamics of neural progenitor proliferation, transcriptional imbalances, and highly specific alterations in neuronal output, leading to precocious excitatory neuron production in 7Dup, which was rescued by restoring physiological GTF2I levels. Transgenic mice with Gtf2i duplication recapitulated progenitor proliferation and neuronal differentiation defects alongside ASD-like behaviors. Consistently, inhibition of lysine demethylase 1 (LSD1), a GTF2I effector, was sufficient to rescue ASD-like phenotypes in transgenic mice, establishing GTF2I-LSD1 axis as a molecular pathway amenable to therapeutic intervention in ASD.

Truncated radial glia as a common precursor in the late corticogenesis of gyrencephalic mammals

The diversity of neural stem cells is a hallmark of the cerebral cortex development in gyrencephalic mammals, such as Primates and Carnivora. Among them, ferrets are a good model for mechanistic studies. However, information on their neural progenitor cells (NPC), termed radial glia (RG), is limited. Here, we surveyed the temporal series of single-cell transcriptomes of progenitors regarding ferret corticogenesis and found a conserved diversity and temporal trajectory between human and ferret NPC, despite the large timescale difference. We found truncated RG (tRG) in ferret cortical development, a progenitor subtype previously described in humans. The combination of in silico and in vivo analyses identified that tRG differentiate into both ependymal and astrogenic cells. Via transcriptomic comparison, we predict that this is also the case in humans. Our findings suggest that tRG plays a role in the formation of adult ventricles, thereby providing the architectural bases for brain expansion.

Unraveling Neuronal Identities Using SIMS: A Deep Learning Label Transfer Tool for Single-Cell RNA Sequencing Analysis

Large single-cell RNA datasets have contributed to unprecedented biological insight. Often, these take the form of cell atlases and serve as a reference for automating cell labeling of newly sequenced samples. Yet, classification algorithms have lacked the capacity to accurately annotate cells, particularly in complex datasets. Here we present SIMS (Scalable, Interpretable Machine Learning for Single-Cell), an end-to-end data-efficient machine learning pipeline for discrete classification of single-cell data that can be applied to new datasets with minimal coding. We benchmarked SIMS against common single-cell label transfer tools and demonstrated that it performs as well or better than state of the art algorithms. We then use SIMS to classify cells in one of the most complex tissues: the brain. We show that SIMS classifies cells of the adult cerebral cortex and hippocampus at a remarkably high accuracy. This accuracy is maintained in trans-sample label transfers of the adult human cerebral cortex. We then apply SIMS to classify cells in the developing brain and demonstrate a high level of accuracy at predicting neuronal subtypes, even in periods of fate refinement, shedding light on genetic changes affecting specific cell types across development. Finally, we apply SIMS to single cell datasets of cortical organoids to predict cell identities and unveil genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. When cell types are obscured by stress signals, label transfer from primary tissue improves the accuracy of cortical organoid annotations, serving as a reliable ground truth. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

Identification of glioblastoma stem cell-associated lncRNAs using single-cell RNA sequencing datasets

Glioblastoma multiforme (GBM) is an aggressive, heterogeneous brain tumor in which glioblastoma stem cells (GSCs) are known culprits of therapy resistance. Long non-coding RNAs (lncRNAs) have been shown to play a critical role in both cancer and normal biology. A few studies have suggested that aberrant expression of lncRNAs is associated with GSCs. However, a comprehensive single-cell analysis of the GSC-associated lncRNA transcriptome has not been carried out. Here, we analyzed recently published single-cell RNA sequencing datasets of adult GBM tumors, GBM organoids, GSC-enriched GBM tumors, and developing human brain samples to identify lncRNAs highly expressed in GSCs. We further revealed that the GSC-specific lncRNAs GIHCG and LINC01563 promote proliferation, migration, and stemness in the GSC population. Together, this study identified a panel of uncharacterized GSC-enriched lncRNAs and set the stage for future in-depth studies to examine their role in GBM pathology and their potential as biomarkers and/or therapeutic targets in GBM.

Single-cell genomics reveals region-specific developmental trajectories underlying neuronal diversity in the human hypothalamus

The development and diversity of neuronal subtypes in the human hypothalamus has been insufficiently characterized. To address this, we integrated transcriptomic data from 241,096 cells (126,840 newly generated) in the prenatal and adult human hypothalamus to reveal a temporal trajectory from proliferative stem cell populations to mature hypothalamic cell types. Iterative clustering of the adult neurons identified 108 robust transcriptionally distinct neuronal subtypes representing 10 hypothalamic nuclei. Pseudotime trajectories provided insights into the genes driving formation of these nuclei. Comparisons to single-cell transcriptomic data from the mouse hypothalamus suggested extensive conservation of neuronal subtypes despite certain differences in species-enriched gene expression. The uniqueness of hypothalamic neuronal lineages was examined developmentally by comparing excitatory lineages present in cortex and inhibitory lineages in ganglionic eminence, revealing both distinct and shared drivers of neuronal maturation across the human forebrain. These results provide a comprehensive transcriptomic view of human hypothalamus development through gestation and adulthood at cellular resolution.

Single-cell analysis reveals specific neuronal transition during mouse corticogenesis

Background: Currently, the mechanism(s) underlying corticogenesis is still under characterization. Methods: We curated the most comprehensive single-cell RNA-seq (scRNA-seq) datasets from mouse and human fetal cortexes for data analysis and confirmed the findings with co-immunostaining experiments. Results: By analyzing the developmental trajectories with scRNA-seq datasets in mice, we identified a specific developmental sub-path contributed by a cell-population expressing both deep- and upper-layer neurons (DLNs and ULNs) specific markers, which occurred on E13.5 but was absent in adults. In this cell-population, the percentages of cells expressing DLN and ULN markers decreased and increased, respectively, during the development suggesting direct neuronal transition (namely D-T-U). Whilst genes significantly highly/uniquely expressed in D-T-U cell population were significantly enriched in PTN/MDK signaling pathways related to cell migration. Both findings were further confirmed by co-immunostaining with DLNs, ULNs and D-T-U specific markers across different timepoints. Furthermore, six genes (co-expressed with D-T-U specific markers in mice) showing a potential opposite temporal expression between human and mouse during fetal cortical development were associated with neuronal migration and cognitive functions. In adult prefrontal cortexes (PFC), D-T-U specific genes were expressed in neurons from different layers between humans and mice. Conclusion: Our study characterizes a specific cell population D-T-U showing direct DLNs to ULNs neuronal transition and migration during fetal cortical development in mice. It is potentially associated with the difference of cortical development in humans and mice.

δ-Catenin controls astrocyte morphogenesis via layer-specific astrocyte-neuron cadherin interactions

Astrocytes control the formation of specific synaptic circuits via cell adhesion and secreted molecules. Astrocyte synaptogenic functions are dependent on the establishment of their complex morphology. However, it is unknown if distinct neuronal cues differentially regulate astrocyte morphogenesis. δ-Catenin was previously thought to be a neuron-specific protein that regulates dendrite morphology. We found δ-catenin is also highly expressed by astrocytes and required both in astrocytes and neurons for astrocyte morphogenesis. δ-Catenin is hypothesized to mediate transcellular interactions through the cadherin family of cell adhesion proteins. We used structural modeling and biochemical analyses to reveal that δ-catenin interacts with the N-cadherin juxtamembrane domain to promote N-cadherin surface expression. An autism-linked δ-catenin point mutation impaired N-cadherin cell surface expression and reduced astrocyte complexity. In the developing mouse cortex, only lower-layer cortical neurons express N-cadherin. Remarkably, when we silenced astrocytic N-cadherin throughout the cortex, only lower-layer astrocyte morphology was disrupted. These findings show that δ-catenin controls astrocyte-neuron cadherin interactions that regulate layer-specific astrocyte morphogenesis.

Pluripotent stem cell-derived neural progenitor cells can be used to model effects of IL-6 on human neurodevelopment

Maternal immune activation (MIA) increases the risks for neurodevelopmental disorders in offspring through inflammatory cytokines, including interleukin-6 (IL-6). We therefore aimed to establish a human two-dimensional (2D) in vitro neural model to investigate the effects of IL-6 exposure on neurodevelopment. IL-6 signal transduction requires two receptors: interleukin-6 signal transducer (IL6ST) and interleukin-6 receptor (IL6R). Prenatally, neural cells lack IL6R, and hence cannot elicit cis IL-6 signaling, but IL6R can be provided by microglia in trans. We demonstrate here that an immortalized human neural progenitor cell (NPC) line, ReNCell CX, expresses IL6ST and elicits both cis and trans IL-6 signaling, limiting its use as a model of MIA. In contrast, induced pluripotent stem cell (iPSC)-derived NPCs only activate the IL-6 cascade in trans. Activation of the trans IL-6 cascade did not result in increased proliferation of iPSC-derived NPCs or ReNCell CX, as has been demonstrated in animal models. iPSC-derived NPCs upregulated NR2F1 expression in response to IL-6 signaling in line with analogous experiments in organoids. Thus, iPSC-derived NPCs can be used to model gene expression changes in response to MIA in 2D cultures.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

Preservation of co-expression defines the primary tissue fidelity of human neural organoids

Human neural organoid models offer an exciting opportunity for studying often inaccessible human-specific brain development; however, it remains unclear how precisely organoids recapitulate fetal/primary tissue biology. Here, we characterize field-wide replicability and biological fidelity through a meta-analysis of single-cell RNA-sequencing data for first and second trimester human primary brain (2.95 million cells, 51 datasets) and neural organoids (1.63 million cells, 130 datasets). We quantify the degree to which primary tissue cell-type marker expression and co-expression are recapitulated in organoids across 12 different protocol types. By quantifying gene-level preservation of primary tissue co-expression, we show neural organoids lie on a spectrum ranging from virtually no signal to co-expression near indistinguishable from primary tissue data, demonstrating high fidelity is within the scope of current methods. Additionally, we show neural organoids preserve the cell-type specific co-expression of developing rather than adult cells, confirming organoids are an appropriate model for primary tissue development. Overall, quantifying the preservation of primary tissue co-expression is a powerful tool for uncovering unifying axes of variation across heterogeneous neural organoid experiments.

NeST: nested hierarchical structure identification in spatial transcriptomic data

Spatial gene expression in tissue is characterized by regions in which particular genes are enriched or depleted. Frequently, these regions contain nested inside them subregions with distinct expression patterns. Segmentation methods in spatial transcriptomic (ST) data extract disjoint regions maximizing similarity over the greatest number of genes, typically on a particular spatial scale, thus lacking the ability to find region-within-region structure. We present NeST, which extracts spatial structure through coexpression hotspots-regions exhibiting localized spatial coexpression of some set of genes. Coexpression hotspots identify structure on any spatial scale, over any possible subset of genes, and are highly explainable. NeST also performs spatial analysis of cell-cell interactions via ligand-receptor, identifying active areas de novo without restriction of cell type or other groupings, in both two and three dimensions. Through application on ST datasets of varying type and resolution, we demonstrate the ability of NeST to reveal a new level of biological structure.

Multi-omic profiling of the developing human cerebral cortex at the single-cell level

The cellular complexity of the human brain is established via dynamic changes in gene expression throughout development that is mediated, in part, by the spatiotemporal activity of cis-regulatory elements (CREs). We simultaneously profiled gene expression and chromatin accessibility in 45,549 cortical nuclei across six broad developmental time points from fetus to adult. We identified cell type-specific domains in which chromatin accessibility is highly correlated with gene expression. Differentiation pseudotime trajectory analysis indicates that chromatin accessibility at CREs precedes transcription and that dynamic changes in chromatin structure play a critical role in neuronal lineage commitment. In addition, we mapped cell type-specific and temporally specific genetic loci implicated in neuropsychiatric traits, including schizophrenia and bipolar disorder. Together, our results describe the complex regulation of cell composition at critical stages in lineage determination and shed light on the impact of spatiotemporal alterations in gene expression on neuropsychiatric disease.

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.

Molecular programs of regional specification and neural stem cell fate progression in macaque telencephalon

During early telencephalic development, intricate processes of regional patterning and neural stem cell (NSC) fate specification take place. However, our understanding of these processes in primates, including both conserved and species-specific features, remains limited. Here, we profiled 761,529 single-cell transcriptomes from multiple regions of the prenatal macaque telencephalon. We deciphered the molecular programs of the early organizing centers and their cross-talk with NSCs, revealing primate-biased galanin-like peptide (GALP) signaling in the anteroventral telencephalon. Regional transcriptomic variations were observed along the frontotemporal axis during early stages of neocortical NSC progression and in neurons and astrocytes. Additionally, we found that genes associated with neuropsychiatric disorders and brain cancer risk might play critical roles in the early telencephalic organizers and during NSC progression.

Single-cell analysis of prenatal and postnatal human cortical development

We analyzed >700,000 single-nucleus RNA sequencing profiles from 106 donors during prenatal and postnatal developmental stages and identified lineage-specific programs that underlie the development of specific subtypes of excitatory cortical neurons, interneurons, glial cell types, and brain vasculature. By leveraging single-nucleus chromatin accessibility data, we delineated enhancer gene regulatory networks and transcription factors that control commitment of specific cortical lineages. By intersecting our results with genetic risk factors for human brain diseases, we identified the cortical cell types and lineages most vulnerable to genetic insults of different brain disorders, especially autism. We find that lineage-specific gene expression programs up-regulated in female cells are especially enriched for the genetic risk factors of autism. Our study captures the molecular progression of cortical lineages across human development.

Systematic investigation of allelic regulatory activity of schizophrenia-associated common variants

Genome-wide association studies (GWASs) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium makes it challenging to discern causal variants. We performed a massively parallel reporter assay (MPRA) on 5,173 fine-mapped schizophrenia GWAS variants in primary human neural progenitors and identified 439 variants with allelic regulatory effects (MPRA-positive variants). Transcription factor binding had modest predictive power, while fine-map posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit expressive quantitative trait loci signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. To predict the combinatorial effect of MPRA-positive variants on gene regulation, we propose an accessibility-by-contact model that combines MPRA-measured allelic activity with neuronal chromatin architecture.

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 or transcript-isoform diversity during human brain development has not been systematically investigated. Here, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 unique isoforms, of which 72.6% are novel (unannotated in Gencode-v33), 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 re-prioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

One-Sentence Summary

A cell-specific atlas of gene isoform expression helps shape our understanding of brain development and disease.

Structured Abstract

INTRODUCTION: The development of the human brain is regulated by precise molecular and genetic mechanisms driving spatio-temporal and cell-type-specific transcript expression programs. Alternative splicing, a major mechanism increasing transcript diversity, is highly prevalent in the human brain, influences many aspects of brain development, and has strong links to neuropsychiatric disorders. Despite this, the cell-type-specific transcript-isoform diversity of the developing human brain has not been systematically investigated.RATIONALE: Understanding splicing patterns and isoform diversity across the developing neocortex has translational relevance and can elucidate genetic risk mechanisms in neurodevelopmental disorders. However, short-read sequencing, the prevalent technology for transcriptome profiling, is not well suited to capturing alternative splicing and isoform diversity. To address this, we employed third-generation long-read sequencing, which enables capture and sequencing of complete individual RNA molecules, to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution.RESULTS: We profiled microdissected GZ and CP regions of post-conception week (PCW) 15-17 human neocortex in bulk and at single-cell resolution across six subjects using high-fidelity long-read sequencing (PacBio IsoSeq). We identified 214,516 unique isoforms, of which 72.6% were novel (unannotated in Gencode), and >7,000 novel exons, expanding the proteome by 92,422 putative proteoforms. We uncovered thousands of isoform switches during cortical neurogenesis predicted to impact RNA regulatory domains or protein structure and implicating previously uncharacterized RNA-binding proteins in cellular identity and neuropsychiatric disease. At the single-cell level, early-stage excitatory neurons exhibited the greatest isoform diversity, and isoform-centric single-cell clustering led to the identification of previously uncharacterized cell states. We systematically assessed the contribution of transcriptomic features, and localized cell and spatio-temporal transcript expression signatures across neuropsychiatric disorders, revealing predominant enrichments in dynamic isoform expression and utilization patterns and that the number and complexity of isoforms per gene is strongly predictive of disease. Leveraging this resource, we re-prioritized thousands of rare de novo risk variants associated with autism spectrum disorders (ASD), intellectual disability (ID), and neurodevelopmental disorders (NDDs), more broadly, to potentially more severe consequences and revealed a larger proportion of cryptic splice variants with the expanded transcriptome annotation provided in this study.CONCLUSION: Our study offers a comprehensive landscape of isoform diversity in the human neocortex during development. This extensive cataloging of novel isoforms and splicing events sheds light on the underlying mechanisms of neurodevelopmental disorders and presents an opportunity to explore rare genetic variants linked to these conditions. The implications of our findings extend beyond fundamental neuroscience, as they provide crucial insights into the molecular basis of developmental brain disorders and pave the way for targeted therapeutic interventions. To facilitate exploration of this dataset we developed an online portal ( https://sciso.gandallab.org/ ).

Tissue morphology influences the temporal program of human brain organoid development

Progression through fate decisions determines cellular composition and tissue architecture, but how that same architecture may impact cell fate is less clear. We took advantage of organoids as a tractable model to interrogate this interaction of form and fate. Screening methodological variations revealed that common protocol adjustments impacted various aspects of morphology, from macrostructure to tissue architecture. We examined the impact of morphological perturbations on cell fate through integrated single nuclear RNA sequencing (snRNA-seq) and spatial transcriptomics. Regardless of the specific protocol, organoids with more complex morphology better mimicked in vivo human fetal brain development. Organoids with perturbed tissue architecture displayed aberrant temporal progression, with cells being intermingled in both space and time. Finally, encapsulation to impart a simplified morphology led to disrupted tissue cytoarchitecture and a similar abnormal maturational timing. These data demonstrate that cells of the developing brain require proper spatial coordinates to undergo correct temporal progression.

Aberrant pace of cortical neuron development in brain organoids from patients with 22q11.2 deletion syndrome and schizophrenia

Adults and children afflicted with the 22q11.2 deletion syndrome (22q11.2DS) exhibit cognitive, social, and emotional impairments, and are at significantly heightened risk for schizophrenia (SCZ). The impact of this deletion on early human brain development, however, has remained unclear. Here we harness organoid models of the developing human cerebral cortex, cultivated from subjects with 22q11.2DS and SCZ, as well as unaffected control samples, to identify cell-type-specific developmental abnormalities arising from this genomic lesion. Leveraging single-cell RNA-sequencing in conjunction with experimental validation, we find that the loss of genes within the 22q11.2 locus leads to a delayed development of cortical neurons. This compromised development was reflected in an elevated proportion of actively proliferating neural progenitor cells, coupled with a decreased fraction of more mature neurons. Furthermore, we identify perturbed molecular imprints linked to neuronal maturation, observe the presence of sparser neurites, and note a blunted amplitude in glutamate-induced Ca2+ transients. The aberrant transcription program underlying impaired development contains molecular signatures significantly enriched in neuropsychiatric genetic liability. MicroRNA profiling and target gene investigation suggest that microRNA dysregulation may drive perturbations of genes governing the pace at which maturation unfolds. Using protein-protein interaction network analysis we define complementary effects stemming from additional genes residing within the deleted locus. Our study uncovers reproducible neurodevelopmental and molecular alterations due to 22q11.2 deletions. These findings have the potential to facilitate disease modeling and promote the pursuit of therapeutic interventions.

LIF signaling regulates outer radial glial to interneuron fate during human cortical development

Radial glial (RG) development is essential for cerebral cortex growth and organization. In humans, the outer radial glia (oRG) subtype is expanded and gives rise to diverse neurons and glia. However, the mechanisms regulating oRG differentiation are unclear. oRG cells express leukemia-inhibitory factor (LIF) receptors during neurogenesis, and consistent with a role in stem cell self-renewal, LIF perturbation impacts oRG proliferation in cortical tissue and organoids. Surprisingly, LIF treatment also increases the production of inhibitory interneurons (INs) in cortical cultures. Comparative transcriptomic analysis identifies that the enhanced IN population resembles INs produced in the caudal ganglionic eminence. To evaluate whether INs could arise from oRGs, we isolated primary oRG cells and cultured them with LIF. We observed the production of INs from oRG cells and an increase in IN abundance following LIF treatment. Our observations suggest that LIF signaling regulates the capacity of oRG cells to generate INs.

Neurodevelopmental deficits and cell-type-specific transcriptomic perturbations in a mouse model of HNRNPU haploinsufficiency

Heterozygous de novo loss-of-function mutations in the gene expression regulator HNRNPU cause an early-onset developmental and epileptic encephalopathy. To gain insight into pathological mechanisms and lay the potential groundwork for developing targeted therapies, we characterized the neurophysiologic and cell-type-specific transcriptomic consequences of a mouse model of HNRNPU haploinsufficiency. Heterozygous mutants demonstrated global developmental delay, impaired ultrasonic vocalizations, cognitive dysfunction and increased seizure susceptibility, thus modeling aspects of the human disease. Single-cell RNA-sequencing of hippocampal and neocortical cells revealed widespread, yet modest, dysregulation of gene expression across mutant neuronal subtypes. We observed an increased burden of differentially-expressed genes in mutant excitatory neurons of the subiculum-a region of the hippocampus implicated in temporal lobe epilepsy. Evaluation of transcriptomic signature reversal as a therapeutic strategy highlights the potential importance of generating cell-type-specific signatures. Overall, this work provides insight into HNRNPU-mediated disease mechanisms and provides a framework for using single-cell RNA-sequencing to study transcriptional regulators implicated in disease.

A Meta-Atlas of the Developing Human Cortex Identifies Modules Driving Cell Subtype Specification

Human brain development requires the generation of hundreds of diverse cell types, a process targeted by recent single-cell transcriptomic profiling efforts. Through a meta-analysis of seven of these published datasets, we have generated 225 meta-modules - gene co-expression networks that can describe mechanisms underlying cortical development. Several meta-modules have potential roles in both establishing and refining cortical cell type identities, and we validated their spatiotemporal expression in primary human cortical tissues. These include meta-module 20, associated with FEZF2+ deep layer neurons. Half of meta-module 20 genes are putative FEZF2 targets, including TSHZ3, a transcription factor associated with neurodevelopmental disorders. Human cortical organoid experiments validated that both factors are necessary for deep layer neuron specification. Importantly, subtle manipulations of these factors drive slight changes in meta-module activity that cascade into strong differences in cell fate - demonstrating how of our meta-atlas can engender further mechanistic analyses of cortical fate specification.

Human microglia maturation is underpinned by specific gene regulatory networks

Microglia phenotypes are highly regulated by the brain environment, but the transcriptional networks that specify the maturation of human microglia are poorly understood. Here, we characterized stage-specific transcriptomes and epigenetic landscapes of fetal and postnatal human microglia and acquired corresponding data in induced pluripotent stem cell (iPSC)-derived microglia, in cerebral organoids, and following engraftment into humanized mice. Parallel development of computational approaches that considered transcription factor (TF) co-occurrence and enhancer activity allowed prediction of shared and state-specific gene regulatory networks associated with fetal and postnatal microglia. Additionally, many features of the human fetal-to-postnatal transition were recapitulated in a time-dependent manner following the engraftment of iPSC cells into humanized mice. These data and accompanying computational approaches will facilitate further efforts to elucidate mechanisms by which human microglia acquire stage- and disease-specific phenotypes.