The role of Tal2 and Tal1 in the differentiation of midbrain GABAergic neuron precursors

Specification of human regional brain lineages using orthogonal gradients of WNT and SHH in organoids

The repertory of neurons generated by progenitor cells depends on their location along antero-posterior and dorso-ventral axes of the neural tube. To understand if recreating those axes was sufficient to specify human brain neuronal diversity, we designed a mesofluidic device termed Duo-MAPS to expose induced pluripotent stem cells (iPSC) to concomitant orthogonal gradients of a posteriorizing and a ventralizing morphogen, activating WNT and SHH signaling, respectively. Comparison of single cell transcriptomes with fetal human brain revealed that Duo-MAPS-patterned organoids generated the major neuronal lineages of the forebrain, midbrain, and hindbrain. Morphogens crosstalk translated into early patterns of gene expression programs predicting the generation of specific brain lineages. Human iPSC lines from six different genetic backgrounds showed substantial differences in response to morphogens, suggesting that interindividual genomic and epigenomic variations could impact brain lineages formation. Morphogen gradients promise to be a key approach to model the brain in its entirety.

Chromatin accessibility during human first-trimester neurodevelopment

The human brain develops through a tightly organized cascade of patterning events, induced by transcription factor expression and changes in chromatin accessibility. Although gene expression across the developing brain has been described at single-cell resolution1, similar atlases of chromatin accessibility have been primarily focused on the forebrain2-4. Here we describe chromatin accessibility and paired gene expression across the entire developing human brain during the first trimester (6-13 weeks after conception). We defined 135 clusters and used multiomic measurements to link candidate cis-regulatory elements to gene expression. The number of accessible regions increased both with age and along neuronal differentiation. Using a convolutional neural network, we identified putative functional transcription factor-binding sites in enhancers characterizing neuronal subtypes. We applied this model to cis-regulatory elements linked to ESRRB to elucidate its activation mechanism in the Purkinje cell lineage. Finally, by linking disease-associated single nucleotide polymorphisms to cis-regulatory elements, we validated putative pathogenic mechanisms in several diseases and identified midbrain-derived GABAergic neurons as being the most vulnerable to major depressive disorder-related mutations. Our findings provide a more detailed view of key gene regulatory mechanisms underlying the emergence of brain cell types during the first trimester and a comprehensive reference for future studies related to human neurodevelopment.

Range of chromatin accessibility configurations are permissive of GABAergic fate acquisition in developing mouse brain

In recent single-cell -omics studies, both the differential activity of transcription factors regulating cell fate determination and differential genome activation have been tested for utility as descriptors of cell types. Naturally, genome accessibility and gene expression are interlinked. To understand the variability in genomic feature activation in the GABAergic neurons of different spatial origins, we have mapped accessible chromatin regions and mRNA expression in single cells derived from the developing mouse central nervous system (CNS). We first defined a reference set of open chromatin regions for scATAC-seq read quantitation across samples, allowing comparison of chromatin accessibility between brain regions and cell types directly. Second, we integrated the scATAC-seq and scRNA-seq data to form a unified resource of transcriptome and chromatin accessibility landscape for the cell types in di- and telencephalon, midbrain and anterior hindbrain of E14.5 mouse embryo. Importantly, we implemented resolution optimization at the clustering, and automatized the cell typing step. We show high level of concordance between the cell clustering based on the chromatin accessibility and the transcriptome in analyzed neuronal lineages, indicating that both genome and transcriptome features can be used for cell type definition. Hierarchical clustering by the similarity in accessible chromatin reveals that the genomic feature activation correlates with neurotransmitter phenotype, selector gene expression, cell differentiation stage and neuromere origins.

Developmental dynamics of the single nucleus regulatory landscape of pig hippocampus

The hippocampus is a brain region associated with memory, learning and spatial navigation, its aging-related dysfunction is a common sign of Alzheimer's disease. Pig is a good model for human neurodegenerative disease, but our understanding of the regulatory program of the pig hippocampus and its cross-species conservation in humans remains limited. Here, we profiled chromatin accessibility in 33,409 high-quality nuclei and gene expression in 8,122 high-quality nuclei of the pig hippocampus at four postnatal stages. We identified 510,908 accessible chromatin regions (ACRs) in 12 major cell types, among which progenitor cells such as neuroblasts and oligodendrocyte progenitor cells showed a dynamic decrease from early to later developmental stages. We revealed significant enrichment of transposable elements in cell type-specific ACRs, particularly in neuroblasts. We identified oligodendrocytes as the most prominent cell type with the greatest number of genes that showed significant changes during the development. We identified ACRs and key transcription factors underlying the trajectory of neurogenesis (such as POU3F3 and EGR1) and oligodendrocyte differentiation (RXRA and FOXO6). We examined 27 Alzheimer's disease-related genes in our data and found that 15 showed cell type-specific activity (TREM2, RIN3 and CLU), and 15 genes displayed age-associated dynamic activity (BIN1, RABEP1 and APOE). We intersected our data with human genome-wide association study results to detect neurological disease-associated cell types. The present study provides a single nucleus-accessible chromatin landscape of the pig hippocampus at different developmental stages and is helpful for the exploration of pigs as a biomedical model in human neurodegenerative diseases.

A lamprey neural cell type atlas illuminates the origins of the vertebrate brain

The vertebrate brain emerged more than ~500 million years ago in common evolutionary ancestors. To systematically trace its cellular and molecular origins, we established a spatially resolved cell type atlas of the entire brain of the sea lamprey-a jawless species whose phylogenetic position affords the reconstruction of ancestral vertebrate traits-based on extensive single-cell RNA-seq and in situ sequencing data. Comparisons of this atlas to neural data from the mouse and other jawed vertebrates unveiled various shared features that enabled the reconstruction of cell types, tissue structures and gene expression programs of the ancestral vertebrate brain. However, our analyses also revealed key tissues and cell types that arose later in evolution. For example, the ancestral brain was probably devoid of cerebellar cell types and oligodendrocytes (myelinating cells); our data suggest that the latter emerged from astrocyte-like evolutionary precursors in the jawed vertebrate lineage. Altogether, our work illuminates the cellular and molecular architecture of the ancestral vertebrate brain and provides a foundation for exploring its diversification during evolution.

The embryonic patterning gene Dbx1 governs the survival of the auditory midbrain via Tcf7l2-Ap2δ transcriptional cascade

At the top of the midbrain is the inferior colliculus (IC), which functions as the major hub for processing auditory information. Despite the functional significance of neurons in the IC, our understanding of their formation is limited. In this study, we identify the embryonic patterning gene Dbx1 as a key molecular player that governs genetic programs for IC survival. We find that Dbx1 plays a critical role in preventing apoptotic cell death in postnatal IC by transcriptionally repressing c-Jun and pro-apoptotic BH3 only factors. Furthermore, by employing combined approaches, we uncover that Tcf7l2 functions downstream of Dbx1. Loss of Tcf7l2 function causes IC phenotypes with striking similarity to those of Dbx1 mutant mice, which include defective embryonic maturation and postnatal deletion of the IC. Finally, we demonstrate that the Dbx1-Tcf7l2 cascade functions upstream of Ap-2δ, which is essential for IC development and survival. Together, these results unravel a novel molecular mechanism for IC maintenance, which is indispensable for normal brain development.

Postsynaptic NMDA Receptor Expression Is Required for Visual Corticocollicular Projection Refinement in the Mouse Superior Colliculus

Efficient sensory processing of spatial information is facilitated through the organization of neuronal connections into topographic maps of space. In integrative sensory centers, converging topographic maps must be aligned to merge spatially congruent information. The superior colliculus (SC) receives topographically ordered visual inputs from retinal ganglion cells (RGCs) in the eye and layer 5 neurons in the primary visual cortex (L5-V1). Previous studies suggest that RGCs instruct the alignment of later-arriving L5-V1 inputs in an activity-dependent manner. However, the molecular mechanisms underlying this remain unclear. Here, we explored the role of NMDA receptors in visual map alignment in the SC using a conditional genetic knockout approach. We leveraged a novel knock-in mouse line that expresses tamoxifen-inducible Cre recombinase under the control of the Tal1 gene (Tal1CreERT2 ), which we show allows for specific recombination in the superficial layers of the SC. We used Tal1CreERT2 mice of either sex to conditionally delete the obligate GluN1 subunit of the NMDA receptor (SC-cKO) during the period of visual map alignment. We observed a significant disruption of L5-V1 axon terminal organization in the SC of SC-cKO mice. Importantly, retinocollicular topography was unaffected in this context, suggesting that alignment is also disrupted. Time-course experiments suggest that NMDA receptors may play a critical role in the refinement of L5-V1 inputs in the SC. Together, these data implicate NMDA receptors as critical mediators of activity-dependent visual map alignment in the SC.SIGNIFICANCE STATEMENT Alignment of topographic inputs is critical for integration of spatially congruent sensory information; however, little is known about the mechanisms underlying this complex process. Here, we took a conditional genetic approach to explore the role of NMDA receptors in the alignment of retinal and cortical visual inputs in the superior colliculus. We characterize a novel mouse line providing spatial and temporal control of recombination in the superior colliculus and reveal a critical role for NMDA expression in visual map alignment. These data support a role for neuronal activity in visual map alignment and provide mechanistic insight into this complex developmental process.

Building thalamic neuronal networks during mouse development

The thalamic nuclear complex contains excitatory projection neurons and inhibitory local neurons, the two cell types driving the main circuits in sensory nuclei. While excitatory neurons are born from progenitors that reside in the proliferative zone of the developing thalamus, inhibitory local neurons are born outside the thalamus and they migrate there during development. In addition to these cell types, which occupy most of the thalamus, there are two small thalamic regions where inhibitory neurons target extra-thalamic regions rather than neighboring neurons, the intergeniculate leaflet and the parahabenular nucleus. Like excitatory thalamic neurons, these inhibitory neurons are derived from progenitors residing in the developing thalamus. The assembly of these circuits follows fine-tuned genetic programs and it is coordinated by extrinsic factors that help the cells find their location, associate with thalamic partners, and establish connections with their corresponding extra-thalamic inputs and outputs. In this review, we bring together what is currently known about the development of the excitatory and inhibitory components of the thalamocortical sensory system, in particular focusing on the visual pathway and thalamic interneurons in mice.

Tal2 is required for generation of GABAergic neurons in the zebrafish midbrain

BACKGROUND

In the zebrafish midbrain, GABAergic neurons develop from precursors located in the nucleus of the medial longitudinal fasciculus (nMLF). However, the precise mechanisms that underline generation of the nMLF GABAergic neuron are poorly understood.

RESULTS

GABAergic neurons in the nMLF co-express transcription factors tal2, gata2a, gata3, and nkx1.2lb. The Nodal-related gene and shh signaling are required for differentiation of nMLF GABAergic neuron precursors. Tal2 is important for nMLF GABAergic neurogenesis. Disruption of Tal2, embryos completely lack the GABA-synthesizing enzyme glutamic acid decarboxylase 67 gene (gad67) expressing cells in the nMLF, and the whole nkx1.2lb expressing cells in the midbrain. Although almost all tal2-expressing cells in the diencephalon and/or nMLF are gata2a- and gata3-positive, simultaneous knockdown of gata2a and gata3 does not affect either tal2 or gad67 expression.

CONCLUSIONS

In the zebrafish midbrain, expression of tal2, gata2a, and/or gata3 is independent of each other. The function of gata2a and gata3 is dispensable for generation of GABAergic neuron in the nMLF. This suggests that the functional connections of the regulatory genes leading to generation of nMLF GABAergic neurons have diverged between mouse and zebrafish.

Neurons gating behavior—developmental, molecular and functional features of neurons in the Substantia Nigra pars reticulata

The Substantia Nigra pars reticulata (SNpr) is the major information output site of the basal ganglia network and instrumental for the activation and adjustment of movement, regulation of the behavioral state and response to reward. Due to both overlapping and unique input and output connections, the SNpr might also have signal integration capacity and contribute to action selection. How the SNpr regulates these multiple functions remains incompletely understood. The SNpr is located in the ventral midbrain and is composed primarily of inhibitory GABAergic projection neurons that are heterogeneous in their properties. In addition, the SNpr contains smaller populations of other neurons, including glutamatergic neurons. Here, we discuss regionalization of the SNpr, in particular the division of the SNpr neurons to anterior (aSNpr) and posterior (pSNpr) subtypes, which display differences in many of their features. We hypothesize that unique developmental and molecular characteristics of the SNpr neuron subtypes correlate with both region-specific connections and notable functional specializations of the SNpr. Variation in both the genetic control of the SNpr neuron development as well as signals regulating cell migration and axon guidance may contribute to the functional diversity of the SNpr neurons. Therefore, insights into the various aspects of differentiation of the SNpr neurons can increase our understanding of fundamental brain functions and their defects in neurological and psychiatric disorders, including movement and mood disorders, as well as epilepsy.

Gata2, Nkx2-2 and Skor2 form a transcription factor network regulating development of a midbrain GABAergic neuron subtype with characteristics of REM-sleep regulatory neurons

The midbrain reticular formation (MRF) is a mosaic of diverse GABAergic and glutamatergic neurons that have been associated with a variety of functions, including sleep regulation. However, the molecular characteristics and development of MRF neurons are poorly understood. As the transcription factor, Gata2 is required for the development of all GABAergic neurons derived from the embryonic mouse midbrain, we hypothesized that the genes expressed downstream of Gata2 could contribute to the diversification of GABAergic neuron subtypes in this brain region. Here, we show that Gata2 is required for the expression of several GABAergic lineage-specific transcription factors, including Nkx2-2 and Skor2, which are co-expressed in a restricted group of post-mitotic GABAergic precursors in the MRF. Both Gata2 and Nkx2-2 function is required for Skor2 expression in GABAergic precursors. In the adult mouse and rat midbrain, Nkx2-2-and Skor2-expressing GABAergic neurons locate at the boundary of the ventrolateral periaqueductal gray and the MRF, an area containing REM-off neurons regulating REM sleep. In addition to the characteristic localization, Skor2+ cells increase their activity upon REM-sleep inhibition, send projections to the dorsolateral pons, a region associated with sleep control, and are responsive to orexins, consistent with the known properties of midbrain REM-off neurons.

Virtual ChIP-seq: predicting transcription factor binding by learning from the transcriptome

Existing methods for computational prediction of transcription factor (TF) binding sites evaluate genomic regions with similarity to known TF sequence preferences. Most TF binding sites, however, do not resemble known TF sequence motifs, and many TFs are not sequence-specific. We developed Virtual ChIP-seq, which predicts binding of individual TFs in new cell types, integrating learned associations with gene expression and binding, TF binding sites from other cell types, and chromatin accessibility data in the new cell type. This approach outperforms methods that predict TF binding solely based on sequence preference, predicting binding for 36 TFs (MCC>0.3).

OUP accepted manuscript

In the last years, many studies were able to identify associations between common genetic variants and complex diseases. However, the mechanistic biological links explaining these associations are still mostly unknown. Common variants are usually associated with a relatively small effect size, suggesting that interactions among multiple variants might be a major genetic component of complex diseases. Hence, elucidating the presence of functional relations among variants may be fundamental to identify putative variants’ interactions. To this aim, we developed Polympact, a web-based resource that allows to explore functional relations among human common variants by exploiting variants’ functional element landscape, their impact on transcription factor binding motifs, and their effect on transcript levels of protein-coding genes. Polympact characterizes over 18 million common variants and allows to explore putative relations by combining clustering analysis and innovative similarity and interaction network models. The properties of the network models were studied and the utility of Polympact was demonstrated by analysing the rich sets of Breast Cancer and Alzheimer's GWAS variants. We identified relations among multiple variants, suggesting putative interactions. Polympact is freely available at bcglab.cibio.unitn.it/polympact.

Satb2 acts as a gatekeeper for major developmental transitions during early vertebrate embryogenesis

Zygotic genome activation (ZGA) initiates regionalized transcription underlying distinct cellular identities. ZGA is dependent upon dynamic chromatin architecture sculpted by conserved DNA-binding proteins. However, the direct mechanistic link between the onset of ZGA and the tissue-specific transcription remains unclear. Here, we have addressed the involvement of chromatin organizer Satb2 in orchestrating both processes during zebrafish embryogenesis. Integrative analysis of transcriptome, genome-wide occupancy and chromatin accessibility reveals contrasting molecular activities of maternally deposited and zygotically synthesized Satb2. Maternal Satb2 prevents premature transcription of zygotic genes by influencing the interplay between the pluripotency factors. By contrast, zygotic Satb2 activates transcription of the same group of genes during neural crest development and organogenesis. Thus, our comparative analysis of maternal versus zygotic function of Satb2 underscores how these antithetical activities are temporally coordinated and functionally implemented highlighting the evolutionary implications of the biphasic and bimodal regulation of landmark developmental transitions by a single determinant.

Activation of the zygotic genome is a critical transition during development, though the link to tissue-specific gene regulation remains unclear. Here the authors demonstrate distinct functions for Satb2 before and after zygotic genome activation, highlighting the temporal coordination of these roles.

New insights into neuropeptides regulation of immune system and hemopoiesis: effects on hematologic malignancies

Several neurotransmitters and neuropeptides were reported to join to or to cooperate with different cells of the immune system, bone marrow, and peripheral cells and numerous data support that neuroactive molecules might control immune system activity and hemopoiesis operating on lymphoid organs, and the primary hematopoietic unit, the hematopoietic niche. Furthermore, many compounds seem to be able to take part to the leukemogenesis and lymphomagenesis process, and in the onset of multiple myeloma. In this review, we will assess the possibility that neurotransmitters and neuropeptides may have a role in the onset of haematological neoplasms, may affect the response to treatment or may represent a useful starting point for a new therapeutic approach. More in vivo investigations are needed to evaluate neuropeptide's role in haematological malignancies and the possible utilization as an antitumor therapeutic target. Comprehending the effect of the pharmacological administration of neuropeptide modulators on hematologic malignancies opens up new possibilities in curing clonal hematologic diseases to achieve more satisfactory outcomes.

Dual midbrain and forebrain origins of thalamic inhibitory interneurons

The ubiquitous presence of inhibitory interneurons in the thalamus of primates contrasts with the sparsity of interneurons reported in mice. Here, we identify a larger than expected complexity and distribution of interneurons across the mouse thalamus, where all thalamic interneurons can be traced back to two developmental programmes: one specified in the midbrain and the other in the forebrain. Interneurons migrate to functionally distinct thalamocortical nuclei depending on their origin: the abundant, midbrain-derived class populates the first and higher order sensory thalamus while the rarer, forebrain-generated class is restricted to some higher order associative regions. We also observe that markers for the midbrain-born class are abundantly expressed throughout the thalamus of the New World monkey marmoset. These data therefore reveal that, despite the broad variability in interneuron density across mammalian species, the blueprint of the ontogenetic organisation of thalamic interneurons of larger-brained mammals exists and can be studied in mice.

A BRD's (BiRD's) eye view of BET and BRPF bromodomains in neurological diseases

Neurological disorders (NLDs) are among the top leading causes for disability worldwide. Dramatic changes in the epigenetic topography of the brain and nervous system have been found in many NLDs. Histone lysine acetylation has prevailed as one of the well characterised epigenetic modifications in these diseases. Two instrumental components of the acetylation machinery are the evolutionarily conserved Bromodomain and PHD finger containing (BRPF) and Bromo and Extra terminal domain (BET) family of proteins, also referred to as acetylation 'readers'. Several reasons, including their distinct mechanisms of modulation of gene expression and their property of being highly tractable small molecule targets, have increased their translational relevance. Thus, compounds which demonstrated promising results in targeting these proteins have advanced to clinical trials. They have been established as key role players in pathologies of cancer, cardiac diseases, renal diseases and rheumatic diseases. In addition, studies implicating the role of these bromodomains in NLDs are gaining pace. In this review, we highlight the findings of these studies, and reason for the plausible roles of all BET and BRPF members in NLDs. A comprehensive understanding of their multifaceted functions would be radical in the development of therapeutic interventions.

The landscape of regulatory genes in brain-wide neuronal phenotypes of a vertebrate brain

Multidimensional landscapes of regulatory genes in neuronal phenotypes at whole-brain levels in the vertebrate remain elusive. We generated single-cell transcriptomes of ~67,000 region- and neurotransmitter/neuromodulator-identifiable cells from larval zebrafish brains. Hierarchical clustering based on effector gene profiles ('terminal features') distinguished major brain cell types. Sister clusters at hierarchical termini displayed similar terminal features. It was further verified by a population-level statistical method. Intriguingly, glutamatergic/GABAergic sister clusters mostly expressed distinct transcription factor (TF) profiles ('convergent pattern'), whereas neuromodulator-type sister clusters predominantly expressed the same TF profiles ('matched pattern'). Interestingly, glutamatergic/GABAergic clusters with similar TF profiles could also display different terminal features ('divergent pattern'). It led us to identify a library of RNA-binding proteins that differentially marked divergent pair clusters, suggesting the post-transcriptional regulation of neuron diversification. Thus, our findings reveal multidimensional landscapes of transcriptional and post-transcriptional regulators in whole-brain neuronal phenotypes in the zebrafish brain.

Human iPSC-Derived Neuronal Cells From CTBP1-Mutated Patients Reveal Altered Expression of Neurodevelopmental Gene Networks

A recurrent de novo mutation in the transcriptional corepressor CTBP1 is associated with neurodevelopmental disabilities in children (Beck et al., 2016, 2019; Sommerville et al., 2017). All reported patients harbor a single recurrent de novo heterozygous missense mutation (p.R342W) within the cofactor recruitment domain of CtBP1. To investigate the transcriptional activity of the pathogenic CTBP1 mutant allele in physiologically relevant human cell models, we generated induced pluripotent stem cells (iPSC) from the dermal fibroblasts derived from patients and normal donors. The transcriptional profiles of the iPSC-derived “early” neurons were determined by RNA-sequencing. Comparison of the RNA-seq data of the neurons from patients and normal donors revealed down regulation of gene networks involved in neurodevelopment, synaptic adhesion and anti-viral (interferon) response. Consistent with the altered gene expression patterns, the patient-derived neurons exhibited morphological and electrophysiological abnormalities, and susceptibility to viral infection. Taken together, our studies using iPSC-derived neuron models provide novel insights into the pathological activities of the CTBP1 p.R342W allele.

ADHD-like behaviors caused by inactivation of a transcription factor controlling the balance of inhibitory and excitatory neuron development in the mouse anterior brainstem

The neural circuits regulating motivation and movement include midbrain dopaminergic neurons and associated inhibitory GABAergic and excitatory glutamatergic neurons in the anterior brainstem. Differentiation of specific subtypes of GABAergic and glutamatergic neurons in the mouse embryonic brainstem is controlled by a transcription factor Tal1. This study characterizes the behavioral and neurochemical changes caused by the absence of Tal1 function. The Tal1^cko mutant mice are hyperactive, impulsive, hypersensitive to reward, have learning deficits and a habituation defect in a novel environment. Only minor changes in their dopaminergic system were detected. Amphetamine induced striatal dopamine release and amphetamine induced place preference were normal in Tal1^cko mice. Increased dopamine signaling failed to stimulate the locomotor activity of the Tal1^cko mice, but instead alleviated their hyperactivity. Altogether, the Tal1^cko mice recapitulate many features of the attention and hyperactivity disorders, suggesting a role for Tal1 regulated developmental pathways and neural structures in the control of motivation and movement.

Expression Quantitative Trait Loci in Equine Skeletal Muscle Reveals Heritable Variation in Metabolism and the Training Responsive Transcriptome

While over ten thousand genetic loci have been associated with phenotypic traits and inherited diseases in genome-wide association studies, in most cases only a relatively small proportion of the trait heritability is explained and biological mechanisms underpinning these traits have not been clearly identified. Expression quantitative trait loci (eQTL) are subsets of genomic loci shown experimentally to influence gene expression. Since gene expression is one of the primary determinants of phenotype, the identification of eQTL may reveal biologically relevant loci and provide functional links between genomic variants, gene expression and ultimately phenotype. Skeletal muscle (gluteus medius) gene expression was quantified by RNA-seq for 111 Thoroughbreds (47 male, 64 female) in race training at a single training establishment sampled at two time-points: at rest (n = 92) and four hours after high-intensity exercise (n = 77); n = 60 were sampled at both time points. Genotypes were generated from the Illumina Equine SNP70 BeadChip. Applying a False Discovery Rate (FDR) corrected P-value threshold (P_FDR < 0.05), association tests identified 3,583 cis-eQTL associated with expression of 1,456 genes at rest; 4,992 cis-eQTL associated with the expression of 1,922 genes post-exercise; 1,703 trans-eQTL associated with 563 genes at rest; and 1,219 trans-eQTL associated with 425 genes post-exercise. The gene with the highest cis-eQTL association at both time-points was the endosome-associated-trafficking regulator 1 gene (ENTR1; Rest: P_FDR = 3.81 × 10^-27, Post-exercise: P_FDR = 1.66 × 10^-24), which has a potential role in the transcriptional regulation of the solute carrier family 2 member 1 glucose transporter protein (SLC2A1). Functional analysis of genes with significant eQTL revealed significant enrichment for cofactor metabolic processes. These results suggest heritable variation in genomic elements such as regulatory sequences (e.g. gene promoters, enhancers, silencers), microRNA and transcription factor genes, which are associated with metabolic function and may have roles in determining end-point muscle and athletic performance phenotypes in Thoroughbred horses. The incorporation of the eQTL identified with genome and transcriptome-wide association may reveal useful biological links between genetic variants and their impact on traits of interest, such as elite racing performance and adaptation to training.

An in-silico approach for discovery of microRNA-TF regulation of DISC1 interactome mediating neuronal migration

Neuronal migration constitutes an important step in corticogenesis; dysregulation of the molecular mechanisms mediating this crucial step in neurodevelopment may result in various neuropsychiatric disorders. By curating experimental data from published literature, we identified eight functional modules involving Disrupted-in-schizophrenia 1 (DISC1) and its interacting proteins that regulate neuronal migration. We then identified miRNAs and transcription factors (TFs) that form functional feedback loops and regulate gene expression of the DISC1 interactome. Using this curated data, we conducted in-silico modeling of the DISC1 interactome involved in neuronal migration and identified the proteins that either facilitate or inhibit neuronal migrational processes. We also studied the effect of perturbation of miRNAs and TFs in feedback loops on the DISC1 interactome. From these analyses, we discovered that STAT3, TCF3, and TAL1 (through feedback loop with miRNAs) play a critical role in the transcriptional control of DISC1 interactome thereby regulating neuronal migration. To the best of our knowledge, regulation of the DISC1 interactome mediating neuronal migration by these TFs has not been previously reported. These potentially important TFs can serve as targets for undertaking validation studies, which in turn can reveal the molecular processes that cause neuronal migration defects underlying neurodevelopmental disorders. This underscores the importance of the use of in-silico techniques in aiding the discovery of mechanistic evidence governing important molecular and cellular processes. The present work is one such step towards the discovery of regulatory factors of the DISC1 interactome that mediates neuronal migration.

Evolution of neuronal types and families

•

Solving nervous system evolution requires cross-species comparison of neuronal types.

•

Neuronal types are commonly defined by their specific structure and function.

•

We provide an operational definition of cell types that allows evolutionary comparison.

•

The identity of neuronal types is best reflected by specifying transcription factors.

•

Families of related neuronal types are conserved across large evolutionary distances.

Major questions in the evolution of neurons and nervous systems remain unsolved, such as the origin of the first neuron, the possible convergent evolution of neuronal phenotypes, and the transition from a relatively simple decentralized nerve net to the complex, centralized nervous systems found in modern bilaterian animals. In recent years, comparative single-cell transcriptomics has opened up new research avenues addressing these issues. Here, we review recent conceptual progress toward an evolutionary definition of cell types, and how it facilitates the identification and large-scale comparison of neuronal types and neuron type families from single-cell data — with the family of GABAergic neurons in distinct parts of the vertebrate forebrain as prime example. We also highlight strategies to infer cell type-specific innovation, so-called apomeres, from single-cell data.

Comprehensive and cell-type-based characterization of the dorsal midbrain during development

The layer structure has been intensively characterized in the developing neocortex and cerebellum based on the various molecular markers. However, as to the developing dorsal midbrain, comprehensive analyses have not been intensely carried out, and thus, the name as well as the definition of each layer is not commonly shared. Here, we redefined the three layers, such as the ventricular zone, intermediate zone and marginal zone, based on various markers for proliferation and differentiation in embryonic dorsal midbrain. Biphasic Ki67 expression defines the classical VZ, in which there is clear separation of the mitotic and interphase zones. Next, we mapped the distribution of immature neurons to the defined layers, based on markers for glutamatergic and GABAergic lineage. Interestingly, Tbr2 and Neurog2 were expressed in the postmitotic neurons. We also report that active (phosphorylated) JNK is a useful marker to demarcate layers during the embryonic stage. Finally, we validated the final arrival layers of the migratory glutamatergic and GABAergic neurons. These results form a foundation for analyses of brain development, especially in the proliferation and migration of excitatory and inhibitory neurons in the dorsal midbrain.

Sequential Role of SOXB2 Factors in GABAergic Neuron Specification of the Dorsal Midbrain

Studies proposed a model for embryonic neurogenesis where the expression levels of the SOXB2 and SOXB1 factors regulate the differentiation status of the neural stem cells. However, the precise role of the SOXB2 genes remains controversial. Therefore, this study aims to investigate the effects of individual deletions of the SOX21 and SOX14 genes during the development of the dorsal midbrain. We show that SOX21 and SOX14 function distinctly during the commitment of the GABAergic lineage. More explicitly, deletion of SOX21 reduced the expression of the GABAergic precursor marker GATA3 and BHLHB5 while the expression of GAD6, which marks GABAergic terminal differentiation, was not affected. In contrast deletion of SOX14 alone was sufficient to inhibit terminal differentiation of the dorsal midbrain GABAergic neurons. Furthermore, we demonstrate through gain-of-function experiments, that despite the homology of SOX21 and SOX14, they have unique gene targets and cannot compensate for the loss of each other. Taken together, these data do not support a pan-neurogenic function for SOXB2 genes in the dorsal midbrain, but instead they influence, sequentially, the specification of GABAergic neurons.

Specification of neurotransmitter identity by Tal1 in thalamic nuclei

BACKGROUND

The neurons contributing to thalamic nuclei are derived from at least two distinct progenitor domains: the caudal (cTH) and rostral (rTH) populations of thalamic progenitors. These neural compartments exhibit unique neurogenic patterns, and the molecular mechanisms underlying the acquisition of neurotransmitter identity remain largely unclear.

RESULTS

T-cell acute lymphocytic leukemia protein 1 (Tal1) was expressed in the early postmitotic cells in the rTH domain, and its expression was maintained in mature thalamic neurons in the ventrolateral geniculate nucleus (vLG) and the intergeniculate leaflet (IGL). To investigate a role of Tal1 in thalamic development, we used a newly generated mouse line driving Cre-mediated recombination in the rTH domain. Conditional deletion of Tal1 did not alter regional patterning in the developing diencephalon. However, in the absence of Tal1, rTH-derived thalamic neurons failed to maintain their postmitotic neuronal features, including neurotransmitter profile. Tal1-deficient thalamic neurons lost their GABAergic markers such as Gad1, Npy, and Penk in IGL/vLG. These defects may be associated at least in part with down-regulation of Nkx2.2, which is known as a critical regulator of rTH-derived GABAergic neurons.

CONCLUSIONS

Our results demonstrate that Tal1 plays an essential role in regulating neurotransmitter phenotype in the developing thalamic nuclei. Developmental Dynamics 246:749-758, 2017. © 2017 Wiley Periodicals, Inc.

Molecular Diversity of Midbrain Development in Mouse, Human, and Stem Cells

Understanding human embryonic ventral midbrain is of major interest for Parkinson's disease. However, the cell types, their gene expression dynamics, and their relationship to commonly used rodent models remain to be defined. We performed single-cell RNA sequencing to examine ventral midbrain development in human and mouse. We found 25 molecularly defined human cell types, including five subtypes of radial glia-like cells and four progenitors. In the mouse, two mature fetal dopaminergic neuron subtypes diversified into five adult classes during postnatal development. Cell types and gene expression were generally conserved across species, but with clear differences in cell proliferation, developmental timing, and dopaminergic neuron development. Additionally, we developed a method to quantitatively assess the fidelity of dopaminergic neurons derived from human pluripotent stem cells, at a single-cell level. Thus, our study provides insight into the molecular programs controlling human midbrain development and provides a foundation for the development of cell replacement therapies.

Expression and Regulation of Tal2 during Neuronal Differentiation in P19 Cells

T-cell acute lymphocytic leukemia 2 (Tal2) is a gene encoding a member of the basic helix-loop-helix transcription factor family, which is essential for the normal development of the mouse brain. We found that Tal2 was induced during neural differentiation in P19 cells, which are pluripotent mouse embryonal carcinoma cells that differentiate into the neural lineage upon both exposure to all-trans retinoic acid (atRA) and the formation of cell aggregation. Tal2 expression during neural differentiation in P19 cells was detected within 3 h after induction with atRA and retinoic acid receptor α (RARα). The atRA-RARα complex is known to bind to a characteristic retinoic acid response element (RARE) located in the promoter of target genes. We found a RARE-like element in the intron of Tal2. We also found a TATA-box-like element in the 5' region. The TATA-box-like element functioned as a core promoter, and TATA- box binding protein bound to this element upstream of Tal2 in P19 cells. The RARE-like element responded to atRA signaling that activated the transcription, and RARα was bound to this element in the intron of Tal2 in P19 cells. Furthermore, the interaction between these elements on Tal2 was confirmed in a chromatin immunoprecipitation assay. Because the neural differentiation of P19 cells mimics in part the development of the nervous system, P19 cells are useful for studying the mechanism underlying the role of Tal2 in neural differentiation. Further work is underway to clarify the function of Tal2 in neural differentiation using the differentiation system of P19 cells.

Reconstructed cell fate-regulatory programs in stem cells reveal hierarchies and key factors of neurogenesis

Cell lineages, which shape the body architecture and specify cell functions, derive from the integration of a plethora of cell intrinsic and extrinsic signals. These signals trigger a multiplicity of decisions at several levels to modulate the activity of dynamic gene regulatory networks (GRNs), which ensure both general and cell-specific functions within a given lineage, thereby establishing cell fates. Significant knowledge about these events and the involved key drivers comes from homogeneous cell differentiation models. Even a single chemical trigger, such as the morphogen all-trans retinoic acid (RA), can induce the complex network of gene-regulatory decisions that matures a stem/precursor cell to a particular step within a given lineage. Here we have dissected the GRNs involved in the RA-induced neuronal or endodermal cell fate specification by integrating dynamic RXRA binding, chromatin accessibility, epigenetic promoter epigenetic status, and the transcriptional activity inferred from RNA polymerase II mapping and transcription profiling. Our data reveal how RA induces a network of transcription factors (TFs), which direct the temporal organization of cognate GRNs, thereby driving neuronal/endodermal cell fate specification. Modeling signal transduction propagation using the reconstructed GRNs indicated critical TFs for neuronal cell fate specification, which were confirmed by CRISPR/Cas9-mediated genome editing. Overall, this study demonstrates that a systems view of cell fate specification combined with computational signal transduction models provides the necessary insight in cellular plasticity for cell fate engineering. The present integrated approach can be used to monitor the in vitro capacity of (engineered) cells/tissues to establish cell lineages for regenerative medicine.

Mesencephalic origin of the rostral Substantia nigra pars reticulata

In embryonic development, the neurons that will constitute a heterogeneous nucleus may have distinct origins. The different components of these populations reach their final location by radial and tangential migrations. The Substantia nigra pars reticulata (SNR) presents a high level of neuronal heterogeneity. It is composed by GABAergic neurons located in the mes-diencephalic basal plate. These inhibitory neurons usually display tangential migrations and it has been already described that the caudal SNR is colonized tangentially from rhombomere 1. Our aim is to unveil the origin of the rostral SNR. We have localized a Nkx6.2 positive ventricular domain located in the alar midbrain. Nkx6.2 derivatives' fate map analysis showed mainly a rostral colonization of this GABAergic neuronal population. We confirmed the mesencephalic origin by the expression of Six3. Both transcription factors are sequentially expressed along the differentiation of these neurons. We demonstrated the origin of the rostral SNR; our data allowed us to postulate that this nucleus is composed by two neuronal populations distributed in opposite gradients with different origins, one from rhombomere 1, caudal to rostral, and the other from the midbrain, rostral to caudal. We can conclude that the SNR has multiple origins and follows complex mechanisms of specification and migration. Our results support vital information for the study of genetic modifications in these extremely complex processes that result in devastating behavioral alterations and predisposition to psychiatric diseases. Understanding the development, molecular identity and functional characteristics of these diverse neuronal populations might lead to better diagnosis and treatment of several forms of neurological and psychiatric disease.

Diversity and development of local inhibitory and excitatory neurons associated with dopaminergic nuclei

For regulation of voluntary movement and motivation the midbrain dopaminergic system receives input from a variety of brain regions. Often this input is mediated by local non-dopaminergic neurons within or closely associated with the dopaminergic nuclei. In addition to the dopaminergic neurons, some of these non-dopaminergic neurons also send functionally important output from the ventral midbrain to forebrain targets. The aim of this review is to introduce subtypes of GABAergic and glutamatergic neurons, which are located in the dopaminergic nuclei or the adjacent brainstem and are important for the regulation of the dopaminergic pathways. In addition, we discuss recent studies beginning to reveal mechanisms for their development, which may hold the key to understanding the diversity of these neurons.

Hairy/Enhancer-of-Split MEGANE and Proneural MASH1 Factors Cooperate Synergistically in Midbrain GABAergic Neurogenesis

GABAergic neurons are the primary inhibitory cell type in the mature brain and their dysfunction is associated with important neurological conditions like schizophrenia and anxiety. We aimed to discover the underlying mechanisms for dorsal/ventral midbrain GABAergic neurogenesis. Previous work by us and others has provided crucial insights into the key function of Mgn and Mash1 genes in determining GABAergic neurotransmitter fate. Induction of dorsal midbrain GABAergic neurons does not take place at any time during development in either of the single mutant mice. However, GABAergic neurons in the ventral midbrain remained unchanged. Thus, the similarities in MB-GABAergic phenotype observed in the Mgn and Mash1 single mutants suggest the existence of other factors that take over the function of MGN and MASH1 in the ventral midbrain or the existence of different molecular mechanisms. We show that this process essentially depends on heterodimers and homodimers formed by MGN and MASH1 and deciphered the in vivo relevance of the interaction by phenotypic analysis of Mgn/Mash1 double knockout and compound mice. Furthermore, the combination of gain- and loss-of-function experiments in the developing midbrain showed co-operative roles for Mgn and Mash1 genes in determining GABAergic identity. Transcription factors belonging to the Enhancer-of-split-related and proneural families have long been believed to counterpart each other's function. This work uncovers a synergistic cooperation between these two families, and provides a novel paradigm for how these two families cooperate for the acquisition of MB-GABAergic neuronal identity. Understanding their molecular mechanisms is essential for cell therapy strategies to amend GABAergic deficits.

Occupancy by key transcription factors is a more accurate predictor of enhancer activity than histone modifications or chromatin accessibility

Background

Regulated gene expression controls organismal development, and variation in regulatory patterns has been implicated in complex traits. Thus accurate prediction of enhancers is important for further understanding of these processes. Genome-wide measurement of epigenetic features, such as histone modifications and occupancy by transcription factors, is improving enhancer predictions, but the contribution of these features to prediction accuracy is not known. Given the importance of the hematopoietic transcription factor TAL1 for erythroid gene activation, we predicted candidate enhancers based on genomic occupancy by TAL1 and measured their activity. Contributions of multiple features to enhancer prediction were evaluated based on the results of these and other studies.

Results

TAL1-bound DNA segments were active enhancers at a high rate both in transient transfections of cultured cells (39 of 79, or 56%) and transgenic mice (43 of 66, or 65%). The level of binding signal for TAL1 or GATA1 did not help distinguish TAL1-bound DNA segments as active versus inactive enhancers, nor did the density of regulation-related histone modifications. A meta-analysis of results from this and other studies (273 tested predicted enhancers) showed that the presence of TAL1, GATA1, EP300, SMAD1, H3K4 methylation, H3K27ac, and CAGE tags at DNase hypersensitive sites gave the most accurate predictors of enhancer activity, with a success rate over 80% and a median threefold increase in activity. Chromatin accessibility assays and the histone modifications H3K4me1 and H3K27ac were sensitive for finding enhancers, but they have high false positive rates unless transcription factor occupancy is also included.

Conclusions

Occupancy by key transcription factors such as TAL1, GATA1, SMAD1, and EP300, along with evidence of transcription, improves the accuracy of enhancer predictions based on epigenetic features.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0009-5) contains supplementary material, which is available to authorized users.

Differentiation and molecular heterogeneity of inhibitory and excitatory neurons associated with midbrain dopaminergic nuclei

Local inhibitory GABAergic and excitatory glutamatergic neurons are important for midbrain dopaminergic and hindbrain serotonergic pathways controlling motivation, mood, and voluntary movements. Such neurons reside both within the dopaminergic nuclei, and in adjacent brain structures, including the rostromedial and laterodorsal tegmental nuclei. Compared to the monoaminergic neurons, the development, heterogeneity, and molecular characteristics of these regulatory neurons are poorly understood. We show here that different GABAergic and glutamatergic subgroups associated with the monoaminergic nuclei express specific transcription factors. These neurons share common origins in the ventrolateral rhombomere 1, where postmitotic selector genes Tal1, Gata2, and Gata3 control the balance between the generation of inhibitory and excitatory neurons. In the absence of Tal1, or both Gata2 and Gata3, the GABAergic precursors adopt glutamatergic fates and populate the glutamatergic nuclei in excessive numbers. Together, our results uncover developmental regulatory mechanisms, molecular characteristics, and heterogeneity of central regulators of monoaminergic circuits.

Dynamic shifts in occupancy by TAL1 are guided by GATA factors and drive large-scale reprogramming of gene expression during hematopoiesis

We used mouse ENCODE data along with complementary data from other laboratories to study the dynamics of occupancy and the role in gene regulation of the transcription factor TAL1, a critical regulator of hematopoiesis, at multiple stages of hematopoietic differentiation. We combined ChIP-seq and RNA-seq data in six mouse cell types representing a progression from multilineage precursors to differentiated erythroblasts and megakaryocytes. We found that sites of occupancy shift dramatically during commitment to the erythroid lineage, vary further during terminal maturation, and are strongly associated with changes in gene expression. In multilineage progenitors, the likely target genes are enriched for hematopoietic growth and functions associated with the mature cells of specific daughter lineages (such as megakaryocytes). In contrast, target genes in erythroblasts are specifically enriched for red cell functions. Furthermore, shifts in TAL1 occupancy during erythroid differentiation are associated with gene repression (dissociation) and induction (co-occupancy with GATA1). Based on both enrichment for transcription factor binding site motifs and co-occupancy determined by ChIP-seq, recruitment by GATA transcription factors appears to be a stronger determinant of TAL1 binding to chromatin than the canonical E-box binding site motif. Studies of additional proteins lead to the model that TAL1 regulates expression after being directed to a distinct subset of genomic binding sites in each cell type via its association with different complexes containing master regulators such as GATA2, ERG, and RUNX1 in multilineage cells and the lineage-specific master regulator GATA1 in erythroblasts.

Mechanisms regulating GABAergic neuron development

Neurons using gamma-aminobutyric acid (GABA) as their neurotransmitter are the main inhibitory neurons in the mature central nervous system (CNS) and show great variation in their form and function. GABAergic neurons are produced in all of the main domains of the CNS, where they develop from discrete regions of the neuroepithelium. Here, we review the gene expression and regulatory mechanisms controlling the main steps of GABAergic neuron development: early patterning of the proliferative neuroepithelium, production of postmitotic neural precursors, establishment of their identity and migration. By comparing the molecular regulation of these events across CNS, we broadly identify three regions utilizing distinct molecular toolkits for GABAergic fate determination: telencephalon-anterior diencephalon (DLX2 type), posterior diencephalon-midbrain (GATA2 type) and hindbrain-spinal cord (PTF1A and TAL1 types). Similarities and differences in the molecular regulatory mechanisms reveal the core determinants of a GABAergic neuron as well as provide insights into generation of the vast diversity of these neurons.