Cdx1 refines positional identity of the vertebrate hindbrain by directly repressing Mafb expression

Loss of the epithelial marker CDX1 predicts poor prognosis in early-stage CRC patients

BACKGROUND

We have previously shown that non-curative chemotherapy imposes fetal conversion and high metastatic capacity to cancer cells. From the set of genes differentially expressed in Chemotherapy Resistant Cells, we obtained a characteristic fetal intestinal cell signature that is present in a group of untreated tumors and is sufficient to predict patient prognosis. A feature of this fetal signature is the loss of CDX1.

METHODS

We have analyzed transcriptomic data in public datasets and performed immunohistochemistry analysis of paraffin embedded tumor samples from two cohorts of colorectal cancer patients.

RESULTS

We demonstrated that low levels of CDX1 are sufficient to identify patients with poorest outcome at the early tumor stages II and III. Presence tumor areas that are negative for CDX1 staining in stage I cancers is associated with tumor relapse.

CONCLUSIONS

Our results reveal the actual possibility of incorporating CDX1 immunostaining as a valuable biomarker for CRC patients.

The control of transitions along the main body axis

Although vertebrates display a large variety of forms and sizes, the mechanisms controlling the layout of the basic body plan are substantially conserved throughout the clade. Following gastrulation, head, trunk, and tail are sequentially generated through the continuous addition of tissue at the caudal embryonic end. Development of each of these major embryonic regions is regulated by a distinct genetic network. The transitions from head-to-trunk and from trunk-to-tail development thus involve major changes in regulatory mechanisms, requiring proper coordination to guarantee smooth progression of embryonic development. In this review, we will discuss the key cellular and embryological events associated with those transitions giving particular attention to their regulation, aiming to provide a cohesive outlook of this important component of vertebrate development.

Enhanced production of mesencephalic dopaminergic neurons from lineage-restricted human undifferentiated stem cells

Current differentiation protocols for generating mesencephalic dopaminergic (mesDA) neurons from human pluripotent stem cells result in grafts containing only a small proportion of mesDA neurons when transplanted in vivo. In this study, we develop lineage-restricted undifferentiated stem cells (LR-USCs) from pluripotent stem cells, which enhances their potential for differentiating into caudal midbrain floor plate progenitors and mesDA neurons. Using a ventral midbrain protocol, 69% of LR-USCs become bona fide caudal midbrain floor plate progenitors, compared to only 25% of human embryonic stem cells (hESCs). Importantly, LR-USCs generate significantly more mesDA neurons under midbrain and hindbrain conditions in vitro and in vivo. We demonstrate that midbrain-patterned LR-USC progenitors transplanted into 6-hydroxydopamine-lesioned rats restore function in a clinically relevant non-pharmacological behavioral test, whereas midbrain-patterned hESC-derived progenitors do not. This strategy demonstrates how lineage restriction can prevent the development of undesirable lineages and enhance the conditions necessary for mesDA neuron generation.

The cardiac neural crest gene MafB ectopically directs CXCR4 expression in the trunk neural crest

The cardiac neural crest is a subpopulation of cells arising from the caudal hindbrain. The delaminated cardiac neural crest cells migrate to the heart using the CXCR/SDF1 chemokine signaling system. These cells contribute to the formation of the cardiovascular system, including the septation of the outflow tract, which is unique to these cells. Here, we investigated the effect of ectopic expression of the cardiac neural crest gene MafB on trunk neural crest cells. First, we found that MafB has the potential to activate its own cis-regulatory element in enteric and trunk neural crest cells but not in cranial neural crest cells. Forced expression of two cardiac neural crest genes, Ets1 and Sox8, together with or without MafB, induced ectopic Sox10E2 enhancer activity in the trunk region. Finally, we uncovered that the expression of MafB, Ets1 and Sox8 can induce ectopic CXCR4 expression in the trunk neural crest cells, resulting in acquisition of responsiveness to the SDF1 signal. These results demonstrate that MafB, Ets1 and Sox8 are critical components for generation of the identity of the cardiac neural crest, especially the cell migration property.

Systematic reconstruction of cellular trajectories across mouse embryogenesis

Mammalian embryogenesis is characterized by rapid cellular proliferation and diversification. Within a few weeks, a single-cell zygote gives rise to millions of cells expressing a panoply of molecular programs. Although intensively studied, a comprehensive delineation of the major cellular trajectories that comprise mammalian development in vivo remains elusive. Here, we set out to integrate several single-cell RNA-sequencing (scRNA-seq) datasets that collectively span mouse gastrulation and organogenesis, supplemented with new profiling of ~150,000 nuclei from approximately embryonic day 8.5 (E8.5) embryos staged in one-somite increments. Overall, we define cell states at each of 19 successive stages spanning E3.5 to E13.5 and heuristically connect them to their pseudoancestors and pseudodescendants. Although constructed through automated procedures, the resulting directed acyclic graph (TOME (trajectories of mammalian embryogenesis)) is largely consistent with our contemporary understanding of mammalian development. We leverage TOME to systematically nominate transcription factors (TFs) as candidate regulators of each cell type's specification, as well as 'cell-type homologs' across vertebrate evolution.

Regulation of axial elongation by Cdx

The primordia of the post-otic mouse embryo forms largely from a bipotential cell population containing neuromesodermal progenitors (NMP) which reside in the tail bud and contribute to the elaboration of the major body axis after gastrulation. The mechanisms by which the NMP population is both maintained and subsequently directed down mesodermal and neural lineages is incompletely understood. The vertebrate transcription factor Cdx2, is essential for axial elongation and has been implicated in maintaining the NMP niche and in specification of NMP derivatives. To better understand the role of the Cdx family in axial elongation, we employed a conditional mutant allele which evokes total loss of Cdx function, and enriched for tail bud progenitors through the use of a Pax2-GFP transgenic reporter. Using this approach, we identified 349 Cdx-dependent genes by RNA sequencing (RNA-seq). From these, Gene Ontology and chromatin immunoprecipitation analysis further revealed a number of putative direct Cdx candidate target genes implicated in axial elongation, including Sp8, Isl1, Evx1, Zic3 and Nr2f1. Additional analysis of available single-cell RNA-seq data from mouse tail buds revealed the co-expression of Sp8, Isl1, Evx1 and Zic3 with Cdx2 in putative NMP cells, while Nr2f1 was excluded from this population. These findings identify a number of novel Cdx targets and provide further insight into the critical roles for Cdx in elaborating the post-otic embryo.

Identification of regulatory elements for MafB expression in the cardiac neural crest

Cardiac neural crest cells arise in the caudal hindbrain and then migrate to the heart through the pharyngeal arches. These cells contribute to the formation of the heart, including the outflow tract, and are unique to this neural crest population. MafB is a transcription factor expressed specifically in early migrating cardiac neural crest cells as well as in rhombomeres (r) 5 and 6. Here, we identified the regulatory region in the chicken genome controlling the expression of endogenous MafB transcripts and used these essential elements to express MafB in the cardiac neural crest in reporter assays. A reporter driven by this regulatory region was employed to trace the migration of these cells into the pharyngeal arches. This regulatory region demonstrated transcriptional activity in the cardiac neural crest but not in other neural crest cell subpopulations, such as the cranial and trunk cells. This study provides insights into the gene regulatory mechanisms that specify cardiac neural crest cells among neural crest cell populations.

Retinoids and developmental neurotoxicity: Utilizing toxicogenomics to enhance adverse outcome pathways and testing strategies

The use of genomic approaches in toxicological studies has greatly increased our ability to define the molecular profiles of environmental chemicals associated with developmental neurotoxicity (DNT). Integration of these approaches with adverse outcome pathways (AOPs), a framework that translates environmental exposures to adverse developmental phenotypes, can potentially inform DNT testing strategies. Here, using retinoic acid (RA) as a case example, we demonstrate that the integration of toxicogenomic profiles into the AOP framework can be used to establish a paradigm for chemical testing. RA is a critical regulatory signaling molecule involved in multiple aspects of mammalian central nervous system (CNS) development, including hindbrain formation/patterning and neuronal differentiation, and imbalances in RA signaling pathways are linked with DNT. While the mechanisms remain unresolved, environmental chemicals can cause DNT by disrupting the RA signaling pathway. First, we reviewed literature evidence of RA and other retinoid exposures and DNT to define a provisional AOP related to imbalances in RA embryonic bioavailability and hindbrain development. Next, by integrating toxicogenomic datasets, we defined a relevant transcriptomic signature associated with RA-induced developmental neurotoxicity (RA-DNT) in human and rodent models that was tested against zebrafish model data, demonstrating potential for integration into an AOP framework. Finally, we demonstrated how these approaches may be systematically utilized to identify chemical hazards by testing the RA-DNT signature against azoles, a proposed class of compounds that alters RA-signaling. The provisional AOP from this study can be expanded in the future to better define DNT biomarkers relevant to RA signaling and toxicity.

Retinoic acid signaling pathways

Retinoic acid (RA), a metabolite of retinol (vitamin A), functions as a ligand for nuclear RA receptors (RARs) that regulate development of chordate animals. RA-RARs can activate or repress transcription of key developmental genes. Genetic studies in mouse and zebrafish embryos that are deficient in RA-generating enzymes or RARs have been instrumental in identifying RA functions, revealing that RA signaling regulates development of many organs and tissues, including the body axis, spinal cord, forelimbs, heart, eye and reproductive tract. An understanding of the normal functions of RA signaling during development will guide efforts for use of RA as a therapeutic agent to improve human health. Here, we provide an overview of RA signaling and highlight its key functions during development.

Maternal Glucocorticoid Elevation and Associated Fetal Thymocyte Apoptosis are Involved in Immune Disorders of Prenatal Caffeine Exposed Offspring Mice

Our previous study showed that prenatal caffeine exposure (PCE) could induce intrauterine growth retardation (IUGR) and glucocorticoid elevation in the fetus. Researchers suggested that IUGR is a risk factor for T helper cell (Th)1/Th2 deviation. However, whether PCE can induce these immune disorders and the underlying mechanisms of that induction remain unknown. This study aimed to observe the effects of PCE on the Th1/Th2 balance in offspring and further explore the developmental origin mechanisms from the perspective of glucocorticoid overexposure-induced thymocyte apoptosis. An IUGR model was established by caffeine administration from gestational day (GD) 9 to GD 18, and the offspring were immunized on postnatal day (PND) 42. The results show that maternal glucocorticoid overexposure increased fetal thymocyte apoptosis by activating both the Fas-mediated and the Bim-regulated apoptotic pathways. After birth, accelerated thymocyte apoptosis and Th1 suppression were also found in the PCE offspring at PND 14 and PND 49. Moreover, the PCE offspring showed immune disorders after immunization, manifesting as increased IgG1/IgG2a ratio and IL-4 production in the serum. In conclusion, PCE could induce fetal overexposure to maternal glucocorticoids and increase thymocyte apoptosis, which could persist into postnatal life and be implicated in Th1 inhibition and further immune disorders.

Segmental arithmetic: summing up the Hox gene regulatory network for hindbrain development in chordates

Organization and development of the early vertebrate hindbrain are controlled by a cascade of regulatory interactions that govern the process of segmentation and patterning along the anterior-posterior axis via Hox genes. These interactions can be assembled into a gene regulatory network that provides a framework to interpret experimental data, generate hypotheses, and identify gaps in our understanding of the progressive process of hindbrain segmentation. The network can be broadly separated into a series of interconnected programs that govern early signaling, segmental subdivision, secondary signaling, segmentation, and ultimately specification of segmental identity. Hox genes play crucial roles in multiple programs within this network. Furthermore, the network reveals properties and principles that are likely to be general to other complex developmental systems. Data from vertebrate and invertebrate chordate models are shedding light on the origin and diversification of the network. Comprehensive cis-regulatory analyses of vertebrate Hox gene regulation have enabled powerful cross-species gene regulatory comparisons. Such an approach in the sea lamprey has revealed that the network mediating segmental Hox expression was present in ancestral vertebrates and has been maintained across diverse vertebrate lineages. Invertebrate chordates lack hindbrain segmentation but exhibit conservation of some aspects of the network, such as a role for retinoic acid in establishing nested Hox expression domains. These comparisons lead to a model in which early vertebrates underwent an elaboration of the network between anterior-posterior patterning and Hox gene expression, leading to the gene-regulatory programs for segmental subdivision and rhombomeric segmentation. WIREs Dev Biol 2017, 6:e286. doi: 10.1002/wdev.286 For further resources related to this article, please visit the WIREs website.

Genetic analysis of the ATG16L1 gene promoter in sporadic Parkinson's disease

Parkinson's disease (PD) is a common and progressive neurodegenerative disease in which the majority of cases arise sporadically. Sporadic PD is caused by the interactions of genetic and environmental factors. To date, genetic causes for sporadic PD remain largely unknown. Autophagy, a highly conserved cellular process, has been implicated in PD pathogenesis. We speculated that genetic variants in autophagy-related genes (ATG) that regulate gene expression may contribute to PD development. In our previous studies, we have identified several functional DNA sequence variants (DSVs) in the ATG5, ATG7 and LC3 genes in sporadic PD patients. In this study, we further genetically and functionally analyzed the promoter of the ATG16L1 gene, a critical gene for autophagosome formation, in groups of sporadic PD patients and ethnic-matched healthy controls. One novel heterozygous DSV, 233251432C>T, was found in one PD patient. Functionally, this DSV did not affect the transcriptional activity of the ATG16L1 gene promoter in human dopaminergic SH-SY5Y cells. Two heterozygous DSVs including one SNP, 233251286G>A (rs539735288) and 233251582C>T, were found only in controls. In addition, five other SNPs were found in both PD patients and controls. Taken together, the data suggested that genetic variants within the ATG16L1 gene promoter were not a risk factor for sporadic PD development.

Targeted resequencing of regulatory regions at schizophrenia risk loci: Role of rare functional variants at chromatin repressive states

There is mounting evidence that regulatory variation plays an important role in genetic risk for schizophrenia. Here, we specifically search for regulatory variants at risk by sequencing promoter regions of twenty-three genes implied in schizophrenia by copy number variant or genome-wide association studies. After strict quality control, a total of 55,206bp per sample were analyzed in 526 schizophrenia cases and 516 controls from Galicia, NW Spain, using the Applied Biosystems SOLiD System. Variants were filtered based on frequency from public databases, chromatin states from the RoadMap Epigenomics Consortium at tissues relevant for schizophrenia, such as fetal brain, mid-frontal lobe, and angular gyrus, and prediction of functionality from RegulomeDB. The proportion of rare variants at polycomb repressive chromatin state at relevant tissues was higher in cases than in controls. The proportion of rare variants with predicted regulatory role was significantly higher in cases than in controls (P=0.0028, OR=1.93, 95% C.I.=1.23-3.04). Combination of information from both sources led to the identification of an excess of carriers of rare variants with predicted regulatory role located at polycomb repressive chromatin state at relevant tissues in cases versus controls (P=0.0016, OR=19.34, 95% C.I.=2.45-2495.26). The variants are located at two genes affected by the 17q12 copy number variant, LHX1 and HNF1B. These data strongly suggest that a specific epigenetic mechanism, chromatin remodeling by histone modification during early development, may be impaired in a subset of schizophrenia patients, in agreement with previous data.

The vertebrate Hox gene regulatory network for hindbrain segmentation: Evolution and diversification: Coupling of a Hox gene regulatory network to hindbrain segmentation is an ancient trait originating at the base of vertebrates

Hindbrain development is orchestrated by a vertebrate gene regulatory network that generates segmental patterning along the anterior-posterior axis via Hox genes. Here, we review analyses of vertebrate and invertebrate chordate models that inform upon the evolutionary origin and diversification of this network. Evidence from the sea lamprey reveals that the hindbrain regulatory network generates rhombomeric compartments with segmental Hox expression and an underlying Hox code. We infer that this basal feature was present in ancestral vertebrates and, as an evolutionarily constrained developmental state, is fundamentally important for patterning of the vertebrate hindbrain across diverse lineages. Despite the common ground plan, vertebrates exhibit neuroanatomical diversity in lineage-specific patterns, with different vertebrates revealing variations of Hox expression in the hindbrain that could underlie this diversification. Invertebrate chordates lack hindbrain segmentation but exhibit some conserved aspects of this network, with retinoic acid signaling playing a role in establishing nested domains of Hox expression.

A direct role for murine Cdx proteins in the trunk neural crest gene regulatory network

Numerous studies in chordates and arthropods currently indicate that Cdx proteins have a major ancestral role in the organization of post-head tissues. In urochordate embryos, Cdx loss-of-function has been shown to impair axial elongation, neural tube (NT) closure and pigment cell development. Intriguingly, in contrast to axial elongation and NT closure, a Cdx role in neural crest (NC)-derived melanocyte/pigment cell development has not been reported in any other chordate species. To address this, we generated a new conditional pan-Cdx functional knockdown mouse model that circumvents Cdx functional redundancy as well as the early embryonic lethality of Cdx mutants. Through directed inhibition in the neuroectoderm, we provide in vivo evidence that murine Cdx proteins impact melanocyte and enteric nervous system development by, at least in part, directly controlling the expression of the key early regulators of NC ontogenesis Pax3,Msx1 and Foxd3 Our work thus reveals a novel role for Cdx proteins at the top of the trunk NC gene regulatory network in the mouse, which appears to have been inherited from their ancestral ortholog.

CDX4 and retinoic acid interact to position the hindbrain-spinal cord transition

The sub-division of the posterior-most territory of the neural plate results in the formation of two distinct neural structures, the hindbrain and the spinal cord. Although many of the molecular signals regulating the development of these individual structures have been elucidated, the mechanisms involved in delineating the boundary between the hindbrain and spinal cord remain elusive. Two molecules, retinoic acid (RA) and the Cdx4 transcription factor have been previously implicated as important regulators of hindbrain and spinal cord development, respectively. Here, we provide evidence that suggests multiple regulatory interactions occur between RA signaling and the Cdx4 transcription factor to establish the anterior-posterior (AP) position of the transition between the hindbrain and spinal cord. Using chemical inhibitors to alter RA concentrations and morpholinos to knock-down Cdx4 function in zebrafish, we show that Cdx4 acts to prevent RA degradation in the presumptive spinal cord domain by suppressing expression of the RA degradation enzyme, Cyp26a1. In the hindbrain, RA signaling modulates its own concentration by activating the expression of cyp26a1 and inhibiting the expansion of cdx4. Therefore, interactions between Cyp26a1 and Cdx4 modulate RA levels along the AP axis to segregate the posterior neural plate into the hindbrain and spinal cord territories.

The orphan GPCR, Gpr161, regulates the retinoic acid and canonical Wnt pathways during neurulation

The vacuolated lens (vl) mouse mutation arose on the C3H/HeSnJ background and results in lethality, neural tube defects (NTDs) and cataracts. The vl phenotypes are due to a deletion/frameshift mutation in the orphan GPCR, Gpr161. A recent study using a null allele demonstrated that Gpr161 functions in primary cilia and represses the Shh pathway. We show the hypomorphic Gpr161(vl) allele does not severely affect the Shh pathway. To identify additional pathways regulated by Gpr161 during neurulation, we took advantage of naturally occurring genetic variation in the mouse. Previously Gpr161(vl-C3H) was crossed to different inbred backgrounds including MOLF/EiJ and the Gpr161(vl) mutant phenotypes were rescued. Five modifiers were mapped (Modvl: Modifier of vl) including Modvl5(MOLF). In this study we demonstrate the Modvl5(MOLF) congenic rescues the Gpr161(vl)-associated lethality and NTDs but not cataracts. Bioinformatics determined the transcription factor, Cdx1, is the only annotated gene within the Modvl5 95% CI co-expressed with Gpr161 during neurulation and not expressed in the eye. Using Cdx1 as an entry point, we identified the retinoid acid (RA) and canonical Wnt pathways as downstream targets of Gpr161. QRT-PCR, ISH and IHC determined that expression of RA and Wnt genes are down-regulated in Gpr161(vl/vl) but rescued by the Modvl5(MOLF) congenic during neurulation. Intraperitoneal RA injection restores expression of canonical Wnt markers and rescues Gpr161(vl/vl) NTDs. These results establish the RA and canonical Wnt as pathways downstream of Gpr161 during neurulation, and suggest that Modvl5(MOLF) bypasses the Gpr161(vl) mutation by restoring the activity of these pathways.

The gene regulatory networks underlying formation of the auditory hindbrain

Development and evolution of auditory hindbrain nuclei are two major unsolved issues in hearing research. Recent characterization of transgenic mice identified the rhombomeric origins of mammalian auditory nuclei and unraveled genes involved in their formation. Here, we provide an overview on these data by assembling them into rhombomere-specific gene regulatory networks (GRNs), as they underlie developmental and evolutionary processes. To explore evolutionary mechanisms, we compare the GRNs operating in the mammalian auditory hindbrain with data available from the inner ear and other vertebrate groups. Finally, we propose that the availability of genomic sequences from all major vertebrate taxa and novel genetic techniques for non-model organisms provide an unprecedented opportunity to investigate development and evolution of the auditory hindbrain by comparative molecular approaches. The dissection of the molecular mechanisms leading to auditory structures will also provide an important framework for auditory processing disorders, a clinical problem difficult to tackle so far. These data will, therefore, foster basic and clinical hearing research alike.

A novel function for Egr4 in posterior hindbrain development

Segmentation of the vertebrate hindbrain is an evolutionarily conserved process. Here, we identify the transcription factor early growth response 4 (egr4) as a novel regulator of posterior hindbrain development in Xenopus. egr4 is specifically and transiently expressed in rhombomeres 5 and 6 (r5/r6), and Egr4 knockdown causes a loss of mafb/kreisler and krox20/egr2 expression in r5/r6 and r5, respectively. This phenotype can be fully rescued by injection of frog or mouse Egr4 mRNA. Moreover Egr4-depleted embryos exhibit a specific loss of the neural crest stream adjacent to r5, and have inner ear defects. While the homeodomain protein vHnf1/Hnf1b directly activates Mafb and Krox20 expression in the mouse hindbrain to specify r5, we show that in Xenopus this process is indirect through the activation of Egr4. We provide evidence that rearrangements in the regulatory sequences around egr4 and mafb genes may account for this difference.

Retinoic acid regulates size, pattern and alignment of tissues at the head-trunk transition

At the head-trunk transition, hindbrain and spinal cord alignment to occipital and vertebral bones is crucial for coherent neural and skeletal system organization. Changes in neural or mesodermal tissue configuration arising from defects in the specification, patterning or relative axial placement of territories can severely compromise their integration and function. Here, we show that coordination of neural and mesodermal tissue at the zebrafish head-trunk transition crucially depends on two novel activities of the signaling factor retinoic acid (RA): one specifying the size and the other specifying the axial position relative to mesodermal structures of the hindbrain territory. These activities are each independent but coordinated with the well-established function of RA in hindbrain patterning. Using neural and mesodermal landmarks we demonstrate that the functions of RA in aligning neural and mesodermal tissues temporally precede the specification of hindbrain and spinal cord territories and the activation of hox transcription. Using cell transplantation assays we show that RA activity in the neuroepithelium regulates hindbrain patterning directly and territory size specification indirectly. This indirect function is partially dependent on Wnts but independent of FGFs. Importantly, RA specifies and patterns the hindbrain territory by antagonizing the activity of the spinal cord specification gene cdx4; loss of Cdx4 rescues the defects associated with the loss of RA, including the reduction in hindbrain size and the loss of posterior rhombomeres. We propose that at the head-trunk transition, RA coordinates specification, patterning and alignment of neural and mesodermal tissues that are essential for the organization and function of the neural and skeletal systems.

Combinatorial analysis of developmental cues efficiently converts human pluripotent stem cells into multiple neuronal subtypes

Specification of cell identity during development depends on exposure of cells to sequences of extrinsic cues delivered at precise times and concentrations. Identification of combinations of patterning molecules that control cell fate is essential for the effective use of human pluripotent stem cells (hPSCs) for basic and translational studies. Here we describe a scalable, automated approach to systematically test the combinatorial actions of small molecules for the targeted differentiation of hPSCs. Applied to the generation of neuronal subtypes, this analysis revealed an unappreciated role for canonical Wnt signaling in specifying motor neuron diversity from hPSCs and allowed us to define rapid (14 days), efficient procedures to generate spinal and cranial motor neurons as well as spinal interneurons and sensory neurons. Our systematic approach to improving hPSC-targeted differentiation should facilitate disease modeling studies and drug screening assays.

An adverse outcome pathway framework for neural tube and axial defects mediated by modulation of retinoic acid homeostasis

Developmental toxicity can be caused through a multitude of mechanisms and can therefore not be captured through a single simple mechanistic paradigm. However, it may be possible to define a selected group of overarching mechanisms that might allow detection of the vast majority of developmental toxicants. Against this background, we have explored the usefulness of retinoic acid mediated regulation of neural tube and axial patterning as a general mechanism that, when perturbed, may result in manifestations of developmental toxicity that may cover a large part of malformations known to occur in experimental animals and in man. Through a literature survey, we have identified key genes in the regulation of retinoic acid homeostasis, as well as marker genes of neural tube and axial patterning, that may be used to detect developmental toxicants in in vitro systems. A retinoic acid-neural tube/axial patterning adverse outcome pathway (RA-NTA AOP) framework was designed. The framework was tested against existing data of flusilazole exposure in the rat whole embryo culture, the zebrafish embryotoxicity test, and the embryonic stem cell test. Flusilazole is known to interact with retinoic acid homeostasis, and induced common and unique NTA marker gene changes in the three test systems. Flusilazole-induced changes were similar in directionality to gene expression responses after retinoic acid exposure. It is suggested that the RA-NTA framework may provide a general tool to define mechanistic pathways and biomarkers of developmental toxicity that may be used in alternative in vitro assays for the detection of embryotoxic compounds.

A novel role of CDX1 in embryonic epicardial development

The molecular mechanism that regulates epicardial development has yet to be understood. In this study, we explored the function of CDX1, a Caudal-related family member, in epicardial epithelial-to-mesenchymal transition (EMT) and in the migration and the differentiation of epicardium-derived progenitors into vascular smooth muscle cells. We detected a transient expression of CDX1 in murine embryonic hearts at 11.5 days post coitum (dpc). Using a doxycycline-inducible CDX1 mouse model, primary epicardium, and ex vivo heart culture, we further demonstrated that ectopic expression of CDX1 promoted epicardial EMT. In addition, a low-dose CDX1 induction led to enhanced migration and differentiation of epicardium-derived cells into α-SMA+ vascular smooth muscles. In contrast, either continued high-level induction of CDX1 or CDX1 deficiency attenuated the ability of epicardium-derived cells to migrate and to mature into smooth muscles induced by TGF-β1. Further RNA-seq analyses showed that CDX1 induction altered the transcript levels of genes involved in neuronal development, angiogenesis, and cell adhesions required for EMT. Our data have revealed a previously undefined role of CDX1 during epicardial development, and suggest that transient expression of CDX1 promotes epicardial EMT, whereas subsequent down-regulation of CDX1 after 11.5 dpc in mice is necessary for further subepicardial invasion of EPDCs and contribution to coronary vascular endothelium or smooth muscle cells.

The Prdm13 histone methyltransferase encoding gene is a Ptf1a-Rbpj downstream target that suppresses glutamatergic and promotes GABAergic neuronal fate in the dorsal neural tube

The basic helix-loop-helix (bHLH) transcriptional activator Ptf1a determines inhibitory GABAergic over excitatory glutamatergic neuronal cell fate in progenitors of the vertebrate dorsal spinal cord, cerebellum and retina. In an in situ hybridization expression survey of PR domain containing genes encoding putative chromatin-remodeling zinc finger transcription factors in Xenopus embryos, we identified Prdm13 as a histone methyltransferase belonging to the Ptf1a synexpression group. Gain and loss of Ptf1a function analyses in both frog and mice indicates that Prdm13 is positively regulated by Ptf1a and likely constitutes a direct transcriptional target. We also showed that this regulation requires the formation of the Ptf1a-Rbp-j complex. Prdm13 knockdown in Xenopus embryos and in Ptf1a overexpressing ectodermal explants lead to an upregulation of Tlx3/Hox11L2, which specifies a glutamatergic lineage and a reduction of the GABAergic neuronal marker Pax2. It also leads to an upregulation of Prdm13 transcription, suggesting an autonegative regulation. Conversely, in animal caps, Prdm13 blocks the ability of the bHLH factor Neurog2 to activate Tlx3. Additional gain of function experiments in the chick neural tube confirm that Prdm13 suppresses Tlx3(+)/glutamatergic and induces Pax2(+)/GABAergic neuronal fate. Thus, Prdm13 is a novel crucial component of the Ptf1a regulatory pathway that, by modulating the transcriptional activity of bHLH factors such as Neurog2, controls the balance between GABAergic and glutamatergic neuronal fate in the dorsal and caudal part of the vertebrate neural tube.

Characterization of a novel transgenic mouse line expressing Cre recombinase under the control of the Cdx2 neural specific enhancer

Several genetically modified mouse models have been generated in order to drive expression of the Cre recombinase in the neuroectoderm. However, none of them specifically targets the posterior neural plate during neurulation. To fill this gap, we have generated a new transgenic mouse line in which Cre expression is controlled by a neural specific enhancer (NSE) from the Caudal-related homeobox 2 (Cdx2) locus. Analyses of Cre activity via breeding with R26R-YFP reporter mice have indicated that the Cdx2NSE-Cre mouse line allows for recombination of LoxP sites in most cells of the posterior neural plate as soon as from the head fold stage. Detailed examination of double-transgenic embryos has revealed that this novel Cre-driver line allows targeting the entire posterior neural tube with an anterior limit in the caudal hindbrain. Of note, the Cdx2NSE regulatory sequences direct Cre expression along the whole dorso-ventral axis (including pre-migratory neural crest cells) and, accordingly, YFP fluorescence has been also observed in multiple non-cranial neural crest derivatives of double-transgenic embryos. Therefore, we believe that the Cdx2NSE-Cre mouse line represents an important novel genetic tool for the study of early events occurring in the caudal neuroectoderm during the formation of both the central and the peripheral nervous systems.

The evolution of human genetic studies of cleft lip and cleft palate

Orofacial clefts (OFCs)--primarily cleft lip and cleft palate--are among the most common birth defects in all populations worldwide, and have notable population, ethnicity, and gender differences in birth prevalence. Interest in these birth defects goes back centuries, as does formal scientific interest; scientists often used OFCs as examples or evidence during paradigm shifts in human genetics, and have also used virtually every new method of human genetic analysis to deepen our understanding of OFC. This review traces the evolution of human genetic investigations of OFC, highlights the specific insights gained about OFC through the years, and culminates in a review of recent key OFC genetic findings resulting from the powerful tools of the genomics era. Notably, OFC represents a major success for genome-wide approaches, and the field is poised for further breakthroughs in the near future.

Caudal-related homeobox (Cdx) protein-dependent integration of canonical Wnt signaling on paired-box 3 (Pax3) neural crest enhancer

One of the earliest events in neural crest development takes place at the neural plate border and consists in the induction of Pax3 expression by posteriorizing Wnt·β-catenin signaling. The molecular mechanism of this regulation is not well understood, but several observations suggest a role for posteriorizing Cdx transcription factors (Cdx1/2/4) in this process. Cdx genes are known as integrators of posteriorizing signals from Wnt, retinoic acid, and FGF pathways. In this work, we report that Wnt-mediated regulation of murine Pax3 expression is indirect and involves Cdx proteins as intermediates. We show that Pax3 transcripts co-localize with Cdx proteins in the posterior neurectoderm and that neural Pax3 expression is reduced in Cdx1-null embryos. Using Wnt3a-treated P19 cells and neural crest-derived Neuro2a cells, we demonstrate that Pax3 expression is induced by the Wnt-Cdx pathway. Co-transfection analyses, electrophoretic mobility shift assays, chromatin immunoprecipitation, and transgenic studies further indicate that Cdx proteins operate via direct binding to an evolutionarily conserved neural crest enhancer of the Pax3 proximal promoter. Taken together, these results suggest a novel neural function for Cdx proteins within the gene regulatory network controlling neural crest development.

TF-Cluster: a pipeline for identifying functionally coordinated transcription factors via network decomposition of the shared coexpression connectivity matrix (SCCM)

BACKGROUND

Identifying the key transcription factors (TFs) controlling a biological process is the first step toward a better understanding of underpinning regulatory mechanisms. However, due to the involvement of a large number of genes and complex interactions in gene regulatory networks, identifying TFs involved in a biological process remains particularly difficult. The challenges include: (1) Most eukaryotic genomes encode thousands of TFs, which are organized in gene families of various sizes and in many cases with poor sequence conservation, making it difficult to recognize TFs for a biological process; (2) Transcription usually involves several hundred genes that generate a combination of intrinsic noise from upstream signaling networks and lead to fluctuations in transcription; (3) A TF can function in different cell types or developmental stages. Currently, the methods available for identifying TFs involved in biological processes are still very scarce, and the development of novel, more powerful methods is desperately needed.

RESULTS

We developed a computational pipeline called TF-Cluster for identifying functionally coordinated TFs in two steps: (1) Construction of a shared coexpression connectivity matrix (SCCM), in which each entry represents the number of shared coexpressed genes between two TFs. This sparse and symmetric matrix embodies a new concept of coexpression networks in which genes are associated in the context of other shared coexpressed genes; (2) Decomposition of the SCCM using a novel heuristic algorithm termed "Triple-Link", which searches the highest connectivity in the SCCM, and then uses two connected TF as a primer for growing a TF cluster with a number of linking criteria. We applied TF-Cluster to microarray data from human stem cells and Arabidopsis roots, and then demonstrated that many of the resulting TF clusters contain functionally coordinated TFs that, based on existing literature, accurately represent a biological process of interest.

CONCLUSIONS

TF-Cluster can be used to identify a set of TFs controlling a biological process of interest from gene expression data. Its high accuracy in recognizing true positive TFs involved in a biological process makes it extremely valuable in building core GRNs controlling a biological process. The pipeline implemented in Perl can be installed in various platforms.