Cell interactions, signals and transcriptional hierarchy governing placode progenitor induction

PMF-GRN: a variational inference approach to single-cell gene regulatory network inference using probabilistic matrix factorization

Inferring gene regulatory networks (GRNs) from single-cell data is challenging due to heuristic limitations. Existing methods also lack estimates of uncertainty. Here we present Probabilistic Matrix Factorization for Gene Regulatory Network Inference (PMF-GRN). Using single-cell expression data, PMF-GRN infers latent factors capturing transcription factor activity and regulatory relationships. Using variational inference allows hyperparameter search for principled model selection and direct comparison to other generative models. We extensively test and benchmark our method using real single-cell datasets and synthetic data. We show that PMF-GRN infers GRNs more accurately than current state-of-the-art single-cell GRN inference methods, offering well-calibrated uncertainty estimates.

Zmym4 is required for early cranial gene expression and craniofacial cartilage formation

Introduction: The Six1 transcription factor plays important roles in the development of cranial sensory organs, and point mutations underlie craniofacial birth defects. Because Six1's transcriptional activity can be modulated by interacting proteins, we previously screened for candidate interactors and identified zinc-finger MYM-containing protein 4 (Zmym4) by its inclusion of a few domains with a bona fide cofactor, Sine oculis binding protein (Sobp). Although Zmym4 has been implicated in regulating early brain development and certain cancers, its role in craniofacial development has not previously been described. Methods: We used co-immunoprecipitation and luciferase-reporter assays in cultured cells to test interactions between Zmym4 and Six1. We used knock-down and overexpression of Zmym4 in embryos to test for its effects on early ectodermal gene expression, neural crest migration and craniofacial cartilage formation. Results: We found no evidence that Zmym4 physically or transcriptionally interacts with Six1 in cultured cells. Nonetheless, knockdown of endogenous Zmym4 in embryos resulted in altered early cranial gene expression, including those expressed in the neural border, neural plate, neural crest and preplacodal ectoderm. Experimentally increasing Zmym4 levels had minor effects on neural border or neural plate genes, but altered the expression of neural crest and preplacodal genes. At larval stages, genes expressed in the otic vesicle and branchial arches showed reduced expression in Zmym4 morphants. Although we did not detect defects in neural crest migration into the branchial arches, loss of Zmym4 resulted in aberrant morphology of several craniofacial cartilages. Discussion: Although Zmym4 does not appear to function as a Six1 transcriptional cofactor, it plays an important role in regulating the expression of embryonic cranial genes in tissues critical for normal craniofacial development.

scRNA-sequencing in chick suggests a probabilistic model for cell fate allocation at the neural plate border

The vertebrate 'neural plate border' is a transient territory located at the edge of the neural plate containing precursors for all ectodermal derivatives: the neural plate, neural crest, placodes and epidermis. Elegant functional experiments in a range of vertebrate models have provided an in-depth understanding of gene regulatory interactions within the ectoderm. However, these experiments conducted at tissue level raise seemingly contradictory models for fate allocation of individual cells. Here, we carry out single cell RNA sequencing of chick ectoderm from primitive streak to neurulation stage, to explore cell state diversity and heterogeneity. We characterise the dynamics of gene modules, allowing us to model the order of molecular events which take place as ectodermal fates segregate. Furthermore, we find that genes previously classified as neural plate border 'specifiers' typically exhibit dynamic expression patterns and are enriched in either neural, neural crest or placodal fates, revealing that the neural plate border should be seen as a heterogeneous ectodermal territory and not a discrete transitional transcriptional state. Analysis of neural, neural crest and placodal markers reveals that individual NPB cells co-express competing transcriptional programmes suggesting that their ultimate identify is not yet fixed. This population of 'border located undecided progenitors' (BLUPs) gradually diminishes as cell fate decisions take place. Considering our findings, we propose a probabilistic model for cell fate choice at the neural plate border. Our data suggest that the probability of a progenitor's daughters to contribute to a given ectodermal derivative is related to the balance of competing transcriptional programmes, which in turn are regulated by the spatiotemporal position of a progenitor.

Anatomical network modules of the human central nervous-craniofacial skeleton system

Anatomical network analysis (AnNA) is a systems biological framework based on network theory that enables anatomical structural analysis by incorporating modularity to model structural complexity. The human brain and facial structures exhibit close structural and functional relationships, suggestive of a co-evolved anatomical network. The present study aimed to analyze the human head as a modular entity that comprises the central nervous system, including the brain, spinal cord, and craniofacial skeleton. An AnNA model was built using 39 anatomical nodes from the brain, spinal cord, and craniofacial skeleton. The linkages were identified using peripheral nerve supply and direct contact between structures. The Spinglass algorithm in the igraph software was applied to construct a network and identify the modules of the central nervous system-craniofacial skeleton anatomical network. Two modules were identified. These comprised an anterior module, which included the forebrain, anterior cranial base, and upper-middle face, and a posterior module, which included the midbrain, hindbrain, mandible, and posterior cranium. These findings may reflect the genetic and signaling networks that drive the mosaic central nervous system and craniofacial development and offer important systems biology perspectives for developmental disorders of craniofacial structures.

Folate Carrier Deficiency Drives Differential Methylation and Enhanced Cellular Potency in the Neural Plate Border

The neural plate border (NPB) of vertebrate embryos segregates from the neural and epidermal regions, and it is comprised of an intermingled group of multipotent progenitor cells. Folate is the precursor of S-adenosylmethionine, the main methyl donor for DNA methylation, and it is critical for embryonic development, including the specification of progenitors which reside in the NPB. Despite the fact that several intersecting signals involved in the specification and territorial restriction of NPB cells are known, the role of epigenetics, particularly DNA methylation, has been a matter of debate. Here, we examined the temporal and spatial distribution of the methyl source and analyzed the abundance of 5mC/5 hmC and their epigenetic writers throughout the segregation of the neural and NPB territories. Reduced representation bisulfite sequencing (RRBS) on Reduced Folate Carrier 1 (RFC1)-deficient embryos leads to the identification of differentially methylated regions (DMRs). In the RFC1-deficient embryos, we identified several DMRs in the Notch1 locus, and the spatiotemporal expression of Notch1 and its downstream target gene Bmp4 were expanded in the NPB. Cell fate analysis on folate deficient embryos revealed a significant increase in the number of cells coexpressing both neural (SOX2) and NPB (PAX7) markers, which may represent an enhancing effect in the cellular potential of those progenitors. Taken together, our findings propose a model where the RFC1 deficiency drives methylation changes in specific genomic regions that are correlated with a dysregulation of pathways involved in early development such as Notch1 and BMP4 signaling. These changes affect the potency of the progenitors residing in the juncture of the neural plate and NPB territories, thus driving them to a primed state.

Repressive Interactions Between Transcription Factors Separate Different Embryonic Ectodermal Domains

The embryonic ectoderm is composed of four domains: neural plate, neural crest, pre-placodal region (PPR) and epidermis. Their formation is initiated during early gastrulation by dorsal-ventral and anterior-posterior gradients of signaling factors that first divide the embryonic ectoderm into neural and non-neural domains. Next, the neural crest and PPR domains arise, either via differential competence of the neural and non-neural ectoderm (binary competence model) or via interactions between the neural and non-neural ectoderm tissues to produce an intermediate neural border zone (NB) (border state model) that subsequently separates into neural crest and PPR. Many previous gain- and loss-of-function experiments demonstrate that numerous TFs are expressed in initially overlapping zones that gradually resolve into patterns that by late neurula stages are characteristic of each of the four domains. Several of these studies suggested that this is accomplished by a combination of repressive TF interactions and competence to respond to local signals. In this study, we ectopically expressed TFs that at neural plate stages are characteristic of one domain in a different domain to test whether they act cell autonomously as repressors. We found that almost all tested TFs caused reduced expression of the other TFs. At gastrulation these effects were strictly within the lineage-labeled cells, indicating that the effects were cell autonomous, i.e., due to TF interactions within individual cells. Analysis of previously published single cell RNAseq datasets showed that at the end of gastrulation, and continuing to neural tube closure stages, many ectodermal cells express TFs characteristic of more than one neural plate stage domain, indicating that different TFs have the opportunity to interact within the same cell. At neurula stages repression was observed both in the lineage-labeled cells and in adjacent cells not bearing detectable lineage label, suggesting that cell-to-cell signaling has begun to contribute to the separation of the domains. Together, these observations directly demonstrate previous suggestions in the literature that the segregation of embryonic ectodermal domains initially involves cell autonomous, repressive TF interactions within an individual cell followed by the subsequent advent of non-cell autonomous signaling to neighbors.

Kallmann syndrome in a patient with Weiss-Kruszka syndrome and a de novo deletion in 9q31.2

Patients with deletions on chromosome 9q31.2 may exhibit delayed puberty, craniofacial phenotype including cleft lip/palate, and olfactory bulb hypoplasia. We report a patient with congenital HH with anosmia (Kallmann syndrome, KS) and a de novo 2.38 Mb heterozygous deletion in 9q31.2. The deletion breakpoints (determined with whole-genome linked-read sequencing) were in the FKTN gene (9:108,331,353) and in a non-coding area (9:110,707,332) (hg19). The deletion encompassed six protein-coding genes (FKTN, ZNF462, TAL2, TMEM38B, RAD23B, and KLF4). ZNF462 haploinsufficiency was consistent with the patient's Weiss-Kruszka syndrome (craniofacial phenotype, developmental delay, and sensorineural hearing loss), but did not explain his KS. In further analyses, he did not carry rare sequence variants in 32 known KS genes in whole-exome sequencing and displayed no aberrant splicing of 15 KS genes that were expressed in peripheral blood leukocyte transcriptome. The deletion was 1.8 Mb upstream of a KS candidate gene locus (PALM2AKAP2) but did not suppress its expression. In conclusion, this is the first report of a patient with Weiss-Kruszka syndrome and KS. We suggest that patients carrying a microdeletion in 9q31.2 should be evaluated for the presence of KS and KS-related features.

Why Does the Face Predict the Brain? Neural Crest Induction, Craniofacial Morphogenesis, and Neural Circuit Development

Mesenchephalic and rhombencephalic neural crest cells generate the craniofacial skeleton, special sensory organs, and subsets of cranial sensory receptor neurons. They do so while preserving the anterior-posterior (A-P) identity of their neural tube origins. This organizational principle is paralleled by central nervous system circuits that receive and process information from facial structures whose A-P identity is in register with that in the brain. Prior to morphogenesis of the face and its circuits, however, neural crest cells act as "inductive ambassadors" from distinct regions of the neural tube to induce differentiation of target craniofacial domains and establish an initial interface between the brain and face. At every site of bilateral, non-axial secondary induction, neural crest constitutes all or some of the mesenchymal compartment for non-axial mesenchymal/epithelial (M/E) interactions. Thus, for epithelial domains in the craniofacial primordia, aortic arches, limbs, the spinal cord, and the forebrain (Fb), neural crest-derived mesenchymal cells establish local sources of inductive signaling molecules that drive morphogenesis and cellular differentiation. This common mechanism for building brains, faces, limbs, and hearts, A-P axis specified, neural crest-mediated M/E induction, coordinates differentiation of distal structures, peripheral neurons that provide their sensory or autonomic innervation in some cases, and central neural circuits that regulate their behavioral functions. The essential role of this neural crest-mediated mechanism identifies it as a prime target for pathogenesis in a broad range of neurodevelopmental disorders. Thus, the face and the brain "predict" one another, and this mutual developmental relationship provides a key target for disruption by developmental pathology.

Building the Border: Development of the Chordate Neural Plate Border Region and Its Derivatives

The paired cranial sensory organs and peripheral nervous system of vertebrates arise from a thin strip of cells immediately adjacent to the developing neural plate. The neural plate border region comprises progenitors for four key populations of cells: neural plate cells, neural crest cells, the cranial placodes, and epidermis. Putative homologues of these neural plate border derivatives can be found in protochordates such as amphioxus and tunicates. In this review, we summarize key signaling pathways and transcription factors that regulate the inductive and patterning events at the neural plate border region that give rise to the neural crest and placodal lineages. Gene regulatory networks driven by signals from WNT, fibroblast growth factor (FGF), and bone morphogenetic protein (BMP) signaling primarily dictate the formation of the crest and placodal lineages. We review these studies and discuss the potential of recent advances in spatio-temporal transcriptomic and epigenomic analyses that would allow a mechanistic understanding of how these signaling pathways and their downstream transcriptional cascades regulate the formation of the neural plate border region.

Insights Into the Early Gene Regulatory Network Controlling Neural Crest and Placode Fate Choices at the Neural Border

The neural crest (NC) cells and cranial placodes are two ectoderm-derived innovations in vertebrates that led to the acquisition of a complex head structure required for a predatory lifestyle. They both originate from the neural border (NB), a portion of the ectoderm located between the neural plate (NP), and the lateral non-neural ectoderm. The NC gives rise to a vast array of tissues and cell types such as peripheral neurons and glial cells, melanocytes, secretory cells, and cranial skeletal and connective cells. Together with cells derived from the cranial placodes, which contribute to sensory organs in the head, the NC also forms the cranial sensory ganglia. Multiple in vivo studies in different model systems have uncovered the signaling pathways and genetic factors that govern the positioning, development, and differentiation of these tissues. In this literature review, we give an overview of NC and placode development, focusing on the early gene regulatory network that controls the formation of the NB during early embryonic stages, and later dictates the choice between the NC and placode progenitor fates.

Accelerated differentiation of human pluripotent stem cells into neural lineages via an early intermediate ectoderm population

Differentiation of human pluripotent stem cells (hPSCs) into ectoderm provides neurons and glia useful for research, disease modeling, drug discovery, and potential cell therapies. In current protocols, hPSCs are traditionally differentiated into an obligate rostro-dorsal ectodermal fate expressing PAX6 after 6 to 12 days in vitro when protected from mesendoderm inducers. This rate-limiting step has performed a long-standing role in hindering the development of rapid differentiation protocols for ectoderm-derived cell types, as any protocol requires 6 to 10 days in vitro to simply initiate. Here, we report efficient differentiation of hPSCs into a naive early ectodermal intermediate within 24 hours using combined inhibition of bone morphogenic protein and fibroblast growth factor signaling. The induced population responds immediately to morphogen gradients to upregulate rostro-caudal neurodevelopmental landmark gene expression in a generally accelerated fashion. This method can serve as a new platform for the development of novel, rapid, and efficient protocols for the manufacture of hPSC-derived neural lineages.

Cell fate decisions during the development of the peripheral nervous system in the vertebrate head

Sensory placodes and neural crest cells are among the key cell populations that facilitated the emergence and diversification of vertebrates throughout evolution. Together, they generate the sensory nervous system in the head: both form the cranial sensory ganglia, while placodal cells make major contributions to the sense organs-the eye, ear and olfactory epithelium. Both are instrumental for integrating craniofacial organs and have been key to drive the concentration of sensory structures in the vertebrate head allowing the emergence of active and predatory life forms. Whereas the gene regulatory networks that control neural crest cell development have been studied extensively, the signals and downstream transcriptional events that regulate placode formation and diversity are only beginning to be uncovered. Both cell populations are derived from the embryonic ectoderm, which also generates the central nervous system and the epidermis, and recent evidence suggests that their initial specification involves a common molecular mechanism before definitive neural, neural crest and placodal lineages are established. In this review, we will first discuss the transcriptional networks that pattern the embryonic ectoderm and establish these three cell fates with emphasis on sensory placodes. Second, we will focus on how sensory placode precursors diversify using the specification of otic-epibranchial progenitors and their segregation as an example.

PRDM1 controls the sequential activation of neural, neural crest and sensory progenitor determinants

During early embryogenesis, the ectoderm is rapidly subdivided into neural, neural crest and sensory progenitors. How the onset of lineage determinants and the loss of pluripotency markers are temporally and spatially coordinated in vivo is still debated. Here, we identify a crucial role for the transcription factor PRDM1 in the orderly transition from epiblast to defined neural lineages in chick. PRDM1 is initially expressed broadly in the entire epiblast, but becomes gradually restricted as cell fates are specified. We find that PRDM1 is required for the loss of some pluripotency markers and the onset of neural, neural crest and sensory progenitor specifier genes. PRDM1 directly activates their expression by binding to their promoter regions and recruiting the histone demethylase Kdm4a to remove repressive histone marks. However, once neural lineage determinants become expressed, they in turn repress PRDM1, whereas prolonged PRDM1 expression inhibits neural, neural crest and sensory progenitor genes, suggesting that its downregulation is necessary for cells to maintain their identity. Therefore, PRDM1 plays multiple roles during ectodermal cell fate allocation.

Transplantation of Neural Tissue: Quail-Chick Chimeras

Tissue transplantation is an important approach in developmental neurobiology to determine cell fate, to uncover inductive interactions required for tissue specification and patterning as well as to establish tissue competence and commitment. Combined with state-of-the-art molecular approaches, transplantation assays have been instrumental for the discovery of gene regulatory networks controlling cell fate choices and how such networks change over time. Avian species are among the favorite model systems for these approaches because of their accessibility and relatively large size. Here we describe two culture techniques used to generate quail-chick chimeras at different embryonic stages and methods to distinguish graft and donor tissue.

Six1 and Irx1 have reciprocal interactions during cranial placode and otic vesicle formation

The specialized sensory organs of the vertebrate head are derived from thickened patches of cells in the ectoderm called cranial sensory placodes. The developmental program that generates these placodes and the genes that are expressed during the process have been studied extensively in a number of animals, yet very little is known about how these genes regulate one another. We previously found via a microarray screen that Six1, a known transcriptional regulator of cranial placode fate, up-regulates Irx1 in ectodermal explants. In this study, we investigated the transcriptional relationship between Six1 and Irx1 and found that they reciprocally regulate each other throughout cranial placode and otic vesicle formation. Although Irx1 expression precedes that of Six1 in the neural border zone, its continued and appropriately patterned expression in the pre-placodal region (PPR) and otic vesicle requires Six1. At early PPR stages, Six1 expands the Irx1 domain, but this activity subsides over time and changes to a predominantly repressive effect. Likewise, Irx1 initially expands Six1 expression in the PPR, but later represses it. We also found that Irx1 and Sox11, a known direct target of Six1, reciprocally affect each other. This work demonstrates that the interactions between Six1 and Irx1 are continuous during PPR and placode development and their transcriptional effects on one another change over developmental time.

A gene regulatory network underlying the formation of pre-placodal ectoderm in Xenopus laevis

BACKGROUND

The neural plate border ectoderm gives rise to key developmental structures during embryogenesis, including the neural crest and the preplacodal ectoderm. Many sensory organs and ganglia of vertebrates develop from cranial placodes, which themselves arise from preplacodal ectoderm, defined by expression of transcription factor Six1 and its coactivator Eya1. Here we elucidate the gene regulatory network underlying the specification of the preplacodal ectoderm in Xenopus, and the functional interactions among transcription factors that give rise to this structure.

RESULTS

To elucidate the gene regulatory network upstream of preplacodal ectoderm formation, we use gain- and loss-of-function studies to explore the role of early ectodermal transcription factors for establishing the preplacodal ectoderm and adjacent ectodermal territories, and the role of Six1 and Eya1 in feedback regulation of these transcription factors. Our findings suggest that transcription factors with expression restricted to ventral (non-neural) ectoderm (AP2, Msx1, FoxI1, Vent2, Dlx3, GATA2) and those restricted to dorsal (neural) ectoderm (Pax3, Hairy2b, Zic1) are required for specification of both preplacodal ectoderm and neural crest in a context-dependent fashion and are cross-regulated by Eya1 and Six1.

CONCLUSION

These findings allow us to elucidate a detailed gene regulatory network at the neural plate border upstream of preplacodal ectoderm formation based on functional interactions between ectodermal transcription factors. We propose a new model to explain the formation of immediately juxtaposed preplacodal ectoderm and neural crest territories at the neural plate border, uniting previous models.

Neural crest and the origin of species-specific pattern

For well over half of the 150 years since the discovery of the neural crest, the special ability of these cells to function as a source of species-specific pattern has been clearly recognized. Initially, this observation arose in association with chimeric transplant experiments among differentially pigmented amphibians, where the neural crest origin for melanocytes had been duly noted. Shortly thereafter, the role of cranial neural crest cells in transmitting species-specific information on size and shape to the pharyngeal arch skeleton as well as in regulating the timing of its differentiation became readily apparent. Since then, what has emerged is a deeper understanding of how the neural crest accomplishes such a presumably difficult mission, and this includes a more complete picture of the molecular and cellular programs whereby neural crest shapes the face of each species. This review covers studies on a broad range of vertebrates and describes neural-crest-mediated mechanisms that endow the craniofacial complex with species-specific pattern. A major focus is on experiments in quail and duck embryos that reveal a hierarchy of cell-autonomous and non-autonomous signaling interactions through which neural crest generates species-specific pattern in the craniofacial integument, skeleton, and musculature. By controlling size and shape throughout the development of these systems, the neural crest underlies the structural and functional integration of the craniofacial complex during evolution.

Neural induction by the node and placode induction by head mesoderm share an initial state resembling neural plate border and ES cells

Around the time of gastrulation in higher vertebrate embryos, inductive interactions direct cells to form central nervous system (neural plate) or sensory placodes. Grafts of different tissues into the periphery of a chicken embryo elicit different responses: Hensen's node induces a neural plate whereas the head mesoderm induces placodes. How different are these processes? Transcriptome analysis in time course reveals that both processes start by induction of a common set of genes, which later diverge. These genes are remarkably similar to those induced by an extraembryonic tissue, the hypoblast, and are normally expressed in the pregastrulation stage epiblast. Explants of this epiblast grown in the absence of further signals develop as neural plate border derivatives and eventually express lens markers. We designate this state as "preborder"; its transcriptome resembles embryonic stem cells. Finally, using sequential transplantation experiments, we show that the node, head mesoderm, and hypoblast are interchangeable to begin any of these inductions while the final outcome depends on the tissue emitting the later signals.

Spatiotemporal expression pattern of Connexin 43 during early chick embryogenesis

During embryogenesis, a single cell develops into new tissues and organs that are made up of a number of different cell types. The assembly of the trigeminal ganglion (cranial nerve V), an important component of the peripheral nervous system, typifies this process. The trigeminal ganglia perform key sensory functions, including sensing pain and touch in the face, and arise from cells of two different progenitor populations, the neural crest and the cranial placodes. One question that remains poorly understood is how these two populations of cells interact with each other during development to form a functional ganglion. Gap junctions are intercellular channels that allow for the passage of small solutes between connected cells and could serve as one potential mechanism by which neural crest and placode cells communicate to create the trigeminal ganglia. To this end, we have generated a comprehensive spatiotemporal expression profile for the gap junction protein Connexin 43, a highly expressed member of the Connexin protein family during development. Our results reveal that Connexin 43 is expressed in the neural folds during neural fold fusion and in premigratory neural crest cells prior to the epithelial-to-mesenchymal transition (EMT), during EMT, and in migratory neural crest cells. During trigeminal gangliogenesis, Connexin 43 is expressed in cranial neural crest cells and the mesenchyme but is strikingly absent in the placode-derived neurons. These data underscore the complexity of bringing two distinct cell populations together to form a new tissue during development and suggest that Connexin 43 may play a key role within neural crest cells during EMT, migration, and trigeminal gangliogenesis.

Retinoic acid is required and Fgf, Wnt, and Bmp signaling inhibit posterior lateral line placode induction in zebrafish

The lateral line system is a mechanosensory systems present in aquatic animals. The anterior and posterior lateral lines develop from anterior and posterior lateral line placodes (aLLp and pLLp), respectively. Although signaling molecules required for the induction of other cranial placodes have been well studied, the molecular mechanisms underlying formation of the lateral line placodes are unknown. In this study we tested the requirement of multiple signaling pathways, such as Wnt, Bmp Fgf, and Retinoic Acid for aLLp and pLLp induction. We determined that aLLp specification requires Fgf signaling, whilst pLLp specification requires retinoic acid which inhibits Fgf signaling. pLLp induction is also independent of Wnt and Bmp activities, even though these pathways limit the boundaries of the pLLp. This is the first report that the aLLp and pLLp depend on different inductive mechanisms and that pLLp induction requires the inhibition of Fgf, Wnt and Bmp signaling.