Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex

Sequential emergence and contraction of epithelial subtypes in the prenatal human choroid plexus revealed by a stem cell model

Despite the major roles of choroid plexus epithelial cells (CPECs) in brain homeostasis and repair, their developmental lineage and diversity remain undefined. In simplified differentiations from human pluripotent stem cells, derived CPECs (dCPECs) display canonical properties and dynamic motile multiciliated phenotypes that interact with Aβ uptake. Single dCPEC transcriptomes over time correlate well with human organoid and fetal CPECs, while pseudotemporal and cell cycle analyses highlight the direct CPEC origin from neuroepithelial cells. In addition, time series analyses define metabolic (type 1) and ciliogenic dCPECs (type 2) at early timepoints, followed by type 1 diversification into anabolic-secretory (type 1a) and catabolic-absorptive subtypes (type 1b) as type 2 cells contract. These temporal patterns are then confirmed in independent derivations and mapped to prenatal stages using human tissues. In addition to defining the prenatal lineage of human CPECs, these findings suggest dynamic models of ChP support for the developing human brain.

Cellular hierarchies of embryonal tumors with multilayered rosettes are shaped by oncogenic microRNAs and receptor-ligand interactions

Embryonal tumor with multilayered rosettes (ETMR) is a pediatric brain tumor with dismal prognosis. Characteristic alterations of the chromosome 19 microRNA cluster (C19MC) are observed in most ETMR; however, the ramifications of C19MC activation and the complex cellular architecture of ETMR remain understudied. Here we analyze 11 ETMR samples from patients using single-cell transcriptomics and multiplexed spatial imaging. We reveal a spatially distinct cellular hierarchy that spans highly proliferative neural stem-like cells and more differentiated neuron-like cells. C19MC is predominantly expressed in stem-like cells and controls a transcriptional network governing stemness and lineage commitment, as resolved by genome-wide analysis of microRNA-mRNA binding. Systematic analysis of receptor-ligand interactions between malignant cell types reveals fibroblast growth factor receptor and Notch signaling as oncogenic pathways that can be successfully targeted in preclinical models and in one patient with ETMR. Our study provides fundamental insights into ETMR pathobiology and a powerful rationale for more effective targeted therapies.

CellWalker2: Multi-omic discovery using hierarchical cell type relationships

Tissues are composed of cells with a wide range of similarities to each other, yet existing methods for single-cell genomics treat cell types as discrete labels. To address this gap, we developed CellWalker2, a graph diffusion-based model for the annotation and mapping of multi-modal data. With our open-source software package, hierarchically related cell types can be probabilistically matched across contexts and used to annotate cells, genomic regions, or gene sets. Additional features include estimating statistical significance and enabling gene expression and chromatin accessibility to be jointly modeled. Through simulation studies, we show that CellWalker2 performs better than existing methods in cell-type annotation and mapping. We then use multi-omics data from the brain and immune system to demonstrate CellWalker2's ability to assign high-resolution cell-type labels to regulatory elements and TFs and to quantify both conserved and divergent cell-type relationships between species.

"Teaching the biology of neurodiversity for social change"

The scientific understanding of autism spectrum disorder (ASD) is historically rooted in the diagnosis of children with divergent behavior and socialization within Western societies. Increased activism and visibility of ASD individuals has helped popularize the neurodiversity movement, which proposes that autism falls within the natural spectrum of human behavior and is not always pathological. The increased interest in neurodiversity has coincided with rapid advances in the field of psychiatric genetics, particularly in understanding the underlying genetic causes of ASD. Identified genes linked to ASD have highlighted mid-fetal cortical development as a nexus when many of these genes are co-expressed. New discoveries in the genetics of ASD present a unique opportunity in teaching developmental biology and especially brain development to engage students with the biology underlying ASD and neurodiversity in general. Informed and inclusive language, engagement with the neurodiversity community, and awareness of the complexity of the issues are good practices in teaching the biology of neurodiversity.

Spatial transcriptomics reveals human cortical layer and area specification

The human cerebral cortex is composed of six layers and dozens of areas that are molecularly and structurally distinct1-4. Although single-cell transcriptomic studies have advanced the molecular characterization of human cortical development, a substantial gap exists owing to the loss of spatial context during cell dissociation5-8. Here we used multiplexed error-robust fluorescence in situ hybridization (MERFISH)9, augmented with deep-learning-based nucleus segmentation, to examine the molecular, cellular and cytoarchitectural development of the human fetal cortex with spatially resolved single-cell resolution. Our extensive spatial atlas, encompassing more than 18 million single cells, spans eight cortical areas across seven developmental time points. We uncovered the early establishment of the six-layer structure, identifiable by the laminar distribution of excitatory neuron subtypes, 3 months before the emergence of cytoarchitectural layers. Notably, we discovered two distinct modes of cortical areal specification during mid-gestation: (1) a continuous, gradual transition observed across most cortical areas along the anterior-posterior axis and (2) a discrete, abrupt boundary specifically identified between the primary (V1) and secondary (V2) visual cortices as early as gestational week 20. This sharp binary transition in V1-V2 neuronal subtypes challenges the notion that mid-gestation cortical arealization involves only gradient-like transitions6,10. Furthermore, integrating single-nucleus RNA sequencing with MERFISH revealed an early upregulation of synaptogenesis in V1-specific layer 4 neurons. Collectively, our findings underscore the crucial role of spatial relationships in determining the molecular specification of cortical layers and areas. This study establishes a spatially resolved single-cell analysis paradigm and paves the way for the construction of a comprehensive developmental atlas of the human brain.

Multimodal analyses reveal genes driving electrophysiological maturation of neurons in the primate prefrontal cortex

The prefrontal cortex (PFC) is critical for myriad high-cognitive functions and is associated with several neuropsychiatric disorders. Here, using Patch-seq and single-nucleus multiomic analyses, we identified genes and regulatory networks governing the maturation of distinct neuronal populations in the PFC of rhesus macaque. We discovered that specific electrophysiological properties exhibited distinct maturational kinetics and identified key genes underlying these properties. We unveiled that RAPGEF4 is important for the maturation of resting membrane potential and inward sodium current in both macaque and human. We demonstrated that knockdown of CHD8, a high-confidence autism risk gene, in human and macaque organotypic slices led to impaired maturation, via downregulation of key genes, including RAPGEF4. Restoring the expression of RAPGEF4 rescued the proper electrophysiological maturation of CHD8-deficient neurons. Our study revealed regulators of neuronal maturation during a critical period of PFC development in primates and implicated such regulators in molecular processes underlying autism.

A human-specific enhancer fine-tunes radial glia potency and corticogenesis

Humans have evolved an extraordinarily expanded and complex cerebral cortex associated with developmental and gene regulatory modifications1-3. Human accelerated regions (HARs) are highly conserved DNA sequences with human-specific nucleotide substitutions. Although there are thousands of annotated HARs, their functional contribution to species-specific cortical development remains largely unknown4,5. HARE5 is a HAR transcriptional enhancer of the WNT signalling receptor Frizzled8 that is active during brain development6. Here, using genome-edited mouse (Mus musculus, Mm) and primate models, we demonstrated that human (Homo sapiens, Hs) HARE5 fine-tunes cortical development and connectivity by controlling the proliferative and neurogenic capacities of neural progenitor cells. Hs-HARE5 knock-in mice have significantly enlarged neocortices, containing more excitatory neurons. By measuring neural dynamics in vivo, we showed that these anatomical features result in increased functional independence between cortical regions. We assessed underlying developmental mechanisms using fixed and live imaging, lineage analysis and single-cell RNA sequencing. We discovered that Hs-HARE5 modifies radial glial cell behaviour, with increased self-renewal at early developmental stages, followed by expanded neurogenic potential. Using genome-edited human and chimpanzee (Pan troglodytes, Pt) neural progenitor cells and cortical organoids, we showed that four human-specific variants of Hs-HARE5 drive increased enhancer activity that promotes progenitor proliferation. Finally, we showed that Hs-HARE5 increased progenitor proliferation by amplifying canonical WNT signalling. These findings illustrate how small changes in regulatory DNA can directly affect critical signalling pathways to modulate brain development. Our study uncovered new functions of HARs as key regulatory elements crucial for the expansion and complexity of the human cerebral cortex.

Patient-derived models of UBA5-associated encephalopathy identify defects in neurodevelopment and highlight potential therapeutic avenues

UBA5 encodes for the E1 enzyme of the UFMylation cascade, which plays an essential role in endoplasmic reticulum (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 from two patients with compound heterozygous variants in UBA5. Both shared the same missense variant, which encodes a hypomorphic allele (p.A371T), along with a nonsense variant (p.G267* or p.A123fs*4). Single-cell RNA sequencing of 100-day organoids identified defects in GABAergic interneuron development. We demonstrated aberrant neuronal firing and reduction in size of patient-derived organoids. Mechanistically, we showed that ER homeostasis is perturbed along with an exacerbated unfolded protein response pathway in engineered U87-MG cells and patient-derived organoids expressing UBA5 pathogenic variants. We also assessed two potential therapeutic modalities that augmented UBA5 protein abundance to rescue aberrant molecular and cellular phenotypes. We assessed SINEUP, a long noncoding RNA that augments translation efficiency, and CRISPRa, a modified CRISPR-Cas9 approach to augment transcription efficiency to increase UBA5 protein production. Our study provides a humanized model that allows further investigations of UBA5 variants in the brain and highlights promising approaches to alleviate cellular aberrations for this rare, developmental disorder.

Deciphering the longitudinal trajectories of glioblastoma ecosystems by integrative single-cell genomics

The evolution of isocitrate dehydrogenase (IDH)-wildtype glioblastoma (GBM) after standard-of-care therapy remains poorly understood. Here we analyzed matched primary and recurrent GBMs from 59 patients using single-nucleus RNA sequencing and bulk DNA sequencing, assessing the longitudinal evolution of the GBM ecosystem across layers of cellular and molecular heterogeneity. The most consistent change was a lower malignant cell fraction at recurrence and a reciprocal increase in glial and neuronal cell types in the tumor microenvironment (TME). The predominant malignant cell state differed between most matched pairs, but no states were exclusive or highly enriched in either time point, nor was there a consistent longitudinal trajectory across the cohort. Nevertheless, specific trajectories were enriched in subsets of patients. Changes in malignant state abundances mirrored changes in TME composition and baseline profiles, reflecting the co-evolution of the GBM ecosystem. Our study provides a blueprint of GBM's diverse longitudinal trajectories and highlights the treatment and TME modifiers that shape them.

The multilayered transcriptional architecture of glioblastoma ecosystems

In isocitrate dehydrogenase wildtype glioblastoma (GBM), cellular heterogeneity across and within tumors may drive therapeutic resistance. Here we analyzed 121 primary and recurrent GBM samples from 59 patients using single-nucleus RNA sequencing and bulk tumor DNA sequencing to characterize GBM transcriptional heterogeneity. First, GBMs can be classified by their broad cellular composition, encompassing malignant and nonmalignant cell types. Second, in each cell type we describe the diversity of cellular states and their pathway activation, particularly an expanded set of malignant cell states, including glial progenitor cell-like, neuronal-like and cilia-like. Third, the remaining variation between GBMs highlights three baseline gene expression programs. These three layers of heterogeneity are interrelated and partially associated with specific genetic aberrations, thereby defining three stereotypic GBM ecosystems. This work provides an unparalleled view of the multilayered transcriptional architecture of GBM. How this architecture evolves during disease progression is addressed in the companion manuscript by Spitzer et al.

DEAD/DEAH-box RNA helicases shape the risk of neurodevelopmental disorders

The DEAD/DEAH-box family of RNA helicases (RHs) is among the most abundant and conserved in eukaryotes. These proteins catalyze the remodeling of RNAs to regulate their splicing, stability, localization, and translation. Rare genetic variants in DEAD/DEAH-box proteins have recently emerged as being associated with neurodevelopmental disorders (NDDs). Analyses in cellular and animal models have uncovered fundamental roles for these proteins during brain development. We discuss the genetic and functional evidence that implicates DEAD/DEAH-box proteins in brain development and NDDs, with a focus on how structural insights from paralogous genes can be leveraged to advance our understanding of the pathogenic mechanisms at play.

Integrated analysis of molecular atlases unveils modules driving developmental cell subtype specification in the human cortex

Human brain development requires generating diverse cell types, a process explored by single-cell transcriptomics. Through parallel meta-analyses of the human cortex in development (seven datasets) and adulthood (16 datasets), we generated over 500 gene co-expression networks that can describe mechanisms of cortical development, centering on peak stages of neurogenesis. These meta-modules show dynamic cell subtype specificities throughout cortical development, with several developmental meta-modules displaying spatiotemporal expression patterns that allude to potential roles in cell fate specification. We validated the expression of these modules in primary human cortical tissues. These include meta-module 20, a module elevated in FEZF2+ deep layer neurons that includes TSHZ3, a transcription factor associated with neurodevelopmental disorders. Human cortical chimeroid experiments validated that both FEZF2 and TSHZ3 are required to drive module 20 activity and deep layer neuron specification but through distinct modalities. These studies demonstrate how meta-atlases can engender further mechanistic analyses of cortical fate specification.

Fate plasticity of interneuron specification

Neuronal subtype generation in the mammalian central nervous system is governed by competing genetic programs. The medial ganglionic eminence (MGE) produces two major cortical interneuron (IN) populations, somatostatin (Sst) and parvalbumin (Pvalb), which develop on different timelines. The extent to which external signals influence these identities remains unclear. Pvalb-positive INs are crucial for cortical circuit regulation but challenging to model in vitro. We grafted mouse MGE progenitors into diverse 2D and 3D co-culture systems, including mouse and human cortical, MGE, and thalamic models. Strikingly, only 3D human corticogenesis models promoted efficient, non-autonomous Pvalb differentiation, characterized by upregulation of Pvalb maturation markers, downregulation of Sst-specific markers, and the formation of perineuronal nets. Additionally, lineage-traced postmitotic Sst-positive INs upregulated Pvalb when grafted onto human cortical models. These findings reveal unexpected fate plasticity in MGE-derived INs, suggesting that their identities can be dynamically shaped by the environment.

CeiTEA: Adaptive Hierarchy of Single Cells with Topological Entropy

Advances in single-cell RNA sequencing (scRNA-seq) enable detailed analysis of cellular heterogeneity, but existing clustering methods often fail to capture the complex hierarchical structures of cell types and subtypes. CeiTEA is introduced, a novel algorithm for adaptive hierarchical clustering based on topological entropy (TE), designed to address this challenge. CeiTEA constructs a multi-nary partition tree that optimally represents relationships and diversity among cell types by minimizing TE. This method combines a bottom-up strategy for hierarchy construction with a top-down strategy for local diversification, facilitating the identification of smaller hierarchical structures within subtrees. CeiTEA is evaluated on both simulated and real-world scRNA-seq datasets, demonstrating superior clustering performance compared to state-of-the-art tools like Louvain, Leiden, K-means, and SEAT. In simulated multi-layer datasets, CeiTEA demonstrated superior performance in retrieving hierarchies with a lower average clustering information distance of 0.15, compared to 0.39 from SEAT and 0.67 from traditional hierarchical clustering methods. On real datasets, the CeiTEA hierarchy reflects the developmental potency of various cell populations, validated by gene ontology enrichment, cell-cell interaction, and pseudo-time analysis. These findings highlight CeiTEA's potential as a powerful tool for understanding complex relationships in single-cell data, with applications in tumor heterogeneity and tissue specification.

Individualized patient tumor organoids faithfully preserve human brain tumor ecosystems and predict patient response to therapy

Tumor organoids are important tools for cancer research, but current models have drawbacks that limit their applications for predicting response to therapy. Here, we developed a fast, efficient, and complex culture system (IPTO, individualized patient tumor organoid) that accurately recapitulates the cellular and molecular pathology of human brain tumors. Patient-derived tumor explants were cultured in induced pluripotent stem cell (iPSC)-derived cerebral organoids, thus enabling culture of a wide range of human tumors in the central nervous system (CNS), including adult, pediatric, and metastatic brain cancers. Histopathological, genomic, epigenomic, and single-cell RNA sequencing (scRNA-seq) analyses demonstrated that the IPTO model recapitulates cellular heterogeneity and molecular features of original tumors. Crucially, we showed that the IPTO model predicts patient-specific drug responses, including resistance mechanisms, in a prospective patient cohort. Collectively, the IPTO model represents a major breakthrough in preclinical modeling of human cancers, which provides a path toward personalized cancer therapy.

Genomic and Transcriptomic Signatures of SETD1A Disruption in Human Excitatory Neuron Development and Psychiatric Disease Risk

Genetic disruption of SETD1A markedly increases the risk for schizophrenia. To elucidate the underlying mechanisms, we generated isogenic organoid models of the developing human cerebral cortex harboring a SETD1A loss-of-function schizophrenia risk mutation. Employing chromatin profiling combined with RNA sequencing, we identified high-confidence SETD1A target genes, analyzed the impact of the mutation on SETD1A binding and transcriptional regulation and validated key findings with orthogonal approaches. Disruption of SETD1A function disturbs the finely tuned temporal gene expression in the excitatory neuron lineage, yielding an aberrant transcriptional program that compromises key regulatory and metabolic pathways essential for neurodevelopmental transitions. Although overall SETD1A binding remains unchanged in mutant neurons, we identified localized alterations in SETD1A binding that correlate with shifts in H3K4me3 levels and gene expression. These changes are enriched at enhancer regions, suggesting that enhancer-regulated genes are especially vulnerable to SETD1A reduction. Notably, target genes with enhancer-bound SETD1A are primarily linked to neuronal functions while those with promoter-bound SETD1A are enriched for basic cellular functions. By mapping the SETD1A binding landscape in excitatory neurons of the human fetal frontal cortex and integrating multimodal neuroimaging and genetic datasets, we demonstrate that the genomic context of SETD1A binding differentially correlates with macroscale brain organization and establish a link between SETD1A-bound enhancers, schizophrenia-associated brain alterations and genetic susceptibility. Our study advances our understanding of the role of SETD1A binding patterns in schizophrenia pathogenesis, offering insights that may guide future therapeutic strategies.

Development of the rodent prefrontal cortex: circuit formation, plasticity, and impacts of early life stress

The prefrontal cortex (PFC), located at the anterior region of the cerebral cortex, is a multimodal association cortex essential for higher-order brain functions, including decision-making, attentional control, memory processing, and regulation of social behavior. Structural, circuit-level, and functional abnormalities in the PFC are often associated with neurodevelopmental disorders. Here, we review recent findings on the postnatal development of the PFC, with a particular emphasis on rodent studies, to elucidate how its structural and circuit properties are established during critical developmental windows and how these processes influence adult behaviors. Recent evidence also highlights the lasting effects of early life stress on the PFC structure, connectivity, and function. We explore potential mechanisms underlying these stress-induced alterations, with a focus on epigenetic regulation and its implications for PFC maturation and neurodevelopmental disorders. By integrating these insights, this review provides an overview of the developmental processes shaping the PFC and their implications for brain health and disease.

Metabolic Atlas of Early Human Cortex Identifies Regulators of Cell Fate Transitions

Characterization of cell type emergence during human cortical development, which enables unique human cognition, has focused primarily on anatomical and transcriptional characterizations. Metabolic processes in the human brain that allow for rapid expansion, but contribute to vulnerability to neurodevelopmental disorders, remain largely unexplored. We performed a variety of metabolic assays in primary tissue and stem cell derived cortical organoids and observed dynamic changes in core metabolic functions, including an unexpected increase in glycolysis during late neurogenesis. By depleting glucose levels in cortical organoids, we increased outer radial glia, astrocytes, and inhibitory neurons. We found the pentose phosphate pathway (PPP) was impacted in these experiments and leveraged pharmacological and genetic manipulations to recapitulate these radial glia cell fate changes. These data identify a new role for the PPP in modulating radial glia cell fate specification and generate a resource for future exploration of additional metabolic pathways in human cortical development.

Massively parallel reporter assay investigates shared genetic variants of eight psychiatric disorders

A meta-genome-wide association study across eight psychiatric disorders has highlighted the genetic architecture of pleiotropy in major psychiatric disorders. However, mechanisms underlying pleiotropic effects of the associated variants remain to be explored. We conducted a massively parallel reporter assay to decode the regulatory logic of variants with pleiotropic and disorder-specific effects. Pleiotropic variants differ from disorder-specific variants by exhibiting chromatin accessibility that extends across diverse cell types in the neuronal lineage and by altering motifs for transcription factors with higher connectivity in protein-protein interaction networks. We mapped pleiotropic and disorder-specific variants to putative target genes using functional genomics approaches and CRISPR perturbation. In vivo CRISPR perturbation of a pleiotropic and a disorder-specific gene suggests that pleiotropy may involve the regulation of genes expressed broadly across neuronal cell types and with higher network connectivity.

CRISPRi-based screens in iAssembloids to elucidate neuron-glia interactions

The complexity of the human brain makes it challenging to understand the molecular mechanisms underlying brain function. Genome-wide association studies have uncovered variants associated with neurological phenotypes. Single-cell transcriptomics have provided descriptions of changes brain cells undergo during disease. However, these approaches do not establish molecular mechanism. To facilitate the scalable interrogation of causal molecular mechanisms in brain cell types, we developed a 3D co-culture system of induced pluripotent stem cell (iPSC)-derived neurons and glia, termed iAssembloids. Using iAssembloids, we ask how glial and neuronal cells interact to control neuronal death and survival. Our CRISPRi-based screens identified that GSK3β inhibits the protective NRF2-mediated oxidative stress response elicited by high neuronal activity. We then investigate the role of APOE-ε4, a risk variant for Alzheimer's disease, on neuronal survival. We find that APOE-ε4-expressing astrocytes may promote neuronal hyperactivity as compared with APOE-ε3-expressing astrocytes. This platform allows for the unbiased identification of mechanisms of neuron-glia cell interactions.

The Principle of Cortical Development and Evolution

Human's robust cognitive abilities, including creativity and language, are made possible, at least in large part, by evolutionary changes made to the cerebral cortex. This paper reviews the biology and evolution of mammalian cortical radial glial cells (primary neural stem cells) and introduces the concept that a genetically step wise process, based on a core molecular pathway already in use, is the evolutionary process that has molded cortical neurogenesis. The core mechanism, which has been identified in our recent studies, is the extracellular signal-regulated kinase (ERK)-bone morphogenic protein 7 (BMP7)-GLI3 repressor form (GLI3R)-sonic hedgehog (SHH) positive feedback loop. Additionally, I propose that the molecular basis for cortical evolutionary dwarfism, exemplified by the lissencephalic mouse which originated from a larger gyrencephalic ancestor, is an increase in SHH signaling in radial glia, that antagonizes ERK-BMP7 signaling. Finally, I propose that: (1) SHH signaling is not a key regulator of primate cortical expansion and folding; (2) human cortical radial glial cells do not generate neocortical interneurons; (3) human-specific genes may not be essential for most cortical expansion. I hope this review assists colleagues in the field, guiding research to address gaps in our understanding of cortical development and evolution.

PTEN mutations impair CSF dynamics and cortical networks by dysregulating periventricular neural progenitors

Enlargement of the cerebrospinal fluid (CSF)-filled brain ventricles (ventriculomegaly) is a defining feature of congenital hydrocephalus (CH) and an under-recognized concomitant of autism. Here, we show that de novo mutations in the autism risk gene PTEN are among the most frequent monogenic causes of CH and primary ventriculomegaly. Mouse Pten-mutant ventriculomegaly results from aqueductal stenosis due to hyperproliferation of periventricular Nkx2.1+ neural progenitor cells (NPCs) and increased CSF production from hyperplastic choroid plexus. Pten-mutant ventriculomegalic cortices exhibit network dysfunction from increased activity of Nkx2.1+ NPC-derived inhibitory interneurons. Raptor deletion or postnatal everolimus treatment corrects ventriculomegaly, rescues cortical deficits and increases survival by antagonizing mTORC1-dependent Nkx2.1+ NPC pathology. Thus, PTEN mutations concurrently alter CSF dynamics and cortical networks by dysregulating Nkx2.1+ NPCs. These results implicate a nonsurgical treatment for CH, demonstrate a genetic association of ventriculomegaly and ASD, and help explain neurodevelopmental phenotypes refractory to CSF shunting in select individuals with CH.

Single-Cell Analysis of Sex and Gender Differences in the Human Brain During Development and Disease

Sex and gender (SG) differences in the human brain are of interest to society and science as numerous processes are impacted by them, including brain development, behavior, and diseases. By collecting publicly available single-cell data from the in-utero to elderly age in healthy, Alzheimer's disease and multiple sclerosis samples, we identified and characterized SG-biased genes in ten brain cell types across 9 age and disease groups. Sex and gender differences in the transcriptome were present throughout the lifespan and across all cell types. Although there was limited overlap among SG-biased genes across different age and disease groups, we observed significant functional overlap. Female-biased genes are consistently enriched for brain-related processes, while male-biased genes are enriched for metabolic pathways. Additionally, mitochondrial genes showed a consistent female bias across cell types. We also found that androgen response elements (not estrogen) were significantly enriched in both male- and female-biased genes, and thymosin hormone targets being consistently enriched only in male-biased genes. We systematically characterised SG differences in brain development and brain-related disorders at a single-cell level, by analysing a total of publicly available 419,885 single nuclei from 161 human brain samples (72 females, 89 males). The significant enrichment of androgen (not estrogen) response elements in both male- and female-biased genes suggests that androgens are important regulators likely establishing these SG differences. Finally, we provide full characterization of SG-biased genes at different thresholds for the scientific community as a web resource.

Cell Type-Agnostic Transcriptomic Signatures Enable Uniform Comparisons of Neurodevelopment

Single-cell transcriptomics has revolutionized our understanding of neurodevelopmental cell identities, yet, predicting a cell type's developmental state from its transcriptome remains a challenge. We perform a meta-analysis of developing human brain datasets comprising over 2.8 million cells, identifying both tissue-level and cell-autonomous predictors of developmental age. While tissue composition predicts age within individual studies, it fails to generalize, whereas specific cell type proportions reliably track developmental time across datasets. Training regularized regression models to infer cell-autonomous maturation, we find that a cell type-agnostic model achieves the highest accuracy (error = 2.6 weeks), robustly capturing developmental dynamics across diverse cell types and datasets. This model generalizes to human neural organoids, accurately predicting normal developmental trajectories (R = 0.91) and disease-induced shifts in vitro. Furthermore, it extends to the developing mouse brain, revealing an accelerated developmental tempo relative to humans. Our work provides a unified framework for comparing neurodevelopment across contexts, model systems, and species.

Multimodal transcriptomics reveal neurogenic aging trajectories and age-related regional inflammation in the dentate gyrus

The mammalian dentate gyrus (DG) is involved in certain forms of learning and memory, and DG dysfunction has been implicated in age-related diseases. Although neurogenic potential is maintained throughout life in the DG as neural stem cells (NSCs) continue to generate new neurons, neurogenesis decreases with advancing age, with implications for age-related cognitive decline and disease. In this study, we used single-cell RNA sequencing to characterize transcriptomic signatures of neurogenic cells and their surrounding DG niche, identifying molecular changes associated with neurogenic aging from the activation of quiescent NSCs to the maturation of fate-committed progeny. By integrating spatial transcriptomics data, we identified the regional invasion of inflammatory cells into the hippocampus with age and show here that early-onset neuroinflammation decreases neurogenic activity. Our data reveal the lifelong molecular dynamics of NSCs and their surrounding neurogenic DG niche with age and provide a powerful resource to understand age-related molecular alterations in the aging hippocampus.

Classifying cell cycle states and a quiescent-like G0 state using single-cell transcriptomics

Single-cell transcriptomics has unveiled a vast landscape of cellular heterogeneity in which the cell cycle is a significant component. We trained a high-resolution cell cycle classifier (ccAFv2) using single cell RNA-seq (scRNA-seq) characterized human neural stem cells. The ccAFv2 classifies six cell cycle states (G1, Late G1, S, S/G2, G2/M, and M/Early G1) and a quiescent-like G0 state (Neural G0), and it incorporates a tunable parameter to filter out less certain classifications. The ccAFv2 classifier performed better than or equivalent to other state-of-the-art methods even while classifying more cell cycle states, including G0. We demonstrate that the ccAFv2 classifier effectively generalizes the S, S/G2, G2/M, and M/Early G1 states across cell types derived from all three germ layers. While the G0, G1, and Late G1 states perform well in neuroepithelial cell types, their accuracy is lower in other cell types. However, misclassifications are confined to the G0, G1, and Late G1 states. We showcased the versatility of ccAFv2 by successfully applying it to classify cells, nuclei, and spatial transcriptomics data in humans and mice, using various normalization methods and gene identifiers. We provide methods to regress the cell cycle expression patterns out of single cell or nuclei data to uncover underlying biological signals. The classifier can be used either as an R package integrated with Seurat or a PyPI package integrated with SCANPY. We proved that ccAFv2 has enhanced accuracy, flexibility, and adaptability across various experimental conditions, establishing ccAFv2 as a powerful tool for dissecting complex biological systems, unraveling cellular heterogeneity, and deciphering the molecular mechanisms by which proliferation and quiescence affect cellular processes.

A subpopulation of cortical neurons altered by mutations in the autism risk gene DDX3X

Cell fate decisions during cortical development sculpt the identity of long-range connections that subserve complex behaviors. These decisions are largely dictated by mutually exclusive transcription factors, including CTIP2/Bcl11b for subcerebral projection neurons and BRN1/Pou3f3 for intra-telencephalic projection neurons. We have recently reported that the balance of cortical CTIP2-expressing neurons is altered in a mouse model of DDX3X syndrome, a female-biased neurodevelopmental disorder associated with intellectual disability, autism spectrum disorder, and significant motor challenges. Here, we studied the developmental dynamics of a subpopulation of cortical neurons co-expressing CTIP2 and BRN1. We found that CTIP2+BRN1+ neurons are born during early phases of neurogenesis like other CTIP2+ neurons, peak in expression during perinatal life, and persist in adult brains. We also found that CTIP2+BRN1+ neurons are excessive in number in prenatal and mature cortical motor areas of Ddx3x mutant mice, translating into altered laminar distribution of subcerebral projection neurons extending axons to the brainstem. These findings underscore the critical role of molecular specification during cortical development in health and disease.

Human neural rosettes secrete bioactive extracellular vesicles enriched in neuronal and glial cellular components

Extracellular vesicles (EVs) play a critical role in the development of neural cells in the central nervous system (CNS). Human neural rosettes (hNRs) are radial cell structures that assemble from induced pluripotent stem cells (hiPSCs) and recapitulate some stages of neural tube morphogenesis. Here we show that hiPSCs and hNRs secrete EVs (hiPSC-EVs and hNR-EVs) with distinctive protein cargoes. Remarkably, hNR-EVs carry neuronal and glial cellular components involved in human CNS development. Importantly, hNR-EVs stimulate stem cells to change their cellular morphology and promote neurite growth in human and murine neurons with a significant dysregulation of SOX2 levels. This transcription factor modulates both neural differentiation and pluripotency. Interestingly, these effects were inhibited by antibodies against an unexpected neuroglial cargo of hNR-EVs: the major proteolipid protein (PLP). These findings show that hNRs secrete bioactive EVs containing neural components and might contribute as trophic factors during human neurodevelopment.

ChromBPNet: bias factorized, base-resolution deep learning models of chromatin accessibility reveal cis-regulatory sequence syntax, transcription factor footprints and regulatory variants

Despite extensive mapping of cis-regulatory elements (cREs) across cellular contexts with chromatin accessibility assays, the sequence syntax and genetic variants that regulate transcription factor (TF) binding and chromatin accessibility at context-specific cREs remain elusive. We introduce ChromBPNet, a deep learning DNA sequence model of base-resolution accessibility profiles that detects, learns and deconvolves assay-specific enzyme biases from regulatory sequence determinants of accessibility, enabling robust discovery of compact TF motif lexicons, cooperative motif syntax and precision footprints across assays and sequencing depths. Extensive benchmarks show that ChromBPNet, despite its lightweight design, is competitive with much larger contemporary models at predicting variant effects on chromatin accessibility, pioneer TF binding and reporter activity across assays, cell contexts and ancestry, while providing interpretation of disrupted regulatory syntax. ChromBPNet also helps prioritize and interpret regulatory variants that influence complex traits and rare diseases, thereby providing a powerful lens to decode regulatory DNA and genetic variation.

MAPbrain: a multi-omics atlas of the primate brain

The brain is the central hub of the entire nervous system. Its development is a lifelong process guided by a genetic blueprint. Understanding how genes influence brain development is critical for deciphering the formation of human cognitive functions and the underlying mechanisms of neurological disorders. Recent advances in multi-omics techniques have now made it possible to explore these aspects comprehensively. However, integrating and analyzing extensive multi-omics data presents significant challenges. Here, we introduced MAPbrain (http://bigdata.ibp.ac.cn/mapBRAIN/), a multi-omics atlas of the primate brain. This repository integrates and normalizes both our own lab's published data and publicly available multi-omics data, encompassing 21 million brain cells from 38 key brain regions and 436 sub-regions across embryonic and adult stages, with 164 time points in humans and non-human primates. MAPbrain offers a unique, robust, and interactive platform that includes transcriptomics, epigenomics, and spatial transcriptomics data, facilitating a comprehensive exploration of brain development. The platform enables the exploration of cell type- and time point-specific markers, gene expression comparison between brain regions and species, joint analyses across transcriptome and epigenome, and navigation of cell types across species, brain regions, and development stages. Additionally, MAPbrain provides an online integration module for users to navigate and analyze their own data within the platform.

Depth-corrected multi-factor dissection of chromatin accessibility for scATAC-seq data with PACS

Single cell ATAC-seq (scATAC-seq) experimental designs have become increasingly complex, with multiple factors that might affect chromatin accessibility, including genotype, cell type, tissue of origin, sample location, batch, etc., whose compound effects are difficult to test by existing methods. In addition, current scATAC-seq data present statistical difficulties due to their sparsity and variations in individual sequence capture. To address these problems, we present a zero-adjusted statistical model, Probability model of Accessible Chromatin of Single cells (PACS), that allows complex hypothesis testing of accessibility-modulating factors while accounting for sparse and incomplete data. For differential accessibility analysis, PACS controls the false positive rate and achieves a 17% to 122% higher power on average than existing tools. We demonstrate the effectiveness of PACS through several analysis tasks, including supervised cell type annotation, compound hypothesis testing, batch effect correction, and spatiotemporal modeling. We apply PACS to datasets from various tissues and show its ability to reveal previously undiscovered insights in scATAC-seq data.

Deficiency of FABP7 Triggers Premature Neural Differentiation in Idiopathic Normocephalic Autism Organoids

Autism spectrum disorder (ASD), which is caused by heterogeneous genetic and environmental factors, is characterized by diverse clinical phenotypes linked to distinct pathological mechanisms. ASD individuals with a shared clinical phenotype might contribute to revealing the molecular mechanism underlying ASD progression. Here, it is generated induced pluripotent stem cell (iPSC)-derived cerebral organoids from normocephalic individuals with ASD in a prospective birth cohort with a shared clinical diagnosis. Multiple cell lines and time series scRNA-seq combined with a histomorphological analysis revealed premature neural differentiation of neural stem cells (NSCs) and decreased expression of Fatty acid binding protein 7 (FABP7) in ASD organoids. It is subsequently revealed alterations in the phosphorylation levels of Mitogen-Activated Protein Kinase Kinase 1/2 (MEK1/2), which are downstream of FABP7, and the regulation of the FABP7/MEK pathway reversed improper neural differentiation in the ASD organoids. Moreover, both Fabp7-knockdown and MEK2-overexpressing mice exhibited repetitive stereotyped behaviors and social defects relevant to autism. This study reveals the role of the FABP7/MEK pathway in abnormal NSC differentiation in normocephalic individuals with ASD, which might provide a promising therapeutic target for ASD treatment.

Understanding ubiquitination in neurodevelopment by integrating insights across space and time

Ubiquitination regulates a myriad of eukaryotic signaling cascades by modifying substrate proteins, thereby determining their functions and fates. In this perspective, we discuss current challenges in investigating the ubiquitin system in the developing brain. We foster the concept that ubiquitination pathways are spatiotemporally regulated and tightly intertwined with molecular and cellular transitions during neurogenesis and neural circuit assembly. Focusing on the neurologically highly relevant class of homologous to E6AP C-terminus (HECT) ubiquitin ligases, we propose cross-disciplinary translational approaches bridging state-of-the-art cell biology, proteomics, biochemistry, structural biology and neuroscience to dissect ubiquitination in neurodevelopment and its specific perturbations in brain diseases. We highlight that a comprehensive understanding of ubiquitin signaling in the brain may reveal new horizons in basic neuroscience and clinical applications.

Protocol for comparative gene expression data analysis between brains and organoids using a cloud-based web app

Here, we present a protocol for using Brain and Organoid Manifold Alignment (BOMA), a cloud-based web app for comparative gene expression data analysis between brains and organoids. We describe steps for performing a global alignment of developmental gene expression data from both brains and organoids. We then detail procedures for investigating both shared and distinctive developmental pathways across brains and organoids by refining alignment locally using manifold learning. This protocol is applicable for working with single-cell and bulk RNA sequencing data. For complete details on the use and execution of this protocol, please refer to He et al.1.

Multi-omics approaches reveal that diffuse midline gliomas present altered DNA replication and are susceptible to replication stress therapy

BACKGROUND

The fatal diffuse midline gliomas (DMG) are characterized by an undruggable H3K27M mutation in H3.1 or H3.3. K27M impairs normal development by stalling differentiation. The identification of targetable pathways remains very poorly explored. Toward this goal, we undertake a multi-omics approach to evaluate replication timing profiles, transcriptomics, and cell cycle features in DMG cells from both H3.1K27M and H3.3K27M subgroups and perform a comparative, integrative data analysis with healthy brain tissue.

RESULTS

DMG cells present differential replication timing in each subgroup, which, in turn, correlates with significant differential gene expression. Differentially expressed genes in S phase are involved in various pathways related to DNA replication. We detect increased expression of DNA replication genes earlier in the cell cycle in DMG cell lines compared to normal brain cells. Furthermore, the distance between origins of replication in DMG cells is smaller than in normal brain cells and their fork speed is slower, a read-out of replication stress. Consistent with these findings, DMG tumors present high replication stress signatures in comparison to normal brain cells. Finally, DMG cells are specifically sensitive to replication stress therapy.

CONCLUSIONS

This whole genome multi-omics approach provides insights into the cell cycle regulation of DMG via the H3K27M mutations and establishes a pharmacologic vulnerability in DNA replication, which resolves a potentially novel therapeutic strategy for this non-curable disease.

Reliability of high-quantity human brain organoids for modeling microcephaly, glioma invasion and drug screening

Brain organoids offer unprecedented insights into brain development and disease modeling and hold promise for drug screening. Significant hindrances, however, are morphological and cellular heterogeneity, inter-organoid size differences, cellular stress, and poor reproducibility. Here, we describe a method that reproducibly generates thousands of organoids across multiple hiPSC lines. These High Quantity brain organoids (Hi-Q brain organoids) exhibit reproducible cytoarchitecture, cell diversity, and functionality, are free from ectopically active cellular stress pathways, and allow cryopreservation and re-culturing. Patient-derived Hi-Q brain organoids recapitulate distinct forms of developmental defects: primary microcephaly due to a mutation in CDK5RAP2 and progeria-associated defects of Cockayne syndrome. Hi-Q brain organoids displayed a reproducible invasion pattern for a given patient-derived glioma cell line. This enabled a medium-throughput drug screen to identify Selumetinib and Fulvestrant, as inhibitors of glioma invasion in vivo. Thus, the Hi-Q approach can easily be adapted to reliably harness brain organoids' application for personalized neurogenetic disease modeling and drug discovery.

Initial regional cytoarchitectonic differences in dorsal and orbitobasal human developing frontal cortex revealed by spatial transcriptomics

Early development of the human fetal cerebral cortex involves a set of precisely coordinated molecular processes that remains rather underexplored. Previous studies indicate that the laminar identity and the molecular specification of cortical neurons driven by genetic programming, as well as associated histogenetic events begin during early fetal development. Our recent study discovered unique regional cytoarchitectonic features in the developing human frontal lobe, including migratory waves of postmitotic neurons in the dorsal frontal cortex and the "double plate" feature in orbitobasal cortex (Kopić et al. in Cells 12:231, 2023). Notably, neurons of these two cytoarchitectonic features typically express deep projection neuron (DPN) markers (TBR1, TLE4, SOX5). This paper aims to conduct an in-depth investigation of these cytoarchitectonic features at the transcriptomic level, whilst preserving spatial information. Here, we employed NanoString GeoMx™ Digital Spatial Profiler (DSP) technology to examine gene expression differences in the transient cortical compartments of the dorsal and ventral regions of the developing frontal lobe, focusing specifically on 15 post-conceptional weeks (PCW), that is a critical period for subplate formation. We identified multiple differentially expressed genes between the transient cellular compartments of the dorsal and orbitobasal regions of the developing human frontal cortex. These new findings additionally confirm that regional patterning and specification of the prospective higher-order association prefrontal cortex emerges early in fetal development, contributing to the highly organized cortical architecture of the human brain.

An ILK/STAT3 pathway controls glioblastoma stem cell plasticity

Glioblastoma (GBM) is driven by malignant neural stem-like cells that display extensive heterogeneity and phenotypic plasticity, which drive tumor progression and therapeutic resistance. Here, we show that the extracellular matrix-cell adhesion protein integrin-linked kinase (ILK) stimulates phenotypic plasticity and mesenchymal-like, invasive behavior in a murine GBM stem cell model. ILK is required for the interconversion of GBM stem cells between malignancy-associated glial-like states, and its loss produces cells that are unresponsive to multiple cell state transition cues. We further show that an ILK/STAT3 signaling pathway controls the plasticity that enables transition of GBM stem cells to an astrocyte-like state in vitro and in vivo. Finally, we find that ILK expression correlates with expression of STAT3-regulated proteins and protein signatures describing astrocyte-like and mesenchymal states in patient tumors. This work identifies ILK as a pivotal regulator of multiple malignancy-associated GBM phenotypes, including phenotypic plasticity and mesenchymal state.

Applying single-cell and single-nucleus genomics to studies of cellular heterogeneity and cell fate transitions in the nervous system

Single-cell and single-nucleus genomic approaches can provide unbiased and multimodal insights. Here, we discuss what constitutes a molecular cell atlas and how to leverage single-cell omics data to generate hypotheses and gain insights into cell transitions in development and disease of the nervous system. We share points of reflection on what to consider during study design and implementation as well as limitations and pitfalls.

Human-specific genetic hallmarks in neocortical development: focus on neural progenitors

The evolutionary expansion of the neocortex in the ape lineage is the basis for the development of higher cognitive abilities. However, the human brain has uniquely increased in size and degree of folding, forming an essential foundation for advanced cognitive functions. This raises the question: what factors distinguish humans from our closest living primate relatives, such as chimpanzees and bonobos, which exhibit comparatively constrained cognitive capabilities? In this review, we focus on recent studies examining (modern) human-specific genetic traits that influence neural progenitor cells, whose behavior and activity are crucial for shaping cortical morphology. We emphasize the role of human-specific genetic modifications in signaling pathways that enhance the abundance of apical and basal progenitors, as well as the importance of basal progenitor metabolism in their proliferation in human. Additionally, we discuss how changes in neuron morphology contribute to the evolution of human cognition and provide our perspective on future directions in the field.

Meta-analysis of single-cell RNA sequencing co-expression in human neural organoids reveals their high variability in recapitulating primary tissue

Human neural organoids offer an exciting opportunity for studying inaccessible human-specific brain development; however, it remains unclear how precisely organoids recapitulate fetal/primary tissue biology. 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 data sets) and neural organoids (1.59 million cells, 173 data sets). We quantify the degree primary tissue cell type marker expression and co-expression are recapitulated in organoids across 10 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 indistinguishable from primary tissue, demonstrating a high degree of variability in biological fidelity among organoid systems. Our preserved co-expression framework provides cell type-specific measures of fidelity applicable to diverse neural organoids, offering a powerful tool for uncovering unifying axes of variation across heterogeneous neural organoid experiments.

Implementation and validation of single-cell genomics experiments in neuroscience

Single-cell or single-nucleus transcriptomics is a powerful tool for identifying cell types and cell states. However, hypotheses derived from these assays, including gene expression information, require validation, and their functional relevance needs to be established. The choice of validation depends on numerous factors. Here, we present types of orthogonal and functional validation experiment to strengthen preliminary findings obtained using single-cell and single-nucleus transcriptomics as well as the challenges and limitations of these approaches.

Hierarchical annotation of eQTLs by H-eQTL enables identification of genes with cell type-divergent regulation

While context-type-specific regulation of genes is largely determined by cis-regulatory regions, attempts to identify cell type-specific eQTLs are complicated by the nested nature of cell types. We present hierarchical eQTL (H-eQTL), a network-based model for hierarchical annotation of bulk-derived eQTLs to levels of a cell type tree using single-cell chromatin accessibility data and no clustering of cells into discrete cell types. Using our model, we annotate bulk-derived eQTLs from the developing brain with high specificity to levels of a cell type hierarchy, which allows sensitive detection of genes with multiple distinct non-coding elements regulating their expression in different cell types.

Molecular signatures of cortical expansion in the human foetal brain

The third trimester of human gestation is characterised by rapid increases in brain volume and cortical surface area. Recent studies have revealed a remarkable molecular diversity across the prenatal cortex but little is known about how this diversity translates into the differential rates of cortical expansion observed during gestation. We present a digital resource, μBrain, to facilitate knowledge translation between molecular and anatomical descriptions of the prenatal brain. Using μBrain, we evaluate the molecular signatures of preferentially-expanded cortical regions, quantified in utero using magnetic resonance imaging. Our findings demonstrate a spatial coupling between areal differences in the timing of neurogenesis and rates of neocortical expansion during gestation. We identify genes, upregulated from mid-gestation, that are highly expressed in rapidly expanding neocortex and implicated in genetic disorders with cognitive sequelae. The μBrain atlas provides a tool to comprehensively map early brain development across domains, model systems and resolution scales.

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.

FGF8-mediated gene regulation affects regional identity in human cerebral organoids

The morphogen FGF8 establishes graded positional cues imparting regional cellular responses via modulation of early target genes. The roles of FGF signaling and its effector genes remain poorly characterized in human experimental models mimicking early fetal telencephalic development. We used hiPSC-derived cerebral organoids as an in vitro platform to investigate the effect of FGF8 signaling on neural identity and differentiation. We found that FGF8 treatment increases cellular heterogeneity, leading to distinct telencephalic and mesencephalic-like domains that co-develop in multi-regional organoids. Within telencephalic regions, FGF8 affects the anteroposterior and dorsoventral identity of neural progenitors and the balance between GABAergic and glutamatergic neurons, thus impacting spontaneous neuronal network activity. Moreover, FGF8 efficiently modulates key regulators responsible for several human neurodevelopmental disorders. Overall, our results show that FGF8 signaling is directly involved in both regional patterning and cellular diversity in human cerebral organoids and in modulating genes associated with normal and pathological neural development.

Astrocyte regional specialization is shaped by postnatal development

Astrocytes are an abundant class of glial cells with critical roles in neural circuit assembly and function. Though many studies have uncovered significant molecular distinctions between astrocytes from different brain regions, how this regionalization unfolds over development is not fully understood. We used single-nucleus RNA sequencing to characterize the molecular diversity of brain cells across six developmental stages and four brain regions in the mouse and marmoset brain. Our analysis of over 170,000 single astrocyte nuclei revealed striking regional heterogeneity among astrocytes, particularly between telencephalic and diencephalic regions, at all developmental time points surveyed in both species. At the stages sampled, most of the region patterning was private to astrocytes and not shared with neurons or other glial types. Though astrocytes were already regionally patterned in late embryonic stages, this region-specific astrocyte gene expression signature changed dramatically over postnatal development, and its composition suggests that regional astrocytes further specialize postnatally to support their local neuronal circuits. Comparing across species, we found divergence in the expression of astrocytic region- and age-differentially expressed genes and the timing of astrocyte maturation relative to birth between mouse and marmoset, as well as hundreds of species differentially expressed genes. Finally, we used expansion microscopy to show that astrocyte morphology is largely conserved across gray matter regions of prefrontal cortex, striatum, and thalamus in the mouse, despite substantial molecular divergence.

The chromodomain protein CDYL confers forebrain identity to human cortical organoids by inhibiting neuronatin

Fate determination of neural stem cells (NSCs) is crucial for cortex development and is closely linked to neurodevelopmental disorders when gene expression networks are disrupted. The transcriptional corepressor chromodomain Y-like (CDYL) is widely expressed across diverse cell populations within the human embryonic cortex. However, its precise role in cortical development remains unclear. Here, we show that CDYL is critical for human cortical neurogenesis and that its deficiency leads to a substantial increase in gamma-aminobutyric acid (GABA)-ergic neurons in cortical organoids. Subsequently, neuronatin (NNAT) is identified as a significant target of CDYL, and its abnormal expression obviously influences the fate commitment of cortical NSCs. Cross-species comparisons of CDYL targets unravel a distinct developmental trajectory between human cortical organoids and the mouse cortex at an analogous stage. Collectively, our data provide insight into the evolutionary roles of CDYL in human cortex development, emphasizing its critical function in maintaining the fate of human cortical NSCs.

Neuronal lineage tracing from progenitors in human cortical organoids reveals mechanisms of neuronal production, diversity, and disease

The contribution of progenitor subtypes to generating the billions of neurons produced during human cortical neurogenesis is not well understood. We developed the cortical organoid lineage-tracing (COR-LT) system for human cortical organoids. Differential fluorescent reporter activation in distinct progenitor cells leads to permanent reporter expression, enabling the progenitor cell lineage of neurons to be determined. Surprisingly, nearly all excitatory neurons produced in cortical organoids were generated indirectly from intermediate progenitor cells. Additionally, neurons of different progenitor lineages were transcriptionally distinct. Isogenic lines made from an autistic individual with and without a likely pathogenic CTNNB1 variant demonstrated that the variant substantially altered the proportion of neurons derived from specific progenitor cell lineages, as well as the lineage-specific transcriptional profiles of these neurons, suggesting a pathogenic mechanism for this mutation. These results suggest individual progenitor subtypes play roles in generating the diverse neurons of the human cerebral cortex.

The Human Accelerated Region HAR202 Controls NPAS3 Expression in the Developing Forebrain Displaying Differential Enhancer Activity Between Modern and Archaic Human Sequences

It has been proposed that the phenotypic differences in cognitive abilities between humans and our closest living relatives, chimpanzees, are largely due to changes in the regulation of neurodevelopmental genes. We have previously found that the neurodevelopmental transcription factor gene NPAS3 accumulates the largest number of human accelerated regions (HARs), suggesting it may play some role in the phenotypic evolution of the human nervous system. In this work, we performed a comparative functional analysis of NPAS3-HAR202 using enhancer reporter assays in transgenic zebrafish and mice. We found that the Homo sapiens HAR202 ortholog failed to drive reporter expression to the zebrafish nervous system, in high contrast to the strong expression displayed by the rest of the vertebrate ortholog sequences tested. Remarkably, the HAR202 ortholog from archaic humans (Neanderthals/Denisovans) also displayed a pan-vertebrate expression pattern, despite the fact that archaic and modern humans have only one nucleotide substitution. Moreover, similar results were found when comparing enhancer activity in transgenic mice, where we observed a loss of activity of the modern human version in the mouse developing brain. To investigate the functional importance of HAR202, we generated mice lacking HAR202 and found a remarkable decrease of Npas3 expression in the forebrain during development. Our results place HAR202 as one of the very few examples of a neurodevelopmental transcriptional enhancer displaying functional evolution in the brain as a result of a fast molecular evolutionary process that specifically occurred in the human lineage.