Spatiotemporal gene expression trajectories reveal developmental hierarchies of the human cortex

The role of histone methyltransferases in neurocognitive disorders associated with brain size abnormalities

Brain size is controlled by several factors during neuronal development, including neural progenitor proliferation, neuronal arborization, gliogenesis, cell death, and synaptogenesis. Multiple neurodevelopmental disorders have co-morbid brain size abnormalities, such as microcephaly and macrocephaly. Mutations in histone methyltransferases that modify histone H3 on Lysine 36 and Lysine 4 (H3K36 and H3K4) have been identified in neurodevelopmental disorders involving both microcephaly and macrocephaly. H3K36 and H3K4 methylation are both associated with transcriptional activation and are proposed to sterically hinder the repressive activity of the Polycomb Repressor Complex 2 (PRC2). During neuronal development, tri-methylation of H3K27 (H3K27me3) by PRC2 leads to genome wide transcriptional repression of genes that regulate cell fate transitions and neuronal arborization. Here we provide a review of neurodevelopmental processes and disorders associated with H3K36 and H3K4 histone methyltransferases, with emphasis on processes that contribute to brain size abnormalities. Additionally, we discuss how the counteracting activities of H3K36 and H3K4 modifying enzymes vs. PRC2 could contribute to brain size abnormalities which is an underexplored mechanism in relation to brain size control.

Identification of brain cell types underlying genetic association with word reading and correlated traits

Neuroimaging studies implicate multiple cortical regions in reading ability/disability. However, the neural cell types integral to the reading process are unknown. To contribute to this gap in knowledge, we integrated genetic results from genome-wide association studies for word reading (n = 5054) with gene expression datasets from adult/fetal human brain. Linkage disequilibrium score regression (LDSC) suggested that variants associated with word reading were enriched in genes expressed in adult excitatory neurons, specifically layer 5 and 6 FEZF2 expressing neurons and intratelencephalic (IT) neurons, which express the marker genes LINC00507, THEMIS, or RORB. Inhibitory neurons (VIP, SST, and PVALB) were also found. This finding was interesting as neurometabolite studies previously implicated excitatory-inhibitory imbalances in the etiology of reading disabilities (RD). We also tested traits that shared genetic etiology with word reading (previously determined by polygenic risk scores): attention-deficit/hyperactivity disorder (ADHD), educational attainment, and cognitive ability. For ADHD, we identified enrichment in L4 IT adult excitatory neurons. For educational attainment and cognitive ability, we confirmed previous studies identifying multiple subclasses of adult cortical excitatory and inhibitory neurons, as well as astrocytes and oligodendrocytes. For educational attainment and cognitive ability, we also identified enrichment in multiple fetal cortical excitatory and inhibitory neurons, intermediate progenitor cells, and radial glial cells. In summary, this study supports a role of excitatory and inhibitory neurons in reading and excitatory neurons in ADHD and contributes new information on fetal cell types enriched in educational attainment and cognitive ability, thereby improving our understanding of the neurobiological basis of reading/correlated traits.

Comprehensive multi-omic profiling of somatic mutations in malformations of cortical development

Malformations of cortical development (MCD) are neurological conditions involving focal disruptions of cortical architecture and cellular organization that arise during embryogenesis, largely from somatic mosaic mutations, and cause intractable epilepsy. Identifying the genetic causes of MCD has been a challenge, as mutations remain at low allelic fractions in brain tissue resected to treat condition-related epilepsy. Here we report a genetic landscape from 283 brain resections, identifying 69 mutated genes through intensive profiling of somatic mutations, combining whole-exome and targeted-amplicon sequencing with functional validation including in utero electroporation of mice and single-nucleus RNA sequencing. Genotype-phenotype correlation analysis elucidated specific MCD gene sets associated with distinct pathophysiological and clinical phenotypes. The unique single-cell level spatiotemporal expression patterns of mutated genes in control and patient brains indicate critical roles in excitatory neurogenic pools during brain development and in promoting neuronal hyperexcitability after birth.

TMEM161B modulates radial glial scaffolding in neocortical development

TMEM161B encodes an evolutionarily conserved widely expressed novel 8-pass transmembrane protein of unknown function in human. Here we identify TMEM161B homozygous hypomorphic missense variants in our recessive polymicrogyria (PMG) cohort. Patients carrying TMEM161B mutations exhibit striking neocortical PMG and intellectual disability. Tmem161b knockout mice fail to develop midline hemispheric cleavage, whereas knock-in of patient mutations and patient-derived brain organoids show defects in apical cell polarity and radial glial scaffolding. We found that TMEM161B modulates actin filopodia, functioning upstream of the Rho-GTPase CDC42. Our data link TMEM161B with human PMG, likely regulating radial glia apical polarity during neocortical development.

Codon modification of Tuba1a alters mRNA levels and causes a severe neurodevelopmental phenotype in mice

The tubulinopathies are an umbrella of rare diseases that result from mutations in tubulin genes and are frequently characterised by severe brain malformations. The characteristics of a given disease reflect the expression pattern of the transcript, the function of a given tubulin gene, and the role microtubules play in a particular cell type. Mouse models have proved to be valuable tools that have provided insight into the molecular and cellular mechanisms that underlie the disease state. In this manuscript we compare two Tuba1a mouse models, both of which express wild-type TUBA1A protein but employ different codon usage. We show that modification of the Tuba1a mRNA sequence results in homozygous lethality and a severe neurodevelopmental phenotype. This is associated with a decrease in the number of post-mitotic neurons, PAX6 positive progenitors, and an increase in the number of apoptotic cells. We attribute this to a decrease in the stability of the modified Tuba1a transcript, and the absence of compensation by the other neurogenic tubulins. Our findings highlight the importance of maintaining the wild-type coding sequence when engineering mouse lines and the impact of synonymous genetic variation.

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

Glioblastoma multiforme (GBM) is an aggressive, heterogeneous grade IV brain tumor. Glioblastoma stem cells (GSCs) initiate the tumor and are known culprits of therapy resistance. Mounting evidence has demonstrated a regulatory role of long non-coding RNAs (lncRNAs) in various biological processes, including pluripotency, differentiation, and tumorigenesis. 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 human GBM tumors, GBM organoids, GSC-enriched GBM tumors, and developing human brains to identify lncRNAs highly expressed in GBM. To categorize GSC populations in the GBM tumors, we used the GSC marker genes SOX2, PROM1, FUT4, and L1CAM. We found three major GSC population clusters: radial glia, oligodendrocyte progenitor cells, and neurons. We found 10â€"100 lncRNAs significantly enriched in different GSC populations. We also validated the level of expression and localization of several GSC-enriched lncRNAs using qRT-PCR, single-molecule RNA FISH, and sub-cellular fractionation. We found that the radial glia GSC-enriched lncRNA PANTR1 is highly expressed in GSC lines and is localized to both the cytoplasmic and nuclear fractions. In contrast, the neuronal GSC-enriched lncRNAs LINC01563 and MALAT1 are highly enriched in the nuclear fraction of GSCs. Together, this study identified a panel of uncharacterized GSC-specific lncRNAs. These findings 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.

Identifying foetal forebrain interneurons as a target for monogenic autism risk factors and the polygenic 16p11.2 microdeletion

BACKGROUND

Autism spectrum condition or 'autism' is associated with numerous genetic risk factors including the polygenic 16p11.2 microdeletion. The balance between excitatory and inhibitory neurons in the cerebral cortex is hypothesised to be critical for the aetiology of autism making improved understanding of how risk factors impact on the development of these cells an important area of research. In the current study we aim to combine bioinformatics analysis of human foetal cerebral cortex gene expression data with anatomical and electrophysiological analysis of a 16p11.2^+/- rat model to investigate how genetic risk factors impact on inhibitory neuron development.

METHODS

We performed bioinformatics analysis of single cell transcriptomes from gestational week (GW) 8-26 human foetal prefrontal cortex and anatomical and electrophysiological analysis of 16p11.2^+/- rat cerebral cortex and hippocampus at post-natal day (P) 21.

RESULTS

We identified a subset of human interneurons (INs) first appearing at GW23 with enriched expression of a large fraction of risk factor transcripts including those expressed from the 16p11.2 locus. This suggests the hypothesis that these foetal INs are vulnerable to mutations causing autism. We investigated this in a rat model of the 16p11.2 microdeletion. We found no change in the numbers or position of either excitatory or inhibitory neurons in the somatosensory cortex or CA1 of 16p11.2^+/- rats but found that CA1 Sst INs were hyperexcitable with an enlarged axon initial segment, which was not the case for CA1 pyramidal cells.

LIMITATIONS

The human foetal gene expression data was acquired from cerebral cortex between gestational week (GW) 8 to 26. We cannot draw inferences about potential vulnerabilities to genetic autism risk factors for cells not present in the developing cerebral cortex at these stages. The analysis 16p11.2^+/- rat phenotypes reported in the current study was restricted to 3-week old (P21) animals around the time of weaning and to a single interneuron cell-type while in human 16p11.2 microdeletion carriers symptoms likely involve multiple cell types and manifest in the first few years of life and on into adulthood.

CONCLUSIONS

We have identified developing interneurons in human foetal cerebral cortex as potentially vulnerable to monogenic autism risk factors and the 16p11.2 microdeletion and report interneuron phenotypes in post-natal 16p11.2^+/- rats.

Single-Nuclei RNA Sequencing of 5 Regions of the Human Prenatal Brain Implicates Developing Neuron Populations in Genetic Risk for Schizophrenia

BACKGROUND

While a variety of evidence supports a prenatal component in schizophrenia, there are few data regarding the cell populations involved. We sought to identify cells of the human prenatal brain mediating genetic risk for schizophrenia by integrating cell-specific gene expression measures generated through single-nuclei RNA sequencing with recent large-scale genome-wide association study (GWAS) and exome sequencing data for the condition.

METHODS

Single-nuclei RNA sequencing was performed on 5 brain regions (frontal cortex, ganglionic eminence, hippocampus, thalamus, and cerebellum) from 3 fetuses from the second trimester of gestation. Enrichment of schizophrenia common variant genetic liability and rare damaging coding variation was assessed in relation to gene expression specificity within each identified cell population.

RESULTS

Common risk variants were prominently enriched within genes with high expression specificity for developing neuron populations within the frontal cortex, ganglionic eminence, and hippocampus. Enrichments were largely independent of genes expressed in neuronal populations of the adult brain that have been implicated in schizophrenia through the same methods. Genes containing an excess of rare damaging variants in schizophrenia had higher expression specificity for developing glutamatergic neurons of the frontal cortex and hippocampus that were also enriched for common variant liability.

CONCLUSIONS

We found evidence for a distinct contribution of prenatal neuronal development to genetic risk for schizophrenia, involving specific populations of developing neurons within the second-trimester fetal brain. Our study significantly advances the understanding of the neurodevelopmental origins of schizophrenia and provides a resource with which to investigate the prenatal antecedents of other psychiatric and neurologic disorders.

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

The Neuroscience Multi-Omic Archive: a BRAIN Initiative resource for single-cell transcriptomic and epigenomic data from the mammalian brain

Scalable technologies to sequence the transcriptomes and epigenomes of single cells are transforming our understanding of cell types and cell states. The Brain Research through Advancing Innovative Neurotechnologies (BRAIN) Initiative Cell Census Network (BICCN) is applying these technologies at unprecedented scale to map the cell types in the mammalian brain. In an effort to increase data FAIRness (Findable, Accessible, Interoperable, Reusable), the NIH has established repositories to make data generated by the BICCN and related BRAIN Initiative projects accessible to the broader research community. Here, we describe the Neuroscience Multi-Omic Archive (NeMO Archive; nemoarchive.org), which serves as the primary repository for genomics data from the BRAIN Initiative. Working closely with other BRAIN Initiative researchers, we have organized these data into a continually expanding, curated repository, which contains transcriptomic and epigenomic data from over 50 million brain cells, including single-cell genomic data from all of the major regions of the adult and prenatal human and mouse brains, as well as substantial single-cell genomic data from non-human primates. We make available several tools for accessing these data, including a searchable web portal, a cloud-computing interface for large-scale data processing (implemented on Terra, terra.bio), and a visualization and analysis platform, NeMO Analytics (nemoanalytics.org).

Molecular subtypes of ALS are associated with differences in patient prognosis

Amyotrophic Lateral Sclerosis (ALS) is a neurodegenerative disease with poorly understood clinical heterogeneity, underscored by significant differences in patient age at onset, symptom progression, therapeutic response, disease duration, and comorbidity presentation. We perform a patient stratification analysis to better understand the variability in ALS pathology, utilizing postmortem frontal and motor cortex transcriptomes derived from 208 patients. Building on the emerging role of transposable element (TE) expression in ALS, we consider locus-specific TEs as distinct molecular features during stratification. Here, we identify three unique molecular subtypes in this ALS cohort, with significant differences in patient survival. These results suggest independent disease mechanisms drive some of the clinical heterogeneity in ALS.

Early Regional Patterning in the Human Prefrontal Cortex Revealed by Laminar Dynamics of Deep Projection Neuron Markers

Early regional patterning and laminar position of cortical projection neurons is determined by activation and deactivation of transcriptional factors (TFs) and RNA binding proteins (RBPs) that regulate spatiotemporal framework of neurogenetic processes (proliferation, migration, aggregation, postmigratory differentiation, molecular identity acquisition, axonal growth, dendritic development, and synaptogenesis) within transient cellular compartments. Deep-layer projection neurons (DPN), subplate (SPN), and Cajal-Retzius neurons (CRN) are early-born cells involved in the establishment of basic laminar and regional cortical architecture; nonetheless, laminar dynamics of their molecular transcriptional markers remain underexplored. Here we aimed to analyze laminar dynamics of DPN markers, i.e., transcription factors TBR1, CTIP2, TLE4, SOX5, and RBP CELF1 on histological serial sections of the human frontal cortex between 7.5-15 postconceptional weeks (PCW) in reference to transient proliferative, migratory, and postmigratory compartments. The subtle signs of regional patterning were seen during the late preplate phase in the pattern of sublaminar organization of TBR1+/Reelin+ CRN and TBR1+ pioneering SPN. During the cortical plate (CP)-formation phase, TBR1+ neurons became radially aligned, forming continuity from a well-developed subventricular zone to CP showing clear lateral to medial regional gradients. The most prominent regional patterning was seen during the subplate formation phase (around 13 PCW) when a unique feature of the orbitobasal frontal cortex displays a "double plate" pattern. In other portions of the frontal cortex (lateral, dorsal, medial) deep portion of CP becomes loose and composed of TBR1+, CTIP2+, TLE4+, and CELF1+ neurons of layer six and later-born SPN, which later become constituents of the expanded SP (around 15 PCW). Overall, TFs and RBPs mark characteristic regional laminar dynamics of DPN, SPN, and CRN subpopulations during remarkably early fetal phases of the highly ordered association cortex development.

A transcription factor atlas of directed differentiation

Transcription factors (TFs) regulate gene programs, thereby controlling diverse cellular processes and cell states. To comprehensively understand TFs and the programs they control, we created a barcoded library of all annotated human TF splice isoforms (>3,500) and applied it to build a TF Atlas charting expression profiles of human embryonic stem cells (hESCs) overexpressing each TF at single-cell resolution. We mapped TF-induced expression profiles to reference cell types and validated candidate TFs for generation of diverse cell types, spanning all three germ layers and trophoblasts. Targeted screens with subsets of the library allowed us to create a tailored cellular disease model and integrate mRNA expression and chromatin accessibility data to identify downstream regulators. Finally, we characterized the effects of combinatorial TF overexpression by developing and validating a strategy for predicting combinations of TFs that produce target expression profiles matching reference cell types to accelerate cellular engineering efforts.

CROCCP2 acts as a human-specific modifier of cilia dynamics and mTOR signaling to promote expansion of cortical progenitors

The primary cilium is a central signaling component during embryonic development. Here we focus on CROCCP2, a hominid-specific gene duplicate from ciliary rootlet coiled coil (CROCC), also known as rootletin, that encodes the major component of the ciliary rootlet. We find that CROCCP2 is highly expressed in the human fetal brain and not in other primate species. CROCCP2 gain of function in the mouse embryonic cortex and human cortical cells and organoids results in decreased ciliogenesis and increased cortical progenitor amplification, particularly basal progenitors. CROCCP2 decreases ciliary dynamics by inhibition of the IFT20 ciliary trafficking protein, which then impacts neurogenesis through increased mTOR signaling. Loss of function of CROCCP2 in human cortical cells and organoids leads to increased ciliogenesis, decreased mTOR signaling, and impaired basal progenitor amplification. These data identify CROCCP2 as a human-specific modifier of cortical neurogenesis that acts through modulation of ciliary dynamics and mTOR signaling.

Neurodevelopmental disorders-high-resolution rethinking of disease modeling

Neurodevelopmental disorders arise due to various risk factors that can perturb different stages of brain development, and a combinatorial impact of these risk factors programs the phenotype in adulthood. While modeling the complete phenotype of a neurodevelopmental disorder is challenging, individual developmental perturbations can be successfully modeled in vivo in animals and in vitro in human cellular models. Nevertheless, our limited knowledge of human brain development restricts modeling strategies and has raised questions of how well a model corresponds to human in vivo brain development. Recent progress in high-resolution analysis of human tissue with single-cell and spatial omics techniques has enhanced our understanding of the complex events that govern the development of the human brain in health and disease. This new knowledge can be utilized to improve modeling of neurodevelopmental disorders and pave the way to more accurately portraying the relevant developmental perturbations in disease models.

Transcriptional networks of transient cell states during human prefrontal cortex development

The human brain is divided into various anatomical regions that control and coordinate unique functions. The prefrontal cortex (PFC) is a large brain region that comprises a range of neuronal and non-neuronal cell types, sharing extensive interconnections with subcortical areas, and plays a critical role in cognition and memory. A timely appearance of distinct cell types through embryonic development is crucial for an anatomically perfect and functional brain. Direct tracing of cell fate development in the human brain is not possible, but single-cell transcriptome sequencing (scRNA-seq) datasets provide the opportunity to dissect cellular heterogeneity and its molecular regulators. Here, using scRNA-seq data of human PFC from fetal stages, we elucidate distinct transient cell states during PFC development and their underlying gene regulatory circuitry. We further identified that distinct intermediate cell states consist of specific gene regulatory modules essential to reach terminal fate using discrete developmental paths. Moreover, using in silico gene knock-out and over-expression analysis, we validated crucial gene regulatory components during the lineage specification of oligodendrocyte progenitor cells. Our study illustrates unique intermediate states and specific gene interaction networks that warrant further investigation for their functional contribution to typical brain development and discusses how this knowledge can be harvested for therapeutic intervention in challenging neurodevelopmental disorders.

Biology-inspired data-driven quality control for scientific discovery in single-cell transcriptomics

BACKGROUND

Quality control (QC) of cells, a critical first step in single-cell RNA sequencing data analysis, has largely relied on arbitrarily fixed data-agnostic thresholds applied to QC metrics such as gene complexity and fraction of reads mapping to mitochondrial genes. The few existing data-driven approaches perform QC at the level of samples or studies without accounting for biological variation.

RESULTS

We first demonstrate that QC metrics vary with both tissue and cell types across technologies, study conditions, and species. We then propose data-driven QC (ddqc), an unsupervised adaptive QC framework to perform flexible and data-driven QC at the level of cell types while retaining critical biological insights and improved power for downstream analysis. ddqc applies an adaptive threshold based on the median absolute deviation on four QC metrics (gene and UMI complexity, fraction of reads mapping to mitochondrial and ribosomal genes). ddqc retains over a third more cells when compared to conventional data-agnostic QC filters. Finally, we show that ddqc recovers biologically meaningful trends in gradation of gene complexity among cell types that can help answer questions of biological interest such as which cell types express the least and most number of transcripts overall, and ribosomal transcripts specifically.

CONCLUSIONS

ddqc retains cell types such as metabolically active parenchymal cells and specialized cells such as neutrophils which are often lost by conventional QC. Taken together, our work proposes a revised paradigm to quality filtering best practices-iterative QC, providing a data-driven QC framework compatible with observed biological diversity.

Temporal and sequential transcriptional dynamics define lineage shifts in corticogenesis

The cerebral cortex contains billions of neurons, and their disorganization or misspecification leads to neurodevelopmental disorders. Understanding how the plethora of projection neuron subtypes are generated by cortical neural stem cells (NSCs) is a major challenge. Here, we focused on elucidating the transcriptional landscape of murine embryonic NSCs, basal progenitors (BPs), and newborn neurons (NBNs) throughout cortical development. We uncover dynamic shifts in transcriptional space over time and heterogeneity within each progenitor population. We identified signature hallmarks of NSC, BP, and NBN clusters and predict active transcriptional nodes and networks that contribute to neural fate specification. We find that the expression of receptors, ligands, and downstream pathway components is highly dynamic over time and throughout the lineage implying differential responsiveness to signals. Thus, we provide an expansive compendium of gene expression during cortical development that will be an invaluable resource for studying neural developmental processes and neurodevelopmental disorders.

An atlas of late prenatal human neurodevelopment resolved by single-nucleus transcriptomics

Late prenatal development of the human neocortex encompasses a critical period of gliogenesis and cortical expansion. However, systematic single-cell analyses to resolve cellular diversity and gliogenic lineages of the third trimester are lacking. Here, we present a comprehensive single-nucleus RNA sequencing atlas of over 200,000 nuclei derived from the proliferative germinal matrix and laminating cortical plate of 15 prenatal, non-pathological postmortem samples from 17 to 41 gestational weeks, and 3 adult controls. This dataset captures prenatal gliogenesis with high temporal resolution and is provided as a resource for further interrogation. Our computational analysis resolves greater complexity of glial progenitors, including transient glial intermediate progenitor cell (gIPC) and nascent astrocyte populations in the third trimester of human gestation. We use lineage trajectory and RNA velocity inference to further characterize specific gIPC subpopulations preceding both oligodendrocyte (gIPC-O) and astrocyte (gIPC-A) lineage differentiation. We infer unique transcriptional drivers and biological pathways associated with each developmental state, validate gIPC-A and gIPC-O presence within the human germinal matrix and cortical plate in situ, and demonstrate gIPC states being recapitulated across adult and pediatric glioblastoma tumors.

RINGs, DUBs and Abnormal Brain Growth-Histone H2A Ubiquitination in Brain Development and Disease

During mammalian neurodevelopment, signaling pathways converge upon transcription factors (TFs) to establish appropriate gene expression programmes leading to the production of distinct neural and glial cell types. This process is partially regulated by the dynamic modulation of chromatin states by epigenetic systems, including the polycomb group (PcG) family of co-repressors. PcG proteins form multi-subunit assemblies that sub-divide into distinct, yet functionally related families. Polycomb repressive complexes 1 and 2 (PRC1 and 2) modify the chemical properties of chromatin by covalently modifying histone tails via H2A ubiquitination (H2AK119ub1) and H3 methylation, respectively. In contrast to the PRCs, the Polycomb repressive deubiquitinase (PR-DUB) complex removes H2AK119ub1 from chromatin through the action of the C-terminal hydrolase BAP1. Genetic screening has identified several PcG mutations that are causally associated with a range of congenital neuropathologies associated with both localised and/or systemic growth abnormalities. As PRC1 and PR-DUB hold opposing functions to control H2AK119ub1 levels across the genome, it is plausible that such neurodevelopmental disorders arise through a common mechanism. In this review, we will focus on advancements regarding the composition and opposing molecular functions of mammalian PRC1 and PR-DUB, and explore how their dysfunction contributes to the emergence of neurodevelopmental disorders.

How mechanisms of stem cell polarity shape the human cerebral cortex

Apical-basal progenitor cell polarity establishes key features of the radial and laminar architecture of the developing human cortex. The unique diversity of cortical stem cell populations and an expansion of progenitor population size in the human cortex have been mirrored by an increase in the complexity of cellular processes that regulate stem cell morphology and behaviour, including their polarity. The study of human cells in primary tissue samples and human stem cell-derived model systems (such as cortical organoids) has provided insight into these processes, revealing that protein complexes regulate progenitor polarity by controlling cell membrane adherence within appropriate cortical niches and are themselves regulated by cytoskeletal proteins, signalling molecules and receptors, and cellular organelles. Studies exploring how cortical stem cell polarity is established and maintained are key for understanding the features of human brain development and have implications for neurological dysfunction.

Thyroid hormone regulators in human cerebral cortex development

Brain development is critically dependent on the timely supply of thyroid hormones. The thyroid hormone transporters are central to the action of thyroid hormones in the brain, facilitating their passage through the blood-brain barrier. Mutations of the monocarboxylate transporter 8 (MCT8) cause the Allan-Herndon-Dudley syndrome, with altered thyroid hormone concentrations in the blood and profound neurological impairment and intellectual deficit. Mouse disease models have revealed interplay between transport, deiodination, and availability of T3 to receptors in specific cells. However, the mouse models are not satisfactory, given the fundamental differences between the mouse and human brains. The goal of the present work is to review human neocortex development in the context of thyroid pathophysiology. Recent developments in single-cell transcriptomic approaches aimed at the human brain make it possible to profile the expression of thyroid hormone regulators in single-cell RNA-Seq datasets of the developing human neocortex. The data provide novel insights into the specific cellular expression of thyroid hormone transporters, deiodinases, and receptors.

Adaptive sequence divergence forged new neurodevelopmental enhancers in humans

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

Global hypo-methylation in a proportion of glioblastoma enriched for an astrocytic signature is associated with increased invasion and altered immune landscape

We describe a subset of glioblastoma, the most prevalent malignant adult brain tumour, harbouring a bias towards hypomethylation at defined differentially methylated regions. This epigenetic signature correlates with an enrichment for an astrocytic gene signature, which together with the identification of enriched predicted binding sites of transcription factors known to cause demethylation and to be involved in astrocytic/glial lineage specification, point to a shared ontogeny between these glioblastomas and astroglial progenitors. At functional level, increased invasiveness, at least in part mediated by SRPX2, and macrophage infiltration characterise this subset of glioblastoma.

Integrated gene analyses of de novo variants from 46,612 trios with autism and developmental disorders

Most genetic studies consider autism spectrum disorder (ASD) and developmental disorder (DD) separately despite overwhelming comorbidity and shared genetic etiology. Here, we analyzed de novo variants (DNVs) from 15,560 ASD (6,557 from SPARK) and 31,052 DD trios independently and also combined as broader neurodevelopmental disorders (NDDs) using three models. We identify 615 NDD candidate genes (false discovery rate [FDR] < 0.05) supported by ≥1 models, including 138 reaching Bonferroni exome-wide significance (P < 3.64e-7) in all models. The genes group into five functional networks associating with different brain developmental lineages based on single-cell nuclei transcriptomic data. We find no evidence for ASD-specific genes in contrast to 18 genes significantly enriched for DD. There are 53 genes that show mutational bias, including enrichments for missense (n = 41) or truncating (n = 12) DNVs. We also find 10 genes with evidence of male- or female-bias enrichment, including 4 X chromosome genes with significant female burden (DDX3X, MECP2, WDR45, and HDAC8). This large-scale integrative analysis identifies candidates and functional subsets of NDD genes.

Human prefrontal cortex gene regulatory dynamics from gestation to adulthood at single-cell resolution

Human brain development is underpinned by cellular and molecular reconfigurations continuing into the third decade of life. To reveal cell dynamics orchestrating neural maturation, we profiled human prefrontal cortex gene expression and chromatin accessibility at single-cell resolution from gestation to adulthood. Integrative analyses define the dynamic trajectories of each cell type, revealing major gene expression reconfiguration at the prenatal-to-postnatal transition in all cell types followed by continuous reconfiguration into adulthood and identifying regulatory networks guiding cellular developmental programs, states, and functions. We uncover links between expression dynamics and developmental milestones, characterize the diverse timing of when cells acquire adult-like states, and identify molecular convergence from distinct developmental origins. We further reveal cellular dynamics and their regulators implicated in neurological disorders. Finally, using this reference, we benchmark cell identities and maturation states in organoid models. Together, this captures the dynamic regulatory landscape of human cortical development.

Aberrant Cortical Layer Development of Brain Organoids Derived from Noonan Syndrome-iPSCs

Noonan syndrome (NS) is a genetic disorder mainly caused by gain-of-function mutations in Src homology region 2-containing protein tyrosine phosphatase 2 (SHP2). Although diverse neurological manifestations are commonly diagnosed in NS patients, the mechanisms as to how SHP2 mutations induce the neurodevelopmental defects associated with NS remain elusive. Here, we report that cortical organoids (NS-COs) derived from NS-induced pluripotent stem cells (iPSCs) exhibit developmental abnormalities, especially in excitatory neurons (ENs). Although NS-COs develop normally in their appearance, single-cell transcriptomic analysis revealed an increase in the EN population and overexpression of cortical layer markers in NS-COs. Surprisingly, the EN subpopulation co-expressing the upper layer marker SATB2 and the deep layer maker CTIP2 was enriched in NS-COs during cortical development. In parallel with the developmental disruptions, NS-COs also exhibited reduced synaptic connectivity. Collectively, our findings suggest that perturbed cortical layer identity and impeded neuronal connectivity contribute to the neurological manifestations of NS.

Single-nucleus transcriptomics of IDH1- and TP53-mutant glioma stem cells displays diversified commitment on invasive cancer progenitors

Glioma is a devastating brain tumor with a high mortality rate attributed to the glioma stem cells (GSCs) possessing high plasticity. Marker mutations in isocitrate dehydrogenase type 1 (IDH1) and tumor protein 53 (TP53) are frequent in gliomas and impact the cell fate decisions. Understanding the GSC heterogeneity within IDH1- and TP53- mutant tumors may elucidate possible treatment targets. Here, we performed single-nucleus transcriptomics of mutant and wild-type glioma samples sorted for Sox2 stem cell marker. For the first time the rare subpopulations of Sox2 + IDH1- and TP53-mutant GSCs were characterized. In general, GSCs contained the heterogeneity root subpopulation resembling active neural stem cells capable of asymmetric division to quiescent and transit amplifying cell branches. Specifically, double-mutant GSCs revealed the commitment on highly invasive oligodendrocyte- and astroglia-like progenitors. Additionally, double-mutant GSCs displayed upregulated markers of collagen synthesis, altered lipogenesis and high migration, while wild-type GSCs expressed genes related to ATP production. Wild-type GSC root population was highly heterogeneous and lacked the signature marker expression, thus glioblastoma treatment should emphasize on establishing differentiation protocol directed against residual GSCs. For the more differentiated IDH1- and TP53-mutant gliomas we suggest therapeutic targeting of migration molecules, such as CD44.

Spatially resolved gene regulatory and disease-related vulnerability map of the adult Macaque cortex

Single cell approaches have increased our knowledge about the cell type composition of the non-human primate (NHP), but a detailed characterization of area-specific regulatory features remains outstanding. We generated single-cell transcriptomic and chromatin accessibility (single-cell ATAC) data of 358,237 cells from prefrontal cortex (PFC), primary motor cortex (M1) and primary visual cortex (V1) of adult female cynomolgus monkey brain, and integrated this dataset with Stereo-seq (spatial enhanced resolution omics-sequencing) of the corresponding cortical areas to assign topographic information to molecular states. We identified area-specific chromatin accessible sites and their targeted genes, including the cell type-specific transcriptional regulatory network associated with excitatory neurons heterogeneity. We reveal calcium ion transport and axon guidance genes related to specialized functions of PFC and M1, identified the similarities and differences between adult macaque and human oligodendrocyte trajectories, and mapped the genetic variants and gene perturbations of human diseases to NHP cortical cells. This resource establishes a transcriptomic and chromatin accessibility combinatory regulatory landscape at a single-cell and spatially resolved resolution in NHP cortex.

Sex-biased effects on hippocampal circuit development by perinatal SERT expression in CA3 pyramidal neurons

Neurodevelopmental disorders ranging from autism to intellectual disability display sex-biased prevalence and phenotypical presentations. Despite increasing knowledge about temporospatial cortical map development and genetic variants linked to neurodevelopmental disorders, when and how sex-biased neural circuit derailment may arise in diseased brain remain unknown. Here, we identify in mice that serotonin uptake transporter (SERT) in non-serotonergic neurons - hippocampal and prefrontal pyramidal neurons - confers sex-biased effects specifically during neural circuit development. A set of gradient-patterned CA3 pyramidal neurons transiently express SERT to clear extracellular serotonin, coinciding with hippocampal synaptic circuit establishment. Ablating pyramidal neuron SERT (SERTPyramidΔ) alters dendritic spine developmental trajectory in the hippocampus, and precipitates sex-biased impairments in long-term activity-dependent hippocampal synaptic plasticity and cognitive behaviors. Transcriptomic analyses identify sex-biased alterations in gene sets associated with autism, dendritic spine structure, synaptic function and male-specific enrichment of dysregulated genes in glial cells in early postnatal SERTPyramidΔ hippocampus. Our data suggest that SERT function in these pyramidal neurons underscores a temporal- and brain region-specific regulation of normal sex-dimorphic circuit development and a source for sex-biased vulnerability to cognitive and behavioral impairments. This article has an associated 'The people behind the papers' interview.

ZEB2 haploinsufficient Mowat-Wilson syndrome induced pluripotent stem cells show disrupted GABAergic transcriptional regulation and function

Mowat-Wilson syndrome (MWS) is a severe neurodevelopmental disorder caused by heterozygous variants in the gene encoding transcription factor ZEB2. Affected individuals present with structural brain abnormalities, speech delay and epilepsy. In mice, conditional loss of Zeb2 causes hippocampal degeneration, altered migration and differentiation of GABAergic interneurons, a heterogeneous population of mainly inhibitory neurons of importance for maintaining normal excitability. To get insights into GABAergic development and function in MWS we investigated ZEB2 haploinsufficient induced pluripotent stem cells (iPSC) of MWS subjects together with iPSC of healthy donors. Analysis of RNA-sequencing data at two time points of GABAergic development revealed an attenuated interneuronal identity in MWS subject derived iPSC with enrichment of differentially expressed genes required for transcriptional regulation, cell fate transition and forebrain patterning. The ZEB2 haploinsufficient neural stem cells (NSCs) showed downregulation of genes required for ventral telencephalon specification, such as FOXG1, accompanied by an impaired migratory capacity. Further differentiation into GABAergic interneuronal cells uncovered upregulation of transcription factors promoting pallial and excitatory neurons whereas cortical markers were downregulated. The differentially expressed genes formed a neural protein-protein network with extensive connections to well-established epilepsy genes. Analysis of electrophysiological properties in ZEB2 haploinsufficient GABAergic cells revealed overt perturbations manifested as impaired firing of repeated action potentials. Our iPSC model of ZEB2 haploinsufficient GABAergic development thus uncovers a dysregulated gene network leading to immature interneurons with mixed identity and altered electrophysiological properties, suggesting mechanisms contributing to the neuropathogenesis and seizures in MWS.

Loss of non-motor kinesin KIF26A causes congenital brain malformations via dysregulated neuronal migration and axonal growth as well as apoptosis

Kinesins are canonical molecular motors but can also function as modulators of intracellular signaling. KIF26A, an unconventional kinesin that lacks motor activity, inhibits growth-factor-receptor-bound protein 2 (GRB2)- and focal adhesion kinase (FAK)-dependent signal transduction, but its functions in the brain have not been characterized. We report a patient cohort with biallelic loss-of-function variants in KIF26A, exhibiting a spectrum of congenital brain malformations. In the developing brain, KIF26A is preferentially expressed during early- and mid-gestation in excitatory neurons. Combining mice and human iPSC-derived organoid models, we discovered that loss of KIF26A causes excitatory neuron-specific defects in radial migration, localization, dendritic and axonal growth, and apoptosis, offering a convincing explanation of the disease etiology in patients. Single-cell RNA sequencing in KIF26A knockout organoids revealed transcriptional changes in MAPK, MYC, and E2F pathways. Our findings illustrate the pathogenesis of KIF26A loss-of-function variants and identify the surprising versatility of this non-motor kinesin.

High-throughput CRISPRi and CRISPRa technologies in 3D genome regulation for neuropsychiatric diseases

Advances in genomics have led to the identification of many risk loci with hundreds of genes and thousands of DNA variants associated with neuropsychiatric disorders. A significant barrier to understanding the genetic underpinnings of complex diseases is the lack of functional characterization of risk genes and variants in biological systems relevant to human health and connecting disease-associated variants to pathological phenotypes. Characterizing gene and DNA variant functions requires genetic perturbations followed by molecular and cellular assays of neurobiological phenotypes. However, generating null or mutant alleles is low throughput, making it impossible to characterize disease-associated variants in large quantities efficiently. CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa) screens can be leveraged to dissect the biological consequences of the tested genes and variants in their native context. Nevertheless, testing non-coding variants associated with complex diseases remains non-trivial. In this review, we first discuss the current challenges of interpreting the function of the non-coding genome and approaches to prioritizing disease-associated variants in the context of the 3D epigenome. Second, we provide a brief overview of high-throughput CRISPRi and CRISPRa screening strategies applicable for characterizing non-coding sequences in appropriate biological systems. Lastly, we discuss the promising prospects of using CRISPR-based technologies to dissect DNA sequences associated with neuropsychiatric diseases.

Two-stage linked component analysis for joint decomposition of multiple biologically related data sets

Integrative analysis of multiple data sets has the potential of fully leveraging the vast amount of high throughput biological data being generated. In particular such analysis will be powerful in making inference from publicly available collections of genetic, transcriptomic and epigenetic data sets which are designed to study shared biological processes, but which vary in their target measurements, biological variation, unwanted noise, and batch variation. Thus, methods that enable the joint analysis of multiple data sets are needed to gain insights into shared biological processes that would otherwise be hidden by unwanted intra-data set variation. Here, we propose a method called two-stage linked component analysis (2s-LCA) to jointly decompose multiple biologically related experimental data sets with biological and technological relationships that can be structured into the decomposition. The consistency of the proposed method is established and its empirical performance is evaluated via simulation studies. We apply 2s-LCA to jointly analyze four data sets focused on human brain development and identify meaningful patterns of gene expression in human neurogenesis that have shared structure across these data sets.

Modeling human telencephalic development and autism-associated SHANK3 deficiency using organoids generated from single neural rosettes

Human telencephalon is an evolutionarily advanced brain structure associated with many uniquely human behaviors and disorders. However, cell lineages and molecular pathways implicated in human telencephalic development remain largely unknown. We produce human telencephalic organoids from stem cell-derived single neural rosettes and investigate telencephalic development under normal and pathological conditions. We show that single neural rosette-derived organoids contain pallial and subpallial neural progenitors, excitatory and inhibitory neurons, as well as macroglial and periendothelial cells, and exhibit predictable organization and cytoarchitecture. We comprehensively characterize the properties of neurons in SNR-derived organoids and identify transcriptional programs associated with the specification of excitatory and inhibitory neural lineages from a common pool of NPs early in telencephalic development. We also demonstrate that neurons in organoids with a hemizygous deletion of an autism- and intellectual disability-associated gene SHANK3 exhibit intrinsic and excitatory synaptic deficits and impaired expression of several clustered protocadherins. Collectively, this study validates SNR-derived organoids as a reliable model for studying human telencephalic cortico-striatal development and identifies intrinsic, synaptic, and clustered protocadherin expression deficits in human telencephalic tissue with SHANK3 hemizygosity.

Orchestrating human neocortex development across the scales; from micro to macro

How our brains have developed to perform the many complex functions that make us human has long remained a question of great interest. Over the last few decades, many scientists from a wide range of fields have tried to answer this question by aiming to uncover the mechanisms that regulate the development of the human neocortex. They have approached this on different scales, focusing microscopically on individual cells all the way up to macroscopically imaging entire brains within living patients. In this review we will summarise these key findings and how they fit together.

Maturation and circuit integration of transplanted human cortical organoids

Self-organizing neural organoids represent a promising in vitro platform with which to model human development and disease1-5. However, organoids lack the connectivity that exists in vivo, which limits maturation and makes integration with other circuits that control behaviour impossible. Here we show that human stem cell-derived cortical organoids transplanted into the somatosensory cortex of newborn athymic rats develop mature cell types that integrate into sensory and motivation-related circuits. MRI reveals post-transplantation organoid growth across multiple stem cell lines and animals, whereas single-nucleus profiling shows progression of corticogenesis and the emergence of activity-dependent transcriptional programs. Indeed, transplanted cortical neurons display more complex morphological, synaptic and intrinsic membrane properties than their in vitro counterparts, which enables the discovery of defects in neurons derived from individuals with Timothy syndrome. Anatomical and functional tracings show that transplanted organoids receive thalamocortical and corticocortical inputs, and in vivo recordings of neural activity demonstrate that these inputs can produce sensory responses in human cells. Finally, cortical organoids extend axons throughout the rat brain and their optogenetic activation can drive reward-seeking behaviour. Thus, transplanted human cortical neurons mature and engage host circuits that control behaviour. We anticipate that this approach will be useful for detecting circuit-level phenotypes in patient-derived cells that cannot otherwise be uncovered.

Characterization of exogenous αSN response genes and their relation to Parkinson’s disease using network analyses

Despite extensive research, the molecular mechanisms underlying the toxicity of αSN in Parkinson’s disease (PD) pathology are still poorly understood. To address this, we used a microarray dataset to identify genes that are induced and differentially expressed after exposure to toxic αSN aggregates, which we call exogenous αSN response (EASR) genes. Using systems biology approaches, we then determined, at multiple levels of analysis, how these EASR genes could be related to PD pathology. A key result was the identification of functional connections between EASR genes and previously identified PD-related genes by employing the proteins’ interactions networks and 9 brain region-specific co-expression networks. In each brain region, co-expression modules of EASR genes were enriched for gene sets whose expression are altered by SARS-CoV-2 infection, leading to the hypothesis that EASR co-expression genes may explain the observed links between COVID-19 and PD. An examination of the expression pattern of EASR genes in different non-neurological healthy brain regions revealed that regions with lower mean expression of the upregulated EASR genes, such as substantia nigra, are more vulnerable to αSN aggregates and lose their neurological functions during PD progression. Gene Set Enrichment Analysis of healthy and PD samples from substantia nigra revealed that a specific co-expression network, “TNF-α signaling via NF-κB”, is an upregulated pathway associated with the PD phenotype. Inhibitors of the “TNF-α signaling via NF-κB” pathway may, therefore, decrease the activity level of this pathway and thereby provide therapeutic benefits for PD patients. We virtually screened FDA-approved drugs against these upregulated genes (NR4A1, DUSP1, and FOS) using docking-based drug discovery and identified several promising drugs. Altogether, our study provides a better understanding of αSN toxicity mechanisms in PD and identifies potential therapeutic targets and small molecules for treatment of PD.

Ensembles of endothelial and mural cells promote angiogenesis in prenatal human brain

Interactions between angiogenesis and neurogenesis regulate embryonic brain development. However, a comprehensive understanding of the stages of vascular cell maturation is lacking, especially in the prenatal human brain. Using fluorescence-activated cell sorting, single-cell transcriptomics, and histological and ultrastructural analyses, we show that an ensemble of endothelial and mural cell subtypes tile the brain vasculature during the second trimester. These vascular cells follow distinct developmental trajectories and utilize diverse signaling mechanisms, including collagen, laminin, and midkine, to facilitate cell-cell communication and maturation. Interestingly, our results reveal that tip cells, a subtype of endothelial cells, are highly enriched near the ventricular zone, the site of active neurogenesis. Consistent with these observations, prenatal vascular cells transplanted into cortical organoids exhibit restricted lineage potential that favors tip cells, promotes neurogenesis, and reduces cellular stress. Together, our results uncover important mechanisms into vascular maturation during this critical period of human brain development.

Cerebrovasculature pumps up progenitors

Fetal human brain stem cell niches that contain multipotent neural progenitors are progressively vascularized during development. Crouch et al. (Crouch et al., 2022) report endothelial and mural lineage trajectories that build developing prenatal vascular in second trimester fetal brain. This cerebral angiogenesis in neural progenitor zones occurs simultaneously with and can promote neurogenesis.

Microfluidics for Neuronal Cell and Circuit Engineering

The widespread adoption of microfluidic devices among the neuroscience and neurobiology communities has enabled addressing a broad range of questions at the molecular, cellular, circuit, and system levels. Here, we review biomedical engineering approaches that harness the power of microfluidics for bottom-up generation of neuronal cell types and for the assembly and analysis of neural circuits. Microfluidics-based approaches are instrumental to generate the knowledge necessary for the derivation of diverse neuronal cell types from human pluripotent stem cells, as they enable the isolation and subsequent examination of individual neurons of interest. Moreover, microfluidic devices allow to engineer neural circuits with specific orientations and directionality by providing control over neuronal cell polarity and permitting the isolation of axons in individual microchannels. Similarly, the use of microfluidic chips enables the construction not only of 2D but also of 3D brain, retinal, and peripheral nervous system model circuits. Such brain-on-a-chip and organoid-on-a-chip technologies are promising platforms for studying these organs as they closely recapitulate some aspects of in vivo biological processes. Microfluidic 3D neuronal models, together with 2D in vitro systems, are widely used in many applications ranging from drug development and toxicology studies to neurological disease modeling and personalized medicine. Altogether, microfluidics provide researchers with powerful systems that complement and partially replace animal models.

Recent Developments in Autism Genetic Research: A Scientometric Review from 2018 to 2022

Genetic research in Autism Spectrum Disorder (ASD) has progressed tremendously in recent decades. Dozens of genetic loci and hundreds of alterations in the genetic sequence, expression, epigenetic transformation, and interactions with other physiological and environmental systems have been found to increase the likelihood of developing ASD. There is therefore a need to represent this wide-ranging yet voluminous body of literature in a systematic manner so that this information can be synthesised and understood at a macro level. Therefore, this study made use of scientometric methods, particularly document co-citation analysis (DCA), to systematically review literature on ASD genetic research from 2018 to 2022. A total of 14,818 articles were extracted from Scopus and analyzed with CiteSpace. An optimized DCA analysis revealed that recent literature on ASD genetic research can be broadly organised into 12 major clusters representing various sub-topics. These clusters are briefly described in the manuscript and potential applications of this study are discussed.

Cerebral Organoids as an Experimental Platform for Human Neurogenomics

The cerebral cortex forms early in development according to a series of heritable neurodevelopmental instructions. Despite deep evolutionary conservation of the cerebral cortex and its foundational six-layered architecture, significant variations in cortical size and folding can be found across mammals, including a disproportionate expansion of the prefrontal cortex in humans. Yet our mechanistic understanding of neurodevelopmental processes is derived overwhelmingly from rodent models, which fail to capture many human-enriched features of cortical development. With the advent of pluripotent stem cells and technologies for differentiating three-dimensional cultures of neural tissue in vitro, cerebral organoids have emerged as an experimental platform that recapitulates several hallmarks of human brain development. In this review, we discuss the merits and limitations of cerebral organoids as experimental models of the developing human brain. We highlight innovations in technology development that seek to increase its fidelity to brain development in vivo and discuss recent efforts to use cerebral organoids to study regeneration and brain evolution as well as to develop neurological and neuropsychiatric disease models.

Identifying enhancer properties associated with genetic risk for complex traits using regulome-wide association studies

Genetic risk for complex traits is strongly enriched in non-coding genomic regions involved in gene regulation, especially enhancers. However, we lack adequate tools to connect the characteristics of these disruptions to genetic risk. Here, we propose RWAS (Regulome Wide Association Study), a new application of the MAGMA software package to identify the characteristics of enhancers that contribute to genetic risk for disease. RWAS involves three steps: (i) assign genotyped SNPs to cell type- or tissue-specific regulatory features (e.g., enhancers); (ii) test associations of each regulatory feature with a trait of interest for which genome-wide association study (GWAS) summary statistics are available; (iii) perform enhancer-set enrichment analyses to identify quantitative or categorical features of regulatory elements that are associated with the trait. These steps are implemented as a novel application of MAGMA, a tool originally developed for gene-based GWAS analyses. Applying RWAS to interrogate genetic risk for schizophrenia, we discovered a class of risk-associated AT-rich enhancers that are active in the developing brain and harbor binding sites for multiple transcription factors with neurodevelopmental functions. RWAS utilizes open-source software, and we provide a comprehensive collection of annotations for tissue-specific enhancer locations and features, including their evolutionary conservation, AT content, and co-localization with binding sites for hundreds of TFs. RWAS will enable researchers to characterize properties of regulatory elements associated with any trait of interest for which GWAS summary statistics are available.

Corticogenesis across species at single-cell resolution

The neocortex (or pallium) consists of diverse cell types that are organized in a highly species-specific manner under strict spatiotemporal control during development. Many of the cell types are present transiently throughout development but contribute to permanent species-specific cortical features that are acquired through evolution. Therefore, capturing cell type-specific biological information has always been an important quest in the field of neurodevelopment. The progress in achieving fine cellular resolution has been slow due to technical challenges. However, with recent advancements in single-cell and multi-omics technologies, many laboratories have begun to successfully interrogate cellular and molecular mechanisms driving corticogenesis at single-cell resolution. In this review, we provide summarized results from many primary publications and several in-depth review articles that utilize or address single-cell genomics techniques to understand important topics, such as cellular and molecular mechanisms governing cortical progenitor proliferation, cell lineage progression, neuronal specification, and arealization, across multiple gyrencephalic (i.e., human and non-human primates) and lissencephalic species (i.e., mouse, reptiles, and songbirds). We also examine findings from recent studies involving epigenomic and posttranscriptional regulation of corticogenesis. In the discussion section, we provide our insights on the challenges the field currently faces as well as promising future applications of single cell technologies.

Pathway-based approach reveals differential sensitivity to E2F1 inhibition in glioblastoma

Analysis of tumor gene expression is an important approach for the classification and identification of therapeutic vulnerabilities. However, targeting glioblastoma (GBM) based on molecular subtyping has not yet translated into successful therapies. Here, we present an integrative approach based on molecular pathways to expose new potentially actionable targets. We used gene set enrichment analysis (GSEA) to conduct an unsupervised clustering analysis to condense the gene expression data from bulk patient samples and patient-derived gliomasphere lines into new gene signatures. We identified key targets that are predicted to be differentially activated between tumors and were functionally validated in a library of gliomasphere cultures. Resultant cluster-specific gene signatures associated not only with hallmarks of cell cycle and stemness gene expression, but also with cell-type specific markers and different cellular states of GBM. Several upstream regulators, such as PIK3R1 and EBF1 were differentially enriched in cells bearing stem cell like signatures and bear further investigation. We identified the transcription factor E2F1 as a key regulator of tumor cell proliferation and self-renewal in only a subset of gliomasphere cultures predicted to be E2F1 signaling dependent. Our in vivo work also validated the functional significance of E2F1 in tumor formation capacity in the predicted samples. E2F1 inhibition also differentially sensitized E2F1-dependent gliomasphere cultures to radiation treatment. Our findings indicate that this novel approach exploring cancer pathways highlights key therapeutic vulnerabilities for targeting GBM.