Laminar and temporal expression dynamics of coding and noncoding RNAs in the mouse neocortex

Transcriptional determinism and stochasticity contribute to the complexity of autism-associated SHANK family genes

Precision of transcription is critical because transcriptional dysregulation is disease causing. Traditional methods of transcriptional profiling are inadequate to elucidate the full spectrum of the transcriptome, particularly for longer and less abundant mRNAs. SHANK3 is one of the most common autism causative genes. Twenty-four Shank3-mutant animal lines have been developed for autism modeling. However, their preclinical validity has been questioned due to incomplete Shank3 transcript structure. We apply an integrative approach combining cDNA-capture and long-read sequencing to profile the SHANK3 transcriptome in humans and mice. We unexpectedly discover an extremely complex SHANK3 transcriptome. Specific SHANK3 transcripts are altered in Shank3-mutant mice and postmortem brain tissues from individuals with autism spectrum disorder. The enhanced SHANK3 transcriptome significantly improves the detection rate for potential deleterious variants from genomics studies of neuropsychiatric disorders. Our findings suggest that both deterministic and stochastic transcription of the genome is associated with SHANK family genes.

Region and layer-specific expression of GABAA receptor isoforms and KCC2 in developing cortex

Introduction

γ-Aminobutyric acid (GABA) type A receptors (GABAARs) are ligand-gated Cl-channels that mediate the bulk of inhibitory neurotransmission in the mature CNS and are targets of many drugs. During cortical development, GABAAR-mediated signals are significantly modulated by changing subunit composition and expression of Cl-transporters as part of developmental processes and early network activity. To date, this developmental evolution has remained understudied, particularly at the level of cortical layer-specific changes. In this study, we characterized the expression of nine major GABAAR subunits and K-Cl transporter 2 (KCC2) in mouse somatosensory cortex from embryonic development to postweaning maturity.

Methods

We evaluated expression of α1-5, β2-3, γ2, and δ GABAAR subunits using immunohistochemistry and Western blot techniques, and expression of KCC2 using immunohistochemistry in cortices from E13.5 to P25 mice.

Results

We found that embryonic cortex expresses mainly α3, α5, β3, and γ2, while expression of α1, α2, α4, β2, δ, and KCC2 begins at later points in development; however, many patterns of nuanced expression can be found in specific lamina, cortical regions, and cells and structures.

Discussion

While the general pattern of expression of each subunit and KCC2 is similar to previous studies, we found a number of unique temporal, regional, and laminar patterns that were previously unknown. These findings provide much needed knowledge of the intricate developmental evolution in GABAAR composition and KCC2 expression to accommodate developmental signals that transition to mature neurotransmission.

Alteration of synaptic protein composition during developmental synapse maturation

The postsynaptic density (PSD) is a collection of specialized proteins assembled beneath the postsynaptic membrane of dendritic spines. The PSD proteome comprises ~1000 proteins, including neurotransmitter receptors, scaffolding proteins and signalling enzymes. Many of these proteins have essential roles in synaptic function and plasticity. During brain development, changes are observed in synapse density and in the stability and shape of spines, reflecting the underlying molecular maturation of synapses. Synaptic protein composition changes in terms of protein abundance and the assembly of protein complexes, supercomplexes and the physical organization of the PSD. Here, we summarize the developmental alterations of postsynaptic protein composition during synapse maturation. We describe major PSD proteins involved in postsynaptic signalling that regulates synaptic plasticity and discuss the effect of altered expression of these proteins during development. We consider the abnormality of synaptic profiles and synaptic protein composition in the brain in neurodevelopmental disorders such as autism spectrum disorders. We also explain differences in synapse development between rodents and primates in terms of synaptic profiles and protein composition. Finally, we introduce recent findings related to synaptic diversity and nanoarchitecture and discuss their impact on future research. Synaptic protein composition can be considered a major determinant and marker of synapse maturation in normality and disease.

A data-driven single-cell and spatial transcriptomic map of the human prefrontal cortex

The molecular organization of the human neocortex historically has been studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally defined spatial domains that move beyond classic cytoarchitecture. We used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex. Integration with paired single-nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we mapped the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains.

Remodeling of the postsynaptic proteome in male mice and marmosets during synapse development

Postsynaptic proteins play crucial roles in synaptic function and plasticity. During brain development, alterations in synaptic number, shape, and stability occur, known as synapse maturation. However, the postsynaptic protein composition changes during development are not fully understood. Here, we show the trajectory of the postsynaptic proteome in developing male mice and common marmosets. Proteomic analysis of mice at 2, 3, 6, and 12 weeks of age shows that proteins involved in synaptogenesis are differentially expressed during this period. Analysis of published transcriptome datasets shows that the changes in postsynaptic protein composition in the mouse brain after 2 weeks of age correlate with gene expression changes. Proteomic analysis of marmosets at 0, 2, 3, 6, and 24 months of age show that the changes in the marmoset brain can be categorized into two parts: the first 2 months and after that. The changes observed in the first 2 months are similar to those in the mouse brain between 2 and 12 weeks of age. The changes observed in marmoset after 2 months old include differential expression of synaptogenesis-related molecules, which hardly overlap with that in mice. Our results provide a comprehensive proteomic resource that underlies developmental synapse maturation in rodents and primates.

Transcriptional Determinism and Stochasticity Contribute to the Complexity of Autism Associated SHANK Family Genes

Precision of transcription is critical because transcriptional dysregulation is disease causing. Traditional methods of transcriptional profiling are inadequate to elucidate the full spectrum of the transcriptome, particularly for longer and less abundant mRNAs. SHANK3 is one of the most common autism causative genes. Twenty-four Shank3 mutant animal lines have been developed for autism modeling. However, their preclinical validity has been questioned due to incomplete Shank3 transcript structure. We applied an integrative approach combining cDNA-capture and long-read sequencing to profile the SHANK3 transcriptome in human and mice. We unexpectedly discovered an extremely complex SHANK3 transcriptome. Specific SHANK3 transcripts were altered in Shank3 mutant mice and postmortem brains tissues from individuals with ASD. The enhanced SHANK3 transcriptome significantly improved the detection rate for potential deleterious variants from genomics studies of neuropsychiatric disorders. Our findings suggest the stochastic transcription of genome associated with SHANK family genes.

MicroRNAs are necessary for the emergence of Purkinje cell identity

Brain computations are dictated by the unique morphology and connectivity of neuronal subtypes, features established by closely timed developmental events. MicroRNAs (miRNAs) are critical for brain development, but current technologies lack the spatiotemporal resolution to determine how miRNAs instruct the steps leading to subtype identity. Here, we developed new tools to tackle this major gap. Fast and reversible miRNA loss-of-function revealed that miRNAs are necessary for cerebellar Purkinje cell (PC) differentiation, which previously appeared miRNA-independent, and resolved distinct miRNA critical windows in PC dendritogenesis and climbing fiber synaptogenesis, key determinants of PC identity. To identify underlying mechanisms, we generated a mouse model, which enables precise mapping of miRNAs and their targets in rare cell types. With PC-specific maps, we found that the PC-enriched miR-206 drives exuberant dendritogenesis and modulates synaptogenesis. Our results showcase vastly improved approaches for dissecting miRNA function and reveal that many critical miRNA mechanisms remain largely unexplored.

Highlights

Fast miRNA loss-of-function with T6B impairs postnatal Purkinje cell developmentReversible T6B reveals critical miRNA windows for dendritogenesis and synaptogenesisConditional Spy3-Ago2 mouse line enables miRNA-target network mapping in rare cellsPurkinje cell-enriched miR-206 regulates its unique dendritic and synaptic morphology.

Gut-derived peptide hormone receptor expression in the developing mouse hypothalamus

OBJECTIVE

In adult organisms, a number of receptors have been identified which modulate metabolic processes related to peptides derived from the intestinal tract. These receptors play significant roles in glucose homeostasis, food intake and energy balance. Here we assess these classical metabolic receptors and their expression as well as their potential role in early development of hypothalamic neuronal circuits.

METHODS

Chow-fed C57BL6/N female mice were mated and hypothalamic tissue was collected from offspring across postnatal development (postnatal day 7-21). Subsequent qPCR and Western Blot analyses were used to determine mRNA and protein changes in gut-derived peptide hormone receptors. Correlations to body weight, blood glucose and circulating leptin levels were analyzed.

RESULTS

We describe the gene expression and dynamic protein regulation of key gut-derived peptide hormone receptors in the early postnatal period of the mouse brain. Specifically, we show changes to Gastric inhibitory polypeptide receptor (GIPR), glucagon-like peptide 1 receptor (GLP1R), and cholecystokinin receptor 2 (CCK2R) in the developing hypothalamus. The changes to GIPR and InsR seem to be strongly negatively correlated with body weight.

CONCLUSIONS

This comprehensive analysis underscores the need to understand the roles of maternal-derived circulating gut hormones and their direct effect on offspring brain development.

Integrated single cell and unsupervised spatial transcriptomic analysis defines molecular anatomy of the human dorsolateral prefrontal cortex

Generation of a molecular neuroanatomical map of the human prefrontal cortex reveals novel spatial domains and cell-cell interactions relevant for psychiatric disease. The molecular organization of the human neocortex has been historically studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally-defined spatial domains that move beyond classic cytoarchitecture. Here we used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex (DLPFC). Integration with paired single nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we map the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains. Finally, we provide resources for the scientific community to explore these integrated spatial and single cell datasets at research.libd.org/spatialDLPFC/.

MicroRNA 3' ends shorten during adolescent brain maturation

MicroRNA (miRNA) dysregulation is well-documented in psychiatric disease, but miRNA dynamics remain poorly understood during adolescent and early adult brain maturation, when symptoms often first appear. Here, we use RNA sequencing to examine miRNAs and their mRNA targets in cortex and hippocampus from early-, mid-, and late-adolescent and adult mice. Furthermore, we use quantitative proteomics by tandem mass tag mass spectrometry (TMT-MS) to examine protein dynamics in cortex from the same subjects. We found that ~25% of miRNAs' 3' ends shorten with age due to increased 3' trimming and decreased U tailing. Particularly, shorter but functionally competent isoforms (isomiRs) of miR-338-3p increase up to 10-fold during adolescence and only in brain. MiRNAs that undergo 3' shortening exhibit stronger negative correlations with targets that decrease with age and stronger positive correlations with targets that increase with age, than miRNAs with stable 3' ends. Increased 3' shortening with age was also observed in available mouse and human miRNA-seq data sets, and stronger correlations between miRNAs that undergo shortening and their mRNA targets were observed in two of the three available data sets. We conclude that age-associated miRNA 3' shortening is a well-conserved feature of postnatal brain maturation.

Molecular and cellular evolution of the primate dorsolateral prefrontal cortex

The granular dorsolateral prefrontal cortex (dlPFC) is an evolutionary specialization of primates that is centrally involved in cognition. We assessed more than 600,000 single-nucleus transcriptomes from adult human, chimpanzee, macaque, and marmoset dlPFC. Although most cell subtypes defined transcriptomically are conserved, we detected several that exist only in a subset of species as well as substantial species-specific molecular differences across homologous neuronal, glial, and non-neural subtypes. The latter are exemplified by human-specific switching between expression of the neuropeptide somatostatin and tyrosine hydroxylase, the rate-limiting enzyme in dopamine production in certain interneurons. The above molecular differences are also illustrated by expression of the neuropsychiatric risk gene FOXP2, which is human-specific in microglia and primate-specific in layer 4 granular neurons. We generated a comprehensive survey of the dlPFC cellular repertoire and its shared and divergent features in anthropoid primates.

Identification of a transcriptional signature found in multiple models of ASD and related disorders

Epigenetic regulation plays a critical role in many neurodevelopmental disorders (NDDs), including autism spectrum disorder (ASD). In particular, many such disorders are the result of mutations in genes that encode chromatin-modifying proteins. However, although these disorders share many features, it is unclear whether they also share gene expression disruptions resulting from the aberrant regulation of chromatin. We examined five chromatin modifiers that are all linked to ASD despite their different roles in regulating chromatin. Specifically, we depleted ASH1L, CHD8, CREBBP, EHMT1, and NSD1 in parallel in a highly controlled neuronal culture system. We then identified sets of shared genes, or transcriptional signatures, that are differentially expressed following loss of multiple ASD-linked chromatin modifiers. We examined the functions of genes within the transcriptional signatures and found an enrichment in many neurotransmitter transport genes and activity-dependent genes. In addition, these genes are enriched for specific chromatin features such as bivalent domains that allow for highly dynamic regulation of gene expression. The down-regulated transcriptional signature is also observed within multiple mouse models of NDDs that result in ASD, but not those only associated with intellectual disability. Finally, the down-regulated transcriptional signature can distinguish between control and idiopathic ASD patient iPSC-derived neurons as well as postmortem tissue, demonstrating that this gene set is relevant to the human disorder. This work identifies a transcriptional signature that is found within many neurodevelopmental syndromes, helping to elucidate the link between epigenetic regulation and the underlying cellular mechanisms that result in ASD.

Mechanisms underlying the role of ankyrin-B in cardiac and neurological health and disease

The ANK2 gene encodes for ankyrin-B (ANKB), one of 3 members of the ankyrin family of proteins, whose name is derived from the Greek word for anchor. ANKB was originally identified in the brain (B denotes “brain”) but has become most widely known for its role in cardiomyocytes as a scaffolding protein for ion channels and transporters, as well as an interacting protein for structural and signaling proteins. Certain loss-of-function ANK2 variants are associated with a primarily cardiac-presenting autosomal-dominant condition with incomplete penetrance and variable expressivity characterized by a predisposition to supraventricular and ventricular arrhythmias, arrhythmogenic cardiomyopathy, congenital and adult-onset structural heart disease, and sudden death. Another independent group of ANK2 variants are associated with increased risk for distinct neurological phenotypes, including epilepsy and autism spectrum disorders. The mechanisms underlying ANKB's roles in cells in health and disease are not fully understood; however, several clues from a range of molecular and cell biological studies have emerged. Notably, ANKB exhibits several isoforms that have different cell-type–, tissue–, and developmental stage– expression profiles. Given the conservation within ankyrins across evolution, model organism studies have enabled the discovery of several ankyrin roles that could shed important light on ANKB protein-protein interactions in heart and brain cells related to the regulation of cellular polarity, organization, calcium homeostasis, and glucose and fat metabolism. Along with this accumulation of evidence suggesting a diversity of important ANKB cellular functions, there is an on-going debate on the role of ANKB in disease. We currently have limited understanding of how these cellular functions link to disease risk. To this end, this review will examine evidence for the cellular roles of ANKB and the potential contribution of ANKB functional variants to disease risk and presentation. This contribution will highlight the impact of ANKB dysfunction on cardiac and neuronal cells and the significance of understanding the role of ANKB variants in disease.

p53-mediated neurodegeneration in the absence of the nuclear protein Akirin2

Proper gene regulation is critical for both neuronal development and maintenance as the brain matures. We previously demonstrated that Akirin2, an essential nuclear protein that interacts with transcription factors and chromatin remodeling complexes, is required for the embryonic formation of the cerebral cortex. Here we show that Akirin2 plays a mechanistically distinct role in maintaining healthy neurons during cortical maturation. Restricting Akirin2 loss to excitatory cortical neurons resulted in progressive neurodegeneration via necroptosis and severe cortical atrophy with age. Comparing transcriptomes from Akirin2-null postnatal neurons and cortical progenitors revealed that targets of the tumor suppressor p53, a regulator of both proliferation and cell death encoded by Trp53, were consistently upregulated. Reduction of Trp53 rescued neurodegeneration in Akirin2-null neurons. These data: (1) implicate Akirin2 as a critical neuronal maintenance protein, (2) identify p53 pathways as mediators of Akirin2 functions, and (3) suggest Akirin2 dysfunction may be relevant to neurodegenerative diseases.

A dominant tubulin mutation causes cerebellar neurodegeneration in a genetic model of tubulinopathy

Mutations in tubulins cause distinct neurodevelopmental and degenerative diseases termed "tubulinopathies"; however, little is known about the functional requirements of tubulins or how mutations cause cell-specific pathologies. Here, we identify a mutation in the gene Tubb4a that causes degeneration of cerebellar granule neurons and myelination defects. We show that the neural phenotypes result from a cell type-specific enrichment of a dominant mutant form of Tubb4a relative to the expression other β-tubulin isotypes. Loss of Tubb4a function does not underlie cellular pathology but is compensated by the transcriptional up-regulation of related tubulin genes in a cell type-specific manner. This work establishes that the expression of a primary tubulin mutation in mature neurons is sufficient to promote cell-autonomous cell death, consistent with a causative association of microtubule dysfunction with neurodegenerative diseases. These studies provide evidence that mutations in tubulins cause specific phenotypes based on expression ratios of tubulin isotype genes.

Identification of cell-type-specific marker genes from co-expression patterns in tissue samples

MOTIVATION

Marker genes, defined as genes that are expressed primarily in a single-cell type, can be identified from the single-cell transcriptome; however, such data are not always available for the many uses of marker genes, such as deconvolution of bulk tissue. Marker genes for a cell type, however, are highly correlated in bulk data, because their expression levels depend primarily on the proportion of that cell type in the samples. Therefore, when many tissue samples are analyzed, it is possible to identify these marker genes from the correlation pattern.

RESULTS

To capitalize on this pattern, we develop a new algorithm to detect marker genes by combining published information about likely marker genes with bulk transcriptome data in the form of a semi-supervised algorithm. The algorithm then exploits the correlation structure of the bulk data to refine the published marker genes by adding or removing genes from the list.

AVAILABILITY AND IMPLEMENTATION

We implement this method as an R package markerpen, hosted on CRAN (https://CRAN.R-project.org/package=markerpen).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

MicroRNAs in the Onset of Schizophrenia

Mounting evidence implicates microRNAs (miRNAs) in the pathology of schizophrenia. These small noncoding RNAs bind to mRNAs containing complementary sequences and promote their degradation and/or inhibit protein synthesis. A single miRNA may have hundreds of targets, and miRNA targets are overrepresented among schizophrenia-risk genes. Although schizophrenia is a neurodevelopmental disorder, symptoms usually do not appear until adolescence, and most patients do not receive a schizophrenia diagnosis until late adolescence or early adulthood. However, few studies have examined miRNAs during this critical period. First, we examine evidence that the miRNA pathway is dynamic throughout adolescence and adulthood and that miRNAs regulate processes critical to late neurodevelopment that are aberrant in patients with schizophrenia. Next, we examine evidence implicating miRNAs in the conversion to psychosis, including a schizophrenia-associated single nucleotide polymorphism in MIR137HG that is among the strongest known predictors of age of onset in patients with schizophrenia. Finally, we examine how hemizygosity for DGCR8, which encodes an obligate component of the complex that synthesizes miRNA precursors, may contribute to the onset of psychosis in patients with 22q11.2 microdeletions and how animal models of this disorder can help us understand the many roles of miRNAs in the onset of schizophrenia.

Do Fragile X Syndrome and Other Intellectual Disorders Converge at Aberrant Pre-mRNA Splicing?

Fragile X Syndrome is a neuro-developmental disorder caused by the silencing of the FMR1 gene, resulting in the loss of its protein product, FMRP. FMRP binds mRNA and represses general translation in the brain. Transcriptome analysis of the Fmr1-deficient mouse hippocampus reveals widespread dysregulation of alternative splicing of pre-mRNAs. Many of these aberrant splicing changes coincide with those found in post-mortem brain tissue from individuals with autism spectrum disorders (ASDs) as well as in mouse models of intellectual disability such as PTEN hamartoma syndrome (PHTS) and Rett Syndrome (RTT). These splicing changes could result from chromatin modifications (e.g., in FXS, RTT) and/or splicing factor alterations (e.g., PTEN, autism). Based on the identities of the RNAs that are mis-spliced in these disorders, it may be that they are at least partly responsible for some shared pathophysiological conditions. The convergence of splicing aberrations among these autism spectrum disorders might be crucial to understanding their underlying cognitive impairments.

Scaling of gene transcriptional gradients with brain size across mouse development

The structure of the adult brain is the result of complex physical mechanisms acting in three-dimensional space through development. Consequently, the brain's spatial embedding plays a key role in its organization, including the gradient-like patterning of gene expression that encodes the molecular underpinning of functional specialization. However, we do not yet understand how changes in brain shape and size that occur across development influence the brain's transcriptional architecture. Here we investigate the spatial embedding of transcriptional patterns of over 1800 genes across seven time points through mouse-brain development using data from the Allen Developing Mouse Brain Atlas. We find that transcriptional similarity decreases exponentially with separation distance across all developmental time points, with a correlation length scale that follows a power-law scaling relationship with a linear dimension of brain size. This scaling suggests that the mouse brain achieves a characteristic balance between local molecular similarity (homogeneous gene expression within a specialized brain area) and longer-range diversity (between functionally specialized brain areas) throughout its development. Extrapolating this mouse developmental scaling relationship to the human cortex yields a prediction consistent with the value measured from microarray data. We introduce a simple model of brain growth as spatially autocorrelated gene-expression gradients that expand through development, which captures key features of the mouse developmental data. Complementing the well-known exponential distance rule for structural connectivity, our findings characterize an analogous exponential distance rule for transcriptional gradients that scales across mouse brain development, providing new understanding of spatial constraints on the brain's molecular patterning.

Hypoxia-Ischemia Induced Age-Dependent Gene Transcription Effects at Two Development Stages in the Neonate Mouse Brain

Human brain lesions in the perinatal period result in life-long neuro-disabilities impairing sensory-motor, cognitive, and behavior functions for years. Topographical aspects of brain lesions depend on gestational age at the time of insult in preterm or term infants and impaired subsequent steps of brain development and maturation. In mice, the Rice-Vannucci procedure of neonate hypoxia-ischemia (HI) was used at 5 days (P5) or P10, mimicking the development of 30 week-gestation fetus/preterm newborn, or full-term infant, respectively. Transcription response to HI was assessed at 3, 6, 12, and 24 h after insult, using micro-array technology. Statistical Pathway and Gene Ontology terms enrichments were investigated using DAVID®, Revigo® and Ingenuity Pathway Analysis (IPA®) to identify a core of transcription response to HI, age-specific regulations, and interactions with spontaneous development. Investigations were based on direction, amplitude, and duration of responses, basal expression, and annotation. Five major points deserve attention; (i) inductions exceeded repressions (60/40%) at both ages, (ii) only 20.3% (393/1938 records) were common to P5 and P10 mice, (iii) at P5, HI effects occurred early and decreased 24 h after insult whereas they were delayed at P10 and increased 24 h after insult, (iv) common responses at P5 and P10 involved inflammation, immunity, apoptosis, and angiogenesis. (v) age-specific effects occurred with higher statistical significance at P5 than at P10. Transient repression of 12 genes encoding cholesterol biosynthesis enzymes was transiently observed 12 h after HI at P5. Synaptogenesis appeared inhibited at P5 while induced at P10, showing reciprocal effects on glutamate receptors. Specific involvement of Il-1 (interleukin-1) implicated in the firing of inflammation was observed at P10. This study pointed out age-differences in HI responses kinetics, e.g., a long-lasting inflammatory response at P10 compared to P5. Whether the specific strong depression of cholesterol biosynthesis genes that could account for white matter-specific vulnerability at P5 or prevent delayed inflammation needs further investigation. Determination of putative involvement of Il-1 and the identification of upstream regulators involved in the delayed inflammation firing at P10 appears promising routes of research in the understandings of age-dependent vulnerabilities in the neonatal brain.

Identification of de novo mutations in prenatal neurodevelopment-associated genes in schizophrenia in two Han Chinese patient-sibling family-based cohorts

Schizophrenia (SCZ) is a severe psychiatric disorder with a strong genetic component. High heritability of SCZ suggests a major role for transmitted genetic variants. Furthermore, SCZ is also associated with a marked reduction in fecundity, leading to the hypothesis that alleles with large effects on risk might often occur de novo. In this study, we conducted whole-genome sequencing for 23 families from two cohorts with unaffected siblings and parents. Two nonsense de novo mutations (DNMs) in GJC1 and HIST1H2AD were identified in SCZ patients. Ten genes (DPYSL2, NBPF1, SDK1, ZNF595, ZNF718, GCNT2, SNX9, AACS, KCNQ1, and MSI2) were found to carry more DNMs in SCZ patients than their unaffected siblings by burden test. Expression analyses indicated that these DNM implicated genes showed significantly higher expression in prefrontal cortex in prenatal stage. The DNM in the GJC1 gene is highly likely a loss function mutation (pLI = 0.94), leading to the dysregulation of ion channel in the glutamatergic excitatory neurons. Analysis of rare variants in independent exome sequencing dataset indicates that GJC1 has significantly more rare variants in SCZ patients than in unaffected controls. Data from genome-wide association studies suggested that common variants in the GJC1 gene may be associated with SCZ and SCZ-related traits. Genes co-expressed with GJC1 are involved in SCZ, SCZ-associated pathways, and drug targets. These evidences suggest that GJC1 may be a risk gene for SCZ and its function may be involved in prenatal and early neurodevelopment, a vulnerable period for developmental disorders such as SCZ.

Spatio-Temporal Roles of ASD-Associated Variants in Human Brain Development

Transcriptional regulation of the genome arguably provides the basis for the anatomical elaboration and dynamic operation of the human brain. It logically follows that genetic variations affecting gene transcription contribute to mental health disorders, including autism spectrum disorder (ASD). A number of recent studies have shown the role of de novo variants (DNVs) in disrupting early neurodevelopment. However, there is limited knowledge concerning the role of inherited variants during the early brain development of ASD. In this study, we investigate the role of rare inherited variations in neurodevelopment. We conducted co-expression network analyses using an anatomically comprehensive atlas of the developing human brain and examined whether rare coding and regulatory variants, identified from our genetic screening of Australian families with ASD, work in different spatio-temporal functions.

Application of Computational Biology to Decode Brain Transcriptomes

The rapid development of high-throughput sequencing technologies has generated massive valuable brain transcriptome atlases, providing great opportunities for systematically investigating gene expression characteristics across various brain regions throughout a series of developmental stages. Recent studies have revealed that the transcriptional architecture is the key to interpreting the molecular mechanisms of brain complexity. However, our knowledge of brain transcriptional characteristics remains very limited. With the immense efforts to generate high-quality brain transcriptome atlases, new computational approaches to analyze these high-dimensional multivariate data are greatly needed. In this review, we summarize some public resources for brain transcriptome atlases and discuss the general computational pipelines that are commonly used in this field, which would aid in making new discoveries in brain development and disorders.

Uncovering the Transcriptional Correlates of Hub Connectivity in Neural Networks

Connections in nervous systems are disproportionately concentrated on a small subset of neural elements that act as network hubs. Hubs have been found across different species and scales ranging from C. elegans to mouse, rat, cat, macaque, and human, suggesting a role for genetic influences. The recent availability of brain-wide gene expression atlases provides new opportunities for mapping the transcriptional correlates of large-scale network-level phenotypes. Here we review studies that use these atlases to investigate gene expression patterns associated with hub connectivity in neural networks and present evidence that some of these patterns are conserved across species and scales.

Dynamic spatiotemporal patterns of alternative splicing of an F-actin scaffold protein, afadin, during murine development

An F-actin scaffold protein, afadin, comprises two splice variants called l-afadin (a long isoform) and s-afadin (a short isoform). It is known that in adult tissues, l-afadin is ubiquitously expressed while s-afadin is restrictedly expressed in brain. In cultured cortical neurons, l-afadin potentiates axonal branching whereas s-afadin blocks axonal branching by functioning as a naturally occurring dominant-negative isoform that forms a heterodimer with l-afadin. However, the temporal and spatial expression pattern of s-afadin during development or across multiple tissues and organs has not been fully understood. In this study, using Western blotting and RT-qPCR techniques and the fluorescent splicing reporters, we examined the detailed expression patterns of l- and s-afadin isoforms across various tissues and cell types, including rat organs at developmental stages, primary cultured neuronal and non-neuronal cells prepared from the developing rat brain, and in neurons in vitro generated from P19 embryonal carcinoma (EC) cells. Both mRNA and protein of s-afadin were abundantly expressed in various regions of rat neuronal tissues, and their expression dynamically changed during development in vivo. The expression of s-afadin was also detected in primary rat cortical neurons, but not in astrocytes or fibroblasts, and the neuronal expression increased during neuronal maturation in vitro. The dynamic alternative splicing event of afadin during development was successfully visualized with the newly developed fluorescent splicing reporter plasmids at a single cell level. Moreover, s-afadin was undetectable in undifferentiated EC cells, and the expression was induced upon differentiation of the cells into neurons. These data suggest that spatiotemporal and dynamic alternative splicing produces differential expression patterns of afadin isoforms and that alternative splicing of afadin is controlled by the neuron-specific regulator(s) whose activity is triggered and dynamically altered during neuronal differentiation and maturation.

A gene expression signature in developing Purkinje cells predicts autism and intellectual disability co-morbidity status

Autism spectrum disorder (ASD) is a complex neurodevelopmental disease whose underpinning molecular mechanisms and neural substrates are subject to intense scrutiny. Interestingly, the cerebellum has emerged as one of the key brain regions affected in ASD. However, the genetic and molecular mechanisms that link the cerebellum to ASD, particularly during development, remain poorly understood. To gain insight into the genetic and molecular mechanisms that might link the cerebellum to ASD, we analysed the transcriptome dynamics of a developing cell population highly enriched for Purkinje cells of the mouse cerebellum across multiple timepoints. We identified a single cluster of genes whose expression is positively correlated with development and which is enriched for genes associated with ASD. This ASD-associated gene cluster was specific to developing Purkinje cells and not detected in the mouse neocortex during the same developmental period, in which we identified a distinct temporally regulated ASD gene module. Furthermore, the composition of ASD risk genes within the two distinct clusters was significantly different in their association with intellectual disability (ID), consistent with the existence of genetically and spatiotemporally distinct endophenotypes of ASD. Together, our findings define a specific cluster of ASD genes that is enriched in developing PCs and predicts co-morbidity status.

A practical guide to linking brain-wide gene expression and neuroimaging data

The recent availability of comprehensive, brain-wide gene expression atlases such as the Allen Human Brain Atlas (AHBA) has opened new opportunities for understanding how spatial variations on molecular scale relate to the macroscopic neuroimaging phenotypes. A rapidly growing body of literature is demonstrating relationships between gene expression and diverse properties of brain structure and function, but approaches for combining expression atlas data with neuroimaging are highly inconsistent, with substantial variations in how the expression data are processed. The degree to which these methodological variations affect findings is unclear. Here, we outline a seven-step analysis pipeline for relating brain-wide transcriptomic and neuroimaging data and compare how different processing choices influence the resulting data. We suggest that studies using the AHBA should work towards a unified data processing pipeline to ensure consistent and reproducible results in this burgeoning field.

The novel lncRNA lnc-NR2F1 is pro-neurogenic and mutated in human neurodevelopmental disorders

Long noncoding RNAs (lncRNAs) have been shown to act as important cell biological regulators including cell fate decisions but are often ignored in human genetics. Combining differential lncRNA expression during neuronal lineage induction with copy number variation morbidity maps of a cohort of children with autism spectrum disorder/intellectual disability versus healthy controls revealed focal genomic mutations affecting several lncRNA candidate loci. Here we find that a t(5:12) chromosomal translocation in a family manifesting neurodevelopmental symptoms disrupts specifically lnc-NR2F1. We further show that lnc-NR2F1 is an evolutionarily conserved lncRNA functionally enhances induced neuronal cell maturation and directly occupies and regulates transcription of neuronal genes including autism-associated genes. Thus, integrating human genetics and functional testing in neuronal lineage induction is a promising approach for discovering candidate lncRNAs involved in neurodevelopmental diseases.

The spatiotemporal expression pattern of microRNAs in the developing mouse nervous system

MicroRNAs (miRNAs) control various biological processes by inducing translational repression and transcript degradation of the target genes. In mammalian development, knowledge of the timing and expression pattern of each miRNA is important to determine and predict its function in vivo So far, no systematic analyses of the spatiotemporal expression pattern of miRNAs during mammalian neurodevelopment have been performed. Here, we isolated total RNAs from the embryonic dorsal forebrain of mice at different developmental stages and subjected these RNAs to microarray analyses. We selected 279 miRNAs that exhibited high signal intensities or ascending or descending expression dynamics. To ascertain the expression patterns of these miRNAs, we used locked nucleic acid (LNA)-modified miRNA probes in in situ hybridization experiments. Multiple miRNAs exhibited spatially restricted/enriched expression in anatomically distinct regions or in specific neuron subtypes in the embryonic brain and spinal cord, such as in the ventricular area, the striatum (and other basal ganglia), hypothalamus, choroid plexus, and the peripheral nervous system. These findings provide new insights into the expression and function of miRNAs during the development of the nervous system and could be used as a resource to facilitate studies in neurodevelopment.

Variation among intact tissue samples reveals the core transcriptional features of human CNS cell classes

It is widely assumed that cells must be physically isolated to study their molecular profiles. However, intact tissue samples naturally exhibit variation in cellular composition, which drives covariation of cell-class-specific molecular features. By analyzing transcriptional covariation in 7,221 intact CNS samples from 840 neurotypical individuals, representing billions of cells, we reveal the core transcriptional identities of major CNS cell classes in humans. By modeling intact CNS transcriptomes as a function of variation in cellular composition, we identify cell-class-specific transcriptional differences in Alzheimer's disease, among brain regions, and between species. Among these, we show that PMP2 is expressed by human but not mouse astrocytes and significantly increases mouse astrocyte size upon ectopic expression in vivo, causing them to more closely resemble their human counterparts. Our work is available as an online resource ( http://oldhamlab.ctec.ucsf.edu/ ) and provides a generalizable strategy for determining the core molecular features of cellular identity in intact biological systems.

Precise temporal regulation of alternative splicing during neural development

Alternative splicing (AS) is one crucial step of gene expression that must be tightly regulated during neurodevelopment. However, the precise timing of developmental splicing switches and the underlying regulatory mechanisms are poorly understood. Here we systematically analyze the temporal regulation of AS in a large number of transcriptome profiles of developing mouse cortices, in vivo purified neuronal subtypes, and neurons differentiated in vitro. Our analysis reveals early-switch and late-switch exons in genes with distinct functions, and these switches accurately define neuronal maturation stages. Integrative modeling suggests that these switches are under direct and combinatorial regulation by distinct sets of neuronal RNA-binding proteins including Nova, Rbfox, Mbnl, and Ptbp. Surprisingly, various neuronal subtypes in the sensory systems lack Nova and/or Rbfox expression. These neurons retain the "immature" splicing program in early-switch exons, affecting numerous synaptic genes. These results provide new insights into the organization and regulation of the neurodevelopmental transcriptome.

Brain Transcriptome Databases: A User's Guide

Transcriptional programs instruct the generation and maintenance of diverse subtypes of neural cells, establishment of distinct brain regions, formation and function of neural circuits, and ultimately behavior. Spatiotemporal and cell type-specific analyses of the transcriptome, the sum total of all RNA transcripts in a cell or an organ, can provide insights into the role of genes in brain development and function, and their potential contribution to disorders of the brain. In the previous decade, advances in sequencing technology and funding from the National Institutes of Health and private foundations for large-scale genomics projects have led to a growing collection of brain transcriptome databases. These valuable resources provide rich and high-quality datasets with spatiotemporal, cell type-specific, and single-cell precision. Most importantly, many of these databases are publicly available via user-friendly web interface, making the information accessible to individual scientists without the need for advanced computational expertise. Here, we highlight key publicly available brain transcriptome databases, summarize the tissue sources and methods used to generate the data, and discuss their utility for neuroscience research.

Single-cell transcriptomics of the developing lateral geniculate nucleus reveals insights into circuit assembly and refinement

Neurons and nonneuronal cells in the developing brain dynamically regulate gene expression as neural connectivity is established. However, the specific gene programs activated in distinct cell populations during the assembly and refinement of many intact neuronal circuits have not been thoroughly characterized. In this study, we take advantage of recent advances in transcriptomic profiling techniques to characterize gene expression in the postnatal developing lateral geniculate nucleus (LGN) at single-cell resolution. Our data reveal that genes involved in brain development are dynamically regulated in all major cell types of the LGN, suggesting that the establishment of neural connectivity depends upon functional collaboration between multiple neuronal and nonneuronal cell types in this brain region.

Coordinated changes in gene expression underlie the early patterning and cell-type specification of the central nervous system. However, much less is known about how such changes contribute to later stages of circuit assembly and refinement. In this study, we employ single-cell RNA sequencing to develop a detailed, whole-transcriptome resource of gene expression across four time points in the developing dorsal lateral geniculate nucleus (LGN), a visual structure in the brain that undergoes a well-characterized program of postnatal circuit development. This approach identifies markers defining the major LGN cell types, including excitatory relay neurons, oligodendrocytes, astrocytes, microglia, and endothelial cells. Most cell types exhibit significant transcriptional changes across development, dynamically expressing genes involved in distinct processes including retinotopic mapping, synaptogenesis, myelination, and synaptic refinement. Our data suggest that genes associated with synapse and circuit development are expressed in a larger proportion of nonneuronal cell types than previously appreciated. Furthermore, we used this single-cell expression atlas to identify the Prkcd-Cre mouse line as a tool for selective manipulation of relay neurons during a late stage of sensory-driven synaptic refinement. This transcriptomic resource provides a cellular map of gene expression across several cell types of the LGN, and offers insight into the molecular mechanisms of circuit development in the postnatal brain.

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.

The BRAIN Initiative Cell Census Consortium: Lessons Learned toward Generating a Comprehensive Brain Cell Atlas

A comprehensive characterization of neuronal cell types, their distributions, and patterns of connectivity is critical for understanding the properties of neural circuits and how they generate behaviors. Here we review the experiences of the BRAIN Initiative Cell Census Consortium, ten pilot projects funded by the U.S. BRAIN Initiative, in developing, validating, and scaling up emerging genomic and anatomical mapping technologies for creating a complete inventory of neuronal cell types and their connections in multiple species and during development. These projects lay the foundation for a larger and longer-term effort to generate whole-brain cell atlases in species including mice and humans.

Maize network analysis revealed gene modules involved in development, nutrients utilization, metabolism, and stress response

BACKGROUND

The advent of big data in biology offers opportunities while poses challenges to derive biological insights. For maize, a large amount of publicly available transcriptome datasets have been generated but a comprehensive analysis is lacking.

RESULTS

We constructed a maize gene co-expression network based on the graphical Gaussian model, using massive RNA-seq data. The network, containing 20,269 genes, assembles into 964 gene modules that function in a variety of plant processes, such as cell organization, the development of inflorescences, ligules and kernels, the uptake and utilization of nutrients (e.g. nitrogen and phosphate), the metabolism of benzoxazionids, oxylipins, flavonoids, and wax, and the response to stresses. Among them, the inflorescences development module is enriched with domestication genes (like ra1, ba1, gt1, tb1, tga1) that control plant architecture and kernel structure, while multiple other modules relate to diverse agronomic traits. Contained within these modules are transcription factors acting as known or potential expression regulators for the genes within the same modules, suggesting them as candidate regulators for related biological processes. A comparison with an established Arabidopsis network revealed conserved gene association patterns for specific modules involved in cell organization, nutrients uptake & utilization, and metabolism. The analysis also identified significant divergences between the two species for modules that orchestrate developmental pathways.

CONCLUSIONS

This network sheds light on how gene modules are organized between different species in the context of evolutionary divergence and highlights modules whose structure and gene content can provide important resources for maize gene functional studies with application potential.

Annotation and cluster analysis of spatiotemporal- and sex-related lncRNA expression in rhesus macaque brain

Long noncoding RNAs (lncRNAs) mediate important epigenetic regulation in a wide range of biological processes and diseases. We applied comprehensive analyses of RNA-seq and CAGE-seq (cap analysis of gene expression and sequencing) to characterize the dynamic changes in lncRNA expression in rhesus macaque (Macaca mulatta) brain in four representative age groups. We identified 18 anatomically diverse lncRNA modules and 14 mRNA modules representing spatial, age, and sex specificities. Spatiotemporal- and sex-biased changes in lncRNA expression were generally higher than those observed in mRNA expression. A negative correlation between lncRNA and mRNA expression in cerebral cortex was observed and functionally validated. Our findings offer a fresh insight into spatial-, age-, and sex-biased changes in lncRNA expression in macaque brain and suggest that the changes represent a previously unappreciated regulatory system that potentially contributes to brain development and aging.

Brain transcriptome atlases: a computational perspective

The immense complexity of the mammalian brain is largely reflected in the underlying molecular signatures of its billions of cells. Brain transcriptome atlases provide valuable insights into gene expression patterns across different brain areas throughout the course of development. Such atlases allow researchers to probe the molecular mechanisms which define neuronal identities, neuroanatomy, and patterns of connectivity. Despite the immense effort put into generating such atlases, to answer fundamental questions in neuroscience, an even greater effort is needed to develop methods to probe the resulting high-dimensional multivariate data. We provide a comprehensive overview of the various computational methods used to analyze brain transcriptome atlases.

Comprehensive transcriptome analysis of neocortical layers in humans, chimpanzees and macaques

While human cognitive abilities are clearly unique, underlying changes in brain organization and function remain unresolved. Here we characterized the transcriptome of the cortical layers and adjacent white matter in the prefrontal cortexes of humans, chimpanzees and rhesus macaques using unsupervised sectioning followed by RNA sequencing. More than 20% of detected genes were expressed predominantly in one layer, yielding 2,320 human layer markers. While the bulk of the layer markers were conserved among species, 376 switched their expression to another layer in humans. By contrast, only 133 of such changes were detected in the chimpanzee brain, suggesting acceleration of cortical reorganization on the human evolutionary lineage. Immunohistochemistry experiments further showed that human-specific expression changes were not limited to neurons but affected a broad spectrum of cortical cell types. Thus, despite apparent histological conservation, human neocortical organization has undergone substantial changes affecting more than 5% of its transcriptome.

MicroRNAs contribute to postnatal development of laminar differences and neuronal subtypes in the rat medial entorhinal cortex

The medial entorhinal cortex (MEC) is important in spatial navigation and memory formation and its layers have distinct neuronal subtypes, connectivity, spatial properties, and disease susceptibility. As little is known about the molecular basis for the development of these laminar differences, we analyzed microRNA (miRNA) and messenger RNA (mRNA) expression differences between rat MEC layer II and layers III–VI during postnatal development. We identified layer and age-specific regulation of gene expression by miRNAs, which included processes related to neuron specialization and locomotor behavior. Further analyses by retrograde labeling and expression profiling of layer II stellate neurons and in situ hybridization revealed that the miRNA most up-regulated in layer II, miR-143, was enriched in stellate neurons, whereas the miRNA most up-regulated in deep layers, miR-219-5p, was expressed in ependymal cells, oligodendrocytes and glia. Bioinformatics analyses of predicted mRNA targets with negatively correlated expression patterns to miR-143 found that miR-143 likely regulates the Lmo4 gene, which is known to influence hippocampal-based spatial learning.

The frontier of RNA metamorphosis and ribosome signature in neocortical development

More than a passive effector of gene expression, mRNA translation (protein synthesis) by the ribosome is a rapidly tunable and dynamic molecular mechanism. Neurodevelopmental disorders are associated with abnormalities in mRNA translation, protein synthesis, and neocortical development; yet, we know little about the molecular mechanisms underlying these abnormalities. Furthermore, our understanding of regulation of the ribosome and mRNA translation during normal brain development is only in its early stages. mRNA translation is emerging as a key driver of the rapid and timed regulation of spatiotemporal gene expression in the developing nervous system, including the neocortex. In this review, we focus on the regulatory role of the ribosome in neocortical development, and construct a current understanding of how ribosomal complex specificity may contribute to the development of the neocortex. We also present a microarray analysis of ribosomal protein-coding mRNAs across the neurogenic phase of neocortical development, in addition to the dynamic enrichment of these mRNAs in actively translating neocortical polysomal ribosomes. Understanding the multivariate control of mRNA translation by ribosomal complex specificity will be critical to reveal the intricate mechanisms of normal brain development and pathologies of neurodevelopmental disorders.

A Novel Analytical Strategy to Identify Fusion Transcripts between Repetitive Elements and Protein Coding-Exons Using RNA-Seq

Repetitive elements (REs) comprise 40-60% of the mammalian genome and have been shown to epigenetically influence the expression of genes through the formation of fusion transcript (FTs). We previously showed that an intracisternal A particle forms an FT with the agouti gene in mice, causing obesity/type 2 diabetes. To determine the frequency of FTs genome-wide, we developed a TopHat-Fusion-based analytical pipeline to identify FTs with high specificity. We applied it to an RNA-seq dataset from the nucleus accumbens (NAc) of mice repeatedly exposed to cocaine. Cocaine was previously shown to increase the expression of certain REs in this brain region. Using this pipeline that can be applied to single- or paired-end reads, we identified 438 genes expressing 813 different FTs in the NAc. Although all types of studied repeats were present in FTs, simple sequence repeats were underrepresented. Most importantly, reverse-transcription and quantitative PCR validated the expression of selected FTs in an independent cohort of animals, which also revealed that some FTs are the prominent isoforms expressed in the NAc by some genes. In other RNA-seq datasets, developmental expression as well as tissue specificity of some FTs differed from their corresponding non-fusion counterparts. Finally, in silico analysis predicted changes in the structure of proteins encoded by some FTs, potentially resulting in gain or loss of function. Collectively, these results indicate the robustness of our pipeline in detecting these new isoforms of genes, which we believe provides a valuable tool to aid in better understanding the broad role of REs in mammalian cellular biology.

RNA-seq analysis of developing olfactory bulb projection neurons

Transmission of olfactory information to higher brain regions is mediated by olfactory bulb (OB) projection neurons, the mitral and tufted cells. Although mitral/tufted cells are often characterized as the OB counterpart of cortical projection neurons (also known as pyramidal neurons), they possess several unique morphological characteristics and project specifically to the olfactory cortices. Moreover, the molecular networks contributing to the generation of mitral/tufted cells during development are largely unknown. To understand the developmental patterns of gene expression in mitral/tufted cells in the OB, we performed transcriptome analyses targeting purified OB projection neurons at different developmental time points with next-generation RNA sequencing (RNA-seq). Through these analyses, we found 1202 protein-coding genes that are temporally differentially-regulated in developing projection neurons. Among them, 388 genes temporally changed their expression level only in projection neurons. The data provide useful resource to study the molecular mechanisms regulating development of mitral/tufted cells. We further compared the gene expression profiles of developing mitral/tufted cells with those of three cortical projection neuron subtypes, subcerebral projection neurons, corticothalamic projection neurons, and callosal projection neurons, and found that the molecular signature of developing olfactory projection neuron bears resemblance to that of subcerebral neurons. We also identified 3422 events that change the ratio of splicing isoforms in mitral/tufted cells during maturation. Interestingly, several genes expressed a novel isoform not previously reported. These results provide us with a broad perspective of the molecular networks underlying the development of OB projection neurons.

A Markov random field-based approach for joint estimation of differentially expressed genes in mouse transcriptome data

The statistical methodology developed in this study was motivated by our interest in studying neurodevelopment using the mouse brain RNA-Seq data set, where gene expression levels were measured in multiple layers in the somatosensory cortex across time in both female and male samples. We aim to identify differentially expressed genes between adjacent time points, which may provide insights on the dynamics of brain development. Because of the extremely small sample size (one male and female at each time point), simple marginal analysis may be underpowered. We propose a Markov random field (MRF)-based approach to capitalizing on the between layers similarity, temporal dependency and the similarity between sex. The model parameters are estimated by an efficient EM algorithm with mean field-like approximation. Simulation results and real data analysis suggest that the proposed model improves the power to detect differentially expressed genes than simple marginal analysis. Our method also reveals biologically interesting results in the mouse brain RNA-Seq data set.

Sequential transcriptional waves direct the differentiation of newborn neurons in the mouse neocortex

During corticogenesis, excitatory neurons are born from progenitors located in the ventricular zone (VZ), from where they migrate to assemble into circuits. How neuronal identity is dynamically specified upon progenitor division is unknown. Here, we study this process using a high-temporal-resolution technology allowing fluorescent tagging of isochronic cohorts of newborn VZ cells. By combining this in vivo approach with single-cell transcriptomics in mice, we identify and functionally characterize neuron-specific primordial transcriptional programs as they dynamically unfold. Our results reveal early transcriptional waves that instruct the sequence and pace of neuronal differentiation events, guiding newborn neurons toward their final fate, and contribute to a road map for the reverse engineering of specific classes of cortical neurons from undifferentiated cells.

Thalamic WNT3 Secretion Spatiotemporally Regulates the Neocortical Ribosome Signature and mRNA Translation to Specify Neocortical Cell Subtypes

Neocortical development requires tightly controlled spatiotemporal gene expression. However, the mechanisms regulating ribosomal complexes and the timed specificity of neocortical mRNA translation are poorly understood. We show that active mRNA translation complexes (polysomes) contain ribosomal protein subsets that undergo dynamic spatiotemporal rearrangements during mouse neocortical development. Ribosomal protein specificity within polysome complexes is regulated by the arrival of in-growing thalamic axons, which secrete the morphogen Wingless-related MMTV (mouse mammary tumor virus) integration site 3 (WNT3). Thalamic WNT3 release during midneurogenesis promotes a change in the levels of Ribosomal protein L7 in polysomes, thereby regulating neocortical translation machinery specificity. Furthermore, we present an RNA sequencing dataset analyzing mRNAs that dynamically associate with polysome complexes as neocortical development progresses, and thus may be regulated spatiotemporally at the level of translation. Thalamic WNT3 regulates neocortical translation of two such mRNAs, Foxp2 and Apc, to promote FOXP2 expression while inhibiting APC expression, thereby driving neocortical neuronal differentiation and suppressing oligodendrocyte maturation, respectively. This mechanism may enable targeted and rapid spatiotemporal control of ribosome composition and selective mRNA translation in complex developing systems like the neocortex.

SIGNIFICANCE STATEMENT

The neocortex is a highly complex circuit generating the most evolutionarily advanced complex cognitive and sensorimotor functions. An intricate progression of molecular and cellular steps during neocortical development determines its structure and function. Our goal is to study the steps regulating spatiotemporal specificity of mRNA translation that govern neocortical development. In this work, we show that the timed secretion of Wingless-related MMTV (mouse mammary tumor virus) integration site 3 (WNT3) by ingrowing axons from the thalamus regulates the combinatorial composition of ribosomal proteins in developing neocortex, which we term the "neocortical ribosome signature." Thalamic WNT3 further regulates the specificity of mRNA translation and development of neurons and oligodendrocytes in the neocortex. This study advances our overall understanding of WNT signaling and the spatiotemporal regulation of mRNA translation in highly complex developing systems.

Global deletion of BCATm increases expression of skeletal muscle genes associated with protein turnover

Consumption of a protein-containing meal by a fasted animal promotes protein accretion in skeletal muscle, in part through leucine stimulation of protein synthesis and indirectly through repression of protein degradation mediated by its metabolite, α-ketoisocaproate. Mice lacking the mitochondrial branched-chain aminotransferase (BCATm/Bcat2), which interconverts leucine and α-ketoisocaproate, exhibit elevated protein turnover. Here, the transcriptomes of gastrocnemius muscle from BCATm knockout (KO) and wild-type mice were compared by next-generation RNA sequencing (RNA-Seq) to identify potential adaptations associated with their persistently altered nutrient signaling. Statistically significant changes in the abundance of 1,486/∼39,010 genes were identified. Bioinformatics analysis of the RNA-Seq data indicated that pathways involved in protein synthesis [eukaryotic initiation factor (eIF)-2, mammalian target of rapamycin, eIF4, and p70S6K pathways including 40S and 60S ribosomal proteins], protein breakdown (e.g., ubiquitin mediated), and muscle degeneration (apoptosis, atrophy, myopathy, and cell death) were upregulated. Also in agreement with our previous observations, the abundance of mRNAs associated with reduced body size, glycemia, plasma insulin, and lipid signaling pathways was altered in BCATm KO mice. Consistently, genes encoding anaerobic and/or oxidative metabolism of carbohydrate, fatty acids, and branched chain amino acids were modestly but systematically reduced. Although there was no indication that muscle fiber type was different between KO and wild-type mice, a difference in the abundance of mRNAs associated with a muscular dystrophy phenotype was observed, consistent with the published exercise intolerance of these mice. The results suggest transcriptional adaptations occur in BCATm KO mice that along with altered nutrient signaling may contribute to their previously reported protein turnover, metabolic and exercise phenotypes.

Transcriptomics of critical period of visual cortical plasticity in mice

In contrast to the prenatal development of the cerebral cortex, when cell production, migration, and layer formation dominate, development after birth involves more subtle processes, such as activity-dependent plasticity that includes refinement of synaptic connectivity by its stabilization and elimination. In the present study, we use RNA-seq with high spatial resolution to examine differential gene expression across layers 2/3, 4, 5, and 6 of the mouse visual cortex before the onset of the critical period of plasticity [postnatal day 5 (P5)], at its peak (P26), and at the mature stage (P180) and compare it with the prefrontal association area. We find that, although genes involved in early developmental events such as cell division, neuronal migration, and axon guidance are still prominent at P5, their expression largely terminates by P26, when synaptic plasticity and associated signaling pathways become enriched. Unexpectedly, the gene expression profile was similar in both areas at this age, suggesting that activity-dependent plasticity between visual and association cortices are subject to the same genetic constraints. Although gene expression changes follow similar paths until P26, we have identified 30 regionally enriched genes that are prominent during the critical period. At P180, we identified several hundred differentially expressed gene isoforms despite subsiding levels of gene expression differences. This result indicates that, once genetic developmental programs cease, the remaining morphogenetic processes may depend on posttranslational events.