Transcriptional architecture of the primate neocortex

Ecological Trait Differences Are Associated with Gene Expression in the Primary Visual Cortex of Primates

Primate species differ drastically from most other mammals in how they visually perceive their environments, which is particularly important for foraging, predator avoidance, and detection of social cues.

BACKGROUND/OBJECTIVES

Although it is well established that primates display diversity in color vision and various ecological specializations, it is not understood how visual system characteristics and ecological adaptations may be associated with gene expression levels within the primary visual cortex (V1).

METHODS

We performed RNA-Seq on V1 tissue samples from 28 individuals, representing 13 species of primates, including hominoids, cercopithecoids, and platyrrhines. We explored trait-dependent differential expression (DE) by contrasting species with differing visual system phenotypes and ecological traits.

RESULTS

Between 4-25% of genes were determined to be differentially expressed in primates that varied in type of color vision (trichromatic or polymorphic di/trichromatic), habitat use (arboreal or terrestrial), group size (large or small), and primary diet (frugivorous, folivorous, or omnivorous).

CONCLUSIONS

Interestingly, our DE analyses revealed that humans and chimpanzees showed the most marked differences between any two species, even though they are only separated by 6-8 million years of independent evolution. These results show a combination of species-specific and trait-dependent differences in the evolution of gene expression in the primate visual cortex.

eLemur: A cellular-resolution 3D atlas of the mouse lemur brain

The gray mouse lemur (Microcebus murinus), one of the smallest living primates, emerges as a promising model organism for neuroscience research. This is due to its genetic similarity to humans, its evolutionary position between rodents and humans, and its primate-like features encapsulated within a rodent-sized brain. Despite its potential, the absence of a comprehensive reference brain atlas impedes the progress of research endeavors in this species, particularly at the microscopic level. Existing references have largely been confined to the macroscopic scale, lacking detailed anatomical information. Here, we present eLemur, a unique resource, comprising a repository of high-resolution brain-wide images immunostained with multiple cell type and structural markers, elucidating the cyto- and chemoarchitecture of the mouse lemur brain. Additionally, it encompasses a segmented two-dimensional reference and 3D anatomical brain atlas delineated into cortical, subcortical, and other vital regions. Furthermore, eLemur includes a comprehensive 3D cell atlas, providing densities and spatial distributions of non-neuronal and neuronal cells across the mouse lemur brain. Accessible via a web-based viewer (https://eeum-brain.com/#/lemurdatasets), the eLemur resource streamlines data sharing and integration, fostering the exploration of different hypotheses and experimental designs using the mouse lemur as a model organism. Moreover, in conjunction with the growing 3D datasets for rodents, nonhuman primates, and humans, our eLemur 3D digital framework enhances the potential for comparative analysis and translation research, facilitating the integration of extensive rodent study data into human studies.

Abnormal structural‒functional coupling patterning in progressive supranuclear palsy is associated with diverse gradients and histological features

The anatomy of the brain supports inherent processes, fostering mental abilities and eventually facilitating adaptive behavior. Recent studies have shown that progressive supranuclear palsy (PSP) is accompanied by alterations in functional and structural networks. However, how the structure and function of PSP coordinates change is not clear, and the relationships between structural‒functional coupling (SFC) and the gradient of hierarchical structure and cellular histology remain largely unknown. Here, we use neuroimaging data from two independent cohorts and a public histological dataset to investigate the relationships among the cellular histology, hierarchical structure, and SFC of PSP patients. We find that the SFC of the entire cortex in PSP is severely disrupted, with higher coupling in the visual network (VN). Moreover, coupling differences in PSP follow a macroscopic organizational principle from unimodal to transmodal gradients. Finally, we elucidate greater laminar differentiation in VN regions sensitive to SFC changes in PSP, which is related mainly to the higher cellular density and smaller size of the internal-granular layer. In conclusion, our findings provide an interpretable framework for understanding SFC changes in PSP and provide new insights into the consistency of structural and functional changes in PSP regarding hierarchical structure and cellular histology.

State-specific Regulation of Electrical Stimulation in the Intralaminar Thalamus of Macaque Monkeys: Network and Transcriptional Insights into Arousal

Long-range thalamocortical communication is central to anesthesia-induced loss of consciousness and its reversal. However, isolating the specific neural networks connecting thalamic nuclei with various cortical regions for state-specific anesthesia regulation is challenging, with the biological underpinnings still largely unknown. Here, simultaneous electroencephalogram-fuctional magnetic resonance imaging (EEG-fMRI) and deep brain stimulation are applied to the intralaminar thalamus in macaques under finely-tuned propofol anesthesia. This approach led to the identification of an intralaminar-driven network responsible for rapid arousal during slow-wave oscillations. A network-based RNA-sequencing analysis is conducted of region-, layer-, and cell-specific gene expression data from independent transcriptomic atlases and identifies 2489 genes preferentially expressed within this arousal network, notably enriched in potassium channels and excitatory, parvalbumin-expressing neurons, and oligodendrocytes. Comparison with human RNA-sequencing data highlights conserved molecular and cellular architectures that enable the matching of homologous genes, protein interactions, and cell types across primates, providing novel insight into network-focused transcriptional signatures of arousal.

Epigenomic landscapes during prefrontal cortex development and aging in rhesus

The prefrontal cortex (PFC) is essential for higher-level cognitive functions. How epigenetic dynamics participates in PFC development and aging is largely unknown. Here, we profiled epigenomic landscapes of rhesus monkey PFCs from prenatal to aging stages. The dynamics of chromatin states, including higher-order chromatin structure, chromatin interaction and histone modifications are coordinated to regulate stage-specific gene transcription, participating in distinct processes of neurodevelopment. Dramatic changes of epigenetic signals occur around the birth stage. Notably, genes involved in neuronal cell differentiation and layer specification are pre-configured by bivalent promoters. We identified a cis-regulatory module and the transcription factors (TFs) associated with basal radial glia development, which was associated with large brain size in primates. These TFs include GLI3, CREB5 and SOX9. Interestingly, the genes associated with the basal radial glia (bRG)-associated cis-element module, such as SRY and SOX9, are enriched in sex differentiation. Schizophrenia-associated single nucleotide polymorphisms are more enriched in super enhancers (SEs) than typical enhancers, suggesting that SEs play an important role in neural network wiring. A cis-regulatory element of DBN1 is identified, which is critical for neuronal cell proliferation and synaptic neuron differentiation. Notably, the loss of distal chromatin interaction and H3K27me3 signal are hallmarks of PFC aging, which are associated with abnormal expression of aging-related genes and transposon activation, respectively. Collectively, our findings shed light on epigenetic mechanisms underlying primate brain development and aging.

Lamination, Borders, and Thalamic Projections of the Primary Visual Cortex in Human, Non-Human Primate, and Rodent Brains

The primary visual cortex (V1) is one of the most studied regions of the brain and is characterized by its specialized and laminated layer 4 in human and non-human primates. However, studies aiming to harmonize the definition of the cortical layers and borders of V1 across rodents and primates are very limited. This article attempts to identify and harmonize the molecular markers and connectional patterns that can consistently link corresponding cortical layers of V1 and borders across mammalian species and ages. V1 in primates has at least two additional and unique layers (L3b2 and L3c) and two sublayers of layer 4 (L4a and L4b) compared to rodent V1. In all species examined, layers 4 and 3b of V1 receive strong inputs from the (dorsal) lateral geniculate nucleus, and V1 is mostly surrounded by the secondary visual cortex except for one location where V1 directly abuts area prostriata. The borders of primate V1 can also be clearly identified at mid-gestational ages using gene markers. In rodents, a novel posteromedial extension of V1 is identified, which expresses V1 marker genes and receives strong inputs from the lateral geniculate nucleus. This V1 extension was labeled as the posterior retrosplenial cortex and medial secondary visual cortex in the literature and brain atlases. Layer 6 of the rodent and primate V1 originates corticothalamic projections to the lateral geniculate, lateral dorsal, and reticular thalamic nuclei and the lateroposterior-pulvinar complex with topographic organization. Finally, the direct geniculo-extrastriate (particularly the strong geniculo-prostriata) projections are probably major contributors to blindsight after V1 lesions. Taken together, compared to rodents, primates, and humans, V1 has at least two unique middle layers, while other layers are comparable across species and display conserved molecular markers and similar connections with the visual thalamus with only subtle differences.

Large-scale coupling of prefrontal activity patterns as a mechanism for cognitive control in health and disease: evidence from rodent models

Cognitive control of behavior is crucial for well-being, as allows subject to adapt to changing environments in a goal-directed way. Changes in cognitive control of behavior is observed during cognitive decline in elderly and in pathological mental conditions. Therefore, the recovery of cognitive control may provide a reliable preventive and therapeutic strategy. However, its neural basis is not completely understood. Cognitive control is supported by the prefrontal cortex, structure that integrates relevant information for the appropriate organization of behavior. At neurophysiological level, it is suggested that cognitive control is supported by local and large-scale synchronization of oscillatory activity patterns and neural spiking activity between the prefrontal cortex and distributed neural networks. In this review, we focus mainly on rodent models approaching the neuronal origin of these prefrontal patterns, and the cognitive and behavioral relevance of its coordination with distributed brain systems. We also examine the relationship between cognitive control and neural activity patterns in the prefrontal cortex, and its role in normal cognitive decline and pathological mental conditions. Finally, based on these body of evidence, we propose a common mechanism that may underlie the impaired cognitive control of behavior.

Selective vulnerability of layer 5a corticostriatal neurons in Huntington's disease

The properties of the cell types that are selectively vulnerable in Huntington's disease (HD) cortex, the nature of somatic CAG expansions of mHTT in these cells, and their importance in CNS circuitry have not been delineated. Here, we employed serial fluorescence-activated nuclear sorting (sFANS), deep molecular profiling, and single-nucleus RNA sequencing (snRNA-seq) of motor-cortex samples from thirteen predominantly early stage, clinically diagnosed HD donors and selected samples from cingulate, visual, insular, and prefrontal cortices to demonstrate loss of layer 5a pyramidal neurons in HD. Extensive mHTT CAG expansions occur in vulnerable layer 5a pyramidal cells, and in Betz cells, layers 6a and 6b neurons that are resilient in HD. Retrograde tracing experiments in macaque brains identify layer 5a neurons as corticostriatal pyramidal cells. We propose that enhanced somatic mHTT CAG expansion and altered synaptic function act together to cause corticostriatal disconnection and selective neuronal vulnerability in HD cerebral cortex.

Heterogeneity of layer 4 in visual areas of rhesus macaque cortex

Recently, single-cell RNA-sequencing (scRNA-seq) has enabled unprecedented insights to the cellular landscape of the brains of many different species, among them the rhesus macaque as a key animal model. Building on previous, broader surveys of the macaque brain, we closely examined five immediately neighboring areas within the visual cortex of the rhesus macaque: V1, V2, V4, MT and TEO. To facilitate this, we first devised a novel pipeline for brain spatial archive - the BrainSPACE - which enabled robust archiving and sampling from the whole unfixed brain. SnRNA-sequencing of ~100,000 nuclei from visual areas V1 and V4 revealed conservation within the GABAergic neuron subtypes, while seven and one distinct principle neuron subtypes were detected in V1 and V4, respectively, all most likely located in layer 4. Moreover, using small molecule fluorescence in situ hybridization, we identified cell type density gradients across V1, V2, V4, MT, and TEO appearing to reflect the visual hierarchy. These findings demonstrate an association between the clear areal specializations among neighboring areas with the hierarchical levels within the visual cortex of the rhesus macaque.

Temporal transcriptomic dynamics in developing macaque neocortex

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

Divergent Evolutionary Rates of Primate Brain Regions as Revealed by Genomics and Transcriptomics

Although the primate brain contains numerous functionally distinct structures that have experienced diverse genetic changes during the course of evolution and development, these changes remain to be explored in detail. Here we utilize two classic metrics from evolutionary biology, the evolutionary rate index (ERI) and the transcriptome age index (TAI), to investigate the evolutionary alterations that have occurred in each area and developmental stage of the primate brain. We observed a higher evolutionary rate for those genes expressed in the non-cortical areas during primate evolution, particularly in human, with the highest rate of evolution being exhibited at brain developmental stages between late infancy and early childhood. Further, the transcriptome age of the non-cortical areas was lower than that of the cerebral cortex, with the youngest age apparent at brain developmental stages between late infancy and early childhood. Our exploration of the evolutionary patterns manifest in each brain area and developmental stage provides important reference points for further research into primate brain evolution.

Gene expression in the primate orbitofrontal cortex related to anxious temperament

Anxiety disorders are among the most prevalent psychiatric disorders, causing significant suffering and disability. Relative to other psychiatric disorders, anxiety disorders tend to emerge early in life, supporting the importance of developmental mechanisms in their emergence and maintenance. Behavioral inhibition (BI) is a temperament that emerges early in life and, when stable and extreme, is linked to an increased risk for the later development of anxiety disorders and other stress-related psychopathology. Understanding the neural systems and molecular mechanisms underlying this dispositional risk could provide insight into treatment targets for anxiety disorders. Nonhuman primates (NHPs) have an anxiety-related temperament, called anxious temperament (AT), that is remarkably similar to BI in humans, facilitating the design of highly translational models for studying the early risk for stress-related psychopathology. Because of the recent evolutionary divergence between humans and NHPs, many of the anxiety-related brain regions that contribute to psychopathology are highly similar in terms of their structure and function, particularly with respect to the prefrontal cortex. The orbitofrontal cortex plays a critical role in the flexible encoding and regulation of threat responses, in part through connections with subcortical structures like the amygdala. Here, we explore individual differences in the transcriptional profile of cells within the region, using laser capture microdissection and single nuclear sequencing, providing insight into the molecules underlying individual differences in AT-related function of the pOFC, with a particular focus on previously implicated cellular systems, including neurotrophins and glucocorticoid signaling.

Selective Vulnerability of Layer 5a Corticostriatal Neurons in Huntington's Disease

The properties of the cell types that are selectively vulnerable in Huntington's disease (HD) cortex, the nature of somatic CAG expansions of mHTT in these cells, and their importance in CNS circuitry have not been delineated. Here we employed serial fluorescence activated nuclear sorting (sFANS), deep molecular profiling, and single nucleus RNA sequencing (snRNAseq) to demonstrate that layer 5a pyramidal neurons are vulnerable in primary motor cortex and other cortical areas of HD donors. Extensive mHTT -CAG expansions occur in vulnerable layer 5a pyramidal cells, and in Betz cells, layer 6a, layer 6b neurons that are resilient in HD. Retrograde tracing experiments in macaque brains identify the vulnerable layer 5a neurons as corticostriatal pyramidal cells. We propose that enhanced somatic mHTT -CAG expansion and altered synaptic function act together to cause corticostriatal disconnection and selective neuronal vulnerability in the HD cerebral cortex.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

A role for cortical interneurons as adversarial discriminators

The brain learns representations of sensory information from experience, but the algorithms by which it does so remain unknown. One popular theory formalizes representations as inferred factors in a generative model of sensory stimuli, meaning that learning must improve this generative model and inference procedure. This framework underlies many classic computational theories of sensory learning, such as Boltzmann machines, the Wake/Sleep algorithm, and a more recent proposal that the brain learns with an adversarial algorithm that compares waking and dreaming activity. However, in order for such theories to provide insights into the cellular mechanisms of sensory learning, they must be first linked to the cell types in the brain that mediate them. In this study, we examine whether a subtype of cortical interneurons might mediate sensory learning by serving as discriminators, a crucial component in an adversarial algorithm for representation learning. We describe how such interneurons would be characterized by a plasticity rule that switches from Hebbian plasticity during waking states to anti-Hebbian plasticity in dreaming states. Evaluating the computational advantages and disadvantages of this algorithm, we find that it excels at learning representations in networks with recurrent connections but scales poorly with network size. This limitation can be partially addressed if the network also oscillates between evoked activity and generative samples on faster timescales. Consequently, we propose that an adversarial algorithm with interneurons as discriminators is a plausible and testable strategy for sensory learning in biological systems.

Cortical microstructural gradients capture memory network reorganization in temporal lobe epilepsy

Temporal lobe epilepsy (TLE), one of the most common pharmaco-resistant epilepsies, is associated with pathology of paralimbic brain regions, particularly in the mesiotemporal lobe. Cognitive dysfunction in TLE is frequent, and particularly affects episodic memory. Crucially, these difficulties challenge the quality of life of patients, sometimes more than seizures, underscoring the need to assess neural processes of cognitive dysfunction in TLE to improve patient management. Our work harnessed a novel conceptual and analytical approach to assess spatial gradients of microstructural differentiation between cortical areas based on high-resolution MRI analysis. Gradients track region-to-region variations in intracortical lamination and myeloarchitecture, serving as a system-level measure of structural and functional reorganization. Comparing cortex-wide microstructural gradients between 21 patients and 35 healthy controls, we observed a reorganization of this gradient in TLE driven by reduced microstructural differentiation between paralimbic cortices and the remaining cortex with marked abnormalities in ipsilateral temporopolar and dorsolateral prefrontal regions. Findings were replicated in an independent cohort. Using an independent post-mortem dataset, we observed that in vivo findings reflected topographical variations in cortical cytoarchitecture. We indeed found that macroscale changes in microstructural differentiation in TLE reflected increased similarity of paralimbic and primary sensory/motor regions. Disease-related transcriptomics could furthermore show specificity of our findings to TLE over other common epilepsy syndromes. Finally, microstructural dedifferentiation was associated with cognitive network reorganization seen during an episodic memory functional MRI paradigm and correlated with interindividual differences in task accuracy. Collectively, our findings showing a pattern of reduced microarchitectural differentiation between paralimbic regions and the remaining cortex provide a structurally-grounded explanation for large-scale functional network reorganization and cognitive dysfunction characteristic of TLE.

Single-cell spatial transcriptome reveals cell-type organization in the macaque cortex

Elucidating the cellular organization of the cerebral cortex is critical for understanding brain structure and function. Using large-scale single-nucleus RNA sequencing and spatial transcriptomic analysis of 143 macaque cortical regions, we obtained a comprehensive atlas of 264 transcriptome-defined cortical cell types and mapped their spatial distribution across the entire cortex. We characterized the cortical layer and region preferences of glutamatergic, GABAergic, and non-neuronal cell types, as well as regional differences in cell-type composition and neighborhood complexity. Notably, we discovered a relationship between the regional distribution of various cell types and the region's hierarchical level in the visual and somatosensory systems. Cross-species comparison of transcriptomic data from human, macaque, and mouse cortices further revealed primate-specific cell types that are enriched in layer 4, with their marker genes expressed in a region-dependent manner. Our data provide a cellular and molecular basis for understanding the evolution, development, aging, and pathogenesis of the primate brain.

Integrative Omics Reveals Rapidly Evolving Regulatory Sequences Driving Primate Brain Evolution

Although the continual expansion of the brain during primate evolution accounts for our enhanced cognitive capabilities, the drivers of brain evolution have scarcely been explored in these ancestral nodes. Here, we performed large-scale comparative genomic, transcriptomic, and epigenomic analyses to investigate the evolutionary alterations acquired by brain genes and provide comprehensive listings of innovatory genetic elements along the evolutionary path from ancestral primates to human. The regulatory sequences associated with brain-expressed genes experienced rapid change, particularly in the ancestor of the Simiiformes. Extensive comparisons of single-cell and bulk transcriptomic data between primate and nonprimate brains revealed that these regulatory sequences may drive the high expression of certain genes in primate brains. Employing in utero electroporation into mouse embryonic cortex, we show that the primate-specific brain-biased gene BMP7 was recruited, probably in the ancestor of the Simiiformes, to regulate neuronal proliferation in the primate ventricular zone. Our study provides a comprehensive listing of genes and regulatory changes along the brain evolution lineage of ancestral primates leading to human. These data should be invaluable for future functional studies that will deepen our understanding not only of the genetic basis of human brain evolution but also of inherited disease.

Brain-wide and cell-specific transcriptomic insights into MRI-derived cortical morphology in macaque monkeys

Integrative analyses of transcriptomic and neuroimaging data have generated a wealth of information about biological pathways underlying regional variability in imaging-derived brain phenotypes in humans, but rarely in nonhuman primates due to the lack of a comprehensive anatomically-defined atlas of brain transcriptomics. Here we generate complementary bulk RNA-sequencing dataset of 819 samples from 110 brain regions and single-nucleus RNA-sequencing dataset, and neuroimaging data from 162 cynomolgus macaques, to examine the link between brain-wide gene expression and regional variation in morphometry. We not only observe global/regional expression profiles of macaque brain comparable to human but unravel a dorsolateral-ventromedial gradient of gene assemblies within the primate frontal lobe. Furthermore, we identify a set of 971 protein-coding and 34 non-coding genes consistently associated with cortical thickness, specially enriched for neurons and oligodendrocytes. These data provide a unique resource to investigate nonhuman primate models of human diseases and probe cross-species evolutionary mechanisms.

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.

Differential gene expression between callosal and ipsilateral projection neurons in the monkey dorsolateral prefrontal and posterior parietal cortices

Reciprocal connections between primate dorsolateral prefrontal (DLPFC) and posterior parietal (PPC) cortices, furnished by subsets of layer 3 pyramidal neurons (PNs), contribute to cognitive processes including working memory (WM). A different subset of layer 3 PNs in each region projects to the homotopic region of the contralateral hemisphere. These ipsilateral (IP) and callosal (CP) projections, respectively, appear to be essential for the maintenance and transfer of information during WM. To determine if IP and CP layer 3 PNs in each region differ in their transcriptomes, fluorescent retrograde tracers were used to label IP and CP layer 3 PNs in the DLPFC and PPC from macaque monkeys. Retrogradely-labeled PNs were captured by laser microdissection and analyzed by RNAseq. Numerous differentially expressed genes (DEGs) were detected between IP and CP neurons in each region and the functional pathways containing many of these DEGs were shared across regions. However, DLPFC and PPC displayed opposite patterns of DEG enrichment between IP and CP neurons. Cross-region analyses indicated that the cortical area targeted by IP or CP layer 3 PNs was a strong correlate of their transcriptome profile. These findings suggest that the transcriptomes of layer 3 PNs reflect regional, projection type and target region specificity.

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.

Identification of visual cortex cell types and species differences using single-cell RNA sequencing

The primate neocortex exerts high cognitive ability and strong information processing capacity. Here, we establish a single-cell RNA sequencing dataset of 133,454 macaque visual cortical cells. It covers major cortical cell classes including 25 excitatory neuron types, 37 inhibitory neuron types and all glial cell types. We identified layer-specific markers including HPCAL1 and NXPH4, and also identified two cell types, an NPY-expressing excitatory neuron type that expresses the dopamine receptor D3 gene; and a primate specific activity-dependent OSTN + sensory neuron type. Comparisons of our dataset with humans and mice show that the gene expression profiles differ between species in relation to genes that are implicated in the synaptic plasticity and neuromodulation of excitatory neurons. The comparisons also revealed that glutamatergic neurons may be more diverse across species than GABAergic neurons and non-neuronal cells. These findings pave the way for understanding how the primary cortex fulfills the high-cognitive functions.

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.

Broad transcriptomic dysregulation occurs across the cerebral cortex in ASD

Neuropsychiatric disorders classically lack defining brain pathologies, but recent work has demonstrated dysregulation at the molecular level, characterized by transcriptomic and epigenetic alterations1-3. In autism spectrum disorder (ASD), this molecular pathology involves the upregulation of microglial, astrocyte and neural-immune genes, the downregulation of synaptic genes, and attenuation of gene-expression gradients in cortex1,2,4-6. However, whether these changes are limited to cortical association regions or are more widespread remains unknown. To address this issue, we performed RNA-sequencing analysis of 725 brain samples spanning 11 cortical areas from 112 post-mortem samples from individuals with ASD and neurotypical controls. We find widespread transcriptomic changes across the cortex in ASD, exhibiting an anterior-to-posterior gradient, with the greatest differences in primary visual cortex, coincident with an attenuation of the typical transcriptomic differences between cortical regions. Single-nucleus RNA-sequencing and methylation profiling demonstrate that this robust molecular signature reflects changes in cell-type-specific gene expression, particularly affecting excitatory neurons and glia. Both rare and common ASD-associated genetic variation converge within a downregulated co-expression module involving synaptic signalling, and common variation alone is enriched within a module of upregulated protein chaperone genes. These results highlight widespread molecular changes across the cerebral cortex in ASD, extending beyond association cortex to broadly involve primary sensory regions.

Genomics, convergent neuroscience and progress in understanding autism spectrum disorder

More than a hundred genes have been identified that, when disrupted, impart large risk for autism spectrum disorder (ASD). Current knowledge about the encoded proteins - although incomplete - points to a very wide range of developmentally dynamic and diverse biological processes. Moreover, the core symptoms of ASD involve distinctly human characteristics, presenting challenges to interpreting evolutionarily distant model systems. Indeed, despite a decade of striking progress in gene discovery, an actionable understanding of pathobiology remains elusive. Increasingly, convergent neuroscience approaches have been recognized as an important complement to traditional uses of genetics to illuminate the biology of human disorders. These methods seek to identify intersection among molecular-level, cellular-level and circuit-level functions across multiple risk genes and have highlighted developing excitatory neurons in the human mid-gestational prefrontal cortex as an important pathobiological nexus in ASD. In addition, neurogenesis, chromatin modification and synaptic function have emerged as key potential mediators of genetic vulnerability. The continued expansion of foundational 'omics' data sets, the application of higher-throughput model systems and incorporating developmental trajectories and sex differences into future analyses will refine and extend these results. Ultimately, a systems-level understanding of ASD genetic risk holds promise for clarifying pathobiology and advancing therapeutics.

New pathogenic insights from large animal models of neurodegenerative diseases

Animal models are essential for investigating the pathogenesis and developing the treatment of human diseases. Identification of genetic mutations responsible for neurodegenerative diseases has enabled the creation of a large number of small animal models that mimic genetic defects found in the affected individuals. Of the current animal models, rodents with genetic modifications are the most commonly used animal models and provided important insights into pathogenesis. However, most of genetically modified rodent models lack overt neurodegeneration, imposing challenges and obstacles in utilizing them to rigorously test the therapeutic effects on neurodegeneration. Recent studies that used CRISPR/Cas9-targeted large animal (pigs and monkeys) have uncovered important pathological events that resemble neurodegeneration in the patient’s brain but could not be produced in small animal models. Here we highlight the unique nature of large animals to model neurodegenerative diseases as well as the limitations and challenges in establishing large animal models of neurodegenerative diseases, with focus on Huntington disease, Amyotrophic lateral sclerosis, and Parkinson diseases. We also discuss how to use the important pathogenic insights from large animal models to make rodent models more capable of recapitulating important pathological features of neurodegenerative diseases.

Differential gene expression in layer 3 pyramidal neurons across 3 regions of the human cortical visual spatial working memory network

Visual spatial working memory (vsWM) is mediated by a distributed cortical network composed of multiple nodes, including primary visual (V1), posterior parietal (PPC), and dorsolateral prefrontal (DLPFC) cortices. Feedforward and feedback information is transferred among these nodes via projections furnished by pyramidal neurons (PNs) located primarily in cortical layer 3. Morphological and electrophysiological differences among layer 3 PNs across these nodes have been reported; however, the transcriptional signatures underlying these differences have not been examined in the human brain. Here we interrogated the transcriptomes of layer 3 PNs from 39 neurotypical human subjects across 3 critical nodes of the vsWM network. Over 8,000 differentially expressed genes were detected, with more than 6,000 transcriptional differences present between layer 3 PNs in V1 and those in PPC and DLPFC. Additionally, over 600 other genes differed in expression along the rostral-to-caudal hierarchy formed by these 3 nodes. Moreover, pathway analysis revealed enrichment of genes in V1 related to circadian rhythms and in DLPFC of genes involved in synaptic plasticity. Overall, these results show robust regional differences in the transcriptome of layer 3 PNs, which likely contribute to regional specialization in their morphological and physiological features and thus in their functional contributions to vsWM.

Evaluation of deep convolutional neural networks for in situ hybridization gene expression image representation

High resolution in situ hybridization (ISH) images of the brain capture spatial gene expression at cellular resolution. These spatial profiles are key to understanding brain organization at the molecular level. Previously, manual qualitative scoring and informatics pipelines have been applied to ISH images to determine expression intensity and pattern. To better capture the complex patterns of gene expression in the human cerebral cortex, we applied a machine learning approach. We propose gene re-identification as a contrastive learning task to compute representations of ISH images. We train our model on an ISH dataset of ~1,000 genes obtained from postmortem samples from 42 individuals. This model reaches a gene re-identification rate of 38.3%, a 13x improvement over random chance. We find that the learned embeddings predict expression intensity and pattern. To test generalization, we generated embeddings in a second dataset that assayed the expression of 78 genes in 53 individuals. In this set of images, 60.2% of genes are re-identified, suggesting the model is robust. Importantly, this dataset assayed expression in individuals diagnosed with schizophrenia. Gene and donor-specific embeddings from the model predict schizophrenia diagnosis at levels similar to that reached with demographic information. Mutations in the most discriminative gene, Sodium Voltage-Gated Channel Beta Subunit 4 (SCN4B), may help understand cardiovascular associations with schizophrenia and its treatment. We have publicly released our source code, embeddings, and models to spur further application to spatial transcriptomics. In summary, we propose and evaluate gene re-identification as a machine learning task to represent ISH gene expression images.

Operating principles of the cerebral cortex as a six-layered network in primates: beyond the classic canonical circuit model

The cerebral cortex performs its computations with many six-layered fundamental units, collectively spreading along the cortical sheet. What is the local network structure and the operating dynamics of such a fundamental unit? Previous investigations of primary sensory areas revealed a classic "canonical" circuit model, leading to an expectation of similar circuit organization and dynamics throughout the cortex. This review clarifies the different circuit dynamics at play in the higher association cortex of primates that implements computation for high-level cognition such as memory and attention. Instead of feedforward processing of response selectivity through Layers 4 to 2/3 that the classic canonical circuit stipulates, memory recall in primates occurs in Layer 5/6 with local backward projection to Layer 2/3, after which the retrieved information is sent back from Layer 6 to lower-level cortical areas for further retrieval of nested associations of target attributes. In this review, a novel "dynamic multimode module (D3M)" in the primate association cortex is proposed, as a new "canonical" circuit model performing this operation.

Relating quantitative 7T MRI across cortical depths to cytoarchitectonics, gene expression and connectomics

Ultra-high field MRI across the depth of the cortex has the potential to provide anatomically precise biomarkers and mechanistic insights into neurodegenerative disease like Huntington's disease that show layer-selective vulnerability. Here we compare multi-parametric mapping (MPM) measures across cortical depths for a 7T 500 μm whole brain acquisition to (a) layer-specific cell measures from the von Economo histology atlas, (b) layer-specific gene expression, using the Allen Human Brain atlas and (c) white matter connections using high-fidelity diffusion tractography, at a 1.3 mm isotropic voxel resolution, from a 300mT/m Connectom MRI system. We show that R2*, but not R1, across cortical depths is highly correlated with layer-specific cell number and layer-specific gene expression. R1- and R2*-weighted connectivity strength of cortico-striatal and intra-hemispheric cortical white matter connections was highly correlated with grey matter R1 and R2* across cortical depths. Limitations of the layer-specific relationships demonstrated are at least in part related to the high cross-correlations of von Economo atlas cell counts and layer-specific gene expression across cortical layers. These findings demonstrate the potential and limitations of combining 7T MPMs, gene expression and white matter connections to provide an anatomically precise framework for tracking neurodegenerative disease.

Comparative cellular analysis of motor cortex in human, marmoset and mouse

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch-seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations.

Primate neuronal connections are sparse in cortex as compared to mouse

Detailing how primate and mouse neurons differ is critical for creating generalized models of how neurons process information. We reconstruct 15,748 synapses in adult Rhesus macaques and mice and ask how connectivity differs on identified cell types in layer 2/3 of primary visual cortex. Primate excitatory and inhibitory neurons receive 2-5 times fewer excitatory and inhibitory synapses than similar mouse neurons. Primate excitatory neurons have lower excitatory-to-inhibitory (E/I) ratios than mouse but similar E/I ratios in inhibitory neurons. In both species, properties of inhibitory axons such as synapse size and frequency are unchanged, and inhibitory innervation of excitatory neurons is local and specific. Using artificial recurrent neural networks (RNNs) optimized for different cognitive tasks, we find that penalizing networks for creating and maintaining synapses, as opposed to neuronal firing, reduces the number of connections per node as the number of nodes increases, similar to primate neurons compared with mice.

Progenitor cell diversity in the developing mouse neocortex

In the mammalian neocortex, projection neuron types are sequentially generated by the same pool of neural progenitors. How neuron type specification is related to developmental timing remains unclear. To determine whether temporal gene expression in neural progenitors correlates with neuron type specification, we performed single-cell RNA sequencing (scRNA-Seq) analysis of the developing mouse neocortex. We uncovered neuroepithelial cell enriched genes such as Hmga2 and Ccnd1 when compared to radial glial cells (RGCs). RGCs display dynamic gene expression over time; for instance, early RGCs express higher levels of Hes5, and late RGCs show higher expression of Pou3f2 Interestingly, intermediate progenitor cell marker gene Eomes coexpresses temporally with known neuronal identity genes at different developmental stages, though mostly in postmitotic cells. Our results delineate neural progenitor cell diversity in the developing mouse neocortex and support that neuronal identity genes are transcriptionally evident in Eomes-positive cells.

Pairwise interactions in gene expression determine a hierarchical transcriptional profile in the human brain

The orchestrated expression of thousands of genes gives rise to the complexity of the human brain. However, the structures governing these myriad gene-gene interactions remain unclear. By analyzing transcription data from more than 2000 sites in six human brains, we found that pairwise interactions between genes, without considering any higher-order interactions, are sufficient to predict the transcriptional pattern of the genome for individual brain regions and the transcriptional profile of the entire brain consisting of more than 200 areas. These findings suggest a quadratic complexity of transcriptional patterns in the human brain, which is much simpler than expected. In addition, using a pairwise interaction model, we revealed that the strength of gene-gene interactions in the human brain gives rise to the nearly maximal number of transcriptional clusters, which may account for the functional and structural richness of the brain.

A comprehensive study of a 29-capsid AAV library in a non-human primate central nervous system

Non-human primates (NHPs) are a preferred animal model for optimizing adeno-associated virus (AAV)-mediated CNS gene delivery protocols before clinical trials. In spite of its inherent appeal, it is challenging to compare different serotypes, delivery routes, and disease indications in a well-powered, comprehensive, multigroup NHP experiment. Here, a multiplex barcode recombinant AAV (rAAV) vector-tracing strategy has been applied to a systemic analysis of 29 distinct, wild-type (WT), AAV natural isolates and engineered capsids in the CNS of eight macaques. The report describes distribution of each capsid in 15 areas of the macaques' CNS after intraparenchymal (putamen) injection, or cerebrospinal fluid (CSF)-mediated administration routes (intracisternal, intrathecal, or intracerebroventricular). To trace the vector biodistribution (viral DNA) and targeted tissues transduction (viral mRNA) of each capsid in each of the analyzed CNS areas, quantitative next-generation sequencing analysis, assisted by the digital-droplet PCR technology, was used. The report describes the most efficient AAV capsid variants targeting specific CNS areas after each route of administration using the direct side-by-side comparison of WT AAV isolates and a new generation of rationally designed capsids. The newly developed bioinformatics and visualization algorithms, applicable to the comparative analysis of several mammalian brain models, have been developed and made available in the public domain.

Prioritizing Genetic Contributors to Cortical Alterations in 22q11.2 Deletion Syndrome Using Imaging Transcriptomics

22q11.2 deletion syndrome (22q11DS) results from a hemizygous deletion that typically spans 46 protein-coding genes and is associated with widespread alterations in brain morphology. The specific genetic mechanisms underlying these alterations remain unclear. In the 22q11.2 ENIGMA Working Group, we characterized cortical alterations in individuals with 22q11DS (n = 232) versus healthy individuals (n = 290) and conducted spatial convergence analyses using gene expression data from the Allen Human Brain Atlas to prioritize individual genes that may contribute to altered surface area (SA) and cortical thickness (CT) in 22q11DS. Total SA was reduced in 22q11DS (Z-score deviance = -1.04), with prominent reductions in midline posterior and lateral association regions. Mean CT was thicker in 22q11DS (Z-score deviance = +0.64), with focal thinning in a subset of regions. Regional expression of DGCR8 was robustly associated with regional severity of SA deviance in 22q11DS; AIFM3 was also associated with SA deviance. Conversely, P2RX6 was associated with CT deviance. Exploratory analysis of gene targets of microRNAs previously identified as down-regulated due to DGCR8 deficiency suggested that DGCR8 haploinsufficiency may contribute to altered corticogenesis in 22q11DS by disrupting cell cycle modulation. These findings demonstrate the utility of combining neuroanatomic and transcriptomic datasets to derive molecular insights into complex, multigene copy number variants.

The impact of (ab)normal maternal environment on cortical development

The cortex in the mammalian brain is the most complex brain region that integrates sensory information and coordinates motor and cognitive processes. To perform such functions, the cortex contains multiple subtypes of neurons that are generated during embryogenesis. Newly born neurons migrate to their proper location in the cortex, grow axons and dendrites, and form neuronal circuits. These developmental processes in the fetal brain are regulated to a large extent by a great variety of factors derived from the mother - starting from simple nutrients as building blocks and ending with hormones. Thus, when the normal maternal environment is disturbed due to maternal infection, stress, malnutrition, or toxic substances, it might have a profound impact on cortical development and the offspring can develop a variety of neurodevelopmental disorders. Here we first describe the major developmental processes which generate neuronal diversity in the cortex. We then review our knowledge of how most common maternal insults affect cortical development, perturb neuronal circuits, and lead to neurodevelopmental disorders. We further present a concept of selective vulnerability of cortical neuronal subtypes to maternal-derived insults, where the vulnerability of cortical neurons and their progenitors to an insult depends on the time (developmental period), place (location in the developing brain), and type (unique features of a cell type and an insult). Finally, we provide evidence for the existence of selective vulnerability during cortical development and identify the most vulnerable neuronal types, stages of differentiation, and developmental time for major maternal-derived insults.

Organisational and neuromodulatory underpinnings of structural-functional connectivity decoupling in patients with Parkinson's disease

Parkinson's dementia is characterised by changes in perception and thought, and preceded by visual dysfunction, making this a useful surrogate for dementia risk. Structural and functional connectivity changes are seen in humans with Parkinson's disease, but the organisational principles are not known. We used resting-state fMRI and diffusion-weighted imaging to examine changes in structural-functional connectivity coupling in patients with Parkinson's disease, and those at risk of dementia. We identified two organisational gradients to structural-functional connectivity decoupling: anterior-to-posterior and unimodal-to-transmodal, with stronger structural-functional connectivity coupling in anterior, unimodal areas and weakened towards posterior, transmodal regions. Next, we related spatial patterns of decoupling to expression of neurotransmitter receptors. We found that dopaminergic and serotonergic transmission relates to decoupling in Parkinson's overall, but instead, serotonergic, cholinergic and noradrenergic transmission relates to decoupling in patients with visual dysfunction. Our findings provide a framework to explain the specific disorders of consciousness in Parkinson's dementia, and the neurotransmitter systems that underlie these.

Distinct Laminar and Cellular Patterns of GABA Neuron Transcript Expression in Monkey Prefrontal and Visual Cortices

The functional output of a cortical region is shaped by its complement of GABA neuron subtypes. GABA-related transcript expression differs substantially between the primate dorsolateral prefrontal cortex (DLPFC) and primary visual (V1) cortices in gray matter homogenates, but the laminar and cellular bases for these differences are unknown. Quantification of levels of GABA-related transcripts in layers 2 and 4 of monkey DLPFC and V1 revealed three distinct expression patterns: 1) transcripts with higher levels in DLPFC and layer 2 [e.g., somatostatin (SST)]; 2) transcripts with higher levels in V1 and layer 4 [e.g., parvalbumin (PV)], and 3) transcripts with similar levels across layers and regions [e.g., glutamic acid decarboxylase (GAD67)]. At the cellular level, these patterns reflected transcript- and cell type-specific differences: the SST pattern primarily reflected differences in the relative proportions of SST mRNA-positive neurons, the PV pattern primarily reflected differences in PV mRNA expression per neuron, and the GAD67 pattern reflected opposed patterns in the relative proportions of GAD67 mRNA-positive neurons and in GAD67 mRNA expression per neuron. These findings suggest that differences in the complement of GABA neuron subtypes and in gene expression levels per neuron contribute to the specialization of inhibitory neurotransmission across cortical circuits.

Transcriptome Response and Spatial Pattern of Gene Expression in the Primate Subventricular Zone Neurogenic Niche After Cerebral Ischemia

The main stem cell niche for neurogenesis in the adult mammalian brain is the subventricular zone (SVZ) that extends along the cerebral lateral ventricles. We aimed at characterizing the initial molecular responses of the macaque monkey SVZ to transient, global cerebral ischemia. We microdissected tissue lining the anterior horn of the lateral ventricle (SVZa) from 7 day post-ischemic and sham-operated monkeys. Transcriptomics shows that in ischemic SVZa, 541 genes were upregulated and 488 genes were down-regulated. The transcription data encompassing the upregulated genes revealed a profile typical for quiescent stem cells and astrocytes. In the primate brain the SVZ is morphologically subdivided in distinct and separate ependymal and subependymal regions. The subependymal contains predominantly neural stem cells (NSC) and differentiated progenitors. To determine in which SVZa region ischemia had evoked transcriptional upregulation, sections through control and ischemic SVZa were analyzed by high-throughput in situ hybridization for a total of 150 upregulated genes shown in the www.monkey-niche.org image database. The majority of the differentially expressed genes mapped to the subependymal layers on the striatal or callosal aspect of the SVZa. Moreover, a substantial number of upregulated genes was expressed in the ependymal layer, implicating a contribution of the ependyma to stem cell biology. The transcriptome analysis yielded several novel gene markers for primate SVZa including the apelin receptor that is strongly expressed in the primate SVZa niche upon ischemic insult.

Perceptual rivalry across animal species

This review in memoriam of Jack Pettigrew provides an overview of past and current research into the phenomenon of multistable perception across multiple animal species. Multistable perception is characterized by two or more perceptual interpretations spontaneously alternating, or rivaling, when animals are exposed to stimuli with inherent sensory ambiguity. There is a wide array of ambiguous stimuli across sensory modalities, ranging from the configural changes observed in simple line drawings, such as the famous Necker cube, to the alternating perception of entire visual scenes that can be instigated by interocular conflict. The latter phenomenon, called binocular rivalry, in particular caught the attention of the late Jack Pettigrew, who combined his interest in the neuronal basis of perception with a unique comparative biological approach that considered ambiguous sensation as a fundamental problem of sensory systems that has shaped the brain throughout evolution. Here, we examine the research findings on visual perceptual alternation and suppression in a wide variety of species including insects, fish, reptiles, and primates. We highlight several interesting commonalities across species and behavioral indicators of perceptual alternation. In addition, we show how the comparative approach provides new avenues for understanding how the brain suppresses opposing sensory signals and generates alternations in perceptual dominance.

Non-Human Primate Blood-Brain Barrier and In Vitro Brain Endothelium: From Transcriptome to the Establishment of a New Model

The non-human primate (NHP)-brain endothelium constitutes an essential alternative to human in the prediction of molecule trafficking across the blood–brain barrier (BBB). This study presents a comparison between the NHP transcriptome of freshly isolated brain microcapillaries and in vitro-selected brain endothelial cells (BECs), focusing on important BBB features, namely tight junctions, receptors mediating transcytosis (RMT), ABC and SLC transporters, given its relevance as an alternative model for the molecule trafficking prediction across the BBB and identification of new brain-specific transport mechanisms. In vitro BECs conserved most of the BBB key elements for barrier integrity and control of molecular trafficking. The function of RMT via the transferrin receptor (TFRC) was characterized in this NHP-BBB model, where both human transferrin and anti-hTFRC antibody showed increased apical-to-basolateral passage in comparison to control molecules. In parallel, eventual BBB-related regional differences were Investig.igated in seven-day in vitro-selected BECs from five brain structures: brainstem, cerebellum, cortex, hippocampus, and striatum. Our analysis retrieved few differences in the brain endothelium across brain regions, suggesting a rather homogeneous BBB function across the brain parenchyma. The presently established NHP-derived BBB model closely mimics the physiological BBB, thus representing a ready-to-use tool for assessment of the penetration of biotherapeutics into the human CNS.

Closing the gap from transcription to the structural connectome enhances the study of connections in the human brain

The brain is composed of a complex web of networks but we have yet to map the structural connections of the human brain in detail. Diffusion MR imaging is a high-throughput method that relies on the principle of diffusion to reconstruct tracts (ie, pathways) across the brain. Although diffusion MR tractography is an exciting method to explore the structural connectivity of the brain in development and across species, the tractography has at times led to questionable interpretations. There are at present few if any alternative methods to trace structural pathways in the human brain. Given these limitations and the potential of diffusion MR imaging to map the human connectome, it is imperative that we develop new approaches to validate neuroimaging techniques. I discuss our recent studies integrating neuroimaging with transcriptional and anatomical variation across humans and other species over the course of development and in adulthood. Developing a novel framework to harness the potential of diffusion MR tractography provides new and exciting opportunities to study the evolution of developmental mechanisms generating variation in connections and bridge the gap between model systems to humans.

Loss of Nuclear TDP-43 Is Associated with Decondensation of LINE Retrotransposons

Loss of the nuclear RNA binding protein TAR DNA binding protein-43 (TDP-43) into cytoplasmic aggregates is the strongest correlate to neurodegeneration in amyotrophic lateral sclerosis and frontotemporal degeneration. The molecular changes associated with the loss of nuclear TDP-43 in human tissues are not entirely known. Using subcellular fractionation andfluorescent-activated cell sorting to enrich for diseased neuronal nuclei without TDP-43 from post-mortem frontotemporal degeneration-amyotro-phic lateral sclerosis (FTD-ALS) human brain, we characterized the effects of TDP-43 loss on the transcriptome and chromatin accessibility. Nuclear TDP-43 loss is associated with gene expression changes that affect RNA processing, nucleocytoplas-mic transport, histone processing, and DNA damage. Loss of nuclear TDP-43 is also associated with chromatin decondensation around long interspersed nuclear elements (LINEs) and increased LINE1 DNA content. Moreover, loss of TDP-43 leads to increased retrotransposition that can be inhibited with antiretro-viral drugs, suggesting that TDP-43 neuropathology is associated with altered chromatin structure including decondensation of LINEs.

Liu et al. fractionated and sorted for diseased neuronal nuclei from post-mortem FTD-ALS human brains and showed that loss of an RNA-binding protein, TDP-43, altered the transcriptome and chromatin accessibility. Their results suggest that loss of nuclear TDP-43 is associated with decondensation of LINE retrotransposons.

Genome-Wide Analysis Identifies NURR1-Controlled Network of New Synapse Formation and Cell Cycle Arrest in Human Neural Stem Cells

Nuclear receptor-related 1 (Nurr1) protein has been identified as an obligatory transcription factor in midbrain dopaminergic neurogenesis, but the global set of human NURR1 target genes remains unexplored. Here, we identified direct gene targets of NURR1 by analyzing genome-wide differential expression of NURR1 together with NURR1 consensus sites in three human neural stem cell (hNSC) lines. Microarray data were validated by quantitative PCR in hNSCs and mouse embryonic brains and through comparison to published human data, including genome-wide association study hits and the BioGPS gene expression atlas. Our analysis identified ~40 NURR1 direct target genes, many of them involved in essential protein modules such as synapse formation, neuronal cell migration during brain development, and cell cycle progression and DNA replication. Specifically, expression of genes related to synapse formation and neuronal cell migration correlated tightly with NURR1 expression, whereas cell cycle progression correlated negatively with it, precisely recapitulating midbrain dopaminergic development. Overall, this systematic examination of NURR1-controlled regulatory networks provides important insights into this protein's biological functions in dopamine-based neurogenesis.

FACS-array-based cell purification yields a specific transcriptome of striatal medium spiny neurons in a murine Huntington disease model

Huntington disease (HD) is a neurodegenerative disorder caused by expanded CAG repeats in the Huntingtin gene. Results from previous studies have suggested that transcriptional dysregulation is one of the key mechanisms underlying striatal medium spiny neuron (MSN) degeneration in HD. However, some of the critical genes involved in HD etiology or pathology could be masked in a common expression profiling assay because of contamination with non-MSN cells. To gain insight into the MSN-specific gene expression changes in presymptomatic R6/2 mice, a common HD mouse model, here we used a transgenic fluorescent protein marker of MSNs for purification via FACS before profiling gene expression with gene microarrays and compared the results of this "FACS-array" with those obtained with homogenized striatal samples (STR-array). We identified hundreds of differentially expressed genes (DEGs) and enhanced detection of MSN-specific DEGs by comparing the results of the FACS-array with those of the STR-array. The gene sets obtained included genes ubiquitously expressed in both MSNs and non-MSN cells of the brain and associated with transcriptional regulation and DNA damage responses. We proposed that the comparative gene expression approach using the FACS-array may be useful for uncovering the gene cascades affected in MSNs during HD pathogenesis.

Measurements of neuronal morphological variation across the rat neocortex

Neuron morphology is highly variable across the mammalian brain. It is thought that these attributes of neuronal cell shape, such as soma surface area and branching frequency, are determined by biological function and information processing. In this study, a large data set of neurons across the rat neocortex were clustered by their anatomical characters for evidence of distinctiveness among neocortical regions and the somatosensory layers. This data set of neuronal morphologies was compiled from 31 different lab sources with a validation procedure so that data records are potentially comparable across research studies. With this large set of heterogeneous data and by clustering analysis, this study shows that neuronal morphological traits overlap among neocortical and somatosensory regions. In the context of past neuroanatomical studies, this result is not congruent with tissue level analysis and strongly suggests further sampling of neuronal data to lessen the effect of confounding factors, such as the influence of different methodologies from use of heterogeneous samples of neuronal data.

Genetic Contributions to Multivariate Data-Driven Brain Networks Constructed via Source-Based Morphometry

Identifying genetic factors underlying neuroanatomical variation has been difficult. Traditional methods have used brain regions from predetermined parcellation schemes as phenotypes for genetic analyses, although these parcellations often do not reflect brain function and/or do not account for covariance between regions. We proposed that network-based phenotypes derived via source-based morphometry (SBM) may provide additional insight into the genetic architecture of neuroanatomy given its data-driven approach and consideration of covariance between voxels. We found that anatomical SBM networks constructed on ~ 20 000 individuals from the UK Biobank were heritable and shared functionally meaningful genetic overlap with each other. We additionally identified 27 unique genetic loci that contributed to one or more SBM networks. Both GWA and genetic correlation results indicated complex patterns of pleiotropy and polygenicity similar to other complex traits. Lastly, we found genetic overlap between a network related to the default mode and schizophrenia, a disorder commonly associated with neuroanatomic alterations.