André M.M. Sousa

Spatio-temporal transcriptome of the human brain

Brain development and function depend on the precise regulation of gene expression. However, our understanding of the complexity and dynamics of the transcriptome of the human brain is incomplete. Here we report the generation and analysis of exon-level transcriptome and associated genotyping data, representing males and females of different ethnicities, from multiple brain regions and neocortical areas of developing and adult post-mortem human brains. We found that 86 per cent of the genes analysed were expressed, and that 90 per cent of these were differentially regulated at the whole-transcript or exon level across brain regions and/or time. The majority of these spatio-temporal differences were detected before birth, with subsequent increases in the similarity among regional transcriptomes. The transcriptome is organized into distinct co-expression networks, and shows sex-biased gene expression and exon usage. We also profiled trajectories of genes associated with neurobiological categories and diseases, and identified associations between single nucleotide polymorphisms and gene expression. This study provides a comprehensive data set on the human brain transcriptome and insights into the transcriptional foundations of human neurodevelopment.

Prolonged myelination in human neocortical evolution

Nerve myelination facilitates saltatory action potential conduction and exhibits spatiotemporal variation during development associated with the acquisition of behavioral and cognitive maturity. Although human cognitive development is unique, it is not known whether the ontogenetic progression of myelination in the human neocortex is evolutionarily exceptional. In this study, we quantified myelinated axon fiber length density and the expression of myelin-related proteins throughout postnatal life in the somatosensory (areas 3b/3a/1/2), motor (area 4), frontopolar (prefrontal area 10), and visual (areas 17/18) neocortex of chimpanzees (N = 20) and humans (N = 33). Our examination revealed that neocortical myelination is developmentally protracted in humans compared with chimpanzees. In chimpanzees, the density of myelinated axons increased steadily until adult-like levels were achieved at approximately the time of sexual maturity. In contrast, humans displayed slower myelination during childhood, characterized by a delayed period of maturation that extended beyond late adolescence. This comparative research contributes evidence crucial to understanding the evolution of human cognition and behavior, which arises from the unfolding of nervous system development within the context of an enriched cultural environment. Perturbations of normal developmental processes and the decreased expression of myelin-related molecules have been related to psychiatric disorders such as schizophrenia. Thus, these species differences suggest that the human-specific shift in the timing of cortical maturation during adolescence may have implications for vulnerability to certain psychiatric disorders.

Zika Virus Disrupts Phospho-TBK1 Localization and Mitosis in Human Neuroepithelial Stem Cells and Radial Glia.

The mechanisms underlying Zika virus (ZIKV)-related microcephaly and other neurodevelopment defects remain poorly understood. Here, we describe the derivation and characterization, including single-cell RNA-seq, of neocortical and spinal cord neuroepithelial stem (NES) cells to model early human neurodevelopment and ZIKV-related neuropathogenesis. By analyzing human NES cells, organotypic fetal brain slices, and a ZIKV-infected micrencephalic brain, we show that ZIKV infects both neocortical and spinal NES cells as well as their fetal homolog, radial glial cells (RGCs), causing disrupted mitoses, supernumerary centrosomes, structural disorganization, and cell death. ZIKV infection of NES cells and RGCs causes centrosomal depletion and mitochondrial sequestration of phospho-TBK1 during mitosis. We also found that nucleoside analogs inhibit ZIKV replication in NES cells, protecting them from ZIKV-induced pTBK1 relocalization and cell death. We established a model system of human neural stem cells to reveal cellular and molecular mechanisms underlying neurodevelopmental defects associated with ZIKV infection and its potential treatment.

Temporal Specification and Bilaterality of Human Neocortical Topographic Gene Expression

Transcriptional events involved in the development of human cerebral neocortex are poorly understood. Here, we analyzed the temporal dynamics and laterality of gene expression in human and macaque monkey neocortex. We found that interareal differences exhibit a temporal hourglass pattern, dividing the human neocortical development into three major phases. The first phase, corresponding to prenatal development, is characterized by the highest number of differential expressed genes among areas and gradient-like expression patterns, including those that are different between human and macaque. The second, preadolescent phase, is characterized by lesser interareal expression differences and by an increased synchronization of areal transcriptomes. During the third phase, from adolescence onward, differential expression among areas increases again driven predominantly by a subset of areas, without obvious gradient-like patterns. Analyses of left-right gene expression revealed population-level global symmetry throughout the fetal and postnatal time span. Thus, human neocortical topographic gene expression is temporally specified and globally symmetric.

Species-Dependent Posttranscriptional Regulation of NOS1 by FMRP in the Developing Cerebral Cortex

Fragile X syndrome (FXS), the leading monogenic cause of intellectual disability and autism, results from loss of function of the RNA-binding protein FMRP. Here, we show that FMRP regulates translation of neuronal nitric oxide synthase 1 (NOS1) in the developing human neocortex. Whereas NOS1 mRNA is widely expressed, NOS1 protein is transiently coexpressed with FMRP during early synaptogenesis in layer- and region-specific pyramidal neurons. These include midfetal layer 5 subcortically projecting neurons arranged into alternating columns in the prospective Brocas area and orofacial motor cortex. Human NOS1 translation is activated by FMRP via interactions with coding region binding motifs absent from mouse Nos1 mRNA, which is expressed in mouse pyramidal neurons, but not efficiently translated. Correspondingly, neocortical NOS1 protein levels are severely reduced in developing human FXS cases, but not FMRP-deficient mice. Thus, alterations in FMRP posttranscriptional regulation of NOS1 in developing neocortical circuits may contribute to cognitive dysfunction in FXS.

Spatially Heterogeneous Choroid Plexus Transcriptomes Encode Positional Identity and Contribute to Regional CSF Production

A sheet of choroid plexus epithelial cells extends into each cerebral ventricle and secretes signaling factors into the CSF. To evaluate whether differences in the CSF proteome across ventricles arise, in part, from regional differences in choroid plexus gene expression, we defined the transcriptome of lateral ventricle (telencephalic) versus fourth ventricle (hindbrain) choroid plexus. We find that positional identities of mouse, macaque, and human choroid plexi derive from gene expression domains that parallel their axial tissues of origin. We then show that molecular heterogeneity between telencephalic and hindbrain choroid plexi contributes to region-specific, age-dependent protein secretion in vitro. Transcriptome analysis of FACS-purified choroid plexus epithelial cells also predicts their cell-type-specific secretome. Spatial domains with distinct protein expression profiles were observed within each choroid plexus. We propose that regional differences between choroid plexi contribute to dynamic signaling gradients across the mammalian cerebroventricular system.

Molecular and cellular reorganization of neural circuits in the human lineage

The makings of the primate brain Although nonhuman primate brains are similar to our own, the disparity between their and our cognitive abilities tells us that surface similarity is not the whole story. Sousa et al. overlaid transcriptome and histological analyses to see what makes human brains different from those of nonhuman primates. Various differentially expressed genes, such as those encoding transcription factors, could alter transcriptional programs. Others were associated with neuromodulatory systems. Furthermore, the dopaminergic interneurons found in the human neocortex were absent from the neocortex of nonhuman African apes. Such differences in neuronal transcriptional programs may underlie a variety of neurodevelopmental disorders. Science, this issue p. 1027 Comparing transcriptome and histology of human and nonhuman primate brains reveals changes that make humans unique. To better understand the molecular and cellular differences in brain organization between human and nonhuman primates, we performed transcriptome sequencing of 16 regions of adult human, chimpanzee, and macaque brains. Integration with human single-cell transcriptomic data revealed global, regional, and cell-type–specific species expression differences in genes representing distinct functional categories. We validated and further characterized the human specificity of genes enriched in distinct cell types through histological and functional analyses, including rare subpallial-derived interneurons expressing dopamine biosynthesis genes enriched in the human striatum and absent in the nonhuman African ape neocortex. Our integrated analysis of the generated data revealed diverse molecular and cellular features of the phylogenetic reorganization of the human brain across multiple levels, with relevance for brain function and disease.

A multiregional proteomic survey of the postnatal human brain

Detailed observations of transcriptional, translational and post-translational events in the human brain are essential to improving our understanding of its development, function and vulnerability to disease. Here, we exploited label-free quantitative tandem mass-spectrometry to create an in-depth proteomic survey of regions of the postnatal human brain, ranging in age from early infancy to adulthood. Integration of protein data with existing matched whole-transcriptome sequencing (RNA-seq) from the BrainSpan project revealed varied patterns of protein–RNA relationships, with generally increased magnitudes of protein abundance differences between brain regions compared to RNA. Many of the differences amplified in protein data were reflective of cytoarchitectural and functional variation between brain regions. Comparing structurally similar cortical regions revealed significant differences in the abundances of receptor-associated and resident plasma membrane proteins that were not readily observed in the RNA expression data.Quantitative mass spectrometry was used to produce a proteomic survey of postnatal human brain regions. Compared to matched RNA-seq, protein levels showed more regional variation, especially for membrane-associated proteins in the neocortex.

XSAnno: a framework for building ortholog models in cross-species transcriptome comparisons

BackgroundThe accurate characterization of RNA transcripts and expression levels across species is critical for understanding transcriptome evolution. As available RNA-seq data accumulate rapidly, there is a great demand for tools that build gene annotations for cross-species RNA-seq analysis. However, prevailing methods of ortholog annotation for RNA-seq analysis between closely-related species do not take inter-species variation in mappability into consideration.ResultsHere we present XSAnno, a computational framework that integrates previous approaches with multiple filters to improve the accuracy of inter-species transcriptome comparisons. The implementation of this approach in comparing RNA-seq data of human, chimpanzee, and rhesus macaque brain transcriptomes has reduced the false discovery of differentially expressed genes, while maintaining a low false negative rate.ConclusionThe present study demonstrates the utility of the XSAnno pipeline in building ortholog annotations and improving the accuracy of cross-species transcriptome comparisons.

Planar cell polarity pathway and development of the human visual cortex

The radial unit hypothesis provides a framework for global (proliferation) and regional (distribution) expansion of the primate cerebral cortex1. Using principal component analysis (PCA), we have identified cortical regions with shared variance in their surface area and cortical thickness, respectively, segmented from magnetic resonance images obtained in 19,171 participants. We then carried out meta-analyses of genome-wide association studies of the first two principal components for each phenotype. For surface area (but not cortical thickness), we have detected strong associations between each of the components and single nucleotide polymorphisms in a number of gene loci. The first (“global”) component was associated mainly with loci on chromosome 17 (8.8×10 -26 ≤ p ≤ 2.3×10-14), including those detected previously as linked with intracranial volume2,3. The second (regional) component captured shared variation in the surface area of the primary and adjacent secondary visual cortices and showed a robust association with polymorphisms in a locus on chromosome 14 containing Disheveled Associated Activator of Morphogenesis 1 (DAAM1; p=3.0×10-32). DAAM1 is a key component in the planar-cell-polarity signaling pathway4,5. In follow-up studies, we have focused on the latter finding and established that: (1) DAAM1 is highly expressed between 12th and 22nd post-conception weeks in the human cerebral cortex; (2) genes co-expressed with DAAM1 in the primary visual cortex are enriched in mitochondria-related pathways; and (3) volume of the lateral geniculate nucleus, which projects to regions of the visual cortex staining for cytochrome oxidase (a mitochondrial enzyme), correlates with the surface area of the visual cortex in major-allele homozygotes but not in carriers of the minor allele. Altogether, we speculate that - in concert with thalamocortical input to cortical subplate - DAAM1 enables migration of neurons to cytochrome-oxidase rich regions of the visual cortex, and, in turn, facilitates regional expansion of this set of cortical regions during development.