Decision by division: making cortical maps

Spatial Single-cell Analysis Decodes Cortical Layer and Area Specification

The human cerebral cortex, pivotal for advanced cognitive functions, is composed of six distinct layers and dozens of functionally specialized areas1,2. The layers and areas are distinguished both molecularly, by diverse neuronal and glial cell subtypes, and structurally, through intricate spatial organization3,4. While single-cell transcriptomics studies have advanced molecular characterization of human cortical development, a critical gap exists due to the loss of spatial context during cell dissociation5,6,7,8. Here, we utilized multiplexed error-robust fluorescence in situ hybridization (MERFISH)9, augmented with deep-learning-based cell segmentation, to examine the molecular, cellular, and cytoarchitectural development of human fetal cortex with spatially resolved single-cell resolution. Our extensive spatial atlas, encompassing 16 million single cells, spans eight cortical areas across four time points in the second and third trimesters. We uncovered an early establishment of the six-layer structure, identifiable in the laminar distribution of excitatory neuronal subtypes by mid-gestation, long before the emergence of cytoarchitectural layers. Notably, while anterior-posterior gradients of neuronal subtypes were generally observed in most cortical areas, a striking exception was the sharp molecular border between primary (V1) and secondary visual cortices (V2) at gestational week 20. Here we discovered an abrupt binary shift in neuronal subtype specification at the earliest stages, challenging the notion that continuous morphogen gradients dictate mid-gestation cortical arealization6,10. Moreover, integrating single-nuclei RNA-sequencing and in situ whole transcriptomics revealed an early upregulation of synaptogenesis in V1-specific Layer 4 neurons, suggesting a role of synaptogenesis in this discrete border formation. Collectively, our findings underscore the crucial role of spatial relationships in determining the molecular specification of cortical layers and areas. This work not only provides a valuable resource for the field, but also establishes a spatially resolved single-cell analysis paradigm that paves the way for a comprehensive developmental atlas of the human brain.

The cytoarchitectonic landscape revealed by deep learning method facilitated precise positioning in mouse neocortex

Neocortex is a complex structure with different cortical sublayers and regions. However, the precise positioning of cortical regions can be challenging due to the absence of distinct landmarks without special preparation. To address this challenge, we developed a cytoarchitectonic landmark identification pipeline. The fluorescence micro-optical sectioning tomography method was employed to image the whole mouse brain stained by general fluorescent nucleotide dye. A fast 3D convolution network was subsequently utilized to segment neuronal somas in entire neocortex. By approach, the cortical cytoarchitectonic profile and the neuronal morphology were analyzed in 3D, eliminating the influence of section angle. And the distribution maps were generated that visualized the number of neurons across diverse morphological types, revealing the cytoarchitectonic landscape which characterizes the landmarks of cortical regions, especially the typical signal pattern of barrel cortex. Furthermore, the cortical regions of various ages were aligned using the generated cytoarchitectonic landmarks suggesting the structural changes of barrel cortex during the aging process. Moreover, we observed the spatiotemporally gradient distributions of spindly neurons, concentrated in the deep layer of primary visual area, with their proportion decreased over time. These findings could improve structural understanding of neocortex, paving the way for further exploration with this method.

Epigenomic tomography for probing spatially defined chromatin state in the brain

Spatially resolved epigenomic profiling is critical for understanding biology in the mammalian brain. Single-cell spatial epigenomic assays were developed recently for this purpose, but they remain costly and labor intensive for examining brain tissues across substantial dimensions and surveying a collection of brain samples. Here, we demonstrate an approach, epigenomic tomography, that maps spatial epigenomes of mouse brain at the scale of centimeters. We individually profiled neuronal and glial fractions of mouse neocortex slices with 0.5 mm thickness. Tri-methylation of histone 3 at lysine 27 (H3K27me3) or acetylation of histone 3 at lysine 27 (H3K27ac) features across these slices were grouped into clusters based on their spatial variation patterns to form epigenomic brain maps. As a proof of principle, our approach reveals striking dynamics in the frontal cortex due to kainic-acid-induced seizure, linked with transmembrane ion transporters, exocytosis of synaptic vesicles, and secretion of neurotransmitters. Epigenomic tomography provides a powerful and cost-effective tool for characterizing brain disorders based on the spatial epigenome.

Massively parallel disruption of enhancers active in human neural stem cells

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

Developmental loss of NMDA receptors results in supernumerary forebrain neurons through delayed maturation of transit-amplifying neuroblasts

Developmental neurogenesis is a tightly regulated spatiotemporal process with its dysregulation implicated in neurodevelopmental disorders. NMDA receptors are glutamate-gated ion channels that are widely expressed in the early nervous system, yet their contribution to neurogenesis is poorly understood. Notably, a variety of mutations in genes encoding NMDA receptor subunits are associated with neurodevelopmental disorders. To rigorously define the role of NMDA receptors in developmental neurogenesis, we used a mutant zebrafish line (grin1-/-) that lacks all NMDA receptors yet survives to 10 days post-fertilization, offering the opportunity to study post-embryonic neurodevelopment in the absence of NMDA receptors. Focusing on the forebrain, we find that these fish have a progressive supernumerary neuron phenotype confined to the telencephalon at the end of embryonic neurogenesis, but which extends to all forebrain regions during postembryonic neurogenesis. This enhanced neuron population does not arise directly from increased numbers or mitotic activity of radial glia cells, the principal neural stem cells. Rather, it stems from a lack of timely maturation of transit-amplifying neuroblasts into post-mitotic neurons, as indicated by a decrease in expression of the ontogenetically-expressed chloride transporter, KCC2. Pharmacological blockade with MK-801 recapitulates the grin1-/- supernumerary neuron phenotype, indicating a requirement for ionotropic signaling. Thus, NMDA receptors are required for suppression of indirect, transit amplifying cell-driven neurogenesis by promoting maturational termination of mitosis. Loss of suppression results in neuronal overpopulation that can fundamentally change brain circuitry and may be a key factor in pathogenesis of neurodevelopmental disorders caused by NMDA receptor dysfunction.

Transcriptional cartography integrates multiscale biology of the human cortex

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Continuity and change in neural plasticity through embryonic morphogenesis, fetal activity-dependent synaptogenesis, and infant memory consolidation

There is an apparent continuity in human neural development that can be traced to venerable themes of vertebrate morphogenesis that have shaped the evolution of the reptilian telencephalon (including both primitive three-layered cortex and basal ganglia) and then the subsequent evolution of the mammalian six-layered neocortex. In this theoretical analysis, we propose that an evolutionary-developmental analysis of these general morphogenetic themes can help to explain the embryonic development of the dual divisions of the limbic system that control the dorsal and ventral networks of the human neocortex. These include the archicortical (dorsal limbic) Papez circuits regulated by the hippocampus that organize spatial, contextual memory, as well as the paleocortical (ventral limbic) circuits that organize object memory. We review evidence that these dorsal and ventral limbic divisions are controlled by the differential actions of brainstem lemnothalamic and midbrain collothalamic arousal control systems, respectively, thereby traversing the vertebrate subcortical neuraxis. These dual control systems are first seen shaping the phyletic morphogenesis of the archicortical and paleocortical foundations of the forebrain in embryogenesis. They then provide dual modes of activity-dependent synaptic organization in the active (lemnothalamic) and quiet (collothalamic) stages of fetal sleep. Finally, these regulatory systems mature to form the major systems of memory consolidation of postnatal development, including the rapid eye movement (lemnothalamic) consolidation of implicit memory and social attachment in the first year, and then-in a subsequent stage-the non-REM (collothalamic) consolidation of explicit memory that is integral to the autonomy and individuation of the second year of life.

Cytoarchitectonic and connection stripes in the dysgranular insular cortex in the macaque monkey

The insula has been classically divided into broad granular, dysgranular, and agranular architectonic sectors. We previously proposed a novel partition, dividing each sector into four to seven sharply delimited architectonic areas, with the dysgranular areas being possibly further subdivided into subtle horizontal partitions or "stripes." In architectonics, discrete subparcellations are prone to subjective variability and need being supported with additional neuroanatomical methods. Here, using a secondary analysis of indirect connectional data in the rhesus macaque monkey, we examined the spatial relationship between the dysgranular architectonic stripes and tract-tracing labeling patterns produced in the insula with injections of neuronal tracers in other cortical regions. The injections consistently produced sharply delimited patches of anterograde and/or retrograde labeling, which formed stripes across consecutive coronal sections of the insula. While the overall pattern of labeling on individual coronal sections varied with the injection site, the boundaries of the patches consistently coincided with architectonic boundaries on an adjacent cyto- (Nissl) and/or myelo- (Gallyas) architectonic section. This overlap supports the existence of a fine dysgranular stripe-like partition of the primate insula, with possibly major implications for interoceptive processing in primates including humans. The modular organization of the insula could underlie a serial stream of integration from a dorsal primary interoceptive cortex toward progressively more ventral egocentric "self-agency" and allocentric "social" dysgranular processing units.

The bidirectional effects between cognitive ability and brain morphology: A life course Mendelian randomization analysis

Introduction

Little is understood about the dynamic interplay between brain morphology and cognitive ability across the life course. Additionally, most existing research has focused on global morphology measures such as estimated total intracranial volume, mean thickness, and total surface area.

Methods

Mendelian randomization was used to estimate the bidirectional effects between cognitive ability, global and regional measures of cortical thickness and surface area, estimated total intracranial volume, total white matter, and the volume of subcortical structures (N=37,864). Analyses were stratified for developmental periods (childhood, early adulthood, mid-to-late adulthood; age range: 8-81 years).

Results

The earliest effects were observed in childhood and early adulthood in the frontoparietal lobes. A bidirectional relationship was identified between higher cognitive ability, larger estimated total intracranial volume (childhood, mid-to-late adulthood) and total surface area (all life stages). A thicker posterior cingulate cortex and a larger surface area in the caudal middle frontal cortex and temporal pole were associated with greater cognitive ability. Contrary, a thicker temporal pole was associated with lower cognitive ability.

Discussion

Stable effects of cognitive ability on brain morphology across the life course suggests that childhood is potentially an important window for intervention.

Differential effects of generalized anxiety and separation anxiety on brain structural development during adolescence

Excessive anxiety is highly prevalent during childhood and adolescence, with detrimental effects on somatic and mental health, and quality of life. Although structural abnormalities in the brain have been found in people with anxiety disorders, whether anxiety affects the brain development of children and adolescents remains unknown. Here, we applied a multivariate approach to two single-site MRI datasets consisting of 733 and 775 participants aged 5-18 years. Using linear support vector regression and cross-validation, brain age is estimated by predicting the chronological age from the features that combine cortical thickness and surface area of 68 brain regions. We found that gray matter can predict the chronological age of children and adolescents with a low mean absolute error. Compared to specific brain network, the whole structural brain measures predicted brain age better. Importantly, adolescents with higher generalized anxiety and those with lower separation anxiety showed lower brain age, indicating a slow development of brain structures. The relationship between anxiety and brain age of youths could also be found in parent-reported separation anxiety. The findings highlight differential effects of different anxiety types on brain structural development and suggest that different types of anxiety during childhood and adolescence should be treated differently.

Schizophrenia Polygenic Risk During Typical Development Reflects Multiscale Cortical Organization

Background

Schizophrenia is widely recognized as a neurodevelopmental disorder. Abnormal cortical development in otherwise typically developing children and adolescents may be revealed using polygenic risk scores for schizophrenia (PRS-SCZ).

Methods

We assessed PRS-SCZ and cortical morphometry in typically developing children and adolescents (3-21 years, 46.8% female) using whole-genome genotyping and T1-weighted magnetic resonance imaging (n = 390) from the PING (Pediatric Imaging, Neurocognition, and Genetics) cohort. We contextualized the findings using 1) age-matched transcriptomics, 2) histologically defined cytoarchitectural types and functionally defined networks, and 3) case-control differences of schizophrenia and other major psychiatric disorders derived from meta-analytic data of 6 ENIGMA (Enhancing Neuro Imaging Genetics through Meta Analysis) working groups, including a total of 12,876 patients and 15,670 control participants.

Results

Higher PRS-SCZ was associated with greater cortical thickness, which was most prominent in areas with heightened gene expression of dendrites and synapses. PRS-SCZ-related increases in vertexwise cortical thickness were mainly distributed in association cortical areas, particularly the ventral attention network, while relatively sparing koniocortical type cortex (i.e., primary sensory areas). The large-scale pattern of cortical thickness increases related to PRS-SCZ mirrored the pattern of cortical thinning in schizophrenia and mood-related psychiatric disorders derived from the ENIGMA consortium. Age group models illustrate a possible trajectory from PRS-SCZ-associated cortical thickness increases in early childhood toward thinning in late adolescence, with the latter resembling the adult brain phenotype of schizophrenia.

Conclusions

Collectively, combining imaging genetics with multiscale mapping, our work provides novel insight into how genetic risk for schizophrenia affects the cortex early in life.

Loss of NMDA receptor function during development results in decreased KCC2 expression and increased neurons in the zebrafish forebrain

Developmental neurogenesis is a tightly regulated spatiotemporal process with its dysregulation implicated in neurodevelopmental disorders. NMDA receptors are glutamate-gated ion channels that are widely expressed in the early nervous system, yet their contribution to neurogenesis is poorly understood. Notably, a variety of mutations in genes encoding NMDA receptor subunits are associated with neurodevelopmental disorders. To rigorously define the role of NMDA receptors in developmental neurogenesis, we used a mutant zebrafish line ( grin1 -/- ) that lacks all NMDA receptors yet survives to 10 days post-fertilization, offering the opportunity to study post-embryonic neurodevelopment in the absence of NMDA receptors. Focusing on the forebrain, we find that these fish have a progressive supernumerary neuron phenotype confined to the telencephalon at the end of embryonic neurogenesis, but which extends to all forebrain regions during postembryonic neurogenesis. This enhanced neuron population does not arise directly from increased numbers or mitotic activity of radial glia cells, the principal neural stem cells. Rather, it stems from a lack of timely maturation of transit-amplifying neuroblasts into post-mitotic neurons, as indicated by a decrease in expression of the ontogenetically-expressed chloride transporter, KCC2. Pharmacological blockade with MK-801 recapitulates the grin1 -/- supernumerary neuron phenotype, indicating a requirement for ionotropic signaling. Thus, NMDA receptors are required for suppression of indirect, transit amplifying cell-driven neurogenesis by promoting maturational termination of mitosis. Loss of suppression results in neuronal overpopulation that can fundamentally change brain circuitry and may be a key factor in pathogenesis of neurodevelopmental disorders caused by NMDA receptor dysfunction.

TempShift Reveals the Sequential Development of Human Neocortex and Skewed Developmental Timing of Down Syndrome Brains

Development is a complex process involving precise regulation. Developmental regulation may vary in tissues and individuals, and is often altered in disorders. Currently, the regulation of developmental timing across neocortical areas and developmental changes in Down syndrome (DS) brains remain unclear. The changes in regulation are often accompanied by changes in the gene expression trajectories, which can be divided into two scenarios: (1) changes of gene expression trajectory shape that reflect changes in cell type composition or altered molecular machinery; (2) temporal shift of gene expression trajectories that indicate different regulation of developmental timing. Therefore, we developed an R package TempShift to separates these two scenarios and demonstrated that TempShift can distinguish temporal shift from different shape (DiffShape) of expression trajectories, and can accurately estimate the time difference between multiple trajectories. We applied TempShift to identify sequential gene expression across 11 neocortical areas, which suggested sequential occurrence of synapse formation and axon guidance, as well as reconstructed interneuron migration pathways within neocortex. Comparison between healthy and DS brains revealed increased microglia, shortened neuronal migration process, and delayed synaptogenesis and myelination in DS. These applications also demonstrate the potential of TempShift in understanding gene expression temporal dynamics during different biological processes.

Two sides of the same coin: distinct neuroanatomical patterns predict crystallized and fluid intelligence in adults

Background

Crystallized intelligence (Gc) and fluid intelligence (Gf) are regarded as distinct intelligence components that statistically correlate with each other. However, the distinct neuroanatomical signatures of Gc and Gf in adults remain contentious.

Methods

Machine learning cross-validated elastic net regression models were performed on the Human Connectome Project Young Adult dataset (N = 1089) to characterize the neuroanatomical patterns of structural magnetic resonance imaging variables that are associated with Gc and Gf. The observed relationships were further examined by linear mixed-effects models. Finally, intraclass correlations were computed to examine the similarity of the neuroanatomical correlates between Gc and Gf.

Results

The results revealed distinct multi-region neuroanatomical patterns predicted Gc and Gf, respectively, which were robust in a held-out test set (R² = 2.40, 1.97%, respectively). The relationship of these regions with Gc and Gf was further supported by the univariate linear mixed effects models. Besides that, Gc and Gf displayed poor neuroanatomical similarity.

Conclusion

These findings provided evidence that distinct machine learning-derived neuroanatomical patterns could predict Gc and Gf in healthy adults, highlighting differential neuroanatomical signatures of different aspects of intelligence.

CellCover Defines Conserved Cell Types and Temporal Progression in scRNA-seq Data across Mammalian Neocortical Development

Accurate identification of cell classes across the tissues of living organisms is central in the analysis of growing atlases of single-cell RNA sequencing (scRNA-seq) data across biomedicine. Such analyses are often based on the existence of highly discriminating "marker genes" for specific cell classes which enables a deeper functional understanding of these classes as well as their identification in new, related datasets. Currently, marker genes are defined by methods that serially assess the level of differential expression (DE) of individual genes across landscapes of diverse cells. This serial approach has been extremely useful, but is limited because it ignores possible redundancy or complementarity across genes, that can only be captured by analyzing several genes at the same time. We wish to identify discriminating panels of genes. To efficiently explore the vast space of possible marker panels, leverage the large number of cells often sequenced, and overcome zero-inflation in scRNA-seq data, we propose viewing panel selection as a variation of the "minimal set-covering problem" in combinatorial optimization which can be solved with integer programming. In this formulation, the covering elements are genes, and the objects to be covered are cells of a particular class, where a cell is covered by a gene if that gene is expressed in that cell. Our method, CellCover, identifies a panel of marker genes in scRNA-seq data that covers one class of cells within a population. We apply this method to generate covering marker gene panels which characterize cells of the developing mouse neocortex as postmitotic neurons are generated from neural progenitor cells (NPCs). We show that CellCover captures cell class-specific signals distinct from those defined by DE methods and that CellCover's compact gene panels can be expanded to explore cell type specific function.Transfer learning experiments exploring these covering panels across in vivo mouse, primate, and human scRNA-seq datasets demonstrate that CellCover identifies markers of conserved cell classes in neurogenesis, as well as markers of temporal progression in the molecular identity of these cell types across development of the mammalian neocortex. The gene covering panels we identify across cell types and developmental time can be freely explored in visualizations across all the public data we use in this report at with NeMo Analytics [1] through https://nemoanalytics.org/p?l=CellCover . The code for CellCover is written in R and the Gurobi R interface and is available at [2].

A human dynein heavy chain mutation impacts cortical progenitor cells causing developmental defects, reduced brain size and altered brain architecture

Dynein heavy chain (DYNC1H1) mutations can either lead to severe cerebral cortical malformations, or alternatively may be associated with the development of spinal muscular atrophy with lower extremity predominance (SMA-LED). To assess the origin of such differences, we studied a new Dync1h1 knock-in mouse carrying the cortical malformation p.Lys3334Asn mutation. Comparing with an existing neurodegenerative Dync1h1 mutant (Legs at odd angles, Loa, +/p.Phe580Tyr), we assessed Dync1h1's roles in cortical progenitor and especially radial glia functions during embryogenesis, and assessed neuronal differentiation. p.Lys3334Asn /+ mice exhibit reduced brain and body size. Embryonic brains show increased and disorganized radial glia: interkinetic nuclear migration occurs in mutants, however there are increased basally positioned cells and abventricular mitoses. The ventricular boundary is disorganized potentially contributing to progenitor mislocalization and death. Morphologies of mitochondria and Golgi apparatus are perturbed in vitro, with different effects also in Loa mice. Perturbations of neuronal migration and layering are also observed in p.Lys3334Asn /+ mutants. Overall, we identify specific developmental effects due to a severe cortical malformation mutation in Dync1h1, highlighting the differences with a mutation known instead to primarily affect motor function.

Larger cerebral cortex is genetically correlated with greater frontal area and dorsal thickness

Human cortical expansion has occurred non-uniformly across the brain. We assessed the genetic architecture of cortical global expansion and regionalization by comparing two sets of genome-wide association studies of 24 cortical regions with and without adjustment for global measures (i.e., total surface area, mean cortical thickness) using a genetically informed parcellation in 32,488 adults. We found 393 and 756 significant loci with and without adjusting for globals, respectively, where 8% and 45% loci were associated with more than one region. Results from analyses without adjustment for globals recovered loci associated with global measures. Genetic factors that contribute to total surface area of the cortex particularly expand anterior/frontal regions, whereas those contributing to thicker cortex predominantly increase dorsal/frontal-parietal thickness. Interactome-based analyses revealed significant genetic overlap of global and dorsolateral prefrontal modules, enriched for neurodevelopmental and immune system pathways. Consideration of global measures is important in understanding the genetic variants underlying cortical morphology.

Divergent topographic projection of cerebral cortical areas to overlapping cerebellar lobules through distinct regions of the pontine nuclei

The massive axonal projection from the cerebrum to the cerebellum through the pontine nuclei supports the cerebrocerebellar coordination of motor and nonmotor functions. However, the cerebrum and cerebellum have distinct patterns of functional localization in their cortices. We addressed this issue by bidirectional neuronal tracing from 22 various locations of the pontine nuclei in the mouse in a comprehensive manner. Cluster analyses of the distribution patterns of labeled cortical pyramidal cells and cerebellar mossy fiber terminals classified all cases into six groups located in six different subareas of the pontine nuclei. The lateral (insular), mediorostral (cingulate and prefrontal), and caudal (visual and auditory) cortical areas of the cerebrum projected to the medial, rostral, and lateral subareas of the pontine nuclei, respectively. These pontine subareas then projected mainly to the crus I, central vermis, and paraflocculus divergently. The central (motor and somatosensory) cortical areas projected to the centrorostral, centrocaudal and caudal subareas of the pontine nuclei, which then projected mainly to the rostral and caudal lobules with a somatotopic arrangement. The results indicate a new pontine nuclei-centric view of the corticopontocerebellar projection: the generally parallel corticopontine projection to pontine nuclei subareas is relayed to the highly divergent pontocerebellar projection terminating in overlapping specific lobules of the cerebellum. Consequently, the mode of the pontine nuclei relay underlies the cerebellar functional organization.

BEAM: a combinatorial recombinase toolbox for binary gene expression and mosaic analysis

Genetic mosaic analysis, in which mutant cells reside intermingled with wild-type cells, is a powerful experimental approach, but has not been widely used in mice because existing genome-based strategies require complicated and protracted breeding schemes. We have developed an alternative approach termed BEAM (for Binary Expression Aleatory Mosaic) that relies on sparse recombinase activation to generate two genetically distinct, non-overlapping populations of cells for comparative analysis. Following delivery of DNA constructs by transfection or viral transduction, combinatorial recombinase activity generates two distinct populations of cells labeled with either green or red fluorescent protein. Any gene of interest can be mis-expressed or deleted in one population for comparison with intermingled control cells. We have extensively optimized and characterized this system both in vitro and in vivo , and demonstrate its power for investigating cell autonomy, identifying temporally or spatially aberrant phenotypes, revealing changes in cell proliferation or death, and controlling for procedural variability.

De novo human brain enhancers created by single-nucleotide mutations

Advanced human cognition is attributed to increased neocortex size and complexity, but the underlying evolutionary and regulatory mechanisms are largely unknown. Using human and macaque embryonic neocortical H3K27ac data coupled with a deep learning model of enhancers, we identified ~4000 enhancer gains in humans, which, per our model, can often be attributed to single-nucleotide essential mutations. Our analyses suggest that functional gains in embryonic brain development are associated with de novo enhancers whose putative target genes exhibit increased expression in progenitor cells and interneurons and partake in critical neural developmental processes. Essential mutations alter enhancer activity through altered binding of key transcription factors (TFs) of embryonic neocortex, including ISL1, POU3F2, PITX1/2, and several SOX TFs, and are associated with central nervous system disorders. Overall, our results suggest that essential mutations lead to gain of embryonic neocortex enhancers, which orchestrate expression of genes involved in critical developmental processes associated with human cognition.

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.

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.

Genetic map of regional sulcal morphology in the human brain from UK biobank data

Genetic associations with macroscopic brain structure can provide insights into brain function and disease. However, specific associations with measures of local brain folding are largely under-explored. Here, we conducted large-scale genome- and exome-wide associations of regional cortical sulcal measures derived from magnetic resonance imaging scans of 40,169 individuals in UK Biobank. We discovered 388 regional brain folding associations across 77 genetic loci, with genes in associated loci enriched for expression in the cerebral cortex, neuronal development processes, and differential regulation during early brain development. We integrated brain eQTLs to refine genes for various loci, implicated several genes involved in neurodevelopmental disorders, and highlighted global genetic correlations with neuropsychiatric phenotypes. We provide an interactive 3D visualisation of our summary associations, emphasising added resolution of regional analyses. Our results offer new insights into the genetic architecture of brain folding and provide a resource for future studies of sulcal morphology in health and disease.

Teleost Fish and Organoids: Alternative Windows Into the Development of Healthy and Diseased Brains

The variety in the display of animals’ cognition, emotions, and behaviors, typical of humans, has its roots within the anterior-most part of the brain: the forebrain, giving rise to the neocortex in mammals. Our understanding of cellular and molecular events instructing the development of this domain and its multiple adaptations within the vertebrate lineage has progressed in the last decade. Expanding and detailing the available knowledge on regionalization, progenitors’ behavior and functional sophistication of the forebrain derivatives is also key to generating informative models to improve our characterization of heterogeneous and mechanistically unexplored cortical malformations. Classical and emerging mammalian models are irreplaceable to accurately elucidate mechanisms of stem cells expansion and impairments of cortex development. Nevertheless, alternative systems, allowing a considerable reduction of the burden associated with animal experimentation, are gaining popularity to dissect basic strategies of neural stem cells biology and morphogenesis in health and disease and to speed up preclinical drug testing. Teleost vertebrates such as zebrafish, showing conserved core programs of forebrain development, together with patients-derived in vitro 2D and 3D models, recapitulating more accurately human neurogenesis, are now accepted within translational workflows spanning from genetic analysis to functional investigation. Here, we review the current knowledge of common and divergent mechanisms shaping the forebrain in vertebrates, and causing cortical malformations in humans. We next address the utility, benefits and limitations of whole-brain/organism-based fish models or neuronal ensembles in vitro for translational research to unravel key genes and pathological mechanisms involved in neurodevelopmental diseases.

Thalamocortical axons control the cytoarchitecture of neocortical layers by area-specific supply of VGF

Neuronal abundance and thickness of each cortical layer are specific to each area, but how this fundamental feature arises during development remains poorly understood. While some of area-specific features are controlled by intrinsic cues such as morphogens and transcription factors, the exact influence and mechanisms of action by cues extrinsic to the cortex, in particular the thalamic axons, have not been fully established. Here, we identify a thalamus-derived factor, VGF, which is indispensable for thalamocortical axons to maintain the proper amount of layer 4 neurons in the mouse sensory cortices. This process is prerequisite for further maturation of the primary somatosensory area, such as barrel field formation instructed by a neuronal activity-dependent mechanism. Our results provide an actual case in which highly site-specific axon projection confers further regional complexity upon the target field through locally secreting signaling molecules from axon terminals.

The ABCD Study: Brain Heterogeneity in Intelligence During a Neurodevelopmental Transition Stage

A complex curvilinear relationship exists between intelligence and age during the neurodevelopment of cortical thickness. To parse out a more fine-grained relationship between intelligence and cortical thickness and surface area, we used a large-scale data set focusing on a critical transition juncture in neurodevelopment in preadolescence. Cortical thickness was derived from T1-weighted structural magnetic resonance images of a large sample of 9- and 11-year-old children from the Adolescent Brain Cognitive Development study. The NIH Toolbox Cognition Battery composite scores, which included fluid, crystallized, and total scores, were used to assess intelligence. Using a double generalized linear model, we assessed the independent association between the mean and dispersion of cortical thickness/surface area and intelligence. Higher intelligence in preadolescents was associated with higher mean cortical thickness in orbitofrontal and primary sensory cortices but with lower thickness in the dorsolateral and medial prefrontal cortex and particularly in the rostral anterior cingulate. The rostral anterior cingulate findings were particularly evident across all subscales of intelligence. Higher intelligence was also associated with greater interindividual similarity in the rostral cingulate. Intelligence during this key transition juncture in preadolescence appears to reflect a dissociation between the cortical development of basic cognitive processes and higher-order executive and motivational processes.

Genetic Influences on Longitudinal Trajectories of Cortical Thickness and Surface Area during the First 2 Years of Life

Genetic influences on cortical thickness (CT) and surface area (SA) are known to vary across the life span. Little is known about the extent to which genetic factors influence CT and SA in infancy and toddlerhood. We performed the first longitudinal assessment of genetic influences on variation in CT and SA in 501 twins who were aged 0-2 years. We observed substantial additive genetic influences on both average CT (0.48 in neonates, 0.37 in 1-year-olds, and 0.44 in 2-year-olds) and total SA (0.59 in neonates, 0.74 in 1-year-olds, and 0.73 in 2-year-olds). In addition, we found strong heritability of the change in average CT (0.49) from neonates to 1-year-olds, but not from 1- to 2-year-olds. Moreover, we found strong genetic correlations for average CT (rG = 0.92) between 1- and 2-year-olds and strong genetic correlations for total SA across all timepoints (rG = 0.96 between neonates and 1-year-olds, rG = 1 between 1- and 2-year-olds). In addition, we found CT and SA are strongly genetic correlated at birth, but weaken over time. Overall, results suggest a dynamic genetic relationship between CT and SA during first 2 years of life and provide novel insights into how genetic influences shape the cortical structure during early brain development.

Structural and Functional Aspects of the Neurodevelopmental Gene NR2F1: From Animal Models to Human Pathology

The assembly and maturation of the mammalian brain result from an intricate cascade of highly coordinated developmental events, such as cell proliferation, migration, and differentiation. Any impairment of this delicate multi-factorial process can lead to complex neurodevelopmental diseases, sharing common pathogenic mechanisms and molecular pathways resulting in multiple clinical signs. A recently described monogenic neurodevelopmental syndrome named Bosch-Boonstra-Schaaf Optic Atrophy Syndrome (BBSOAS) is caused by NR2F1 haploinsufficiency. The NR2F1 gene, coding for a transcriptional regulator belonging to the steroid/thyroid hormone receptor superfamily, is known to play key roles in several brain developmental processes, from proliferation and differentiation of neural progenitors to migration and identity acquisition of neocortical neurons. In a clinical context, the disruption of these cellular processes could underlie the pathogenesis of several symptoms affecting BBSOAS patients, such as intellectual disability, visual impairment, epilepsy, and autistic traits. In this review, we will introduce NR2F1 protein structure, molecular functioning, and expression profile in the developing mouse brain. Then, we will focus on Nr2f1 several functions during cortical development, from neocortical area and cell-type specification to maturation of network activity, hippocampal development governing learning behaviors, assembly of the visual system, and finally establishment of cortico-spinal descending tracts regulating motor execution. Whenever possible, we will link experimental findings in animal or cellular models to corresponding features of the human pathology. Finally, we will highlight some of the unresolved questions on the diverse functions played by Nr2f1 during brain development, in order to propose future research directions. All in all, we believe that understanding BBSOAS mechanisms will contribute to further unveiling pathophysiological mechanisms shared by several neurodevelopmental disorders and eventually lead to effective treatments.

Primary Cilia and Centrosomes in Neocortex Development

During mammalian brain development, neural stem and progenitor cells generate the neurons for the six-layered neocortex. The proliferative capacity of the different types of progenitor cells within the germinal zones of the developing neocortex is a major determinant for the number of neurons generated. Furthermore, the various modes of progenitor cell divisions, for which the orientation of the mitotic spindle of progenitor cells has a pivotal role, are a key parameter to ensure the appropriate size and proper cytoarchitecture of the neocortex. Here, we review the roles of primary cilia and centrosomes of progenitor cells in these processes during neocortical development. We specifically focus on the apical progenitor cells in the ventricular zone. In particular, we address the alternating, dual role of the mother centriole (i) as a component of one of the spindle poles during mitosis, and (ii) as the basal body of the primary cilium in interphase, which is pivotal for the fate of apical progenitor cells and their proliferative capacity. We also discuss the interactions of these organelles with the microtubule and actin cytoskeleton, and with junctional complexes. Centriolar appendages have a specific role in this interaction with the cell cortex and the plasma membrane. Another topic of this review is the specific molecular composition of the ciliary membrane and the membrane vesicle traffic to the primary cilium of apical progenitors, which underlie the ciliary signaling during neocortical development; this signaling itself, however, is not covered in depth here. We also discuss the recently emerging evidence regarding the composition and roles of primary cilia and centrosomes in basal progenitors, a class of progenitors thought to be of particular importance for neocortex expansion in development and evolution. While the tight interplay between primary cilia and centrosomes makes it difficult to allocate independent roles to either organelle, mutations in genes encoding ciliary and/or centrosome proteins indicate that both are necessary for the formation of a properly sized and functioning neocortex during development. Human neocortical malformations, like microcephaly, underpin the importance of primary cilia/centrosome-related processes in neocortical development and provide fundamental insight into the underlying mechanisms involved.

Enhanced phosphorylation of S6 protein in mouse cortical layer V and subplate neurons

The mammalian neocortex is composed of six major layers of neurons. Each group of neurons in the cortical layers has distinct characteristics based on the expression of specific genes and connectivity patterns of neural circuits. Neuronal subtype transition and regional identity acquisition are established by temporal cues and interaction between several transcription factors during neurogenesis. The impairment of cortical lamination or neural circuits results in a wide range of neurodevelopmental disorders such as autism, schizophrenia, and certain forms of childhood epilepsy. Despite continuous efforts to classify neurons with the aid of genetic and epigenetic analyses, the neuron-specific properties associated with post-transcriptional modification remain unclear. In the present study, the distribution of phosphorylated S6-positive layers across the neocortex was examined using several layer markers. The development of pS6 S235/236 layers in layer V and the subplate was spatiotemporally regulated in the mouse brain. In addition, enhanced phosphorylation of ribosomal protein S6 in Ctip2-positive layer V neurons in vivo was sustained under in-vitro conditions using a culture of primary cortical neurons.

Distinct Regionalization Patterns of Cortical Morphology are Associated with Cognitive Performance Across Different Domains

Cognitive performance in children is predictive of academic and social outcomes; therefore, understanding neurobiological mechanisms underlying individual differences in cognition during development may be important for improving quality of life. The belief that a single, psychological construct underlies many cognitive processes is pervasive throughout society. However, it is unclear if there is a consistent neural substrate underlying many cognitive processes. Here, we show that a distributed configuration of cortical surface area and apparent thickness, when controlling for global imaging measures, is differentially associated with cognitive performance on different types of tasks in a large sample (N = 10 145) of 9-11-year-old children from the Adolescent Brain and Cognitive DevelopmentSM (ABCD) study. The minimal overlap in these regionalization patterns of association has implications for competing theories about developing intellectual functions. Surprisingly, not controlling for sociodemographic factors increased the similarity between these regionalization patterns. This highlights the importance of understanding the shared variance between sociodemographic factors, cognition and brain structure, particularly with a population-based sample such as ABCD.

Disruption of Transient SERT Expression in Thalamic Glutamatergic Neurons Alters Trajectory of Postnatal Interneuron Development in the Mouse Cortex

In mice, terminal differentiation of subpopulations of interneurons occurs in late postnatal stages, paralleling the emergence of the adult cortical architecture. Here, we investigated the effects of altered initial cortical architecture on later interneuron development. We identified that a class of somatostatin (SOM)-expressing GABAergic interneurons undergoes terminal differentiation between 2nd and 3rd postnatal week in the mouse somatosensory barrel cortex and upregulates Reelin expression during neurite outgrowth. Our previous work demonstrated that transient expression (E15-P10) of serotonin uptake transporter (SERT) in thalamocortical projection neurons regulates barrel elaboration during cortical map establishment. We show here that in thalamic neuron SERT knockout mice, these SOM-expressing interneurons develop at the right time, reach correct positions and express correct neurochemical markers, but only 70% of the neurons remain in the adult barrel cortex. Moreover, those neurons that remain display altered dendritic patterning. Our data indicate that a precise architecture at the cortical destination is not essential for specifying late-developing interneuron identities, their cortical deposition, and spatial organization, but dictates their number and dendritic structure ultimately integrated into the cortex. Our study illuminates how disruption of temporal-specific SERT function and related key regulators during cortical map establishment can alter interneuron development trajectory that persists to adult central nervous system.

Cellular and Widefield Imaging of Sound Frequency Organization in Primary and Higher Order Fields of the Mouse Auditory Cortex

The mouse auditory cortex (ACtx) contains two core fields-primary auditory cortex (A1) and anterior auditory field (AAF)-arranged in a mirror reversal tonotopic gradient. The best frequency (BF) organization and naming scheme for additional higher order fields remain a matter of debate, as does the correspondence between smoothly varying global tonotopy and heterogeneity in local cellular tuning. Here, we performed chronic widefield and two-photon calcium imaging from the ACtx of awake Thy1-GCaMP6s reporter mice. Data-driven parcellation of widefield maps identified five fields, including a previously unidentified area at the ventral posterior extreme of the ACtx (VPAF) and a tonotopically organized suprarhinal auditory field (SRAF) that extended laterally as far as ectorhinal cortex. Widefield maps were stable over time, where single pixel BFs fluctuated by less than 0.5 octaves throughout a 1-month imaging period. After accounting for neuropil signal and frequency tuning strength, BF organization in neighboring layer 2/3 neurons was intermediate to the heterogeneous salt and pepper organization and the highly precise local organization that have each been described in prior studies. Multiscale imaging data suggest there is no ultrasonic field or secondary auditory cortex in the mouse. Instead, VPAF and a dorsal posterior (DP) field emerged as the strongest candidates for higher order auditory areas.

Heritable functional architecture in human visual cortex

•

We analyzed retinotopic maps from monozygotic and dizygotic twin pairs.

•

Visual field maps in V1-V3 are more similar in monozygotic twins.

•

Heritability is greater in V1 and V3 for polar angle and population receptive field sizes.

•

Eccentricity maps show lesser degree of heritability.

•

Further evidence for link between cortical morphology and topology of retinotopic maps.

How much of the functional organization of our visual system is inherited? Here we tested the heritability of retinotopic maps in human visual cortex using functional magnetic resonance imaging. We demonstrate that retinotopic organization shows a closer correspondence in monozygotic (MZ) compared to dizygotic (DZ) twin pairs, suggesting a partial genetic determination. Using population receptive field (pRF) analysis to examine the preferred spatial location and selectivity of these neuronal populations, we estimate a heritability around 10–20% for polar angle preferences and spatial selectivity, as quantified by pRF size, in extrastriate areas V2 and V3. Our findings are consistent with heritability in both the macroscopic arrangement of visual regions and stimulus tuning properties of visual cortex. This could constitute a neural substrate for variations in a range of perceptual effects, which themselves have been found to be at least partially genetically determined. These findings also add convergent evidence for the hypothesis that functional map topology is linked with cortical morphology.

Creative Destruction: A Basic Computational Model of Cortical Layer Formation

One of the most characteristic properties of many vertebrate neural systems is the layered organization of different cell types. This cytoarchitecture exists in the cortex, the retina, the hippocampus, and many other parts of the central nervous system. The developmental mechanisms of neural layer formation have been subject to substantial experimental efforts. Here, we provide a general computational model for cortical layer formation in 3D physical space. We show that this multiscale, agent-based model, comprising two distinct stages of apoptosis, can account for the wide range of neuronal numbers encountered in different cortical areas and species. Our results demonstrate the phenotypic richness of a basic state diagram structure. Importantly, apoptosis allows for changing the thickness of one layer without automatically affecting other layers. Therefore, apoptosis increases the flexibility for evolutionary change in layer architecture. Notably, slightly changed gene regulatory dynamics recapitulate the characteristic properties observed in neurodevelopmental diseases. Overall, we propose a novel computational model using gene-type rules, exhibiting many characteristics of normal and pathological cortical development.

Early Forebrain Neurons and Scaffold Fibers in Human Embryos

Neural progenitor proliferation, neuronal migration, areal organization, and pioneer axon wiring are critical events during early forebrain development, yet remain incompletely understood, especially in human. Here, we studied forebrain development in human embryos aged 5 to 8 postconceptional weeks (WPC5-8), stages that correspond to the neuroepithelium/early marginal zone (WPC5), telencephalic preplate (WPC6 & 7), and incipient cortical plate (WPC8). We show that early telencephalic neurons are formed at the neuroepithelial stage; the most precocious ones originate from local telencephalic neuroepithelium and possibly from the olfactory placode. At the preplate stage, forebrain organization is quite similar in human and mouse in terms of areal organization and of differentiation of Cajal-Retzius cells, pioneer neurons, and axons. Like in mice, axons from pioneer neurons in prethalamus, ventral telencephalon, and cortical preplate cross the diencephalon-telencephalon junction and the pallial-subpallial boundary, forming scaffolds that could guide thalamic and cortical axons at later stages. In accord with this model, at the early cortical plate stage, corticofugal axons run in ventral telencephalon in close contact with scaffold neurons, which express CELSR3 and FZD3, two molecules that regulates formation of similar scaffolds in mice.

Neocortical Anatomy in the South American Plains Vizcacha, Lagostomus maximus, Reveals Different Strategies in Encephalic Development among Hystricomorpha and Myomorpha Rodents

Depending on the presence or absence of sulci and convolutions, the brains of mammals are classified as gyrencephalic or lissencephalic. We analyzed the encephalic anatomy of the hystricomorph rodent Lagostomus maximus in comparison with other evolutionarily related species. The encephalization quotient (EQ), gyrencephaly index (GI), and minimum cortical thickness (MCT) were calculated for the plains vizcacha as well as for other myomorph and hystricomorph rodents. The vizcacha showed a gyrencephalic brain with a sagittal longitudinal fissure that divides both hemispheres, and 3 pairs of sulci with bilateral symmetry; that is, lateral-rostral, intraparietal, and transverse sulci. The EQ had one of the lowest values among Hystricomorpha, while GI was one of the highest. Besides, the MCT was close to the mean value for the suborder. The comparison of EQ, GI, and MCT values between hystricomorph and myomorph species allowed the detection of significant variations. Both EQ and GI showed a significant increase in Hystricomorpha compared to Myomorpha, whereas a Pearson's analysis between EQ and GI depicted an inverse correlation pattern for Hystricomorpha. Furthermore, the ratio between MCT and GI also showed a negative correlation for Hystricomorpha and Myomorpha. Our phylogenetic analyses showed that Hystricomorpha and Myomorpha do not differ in their allometric patterning between the brain and body mass, GI and brain mass, and MCT and GI. In conclusion, gyrencephalic neuroanatomy in the vizcacha could have developed from the balance between the brain size, the presence of invaginations, and the cortical thickness, which resulted in a mixed encephalization strategy for the species. Gyrencephaly in the vizcacha, as well as in other Hystricomorpha, advocates in favor of the proposal that in the more recently evolved Myomorpha lissencephaly would have arisen from a phenotype reversal process.

Enrichment of Disease-Associated Genes in Cortical Areas Defined by Transcriptome-Based Parcellation

BACKGROUND

Parcellation of the cerebral cortex serves the investigation of the emergence of uniquely human brain functions and disorders. Transcriptome data enable the characterization of the molecular properties of cortical areas in unprecedented detail. Previously, we predicted the expression of 18,686 genes in the entire human brain based on microarray data. Here, we employed these data to parcellate the cortex and study the regional enrichment of disease-associated genes.

METHODS

We performed agglomerative hierarchical clustering based on normalized transcriptome data to delineate areas with distinct gene expression profiles. Subsequently, we tested these profiles for the enrichment of gene sets associated with brain disorders by genome-wide association studies and expert-curated databases using gene set enrichment analysis.

RESULTS

Transcriptome-based parcellation identified borders in line with major anatomical landmarks and the functional differentiation of primary motor, somatosensory, visual, and auditory areas. Gene set enrichment analysis based on curated databases suggested new roles of specific areas in psychiatric and neurological disorders while reproducing well-established links for movement and neurodegenerative disorders, for example, amyotrophic lateral sclerosis (motor cortex) and Alzheimer's disease (entorhinal cortex). Meanwhile, gene sets derived from genome-wide association studies on psychiatric disorders exhibited similar enrichment patterns driven by pleiotropic genes expressed in the posterior fusiform gyrus and inferior parietal lobule.

CONCLUSIONS

The identified enrichment patterns suggest the vulnerability of specific cortical areas to various influences that might alter the risk of developing one or several brain disorders. For several diseases, specific genes were highlighted, which could lead to the discovery of novel disease mechanisms and urgently needed treatments.

Cannabinoid Type 1 Receptor is Undetectable in Rodent and Primate Cerebral Neural Stem Cells but Participates in Radial Neuronal Migration

Cannabinoid type 1 receptor (CB₁R) is expressed and participates in several aspects of cerebral cortex embryonic development as demonstrated with whole-transcriptome mRNA sequencing and other contemporary methods. However, the cellular location of CB₁R, which helps to specify molecular mechanisms, remains to be documented. Using three-dimensional (3D) electron microscopic reconstruction, we examined CB₁R immunolabeling in proliferating neural stem cells (NSCs) and migrating neurons in the embryonic mouse (Mus musculus) and rhesus macaque (Macaca mulatta) cerebral cortex. We found that the mitotic and postmitotic ventricular and subventricular zone (VZ and SVZ) cells are immunonegative in both species while radially migrating neurons in the intermediate zone (IZ) and cortical plate (CP) contain CB₁R-positive intracellular vesicles. CB₁R immunolabeling was more numerous and more extensive in monkeys compared to mice. In CB₁R-knock out mice, projection neurons in the IZ show migration abnormalities such as an increased number of lateral processes. Thus, in radially migrating neurons CB₁R provides a molecular substrate for the regulation of cell movement. Undetectable level of CB₁R in VZ/SVZ cells indicates that previously suggested direct CB₁R-transmitted regulation of cellular proliferation and fate determination demands rigorous re-examination. More abundant CB₁R expression in monkey compared to mouse suggests that therapeutic or recreational cannabis use may be more distressing for immature primate neurons than inferred from experiments with rodents.

Continuity and Discontinuity in Human Cortical Development and Change From Embryonic Stages to Old Age

The human cerebral cortex is highly regionalized, and this feature emerges from morphometric gradients in the cerebral vesicles during embryonic development. We tested if this principle of regionalization could be traced from the embryonic development to the human life span. Data-driven fuzzy clustering was used to identify regions of coordinated longitudinal development of cortical surface area (SA) and thickness (CT) (n = 301, 4-12 years). The principal divide for the developmental SA clusters extended from the inferior-posterior to the superior-anterior cortex, corresponding to the major embryonic morphometric anterior-posterior (AP) gradient. Embryonic factors showing a clear AP gradient were identified, and we found significant differences in gene expression of these factors between the anterior and posterior clusters. Further, each identified developmental SA and CT clusters showed distinguishable life span trajectories in a larger longitudinal dataset (4-88 years, 1633 observations), and the SA and CT clusters showed differential relationships to cognitive functions. This means that regions that developed together in childhood also changed together throughout life, demonstrating continuity in regionalization of cortical changes. The AP divide in SA development also characterized genetic patterning obtained in an adult twin sample. In conclusion, the development of cortical regionalization is a continuous process from the embryonic stage throughout life.

Temporal Embryonic Origin Critically Determines Cellular Physiology in the Dentate Gyrus

The dentate gyrus, the entry gate to the hippocampus, comprises 3 types of glutamatergic cells, the granule, the mossy and the semilunar granule cells. Whereas accumulating evidence indicates that specification of subclasses of neocortical neurons starts at the time of their final mitotic divisions, when cellular diversity is specified in the Dentate Gyrus remains largely unknown. Here we show that semilunar cells, like mossy cells, originate from the earliest stages of developmental neurogenesis and that early born neurons form age-matched circuits with each other. Besides morphology, adult semilunar cells display characteristic electrophysiological features that differ from most neurons but are shared among early born granule cells. Therefore, an early birthdate specifies adult granule cell physiology and connectivity whereas additional factors may combine to produce morphological identity.

Region and Cell Type Distribution of TCF4 in the Postnatal Mouse Brain

Transcription factor 4 is a class I basic helix-loop-helix transcription factor regulating gene expression. Altered TCF4 gene expression has been linked to non-syndromic intellectual disability, schizophrenia, and a severe neurodevelopmental disorder known as Pitt-Hopkins syndrome. An understanding of the cell types expressing TCF4 protein in the mouse brain is needed to help identify potential pathophysiological mechanisms and targets for therapeutic delivery in TCF4-linked disorders. Here we developed a novel green fluorescent protein reporter mouse to visualize TCF4-expressing cells throughout the brain. Using this TCF4 reporter mouse, we observed prominent expression of TCF4 in the pallial region and cerebellum of the postnatal brain. At the cellular level, both glutamatergic and GABAergic neurons express TCF4 in the cortex and hippocampus, while only a subset of GABAergic interneurons express TCF4 in the striatum. Among glial cell groups, TCF4 is present in astrocytes and immature and mature oligodendrocytes. In the cerebellum, cells in the granule and molecular layer express TCF4. Our findings greatly extend our knowledge of the spatiotemporal and cell type-specific expression patterns of TCF4 in the brain, and hence, lay the groundwork to better understand TCF4-linked neurological disorders. Any effort to restore TCF4 functions through small molecule or genetic therapies should target these brain regions and cell groups to best recapitulate TCF4 expression patterns.

Environmental Factors Shape Sleep EEG Connectivity During Early Adolescence

Quantifying the degree to which genetic and environmental factors shape brain network connectivity is critical to furthering our understanding of the developing human brain. Sleep, a state of sensory disengagement, provides a unique opportunity to study brain network activity noninvasively by means of sleep electroencephalography (EEG) coherence. We conducted a high-density sleep EEG study in monozygotic (MZ; n = 38; mean age = 12.46; 20 females) and dizygotic (DZ; n = 24; mean age = 12.50; 12 females) twins to assess the heritability of sleep EEG coherence in early adolescence—a period of significant brain rewiring. Structural equation modeling was used to estimate three latent factors: genes, environmental factors shared between twins and environmental factors unique to each twin. We found a strong contribution of unique environmental factors (66% of the variance) and moderate genetic influence (19% of the variance) on sleep EEG coherence across frequencies and sleep states. An exception to this was sleep spindle activity, an index of the thalamocortical network, which showed on average a genetic contribution of 48% across connections. Furthermore, we observed high intraindividual stability of coherence across two consecutive nights suggesting that despite only a modest genetic contribution, sleep EEG coherence is like a trait. Our findings in adolescent humans are in line with earlier findings in animals that show the primordial cerebral map and its connections are plastic and it is through interaction with the environment that the pattern of brain network connectivity is shaped. Therefore, even in twins living together, small differences in the environment may cascade into meaningful differences in brain connectivity.

Dynamic N6-methyladenosine RNA methylation in brain and diseases

N6-methyladenosine (m⁶A) is a dynamic RNA modification that regulates various aspects of RNA metabolism and has been implicated in many biological processes and transitions. m⁶A is highly abundant in the brain; however, only recently has the role of m⁶A in brain development been a focus. The machinery that controls m⁶A is critically important for proper neurodevelopment, and the precise mechanisms by which m⁶A regulates these processes are starting to emerge. However, the role of m⁶A in neurodegenerative and neuropsychiatric diseases still requires much elucidation. This review discusses and summarizes the current body of knowledge surrounding the function of the m⁶A modification in regulating normal brain development, neurodegenerative diseases and outlines possible future directions.

Cortical Malformations: Lessons in Human Brain Development

Creating a functional cerebral cortex requires a series of complex and well-coordinated developmental steps. These steps have evolved across species with the emergence of cortical gyrification and coincided with more complex behaviors. The presence of diverse progenitor cells, a protracted timeline for neuronal migration and maturation, and diverse neuronal types are developmental features that have emerged in the gyrated cortex. These factors could explain how the human brain has expanded in size and complexity. However, their complex nature also renders new avenues of vulnerability by providing additional cell types that could contribute to disease and longer time windows that could impact the composition and organization of the cortical circuit. We aim to discuss the unique developmental steps observed in human corticogenesis and propose how disruption of these species-unique processes could lead to malformations of cortical development.

Differential Contribution of Cortical Thickness, Surface Area, and Gyrification to Fluid and Crystallized Intelligence

Human intelligence can be broadly subdivided into fluid (gf) and crystallized (gc) intelligence, each tapping into distinct cognitive abilities. Although neuroanatomical correlates of intelligence have been previously studied, differential contribution of cortical morphologies to gf and gc has not been fully delineated. Here, we tried to disentangle the contribution of cortical thickness, cortical surface area, and cortical gyrification to gf and gc in a large sample of healthy young subjects (n = 740, Human Connectome Project) with high-resolution MRIs, followed by replication in a separate data set with distinct cognitive measures indexing gf and gc. We found that while gyrification in distributed cortical regions had positive association with both gf and gc, surface area and thickness showed more regional associations. Specifically, higher performance in gf was associated with cortical expansion in regions related to working memory, attention, and visuo-spatial processing, while gc was associated with thinner cortex as well as higher cortical surface area in language-related networks. We discuss the results in a framework where "horizontal" cortical expansion enables higher resource allocation, computational capacity, and functional specificity relevant to gf and gc, while lower cortical thickness possibly reflects cortical pruning facilitating "vertical" intracolumnar efficiency in knowledge-based tasks relevant mostly to gc.

Cell-Specific Gene-Expression Profiles and Cortical Thickness in the Human Brain

Neurobiological underpinnings of cortical thickness in the human brain are largely unknown. Here we use cell-type-specific gene markers to evaluate the contribution of 9 neural cell-types in explaining inter-regional variations in cortical thickness and age-related cortical thinning in the adolescent brain. Gene-expression data were derived from the Allen Human Brain Atlas (and validated using the BrainSpan Atlas). Values of cortical thickness/thinning were obtained with magnetic resonance imaging in a sample of 987 adolescents. We show that inter-regional profiles in cortical thickness relate to those in the expression of genes marking CA1 pyramidal cells, astrocytes, and microglia; taken together, the 3 cell types explain 70% of regional variation in cortical thickness. We also show that inter-regional profiles in cortical thinning relate to those in the expression of genes marking CA1 and S1 pyramidal cells, astrocytes and microglia. Using Gene Ontology analysis, we demonstrate that the difference in the contribution of CA1 and S1 pyramidal cells may relate to biological processes such as neuronal plasticity and potassium channel activity, respectively. This "virtual histology" approach (scripts provided) can be used to examine neurobiological underpinnings of cortical profiles associated with development, aging, and various disorders.

Brain Development and Stochastic Processes During Prenatal and Early Life: You Can't Lose It if You've Never Had It; But It's Better to Have It and Lose It, Than Never to Have Had It at All

Brain development, although largely driven by genetic processes, also is influenced by environmental factors. However, there has been little discussion in the psychiatric literature on the role of stochastic, or chance, events that take place during neurodevelopment. Studies suggest that the brain capitalizes on and regulates the extent of stochastic processes during development. Furthermore, because neurodevelopment is influenced by environmental factors, there is emerging evidence that fostering those positive environmental factors during prenatal and early life could optimize neurodevelopment and provide greater resilience, including those potentially resulting from stochastic processes. Evidence for the role of environmental factors in optimizing early brain development is supported by work in large population-based studies of child development, randomized control trials in high-risk populations, and early-life adoption studies. The public health message is that creating an environment that fosters optimal brain development during prenatal and early life could prevent psychopathology and provide the developing brain the best chance against negative stochastic processes and potential stressors that are inevitable later in life.

New insights into the development of the human cerebral cortex

The cerebral cortex constitutes more than half the volume of the human brain and is presumed to be responsible for the neuronal computations underlying complex phenomena, such as perception, thought, language, attention, episodic memory and voluntary movement. Rodent models are extremely valuable for the investigation of brain development, but cannot provide insight into aspects that are unique or highly derived in humans. Many human psychiatric and neurological conditions have developmental origins but cannot be studied adequately in animal models. The human cerebral cortex has some unique genetic, molecular, cellular and anatomical features, which need to be further explored. The Anatomical Society devoted its summer meeting to the topic of Human Brain Development in June 2018 to tackle these important issues. The meeting was organized by Gavin Clowry (Newcastle University) and Zoltán Molnár (University of Oxford), and held at St John's College, Oxford. The participants provided a broad overview of the structure of the human brain in the context of scaling relationships across the brains of mammals, conserved principles and recent changes in the human lineage. Speakers considered how neuronal progenitors diversified in human to generate an increasing variety of cortical neurons. The formation of the earliest cortical circuits of the earliest generated neurons in the subplate was discussed together with their involvement in neurodevelopmental pathologies. Gene expression networks and susceptibility genes associated to neurodevelopmental diseases were discussed and compared with the networks that can be identified in organoids developed from induced pluripotent stem cells that recapitulate some aspects of in vivo development. New views were discussed on the specification of glutamatergic pyramidal and γ-aminobutyric acid (GABA)ergic interneurons. With the advancement of various in vivo imaging methods, the histopathological observations can be now linked to in vivo normal conditions and to various diseases. Our review gives a general evaluation of the exciting new developments in these areas. The human cortex has a much enlarged association cortex with greater interconnectivity of cortical areas with each other and with an expanded thalamus. The human cortex has relative enlargement of the upper layers, enhanced diversity and function of inhibitory interneurons and a highly expanded transient subplate layer during development. Here we highlight recent studies that address how these differences emerge during development focusing on diverse facets of our evolution.

DOT1L promotes progenitor proliferation and primes neuronal layer identity in the developing cerebral cortex

Cortical development is controlled by transcriptional programs, which are orchestrated by transcription factors. Yet, stable inheritance of spatio-temporal activity of factors influencing cell fate and localization in different layers is only partly understood. Here we find that deletion of Dot1l in the murine telencephalon leads to cortical layering defects, indicating DOT1L activity and chromatin methylation at H3K79 impact on the cell cycle, and influence transcriptional programs conferring upper layer identity in early progenitors. Specifically, DOT1L prevents premature differentiation by increasing expression of genes that regulate asymmetric cell division (Vangl2, Cenpj). Loss of DOT1L results in reduced numbers of progenitors expressing genes including SoxB1 gene family members. Loss of DOT1L also leads to altered cortical distribution of deep layer neurons that express either TBR1, CTIP2 or SOX5, and less activation of transcriptional programs that are characteristic for upper layer neurons (Satb2, Pou3f3, Cux2, SoxC family members). Data from three different mouse models suggest that DOT1L balances transcriptional programs necessary for proper neuronal composition and distribution in the six cortical layers. Furthermore, because loss of DOT1L in the pre-neurogenic phase of development impairs specifically generation of SATB2-expressing upper layer neurons, our data suggest that DOT1L primes upper layer identity in cortical progenitors.