Repressive interactions in gene regulatory networks: When you have no other choice

Activation-derepression synergy enables a bHLH network to coordinate a signal-specific fate response

Stem cells integrate multiple environmental signals to activate appropriate fate programs. To ensure coherent responses, alternative fates must be concomitantly inactivated. However, mechanisms that coordinate fates in a signal-specific manner are not fully understood. Here, we investigate the role of a network of basic-helix-loop-helix (bHLH) transcription factors in neural stem cells, which integrate leukemia inhibitory factor (LIF) and bone morphogenetic protein (BMP) signaling to synergistically induce glial fibrillary acidic protein (GFAP), a key astrocyte-fate determinant. Using quantitative RNA-fluorescence in situ hybridization (FISH) and ectopic expression, we find that multiple bHLHs that promote alternative fates also repress GFAP but are all suppressed by BMP and, to a lesser extent, LIF. Mathematical modeling shows that synergy arises from this coordinated derepression of GFAP combined with its activation by LIF signaling. Finally, we determine how coordinated and tunable derepression results from extensive cross-regulation among bHLHs. Activation-derepression synergy could be broadly utilized to couple signaling and fate, particularly across the numerous developmental systems regulated by bHLH factors.

Lhx4 surpasses its paralog Lhx3 in promoting the differentiation of spinal V2a interneurons

Paralog factors are considered to ensure the robustness of biological processes by providing redundant activity in cells where they are co-expressed. However, the specific contribution of each factor is frequently underestimated. In the developing spinal cord, multiple families of transcription factors successively contribute to differentiate an initially homogenous population of neural progenitors into a myriad of neuronal subsets with distinct molecular, morphological, and functional characteristics. The LIM-homeodomain transcription factors Lhx3, Lhx4, Isl1 and Isl2 promote the segregation and differentiation of spinal motor neurons and V2 interneurons. Based on their high sequence identity and their similar distribution, the Lhx3 and Lhx4 paralogs are considered to contribute similarly to these processes. However, the specific contribution of Lhx4 has never been studied. Here, we provide evidence that Lhx3 and Lhx4 are present in the same cell populations during spinal cord development. Similarly to Lhx3, Lhx4 can form multiproteic complexes with Isl1 or Isl2 and the nuclear LIM interactor NLI. Lhx4 can stimulate a V2-specific enhancer more efficiently than Lhx3 and surpasses Lhx3 in promoting the differentiation of V2a interneurons in chicken embryo electroporation experiments. Finally, Lhx4 inactivation in mice results in alterations of differentiation of the V2a subpopulation, but not of motor neuron production, suggesting that Lhx4 plays unique roles in V2a differentiation that are not compensated by the presence of Lhx3. Thus, Lhx4 could be the major LIM-HD factor involved in V2a interneuron differentiation during spinal cord development and should be considered for in vitro differentiation of spinal neuronal populations.

Evolution and development of extraocular motor neurons, nerves and muscles in vertebrates

The purpose of this review is to analyze the origin of ocular motor neurons, define the pattern of innervation of nerve fibers that project to the extraocular eye muscles (EOMs), describe congenital disorders that alter the development of ocular motor neurons, and provide an overview of vestibular pathway inputs to ocular motor nuclei. Six eye muscles are innervated by axons of three ocular motor neurons, the oculomotor (CNIII), trochlear (CNIV), and abducens (CNVI) neurons. Ocular motor neurons (CNIII) originate in the midbrain and innervate the ipsilateral orbit, except for the superior rectus and the levator palpebrae, which are contralaterally innervated. Trochlear motor neurons (CNIV) originate at the midbrain-hindbrain junction and innervate the contralateral superior oblique muscle. Abducens motor neurons (CNVI) originate variously in the hindbrain of rhombomeres r4-6 that innervate the posterior (or lateral) rectus muscle and innervate the retractor bulbi. Genes allow a distinction between special somatic (CNIII, IV) and somatic (CNVI) ocular motor neurons. Development of ocular motor neurons and their axonal projections to the EOMs may be derailed by various genetic causes, resulting in the congenital cranial dysinnervation disorders. The ocular motor neurons innervate EOMs while the vestibular nuclei connect with the midbrain-brainstem motor neurons.

Gene networks and the evolution of olfactory organs, eyes, hair cells and motoneurons: a view encompassing lancelets, tunicates and vertebrates

Key developmental pathways and gene networks underlie the formation of sensory cell types and structures involved in chemosensation, vision and mechanosensation, and of the efferents these sensory inputs can activate. We describe similarities and differences in these pathways and gene networks in selected species of the three main chordate groups, lancelets, tunicates, and vertebrates, leading to divergent development of olfactory receptors, eyes, hair cells and motoneurons. The lack of appropriately posited expression of certain transcription factors in lancelets and tunicates prevents them from developing vertebrate-like olfactory receptors and eyes, although they generate alternative structures for chemosensation and vision. Lancelets and tunicates lack mechanosensory cells associated with the sensation of acoustic stimuli, but have gravisensitive organs and ciliated epidermal sensory cells that may (and in some cases clearly do) provide mechanosensation and thus the capacity to respond to movement relative to surrounding water. Although functionally analogous to the vertebrate vestibular apparatus and lateral line, homology is questionable due to differences in the expression of the key transcription factors Neurog and Atoh1/7, on which development of vertebrate hair cells depends. The vertebrate hair cell-bearing inner ear and lateral line thus likely represent major evolutionary advances specific to vertebrates. Motoneurons develop in vertebrates under the control of the ventral signaling molecule hedgehog/sonic hedgehog (Hh,Shh), against an opposing inhibitory effect mediated by dorsal signaling molecules. Many elements of Shh-signaling and downstream genes involved in specifying and differentiating motoneurons are also exhibited by lancelets and tunicates, but the repertoire of MNs in vertebrates is broader, indicating greater diversity in motoneuron differentiation programs.

From signalling to form: the coordination of neural tube patterning

The development of the vertebrate spinal cord involves the formation of the neural tube and the generation of multiple distinct cell types. The process starts during gastrulation, combining axial elongation with specification of neural cells and the formation of the neuroepithelium. Tissue movements produce the neural tube which is then exposed to signals that provide patterning information to neural progenitors. The intracellular response to these signals, via a gene regulatory network, governs the spatial and temporal differentiation of progenitors into specific cell types, facilitating the assembly of functional neuronal circuits. The interplay between the gene regulatory network, cell movement, and tissue mechanics generates the conserved neural tube pattern observed across species. In this review we offer an overview of the molecular and cellular processes governing the formation and patterning of the neural tube, highlighting how the remarkable complexity and precision of vertebrate nervous system arises. We argue that a multidisciplinary and multiscale understanding of the neural tube development, paired with the study of species-specific strategies, will be crucial to tackle the open questions.

Molecular analyses of zebrafish V0v spinal interneurons and identification of transcriptional regulators downstream of Evx1 and Evx2 in these cells

BACKGROUND

V0v spinal interneurons are highly conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells.

METHODS

To identify candidate members of V0v gene regulatory networks, we FAC-sorted wild-type and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

RESULTS

Our data reveal two molecularly distinct subtypes of zebrafish V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuron expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression.

CONCLUSIONS

This study identifies two molecularly distinct subsets of zebrafish V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of these cells. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

The Drosophila ecdysone receptor promotes or suppresses proliferation according to ligand level

The steroid hormone 20-hydroxy-ecdysone (20E) promotes proliferation in Drosophila wing precursors at low titer but triggers proliferation arrest at high doses. Remarkably, wing precursors proliferate normally in the complete absence of the 20E receptor, suggesting that low-level 20E promotes proliferation by overriding the default anti-proliferative activity of the receptor. By contrast, 20E needs its receptor to arrest proliferation. Dose-response RNA sequencing (RNA-seq) analysis of ex vivo cultured wing precursors identifies genes that are quantitatively activated by 20E across the physiological range, likely comprising positive modulators of proliferation and other genes that are only activated at high doses. We suggest that some of these "high-threshold" genes dominantly suppress the activity of the pro-proliferation genes. We then show mathematically and with synthetic reporters that combinations of basic regulatory elements can recapitulate the behavior of both types of target genes. Thus, a relatively simple genetic circuit can account for the bimodal activity of this hormone.

Molecular Analyses of V0v Spinal Interneurons and Identification of Transcriptional Regulators Downstream of Evx1 and Evx2 in these Cells

Background

V0v spinal interneurons are highly conserved, glutamatergic, commissural neurons that function in locomotor circuits. We have previously shown that Evx1 and Evx2 are required to specify the neurotransmitter phenotype of these cells. However, we still know very little about the gene regulatory networks that act downstream of these transcription factors in V0v cells.

Methods

To identify candidate members of V0v gene regulatory networks, we FAC-sorted WT and evx1;evx2 double mutant zebrafish V0v spinal interneurons and expression-profiled them using microarrays and single cell RNA-seq. We also used in situ hybridization to compare expression of a subset of candidate genes in evx1;evx2 double mutants and wild-type siblings.

Results

Our data reveal two molecularly distinct subtypes of V0v spinal interneurons at 48 h and suggest that, by this stage of development, evx1;evx2 double mutant cells transfate into either inhibitory spinal interneurons, or motoneurons. Our results also identify 25 transcriptional regulator genes that require Evx1/2 for their expression in V0v interneurons, plus a further 11 transcriptional regulator genes that are repressed in V0v interneurons by Evx1/2. Two of the latter genes are hmx2 and hmx3a. Intriguingly, we show that Hmx2/3a, repress dI2 interneuronal expression of skor1a and nefma, two genes that require Evx1/2 for their expression in V0v interneurons. This suggests that Evx1/2 might regulate skor1a and nefma expression in V0v interneurons by repressing Hmx2/3a expression.

Conclusions

This study identifies two molecularly distinct subsets of V0v spinal interneurons, as well as multiple transcriptional regulators that are strong candidates for acting downstream of Evx1/2 to specify the essential functional characteristics of these cells. Our data further suggest that in the absence of both Evx1 and Evx2, V0v spinal interneurons initially change their neurotransmitter phenotypes from excitatory to inhibitory and then, later, start to express markers of distinct types of inhibitory spinal interneurons, or motoneurons. Taken together, our findings significantly increase our knowledge of V0v and spinal development and move us closer towards the essential goal of identifying the complete gene regulatory networks that specify this crucial cell type.

Mammalian olfactory cortex neurons retain molecular signatures of ancestral cell types

The cerebral cortex diversified extensively during vertebrate evolution. Intriguingly, the three-layered mammalian olfactory cortex resembles the cortical cytoarchitecture of non-mammals yet evolved alongside the six-layered neocortex, enabling unique comparisons for investigating cortical neuron diversification. We performed single-nucleus multiome sequencing across mouse three- to six-layered cortices and compared neuron types across mice, reptiles and salamander. We identified neurons that are olfactory cortex-specific or conserved across mouse cortical areas. However, transcriptomically similar neurons exhibited area-specific epigenetic states. Additionally, the olfactory cortex showed transcriptomic divergence between lab and wild-derived mice, suggesting enhanced circuit plasticity through adult immature neurons. Finally, olfactory cortex neurons displayed marked transcriptomic similarities to reptile and salamander neurons. Together, these data indicate that the mammalian olfactory cortex retains molecular signatures representative of ancestral cortical traits.

Transcriptome profiling of the Caenorhabditis elegans intestine reveals that ELT-2 negatively and positively regulates intestinal gene expression within the context of a gene regulatory network

ELT-2 is the major transcription factor (TF) required for Caenorhabditis elegans intestinal development. ELT-2 expression initiates in embryos to promote development and then persists after hatching through the larval and adult stages. Though the sites of ELT-2 binding are characterized and the transcriptional changes that result from ELT-2 depletion are known, an intestine-specific transcriptome profile spanning developmental time has been missing. We generated this dataset by performing Fluorescence Activated Cell Sorting on intestine cells at distinct developmental stages. We analyzed this dataset in conjunction with previously conducted ELT-2 studies to evaluate the role of ELT-2 in directing the intestinal gene regulatory network through development. We found that only 33% of intestine-enriched genes in the embryo were direct targets of ELT-2 but that number increased to 75% by the L3 stage. This suggests additional TFs promote intestinal transcription especially in the embryo. Furthermore, only half of ELT-2's direct target genes were dependent on ELT-2 for their proper expression levels, and an equal proportion of those responded to elt-2 depletion with over-expression as with under-expression. That is, ELT-2 can either activate or repress direct target genes. Additionally, we observed that ELT-2 repressed its own promoter, implicating new models for its autoregulation. Together, our results illustrate that ELT-2 impacts roughly 20-50% of intestine-specific genes, that ELT-2 both positively and negatively controls its direct targets, and that the current model of the intestinal regulatory network is incomplete as the factors responsible for directing the expression of many intestinal genes remain unknown.

The time integral of BMP signaling determines fate in a stem cell model for early human development

How paracrine signals are interpreted to yield multiple cell fate decisions in a dynamic context during human development in vivo and in vitro remains poorly understood. Here we report an automated tracking method to follow signaling histories linked to cell fate in large numbers of human pluripotent stem cells (hPSCs). Using an unbiased statistical approach, we discovered that measured BMP signaling history correlates strongly with fate in individual cells. We found that BMP response in hPSCs varies more strongly in the duration of signaling than the level. However, we discovered that both the level and duration of signaling activity control cell fate choices only by changing the time integral of signaling and that duration and level are therefore interchangeable in this context. In a stem cell model for patterning of the human embryo, we showed that signaling histories predict the fate pattern and that the integral model correctly predicts changes in cell fate domains when signaling is perturbed. Using an RNA-seq screen we then found that mechanistically, BMP signaling is integrated by SOX2.

Developmental cell fate choice in neural tube progenitors employs two distinct cis-regulatory strategies

In many developing tissues, the patterns of gene expression that assign cell fate are organized by graded secreted signals. Cis-regulatory elements (CREs) interpret these signals to control gene expression, but how this is accomplished remains poorly understood. In the neural tube, a gradient of the morphogen sonic hedgehog (Shh) patterns neural progenitors. We identify two distinct ways in which CREs translate graded Shh into differential gene expression in mouse neural progenitors. In most progenitors, a common set of CREs control gene activity by integrating cell-type-specific inputs. By contrast, the most ventral progenitors use a unique set of CREs, established by the pioneer factor FOXA2. This parallels the role of FOXA2 in endoderm, where FOXA2 binds some of the same sites. Together, the data identify distinct cis-regulatory strategies for the interpretation of morphogen signaling and raise the possibility of an evolutionarily conserved role for FOXA2 across tissues.

Learning in Transcriptional Network Models: Computational Discovery of Pathway-Level Memory and Effective Interventions

Trainability, in any substrate, refers to the ability to change future behavior based on past experiences. An understanding of such capacity within biological cells and tissues would enable a particularly powerful set of methods for prediction and control of their behavior through specific patterns of stimuli. This top-down mode of control (as an alternative to bottom-up modification of hardware) has been extensively exploited by computer science and the behavioral sciences; in biology however, it is usually reserved for organism-level behavior in animals with brains, such as training animals towards a desired response. Exciting work in the field of basal cognition has begun to reveal degrees and forms of unconventional memory in non-neural tissues and even in subcellular biochemical dynamics. Here, we characterize biological gene regulatory circuit models and protein pathways and find them capable of several different kinds of memory. We extend prior results on learning in binary transcriptional networks to continuous models and identify specific interventions (regimes of stimulation, as opposed to network rewiring) that abolish undesirable network behavior such as drug pharmacoresistance and drug sensitization. We also explore the stability of created memories by assessing their long-term behavior and find that most memories do not decay over long time periods. Additionally, we find that the memory properties are quite robust to noise; surprisingly, in many cases noise actually increases memory potential. We examine various network properties associated with these behaviors and find that no one network property is indicative of memory. Random networks do not show similar memory behavior as models of biological processes, indicating that generic network dynamics are not solely responsible for trainability. Rational control of dynamic pathway function using stimuli derived from computational models opens the door to empirical studies of proto-cognitive capacities in unconventional embodiments and suggests numerous possible applications in biomedicine, where behavior shaping of pathway responses stand as a potential alternative to gene therapy.

Molecular Organization and Patterning of the Medulla Oblongata in Health and Disease

The medulla oblongata, located in the hindbrain between the pons and the spinal cord, is an important relay center for critical sensory, proprioceptive, and motoric information. It is an evolutionarily highly conserved brain region, both structural and functional, and consists of a multitude of nuclei all involved in different aspects of basic but vital functions. Understanding the functional anatomy and developmental program of this structure can help elucidate potential role(s) of the medulla in neurological disorders. Here, we have described the early molecular patterning of the medulla during murine development, from the fundamental units that structure the very early medullary region into 5 rhombomeres (r7–r11) and 13 different longitudinal progenitor domains, to the neuronal clusters derived from these progenitors that ultimately make-up the different medullary nuclei. By doing so, we developed a schematic overview that can be used to predict the cell-fate of a progenitor group, or pinpoint the progenitor domain of origin of medullary nuclei. This schematic overview can further be used to help in the explanation of medulla-related symptoms of neurodevelopmental disorders, e.g., congenital central hypoventilation syndrome, Wold–Hirschhorn syndrome, Rett syndrome, and Pitt–Hopkins syndrome. Based on the genetic defects seen in these syndromes, we can use our model to predict which medullary nuclei might be affected, which can be used to quickly direct the research into these diseases to the likely affected nuclei.

Lessons from a transcription factor: Alx1 provides insights into gene regulatory networks, cellular reprogramming, and cell type evolution

The skeleton-forming cells of sea urchins and other echinoderms have been studied by developmental biologists as models of cell specification and morphogenesis for many decades. The gene regulatory network (GRN) deployed in the embryonic skeletogenic cells of euechinoid sea urchins is one of the best understood in any developing animal. Recent comparative studies have leveraged the information contained in this GRN, bringing renewed attention to the diverse patterns of skeletogenesis within the phylum and the evolutionary basis for this diversity. The homeodomain-containing transcription factor, Alx1, was originally shown to be a core component of the skeletogenic GRN of the sea urchin embryo. Alx1 has since been found to be key regulator of skeletal cell identity throughout the phylum. As such, Alx1 is currently serving as a lens through which multiple developmental processes are being investigated. These include not only GRN organization and evolution, but also cell reprogramming, cell type evolution, and the gene regulatory control of morphogenesis. This review summarizes our current state of knowledge concerning Alx1 and highlights the insights it is yielding into these important developmental and evolutionary processes.

The developmental hourglass model is applicable to the spinal cord based on single-cell transcriptomes and non-conserved cis-regulatory elements

The developmental hourglass model predicts that embryonic morphology is most conserved at the mid-embryonic stage and diverges at the early and late stages. To date, this model has been verified by examining the anatomical features or gene expression profiles at the whole embryonic level. Here, by data mining approach utilizing multiple genomic and transcriptomic datasets from different species in combination, and by experimental validation, we demonstrate that the hourglass model is also applicable to a reduced element, the spinal cord. In the middle of spinal cord development, dorsoventrally arrayed neuronal progenitor domains are established, which are conserved among vertebrates. By comparing the publicly available single-cell transcriptome datasets of mice and zebrafish, we found that ventral subpopulations of post-mitotic spinal neurons display divergent molecular profiles. We also detected the non-conservation of cis-regulatory elements located around the progenitor fate determinants, indicating that the cis-regulatory elements contributing to the progenitor specification are evolvable. These results demonstrate that, despite the conservation of the progenitor domains, the processes before and after the progenitor domain specification diverged. This study will be helpful to understand the molecular basis of the developmental hourglass model.

Single-cell transcriptome profiling of the human developing spinal cord reveals a conserved genetic programme with human-specific features

The spinal cord receives input from peripheral sensory neurons and controls motor output by regulating muscle innervating motor neurons. These functions are carried out by neural circuits comprising molecularly distinct neuronal subtypes generated in a characteristic spatiotemporal arrangement from progenitors in the embryonic neural tube. To gain insight into the diversity and complexity of cells in the developing human neural tube, we used single-cell mRNA sequencing to profile cervical and thoracic regions in four human embryos of Carnegie stages (CS) CS12, CS14, CS17 and CS19 from gestational weeks 4-7. Analysis of progenitor and neuronal populations from the neural tube and dorsal root ganglia identified dozens of distinct cell types and facilitated the reconstruction of the differentiation pathways of specific neuronal subtypes. Comparison with mouse revealed overall similarity of mammalian neural tube development while highlighting some human-specific features. These data provide a catalogue of gene expression and cell type identity in the human neural tube that will support future studies of sensory and motor control systems. The data can be explored at https://shiny.crick.ac.uk/scviewer/neuraltube/.

Cross-species comparisons and in vitro models to study tempo in development and homeostasis

Time is inherent to biological processes. It determines the order of events and the speed at which they take place. However, we still need to refine approaches to measure the course of time in biological systems and understand what controls the pace of development. Here, we argue that the comparison of biological processes across species provides molecular insight into the timekeeping mechanisms in biology. We discuss recent findings and the open questions in the field and highlight the use of in vitro systems as tools to investigate cell-autonomous control as well as the coordination of temporal mechanisms within tissues. Further, we discuss the relevance of studying tempo for tissue transplantation, homeostasis and lifespan.

MicroRNA governs bistable cell differentiation and lineage segregation via a noncanonical feedback

Positive feedback driven by transcriptional regulation has long been considered a key mechanism underlying cell lineage segregation during embryogenesis. Using the developing spinal cord as a paradigm, we found that canonical, transcription-driven feedback cannot explain robust lineage segregation of motor neuron subtypes marked by two cardinal factors, Hoxa5 and Hoxc8. We propose a feedback mechanism involving elementary microRNA-mRNA reaction circuits that differ from known feedback loop-like structures. Strikingly, we show that a wide range of biologically plausible post-transcriptional regulatory parameters are sufficient to generate bistable switches, a hallmark of positive feedback. Through mathematical analysis, we explain intuitively the hidden source of this feedback. Using embryonic stem cell differentiation and mouse genetics, we corroborate that microRNA-mRNA circuits govern tissue boundaries and hysteresis upon motor neuron differentiation with respect to transient morphogen signals. Our findings reveal a previously underappreciated feedback mechanism that may have widespread functions in cell fate decisions and tissue patterning.